- Open Access
The phospholipid transporter PITPNC1 links KRAS to MYC to prevent autophagy in lung and pancreatic cancer
Molecular Cancer volume 22, Article number: 86 (2023)
The discovery of functionally relevant KRAS effectors in lung and pancreatic ductal adenocarcinoma (LUAD and PDAC) may yield novel molecular targets or mechanisms amenable to inhibition strategies. Phospholipids availability has been appreciated as a mechanism to modulate KRAS oncogenic potential. Thus, phospholipid transporters may play a functional role in KRAS-driven oncogenesis. Here, we identified and systematically studied the phospholipid transporter PITPNC1 and its controlled network in LUAD and PDAC.
Genetic modulation of KRAS expression as well as pharmacological inhibition of canonical effectors was completed. PITPNC1 genetic depletion was performed in in vitro and in vivo LUAD and PDAC models. PITPNC1-deficient cells were RNA sequenced, and Gene Ontology and enrichment analyses were applied to the output data. Protein-based biochemical and subcellular localization assays were run to investigate PITPNC1-regulated pathways. A drug repurposing approach was used to predict surrogate PITPNC1 inhibitors that were tested in combination with KRASG12C inhibitors in 2D, 3D, and in vivo models.
PITPNC1 was increased in human LUAD and PDAC, and associated with poor patients’ survival. PITPNC1 was regulated by KRAS through MEK1/2 and JNK1/2. Functional experiments showed PITPNC1 requirement for cell proliferation, cell cycle progression and tumour growth. Furthermore, PITPNC1 overexpression enhanced lung colonization and liver metastasis. PITPNC1 regulated a transcriptional signature which highly overlapped with that of KRAS, and controlled mTOR localization via enhanced MYC protein stability to prevent autophagy. JAK2 inhibitors were predicted as putative PITPNC1 inhibitors with antiproliferative effect and their combination with KRASG12C inhibitors elicited a substantial anti-tumour effect in LUAD and PDAC.
Our data highlight the functional and clinical relevance of PITPNC1 in LUAD and PDAC. Moreover, PITPNC1 constitutes a new mechanism linking KRAS to MYC, and controls a druggable transcriptional network for combinatorial treatments.
KRAS is a driver oncogene in several epithelial tumours. In lung and pancreatic ductal adenocarcinoma (LUAD and PDAC), KRAS mutation frequency reaches ~ 25 and 90% of cases respectively [1, 2]. These cancers share a notable dependency on aberrant KRAS expression through activation of canonical ‘‘proximal’’ effectors, mainly the RAF-MEK-ERK and the PI3K-AKT-mTOR pathways [3, 4]. Consequently, BRAF, MEK1/2, PI3K or mTOR inhibitors were developed and progressed to clinical trials. However, such inhibitors have had limited or no impact on cancers bearing KRAS mutations. Activation of KRAS canonical pathways translates into transcriptomic changes that involve far many ‘‘distal’’ effectors, whose defined role in oncogenesis is less understood. The discovery of transcriptomic changes involving functional effectors in KRAS-driven cancers is thus a relevant question in the field that may lead to the identification of molecular targets for novel therapeutic strategies. Our group previously reported a KRAS signature upregulated across KRAS-driven cancers . Follow up studies confirmed the relevant role of a member of the signature, the transcription factor FOSL1, in LUAD and PDAC  as well as in cholangiocarcinoma . These data suggest that additional genes pertaining to the KRAS signature could have a functional role in mutant (mut) KRAS cancers.
Traditionally regarded as integral components of biological membranes, phospholipids play a relevant role as signalling elements in normal homeostasis and cancer . For instance, phosphatidic acid (PA) functions as a docking site for the selective recruitment of effector proteins to local cell membrane compartments that are involved in transducing signals . Also, phosphatidylinositol (PI) has a central role in the regulation of PI3K-mediated oncogenesis as a precursor to phosphatidylinositol (3,4,5)-trisphosphate (PIP3) . In the KRAS setting, mut Kras increases the levels of phospholipids (i.e., phosphoinositide derivatives) to foster oncogenesis . Conversely, depletion of phospholipid (i.e., phosphatidylserine, PS) trafficking or membrane content compromises KRAS membrane association and subsequent activation [11,12,13]. Hence, it is plausible that proteins regulating the distribution and availability of phospholipids may contribute to the KRAS oncogenic phenotype. A recent study describing the relevance of the PS transport proteins ORP5 and ORP8 to support KRAS oncogenic function in PDAC  strongly favours this idea and warrants further investigation of additional members within the phospholipid transport machinery.
The phosphatidylinositol transfer protein (PITP) family participates in phospholipid transport between cell membranes . Class I PITPs, PITPα (PITPNA) and PITPβ (PITPNB) bind PI or phosphatidylcholine (PC). PITPNC1, which is classified as a Class II PITP, binds PI and PA instead of PC . PITPNC1 was originally reported as a gene amplified in human breast cancer and over-expressed in breast, colon and melanoma metastasis, where it fosters the pro-metastatic phenotype via secretion of pro-invasive and pro-angiogenic mediators . Subsequent studies characterizing PITPNC1 expression in human specimens revealed its association with advanced clinical stage and poor prognosis in gastric cancer , and with radio-resistance in rectal cancer . Nonetheless, the information about PITPNC1 in cancer is far from being completed, furthermore in the context of dominant oncogenes such as KRAS.
In this study, we uncover PITPNC1 as a KRAS-dependent gene with functional implications in LUAD and PDAC, in part by the unanticipated regulation of autophagy by modulating mTOR localization via cMYC (hereinafter referred to as MYC). Most importantly, PITPNC1 controls a druggable transcriptome that offers opportunities for therapeutic intervention in both cancers.
Materials and methods
Human LUAD cell lines (wild-type -wt- KRAS: H2126, H1568; mut KRAS: A549, H23, H2009, H1792, H358, H2347), PDAC cell lines (mut KRAS: PATU8902, HPAFII, PANC1, MiaPaca2) or normal human pancreatic duct epithelial (HPDE) cells, H6C7, were used. LUAD and PDAC cells were grown in fully supplemented RPMI1640 or DMEM media (Gibco) respectively. H6C7 were grown in serum-free Keratinocyte medium (Gibco). Cell lines were authenticated by the Genomics Unit at CIMA, using Short Tandem Repeat profiling (AmpFLSTR® Identifiler® Plus PCR Amplification Kit). Hras−/−; Nras−/−; Kraslox/lox; RERTnert/ert (hereinafter referred to as Kraslox/lox) mouse embryo fibroblasts (MEFs) were grown in fully supplemented DMEM (Gibco) and have been described previously . Mycoplasma test was performed in all cell lines every other week using the MycoAlert Mycoplasma Detection Kit (LONZA). Only mycoplasma-free cell lines were used.
Specific shRNAs oligonucleotides against PITPNC1 (sh6: TRCN0000059479, sh7: TRCN0000059481), MYC (sh42 TRCN0000039642, sh89 TRCN0000010389) were annealed and cloned into a pLKO.1 lentiviral vector (Addgene #10,878). The tet-pLKO-shKRAS (TRCN0000033260) and shRNA against green fluorescence protein (GFP) were already published . For CRISPR knockout experiments, specific sgRNAs against LKB1 (GTACTCCATCACCATATACG and CTTCAAGGTGGACATCTGGT), Lkb1 (GACACTCAAGATCTCCGACCT), Keap1 (GGGTTCGGTTACCGTCCTGCG), and Trp53 (GAACAGATCGTCCATGCAGTG) were annealed and cloned into lentiCRISPRv2 (Addgene #52,961) or lentiCRISPRv2-mCherry (Addgene #99,154) plasmids. Wt (4B) and mut (G12D, G12C and G12V) KRAS cDNAs were from the RAS Initiative (Addgene, Kit #1,000,000,089) and were described previously . The pBABE-puromycin PITPNC1-flagged plasmid was described earlier . For MEFs experiments, KRASG12C, KRASG12D, KRASG12V, KRASG12R, KRASG12S, KRASG13D and KRASQ61H retroviral plasmids were created by point mutagenesis from pBABE HA-tagged KRASWT plasmid (provided by Channing Der, Addgene plasmid # 75,282). A MYC cDNA in a pDONR221 vector was provided by Alejandro Sweet-Cordero (University of California San Francisco, USA) and cloned into a pLenti6/V5-DEST using the Gateway system (Thermofisher).
Virus production and infection
KRASMUT-expressing retroviruses were generated by co-transfection of pBABE plasmids together with pAmpho plasmid into HEK293T cells using FuGENE HD Transfection Reagent (Promega). The retroviruses were transduced into Kraslox/lox MEFs followed by 2 weeks of puromycin selection (1 μg/mL) in complete DMEM medium. To obtain Kras null-KRASMUT clones, cells were then cultured in the presence of 4-hydroxytamoxifen (4OHT) (Sigma, 600 nM) for another 2 weeks in order to achieve complete deletion of endogenous Kras allele. The remaining retrovirus and lentivirus were produced as previously described , then filtered and applied directly to cells for infection at a MOI lower than 1. Selection was done with puromycin, neomycin, or blasticidin (Sigma).
Quantitative PCR (qPCR) analysis
mRNA analysis was done as published earlier . Primers sequences are listed as Suppl. Material.
Immunoblot analysis was performed as previously described . Antibody information is found as Suppl. Material.
Cell proliferation assay
Cell proliferation assays were done as described earlier . A detailed explanation is found as Suppl. Material.
Clonogenic assays were performed as previously described . A detailed explanation is found as Suppl. Material.
Drug combination studies in vitro
Cell lines were plated at a density ranging from 300 to 2,000 cells in 96-well plates, treated on the following day with single drugs or combination, and cultured for 5 days. At day 5, cells were fixed with 4% formaldehyde (Panreac) for 15 min at RT, stained with crystal violet solution (Sigma-Aldrich) (1% crystal violet in H2O) for 15 min. Relative growth was quantified by measuring absorbance at 570 nm after crystal violet dissolution with 10% acetic acid. SynergyFinder software (https://synergyfinder.fimm.fi/) was used to determine the potential synergism of two drugs. Bliss score values > 10 were considered synergistic.
Long-term drug combination assays
Long-term drug combination assays were performed and analysed as previously described . A detailed explanation is found as Suppl. Material.
3D culture assays
3D culture assays were done as described previously . A detailed explanation is found as Suppl. Material.
Cell cycle and apoptosis assays
Cell cycle and apoptosis assays were performed as previously described . A detailed explanation is found as Suppl. Material.
The Connectivity Map (https://clue.io/) was used to predict genes or pharmacological compounds able to phenocopy a PITPNC1-knockdown signature induced by two specific shRNAs.
SP600125 (JNKi) and hydro-chloroquine (C6628) were from Sigma; BIX02189 (MEK5i) was from Tocris; Trametinib (MEKi), GSK2126458 (PI3Ki), Fedratinib (JAK2i), BI2536 (PLK1i), MG-132 and Sotorasib (KRASG12Ci) were from MedChemExpress.
All experiments in mice were performed according to the institutional Animal Care Committee (CEEA) of the University of Navarra under the protocols CEEA #057–18 approved by the regional Government of Navarra. A detailed description of the mouse experiments is provided as Suppl. Material.
Microscope image acquisition
Human tissue slides were scanned at 40 × magnification and digitalized using the Aperio Scan-Scope XT Slide Scanner (Aperio Technologies). Mouse slides were scanned at 40 × magnification and acquired with Aperio CS2 Leica Biosystems.
Immunofluorescence to determine mTOR localization was done using anti-mTOR (#2983 7C10 CST) and AlexaFluor 488-conjugated anti-LAMP1 antibody (clone H4A3, 328,609 Biolegend). A detailed explanation of the protocol and analysis is provided as Suppl. Material.
RNA sequencing (RNAseq) and analysis
Low full-length RNA seq libraries were prepared by adapting the Smart-seq 2 protocol to 1 ng of RNA . RNA sequencing data analysis was performed using the following workflow: (1) the quality of the samples was verified using FastQC software (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/); (2) the alignment of reads to the human genome (hg38) was performed using STAR ; (3) gene expression quantification using read counts of exonic gene regions was carried out with featureCounts ; (4) the gene annotation reference was Gencode v35 ; and (5) differential expression statistical analysis was performed using R/Bioconductor . Data are publicly available in GEO database with the accession number GSE205767. A detailed description of the RNA-seq analysis is provided as Suppl. Material.
Copy number variation data from the The Cancer Genome Atlas (TCGA) LUAD dataset  were downloaded and analyzed with GISTIC2 .
The association between loss of heterozygosity at the KRAS locus and PITPNC1 gene expression in mut KRAS LUAD was assessed in R using the ploidy estimated by ABSOLUTE and RNAseq batch corrected matrix published in Hoadley et al . The statistical significance was measured by a Wilcoxon’s test.
PITPNC1 expression profile
The expression profile of PITPNC1 gene was studied in the TCGA LUAD transcriptome dataset (downloaded from https://tcga-data.nci.nih.gov/tcga/tcgaHome2.jsp) to assess differential expression in patients. Data processing and statistical analyses were performed as previously described .
Survival analyses were conducted on both PITPNC1 gene and a PITPNC1 gene set in the TCGA  and Shedden et al  data sets and in the International Cancer Genome Consortium (ICGC)  and TCGA PDAC data sets. Log-rank test was used to calculate differences in Kaplan–Meier curves . For gene set studies, a summation of all the genes for a particular sample was calculated as previously described . Survival analyses were done with R  and p-values < 0.05 were statistically significant.
PITPNC1 gene set enrichment analyses
GSEA analyses were performed using the dPITPNC1 gene set. Enrichment was studied in samples with inhibition of KRAS expression (GSE196596 and GSE103021), in LUAD patients with KRAS mutation (GSE36133, GSE12667, GSE31210, and GSE26939), in PDAC patients (GSE15471 and GSE16515) or in PDAC models (GSE32277). Data processing of each experiment was previously described . FDR < 0.05 was considered statistically significant and the results were represented using GSEA Multi-sample Running Enrichment plots (https://github.com/GryderArt/VisualizeRNAseq).
Sample size was chosen using www.biomath.info/power/ttest.htm or based on similar experiments previously published by the authors. For comparison of two groups, normality (Shapiro–Wilk test) and variance (Levene test) was assessed. Groups with normal distribution followed a t-test. Non-normal samples were analysed using a Mann–Whitney test (equal variances) or a Median test (unequal variances). For comparison of more than two groups, a residual test was performed to study normality and Levene test assessed homoscedasticity. ANOVA, Brown Forsythe, Kruskal Wallis or Median tests were performed depending on data distribution. A post-hoc test (Dunnet or Bonferroni) explored paired comparisons. All analyses were two-tailed. Error bars correspond to either standard deviation (S.D., n < 8) or standard error of the mean (S.E.M, n ≥ 8). Survival analyses were done using the Log-rank test. Statistical analyses were done with GraphPad software v8.
PITPNC1 is regulated by KRAS oncogene and predicts poor survival in LUAD and PDAC
To identify novel effectors with relevance to KRAS oncogenesis, we used a gene signature derived from experimental models expressing mut KRAS (n = 41 genes) reported by our group . We queried this gene signature against TCGA data set by comparing gene expression profiles of wt and mut KRAS LUAD patients. Among the 15 differentially expressed genes, SPRY4, DUSP6, CCND1, PHLDA1, DUSP4, and PITPNC1 were the most robustly upregulated in KRAS-mutated patients (p < 0.0001) (Fig. 1A). We focused our attention on the phosphatidylinositol cytoplasmic transfer protein 1 (PITPNC1) which, unlike the other genes, had not been previously linked to KRAS oncogene biology. PITPNC1 was also upregulated when compared to normal lung tissue, indirectly suggesting a link to the tumour phenotype (Fig. 1B). PITPNC1 mRNA increase was dependent on KRAS oncogene expression, since the presence of other dominant oncogenic drivers (e.g. BRAF or EGFR) did not affect PITPNC1 transcript levels (Suppl. Figure 1A and B). PITPNC1 upregulation was not due to differential PITPNC1 amplification in mut vs wt KRAS LUAD patients either (p = 0.815) (Fig. 1C). Thus, the differential transcriptional regulation of PITPNC1 may be a consequence of aberrant KRAS activation. We further tested PITPNC1’s clinical role in human cancer by performing survival analysis in LUAD patients. High PITPNC1 expression was associated with poor overall survival in mut KRAS patients but not in wt (Fig. 1D). Since PITPNC1 was part of a mut KRAS signature that included genes with a role in LUAD and PDAC, we studied human PDAC specimens. Notably, high PITPNC1 was also a worse prognosis marker in PDAC (Fig. 1E).
The clinical data led us to test the connection between KRAS and PITPNC1 via genetic gain- and loss-of-function experiments in lung and pancreas cellular models. Overexpression of mut KRAS (G12D, G12C and G12V) in wt KRAS LUAD cells (H2126 and H1568) and in immortalized normal human pancreatic duct epithelial cells (H6C7) increased PITPNC1 protein and mRNA levels (Fig. 1F). Such PITPNC1 upregulation was also observed in human LUAD patients with different KRAS mutations (Suppl. Figure 1C). Conversely, KRAS inhibition in LUAD (A549, H2009, H1792) and PDAC (PATU8902, HPAFII) cells using a specific shRNA decreased PITPNC1 protein (Fig. 1G, and Suppl. Figure 1D and E). PITPNC1 was consistently downregulated upon inactivation of MEK1/2 and JNK1/2 in both LUAD (A549, H2009) and PDAC (HPAFII) cells (Fig. 1H and Suppl. Figure 1F), indicating a regulation by KRAS through different effector pathways. Notably, PITPNC1 was the unique member of the PITP family controlled by KRAS, as the expression of PITPNA, PITPNB, PITPNM1, PITPMN2 and PITPNM3 did not change upon KRAS genetic modulation (Suppl. Figure 1G-I).
In addition to KRAS activating mutations , an imbalance between wt and mut KRAS alleles can influence cancer cells’ fitness, expression profile and therapy response in LUAD and PDAC [37,38,39,40]. Thus, we investigated PITPNC1 levels in relationship to KRAS dosage. First, we used Kraslox/lox MEFs expressing different KRAS mutations to study PITPNC1 in the context of loss-of-heterozygosity (LOH) . Similar to human cell lines, exogenous expression of the various KRAS mutations increased PITPNC1 expression (Suppl. Figure 1J). Notably, Cre-excision of the wt allele in MEFs via 4-OHT treatment reduced PITPNC1 levels in all mutants but G12V (Fig. 1I). Such decrease was also observed in mut KRAS human samples of the LUAD TCGA data set (p = 0.047) (Fig. 1J). Second, we assessed the impact of mut KRAS amplification on PITPNC1 expression in the LUAD data set. However, no significant differences were found (Suppl. Figure 1K). These data may indicate that PITPNC1 represents a functional node downstream of KRAS integrating signals from receptor tyrosine kinases which become activated upon mut KRAS expression and require wt KRAS for downstream signalling.
Given the relevance of concurrent mutations in mut KRAS LUAD prognosis and response to therapy , we explored the association of PITPNC1 expression with prevalently mutated tumour suppressor genes (TSGs). LKB1 mutations were mostly found in mut KRAS with high PITPNC1 expression (p = 0.005) while ARID1A mutations appeared mostly in low PITPNC1-expressing tumours (p = 0.0208) (Suppl. Figure 2A). These results led us to test the impact of LKB1 mutations on PITPNC1 expression. CRISPR/Cas9-based LKB1 knockout in KRAS-mutated LUAD cells (H2009) enhanced PITPNC1 expression (Suppl. Figure 2B). This finding was further recapitulated in mouse LUAD cell lines driven by mut Kras (KLA and LKR10) upon LKB1 abrogation with specific sgRNAs (Suppl. Figure 2C). Thus, PITPNC1 is regulated by KRAS through MEK1/2 and JNK1/2 signalling pathways, and its expression may be exacerbated by LKB1 loss.
PITPNC1 inhibition reduces cell proliferation in vitro and impairs tumour growth in vivo in LUAD and PDAC
To characterise the functional role of PITPNC1, genetic depletion using two independent shRNAs, one of which had been previously validated via rescue experiments , was carried out in a panel of LUAD (n = 6) and PDAC (n = 4) cell lines (Fig. 2A). PITPNC1 inhibition consistently reduced cell proliferation of all cell lines (Fig. 2B). Likewise, a decreased colony-forming capacity was also observed in both tumour types (Fig. 2C). However, we did not find a consistent effect on apoptosis in PITPNC1-depleted cells (Suppl. Figure 3A and B).
Next, we investigated if PITPNC1 is necessary for KRAS-driven tumourigenesis in vivo. First, LUAD cell lines infected with PITPNC1 shRNAs were subcutaneously injected in immunocompromised mice. PITPNC1 knocked-down cells generated tumours of a smaller volume and weight than controls (Fig. 2D-F). PITPNC1 abrogation in PDAC cells also impaired tumour growth, yielding lighter tumours (Fig. 2G-I). The effect of PITPNC1 loss in vivo was related to decreased tumour proliferation and enhanced cytotoxic activity in these models (Fig. 2J and K, and Suppl. Figure 3C).
Complementary to PITPNC1 inhibition experiments, the effect of its overexpression was also assessed in mut KRAS LUAD cell lines (A549 and H358) (Suppl. Figure 3D). No effect on colony formation was observed (Suppl. Figure 3E). Moreover, exogenous PITPNC1 did not confer a growth advantage in vivo when cells were injected subcutaneously in immunodeficient mice (Suppl. Figure 3F-K).
Mut KRAS LUAD harbouring inactivating LKB1 mutations display poor prognosis, in part due to an enhanced metastatic potential [41, 42]. Since increased PITPNC1 expression was observed upon LKB1 loss, we explored PITPNC1 overexpression in the metastatic setting. A549 cells were first constructed to express luciferase and transduced with a PITPNC1-expressing or a control vector. We used a mouse model of lung colonization in vivo where cancer cells initially seed in the lungs after intravenous injection (~ 10–15 min post-injection) (Fig. 2L). Subsequent bioluminescence monitoring revealed that, on week 1, PITPNC1-expressing cells colonized the lung more efficiently. Interestingly, while the bioluminescence signal of the two groups become closer by week 2, it increased in cancer cells over-expressing PITPNC1 at week 3 and 4 (Fig. 2M). Macroscopic and microscopic analysis of tissues at endpoint revealed a higher number of metastatic liver foci in the group of mice injected with PITPNC1-overexpressing cells while the lung tumor burden was similar (Fig. 2N-Q), suggesting that secondary metastasis to the liver contribute to the distinct bioluminescence signal. No differences in the migratory capacity in vitro or the metastatic tropism in vivo of control and PITPNC1-overexpressing cells were detected that could explain these findings (Suppl. Figure 3L-N), suggesting the involvement of heterotypic interactions as described in breast cancer . Thus, PITPNC1 upregulation contributes to the metastatic phenotype of mut KRAS LUAD.
A PITPNC1 gene signature features KRAS-regulated genes and predicts poor survival in LUAD and PDAC
To get a better understanding of PITPNC1 as a KRAS effector, we interrogated the transcriptome of KRAS-mutated LUAD cells (A549) after PITPNC1 inhibition with two shRNAs. A total of 429 genes were found differentially expressed (logFC ± 1, B > 0) with regard to control cells (Fig. 3A). The downregulated PITPNC1 gene signature (dPITPNC1 GS; n = 233 genes), which a priori would contain transcriptional targets whose overexpression fosters the oncogenic phenotype, was used. This signature was queried against two independent data sets where genetic or pharmacological blockade of KRAS, via a tet-inducible KRAS shRNA or the KRASG12C ARS160 inhibitor respectively, was carried out. A consistent enrichment of the dPITPNC1 GS in genes repressed upon KRAS inhibition was found (Fig. 3B). To expand these findings to the pancreas setting, we took advantage of gene expression data from cancer cell lines (iKrasC) and xenograft tumours (iKrasT) derived from an inducible genetically engineered mouse (GEM) model of Kras-driven PDAC in which doxycycline administration activates expression of a mut Kras allele  (Fig. 3C). In both data sets, a large overlap of the dPITPNC1 GS with genes decreased after oncogenic KRAS inactivation was found, suggesting that multiple PITPNC1-regulated genes are part of the KRAS signalling pathway.
To test the PITPNC1-regulated genes in a more clinically relevant setting, we performed GSEA using human LUAD data sets (n = 4) with information on the KRAS mutational status. A general enrichment of the dPITPNC1 GS was found in LUAD tumours harbouring KRAS mutations compared to those with native alleles (Fig. 3D). Likewise, we found a strong enrichment of the dPITPNC1 GS in human PDAC samples with regard to normal pancreas in two data sets (Fig. 3E). Additional analysis of genes whose expression was diminished in response to PITPNC1 were recurrently present in the leading edges of the previously investigated data sets was done by qPCR. A dramatic reduction in mRNA expression was detected for all genes (Suppl. Figure 4A), validating the RNAseq data.
We next explored the clinical relevance of the dPITPNC1 GS. We observed that high dPITPNC1 GS levels were associated with the LUAD and PDAC patient subgroup with the worst prognosis (Fig. 3F-I). Analysis of the signature in the context of tumour stage revealed no differences in either tumour type (Suppl. Figure 4B and C). Likewise, no significant changes in two of the main patient subgroups of LUAD (mut KRAS/P53-mutated and mut KRAS/LKB1-mutated) and PDAC (classical and basal), which display differential outcome, response to therapy and gene expression profiles [41, 44, 45], were found (Suppl. Figure 4D and E). Collectively, these results indicate that PITPNC1 controls the expression of a gene signature with clinical implications for KRAS-mutated tumours.
PITPNC1 loss induces a G1 phase arrest and MYC downregulation
To expand our understanding of PITPNC1’s functional role in KRAS-driven oncogenesis, we performed Gene Ontology analysis to infer the biological pathways (BP) related to PITPNC1-regulated genes. First, the dPITPNC1 GS was used as input. The top BP included general cell cycle, sodium ion transmembrane transport, regulation of hormone levels or establishment/maintenance of cell polarity (Fig. 4A). These findings prompted us to inquiry about the impact of PITPNC1 loss on the cell cycle. We found a consistent G1 arrest and S phase decrease across all LUAD cell lines (Fig. 4B and Suppl. Figure 5A). These observations were extended to the PDAC setting (Fig. 4C), suggesting the regulation of common cellular mechanisms across mut KRAS tumours.
BP analysis also featured a MYC active pathway, which led us to test MYC expression in mut KRAS LUAD and PDAC cells with depleted PITPNC1. We found an overt MYC downregulation across all cell lines studied, which mainly occurred at the protein level (Fig. 4D-F), positioning MYC downstream of PITPNC1 and providing a direct link to the KRAS pathway. Such MYC downregulation was recapitulated upon PITPNC1 inhibition in vivo (Suppl. Figure 5B). MYC cooperates with oncogenic RAS to regulate G1 to S phase transition of cell cycle , a phenotype observed in PITPNC1-depleted cells. Indeed, MYC inhibition using two specific shRNAs revealed a G1 arrest similar to that found in cells with PITPNC1 loss (Suppl. Figure 5C). This mechanism involves repression of various cyclin kinase inhibitors, such as CDKN1B (p27) and CDKN1C (p57), and activation of E2F transcription factors , a link sustained in our experimental models (Suppl. Figure 5D). These observations led us to investigate the molecular consequences of PITPNC1 loss on the cell cycle. Detailed analysis of transcriptomics data showed upregulation of p27 and p57, and downregulation of E2F1 (Suppl. Figure 5E). These results were validated using qPCR and Western blot analyses in independent samples (Fig. 4G and H, and Suppl. Figure 5F and G).
To investigate if PITPNC1 regulates cell cycle through MYC, exogenous MYC was overexpressed in PITPNC1-depleted cells. No rescue of the proliferative phenotype was found, most likely because MYC levels were still low even in the overexpressing cells (Suppl. Figure 5H and I). Notably, blocking the proteasome activity with the specific inhibitor MG132 rescued MYC expression, suggesting post-translational regulatory mechanisms (Fig. 4I). This prompted us to scan the PITPNC1-knockdown RNAseq data for potential kinases involved in MYC protein regulation, and found downregulation of AURKA and PLK1 (Suppl. Figure 5J), two kinases previously reported to stabilize MYC protein via direct phosphorylation [48, 49]. Notably, only PLK1 was consistently decreased across the various LUAD and PDAC cell lines upon PITPNC1 loss, with an expression pattern mimicking that of MYC protein (Fig. 4J and K). We next tested the possibility that PLK1 regulates MYC protein. Using the PLK1 inhibitor BI-2536, we found reduced MYC protein expression that is rescued by proteasome inhibition (Fig. 4L). Thus, PITPNC1 may be regulating MYC protein expression in part by PLK1. Taken together, these results suggest that PITPNC1 represents a functional link that connects oncogenic KRAS to MYC.
PITPNC1 controls mTOR localization via MYC to prevent autophagy
To complement the previous findings, we explored those genes upregulated upon PITPNC1 abrogation (i.e., uPITPNC1 GS). The top 5 BPs of the GO analysis involved P53 transcriptional gene network, regulation of mTORC1 signalling, antigen processing and presentation of endogenous peptide antigen via MHC class I via ER pathway, natural killer cell-mediated toxicity, and genotoxicity pathway (Fig. 5A). We focused on mTOR as it is an effector of the PI3K pathway that can function within the KRAS signalling network. Enriched genes in the regulation of mTORC1 signalling BP feature included CASTOR1, RRAGD, SESN1, SESN3, and GPR137C. Upregulation of SESN1, SESN2 and SESN3 was validated at the mRNA level using qPCR and the results confirmed in additional PITPNC1-depleted cells (H2009 and HPAFII) (Fig. 5B and C, and Suppl. Figure 6A). These results suggested that activation of the mTOR pathway is altered upon PITPNC1 inhibition.
SESTRINS (SESN1-3) inactivate GATOR2 to inhibit mTOR activity, constraining the localization of mTOR to the lysosome where it gets activated [50, 51]. Thus, we investigated mTOR localization in response to PITPNC1 abrogation by immunofluorescence. A549, H2009 and HPAFII control cells showed mTOR activation, as inferred from the overlapping signal with the lysosome marker LAMP1. However, this colocalization was impaired when PITPNC1 was inhibited (Fig. 5D-G, and Suppl. Figure 6B and C). This modification occurred without changes in mTOR protein abundance (Suppl. Figure 6D), suggesting that PITPNC1 controls mTOR lysosomal recruitment.
A close visualization of the lysosomes in PITPNC1-depleted cells revealed increased number and size compared to PITPNC1-proficient ones (Fig. 5H and I, and Suppl. Figure 6E). Expansion of the lysosomal compartment or lysosomal biogenesis has been related to enhanced autophagy . mTOR functions as a counter-regulator of autophagy [53, 54]. Thus, we analysed the level of the autophagy marker LC3-II, which tightly correlates with the number of autophagosomes/autophagolysomes . Increased LC3-II was observed in cell lines lacking PITPNC1 (Fig. 5J). In keeping with autophagy induction, downregulation of gene signature featuring autophagy and lysosome biogenesis  was also found in PITPNC1-inhibited cells (Suppl. Figure 6F).
Enhanced LC3-II expression could indicate either upregulation of autophagic flux (i.e., autophagosome formation) or blockade of autophagic degradation . To confirm the underlying mechanism, we compared changes in LC3-II under the presence of the lysosomal protease inhibitor hydrochloroquine, which accumulates within lysosomes leading to lysosome neutralization and the inhibition of autophagic flux/autophagosome formation . Hydroxichloroquine treatment elicited a further accumulation of LC3-II (Suppl. Figure 6G and H), indicating that PITPNC1 inhibition enhances autophagic flux. This mechanism occurred without activation changes in S6K and 4EBP1 (Suppl. Figure 6I).
MYC suppresses autophagy in B cell lymphomas by antagonizing the function of TFEB transcription factors , raising the possibility that PITPNC1 could control autophagy through MYC in LUAD and PDAC. To address this possibility, we first tested if MYC regulated autophagy in our experimental systems. MYC inhibition by specific shRNAs induced LC3II/I ratio in LUAD and PDAC cell lines (Fig. 5K). Autophagy induction was associated with reduced mTOR localization to lysosomes (Suppl. Figure 7A-D). This was associated with an increase in number and size of lysosomes (Suppl. Figure 7E and F). To define how MYC regulates autophagy, we tested if MYC could be transcriptionally controlling the negative regulators of mTOR localization, SESTRIN1-3, which are downregulated after PITPNC1 inhibition. qPCR analysis of MYC-depleted cells showed that MYC inhibition significantly enhanced their expression (Suppl. Figure 7G and H), positioning SESTRINS downstream of MYC. Taken together, these observations suggest that PITPNC1 controls mTOR activity via MYC to prevent autophagy. A proposed model for the role of PITPNC1 in LUAD AND PDAC KRAS-driven tumours is depicted in Fig. 5L.
JAK2 inhibitors reverse the expression of a PITPNC1-regulated transcriptome and synergize with Sotorasib
Given the lack of pharmacological tools to inhibit PITPNC1 and aiming to increase the translational value of our findings, we followed a drug repurposing strategy that predicts compounds capable of reversing the expression profile of the PITPNC1-regulated transcriptome. The top 200 up and down differentially expressed genes obtained after PITPNC1 knockdown (logFC ± 1, B > 0) were used as input and a repurposing score > 90 was used as cut-off. The top 5 drug families predicted to reverse the PITPNC1 transcriptome were JAK, HDAC, DNA synthesis, bromodomain and DNA dependent protein kinase inhibitors (i). Additional drug families scoring in this analysis were PI3Ki, mTORi or MEKi, known downstream effectors of KRAS oncogene (Fig. 6A). The same drug repurposing approach applied to a KRAS-dependent transcriptome uncovered common drug families (Suppl. Figure 8A), consistent with the overlap of PITPNC1- and KRAS-regulated genes. JAK inhibitors, particularly those against JAK2, scored highest in both repurposing studies and were selected for downstream analyses.
To investigate if JAK2 inhibitors would phenocopy to some extent the effect caused by PITPNC1 abrogation, LUAD and PDAC cell lines were treated with Fedratinib, a highly specific JAK2i. A gradual decrease of cell proliferation was observed with increasing concentrations of Fedratinib in all cell lines (Fig. 6B). IC50 values ranged from 1.2 to 2 µM. Of note, Fedratinib treatment induced MYC depletion, LC3 upregulation or both, partially recapitulating the PITPNC1 inhibition phenotype (Fig. 6C). Thus, JAK2 inhibitors could function as a surrogate tool of PITPNC1 depletion.
KRASG12C inhibitors have emerged as promising targeted agents for mut KRAS-driven tumours [59, 60], albeit the clinical data suggest that combinatorial strategies may be required for more durable antitumour responses [61,62,63]. A rational concept for combination therapies builds on maximal driver pathway inhibition . Thus, we tested the combination of the KRASG12Ci Sotorasib with Fedratinib, both approved by the Food and Drug Administration (FDA) [65, 66]. To do this, human KRASG12C LUAD and PDAC cell lines (n = 5) were treated alone or in combination with various concentrations < IC25 of each drug. A larger antiproliferative phenotype was elicited by the drug combination compared to individual drugs (Fig. 6D and Suppl. Figure 8B). SynergyFinder revealed that the dual treatment was synergistic (Suppl. Figure 8C). The drug combination was also studied in 3D organoid cultures, given the enhanced KRAS oncogene dependence observed compared to 2D cultures [67, 68]. Likewise, the dual treatment effected cell proliferation more largely than single drugs (Fig. 6E and F, Suppl. Figure 8D). Molecular characterization of downstream targets revealed activation of STAT3 in response to KRASG12C inhibition that was suppressed in the drug combination. Furthermore, specific apoptosis induction was also observed in the combined treatment, suggesting a cytotoxic effect (Fig. 6G).
Next, we assessed the effect of the Fedratinib-Sotorasib combination in the context of adaptive resistance. First, using a 10 day-treatment colony forming assay, where early adaptive, non-genetic resistant mechanisms are likely to arise in response to individual drugs, we found that both drugs restrict cell proliferation to a greater extent than drugs alone (Fig. 6H). Second, the dual combination was tested in cell lines that had been made resistant to Sotorasib through gradual treatment with increasing drug concentrations for over 1 month (H358SR and H23SR) (Macaya and Roman et al., under review). The combined treatment yielded an antiproliferative response that was higher than each drug in the two resistant cells and also synergistic (Fig. 6I and J, and Suppl. Figure 8E). Thus, a JAK2i potentiates KRASG12Ci’s effect in both treatment-naïve and KRASG12Ci-resistant cells.
Anti-tumour activity of combined JAK2 and KRASG12C inhibitors in vivo
To investigate the impact of the dual combination in a more physiologically representative system, we generated cell-derived xenograft (CDX) models from LUAD (H358) and PDAC (MiaPaca2) cell lines in immunodeficient mice (n = 10–12 group). Fedratinib was administered at 60 mg/kg (twice daily) and Sotorasib at suboptimal concentrations 10 mg/kg (daily) as described earlier . Treatment started when the average volume of H358- and MiaPaca2-derived tumours reached 80 or 90 mm3 respectively. In the H358 CDX model, tumours treated with Fedratinib displayed growth kinetics similar to the vehicle-treated group, while those treated with Sotorasib were much lower, with a slight volume increase over time. Notably, concomitant drug administration led to generalized tumour regressions (Fig. 7A-D). Analysis of cell proliferation (pH3) and apoptosis (CC3) 7 days post-treatment revealed a notable decrease and increase respectively upon combined treatment (Fig. 7E and F, and Suppl. Figure 9A). In the MiaPaca2 CDX model, Fedratinib had no impact on tumour growth as compared to control mice whereas Sotorasib delayed tumour growth. The drug combination had a more profound impact on tumour growth, with tumour volume barely changing from the start point. The overt effect of the drug combination on the tumour volume translated into significantly smaller tumours (Fig. 7G-J). Histological analysis 7 days post-treatment showed increased apoptosis but no changes on tumour proliferation after the dual treatment (Fig. 7K and L, and Suppl. Figure 9B). More importantly, the effect of the drug combination had no consequences on the mice weight (Suppl. Figure 9C and D).
Collectively, these results nominate the combination of JAK2 and KRASG12C inhibitors as a potential strategy to treat KRAS-driven LUAD and PDAC harboring G12C mutation.
Through a series of clinical, cellular, molecular, and in vivo analyses, our study provides the first evidence that the phospholipid transporter PITPNC1 stands as a functional KRAS effector in LUAD and PDAC, further confirming the critical role of the phospholipid transport machinery in cancers driven by the KRAS oncogene.
PITPNC1 functions as a phospholipid transporter that was originally linked to the metastatic phenotype in breast cancer . However, while subsequent studies using PITPNC1 staining of primary gastric and rectal cancer revealed an association with clinical stage and poor prognosis, and radio-resistance respectively , the involvement of PITPNC1 in the regulation of additional cellular and molecular mechanisms in cancer was still unknown. We provide data indicating that PITPNC1 consistently controls the cell cycle of LUAD and PDAC driven by oncogenic KRAS, in part due to the regulation of positive (E2F1) and negative (P27, P57) cell cycle modulators by the transcription factor MYC. MYC protein expression consistently decreases in the PITPNC1 loss phenotype, a mechanism that involves PLK1 downregulation, a kinase that stabilizes MYC protein via direct phosphorylation . Indeed, PLK1 or proteasome pharmacological inhibition rescues MYC expression. In addition to MYC protein stabilization via ERK1/2 phosphorylation [69,70,71] under normal conditions, or through ERK5 when cancer cells are treated with MEK1/2 or ERK1/2 inhibitors [72, 73], our data provides a new mechanism for MYC protein regulation by KRAS oncogene.
We find that PITPNC1 loss also exacerbates the antiproliferative effect by inducing autophagy across mut KRAS LUAD and PDAC. This function depends on the unexpected regulation of mTOR activity via control of its cellular localization, as PITPNC1-deficient mut KRAS cells display a dramatic mTOR delocalization from lysosomes. This mechanism seems to be dependent on MYC activity, as MYC depletion phenocopies the autophagy induction elicited by PITPNC1 inhibition. A plausible explanation to mTOR delocalization may be the consistent upregulation of SESTRIN1-3 observed in the RNAseq and qPCR data. SESTRINS have been reported to control mTORC1 localization and subsequent activity via negatively regulation of the amino acid sensing pathway upstream of mTORC1 through GATOR2 binding , a mechanism that is partially phenocopied by PITPNC1 loss. Alternatively, PITPNC1’s canonical role regulating PA availability as a lipid transporter could also explain mTOR activation. Indeed, phospholipase D-dependent accumulation of cellular PA after mitogenic stimulation is required for phosphorylation of mTOR downstream effectors S6K and 4EBP1 . While we cannot completely rule out that the canonical function of PITPNC1 may also influence mTOR activity, we have not detected consistent activation of S6K and 4EBP1 in PITPNC1-depleted conditions that would support such activation mechanism. In any case, our data highlight a novel mechanism whereby KRAS can regulate mTOR activation via MYC, what adds to the previous reports showing that mTOR activation by the KRAS pathway occurred via phosphorylation by the ERK-RSK axis [75,76,77].
The cellular and molecular data highlighting the contribution of PITPNC1 to KRAS oncogenesis are supported by complementary analyses using LUAD and PDAC clinical data. First, PITPNC1 is uncovered as a marker of poor prognosis in both tumour types. In tune with these data, a PITPNC1-downregulated gene signature showed an inverse correlation with overall survival. The consistent results across tumours suggest that a PITPNC1-regulated network represents a relevantly common signalling node in KRAS oncogenesis with clinical implications that could be exploited for therapeutic purposes.
Studies focused on PITPNC1 have been limited by the lack of pharmacological inhibitors, which precluded exploring its role as a molecular target in cancer or other diseases. Our study supports the use of gene expression-based drug repurposing to infer drugs that could recapitulate the gene expression network controlled by a particular gene. This approach unveiled the family of JAK inhibitors in general and JAK2i in particular among the top predicted drugs in addition to PI3K, mTOR, or MEK1/2 inhibitors, known downstream effectors of KRAS oncogene. The fact that mTOR inhibitors are predicted to control a PITPNC1 signature further supports the link between mTOR and the phospholipid transporter. Of note, drug predictions largely overlapped with those obtained using a KRAS-regulated gene signature, suggesting that JAK inhibitors may indeed function in mut KRAS cancer. In keeping with this possibility, LUAD and PDAC lines are sensitive to the JAK2i Fedratinib at the low micromolar range. Although JAK2 inhibitors have only been approved for the treatment of myeloproliferative disorders , several preclinical studies have shown the implication of the JAK/STAT pathway in solid tumours [79,80,81]. Corcoran et al. demonstrated how activation of the JAK/STAT pathway is critical for the maintenance and development of PDAC . Furthermore, the relevance of this pathway was also reported in NSCLC , pulmonary fibrosis , or colorectal cancer [85, 86]. Collectively, these findings nominate JAK2 inhibitors as potential drugs for the treatment of KRAS-driven cancers.
Given that single drugs elicit limited antitumour responses in KRAS-mutated cancers, in part due to resistance mechanisms, this study proposes the combination of JAK2i with the recently FDA-approved KRASG12Ci Sotorasib based on a notable synergistic effect in vitro and in vivo. Testing this drug combination is mechanistically supported by the fact that the JAK/STAT pathway can act as a mechanism of compensation to MAPK pathway inhibition treatment . Furthermore, concomitant JAK2 and MEK1/2 inhibition reprograms the cancer-associated fibroblast (CAF) and immune microenvironment to overcome resistance to anti-PD-1 therapy in PDAC . These observations warrant the exploration of the JAK2i-KRASG12Ci drug combination alone or in the context of immune checkpoint inhibitors. The fact that both drugs are already approved for clinical applications and known to be well-tolerated may facilitate their progression to clinical trials. Lastly, while our observations are restricted to combinations using a KRASG12Ci, it is tempting to speculate that upcoming KRAS inhibitors targeting alternative mutations (e.g. G12D or G12V) may also synergize with JAK2 inhibitors, opening a new avenue for the treatment of a larger fraction of LUAD and PDAC patients.
In conclusion, our work uncovers the phospholipid transporter PITPNC1 as a KRAS effector that controls central transcriptional and signalling nodes (i.e. MYC and mTOR) and unveils novel therapeutic strategies for KRAS-driven tumours in the context of a PITPNC1-regulated transcriptional network.
Availability of data and materials
The datasets generated and/or analysed during the current study are available through the referenced publications or at GEO website as described in the Methods section.
Eukaryotic translation initiation factor 4E binding protein 1
AT-rich interaction domain 1A
Aurora kinase A
V-raf murine marcoma viral oncogene homolog B1
Cytosolic arginine sensor for MTORC1 subunit 1
- dPITPNC1 GS:
Downregulated PITPNC1 gene signature
Dual specificity phosphatase 4
Dual specificity phosphatase 6
E2F transcription factor 1
Epidermal growth ractor receptor
GATOR complex protein 2
Genetically engineered mouse
International Cancer Genome Consortium
Janus kinase 2
C-Jun N-terminal kinase
Kelch-like ECH associated protein 1
Kirsten rat sarcoma viral oncogene homologue
Lysosomal associated membrane protein 1 or CD107a
Microtubule associated protein 1 light chain 3 alpha/beta
Liver kinase B1
Mitogen-activated protein kinase
Mitogen-activated protein kinase
Mammalian target of rapamycin kinase
CMYC proto-oncogene, bHLH transcription factor
Oxysterol Binding Protein Like 5
Oxysterol Binding Protein Like 8
Cyclin dependent kinase inhibitor 1B or CDKN1B
Tumour protein 53 (TP53)
Cyclin dependent kinase inhibitor 1C or CDKN1C
Pancreatic ductal adenocarcinoma
Pleckstrin homology-like domain family A member
Phosphatidyl inositol 3-kinase
Phosphatidylinositol transfer protein
Phosphatidylinositol transfer protein α
Phosphatidylinositol transfer protein β
Phosphatidylinositol transfer protein cytoplasmatic 1
Phosphatidylinositol Transfer Protein Membrane Associated 1
Phosphatidylinositol Transfer Protein Membrane Associated 2
Phosphatidylinositol Transfer Protein Membrane Associated 3
Polo-like kinase 1
Ribosomal protein S6 kinase
Signal transducer and activator of transcription
The Cancer Genome Atlas
Tumour suppressor genes
- uPITPNC1 GS:
Upregulated PITPNC1 gene signature
Cox AD, Fesik SW, Kimmelman AC, Luo J, Der CJ. Drugging the undruggable RAS: Mission possible? Nat Rev Drug Discov. 2014;13:828–51.
Johnson C, Burkhart DL, Haigis KM. Classification of KRAS-activating mutations and the implications for therapeutic intervention. Cancer Discov. 2022;12:913–23.
Drosten M, Barbacid M. Targeting the MAPK pathway in KRAS-driven tumors. Cancer Cell. 2020;37:543–50.
Stephen AG, Esposito D, Bagni RK, McCormick F. Dragging ras back in the ring. Cancer Cell. 2014;25:272–81.
Vallejo A, et al. An integrative approach unveils FOSL1 as an oncogene vulnerability in KRAS-driven lung and pancreatic cancer. Nat Commun. 2017;8:14294.
Vallejo A, et al. FOSL1 promotes cholangiocarcinoma via transcriptional effectors that could be therapeutically targeted. J Hepatol. 2021;75:363–76.
Sunshine H, Iruela-Arispe ML. Membrane lipids and cell signaling. Curr Opin Lipidol. 2017;28:408–13.
Zhao C, Du G, Skowronek K, Frohman MA, Bar-Sagi D. Phospholipase D2-generated phosphatidic acid couples EGFR stimulation to Ras activation by Sos. Nat Cell Biol. 2007;9:707–12.
Wong K-K, Engelman JA, Cantley LC. Targeting the PI3K signaling pathway in cancer. Curr Opin Genet Dev. 2010;20:87.
Valitutti S, Cucchi P, Colletta G, Di Filippo C, Corda D. Transformation by the K-RAS oncogene correlates with increases in phospholipase A2 activity, glycerophosphoinositol production and phosphoinositide synthesis in thyroid cells. Cell Signal. 1991;3:321–32.
Cho K, et al. Staurosporines disrupt phosphatidylserine trafficking and mislocalize Ras proteins. J Biol Chem. 2012;287:43573–84.
Zhou Y, et al. Membrane potential modulates plasma membrane phospholipid dynamics and K-Ras signaling. Science. 2015;349:873–6.
Kattan WE, et al. Targeting plasma membrane phosphatidylserine content to inhibit oncogenic KRAS function. Life Sci Alliance. 2019;2:e201900431.
Kattan WE, et al. Components of the phosphatidylserine endoplasmic reticulum to plasma membrane transport mechanism as targets for KRAS inhibition in pancreatic cancer. Proc Natl Acad Sci USA. 2021;118:e2114126118.
Cockcroft S. The diverse functions of phosphatidylinositol transfer proteins. Curr Top Microbiol Immunol. 2012;362:185–208.
Garner K, et al. Phosphatidylinositol transfer protein, cytoplasmic 1 (PITPNC1) binds and transfers phosphatidic acid. J Biol Chem. 2012;287:32263–76.
Halberg N, Sengelaub CA, Navrazhina K, Molina H, Uryu K. PITPNC1 recruits RAB1B to the Golgi network to drive malignant secretion. Cancer Cell. 2016;29:339–53.
Tan Y, et al. Adipocytes fuel gastric cancer omental metastasis via PITPNC1-mediated fatty acid metabolic reprogramming. Theranostics. 2018;8:5452–68.
Tan Y, et al. PITPNC1 fuels radioresistance of rectal cancer by inhibiting reactive oxygen species production. Ann Transl Med. 2020;8:126.
Drosten M, et al. Genetic analysis of Ras signalling pathways in cell proliferation, migration and survival. EMBO J. 2010;29:1091–104.
Vicent S, et al. Wilms tumor 1 (WT1) regulates KRAS-driven oncogenesis and senescence in mouse and human models. J Clin Invest. 2010;120:3940–52.
Valencia K, et al. The Mir181ab1 cluster promotes KRAS-driven oncogenesis and progression in lung and pancreas. J Clin Invest. 2020;130:1879–95.
Picelli S, et al. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc. 2014;9:171–81.
Dobin A, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Harrow J, et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 2012;22:1760–74.
Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102:15545–50.
Mermel CH, et al. GISTIC20 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011;12:R41.
Hoadley KA, et al. Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer. Cell. 2018;173:291-304.e6.
Director’s Challenge Consortium for the Molecular Classification of Lung Adenocarcinoma, et al. Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat Med. 2008;14:822–7.
Bailey P, et al. Genomic analyses identify molecular subtypes of pancreatic cancer. Nature. 2016;531:47–52.
Clark TG, Bradburn MJ, Love SB, Altman DG. Survival analysis part I: basic concepts and first analyses. Br J Cancer. 2003;89:232–8.
Vicent S, et al. Cross-species functional analysis of cancer-associated fibroblasts identifies a critical role for CLCF1 and IL-6 in non-small cell lung cancer in vivo. Can Res. 2012;72:5744–56.
Vallejo A, et al. An integrative approach unveils FOSL1 as an oncogene vulnerability in KRAS-driven lung and pancreatic cancer. Nat Commun. 2017;8:14294.
Ricciuti B, et al. Comparative analysis and isoform-specific therapeutic vulnerabilities of KRAS mutations in non-small cell lung cancer. Clin Cancer Res. 2022;28:1640–50.
Zhang Z, et al. Wildtype Kras2 can inhibit lung carcinogenesis in mice. Nat Genet. 2001;29:25–33.
To MD, Rosario RD, Westcott PMK, Banta KL, Balmain A. Interactions between wild-type and mutant Ras genes in lung and skin carcinogenesis. Oncogene. 2013;32:4028–33.
Ambrogio C, et al. KRAS dimerization impacts MEK inhibitor sensitivity and oncogenic activity of mutant KRAS. Cell. 2018;172:857-868.e15.
Mueller S, et al. Evolutionary routes and KRAS dosage define pancreatic cancer phenotypes. Nature. 2018;554:62–8.
Skoulidis F, et al. Co-occurring genomic alterations define major subsets of KRAS - mutant lung adenocarcinoma with distinct biology, immune profiles, and therapeutic vulnerabilities. Cancer Discov. 2015;5:860–77.
Ji H, et al. LKB1 modulates lung cancer differentiation and metastasis. Nature. 2007;448:807–10.
Kapoor A, et al. Yap1 activation enables bypass of oncogenic Kras addiction in pancreatic cancer. Cell. 2014;158:185–97.
Chan-Seng-Yue M, et al. Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution. Nat Genet. 2020;52:231–40.
Aung KL, et al. Genomics-driven precision medicine for advanced pancreatic cancer: early results from the COMPASS trial. Clin Cancer Res. 2018;24:1344–54.
Leone G, DeGregori J, Sears R, Jakoi L, Nevins JR. Myc and Ras collaborate in inducing accumulation of active cyclin E/Cdk2 and E2F. Nature. 1997;387:422–6.
Bretones G, Delgado MD, León J. Myc and cell cycle control. Biochim Biophys Acta. 2015;1849:506–16.
Dauch D, et al. A MYC-aurora kinase A protein complex represents an actionable drug target in p53-altered liver cancer. Nat Med. 2016;22:744–53.
Ren Y, et al. PLK1 stabilizes a MYC-dependent kinase network in aggressive B cell lymphomas. J Clin Invest. 2018;128:5517–30.
Chantranupong L, et al. The Sestrins interact with GATOR2 to negatively regulate the amino-acid-sensing pathway upstream of mTORC1. Cell Rep. 2014;9:1–8.
Chantranupong L, et al. The CASTOR Proteins Are Arginine Sensors for the mTORC1 Pathway. Cell. 2016;165:153–64.
Yim WW-Y, Mizushima N. Lysosome biology in autophagy Cell Discov. 2020;6:6.
Rabanal-Ruiz Y, Otten EG, Korolchuk VI. mTORC1 as the main gateway to autophagy. Essays Biochem. 2017;61:565–84.
Deleyto-Seldas N, Efeyan A. The mTOR-autophagy axis and the control of metabolism. Front Cell Dev Biol. 2021;9:655731.
Klionsky DJ, et al. Guidelines for the use and interpretation of assays for monitoring autophagy (4th edition). Autophagy. 2021;17:1–382.
Ravichandran M, et al. Coordinated transcriptional and catabolic programs support iron-dependent adaptation to RAS-MAPK pathway inhibition in pancreatic cancer. Cancer Discov. 2022;12:2198–219.
Mizushima N, Yoshimori T. How to Interpret LC3 Immunoblotting. Autophagy. 2007;3:542–5.
Fernandez MR, et al. Disrupting the MYC-TFEB circuit impairs amino acid homeostasis and provokes metabolic anergy. Cancer Res. 2022;82:1234–50.
Canon J, et al. The clinical KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. Nature. 2019;575:217–23.
Hallin J, et al. The KRASG12C inhibitor, MRTX849, provides insight toward therapeutic susceptibility of KRAS mutant cancers in mouse models and patients. Cancer Discov. 2020;10:54–71.
Hong DS, et al. KRASG12C inhibition with sotorasib in advanced solid tumors. N Engl J Med. 2020;383:1207–17.
Zhao Y, et al. Diverse alterations associated with resistance to KRAS(G12C) inhibition. Nature. 2021;599:679–83.
Awad MM, et al. Acquired resistance to KRASG12C inhibition in cancer. N Engl J Med. 2021;384:2382–93.
Boshuizen J, Peeper DS. Rational cancer treatment combinations: an urgent clinical need. Mol Cell. 2020;78:1002–18.
Nakajima EC, et al. FDA Approval summary: sotorasib for KRAS G12C-mutated metastatic NSCLC. Clin Cancer Res. 2022;28:1482–1486.
Talpaz M, Kiladjian J-J. Fedratinib, a newly approved treatment for patients with myeloproliferative neoplasm-associated myelofibrosis. Leukemia. 2021;35:1–17.
Fujita-Sato S, et al. Enhanced MET translation and signaling sustains K-Ras driven proliferation under anchorage-independent growth conditions. Cancer Res. 2015;75:2851–62.
Patricelli MP, et al. Selective inhibition of oncogenic KRAS output with small molecules targeting the inactive state. Cancer Discov. 2016;6:316–29.
Hayes TK, et al. Long-term ERK Inhibition in KRAS-mutant pancreatic cancer is associated with MYC degradation and senescence-like growth suppression. Cancer Cell. 2016;29:75–89.
Lee T, Yao G, Nevins J, You L. Sensing and integration of Erk and PI3K signals by Myc. PLoS Comput Biol. 2008;4:e1000013.
Sears R, Leone G, DeGregori J, Nevins JR. Ras enhances Myc protein stability. Mol Cell. 1999;3:169–79.
Vaseva AV, et al. KRAS Suppression-induced degradation of MYC is antagonized by a MEK5-ERK5 compensatory mechanism. Cancer Cell. 2018;34:807-822.e7.
Bilal F, et al. The Transcription Factor SLUG Uncouples Pancreatic Cancer Progression from the RAF-MEK1/2-ERK1/2 Pathway. Cancer Res. 2021;81:3849–61.
Fang Y, Vilella-Bach M, Bachmann R, Flanigan A, Chen J. Phosphatidic acid-mediated mitogenic activation of mTOR Signaling. Science. 2001;294:1942–5.
Roux PP, Ballif BA, Anjum R, Gygi SP, Blenis J. Tumor-promoting phorbol esters and activated Ras inactivate the tuberous sclerosis tumor suppressor complex via p90 ribosomal S6 kinase. Proc Natl Acad Sci U S A. 2004;101:13489–94.
Carriere A, et al. ERK1/2 phosphorylate Raptor to promote Ras-dependent activation of mTOR complex 1 (mTORC1). J Biol Chem. 2011;286:567–77.
Foster KG, et al. Regulation of mTOR complex 1 (mTORC1) by raptor Ser863 and multisite phosphorylation. J Biol Chem. 2010;285:80–94.
Talpaz M, Kiladjian J-J. Fedratinib, a newly approved treatment for patients with myeloproliferative neoplasm-associated myelofibrosis. Leukemia. 2021;35:1–17.
Qureshy Z, Johnson DE, Grandis JR. Targeting the JAK/STAT pathway in solid tumors. J Cancer Metastasis Treat. 2020;6:27.
Lokau J, Schoeder V, Haybaeck J, Garbers C. Jak-Stat Signaling Induced by Interleukin-6 Family Cytokines in Hepatocellular Carcinoma. Cancers (Basel). 2019;11:1704.
Khanna P, Lee JS, Sereemaspun A, Lee H, Baeg GH. GRAMD1B regulates cell migration in breast cancer cells through JAK/STAT and Akt signaling. Sci Rep. 2018;8:9511.
Corcoran RB, et al. STAT3 plays a critical role in KRAS-induced pancreatic tumorigenesis. Cancer Res. 2011;71:5020–9.
Mohrherr J, et al. JAK-STAT inhibition impairs K-RAS-driven lung adenocarcinoma progression. Int J Cancer. 2019;145:3376–88.
Ruan H, et al. Fedratinib attenuates bleomycin-induced pulmonary fibrosis via the JAK2/STAT3 and TGF-β1 signaling pathway. Molecules. 2021;26:4491.
Yue Y, et al. Identification of key genes involved in JAK/STAT pathway in colorectal cancer. Mol Immunol. 2020;128:287–97.
Park S-Y, et al. The JAK2/STAT3/CCND2 Axis promotes colorectal Cancer stem cell persistence and radioresistance. J Exp Clin Cancer Res. 2019;38:399.
Nagathihalli NS, et al. Inverse Correlation of STAT3 and MEK Signaling Mediates Resistance to RAS Pathway Inhibition in Pancreatic Cancer. Cancer Res. 2018;78:6235–46.
Datta J, et al. Combined Mitogen-Activated Protein Kinase/Extracellular Signal-Regulated Kinase and Signal Transducer and Activator of Transcription 3 Inhibition Uncovers Stromal Plasticity by Enriching for Cancer-Associated Fibroblasts With Mesenchymal Stem Cell-Like Features to Overcome Immunotherapy Resistance in Pancreatic Cancer. Gastroenterology. 2022;163:1593–612.
We thank members of the Solid Tumours Program for insightful comments. We also thank members of the Morphology, Genomics, and Cytometry Core Services for assistance in IHC, RNAseq and FACS work. We thank Sohail Tavazoie for developing and providing critical reagents to study PITPNC1.
R.E-C. is funded by a donation from Maria Eugenia Burgos de la Iglesia’s family. O.E. was funded by FSE; MINECO; FJCI-2017–34233. C.A. is supported by the Giovanni Armenise–Harvard Foundation, the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. ) and AIRC under IG 2021—ID. 25737 project (to C.A.). F.L. was funded by the Gobierno de Navarra (Ref. 34/2021), the Cancer Research Thematic Network of the Instituto de Salud Carlos III (RTICC RD12/0036/0066), PID2021-122638OB-I00 MCIN/AEI/10.13039/501100011033/ FEDER, UE and by FEDER “Una manera de hacer Europa”. S.V. was supported by Ministerio de Ciencia, Innovación y Universidades, Convocatoria 2019 para incentivar la Incorporación Estable de Doctores (IED2019-001007-I), by FEDER/Ministerio de Ciencia, Innovación y Universidades—Agencia Estatal de Investigación (SAF2017-89944-R and PID2020‐116344‐RB‐100), by Proyecto/MCIN/AEI /10.13039/501100011033/ FEDER, UE and European Union NextGenerationEU/ PRTR, by Fundació La Marató de TV3 (474/C/2019), and by the Scientific Foundation of Spanish Association Against Cancer—Strategic Projects 2020 (PROYE20029VICE).
Ethics approval and consent to participate
All experiments in mice were performed following ARRIVE guidelines and in concordance with the Institutional Animal Care and Use Committee (IACUC) of the University of Navarra (CEEA) under the protocol CEEA #057–18 approved by the regional Government of Navarra.
Consent for publication
F.L and S.V. receive research funding from Roche. C.A. and S.V. hold research grants by Revolution Medicines. S.V. is a consultant for LiberaBio. None of the disclosed sources are related to this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Submission declaration. The work described herein has not been published previously and is not under consideration for publication elsewhere.
Additional file 1: Suppl. Figure 1.
A. PITPNC1 expression levels (log2) in TCGA LUAD patient’s database. Mut EGFR: mutant EGFR, wt: wild type EGFR, N: Normal tissue. Mut vs wt (p=0.014), Mut vs N (p=0.021). B. PITPNC1 expression levels (log2) in TCGA LUAD patient’s database. Mut BRAF: mutant BRAF, wt: wild type BRAF, N: Normal tissue. Mut vs wt (p=0.275), Mut vs N (p=0.336). C. PITPNC1 expression levels (log2) in TCGA LUAD patient’s database with different KRAS point mutations. N: Normal tissue D. Western blot of PITPNC1 and KRAS expression in H1792 cells, expressing a control (GFPsh) or a tet-inducible KRAS shRNA (KRASsh) (activated by 1 μg/ml doxycycline). Twenty μg of protein were loaded per sample. β-TUBULIN was used as loading markes E. PITPNC1 mRNA expression in H2009 and H1792 cells expressing a control (GFPsh) or a tet-inducible KRAS shRNA (KRASsh) (activated by 1 μg/ml doxycycline) (Mann-Whitney or unpaired t-test). F. Western blot of pERK1/2, ERK1/2, pAKT, AKT, p-cJUN, cJUN, pERK5 and ERK5 in A549, H2009 and HPAFII treated with pharmacologic inhibitors: trametinib (MEKi, 0.5 μmol/L), BIX02189 (MEK5i, 10 μmol/L) or GSK2126458 (PI3Ki, 0.1 μmol/L) for 24 h, and SP600125 (JNKi, 10 μmol/L) for 2 h. Twenty μg of protein were loaded per sample. G. PITPNA, PITPNB, PITPNM1, PITPNM2, PITPNM3 mRNA expression in H2126 cells overexpressing a mutant (KRASG12D) or a wild type form (KRAS4B) of KRAS compared to the control (LacZ) (Dunnett’s multiple comparison test). H and I. PITPNA, PITPNB, PITPNM1, PITPNM2, PITPNM3 mRNA expression in H2009 (G) and H1792 (H) cells expressing a control (GFPsh) or an inducible KRAS shRNA (KRASsh) (activated by 1 μg/ml doxycycline) (Mann-Whitney or unpaired t-test). J. Western blot of PITPNC1 and KRAS expression in Kraslox/lox MEFs transduced with different human HA-tagged KRAS mutants (G12C, G12D, G12V, G12R, G12S, G13D and Q61H). K. PITPNC1 expression levels (log2) in TCGA LUAD patient’s database. Mut KRAS: mutant KRAS. Mut KRAS Amp: mutant KRAS amplification, N: Normal tissue. Mut KRAS vs Mut KRAS Amp (n.s). Suppl. Figure 2. A. PITPNC1 expression levels (log2) in TCGA LUAD patient’s database in presence of a panel of co-occurrence mutations. B. Western blot of LKB1 and PITPNC1 expression in H2009 cells expressing a control (Control) or LKB1 sgRNAs (LKB1 sgRNA1 or LKB1 sgRNA2). Twenty μg of protein were loaded per sample. β-TUBULIN was used as loading marker. C. Western blot of PITPNC1, p53, Keap1 or Lkb1 expression in KLA and LKR10 cells, expressing a control (Control) or p53, Lkb1 or Keap1 sgRNAs (p53 sgRNA, Lkb1 sgRNA or Keap sgRNA). Twenty μg of protein were loaded per sample. HSP90 and β-TUBULIN were used as loading markers. Suppl. Figure 3. A. Apoptosis analysis by Annexin V/7AAD labelling in the human LUAD A549, H2009, H1792 and H358 cell lines after PITPNC1 knockdown with a specific shRNA (sh6 or sh7) compared to the control (GFPsh). (Tukey’s multiple comparison test). B. Representative flow cytometry images depicting gating strategy in A. C. Representative images of A549- and PATU8902-derived xenografts (GFPsh, PITPNC1 sh6 and PITPNC1 sh7) stained for phospho-histone 3 (pH3) or cleaved caspase 3 (CC3). D. Western blot of PITPNC1 in A549 and H358 cells transfected with a control (pBabe) or PITPNC1 cDNA. Twenty μg of protein were loaded per sample. β-TUBULIN was used as loading marker E. Representative images of clonogenic ability of A549 and H358 PITPNC1-overexpressing cells compared with the control. F. Tumour volume (mm3) of A549 derived xenografts, (n=8), (Bonferroni’s multiple comparison test). G. Representative images of tumours from F. H. Tumor weight (g) of tumours in F (n=8) (Mann-Whitney or Unpaired T-test). I. Tumour volume (mm3) of H358 derived xenografts, (n=8), (Bonferroni’s multiple comparison test). J. Representative images of tumours from I. K. Tumour weight (g) of tumours in I (n=8) (Mann-Whitney or Unpaired T-test). L. Migration assay experiment in LacZ and PITPNC1-overexpressing A549 and H2009 cell lines. M. Representative bioluminescence images of a metastasis assay via intracardiac injection of LacZ and PITPNC1-overexpressing A549 cells. N. Ex vivo analysis of bioluminescence in lung, liver and kidney. Suppl. Figure 4. A. PHLDA2, GJB2, GPX2, BIRC5, PITPNC1, RASSF6 and ARK1B10 mRNA expression in A549 cells expressing a control (GFPsh) or a PITPNC1 shRNA (sh6 or sh7) (Dunnett´s multiple comparison test). B. Box plot comparing early (I-II) and advanced (III-IV) mut KRAS LUAD according to the expression of the dPITPNC1 gene signature (GS). C. Box plot comparing localized and advanced (locally advanced and metastatic) PDAC based on the expression of the dPITPNC1 gene signature (GS). D. Box plot comparing P53-mutated and LKB1-mutated LUAD patients with mut KRAS according to the expression of the dPITPNC1 gene signature (GS). E. Box plot comparing classical and basal PDAC patients based on the expression of the dPITPNC1 gene signature (GS). Suppl. Figure 5. A. Cell cycle analysis by EdU labelling in the human LUAD H1792 and H358 cell lines after PITPNC1 knockdown with a specific shRNA (sh6 or sh7) compared to the control (GFPsh) (Bonferroni´s multiple comparison test). B. Western blot of MYC expression in A549- and PATU8902-xenografts tumours. Twenty μg of protein were loaded per sample. HSP90 and GAPDH were used as loading markers. C. Cell cycle analysis by EdU labelling in A549 and Miapaca2 cell lines after PITPNC1 knockdown with a specific shRNA (sh6 or sh7) compared to the control (GFPsh) (Dunnet´s multiple comparison test). D. Western blot of MYC, E2F1 and P27 in A549, H2009, H1792 and Miapaca2 cells, expressing a control (GFPsh) or a shRNA against PITPNC1 (sh6 or sh7). Twenty μg of protein were loaded per sample. β-TUBULIN was used as loading control. E. mRNA analysis of RNAseq data of E2F1, P27 and P57 in A549 cells expressing a control (GFPsh) or a shRNA against PITPNC1 (sh6 or sh7). F and G. QPCR analysis of P27 (F) and P57 (G) mRNA expression in A549, H2009, H1792 and PATU8902 cells expressing a control (GFPsh) or a PITPNC1 shRNA (sh6 or sh7) (Dunnett´s multiple comparison test). H. Cell proliferation assay in A549 and H2009 cells expressing exogenous LacZ or MYC- and submitted to inhibition of PITPNC1 by specific shRNAs. I. Western blot of MYC, PITPNC1 in A549 and H2009 cells expressing exogenous LacZ or MYC- and submitted to inhibition of PITPNC1 by specific shRNAS. Twenty μg of protein were loaded per sample. HSP90 was used as loading control. J. RNAseq data of AURKA and PLK1 in A549 expressing a control (GFP) or two PITPNC1 shRNAs. Suppl. Figure 6. A. SESN1, SESN2 and SESN3 expression levels in H2009 cell line were measured by qPCR. Cells were virally infected to express a control (GFPsh) or a PITPCN1 shRNA (sh6 and sh7) (Dunnet´s multiple comparison test). GAPDH was used as housekeeping gene. B. mTOR/LAMP1 colocalization analysis by immunofluorescence in H2009 PITPNC1-depleted cells. C. Quantification of mTOR/LAMP1 Mander´s overlap coefficient (MOC) in H2009 of B (Dunnett’s multiple comparison test). D. Western blot of mTOR in A549, H2009 and H1792 cell lines expressing a control (GFP) or two PITPNC1 shRNAs. Twenty μg of protein were loaded per sample and HSP90 was used as loading control. E. Lysosomes per cell and average lysosomes size in H2009 of B (Dunn’s multiple comparison test). Lysosomes per cell and average lysosomes size in H2009 of B (Dunn’s multiple comparison test). F. Heatmap of autophagy and lysosome biogenesis genes upregulated upon PITPNC1 inhibition in A549 cells (data from RNAseq analysis). G and H. Western blot of protein level of LC3-I and LC3-II in A549, H2009 (G) and HPAFII (H) virally infected to express a shRNA control (C) or two PITPNC1 shRNAs (sh6 and sh7) and treated with or without hydroxychloroquine (CQ) (60 μM) for 6 h. Twenty μg of protein were loaded per sample and HSP90 was used as loading control. I. Western blot of mTOR signalling pathway (mTOR, 4EBP1 and S6K) in A549, H2009 and H1792 cell lines in which PITPNC1 was inhibited with specific shRNAs (sh6 and sh7). Twenty μg of protein were loaded per sample and HSP90 was used as loading control. Suppl. Figure 7. A. mTOR/LAMP1 colocalization analysis by immunofluorescence in A549 MYC-depleted cells. B. Quantification of mTOR/LAMP1 Mander´s overlap coefficient (MOC) in A549 of A (Dunnett’s multiple comparison test). C. mTOR/LAMP1 colocalization analysis by immunofluorescence in A549 MYC-depleted cells. D. Quantification of mTOR/LAMP1 Mander´s overlap coefficient (MOC) in MiaPaca2 MYC-depleted cells (unpaired t-test). E. Lysosomes per cell in A549 and MiaPaca2 of A and C respectively (Dunn’s multiple comparison test and unpaired t-test respectively). F. Average lysosomes size in A549 and MiaPaca2 of A and C respectively (Dunn’s multiple comparison test and unpaired t-test respectively). G. SESN1, SESN2 and SESN3 expression levels in A549 cell line were measured by qPCR. Cells were virally infected to express a control (GFPsh) or a MYC shRNA (sh42 and sh89) (Dunnet´s multiple comparison test). GAPDH was used as housekeeping gene. H. SESN1, SESN2 and SESN3 expression levels in Miapaca2 cell line were measured by qPCR. Cells were virally infected to express a control (GFPsh) or a MYC shRNA (sh42 and sh89) (Dunnet´s multiple comparison test). GAPDH was used as housekeeping gene. Suppl. Figure 8. A. Connectivity Map (CMap) analysis for dKRAS GS H358 transcriptomics. Perturbagen classes with mean connectivity scores >90% are displayed. Each dot represents an individual drug included in the specific class. B. Representative image of crystal violet stained plates for drug combination experiment in G12C cell lines. C. Synergistic score (Bliss score) heatmaps of H1792, H2030, H358, H23 and MiaPaca2 treated for 5 days as indicated D. Representative image of 3D proliferation assay in H358 and H1792 cell lines treated with DMSO (Ctrl) Soto (60 nM) Fedra (1 μM) or both (Combo) for 5 days. E. Representative image of crystal violet stained plates for drug combination experiment in Soto-resistant (KR) G12C cell lines H358 and H23. Suppl. Figure 9. A and B. Representative images of H358- and Miapaca2-derived xenografts stained for phospho-histone 3 (pH3) or cleaved caspase 3 (CC3). C and D. Mouse weight change upon different treatments. S= start of the experiment; E= end of the experiment (Mann-Whitney or unpaired t-test).
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Entrialgo-Cadierno, R., Cueto-Ureña, C., Welch, C. et al. The phospholipid transporter PITPNC1 links KRAS to MYC to prevent autophagy in lung and pancreatic cancer. Mol Cancer 22, 86 (2023). https://doi.org/10.1186/s12943-023-01788-w