Lysine-specific demethylase 1 (LSD1) is a histone demethylase that contributes to the etiology of oral squamous cell carcinoma (OSCC) in part by promoting cancer stem cell phenotypes. The molecular signals regulated by LSD1, or acting with LSD1, are poorly understood, particularly in the development of OSSC. In this study, we show that conditional deletion of the Lsd1 gene or pharmacologic inhibition of LSD1 in the tongue epithelium leads to reduced development of OSCC following exposure to the tobacco carcinogen 4NQO. LSD1 inhibition attenuated proliferation and clonogenic survival and showed an additive effect when combined with the YAP inhibitor Verteporfin. Interestingly, LSD1 inhibition upregulated the expression of PD-L1, leading to immune checkpoint inhibitor therapy responses.

Implications:

Collectively, our studies reveal a critical role for LSD1 in OSCC development and identification of tumor growth targeting strategies that can be combined with LSD1 inhibition for improved therapeutic application.

Recent statistics show that oral cavity-related cancers account for 53,000 cases and 10,860 deaths in the United States annually (1, 2). The most common type of head and neck cancer is oral squamous cell carcinoma (OSCC). Treatment for OSCC frequently involves a combination of surgery, radiotherapy, and chemotherapy. Resistance to therapy complicates treatment effectiveness and the 5-year survival rate remains at ∼65% (3). In particular, epigenetic modifications can result in acquired resistance to treatment, leading to poor outcomes among patients with OSCC. However, the role of epigenetic mediators is not well understood, hindering further therapeutic development. One such aberrantly upregulated epigenetic regulator is Lysine-specific demethylase 1 (LSD1), which is encoded by the KDM1A gene. LSD1 is a flavin adenine dinucleotide—dependent amine oxidase. LSD1 also promotes cancer stem cells and chemoresistance (4, 5). We also showed that LSD1 inhibition attenuates OSCC and key signaling pathways, including EGFR signaling. EGFR is also implicated in OSCC progression and therapeutic resistance. The composition of the tumor microenvironment could attenuate immunosurveillance (6). Consistent with immune altering roles, ablation of LSD1 promotes antitumor immunity in combination with checkpoint blockade in melanoma (2).

LSD1 has also been shown to influence the Hippo signaling effector YAP (7). YAP is a transcriptional effector that has been shown to contribute to OSCC development, regulating pro-tumorigenic gene expression (8–10). While connections between LSD1 and pro-tumorigenic signaling effectors such as EGFR and YAP have been identified, the relationship between LSD1 and YAP in OSCC development and potential therapeutic application is not well understood. To gain insight into the roles of LSD1 in OSCC development, we examined LSD1 function using the tobacco smoke toxin, 4 Nitroquanolone-1-oxide (4NQO), model of OSCC. This model captures pathologic changes similar to those found in human OSCC (11–15).

Here, we evaluated the role of LSD1 in the genetic model in combination with 4NQO and determined if the topical application of the LSD1 inhibitor could attenuate its progression to OSCC. We evaluated the topical application of the LSD1 inhibitor in OSCC mice and further report that conditional deletion of Lsd1 in expanding oral epithelium reverses tumorigenic phenotypes in 4NQO treated mice. LSD1-inhibited epithelium revealed signaling networks impacted by LSD1 in tumor development, including YAP/TAZ and PD-L1 pathways, which play a critical role in OSCC. Combination therapies targeting these pathways showed improved targeting of 4NQO-induced OSCC, suggesting potential co-inhibition of these signals in OSCC therapy.

4NQO mouse model

All the experiments were performed with prior approval of the Institutional Animal Care and Use Committee, Boston University Medical Campus, Boston University. C57BL/6J mice were fed with 4NQO (100 μg/mL in propylene glycol and water) in their drinking water for 16 weeks, followed by regular water for the remainder of the experiment using a well-established protocol (16). The progressive changes at different stages of OSCC were determined at specific time points. To evaluate if LSD1 expression progressively increases in our models as in human OSCC, we performed an immunostaining analysis of tongue OSCC. All time points are hereafter indicated as weeks post-4NQO treatment.

LSD1 knockout mice

Lsd1-deficient mice were generated by crossing conditional floxed mice (loxP-Lsd1-loxP; ref. 17) with K14 promoter-driven tamoxifen-inducible Cre mice (K14CreERT, Jackson Laboratory, Stock No. 005107). Male and female mice were fed 4NQO in their drinking water. Given the early increase in LSD1 expression post-4NQO treatment (Supplementary Fig. S1), we treated mice at week 8 with tamoxifen to induce Lsd1 deletion. After 8 weeks of initiating 4NQO treatment, mice were randomized into two groups: half were treated with vehicle and half with tamoxifen. Tamoxifen was applied topically onto the tongue (75 mg/kg body weight, every 24 hours for 5 consecutive days, as per standard protocol; ref. 18) for tongue epithelium-specific deletion of Lsd1. K14 promoter-driven conditional Lsd1-floxed mice treated with vehicle were designated Lsd1WT/WT, whereas tamoxifen-treated mice with deleted Lsd1 were designated Lsd1−/− for homozygotes. The immunostaining and hematoxylin and eosin (H&E) staining were performed on n = 5–8/group and 4 sections/mouse.

RNA sequencing analysis

