The oncogenic axis YAP/MYC/EZH2 impairs PTEN tumor suppression activity enhancing lung tumorigenicity

doi:10.21203/rs.3.rs-5032978/v1

Download PDF

Article

The oncogenic axis YAP/MYC/EZH2 impairs PTEN tumor suppression activity enhancing lung tumorigenicity

https://doi.org/10.21203/rs.3.rs-5032978/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

The tumor suppressor PTEN (phosphatase and tensin homolog deleted in chromosome 10) is genetically deleted or downregulated in many cancer types. Loss of PTEN protein expression is frequently found in lung cancer while genetic alterations are less abundant. PTEN expression is regulated at multiple genetic and epigenetic levels and even partial reduction of its expression increases cancer occurrence.

We show that YAP and TAZ cooperate with EZH2, and MYC to transcriptionally repress onco-suppressor genes, including PTEN, in Non-small cell lung cancer (NSCLC) cells. YAP/TAZ-EZH2-MYC transcriptional regulators form a nuclear complex that represses PTEN transcription, while their combinatorial targeting restores PTEN expression, attenuates NSCLC cell growth, and prevents compensatory responses induced by single treatments. Datasets analysis of NSCLC patients revealed that PTEN expression is negatively correlated to YAP/TAZ, EZH2 and MYC and that low expression of PTEN is predictive of poor prognosis, especially at earlier stages of the disease. These findings highlight the repressive role of the YAP/TAZ-EZH2-MYC axis on tumor suppressor genes and offer a potential therapeutic strategy for lung cancer patients with low PTEN levels.

Biological sciences/Cancer/Lung cancer/Non-small-cell lung cancer

Biological sciences/Molecular biology/Epigenetics/Gene silencing

YAP

MYC

PRC2

NSCLC

tazemetostat

Lung cancer is the leading cause of cancer-related mortality in the world. It is classified into two major subtypes exhibiting different origin tissue, histology and response to therapy: small cell lung cancer (SCLC) accounts for 10–15% of all lung cancers while non-small cell lung cancer (NSCLC) for 80–85% of cases (https://cancerstatisticscenter.cancer.org). NSCLC is further categorized into lung adenocarcinoma (LUAD), the most common subtype arising from secretory cells of distal airways, and lung squamous cell carcinoma (LUSC) that, together with large cell carcinoma (LCC), originates from more proximal airways (1–4). The high mortality rate of lung cancer, despite its relatively low incidence rate, is frequently associated with late detection at advanced stages. Nevertheless, in recent years, the death rate is progressively declining due to powerful early detection, reduced smoking, and improved treatments using target- and immuno-therapy in addition to traditional chemo- and radio-therapy (https://cancerstatisticscenter.cancer.org) (4–6). Yet, not all patients are eligible for target- and immuno-therapy, and even the eligible ones, frequently develop resistance over time. It is, therefore, important to identify potent targets and to define the underlying mechanisms of resistance to commonly used therapies.

PTEN is a key onco-suppressor, which negatively regulates the PI3K/AKT/mTOR oncogenic pathway, and is involved in DNA repair. PTEN downregulation is an early event in lung tumorigenesis and it correlates with poor overall survival, progression-free survival and disease-free survival (7). Molecularly, the reduction of PTEN expression induces the aberrant activation of the PI3K/AKT/mTOR pathway which in turn controls a plethora of oncogenic events, including cell growth, stemness, proliferation and survival, protein synthesis, energy metabolism, motility and cellular architecture (8–12). Moreover, PTEN loss or downregulation confers resistance to radiotherapy, chemotherapy, target therapy and immunotherapy in some settings (13–15). According to the Cancer Genome Atlas (TCGA), PTEN genetic mutations in LUSC are found in approximately 15% of cases, but only in ~ 2–3% of the LUAD patients. The low frequency of PTEN genetic alteration in LUAD cannot explain the higher frequency of downregulation of PTEN protein, which was shown to be dependent on multiple genetic and epigenetic mechanisms as well as protein-protein interactions (8). In the present work, we show that YAP and TAZ cooperate with EZH2 and MYC to transcriptionally repress PTEN in LUAD cells. These transcriptional regulators maintain low levels of PTEN.

YAP and TAZ are widely known as transcriptional coactivators of different oncogenes, mainly those who are related to cell cycle, stemness and epithelial-mesenchymal transition (EMT). Genes commonly activated by those proteins in different cancer types are known as the oncogenic YAP- and TAZ-mediated gene signature (16). In this work, we focus on the role of YAP and TAZ as transcriptional repressors of onco-suppressor genes. We observed that pharmacological inhibition of EZH2 or MYC alone inefficiently upregulates PTEN expression, possibly through activation of oncogenic pathways including the YAP/TAZ signature. However, the concomitant inhibition of YAP/TAZ together with EZH2 and MYC robustly derepressed the expression of PTEN and other onco-suppressor genes and suppressed the hyperactivation of the YAP and TAZ signature. Our study identified a novel oncogenic axis that can be targeted in lung cancer and highlights the suppression activity of YAP/TAZ on PTEN transcription.

YAP, PRC2 and MYC contribute to the transcriptional repression of tumor-suppressor genes in NSCLC

Previous studies showed that YAP and EZH2 synergize to repress the transcription of lineage-specific genes and onco-suppressor genes in NSCLC, Schwann cells, trophoblast and intestine, and that YY1 (Yin Yang 1) transcription factor may bridge YAP and EZH2 at the promoter of cell-cycle regulators (17–20). MYC was also proposed to cooperate with EZH2 or YY1 to repress the transcription of onco-suppressor genes in cancer (21–23) while cooperating with YAP/TAZ mainly to enhance the transcription of oncogenes (24)].

