1. Autophagy is structurally and functionally impaired in HPV + HNC cell lines
To assess the differences in autophagy activation between HPV + and HPV- HNC cells, we examined the number of autophagic vacuoles (AVs) (which include autophagosomes, autolysosomes and amphisomes) by electron microscopy (EM) in 2 HPV- and 2 HPV + HNC cell lines (Fig. 1A). HPV + HNC cell lines showed a reduced number and size of AVs compared to HPV- ones, suggesting a possible inhibition of the pathway (Fig. 1B, S1A).
We then assessed the levels of autophagy markers in HNC cell lines at basal conditions. HPV + HNC cell lines showed lower levels of the lipidated form of LC3 (LC3-II) compared to HPV- HNC cell lines (Fig. 1C, 1D), while other autophagy markers like Beclin or p62 were only respectively slightly or not affected (Figure S1B, S1C). We also measured the activation of Akt and mTORC, described to inhibit autophagy (33), by detecting Akt and p70 S6 kinase (S6K) phosphorylation by western blot in HNC cell lines (Figure S1D, S1E). While Akt phosphorlylation at S473 seems not to be related to HPV status, p70 (S6K) phosphorylation can be detected in the UD-SCC-2 and UPCI-SCC-90 HPV + HNC cell lines, suggesting that mTORC activity could promote autophagy inhibition in some HPV + HNC cell lines.
To functionally characterize autophagy status, we monitored the autophagic flux by using a tandem tagged LC3-GFP-RFP fusion protein (34). GFP is quenched in acidic pH such that only red dots can be observed in functional autolysosomes, while in neutral pH red and green colocalized dots can be detected. The HPV- UM-SCC-19 and the HPV + UD-SCC-2 HNC cell lines transduced with the lentiviral LC3-GFP-RFP vector were treated with the late stage autophagic inhibitor Bafilomycin A1, the mTORC1 inhibitor Torin-1, or the vehicle (DMSO) and observed by confocal microscopy (Fig. 1E). Quantification of colocalized GFP and RFP dots showed a significant increase upon autophagy inhibition and a significant decrease upon autophagy activation with respect to the vehicle only in the HPV- HNC cell lines, as expected in cells in which autophagy is functional (Fig. 1F, left). However, the levels of colocalized dots in the HPV + HNC cell lines do not vary, or increase, with respect to the vehicle, upon pharmacological induction or inhibition of autophagy (Fig. 1F, right). Taken together, these results suggest a global impairment in autophagic functions in HPV + HNC cell lines.
2. Hpv16 Oncoproteins Transcriptionally Inhibit Autophagy Gene Expression
To test whether HPV16 oncoproteins are involved in autophagy inhibition in HPV + HNC cells, we knocked down E6/E7 by RNAi in UD-SCC-2 cells. We then stained Lysosome-associated membrane protein 2 (LAMP2) (lysosomal) and LC3 (autophagosomal) proteins, revealing an increase in colocalization upon reduction of E6/E7, suggesting a reactivation of the pathway (Fig. 2A, 2B). This result was confirmed also by a significant increase in the LC3II/LC3I ratio upon siE6/E7 (Fig. 2C, 2D). Interestingly, downregulation of E6 and E7 was able to induce the mRNA upregulation of some autophagic and lysosomal genes (Fig. 2E).
To investigate the role of E6 and E7 in mRNA regulation of autophagy and lysosomal biogenesis, we transduced primary human keratinocytes (PHK) with HPV16 E6 and E7 (Fig. 2F, S2A), followed by RNA sequencing. Gene ontology analysis on genes downregulated by E6/E7 indicated genes belonging to lysosomal and vacuolar cellular compartments as some of the most relevant genes affected by the two oncoproteins (Fig. 2G), in addition, the pathway analysis highlighted lysosome and autophagy pathway among the top regulated hits (Fig. 2H). Moreover, only HPV16 E7 alone, and not E6, seems to induce mRNA downregulation of autophagy and lysosomal biogenesis (Figure S2D, S2E). These results were also confirmed by qPCR analysis (Figure S2B, S2C).
Next, we wanted to confirm these results also in HNC cell lines by monitoring the mRNA expression of autophagy and lysosomal genes in two HPV- (UM-SCC-4, UM-SCC-19) and two HPV+ (UD-SCC-2, UPCI-SCC-152) HNC cell lines at basal conditions. Figure 2I shows that many autophagy (the LC3 gene MAP1LC3B, the p62 gene SQSTM1, BECN1, GABARAP, GABARAPL1) and, more strongly, lysosomal (NEU1, CTSA, CTSB, CTSD, ATP6V0E1, LAMP1, LAMP2) genes are differentially expressed between HPV- and HPV + cell lines. These data suggest that HPV16 oncoproteins, and in particular E7, transcriptionally downregulate autophagy and lysosomal biogenesis.
3. c-MYC competes with TFEB and TFE3 for the binding to autophagy and lysosomal genes in HPV + HNC cells
As transcriptional repression of autophagy and lysosomal biogenesis occurred in HPV + cells, we investigated the molecular mechanism underlying this process. We evaluated the mRNA and protein expression (Figure S3A-C, Fig. 3A) of the most important transcription factors involved in lysosomal biogenesis: TFEB and TFE3, as well as c-MYC, recently found to repress autophagy (35). Western blots performed on HNC cell lines did not show any significant difference in the expression of these genes between HPV- and HPV + cell lines (Fig. 3A).
Thus, we focused on the differential activation of the transcription factors, and on their binding to the promoters of autophagy genes. We performed chromatin immunoprecipitation of TFEB, TFE3, and c-MYC in the UM-SCC-4 and UD-SCC-2 cell lines, followed by qPCR analysis of autophagy and lysosomal genes. We found that the binding of c-MYC to the promoter of LC3 gene (MAP1LC3B) was significantly higher in the HPV + cell line (UD-SCC-2) compared to the HPV- cell line (UM-SCC-4), and this corresponded to a significant decrease in the binding of TFEB and TFE3 to the same promoter (Fig. 3B). The same result was obtained also for GLA and COLN1 promoters (Fig. 3C, 3D), while no significant differences could be appreciated for the promoter of p62 gene (SQSTM1) (Fig. 3E). This is in line with Annunziata et al. (35) study, which showed that c-MYC competes with TFE3 and TFEB, repressing lysosomal and autophagic transcription. Our data suggest that HPV + HNC cells present a more active c-MYC compared to HPV- cells.
4. C-myc Stability Increases In Hpv + hnc Cells, Affecting The Autophagy Pathway
To assess whether c-MYC is involved in the inhibition of autophagy, we performed c-MYC knock-down in HPV + HNC cell lines. c-MYC downregulation affected both protein expression of LC3 and the mRNA levels of different autophagy genes (Figs. 4A, 4B, S4A, S4B, S4D). Moreover, c-MYC overexpression UD-SCC-2 and UPCI-SCC-152 cell lines was able to affect LC3 protein levels and to induce the downregulation of autophagy genes, confirming the involvement of c-MYC in autophagy repression (Figs. 4C, S4C, 4D).
Since no significant difference in the mRNA expression of MYC between HPV + and HPV- HNC cell lines was detected (Figure S3C), we evaluated the stability of c-MYC protein by treating HNC cell lines with cycloheximide (CHX), a widely used inhibitor of protein synthesis, at different time points. Since c-MYC is very unstable, we used 5 minutes as the first time point and 90 minutes as the last. Western blot analysis revealed a higher stability of c-MYC in HPV + with respect to HPV- HNC cell lines (Fig. 4E, S4E, S4F). We calculated the half-life of c-MYC in 3 HPV- and 3 HPV + HNC cell lines (Fig. 4F) and we observed a significant higher half-life in the HPV + HNC cell lines, revealing a higher stability of c-MYC in these cell lines (Fig. 4G). Moreover, E6/E7 silencing in the HPV + UD-SCC-2 cell line induced the downregulation of c-MYC protein level (Figure S4G), suggesting that the stability of c-MYC relies on E6/E7 oncoproteins.
Since c-MYC stability is affected by its phosphorylation, we tried to assess whether the phosphatase inhibitor was responsible for c-MYC stabilization in the UD-SCC-2 cell line. Silencing of CIP2A by shRNA induced a strong downregulation of c-MYC protein levels (Fig. 4H), which was specifically imputable to the reduction of c-MYC half-life (Fig. 4I). These results confirmed our hypothesis that, also in our system, c-MYC protein stability is regulated by its phosphorylation controlled by CIP2A expression.
5. Cip2a Inhibits Autophagy And Lysosomal Biogenesis In Hnc
As CIP2A affected c-MYC stability, we investigated the role of CIP2A mRNA expression in HNC tumors. We took advantage of the Pancancer Dataset of Head and Neck cancer tumors from the Cancer Genome Atlas (TCGA) and assessed the correlation between CIP2A mRNA and autophagy or lysosomal genes expression. We found that many genes involved in autophagy and lysosomal biogenesis inversely correlated with CIP2A expression (Fig. 5A, S5A, S5B, S5C). We then took the top 100 genes that inversely correlated with CIP2A expression and performed a gene ontology analysis, looking at cellular components, biological processes and pathways enriched (Figs. 5B, 5C, 5D). As shown in Fig. 5B-D, among the most enriched hits we found lysosome related genes, macroautophagy, mitophagy and selective autophagy.
We then decided to assess whether CIP2A expression affects autophagy in HPV + HNC cell lines. Upon CIP2A silencing in HPV + HNC cell lines we showed the upregulation of LC3-II (Fig. 5E and 5F), and also the increased mRNA LC3 gene expression (Fig. 5G). Moreover, the overexpression of CIP2A was able to reduce the colocalization of LC3 with LAMP2 (Figs. 5H, 5I). Taken together these results evidence a novel and important role of CIP2A expression in autophagy gene expression in HNC.
6. Hpv16 E6/e7 Induce Cip2a Overexpression In Hnc
After demonstrating the importance of CIP2A in c-MYC stability and autophagy in HNC, we evaluated the expression of CIP2A in HPV + vs HPV- HNC cell lines. Upon assessing mRNA expression of CIP2A in the Pancancer dataset, we observed significant higher expression levels in HPV + compared to HPV- HNC samples. (Fig. 6A). Interestingly, the same result was confirmed also in primary tumors obtained from IEO hospital (Fig. 6B), and in HNC cell lines (Fig. 6C).
We then evaluated the possible regulation of HPV16 E6/E7 on CIP2A expression, by transducing PHK with E6/E7 (Fig. 6D). Western blot analysis showed the upregulation of CIP2A in the presence of E6 and E7 proteins, confirmed by higher mRNA levels of CIP2A (Figure S6D).
To strengthen our observations, we checked CIP2A protein levels in HPV + and HPV- HNC cell lines, confirming significant higher levels of CIP2A in the HPV + subgroup (Fig. 6E, 6F), and that this upregulation can be significantly reverted upon knock-down of E6/E7 (Fig. 6G, 6H).
These results are in line with published literature showing an overexpression of CIP2A by high risk HPV due to E7 inhibition of pRb in other systems (36–38). Along the same line, we also observed the overexpression of E2F1 (regulated by pRb) both in TCGA cases and in HNC cell lines (Figure S6A, S6C). Moreover, gene set enrichment analysis (GSEA) performed on RNAseq data from PHK transduced with E6/E7 showed a significant enrichment in E2F1 pathways and E2F targets, as well as MYC upregulation, as expected (S6E, S6F, S6G). The same analysis did not show any significant enrichment in mTORC signaling (Figure S6H, S6I). All these data suggest that HPV indirectly controls CIP2A expression activating E2F1 signaling upon pRb inhibition.
7. Cip2a And C-myc Downregulation Confer Resistance To Cisplatin In Hpv + hnc Cell Lines
Autophagy has a fundamental role in resistance to chemotherapy (39–41). Recent findings suggest that inhibition of autophagy sensitizes cells to Cisplatin treatment (42, 43). It has been described that HPV + HNC cells are more sensitive to Cisplatin than HPV- cells (44). To understand the role of c-MYC and CIP2A inhibition in chemosensitivity, we evaluated Cisplatin IC50 in the HPV + UD-SCC-2 cell line upon knock-down of c-MYC or CIP2A. Figure 7A and 7B show that the downregulation of both c-MYC and CIP2A induced resistance to Cisplatin, with almost a doubling of the IC50 values in UD-SCC-2 cells.
To confirm these results, we downregulated c-MYC pharmacologically, by treating UD-SCC-2 and UPCI-SCC-152 cells with the Brd4 inhibitor JQ1, already known to activate autophagy (45, 46). We also tested the effects of the PP2A activator DT-061, to abrogate the impact of CIP2A overexpression on c-MYC (47). Both treatments were able to downregulate c-MYC levels (Figure S7C, S7D). JQ1 efficiently downregulated mRNA MYC expression and to upregulate autophagy in UD-SCC-2 and UPCI-SCC-152 cell lines (S7A, S7B). We then assessed the IC50 of Cisplatin in the presence of 1µM JQ1, 1µM DT-061 or vehicle (DMSO). Figure 7C and 7D showed an increased resistance to Cisplatin in comparison to the control in UD-SCC-2 cells. These results were confirmed also in the UPCI-SCC-152 cell line (Fig. 7E, 7F).
Taken together these results suggest the importance of c-MYC activity in chemosensitivity in HPV + HNC cells.