RNA sequencing (RNA-seq) was performed with 100 ng of total RNA using BGISEQ-500 (Beijing Genomic Institution, https://www.genomics.cn/, BGI) and published in other studies (19, 20). Raw FASTQ sequencing reads were mapped against the reference genome of Mus musculus (mm10). Total read counts and normalized transcripts per million were obtained. Before, differential expression, batch effects or sample heterogeneity was tested using iSeqQC (21). Differential gene expression was tested between Lsd1−/− and Lsd1WT/WT samples using DESeq2 package in R/Bioconductor. To evaluate the biological processes implicated in LSD1 deleted tongue, functional annotations were performed on the differentially expressed genes using DAVID (22). Finally, to characterize the modulation of the immune cell after the LSD1 deletion, we used the deconvolution methods as used in the reported study (23). The gene signatures for 29 different immune cell types were obtained from the study (24) and supervised gene set enrichment analysis using GSEA-v4.0 (25) was performed.

In vitro assays

The cell line authentication was performed by Genetica, Inc and routine Mycoplasma testing was performed using Mycoplasma Detection Kit from Southern Biotech (Catalog No. 13100–01). Low passage number HSC-3 and CAL-27 cells were grown in 100-mm3 plates for 48 to 72 hours or at confluency less than 90%. Next, proliferation assays were performed by plating HSC-3 and CAL-27 cells (20,000 cells per well, 6-replicates per treatment) treated for 48 hours with either vehicle control or with SP2509, then evaluated by CyQUANT Assay. Clonogenic survival was evaluated after 14 and 21 days (26). Most of the inhibitors used in our in vitro and in vivo studies are being used for clinical applications. SP2509 (MedChemExpress) is a reversible LSD1 enzyme activity inhibitor (27) and has similar properties to clinical compound SP-2577 (Seclidemstat; ref. 27). Verteporfin (Sellekchem, Inc) is a YAP inhibitor (28–30). EGFR kinase inhibitor Erolitinib (MedChemExpress; ref. 31) and TGFβ receptor type I (TGFβRI) kinase inhibitor Galunisertib (LY2157299; Sellekchem, Inc; ref. 32) are also currently being tested in clinical trials.

Mouse treatment with SP2509 and YAP inhibitor combination

Two studies were performed simultaneously with the different combinations. The first study was performed with 40 total C57BL/6 male and female mice with 4NQO for 18 weeks. They were randomized into 4 groups (n = 10/condition) and started treatment with: (i) vehicle, (ii) SP2509, (iii) Verteporfin, and (iv) SP2509 + Verteporfin. The route of drug administration for SP2509 and Verteporfin in 25-μL corn oil, 5% DMSO were topical applications onto the tongue after brief anesthesia with local distribution. The naïve group (n = 5) was without treatment of 4NQO as a control. We optimized a dose of SP2509 to 30 mg/kg and Verteporfin 2.5 mg/kg three times a week for 5 weeks by local application in DMSO and corn oil.

Mouse OSCC treatment with SP2509 and anti–PD-1 and anti–PD-L1 neutralizing antibody

In the second study, a total 60 C57BL/6 male and female mice treated with 4NQO for 18 weeks were randomized into different groups: (i) IgG isotype control antibody, (ii) SP2509, (iii) anti–PD-1 antibody, 4) SP2509+ anti–PD-1 antibody, 5) anti–PD-L11 antibody, and 6) SP2509+ anti–PD-L1 antibody. Anti-mouse PD-1 (CD279; Clone: 29F.1A12; Bio X Cell; ref. 33) and anti-mouse PD-l1 (B7-H1; Clone: 10F.9G2; Bio X Cell; ref. 34) were injected intraperitoneally in the respective groups at 10 mg/kg (100 μL buffer). SP2509 in 25-μL corn oil, 5% DMSO was applied topically onto the tongue.

RNA extraction and analysis

Total RNA was extracted by the TRizol protocol according to the manufacturer's instructions (Qiagen). RT-qPCR analysis was performed using TaqMan gene expression assays from Life Technologies, according to the standard protocol. The gene expression was normalized using three housekeeping genes: RSP18, GAPDH, and β-actin.

Pathologic characterization

Tongue sections were stained with H&E and evaluated for pathology by a Board Certified pathologist.

Clinical correlation of KDM1A expression