Considering the oncogenic impact of PTEN loss on NSCLC tumorgenicity and its possible downregulation by transcriptional and post-transcriptional mechanisms, we examined its promoter region using the Cistrome Data Browser and found a clear co-localization of TEAD4, YY1 and MYC onto its promoter in the A549 NSCLC cell line (Fig. 1A). TEAD1-4 transcription factors are known to bind YAP/TAZ on both activated and repressed genes promoters. Therefore, despite the absence of YAP track in the Cistrome DB of A549 cells, we assume that YAP binds PTEN promoter since a YAP peak overlapping with TEAD4, YY1 and MYC is present in the mesothelioma cell line H2052 (Fig. 1A). The presence of YAP, MYC and YY1 on PTEN promoter was validated by chromatin immunoprecipitation (ChIP) in the H1299 NSCLC cell line (Fig. 1B). Furthermore, co-immunoprecipitation (Co-IP) and immunofluorescence (IF) experiments in H1299 cells revealed physical interaction and cellular co-localization of endogenous YAP with TEAD1, pMYC, EZH2, YY1 as well as with LaminA/C, a structural protein associated with the inner nuclear membrane where the epigenetically repressed chromatin is clustered (25) (26) (Fig. 1C-J, S8). Interestingly, pMYC also co-immuno-precipitates, co-localizes and interacts in the nucleus with Lamin A/C (Fig. 1C, F, S8). We obtained similar patterns of YAP distribution and co-localization with Lamin A/C, pMYC, YY1 and H3K27me3 in H1975 cells (Fig. S1). Notably, H3K27me3 is enriched at the nuclear lamina (Fig. S1D) with a ring-like pattern, similar to YAP (Fig. S1A).

Proximity ligation assay (PLA) in H1299 confirmed in situ the interaction of YAP with EZH2, YY1, pMYC, and lamin A/C as well as the interaction between pMYC and lamin A/C (Fig. 1K).

By matching published lists of genes upregulated in YAP- or TAZ-depleted cells (GSE151200) (27), which are possibly target genes transcriptionally repressed by YAP or TAZ, with those that are bound by MYC (ChIP-Atlas) and either EZH2 (28) or H3K27me3 (GSE29611) in A549 cells, we obtained two distinct lists of genes (Table S1). These lists were obtained by applying pval < 0.05 and padj < 0.3. The heatmaps and Venn diagram shown in Fig.S2A, B show that TAZ/EZH2/MYC co-represses a higher number of genes (N = 165) compared to YAP/EZH2/MYC (N = 50 genes). Moreover, 20 genes are commonly repressed by both YAP and TAZ, indicating that almost half of YAP repressed genes (40%) are also repressed by TAZ. However, only 12% of genes repressed by TAZ are also repressed by YAP (Fig.S2A, B, Table S1). Genes set enrichment analysis (GSEA) of genes repressed by either YAP or TAZ and co-bound by MYC and EZH2/H3k27me3 highlights apical junction, apical surface, EMT, genes downregulated by KRAS activation and IL2-STAT5 signaling, as the most significant pathways affected by these proteins (Fig. S2C, Table S1, sheet 3). Collectively, these data suggest that the transcriptional repression mediated by YAP, TAZ, MYC and EZH2 may be a general mechanism, and supports the emerging notion that YAP and TAZ may have both overlapping and distinct functions not only in the transcriptional activation but also in the transcriptional repression of different subsets of genes (29).

Oncogenic MYC cooperates with YAP and TAZ to repress transcription of PTEN

We have noticed that PTEN was not present in the heatmaps of Fig. S2A and Table S1, possibly because the analysis included upregulated genes in response to single KD of either YAP or TAZ. Indeed, when YAP and TAZ were concomitantly depleted in two NSCLC cell lines, H1299 and H1975, a strong upregulation of PTEN was observed, concurrent with a decrease in CYCD1 (cyclin D1) transcript, a cell cycle regulator that is negatively modulated by PTEN (Fig. 2A, B, S9A) (30). Importantly, the expression of CYCD1 and PTEN are anti-correlated in NSCLC patients (Fig. 2C). KD of PTEN robustly upregulated CYCD1 and rescued CYCD1 levels in siYAP/TAZ double KD cells (Fig. 2D, S9B). Conversely, PTEN overexpression reduced CYCD1 levels and increased the fraction of cells at the G1/S transition in two other NSCLC cell lines, H1975 and H358, with no evident effects in H1299 cells (Fig. S2D, E, S12B) (30).

To evaluate the impact of MYC on PTEN expression, we downregulated its expression by siRNA in H1299 and H1975 cells, and observed a significant upregulation of PTEN only in H1975 cells (Fig. 2E, F, S9C). However, colony formation was inhibited in both cell lines upon MYC depletion (Fig. 2G), suggesting that MYC elicits its oncogenic function through several mechanisms beyond PTEN downregulation. Upon MYC depletion, we also observed a cell line-specific de-repression of TGFBR2 and p21, previously shown to be repressed by YAP, TAZ and EZH2, implying a similar repression mechanism for those genes (Fig. 2E, F, S9C) (18). Although depletion of MYC alone did not affect PTEN expression in H1299 cells, we found that MYC overexpression in YAP/TAZ depleted H1299 cells could rescue PTEN, p21 and TGFBR2 repression, as well as a partial rescue of the colony-forming ability (Fig. 2H, I, S9D) implying a role for MYC in PTEN repression also in H1299 cells. In LUAD patients, upregulation of EZH2, MYC and YY1 was associated with downregulation of PTEN expression, both at the protein and the RNA level. This association was specifically observed in patients harboring somatic alterations in the Hippo pathway and MYC-related pathways (Fig. S3A, B). Moreover, the expression level of PTEN protein negatively correlates with its putative transcriptional repressors TAZ, EZH2 and MYC while it positively correlates with pYAPS127 (the YAP phosphorylated form that is retained in the cytoplasm), as well as with the onco-suppressive transcripts TGFBR2, CDKN1A (p21) and CDKN2B (p15), which were previously shown to be repressed by YAP/TAZ and EZH2 (Fig. S3C) (18). Finally, Kaplan Mayer survival analysis shows a better prognosis for LUAD patients expressing high levels of PTEN, especially in the earlier stages of the disease, supporting its role as an early prognostic onco-suppressor (Fig. S3D).

EZH2 and MYC inhibition activates YAP and TAZ targets

The oncogenic role of EZH2 and MYC in the transcriptional repression of onco-suppressors suggests that their inhibition would reactivate those onco-suppressors. Several specific inhibitors of EZH2, which catalyzes tri-methylation of histone H3 at Lys 27 (H3K27me3), have been developed, and Tazemetostat has been approved by the FDA for epithelioid sarcoma (31), diffuse large B-cell lymphoma and relapsed or refractory follicular lymphoma (32). While Tazemetostat is effective for the above cancer types, it is ineffective as monotherapy for other malignancies, especially for solid cancers, where cells show either innate or acquired resistance to this compound (33). Huang and coworkers showed that Tazemetostat treatment of solid tumors induced a decrease in global H3K27me3, as expected, which was accompanied by an increase in total H3K27ac. It was, therefore, proposed that an oncogenic transcriptional reprogramming mediated by the interaction of the epigenetic protein MLL1 with the p300/CBP complex onto several oncogenes induces global H3K27ac upregulation upon H3K27me loss, driving the hyper activation of genes belonging to oncogenic pathways that counteract the de-repression of onco-suppressor pathways (33).

We previously showed that H1299 and H1975 cells are resistant to Tazemetostat, with an IC50 of 24 µM and 36 µM, respectively (18). We also observed a reduced number of colonies in H1299 but not in H1975 or h358 cells at micromolar doses of Tazemetostat (Fig. S4A). In these cells, depletion of EZH2 by siRNA or treatment with 2µM Tazemetostat increased the expression of YAP and TAZ and their transcriptional targets CTGF, ANKRD1 and CYR61 (Fig. 3A-C, S10A). These effects were also observed in H1299 and H1975 cells upon treatment with a lower dosage of Tazemetostat (1µM, Fig. S4B, S12C, D). Moreover, Tazemetostat treatment induces a dose-dependent increase of global H3K27ac that counteracts the H3K27me3 decrease (Fig. 3D, S10B, C). This epigenetic switch was observed previously in other Tazemetostat-resistant cells but not in sensitive cell lines (33). We found that global levels of YAP and TAZ are increased upon Tazemetostat treatment in a dose-dependent fashion, especially at high doses and more strongly in H1975 cells (Fig. 3D, S10B, C) which are more resistant to Tazemetostat than H1299 (Fig. S4A) (18). In accordance, both in lung adenocarcinoma (LUAD) and in lung squamous cell carcinoma (LUSC) patients, the EZH2 mRNA has a negative correlation with CTGF, ANKRD1 and CYR61 transcripts (Fig. 3E).

Based on these results, we speculated that the upregulation of YAP- and TAZ- oncogenic signature may associate with the acquired resistance of H1299 and more strongly of H1975 cells to Tazemetostat treatment. Strikingly, the knockdown of YAP and TAZ increased the sensitivity of H1975 cells to very low doses of Tazemetostat (Fig. S4C) compared to Tazemetostat alone, which did not affect colony formation even at higher doses (up to 8µM, Fig. S4A). Interestingly, MYC depletion or inhibition by MYCi975 also increased the expression of YAP and TAZ targets in H1299, H1975 and H358 cells (Fig. 3F-G, S4D-E, S10D-E, S12E, S13A) consistent with previous reports in liver and breast cancer models (24, 34). Furthermore, the three cell lines were resistant to micromolar doses of MYCi975 as shown by the colony formation assay (Fig. S4F). In contrast, treatment with JQ1, a BRD4 inhibitor that indirectly inhibits oncogenic MYC function (35), induced opposite effects in different cell lines, increasing the expression of total YAP, TAZ, CTGF and ANKRD1 in H1299 but decreasing them in H1975 cells (Fig. S4G, S13B). This might be related to the broad effects of JQ1, which not only inhibits MYC but also other oncogenic pathways through epigenetic mechanisms. Notably, it was previously shown that combination of Tazemetostat with JQ1 inhibits MLL1-mediated H3K27ac accumulation onto oncogenes, and reverses the resistance to EZH2 inhibitors in different cellular models (33). However, considering the different effects of JQ1 treatment in our models, we decided to inhibit MYC directly with MYCi975 and hypothesized that combined inhibition of YAP/TAZ together with MYC and EZH2 may be the best strategy to prevent activation of YAP/TAZ-mediated oncogenic pathways and concurrently de-repressing onco-suppressive pathways.

Concomitant inhibition of EZH2, MYC and YAP/TAZ maximizes the de-repression of onco-suppressors while preventing the activation of oncogenes