Pearson rho and Spearman rho correlation of KDM1A with other genes were evaluated using Xena UCSC genome browser (https://xena.ucsc.edu) to evaluate clinical databases such as the Therapeutically Applicable Research to Generate Effective Treatments (TARGET; 5968 samples). Pearson rho and Spearman's rank rho correlation were considered very high and positive for value 1, whereas the negative correlation was –1.

Genetic deletion of LSD1 attenuates OSCC

To gain insight into the dynamics of LSD1 in OSCC, we crossed Lsd1-floxed mice obtained from Stuart Orkin laboratory (Mass General Hospital, Boston; ref. 17) with mice expressing Cre recombinase from the Krt14 promoter tumor studies (35), allowing us to conditionally delete the Lsd1 gene in Krt14-expressing cells that expand in OSCC. We treated these animals and control animals with the tobacco carcinogen 4NQO (experimental outline is depicted in Fig. 1A), which induced OSCC that mirrors the progressive onset of the human disease (11–15). Lsd1−/− mice sacrificed at 24 weeks post-4NQO exposure showed a total absence or less severe gross pathologic development of OSCC compared with control Lsd1WT/WT mice (Fig. 1B). Gene expression changes identified by RT-qPCR demonstrated reduced expression of OSCC related genes that are known to be regulated by LSD1 (7), including reduced expression of Ccn1/Ctgf, Yap, Wwtr1/Taz, Egfr, as well as Lsd1 itself (Fig. 1C). The survival of Lsd1−/− mice was reduced compared with Lsd1WT/WT mice (Supplementary Fig. S2); however, this was not statistically signficant. The volume of tumors in the tongues of Lsd1−/− mice was reduced compared with Lsd1WT/WT mice (Fig. 1D) and H&E staining of tissues isolated from the tongues showed reduced OSCC development in Lsd1−/− mice at week 25 post-4NQO treatment compared with Lsd1WT/WT mice (Fig. 1E). Immunostaining analysis showed that LSD1 expression was reduced, whereas H3K4me1, which is the substrate of LSD1 enzyme, was increased (Fig. 1F).

Figure 1.

Deletion of Lsd1 in Krt14-expressing cells results in reduced OSCC following 4NQO exposure. A, Experimental design to evaluate the role of LSD1 in 4NQO mouse model. B, Tongue epithelial cells—specific deletion of Lsd1 shows the reduction of the gross pathologic lesion on the tongue. Tongue histopathologic changes evaluated after K14 promoter-specific deletion of LSD1 show the reduction of the pathologic lesion on the tongue in the drinking water starts at week 24. C, LSD1 deletion results in loss of Lsd1 expression and reduced expression of OSCC-related genes, including Ctgf, Yap, Egfr, and Wwtr1/Taz. Tumor volume measurements (D), H&E staining of tongue tissue (E), and immunostaining of tongue tissues at 25 weeks post-4NQO start (F). Analysis by one way-ANOVA. *, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001. Each data point represents the average gene expression represented from individual mice in each graph.

Figure 1.

Deletion of Lsd1 in Krt14-expressing cells results in reduced OSCC following 4NQO exposure. A, Experimental design to evaluate the role of LSD1 in 4NQO mouse model. B, Tongue epithelial cells—specific deletion of Lsd1 shows the reduction of the gross pathologic lesion on the tongue. Tongue histopathologic changes evaluated after K14 promoter-specific deletion of LSD1 show the reduction of the pathologic lesion on the tongue in the drinking water starts at week 24. C, LSD1 deletion results in loss of Lsd1 expression and reduced expression of OSCC-related genes, including Ctgf, Yap, Egfr, and Wwtr1/Taz. Tumor volume measurements (D), H&E staining of tongue tissue (E), and immunostaining of tongue tissues at 25 weeks post-4NQO start (F). Analysis by one way-ANOVA. *, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001. Each data point represents the average gene expression represented from individual mice in each graph.

Close modal

Genetic deletion of LSD1 broadly affects various signaling pathways and cell types

RNA-seq analysis of Lsd1−/− mouse tongues subjected to 4NQO showed differential expression of gene encoding factors associated with several oncogenic signaling pathways compared with Lsd1WT/WT mice tongues (Fig. 2AC). Pathway analysis of these differentially expressed genes suggested broad effects on the tumor microenvironment with a reduction in EGFR activity, JNK signaling, IL6 production, inflammatory responses, innate immune responses, and proteolysis. A comparative analysis Lsd1−/− gene expression with published signature from 29 different cell types (24) suggested that Lsd1−/− mouse tongues exhibit elevated neutrophils, T-regulatory cells, memory cells, and activated dendritic cells, along with a reduction in M2 macrophages, T helper cells, gamma-delta T cells, and eosinophils. Overall, candidate-based and RNA-seq analysis further suggested that Lsd1 deletion could affect YAP/TAZ-mediated effects, immune regulation, and immune checkpoints.

Figure 2.

RNA-seq analysis after deletion of Lsd1 in Krt14-expressing cells regulates oncogenic genes, signaling pathways and cell types in OSCC following 4NQO exposure. A, Heat map showing differences in the genes in 4NQO mouse model. B, Differential regulation of oncogenic pathways. C, Immune cell types.

Figure 2.

RNA-seq analysis after deletion of Lsd1 in Krt14-expressing cells regulates oncogenic genes, signaling pathways and cell types in OSCC following 4NQO exposure. A, Heat map showing differences in the genes in 4NQO mouse model. B, Differential regulation of oncogenic pathways. C, Immune cell types.

Close modal

Therapeutic LSD1 inhibitor application shows additive effects with YAP/TAZ inhibition in OSCC

Reversible inhibitors of Lsd1 have shown the ability to reduce the growth of several tumor cell models (27). Accordingly, we found that treatment of human OSCC HSC-3 cells with the LSD1 inhibitor SP2509 showed a dose-dependent reduction in proliferation (Fig. 3A) and clonogenic survival (Fig. 3B and C). Our observations suggested that LSD1 intersects with several signaling pathways in OSCC development (7), which prompted us to further explore these relationships for potential therapeutic applications. Specifically, we tested the effects of SP2509 together with inhibitors of YAP (Verteporfin), TGFβ (Galunasirtib), and EGFR (Erlotinib). Verteporfin and Erlotinib showed dose-dependent inhibition of HSC3 proliferation (Fig. 3DE) and Verteporfin showed more efficiency. Galunasirtib did not inhibit proliferation (Fig. 3F). We initially evaluated the combination of these inhibitors using proliferation assays with human HSC3 cells and found the additive effects of SP2509 and Verteporfin at 10:1 or 5:1 ratios showed more than 90% inhibition of HSC3 cell proliferation (Fig. 3G).

Figure 3.

LSD1 inhibitor SP2509 prevents proliferation alone or in combination. A, SP2509 shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. B and C, SP2509 shows dose-dependent inhibition of HSC3 and CAL27 cells clonogenic survival at days 14 and 21. D, YAP inhibitor (Verteporfin) shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. E, EGF inhibitor (erlotinib) shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. TGFβ1 inhibitor (galunasertib) does not show inhibition of HSC3 cell proliferation (F) and combination of SP2509 with Verteporfin shows efficient inhibition of cell proliferation in combination with SP2509: Verteporfin (10:1 and 5:1 ratio; G). Analysis by one way-ANOVA: *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 3.

LSD1 inhibitor SP2509 prevents proliferation alone or in combination. A, SP2509 shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. B and C, SP2509 shows dose-dependent inhibition of HSC3 and CAL27 cells clonogenic survival at days 14 and 21. D, YAP inhibitor (Verteporfin) shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. E, EGF inhibitor (erlotinib) shows dose-dependent inhibition of HSC3 cells proliferation at days 3 and 7. TGFβ1 inhibitor (galunasertib) does not show inhibition of HSC3 cell proliferation (F) and combination of SP2509 with Verteporfin shows efficient inhibition of cell proliferation in combination with SP2509: Verteporfin (10:1 and 5:1 ratio; G). Analysis by one way-ANOVA: *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

Our in vitro observations prompted us to test the combination of these inhibitors in vivo. For this, we first titrated the dose of Verteporfin for the local application onto the mouse tongue to any adverse effects. Our study showed that 2.5 to 3 mg/kg of Verteporfin was safe and did not induce any adverse effects. SP2509 25 mg/kg + 2.5 mg/kg Verteporfin in DMSO and corn oil emulsion was therefore tested in combination for in vivo studies using the 4NQO-induced model of OSCC progression (experimental outline depicted in Fig. 4A). Topical application of SP2509 alone inhibited pathologic changes, including reduced papillomatous lesions, SCC, and inflammatory cells on the tongue, as well as reduced OSCC growth, but the combination with Verteporfin showed striking additive effects compared with a single drug alone (Fig. 4BD). Analysis of lysates collected from tissues treated with these inhibitors showed reduced expression of Lsd1 (Fig. 4E).

Figure 4.

SP2509 in combination with YAP inhibitor shows efficient OSCC inhibition. A, Experimental design to evaluate the role LSD1 inhibitor and combination in 4NQO mouse model. B and C, The gross pathologic phenotype of tongue treated with SP2509, Verteporfin, or their combination showing SP2509 + Verteporfin inhibits phenotypic and pathologic changes on the tongue, D, Tongue tumor volume of mice treated with SP2509, Verteporfin, or their combination showing SP2509 + Verteporfin inhibits tumor volume, E, The qRT-PCR analysis Lsd1 and PD-1 expression in various groups. Analysis by one way-ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Figure 4.

SP2509 in combination with YAP inhibitor shows efficient OSCC inhibition. A, Experimental design to evaluate the role LSD1 inhibitor and combination in 4NQO mouse model. B and C, The gross pathologic phenotype of tongue treated with SP2509, Verteporfin, or their combination showing SP2509 + Verteporfin inhibits phenotypic and pathologic changes on the tongue, D, Tongue tumor volume of mice treated with SP2509, Verteporfin, or their combination showing SP2509 + Verteporfin inhibits tumor volume, E, The qRT-PCR analysis Lsd1 and PD-1 expression in various groups. Analysis by one way-ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Close modal

LSD1 inhibition promotes PD-L1

To evaluate the relationship between LSD1 and PD-L1, analysis of gene expression data available from The Cancer Genome Atlas (TCGA) studies showed an inverse correlation between Lsd1 and PD-L1 (Fig. 5A). Consistent with these observations, we observed an increase in the levels of Pd-l1 in SP2509 and SP2509/Verteporfin-treated tissues (Fig. 5B). Further, Lsd−/− mouse tongues also showed increased expression of Pd-l1 compared with LsdWT/WT mouse tongues (Fig. 5C). Finally, immunostaining analysis also showed that LSD1 pharmacologic inhibition (Fig. 5D) of LSD1 knockout (Fig. 5E) increased the levels of PD-L1 expression. Overall, above data indicates that Pd-L1 is repressed by LSD1 and that LSD11 deletion promotes the expression of PD-L1 (Fig. 5AE). These observations suggested the possibility that Lsd1 inhibition may promote efficacy towards anti–PD-1 and anti–PD-L1 immunotherapies.

Figure 5.

LSD1 inhibition promotes PD-L1 and has inverse relationship. A, Lsd1 and Pd-l1 have significantly higher negative correlation in TCGA data. B, The qRT-PCR analysis Lsd1 and Pd-1 expression in various groups treated with SP2509 or combination, C, PD-L1 expression in Lsd1−/− mice tongue compared with Lsd1WT/WT mice tongue. Immunostaining analysis of PD-L1 expression (D) and its quantification (E). Analysis by one-way ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Figure 5.

LSD1 inhibition promotes PD-L1 and has inverse relationship. A, Lsd1 and Pd-l1 have significantly higher negative correlation in TCGA data. B, The qRT-PCR analysis Lsd1 and Pd-1 expression in various groups treated with SP2509 or combination, C, PD-L1 expression in Lsd1−/− mice tongue compared with Lsd1WT/WT mice tongue. Immunostaining analysis of PD-L1 expression (D) and its quantification (E). Analysis by one-way ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Close modal

Therapeutic LSD1 inhibitor application shows additive effects with immune checkpoint antibodies in OSCC

To evaluate potential additive effects of the Lsd1 inhibitor SP2509 with immune checkpoint inhibitors, we treated mice exposed to 4NQO with SP2509 alone and in combination with anti–PD-1 and anti–PD-L1 (experimental outline is depicted in Fig. 6A). SP2509 was topically applied onto the tongue of developing tumors, and the mice were simultaneously administered anti–PD-1 or anti–PD-L1 neutralizing antibodies by intraperitoneal injection, a common administration route for these antibodies (33, 34). Monitoring of tumors showed reduced growth with SP2509, anti–PD-1, or anti–PD-L1 treatment alone and demonstrated a more effective reduction of tumor growth with combined SP2509 and anti–PD-1 treatment (Fig 6BD). The combination of SP2509 reduced gross pathologic and papillomatous lesions, normal basal epithelial morphology and inflammatory lesions, and Lsd1 expression (Fig. 6E). Immunostaining analysis also showed that LSD1 inhibition by SP2509 promoted the presence of CD8+ cells, suggesting elevated T-cell infiltration. Correlation analysis of LSD1/KDM1A in TCGA datasets (Supplementary Table S1) showed that Pearson correlation of LSD1 was moderate to high with EGFR (0.42), YAP (0.48), and CTGF (0.41). Spearman's rank rho was also moderate to the high of LSD1 with EGFR (0.43), YAP (0.51), and CTGF (0.31). Thus, we proposed a model from our studies that LSD1-driven feed-forward loop in OSCC (Fig. 6G).

Figure 6.

LSD1 regulates tumor immune microenvironment and sensitizes to anti–PD-1 and anti–PD-L1 antibody immunotherapy. A, Experimental design to evaluate the role LSD1 inhibitor and combination with immune checkpoint antibodies in 4NQO mouse model. B and C, The gross pathologic phenotype, H&E staining of tongue treated with SP2509, alone or in combination with anti–PD-1 and anti–PD-L1 neutralizing antibody, shows the combination shows more favorable response than a single drug alone to treat OSCC. D, Tongue tumor volume of mice topically treated with SP2509, alone or in combination with anti–PD-1 and anti–PD-L1 neutralizing antibody, shows SP2509+anti–PD-1 antibody treatment inhibits OSCC and inhibits tumor volume, E, The qRT-PCR analysis shows the SP2509 alone or and in combination with anti–PD-1 or anti–PD-L1 attenuates Lsd1 expression. F, Immunostaining analysis of tongue tissue sections with anti-LSD1 and anti-CD8 antibody with different treatments. G, The potential mechanism showing that LSD1 and YAP/TAZ coregulate PD-L1 as a feed-forward loop in OSCC; the disruption of this loop is critical for therapy. Analysis by one way-ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Figure 6.

LSD1 regulates tumor immune microenvironment and sensitizes to anti–PD-1 and anti–PD-L1 antibody immunotherapy. A, Experimental design to evaluate the role LSD1 inhibitor and combination with immune checkpoint antibodies in 4NQO mouse model. B and C, The gross pathologic phenotype, H&E staining of tongue treated with SP2509, alone or in combination with anti–PD-1 and anti–PD-L1 neutralizing antibody, shows the combination shows more favorable response than a single drug alone to treat OSCC. D, Tongue tumor volume of mice topically treated with SP2509, alone or in combination with anti–PD-1 and anti–PD-L1 neutralizing antibody, shows SP2509+anti–PD-1 antibody treatment inhibits OSCC and inhibits tumor volume, E, The qRT-PCR analysis shows the SP2509 alone or and in combination with anti–PD-1 or anti–PD-L1 attenuates Lsd1 expression. F, Immunostaining analysis of tongue tissue sections with anti-LSD1 and anti-CD8 antibody with different treatments. G, The potential mechanism showing that LSD1 and YAP/TAZ coregulate PD-L1 as a feed-forward loop in OSCC; the disruption of this loop is critical for therapy. Analysis by one way-ANOVA with post hoc Bonferroni correction: *, P < 0.05; **, P < 0.01; ***, P < 0.001. Each data point represents average gene expression represented from individual mice in each graph.

Close modal

Loss of tongue tumors by LSD1 deletion, the most striking finding, prompted us to evaluate gene expression changes and the potential therapeutic potential of pharmacologic inhibition of LSD1. LSD1 deletion reduced expression of OSCC-related genes, including Ctgf, Yap, Wwtr1/Taz, and Egfr. LSD1 has been shown to collaborate with YAP (7), a Hippo signaling effector that promotes phenotypes associated with OSCC (8) and other cancers (36, 37). Studies have further shown that both LSD1 and YAP localize to the NuRD chromatin remodeling complex and LSD1 is recruited to the chromatin by a YAP/TAZ-TEAD complex (38, 39). However, the regulatory relationship between LSD1 and YAP is unknown. Cross-talk between LSD1 and nuclear YAP/TAZ signals may affect common downstream transcriptional targets, such as CTGF and the PD-1/PD-L1 axis. PD-L1 expression directly protects tumors from cytolytic T cell killing (40), so the expression of PD-L1 within a tumor cell or tumor-infiltrating cells explains the therapeutic effect of PD-1/PD-L1 blocking antibodies (41–44). LSD1 has been reported to contribute to cancer initiation, progression, and relapse through multiple mechanisms (45–47). Our observations suggest that LSD1 has a specific role in reprogramming the immune microenvironment during OSCC. A previous study showed that LSD1 promotes the progression of a type of medulloblastoma (48), while the ablation of LSD1 promotes antitumor immunity in combination with checkpoint blockade in melanoma (2). Both tumor- and host-derived PD-L1 can play critical roles in immunosuppression (43) and our findings indicate that LSD1 deletion in K14 positive cells in 4NQO mouse could be direct or indirect PD-L1 expression on K14 positive and other cell types, and possibly not restricted to LSD1−/− cells.

Consistent with our study, the histone demethylase KDM4A has been shown to promote antitumor immunity and sensitivity to PD-1 antibody in OSCC growth and metastasis (49). Epigenetic therapies have primarily focused on enzymes that regulate histone acetylation and DNA methylation, with several DNA methyltransferase and histone deacetylase inhibitors now in clinical use and approved by the FDA to treat leukemias and other cancers. Histone lysine demethylases remain a relatively untapped source of potential “druggable targets” and, thus, the LSD1 signaling network may offer opportunities for exploiting immunotherapy. Indeed, we found that SP2509 and Verteporfin promote the expression of PD-L1, which could provide a more sustained effect for immunotherapy and combination therapies. Ablation of LSD1 promotes PD-L1 expression, and antitumor immunity is shown in breast cancer (2). Thus, our study validates this observation in OSCC and further identifies that LSD1 inhibition promotes PD-L1 expression. Our data also suggests that LSD1 inhibition by SP2509 promotes CD8+ T-cell infiltration. Further, our analyses also showed that LSD1 correlates with EGFR, YAP/TAZ, and PD-L1 expression and likely participates in a network of signals induced by these genes in specific clinical cancers (as shown in Fig. 6F).

This study has translational significance by leveraging the fact that compounds targeting LSD1 regulated genes/pathways are extensively studied and are in clinical trials. SP2509 and SP-2577 (Seclidemstat) have similar properties (27) and Seclidemstat is in clinical trials in patients with relapsed or refractory Ewing sarcoma (NCT03600649), advanced solid tumors (NCT03895684), and gynecologic cancers (NCT04611139). Verteporfin (Visudyne) disrupts YAP/TAZ-TEAD interactions and is used in combination with other drugs (28–30). It is also being tested in clinical trials for various other indications (NCT04590664, NCT03033225). Our study showed the effect and mechanism for future application in OSCC. Finally, there is no ongoing clinical trial for OSCC testing LSD1 inhibitors or their combination with YAP/TAZ-TEAD interactions as pharmacologic inhibitors or PD-1/PD-L1 checkpoint antibodies, so our work offers translational directions. Thus, our studies have identified a basic mechanism of LSD1 in OSCC. Overall, our findings could have translational importance in future clinical studies for smoking-induced oral cancer.

G.J. Hanna reports grants from Bristol-Myers Squibb outside the submitted work. No disclosures were reported by the other authors.

T. Alhousami: Conceptualization, data curation, formal analysis. M. Diny: Conceptualization, data curation, formal analysis. F. Ali: Conceptualization, data curation, formal analysis. J. Shin: Conceptualization, data curation, formal analysis. G. Kumar: Conceptualization, resources, data curation, software, formal analysis, investigation. V. Kumar: Software, formal analysis, visualization. J.D. Campbell: Software, formal analysis. V. Noonan: Data curation, writing–review and editing. G.J. Hanna: Conceptualization, writing–review and editing. G.V. Denis: Writing–review and editing. S. Monti: Funding acquisition, writing–review and editing. M.A. Kukuruzinska: Resources, writing–review and editing. X. Varelas: Funding acquisition, writing–review and editing. M.V. Bais: Conceptualization, resources, formal analysis, supervision, funding acquisition, validation, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing.

The authors acknowledge M.V. Bais was supported by NIH/NIDCR grants R21DE026892 and R03DE025274. M.A. Kukuruzinska, S. Monti, and X. Varelas were supported by NIH/NIDCR grant R01DE030350. X. Varelas was funded by NIH/NHLBI R01HL124392 and an American Cancer Society - Ellison New England Research Scholar Grant (RSG-17–138–01-CSM). The authors thank CTSA grant U54-TR001012 for their support of pilot studies.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Siegel
RL
,
Miller
KD
,
Jemal
A
.
Cancer statistics, 2019
.
CA Cancer J Clin
2019
;
69
:
7
34
.
2.
Sheng
W
,
LaFleur
MW
,
Nguyen
TH
,
Chen
S
,
Chakravarthy
A
,
Conway
JR
, et al
.
LSD1 ablation stimulates antitumor immunity and enables checkpoint blockade
.
Cell
2018
;
174
:
549
63
.
3.
Rothenberg
SM
,
Ellisen
LW
.
The molecular pathogenesis of head and neck squamous cell carcinoma
.
J Clin Invest
2012
;
122
:
1951
7
.
4.
Verigos
J
,
Karakaidos
P
,
Kordias
D
,
Papoudou-Bai
A
,
Evangelou
Z
,
Harissis
HV
, et al
.
The histone demethylase LSD1/KappaDM1A mediates chemoresistance in breast cancer via regulation of a stem cell program
.
Cancers (Basel)
2019
;
11
:
1585
.
5.
Lei
ZJ
,
Wang
J
,
Xiao
HL
,
Guo
Y
,
Wang
T
,
Li
Q
, et al
.
Lysine-specific demethylase 1 promotes the stemness and chemoresistance of Lgr5(+) liver cancer initiating cells by suppressing negative regulators of beta-catenin signaling
.
Oncogene
2015
;
34
:
3188
98
.
6.
Pitt
JM
,
Marabelle
A
,
Eggermont
A
,
Soria
JC
,
Kroemer
G
,
Zitvogel
L
.
Targeting the tumor microenvironment: removing obstruction to anticancer immune responses and immunotherapy
.
Ann Oncol
2016
;
27
:
1482
92
.
7.
Alsaqer
SF
,
Tashkandi
MM
,
Kartha
VK
,
Yang
YT
,
Alkheriji
Y
,
Salama
A
, et al
.
Inhibition of LSD1 epigenetically attenuates oral cancer growth and metastasis
.
Oncotarget
2017
;
8
:
73372
86
.
8.
Hiemer
SE
,
Zhang
L
,
Kartha
VK
,
Packer
TS
,
Almershed
M
,
Noonan
V
, et al
.
A YAP/TAZ-regulated molecular signature is associated with oral squamous cell carcinoma
.
Mol Cancer Res
2015
;
13
:
957
68
.
9.
Hiemer
SE
,
Szymaniak
AD
,
Varelas
X
.
The transcriptional regulators TAZ and YAP direct transforming growth factor beta-induced tumorigenic phenotypes in breast cancer cells
.
J Biol Chem
2014
;
289
:
13461
74
.
10.
Mohammad
HP
,
Smitheman
KN
,
Kamat
CD
,
Soong
D
,
Federowicz
KE
,
Van Aller
GS
, et al
.
A DNA hypomethylation signature predicts antitumor activity of LSD1 inhibitors in SCLC
.
Cancer Cell
2015
;
28
:
57
69
.
11.
Schoop
RA
,
Noteborn
MH
,
Baatenburg de Jong
RJ
.
A mouse model for oral squamous cell carcinoma
.
J Mol Histol
2009
;
40
:
177
81
.
12.
Suzuki
R
,
Kohno
H
,
Suzui
M
,
Yoshimi
N
,
Tsuda
H
,
Wakabayashi
K
, et al
.
An animal model for the rapid induction of tongue neoplasms in human c-Ha-ras proto-oncogene transgenic rats by 4-Nitroquinoline 1-oxide: its potential use for preclinical chemoprevention studies
.
Carcinogenesis
2006
;
27
:
619
30
.
13.
Czerninski
R
,
Amornphimoltham
P
,
Patel
V
,
Molinolo
AA
,
Gutkind
JS
.
Targeting mammalian target of rapamycin by rapamycin prevents tumor progression in an oral-specific chemical carcinogenesis model
.
Cancer Prev Res
2009
;
2
:
27
36
.
14.
Wang
Z
,
Wu
VH
,
Allevato
MM
,
Gilardi
M
,
He
Y
,
Callejas-Valera
LJ
, et al
.
Syngeneic animal models of tobacco-associated oral cancer reveal the activity of in situ anti–CTLA-4
.
Nat Commun
2019
;
10
:
5546
.
15.
Tang
XH
,
Osei-Sarfo
K
,
Urvalek
AM
,
Zhang
T
,
Scognamiglio
T
,
Gudas
LJ
.
Combination of bexarotene and the retinoid CD1530 reduces murine oral-cavity carcinogenesis induced by the carcinogen 4-Nitroquinoline 1-oxide
.
Proc Nat Acad Sci USA
2014
;
111
:
8907
12
.
16.
Mahjour
F
,
Dambal
V
,
Shrestha
N
,
Singh
V
,
Noonan
V
,
Kantarci
A
, et al
.
Mechanism for oral tumor cell lysyl oxidase like-2 in cancer development: synergy with PDGF-AB
.
Oncogenesis
2019
;
8
:
34
.
17.
Kerenyi
MA
,
Shao
Z
,
Hsu
YJ
,
Guo
G
,
Luc
S
,
O'Brien
K
, et al
.
Histone demethylase Lsd1 represses hematopoietic stem and progenitor cell signatures during blood cell maturation
.
Elife
2013
;
2
:
e00633
.
18.
Madisen
L
,
Zwingman
TA
,
Sunkin
SM
,
Oh
SW
,
Zariwala
HA
,
Gu
H
, et al
.
A robust and high-throughput Cre reporting and characterization system for the whole mouse brain
.
Nat Neurosci
2010
;
13
:
133
40
.
19.
Xu
J
,
Xu
X
,
Wang
B
,
Ma
Y
,
Zhang
L
,
Xu
H
, et al
.
Nuclear carbonic anhydrase 6B associates with PRMT5 to epigenetically promote IL12 expression in innate response
.
Proc Natl Acad Sci U S A
2017
;
114
:
8620
5
.
20.
Huang
Y
,
Xu
Z
,
Xiong
S
,
Sun
F
,
Qin
G
,
Hu
G
, et al
.
Repopulated microglia are solely derived from the proliferation of residual microglia after acute depletion
.
Nat Neurosci
2018
;
21
:
530
40
.
21.
Kumar
G
,
Ertel
A
,
Feldman
G
,
Kupper
J
,
Fortina
P
.
iSeqQC: a tool for expression-based quality control in RNA sequencing
.
BMC Bioinf
2020
;
21
:
56
.
22.
Huang da
W
,
Sherman
BT
,
Lempicki
RA
.
Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources
.
Nat Protoc
2009
;
4
:
44
57
.
23.
Mandal
R
,
Senbabaoglu
Y
,
Desrichard
A
,
Havel
JJ
,
Dalin
MG
,
Riaz
N
, et al
.
The head and neck cancer immune landscape and its immunotherapeutic implications
.
JCI Insight
2016
;
1
:
e89829
.
24.
Bindea
G
,
Mlecnik
B
,
Tosolini
M
,
Kirilovsky
A
,
Waldner
M
,
Obenauf
AC
, et al
.
Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer
.
Immunity
2013
;
39
:
782
95
.
25.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
,
Mukherjee
S
,
Ebert
BL
,
Gillette
MA
, et al
.
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci USA
2005
;
102
:
15545
50
.
26.
Bais
MV
,
Ozdener
GB
,
Sonenshein
GE
,
Trackman
PC
.
Effects of tumor-suppressor lysyl oxidase propeptide on prostate cancer xenograft growth and its direct interactions with DNA repair pathways
.
Oncogene
2015
;
34
:
1928
37
.
27.
Welch
D
,
Kahen
E
,
Fridley
B
,
Brohl
AS
,
Cubitt
CL
,
Reed
DR
.
Small molecule inhibition of lysine-specific demethylase 1 (LSD1) and histone deacetylase (HDAC) alone and in combination in Ewing sarcoma cell lines
.
PLoS One
2019
;
14
:
e0222228
.
28.
Morishita
T
,
Hayakawa
F
,
Sugimoto
K
,
Iwase
M
,
Yamamoto
H
,
Hirano
D
, et al
.
The photosensitizer verteporfin has light-independent antileukemic activity for Ph-positive acute lymphoblastic leukemia and synergistically works with dasatinib
.
Oncotarget
2016
;
7
:
56241
52
.
29.
Al-Moujahed
A
,
Brodowska
K
,
Stryjewski
TP
,
Efstathiou
NE
,
Vasilikos
I
,
Cichy
J
, et al
.
Verteporfin inhibits growth of human glioma in vitro without light activation
.
Sci Rep
2017
;
7
:
7602
.
30.
Gibault
F
,
Corvaisier
M
,
Bailly
F
,
Huet
G
,
Melnyk
P
,
Cotelle
P
.
Non-photoinduced biological properties of verteporfin
.
Curr Med Chem
2016
;
23
:
1171
84
.
31.
Molina-Arcas
M
,
Moore
C
,
Rana
S
,
van Maldegem
F
,
Mugarza
E
,
Romero-Clavijo
P
, et al
.
Development of combination therapies to maximize the impact of KRAS-G12C inhibitors in lung cancer
.
Sci Transl Med
2019
;
11
:
eaaw7999
.
32.
Bueno
L
,
de Alwis
DP
,
Pitou
C
,
Yingling
J
,
Lahn
M
,
Glatt
S
, et al
.
Semi-mechanistic modelling of the tumor growth inhibitory effects of LY2157299, a new type I receptor TGFβ kinase antagonist, in mice
.
Eur J Cancer
2008
;
44
:
142
50
.
33.
Gordon
SR
,
Maute
RL
,
Dulken
BW
,
Hutter
G
,
George
BM
,
McCracken
MN
, et al
.
PD-1 expression by tumor-associated macrophages inhibits phagocytosis and tumor immunity
.
Nature
2017
;
545
:
495
9
.
34.
Stathopoulou
C
,
Gangaplara
A
,
Mallett
G
,
Flomerfelt
FA
,
Liniany
LP
,
Knight
D
, et al
.
PD-1 inhibitory receptor downregulates asparaginyl endopeptidase and maintains foxp3 transcription factor stability in induced regulatory T cells
.
Immunity
2018
;
49
:
247
63
.
35.
Caulin
C
,
Nguyen
T
,
Longley
MA
,
Zhou
Z
,
Wang
XJ
,
Roop
DR
.
Inducible activation of oncogenic K-ras results in tumor formation in the oral cavity
.
Cancer Res
2004
;
64
:
5054
8
.
36.
Wang
G
,
Lu
X
,
Dey
P
,
Deng
P
,
Wu
CC
,
Jiang
S
, et al
.
Targeting YAP-dependent MDSC infiltration impairs tumor progression
.
Cancer Discov
2016
;
6
:
80
95
.
37.
Zanconato
F
,
Cordenonsi
M
,
Piccolo
S
.
YAP/TAZ at the roots of cancer
.
Cancer Cell
2016
;
29
:
783
803
.
38.
Wang
Y
,
Zhang
H
,
Chen
Y
,
Sun
Y
,
Yang
F
,
Yu
W
, et al
.
LSD1 is a subunit of the NuRD complex and targets the metastasis programs in breast cancer
.
Cell
2009
;
138
:
660
72
.
39.
Kim
M
,
Kim
T
,
Johnson
RL
,
Lim
DS
.
Transcriptional co-repressor function of the hippo pathway transducers YAP and TAZ
.
Cell Rep
2015
;
11
:
270
82
.
40.
Lau
J
,
Cheung
J
,
Navarro
A
,
Lianoglou
S
,
Haley
B
,
Totpal
K
, et al
.
Tumor and host cell PD-L1 is required to mediate suppression of antitumor immunity in mice
.
Nat Commun
2017
;
8
:
14572
.
41.
Kleinovink
JW
,
van Hall
T
,
Ossendorp
F
,
Fransen
MF
.
PD-L1 immune suppression in cancer: tumor cells or host cells?
Oncoimmunology
2017
;
6
:
e1325982
.
42.
Noguchi
T
,
Ward
JP
,
Gubin
MM
,
Arthur
CD
,
Lee
SH
,
Hundal
J
, et al
.
Temporally distinct PD-L1 expression by tumor and host cells contributes to immune escape
.
Cancer Immunol Res
2017
;
5
:
106
17
.
43.
Juneja
VR
,
McGuire
KA
,
Manguso
RT
,
LaFleur
MW
,
Collins
N
,
Haining
WN
, et al
.
PD-L1 on tumor cells is sufficient for immune evasion in immunogenic tumors and inhibits CD8 T-cell cytotoxicity
.
J Exp Med
2017
;
214
:
895
904
.
44.
Kleinovink
JW
,
Marijt
KA
,
Schoonderwoerd
MJA
,
van Hall
T
,
Ossendorp
F
,
Fransen
MF
.
PD-L1 expression on malignant cells is no prerequisite for checkpoint therapy
.
Oncoimmunology
2017
;
6
:
e1294299
.
45.
Zhang
X
,
Lu
F
,
Wang
J
,
Yin
F
,
Xu
Z
,
Qi
D
, et al
.
Pluripotent stem cell protein Sox2 confers sensitivity to LSD1 inhibition in cancer cells
.
Cell Rep
2013
;
5
:
445
57
.
46.
Huang
J
,
Sengupta
R
,
Espejo
AB
,
Lee
MG
,
Dorsey
JA
,
Richter
M
, et al
.
p53 is regulated by the lysine demethylase LSD1
.
Nature
2007
;
449
:
105
8
.
47.
Wang
J
,
Hevi
S
,
Kurash
JK
,
Lei
H
,
Gay
F
,
Bajko
J
, et al
.
The lysine demethylase LSD1 (KDM1) is required for maintenance of global DNA methylation
.
Nat Genet
2009
;
41
:
125
9
.
48.
Lee
C
,
Rudneva
VA
,
Erkek
S
,
Zapatka
M
,
Chau
LQ
,
Tacheva-Grigorova
SK
, et al
.
Lsd1 as a therapeutic target in Gfi1-activated medulloblastoma
.
Nat Commun
2019
;
10
:
332
.
49.
Zhang
W
,
Liu
W
,
Jia
L
,
Chen
D
,
Chang
I
,
Lake
M
, et al
.
Targeting KDM4A epigenetically activates tumor cell–intrinsic immunity by inducing DNA replication stress
.
Mol Cell
2021
;
81
:
2148
2165.e9
.

Supplementary data