We next depleted EZH2, MYC and YAP/TAZ by RNA interference (RNAi) either alone or in combination. We found that the triple knock-down (KD) of EZH2, MYC and YAP/TAZ led to stronger de-repression of onco-suppressor genes (TGFBR2, PTEN, p21) while reducing the activation of YAP and TAZ targets, which was observed upon double EZH2/MYC KD (Fig. 4A, S11A). Consistent with these results, we found that triple KD increased the binding of active RNAPol2 onto the PTEN promoter but decreased the binding onto the CTGF and CYR61 promoters (Fig. 4B). Functionally, the triple KD had stronger anti-tumorigenic effects compared to single or double KD in NSCLC cell lines, as shown by the reduced colony formation potential (Fig. 4C) and by the apoptotic phenotype demonstrated by AnnexinV/PI staining and by cleaved casp3 (Fig. 4AD-E, S5A, S11B). Conversely, KD of PTEN reduced the percentage of apoptotic cells already 24h post-transfection (Fig. S5B) and the level of cleaved PARP/casp3 (Fig. S5C, S14). These results suggest that the increased apoptosis observed in response to triple KD is possibly mediated, at least in part, by the increase in PTEN and the concurrent decrease of YAP and TAZ. Moreover, PTEN KD increased colony formation in H1299, H1975 and H358 cells (Fig. S5D). To further confirm the KD results, we treated H1299 and H1975 cells with low doses of MYCi975 (36), Tazemetostat and Dasatinib (an indirect YAP and TAZ inhibitor) (37). Low doses of the combined three inhibitors reduced the viability of both cell lines and de-repressed onco-suppressor genes while suppressing reactivation of YAP and TAZ oncogenes (Fig. S6A, B, S15) consistent with the KD experiments.

Knockdown of PTEN restores the tumorigenic properties of MYC, EZH2 and YAP/TAZ depleted NSCLC cells.

To further demonstrate the link between the three transcriptional regulators: MYC, EZH2 and YAP/TAZ and the tumor suppressor PTEN, we downloaded a list of transcripts and proteins positively or negatively correlated to PTEN in LUAD patients from the lung adenocarcinoma TCGA, Firehose Legacy study (Table S2). Gene set enrichment analysis (GSEA) of anti-correlated proteins highlighted E2F targets, mTORC1 signaling, and MYC targets as significantly anti-correlated to PTEN, consistent with the roles of PTEN in cell proliferation, cell growth, and cell metabolism (Table S2, Fig. 5A).

To demonstrate that inhibition of PTEN expression through the concerted function of MYC, EZH2, and YAP/TAZ contributes to NSCLC tumorigenicity, we knocked down PTEN together with MYC, EZH2, and YAP/TAZ (defined as “MEY”). MEY KD induced PTEN and p21 reactivation and concurrently inhibited the PI3K/AKT/mTOR pathway as shown by reduced pAKT levels (Fig. 5B, S12A). These effects were accompanied by reduced mitochondrial respiration (Fig. 5C), which was previously shown to be essential for early events of lung tumorigenesis (38). Strikingly, the concomitant KD of PTEN could partially rescue these cellular changes (Fig. 5B, C, S7, S12A).

Altogether, these observations are summarized in Fig. 5D: in cancer context, onco-suppressors are repressed while oncogenes are hyper-activated due to the aberrant functions of EZH2, MYC and YAP/TAZ. Treatment with Tazemetostat (or with MYC inhibitors) as a single agent, reactivates onco-suppressors on one hand but concomitantly may induce stronger activation of oncogenic signaling, including the activation of YAP/TAZ signature, as a mechanism of acquired resistance. Targeting of YAP/TAZ oncogenic signature may overcome the acquired resistance to Tazemetostat and/or MYC inhibitors. Combined targeting of EZH2, MYC and YAP/TAZ can, therefore, effectively de-repress expression of tumor suppressors and concurrently repress the expression of oncogenic pathways, to ultimately inhibit tumor growth.

In this study, we highlight the role of YAP/TAZ not only as transcriptional co-activators of oncogenes but also as transcriptional co-repressors of onco-suppressor genes. Previous studies demonstrated that YAP can function as a transcriptional repressor in breast (39), Schwann cells, trophoblast and intestine, through different mechanisms (17–20). In breast, YAP was found at the promoter of a subset of genes, including those of the 1q21.3 locus. In that context, a complex of YAP, TEAD4, AIB1 and ANCO1 functions as a transcriptional repressor in normal breast cells, while during tumorigenesis, the loss of ANCO1 leads to a switch from transcriptional repression toward transcriptional activation (39). Furthermore, the Nucleosome Remodeling Deacetylase (NuRD) complex mediates the transcriptional repressive function of YAP and TAZ on several target genes, including DDIT4, a negative regulator of mTORC1 (40). Our work provided a novel mechanism mediated by YAP and TAZ to sustain the mTOR pathway, specifically through PTEN downregulation directly at the transcriptional level, in addition to previous works highlighting DDIT4 repression (40) or the upregulation of miR29 that leads to post-transcriptional repression of PTEN (41). Our newly characterized repression mechanism introduces new players, the EZH2 methyltransferase and the MYC transcription factor. EZH2 was found to co-repress onco-suppressor genes with YAP also in Schwann cells (19) and in breast cancer cells (42) as well as inflammation-related genes in intestinal cells (20). In the study by Wang and co-worker, YAP was shown to recruit EZH2 onto GDF15 promoter, resulting in its repression and metastatization of breast cancer cells (42), suggesting that this mechanism of transcriptional repression might be conserved in different cell types. The discovery of the EZH2/MYC/YAP/TAZ oncogenic axis provides the rationale for the inhibition of these proteins, and we found in vitro that the best functional effect can be obtained through the simultaneous inhibition of all these players.

Cell lines used for in vitro experiments

NCI-H1299 (RRID:CVCL0060), NCI-H1975 (RRID:CVCL_1511), NCI-H358 (RRID:CVCL_1559) were purchased from the American Type Culture Collection (ATCC, Manassas, VA) and routinely tested by PCR for mycoplasma contamination by using the following primers: Myco_fw1: 5′-ACACCATGGGAGCTGGTAAT-3′, Myco_rev1: 5′-CTTCATCGACTTTCAGACCCAAGGCA-3′. Cells were grown in RPMI medium (Invitrogen, Carlsbad, CA) supplemented with 10% fetal bovine serum and Pen/Strep antibiotic at 37 °C in a balanced air humidified incubator with 5% CO2. Lipofectamine RNAimax (Invitrogen) was used in accordance with the manufacturer's instruction for transfection with siRNAs. SiRNAs were used at the final amount of 300 pmol in 60 mm dish. List of siRNAs used for functional in vitro experiments is given below. Cells were collected 48–72 h post transfection for subsequent analyses. All cell lines have been authenticated in the past three years.

Transient transfection with plasmids

The plasmid for MYC overexpression was obtained by cloning the human MYC cDNA into pCDNA3 backbone and was a kind gift from the laboratory of Dr. Maurizio Fanciulli (IRCCS Regina Elena, Rome Italy). pCDNA empty vector was used as control. Cells in suspension were transfected with siGFP or siYAP/TAZ using Lipofectamine RNAimax. 24h post-transfection, cells in adhesion were transiently transfected with pCDNA3 empty vector of with pCDNA3-MYC using Lipofectamine 2000 according to manufacturer instructions. Cells were collected 48h after plasmid transfection for subsequent analyses. Plasmid for PTEN overexpression was obtained by cloning the PTEN wt cDNA into pEGFP-C2 backbone and was a kind gift of Dr. Fabiana Conciatori (IRCCS Regina Elena, Rome Italy). Empty vector was used as negative control.

Clonogenic assay

Cells were transfected as indicated above and 48 h–72 h later they were detached and seeded at 500–1000 cells/well into 6-well or 12-well dishes. Fresh medium was added every 4 days. After 7–14 days, colonies were stained with crystal violet and counted.

For colony assay upon treatment, cells were seeded in triplicate at a density of 1000 cells/well in a six-well multiwell with or without Tazemetostat or MYCi975 at the indicated doses. Every 3-4 days, fresh medium was added with or without Tazemetostat or MYCi975, then after 10-14 days the colonies were stained with crystal violet and counted. Two-tailed t-test analysis was applied to calculate the P values

Analysis of cell cycle profile

For cell cycle analysis, cells were collected 48-72h after PTEN overexpression and fixed in 70% ethanol overnight. Fixed cells were washed and treated with RNase at 1 mg/ml final concentration for 30 min at 37°C before adding 5 mg/ml PI and analyzed with the Attune flow cytometer (Invitrogen) and FlowJo software.

Annexin V/PI assay

Analysis of apoptotic cells upon different interference and treatment conditions was performed by using the Annexin V-FITC Apoptosis Kit 300 (Invitrogen #BMS500FI-300) following manufacturer instructions. Stained cells were analyzed with a Guava Easycyte 8HT flow cytometer equipped with Guava Soft 2.1 (Millipore) according to the manufacturer’s instructions.

Pharmacological treatment and chemical reagents

Dasatinib (#S1021), Tazemetostat (EPZ-6438 #S7128), JQ1 (#S7110) and MYCi975 (#S8906) were obtained from Selleck Chemicals (Houston, TX).

Protein extracts and western blot analysis

For the preparation of whole-cell lysates, cells were lysed in buffer with 50 mM Tris–HCl pH 7.6, 0.15 M NaCl, 5 mM EDTA, 1% Triton X-100 and fresh protease inhibitors. Extracts were sonicated for 10 + 15 s at 80% amplitude and centrifuged at 12 000~ rpm for 10 min to remove cell debris. Protein concentrations were determined by colorimetric assay (Bio-Rad). Antibodies used for Western Blotting are listed in Supplementary Materials. Primary antibodies were diluted 1:1000, while secondary antibodies were diluted 1:10000 in Tris-Buffered Saline (TBS) with 5% BSA and 0.1% Tween® 20 detergent.

Transcript analysis

Total RNA was extracted using TRIzol (Ambion) according to the manufacturer's recommendations. For gene transcript analysis, 1γ RNA was retrotranscribed using M-MLV reverse transcriptase (Invitrogen) following the manufacturer instructions. Real time PCR was performed into a final volume of 10 μl using Sybr Green PCR master mix, and normalized on GAPDH. All the real-time PCR assays were performed by using an Applied Biosystems® 7500 fast or a StepOne Real-Time PCR Instrument. Each analysis was performed at least on three independent biological replicates and two-tailed t-test analysis was applied to calculate the P values. List of primers used for Real time PCR is given in Supplementary Materials.

Chromatin Immunoprecipitation (ChIp)

ChIP-qPCR was performed as described previously (43). Briefly, cells were fixed with 1% formaldehyde (Sigma) in culture medium for 10 min at room temperature, and chromatin from lysed nuclei was sheared to 200–600 bp fragments using a Bioruptor sonicator (Diagenode). 100 μg of sheared chromatin and 5 μg of antibody plus 40 μl of magnetic beads (Dynabeads® Protein G 10004D ThermoFisher) were used for each Ip. Quantitative real-time PCR was carried out with Applied Biosystems® 7500 fast or StepOne Real-Time PCR Instruments. Each sample was analyzed in triplicate. The amount of immunoprecipitated DNA in each sample was determined as the fraction of the input (or mock) [amplification efficiency (Ct INPUT–Ct ChIP)], and normalized to the IgG control.

Co-immunoprecipitation (CoIP)

Co-Ip was performed as described previously (44).

Mitochondrial respiration analysis

For mitochondrial respiration analysis, cells were collected 48 h after interference and seeded in a Seahorse XFp PDL Miniplate (Agilent) to grow overnight in fresh medium.The next day, after a 30-min calibration of the XF sensor with the preincubated sensor cartridge, the cell plate was loaded into the analyzer, and Oxygen Consumption Rate (OCR) was analyzed under basal conditions. Net mitochondrial oxygen consumption was then analyzed by sequential injection of the complex inhibitors oligomycin (1 µM), FCCP (1.2 µM), and the mixture of rotenone and antimycin A (0.5 μM each). Cells were stained with DAPI and the cell count of each well was determined by imaging the cells using Cell Profiler. Three replicates of each sample were analyzed Two-tailed t-test analysis was applied to calculate the P values.

Immunofluorescence

Cells were seeded overnight onto cover glasses in multi-well at different confluence.

The next day cells were washed twice in Tris-Buffered Saline (TBS), fixed in 4% paraformaldehyde (PFA) in TBS for 10 min at RT and washed three times in TBS for 5 min at RT. Then cells were permeabilized with 0.1% triton-X 100 in TBS for 10 min at RT, washed three times in TBS for 5 min at RT, blocked in 3% BSA in TBS-T (0.1% tween) for 45’ at RT. After blocking, cells were incubated with primary antibody overnight in humidity chamber at 4°C. Primary antibody dilution ranges from 1:25 to 1:500 depending on the antibody. The next day, cells were washed three times in TBS-T for 5’ at RT and incubated with secondary antibody (1:500 dilution) for 1h at RT in the dark, then washed three times in TBS-T 5’ at RT in the dark and incubated with Hoechst 1:2000 in TBS for 10’ at RT in the dark. Finally, cover glasses with cells were washed three times in TBS-T for 5’ at RT in the dark, mounted with 10 microliters of 70% glycerol onto microscope slides and analyzed when dried. PLA assay was performed using a Duolink in situ PLA kit (Sigma-Aldrich) and following manufacturing instructions.

Confocal Microscopy and Image Analysis.

Fixed and stained cells were imaged using a Zeiss 710 laser-scanning microscope under 63× Plan Apochromat (1.4 numeric aperture) oil immersion lens with Zen 2008 imaging software (Zeiss). Stacks were collected in 0.5-μm Z-section increments, and a representative stack is shown for each image.

Opera phenix plus platform

For immunofluorescence experiments in H1975 cells, 500-100 cells/well were plated in 96-well phenoplates (Revvity). Images were acquired using the Opera Phenix Plus (Revvity) platform with a 20x water objective. The following excitation laser were used: 375 nm for the identification of nuclei stained with Hoechst 33342 (Sigma), 488nm for the identification of YAP (rabbit), pMYC, YY1 and H3K27me3, incubated with anti-rabbit secondary antibody conjugated with Alexa Fluor 488 (Invitrogen), 640nm for the identification of YAP (mouse) and LaminA/C incubated with anti-mouse secondary antibody conjugated with Alexa Fluor 647 (Invitrogen). The Intensity of fluorescence and their co-localization were analyzed using Harmony software (Revvity).

Bioinformatic analysis

Targets of MYC in A549 cell line were identified through the online database ChiP-Atlas, (https://chip-atlas.org/). A hierarchical clustering analysis was performed using the targets modulated after silencing of YAP and TAZ. The RNAseq data used for this analysis was obtained from the Gene Expression Omnibus (GEO) with the accession number GSE151200. Genes obtained from GSE151200 were matched with MYC targets, (ChiP-Atlas), and also with genes bound by either EZH2 (28) or H3K27me3 (GSE29611) in A549 NSCLC cells. To gain insights into the biological pathways associated an enrichment pathway analysis was conducted. ShinyGO 0.77, a web-based tool (http://bioinformatics.sdstate.edu/go/), was utilized for this purpose.

Curves of OS

Curves of overall survival in lung cancer patients were evaluated by Kaplan–Meier method by the Kaplan Mayer plotter tool (kmplot.com). Curves of patients with high and low signals were established by median value. Statistical significance between curves was assessed by using the logrank test.

Analysis of differentially expressed genes

Transcripts or proteins differentially expressed in tumoral vs non-tumoral tissue in NSCLC patients were obtained from the UALCAN database (https://ualcan.path.uab.edu/analysis.html). Pathway alteration was measured through systematic evaluation of pathway-level somatic alterations by small mutation or copy number alteration across tumors with combined proteomic, whole-exome, and CNA data, involving key pathways and genes previously annotated across multiple cancer types based on domain knowledge (pathway annotations from PMID: 33568653).

Correlation analyses

Correlation analyses between proteins and transcripts were performed through the cBioportal Database (https://www.cbioportal.org).

Data availability

The gene expression and Chip-Seq data analyzed during the current study are available in the Gene Expression Omnibus (GEO) repository under accession number GSE151200 (genes upregulated upon YAP or TAZ interference in A549 cells), GSE29611 (genes enriched of H3K27me3 in A549 cells) and Chip-Atlas (MYC targets in A549 cells: https://chip-atlas.org/).

Acknowledgements

We thank dott. Fabiana Conciatori for kindly providing PTEN overexpression plasmid.

GB work is supported by the MUR-PNRR M4C2I1.3 PE6 project PE00000019 Heal italia (CUP H83C22000550006) and PNC D3-4H Health Digital Driven Diagnostics, prognostics and therapeutics for sustainable Health care.

Fundings

GB received fundings from the MUR-PNRR M4C2I1.3 PE6 project PE00000019 Heal Italia (CUP H83C22000550006) and PNC D3-4H Health Digital Driven Diagnostics, prognostics and therapeutics for sustainable Health care.

Author contribution

FLS designed the research study and wrote the manuscript. FLS and CT performed the majority of the experiments. BM performed some of the WB and colony assays. AS performed all the bioinformatic analyses. CP performed the Opera Phoenix analysis. SS discussed some research results. FLS, GB and SL interpreted and discussed the research results.

All authors read and approved the final manuscript.

Correspondence should be addressed to [email protected]

Conflict of interest

All the authors declare no conflict of interest.

Relli V, Trerotola M, Guerra E, Alberti S. Abandoning the Notion of Non-Small Cell Lung Cancer. Trends Mol Med. 2019;25(7):585-94.
Devarakonda S, Morgensztern D, Govindan R. Genomic alterations in lung adenocarcinoma. Lancet Oncol. 2015;16(7):e342-51.
Sainz de Aja J, Dost AFM, Kim CF. Alveolar progenitor cells and the origin of lung cancer. J Intern Med. 2021;289(5):629-35.
Leiter A, Veluswamy RR, Wisnivesky JP. The global burden of lung cancer: current status and future trends. Nat Rev Clin Oncol. 2023;20(9):624-39.
Meng G, Liu X, Ma T, Lv D, Sun G. Predictive value of tumor mutational burden for immunotherapy in non-small cell lung cancer: A systematic review and meta-analysis. PLoS One. 2022;17(2):e0263629.
Mithoowani H, Febbraro M. Non-Small-Cell Lung Cancer in 2022: A Review for General Practitioners in Oncology. Curr Oncol. 2022;29(3):1828-39.
Xiao J, Hu CP, He BX, Chen X, Lu XX, Xie MX, et al. PTEN expression is a prognostic marker for patients with non-small cell lung cancer: a systematic review and meta-analysis of the literature. Oncotarget. 2016;7(36):57832-40.
Gkountakos A, Sartori G, Falcone I, Piro G, Ciuffreda L, Carbone C, et al. PTEN in Lung Cancer: Dealing with the Problem, Building on New Knowledge and Turning the Game Around. Cancers (Basel). 2019;11(8).
Carnero A, Blanco-Aparicio C, Renner O, Link W, Leal JF. The PTEN/PI3K/AKT signalling pathway in cancer, therapeutic implications. Curr Cancer Drug Targets. 2008;8(3):187-98.
Worby CA, Dixon JE. PTEN. Annu Rev Biochem. 2014;83:641-69.
Goo CK, Lim HY, Ho QS, Too HP, Clement MV, Wong KP. PTEN/Akt signaling controls mitochondrial respiratory capacity through 4E-BP1. PLoS One. 2012;7(9):e45806.
Sanaei MJ, Razi S, Pourbagheri-Sigaroodi A, Bashash D. The PI3K/Akt/mTOR pathway in lung cancer; oncogenic alterations, therapeutic opportunities, challenges, and a glance at the application of nanoparticles. Transl Oncol. 2022;18:101364.
Exposito F, Redrado M, Houry M, Hastings K, Molero-Abraham M, Lozano T, et al. PTEN Loss Confers Resistance to Anti-PD-1 Therapy in Non-Small Cell Lung Cancer by Increasing Tumor Infiltration of Regulatory T Cells. Cancer Res. 2023;83(15):2513-26.
Fischer T, Hartmann O, Reissland M, Prieto-Garcia C, Klann K, Pahor N, et al. PTEN mutant non-small cell lung cancer require ATM to suppress pro-apoptotic signalling and evade radiotherapy. Cell Biosci. 2022;12(1):50.
Ferrara MG, Martini M, D'Argento E, Forcella C, Vita E, Di Noia V, et al. PTEN Loss as a Predictor of Tumor Heterogeneity and Poor Prognosis in Patients With EGFR-mutant Advanced Non-small-cell Lung Cancer Receiving Tyrosine Kinase Inhibitors. Clin Lung Cancer. 2021;22(4):351-60.
Cordenonsi M, Zanconato F, Azzolin L, Forcato M, Rosato A, Frasson C, et al. The Hippo transducer TAZ confers cancer stem cell-related traits on breast cancer cells. Cell. 2011;147(4):759-72.
Meinhardt G, Haider S, Kunihs V, Saleh L, Pollheimer J, Fiala C, et al. Pivotal role of the transcriptional co-activator YAP in trophoblast stemness of the developing human placenta. Proc Natl Acad Sci U S A. 2020;117(24):13562-70.
Lo Sardo F, Pulito C, Sacconi A, Korita E, Sudol M, Strano S, et al. YAP/TAZ and EZH2 synergize to impair tumor suppressor activity of TGFBR2 in non-small cell lung cancer. Cancer Lett. 2021;500:51-63.
Hoxha S, Shepard A, Troutman S, Diao H, Doherty JR, Janiszewska M, et al. YAP-Mediated Recruitment of YY1 and EZH2 Represses Transcription of Key Cell-Cycle Regulators. Cancer Res. 2020;80(12):2512-22.
Zhu H, Lu J, Fu M, Chen P, Yu Y, Chen M, et al. YAP represses intestinal inflammation through epigenetic silencing of JMJD3. Clin Epigenetics. 2024;16(1):14.
Barton K, Margolis D. Selective targeting of the repressive transcription factors YY1 and cMyc to disrupt quiescent human immunodeficiency viruses. AIDS Res Hum Retroviruses. 2013;29(2):289-98.
de Nigris F, Botti C, Rossiello R, Crimi E, Sica V, Napoli C. Cooperation between Myc and YY1 provides novel silencing transcriptional targets of alpha3beta1-integrin in tumour cells. Oncogene. 2007;26(3):382-94.
Kuser-Abali G, Alptekin A, Cinar B. Overexpression of MYC and EZH2 cooperates to epigenetically silence MST1 expression. Epigenetics. 2014;9(4):634-43.
Croci O, De Fazio S, Biagioni F, Donato E, Caganova M, Curti L, et al. Transcriptional integration of mitogenic and mechanical signals by Myc and YAP. Genes Dev. 2017;31(20):2017-22.
Cesarini E, Mozzetta C, Marullo F, Gregoretti F, Gargiulo A, Columbaro M, et al. Lamin A/C sustains PcG protein architecture, maintaining transcriptional repression at target genes. J Cell Biol. 2015;211(3):533-51.
Mustafin RN, Khusnutdinova EK. The relationship of lamins with epigenetic factors during aging. Vavilovskii Zhurnal Genet Selektsii. 2022;26(1):40-9.
Shreberk-Shaked M, Dassa B, Sinha S, Di Agostino S, Azuri I, Mukherjee S, et al. A Division of Labor between YAP and TAZ in Non-Small Cell Lung Cancer. Cancer Res. 2020;80(19):4145-57.
Gallardo A, Molina A, Asenjo HG, Lopez-Onieva L, Martorell-Marugán J, Espinosa-Martinez M, et al. EZH2 endorses cell plasticity to non-small cell lung cancer cells facilitating mesenchymal to epithelial transition and tumour colonization. Oncogene. 2022;41(28):3611-24.
Lo Sardo F, Canu V, Maugeri-Saccà M, Strano S, Blandino G. YAP and TAZ: Monocorial and bicorial transcriptional co-activators in human cancers. Biochim Biophys Acta Rev Cancer. 2022;1877(4):188756.
Chung JH, Ostrowski MC, Romigh T, Minaguchi T, Waite KA, Eng C. The ERK1/2 pathway modulates nuclear PTEN-mediated cell cycle arrest by cyclin D1 transcriptional regulation. Hum Mol Genet. 2006;15(17):2553-9.
Gounder M, Schöffski P, Jones RL, Agulnik M, Cote GM, Villalobos VM, et al. Tazemetostat in advanced epithelioid sarcoma with loss of INI1/SMARCB1: an international, open-label, phase 2 basket study. Lancet Oncol. 2020;21(11):1423-32.
Morschhauser F, Tilly H, Chaidos A, McKay P, Phillips T, Assouline S, et al. Tazemetostat for patients with relapsed or refractory follicular lymphoma: an open-label, single-arm, multicentre, phase 2 trial. Lancet Oncol. 2020;21(11):1433-42.
Huang X, Yan J, Zhang M, Wang Y, Chen Y, Fu X, et al. Targeting Epigenetic Crosstalk as a Therapeutic Strategy for EZH2-Aberrant Solid Tumors. Cell. 2018;175(1):186-99.e19.
Elster D, Jaenicke LA, Eilers M, von Eyss B. TEAD activity is restrained by MYC and stratifies human breast cancer subtypes. Cell Cycle. 2016;15(19):2551-6.
Pang Y, Bai G, Zhao J, Wei X, Li R, Li J, et al. The BRD4 inhibitor JQ1 suppresses tumor growth by reducing c-Myc expression in endometrial cancer. J Transl Med. 2022;20(1):336.
Truica MI, Burns MC, Han H, Abdulkadir SA. Turning Up the Heat on MYC: Progress in Small-Molecule Inhibitors. Cancer Res. 2021;81(2):248-53.
Oku Y, Nishiya N, Shito T, Yamamoto R, Yamamoto Y, Oyama C, et al. Small molecules inhibiting the nuclear localization of YAP/TAZ for chemotherapeutics and chemosensitizers against breast cancers. FEBS Open Bio. 2015;5:542-9.
Bennett NK, Nakaoka HJ, Laurent D, Okimoto RA, Sei Y, Horvai AE, et al. Primary and metastatic tumors exhibit systems-level differences in dependence on mitochondrial respiratory function. PLoS Biol. 2022;20(9):e3001753.
Kushner MH, Ory V, Graham GT, Sharif GM, Kietzman WB, Thevissen S, et al. Loss of ANCO1 repression at AIB1/YAP targets drives breast cancer progression. EMBO Rep. 2020;21(1):e48741.
Kim M, Kim T, Johnson RL, Lim DS. Transcriptional co-repressor function of the hippo pathway transducers YAP and TAZ. Cell Rep. 2015;11(2):270-82.
Tumaneng K, Schlegelmilch K, Russell RC, Yimlamai D, Basnet H, Mahadevan N, et al. YAP mediates crosstalk between the Hippo and PI(3)K–TOR pathways by suppressing PTEN via miR-29. Nat Cell Biol. 2012;14(12):1322-9.
Wang T, Mao B, Cheng C, Zou Z, Gao J, Yang Y, et al. YAP promotes breast cancer metastasis by repressing growth differentiation factor-15. Biochim Biophys Acta Mol Basis Dis. 2018;1864(5 Pt A):1744-53.
Zanconato F, Forcato M, Battilana G, Azzolin L, Quaranta E, Bodega B, et al. Genome-wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growth. Nat Cell Biol. 2015;17(9):1218-27.
Lo Sardo F. Co-Immunoprecipitation (Co-Ip) in Mammalian Cells. Methods Mol Biol. 2023;2655:67-77.
Gertz J, Savic D, Varley KE, Partridge EC, Safi A, Jain P, et al. Distinct properties of cell-type-specific and shared transcription factor binding sites. Mol Cell. 2013;52(1):25-36.
Stein C, Bardet AF, Roma G, Bergling S, Clay I, Ruchti A, et al. YAP1 Exerts Its Transcriptional Control via TEAD-Mediated Activation of Enhancers. PLoS Genet. 2015;11(8):e1005465.

(Not answered)

Download PDF

Editorial decision: revise
13 Sep, 2024
Review #2 received at journal
07 Sep, 2024
Reviewer #2 agreed at journal
07 Sep, 2024
Review #1 received at journal
07 Sep, 2024
Reviewers invited by journal
07 Sep, 2024
Reviewer #1 agreed at journal
07 Sep, 2024
Submission checks completed at journal
06 Sep, 2024
First submitted to journal
05 Sep, 2024
Unknown event
05 Sep, 2024
Editor assigned by journal
04 Sep, 2024

You are reading this latest preprint version

The oncogenic axis YAP/MYC/EZH2 impairs PTEN tumor suppression activity enhancing lung tumorigenicity

Status:

Version 1

Abstract

Figures

Introduction

Results

YAP, PRC2 and MYC contribute to the transcriptional repression of tumor-suppressor genes in NSCLC

Oncogenic MYC cooperates with YAP and TAZ to repress transcription of PTEN

EZH2 and MYC inhibition activates YAP and TAZ targets

Discussion

Materials and Methods

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1