Article Text

Original research
SP140 inhibits STAT1 signaling, induces IFN-γ in tumor-associated macrophages, and is a predictive biomarker of immunotherapy response
  1. Kranthi Kiran Kishore Tanagala1,2,
  2. Joshua Morin-Baxter1,2,
  3. Richard Carvajal3,
  4. Maryum Cheema1,2,
  5. Sunil Dubey1,2,
  6. Hiroshi Nakagawa4,
  7. Angela Yoon5,
  8. Yi-Shing L Cheng6,
  9. Alison Taylor7,
  10. Jeffrey Nickerson8,
  11. Akiva Mintz1,9 and
  12. Fatemeh Momen-Heravi1,2,10
  1. 1Section of Oral, Diagnostic and Rehabilitation Sciences, Columbia University College of Dental Medicine, Columbia University Irving Medical Center, New York City, New York, USA
  2. 2Cancer Biology and Immunology Laboratory, College of Dental Medicine, Columbia University Irving Medical Center, New York, New York, USA
  3. 3Herbert Irving Comprehensive Cancer Center, Columbia University Medical Center, New York City, New York, USA
  4. 4Division of Digestive and Liver Diseases, Department of Medicine, Vagelos College of Physicians and Surgeons, Herbert Irving Comprehensive Cancer Center, New York, NY, USA
  5. 5Department of Stomatology, Division of Diagnostic science and Services, College of Dental Medicine, Medical University of South Carolina, Charleston, South Carolina, USA
  6. 6Department of Diagnostic Sciences, Texas A&M University System, Dallas, Texas, USA
  7. 7Department of Pathology and Cell Biology, Vagelos College of Physicians and Surgeons, Columbia University Irving Medical Center, New York, NY, 10032
  8. 8Division of Genes & Development, Department of Pediatrics, University of Massachusetts Chan Medical School, Worcester, Massachusetts, USA
  9. 9Department of Radiology, Columbia University Medical Center, New York, NY, 10032
  10. 10Herbert Irving Comprehensive Cancer Center, Columbia University Irving Medical Center, New York, NY, USA
  1. Correspondence to Dr Fatemeh Momen-Heravi; f.m.heravi{at}gmail.com; fm2540{at}cumc.columbia.edu

Abstract

Background Understanding the role and potential therapeutic targeting of tumor-associated macrophages (TAMs) is crucial to developing new biomarkers and therapeutic strategies for cancer immunotherapies. The epigenetic reader SP140 has emerged as a master regulator of macrophage transcriptional programs; however, its role in the signaling of TAMs and response to immunotherapy has not been investigated.

Methods We evaluated the correlation between SP140 expression in head and neck squamous cell carcinoma (HNSCC) TAMs and clinical outcomes. We also used complementary bioinformatics and experimental approaches to study the association of SP140 expression with tumor mutation burden, patient survival, immunogenic signature of tumors, and signaling of TAMs. SP140 overexpression or knockdown was implemented to identify the role of SP140 in downstream signaling and production of inflammatory cytokine and chemokines. Chromatin immunoprecipitation and analysis of assay of transposase accessible chromatin sequencing data were used to demonstrate the direct binding of SP140 on the promoters of STAT1. Finally, correlation of SP140 with immune cell infiltrates and response to immune-checkpoint blockade in independent cohorts of HNSCC, metastatic melanoma, and melanoma was assessed.

Results We found that SP140 is highly expressed in TAMs across many cancer types, including HNSCCs. Interestingly, higher expression of SP140 in the tumors was associated with higher tumor mutation burden, improved survival, and a favorable response to immunotherapy. Tumors with high SP140 expression showed enrichment of inflammatory response and interferon-gamma (IFN-γ) pathways in both pan-cancer analysis and HNSCC-specific analysis. Mechanistically, SP140 negatively regulates transcription and phosphorylation of STAT1 and induces IFN-γ signaling. Activating SP140 in macrophages and TAMs induced the proinflammatory macrophage phenotype, increased the antitumor activity of macrophages, and increased the production of IFN-γ and antitumor cytokines and chemokines including interleukin-12 and CXCL10. SP140 expression provided higher sensitivity and specificity to predict antiprogrammed cell death protein 1 immunotherapy response compared with programmed death-ligand 1 in HNSCCs and lung cancer. In metastatic melanoma, higher levels of SP140 were associated with a durable response to immunotherapy, higher immune score estimates, high infiltrations of CD8+ T cells, and inflammatory TAMs.

Conclusions Our findings suggest that SP140 could serve as both a therapeutic target and a biomarker to identify immunotherapy responders.

  • head and neck neoplasms
  • biomarkers, tumor
  • immunotherapy
  • macrophages

Data availability statement

Data are available in a public, open access repository. Data are available in a public, open-access repository from the following sources: https://gdc.cancer.gov/node/977 ; https://www.cbioportal.org/datasets; The Gene Expression Omnibus under accession number GSE121810, GSE181815, GSE89177, GSM1625982, GSM1625981.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Immune-checkpoint blockade is beneficial in a subset of patients. Combination of immune-checkpoint blockade with other therapies to improve clinical outcomes and targeting other immune cells such as tumor-associated macrophages (TAMs) are ongoing. There is an urgent need to investigate prognostic and predictive biomarkers for patient selection to improve immunotherapy response in different tumors.

WHAT THIS STUDY ADDS

  • In both pan-cancer and head and neck cancer analysis, high expression of SP140 in the tumors was associated with higher tumor mutation burden, improved survival, and a favorable response to immunotherapy. Mechanistically, SP140 negatively regulates transcription and phosphorylation of STAT1 and induces interferon gamma (IFN-γ) signaling. Activating SP140 in macrophages and TAMs induced the proinflammatory macrophage phenotype, increased the antitumor activity of macrophages, and increased the production of IFN-γ and antitumor cytokines and chemokines including interleukin-12 and CXCL10.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Our findings prove for the first time that SP140 potentially serves as a prognostic marker for immunotherapy and a therapeutic target for sensitizing tumors to immunotherapy.

Introduction

Changes in the tumor microenvironment (TME) and especially changes in tumor-associated immune cells play an important role in both the development and progression of solid tumors including head and neck squamous cell carcinoma (HNSCC).1 2 In the presence of tumor-cell antigens, the immune system triggers a natural T-cell response via cytotoxic T lymphocytes. In contrast, tumor cells promote immunosuppression by various tumor immune evasion mechanisms including overexpression of programmed death-ligand 1 (PD-L1) in tumor cells, immune infiltration by programmed cell death protein 1 (PD-1) T lymphocytes, and recruitment and activation of macrophages.3

Macrophages associated with the TME are referred to as tumor-associated macrophages (TAMs). They usually reside in close proximity to tumor cells and may represent up to 50% of the tumor mass.4 Historically, macrophages were phenotypically classified into two types: proinflammatory phenotype (M1) and anti-inflammatory phenotype (M2). M1 macrophages demonstrate antitumor functions, produce proinflammatory cytokines, phagocytize microbes, and initiate an immune response. On the other hand, tumor cells activate or switch the macrophages to an M2-like phenotype via several pathways such as CCL-2, interleukin (IL)-6, CSF-1, PD-1/PD-L1, CD47/SIRPα, which support tumor growth.5 Recent reports suggest that a binary classification of M1 and M2 phenotypes cannot capture the complexity of macrophage signaling and phenotypes of the TME. Most TAMs demonstrate a mixture of both phenotypes with a dominance of anti-inflammatory and protumorigenic activity.6 7 A subset of TAMs in the TME express immune-checkpoint modulators including PD-L1, which directly suppress the T-cell activation and have been reported to be associated with a better response to immune-checkpoint blockade.8 However, the complete network of macrophage signaling responsible for the immunosuppressive tumor phenotype, potential molecules to target these pathways, and markers of tumor vulnerabilities in TAMs for patient selection have not been completely identified.

Speckled protein 140 (SP140) is a nuclear protein that belongs to the SP family, consisting of SP100, SP110, and SP140L; SP family proteins have high sequence homology with the autoimmune regulator and are highly expressed in innate and adaptive immune cells with SP140 being the immune restricted.9–11 SP140 is interferon gamma (IFN-γ) inducible, and variant forms or differential expression of SP140 have been associated with immune disorders such as Crohn’s disease,12 multiple sclerosis,13 14 inflammatory bowel disease,12 and biliary cirrhosis.15 Previous reports demonstrated that SP140 is a transcriptional repressor that plays a pivotal role in the transcriptional programming of macrophages by suppressing lineage-inappropriate genes and is essential for lipopolysaccharide (LPS)-induced macrophage transcriptional programs.12 However, SP140 roles in the pathogenesis of cancer, response to immunotherapy, as well as TAMs are not investigated.

In the present study, we investigated the role of SP140 and its downstream signaling in HNSCC and macrophages. High levels of SP140 augmented the antitumorigenic macrophage phenotype by inhibiting STAT1 at the transcription and post-transcription level and induction of IFN-γ-related genes. High levels of SP140 were associated with high immune infiltrate in tumors in HNSCC and pan-cancer studies. In patients who underwent immunotherapy with immune-checkpoint inhibitors, high levels of SP140 were associated with response to immunotherapy and increased overall survival in HNSCCs, lung cancer, and metastatic melanoma. SP140 expression in the tumor conferred higher diagnostic accuracy compared with PD-L1 for predicting response to anti-PD-1 immunotherapy.

Materials and methods

Clinical samples

HNSCC samples were collected at the time of surgical resection and reviewed by pathologists at Columbia University. After collection from surgical resection, all samples were embedded in formalin or snap-frozen. The HNSCC cohort used in this study included an untreated cohort of 48 human papillomavirus (HPV)-negative HNSCC and 20 age-matched and sex-matched controls. For the head and neck sample cohort, all samples were histologically confirmed with 80% tumor contact in each sample. The cohort demographics were the following: 65% men, 35% women; mean age of 68 (SD 3.2) years. Sixty-five percent of patients were stage I and II, and 35% were stage III and IV HNSCC.

For immunotherapy samples, pretreatment tumor tissue of 21 patients with HNSCC and lung cancer who received anti-PD-1 immunotherapy was obtained from Columbia Biobank. This cohort included 12 responders and 9 non-responders. The cohort demographics were 65% men, 35% women; mean age of 66 (SD 4.21) years. We used the Response Evaluation Criteria in Solid Tumors V.1.116 to retrospectively identify the treatment response to immune-checkpoint blockade. As previously described,17 clinical benefits (responders) were defined as having experienced partial response or stable disease lasting at least 6 months, whereas non-clinical benefits (non-responders) were defined as a primary progressive disease or stable disease lasting less than 6 months.

Data and bioinformatic analysis

Our analysis included the The Cancer Genome Atlas (TCGA) pan-cancer atlas dataset, which contains 10 522 tumors across 33 cancer types, including 528 HNSCC samples. The following data are available at gdc.cancer.gov/node/977: cancer type, mutation status, clinical data, HPV stats, reverse phase protein array, and mRNA gene expression. Immunogenomic studies of clinical datasets of metastatic melanoma treated with anti-PD-1 (n=38)18 and metastatic melanoma treated with CTLA-4 (n=110)19 were downloaded from cBioportal.20

Estimation of individual immune subtype fractions by xCell in TCGA samples was performed using xcell.ucsf.edu21 and TIMER.22 For the correlation analysis, Pearson’s correlation coefficient was calculated with a Bonferroni-corrected p value of ≤0.05, which was considered statistically significant. We generated survival curves of HNSCC cases in the TCGA cohort according to the expression status of the SP140 gene. A group cut-off of ‘median’ was identified, and the Kaplan-Meier curve was plotted. Immune Subtype Classifier in R was used for the classification of pan-cancer immune subtypes (five-gene signature)23 to identify the correlation of SP140 with pan-cancer immune subtypes. CIBERSORT and TIMER were used for the immune cell analysis of the immunogenomic studies.22 24

Gene set enrichment analysis (GSEA)

GSEA analysis was performed to identify pathways upregulated in SP140 high tumors compared with tumors with low expression of SP140. A total of 528 HNSCC TCGA samples were grouped into SP140 high tumors and SP140 low tumors based on the median expression value of SP140 in the entire cohort. Expression with a Bonferroni corrected p value of <0.01 was analyzed by the GSEA preranked algorithm.25 For each gene set, the normalized enrichment score, p value, and false discovery rate (FDR) q values were calculated based on the Cancer Hallmark pathway and Gene Ontology Pathways.

Chromatin immunoprecipitation (ChIP)

After being transfection with SP140 small interfering RNA (siRNA) for 48 hours, cells were fixed with 1% (v/v) formaldehyde for 7 min, and cross-linking was stopped with 0.125 M glycine. Chromatin was isolated by using the Zymo kit (D5209), according to the manufacturer’s instructions. Sheared chromatin was conjugated with ChIP grade SP140 and control IgG antibodies overnight at 4°C while rotating. ZymoMag Protein A beads (15 µL) were added to ChIP reactions and incubated for 1 hour at 4°C while rotating. Tubes were placed in a magnetic field. Protein A bead/antibody/chromatin complexes were washed with 1 mL of chromatin wash buffer, and the beads were diluted in 500 µL of elution buffer. After washing, elution, and de-cross-linking, the ChIP DNA was detected by PCR. The PCR primers used for this assay are listed in online supplemental table 1.

Supplemental material

Cell culture

Human THP-1 monocytes were cultured in Roswell Park Memorial Institute (RPMI) medium with 10% FBS, 1% penicillin, and streptomycin. Cells were incubated at 37°C in a 5% v/v CO2 atmosphere. For monocyte–macrophage differentiation, cells were seeded at a density of 2.5×105 cells/mL, and macrophage differentiation was induced by exposing the cells to 100 nM/mL phorbol-12-myristate-13-acetate (Invivogen) for 24 hours and replaced with fresh media for 24 hours. Head and neck cancer cell lines (Cal27, FaDu) were grown in complete Dulbecco's Modified Eagle Medium (DMEM) media (Gibco) supplemented with 10% heat-inactivated fetal bovine serum (FBS) (Corning), 1% penicillin, and streptomycin (Corning). SCC7, a murine head and neck carcinoma cell line, was cultured in DMEM media supplemented with 10% FBS, 2 mM L-glutamine, and 100 U/mL penicillin/streptomycin.

In vivo syngeneic tumor model and TAMs’ isolation

Eight-week C3H/HeOuJ female mice were obtained from The Jackson Laboratory (n=8). 1×105 SCC7 cells were inoculated subcutaneously and mice were sacrificed on day 21. Tumors were dissociated using a mouse tumor dissociation kit (Miltenyi Biotec) with the gentleMACS Dissociator according to the manufacturer’s protocol (Miltenyi Biotec). CD11b magnetic beads (Miltenyi) were used for isolation of TAMs after tumor dissociation as recommended by the manufacturer. TAMs were cultured in RPMI 1640 culture medium (Gibco) supplemented with 10% heat-inactivated FBS (Corning), 1% penicillin, and streptomycin (Corning) and were used for ex vivo experiments. This study was approved by the Columbia University IACUC committee (AC-AABS4603).

Flow cytometry

A flow cytometry panel consisting of human CD80 (PE, BioLegend), human CD86 (BV711, BioLegend), and human CD206 (PerCP-Cy5.5, BioLegend) was used for the characterization of macrophages. For intracellular staining of phospho-STAT1 (pSTAT1ty701), cells were prepared using the Fix and Perm kit as recommended by the manufacturer (Invitrogen). Antibody-capture beads (CompBeads, BD Biosciences) were used as single-color compensation controls. Cytometer calibration was performed daily by using rainbow fluorescent particles (BD Biosciences) after acquiring unstained and single-color control samples to calculate the compensation matrix. Data were analyzed using FCS Express software.

Cell viability assay

Cell viability was assessed using a Cell Viability Assay Kit (Fluorometric-Blue) (ab112120) following the manufacturer’s recommendation. The same amount of dye-loading solution (100 μL/well) was added to the wells. Cell were incubated at room temperature for 1 hour and the fluorescent intensity was quantified by excitation and emission of 405 and 460, respectively.

siRNA transfection

Human THP-1 and peripheral blood mononuclear cells (PBMC) naïve macrophages were transfected with SP140 siRNA for 6 hours in Opti-MEM using lipofectamine RNAimax (Thermo Fisher). The transfection procedure was followed as per manufacturer instructions. After transfection, Opti-MEM media were replaced with fresh RPMI complete media. Cells were collected for RNA and protein expression analysis at 24 and 48 hours, respectively.

RNA isolation and quantitative real time-PCR

RNA isolation from the formalin-fixed paraffin-embedded (FFPE) tissue section was performed using a quick DNA/RNA FFPE kit (Zymo Research) according to the manufacturer’s instructions. In cell culture experiments, cells were lysed with TRI reagent (Zymo Research), and RNA was extracted with Direct-zol RNA Miniprep plus (Zymo Research) according to the manufacturer’s instructions. The quantity and quality of the RNA were determined by NanoDrop 1000 (260:280 and 260:230 ratios). Typically, 500–1 μg of total RNA was reverse transcribed using iScript reverse transcription supermix (Bio-Rad) using a T100 PCR Thermal Cycler (Bio-Rad). The amplification protocol included incubations at 94°C for 15 s, 63°C for 30 s, and 72°C for 60 s. The cDNA generated was used for quantitative PCR (qPCR) on a CFX96 Real-Time thermocycler (Bio-Rad) as per the manufacturer’s instructions using iTaq Universal SYBR Green Supermix (Bio-Rad). Target gene expression was normalized to the expression of the housekeeping gene 18S. Relative gene expression was calculated using the standard 2-ΔΔCT method. The primers used for this study are presented in online supplemental table 1.

SP140 CRISPR/dCas9 lentiviral activation particle transduction

Cloning the SP140 overexpressed into macrophages (THP1-derived or TAMs isolated from syngeneic mouse tumors) was performed using SP140 CRISPR/dCas9 lentiviral activation transduction particles, a synergistic activation mediator transcription activation system for the transcriptional activation of endogenous genes, according to the manufacturer’s protocol (Santa Cruz Biotechnology). Briefly, cells were plated onto six-well plates 24 hours before transfection. The culture media was replaced with RPMI-1640 medium plus 10% fetal bovine serum (FBS) and 5 µg/mL polybrene (Thermo Fisher), then transduced with SP140 lentiviral activation transduction particles (Santa Cruz Biotechnology). Forty-eight hours after transduction, the cells were selected by the addition of 10 mg/mL puromycin dihydrochloride for at least 7 days. The mock-transfected cells were transduced with control lentiviral activation particles and went through the same selection protocol.

Assay of transposase accessible chromatin sequencing (ATAC-seq) data processing and analysis

ATAC-seq reads of macrophages treated with LPS (0 hour), macrophages treated with LPS (4 hours), SP140 KO macrophages (0 hour), and SP140 KO macrophages (4 hours) were obtained from the Gene Expression Omnibus (GSE89177).26 In addition, H3K27me3 ChIP-seq reads for macrophages treated with GM-CSF+IFN-γ (GSM1625982), along with ChIP-seq reads for baseline macrophages (GSM1625981), were obtained.26 Three sets of reads were available for each of the H3K27me3 conditions and, due to low read counts, the three reads in each condition were merged as BAM files after alignment. ATAC-seq data were aligned via Bowtie2 V.2.427 using alignment parameters—very-sensitive and −k 10 per Harvard FAS Informatics guidelines to an hg38 index provided by the Bowtie2 package (GRCh38 no-alt analysis set). H3K27me3 data were aligned using default package parameters. After coordinate sorting and indexing (and merging in the case of the H3Kme3 data) via Samtools V.1.11,28 BAM files were imported as GeneRegionTracks using Gviz V.1.34.129 in R V.4.0.2 and mapped against appropriate chromosomal regions (histogram plot, with a sliding window of an arbitrary size determined by package).

Immunofluorescence

HNSCC FFPE tissue array sections were deparaffinized and washed with xylene and graded alcohol, followed by antigen retrieval by using 10 mM citrate buffer (pH 6.0). Non-specific antibody binding was reduced by incubating the sections with 3% bovine serum albumin (BSA) for 30 min at room temperature. The sections were then incubated with anti-SP140 primary antibody (MyBioSource) at 4°C overnight. Samples were washed with phosphate-buffered saline (PBST) (Triton) three times and incubated with Alexa Fluor 488-conjugated goat anti-rabbit IgG (Invitrogen) secondary antibody (Green) and PE/Cyanine7 anti-human CD68 antibody (orange) for 2 hours under dark and subsequently washed with PBST. Slides were incubated with DAPI (1:1000) (Sigma) for 5 min at room temperature in the dark, protected from light. Tissue sections were washed with PBS three times and mounted with Prolong Diamond Antifade Mountant (Invitrogen). After 48 hours, tissue sections were digitalized under an automated bright-field and fluorescence scanner at 20× (Leica Aperio VERSA 8 Scanning System). The QuPath open-source software platform30 was used to identify and quantify TAMs (CD68+) and colocalization of SP140 with TAMs.

Cytokine analysis

The LegendPlex Human Proinflammatory Chemokine Panel and Cytokine Panels were used according to the manufacturer’s instructions. Briefly, the assay was carried out in V-bottom 96-well plates. Mixed beads (25 µL), streptavidin-phycoerythrin (PE), and detection antibodies were used. Cell culture supernatant was added to mixed beads and incubated for 2 hours; wells were washed twice with wash buffer. Samples and standards were incubated with detection antibodies for 1 hour. Streptavidin–PE was then added and the plate was incubated for 30 min. Finally, beads were washed and resuspended in 200 µL of wash buffer, and data were acquired on a four-laser BD LSR Fortessa X20. All incubation steps were carried out in the dark at room temperature, on an orbital shaker set at 600–800 RPM. Data were analyzed using the LegendPlex data analysis software as instructed by the manufacturer.

Immunoblotting

Whole-cell lysates were prepared on ice using RIPA lysis buffer solution (Thermo Fisher), 1% protease, and phosphatase inhibitor cocktail (Invitrogen). Protein concentrations were determined with the Protein BCA Assay (Thermo Fisher). Lysate protein was subjected to 10% sodium dodecyl-sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a nitrocellulose membrane. After blotting, membranes were probed with a primary antibody at 4°C overnight. The following primary antibodies were used: anti-SP140 (MyBioSource), anti-STAT5a (Invitrogen), anti-phospho-STAT5aY694 (Invitrogen), anti-STAT1 (Invitrogen), and anti-phospho-STAT1Y701 (p-STAT1Y701) (Cell Signaling Technology). After washing three times with Tris-buffered saline with 0.1% Tween 20 Detergent (TBST), membranes were incubated with horseradish peroxidase (HRP)-conjugated secondary goat anti-rabbit (Bio-Rad) or anti-mouse (Bio-Rad) antibodies as indicated for 2 hours at room temperature. Blots were developed with the Clarity Max Western ECL detection system (Bio-Rad) according to the manufacturer’s instructions, and images were captured using an iBrightCL1000 (Invitrogen).

Statistical analysis

For all experiments, the number of technical and/or biological replicates is reported in the figure legends or text. To identify mRNA expression patterns that significantly correlate with SP140 expression, a Pearson’s correlation was performed, and p values were corrected using Bonferroni adjustment. For experimental data, parametric Student’s t-test or Mann-Whitney U test was performed to compare two groups based on the underlining distribution of data. Survival curves were calculated using the Kaplan-Meier method and compared using the log-rank test. Data are demonstrated as mean±SEM. Two-sided p values less than 0.05 were considered statistically significant. In GSEA, the p values were calculated by the permutation test. Data were analyzed using R V.4.0.3 and GraphPad Prism V.8 (GraphPad Software). In all cases, ****p<0.0001, ***p<0.001, **p<0.01, and *p<0.05.

Results

SP140 expression in HNSCC clinical samples

We first analyzed the expression levels of SP140 in different types of cancers. Overall, increased expression of SP140 was seen in many cancers in the TCGA, including HNSCC (online supplemental figure 1).

In the TCGA dataset, patients with HNSCC showed overall higher expression of SP140 in tumor tissue compared with normal adjacent tissue with some heterogeneity in the expression (n=44 normal, n=520 primary tumor, p<0.0001) (figure 1A). Since mutation in TP53 is the most common molecular driver of HNSCC, we assessed the SP140 expression levels in relation to TP53 mutation status. SP140 was upregulated in both TP53 mutant (n=327) and non-mutant (n=175) patients with HNSCC, with a significant increase in the TP53 non-mutant (p<0.0001) (figure 1B). Since a subset of HNSCCs is HPV+, we next compared SP140 expression in HPV+ and HPV tumors. Although the level of SP140 was higher than in normal tissue for both HPV+ and HPV HNSCC, the level of SP140 was significantly higher in HPV+ tumors (n=80) compared with HPV tumors (n=434) (p<0.0001) (figure 1C).

Figure 1

SP140 expressed heterogeneously in HNSCCs and its high expression in TAMs is associated with a favorable survival and high tumor mutation burden. (A) SP140 expression in primary HNSCC tumors and normal adjacent tissue in TCGA data (n=520 HNSCC, n=44 normal). (B) SP140 expression versus TP53 mutation status in patient with HNSCC tumors in the TCGA dataset (n=327, TP53 mutant, and n=175, TP53 non-mutant). (C) SP140 expression in primary HPV+ HNSCC tumors, HPV− HNSCC tumors, and normal adjacent tissue in the TCGA data (n=44, normal; n=80, HPV+ HNSC, and n=434, HPV− HNSCC). (D) Representative images of SP140 expression in HNSCC TAMs (CD68+) of normal tissues and HNSCC tumor tissues. (E) Mean fluorescent intensity was calculated in TAMs (CD68+) by quantification of SP140 expression using QuPath software and compared between controls, patients with HNSCC with unfavorable 5-year survival, and patients with HNSCC with favorable 5-year survival. (F) Correlation of SP140 expression and number of non-silent mutations in 497 HNSCCs in the TCGA dataset. (G) Correlation of SP140 expression and number of non-silent mutations in over 9766 cancers in the TCGA dataset. (A–C) Box and whisker plots, with whiskers representing minimum and maximum for the panel. (E) Scatter plot (median and 95% CI). P values (A–C) were calculated by Student’s t-test and corrected for multiple comparison by Bonferroni correction. The correlation between SP140 and mutation burden was quantified with Pearson’s correlation coefficient and corrected for multiple comparison. ***P< 0.001, ****P< 0.0001 after correction for multiple comparisons. The scale bars indicate 100 μm in the left images and 30 μm in the right images. HNSCC, head and neck squamous cell carcinoma; HPV, human papillomavirus; RSEM, RNA-Seq by Expectation-Maximization; SP140, speckled protein 140; TAM, tumor-associated macrophage; TCGA, The Cancer Genome Atlas.

SP140 is described as an immune-restricted gene.26 Consistently, qPCR analysis of HNSCC lines (FaDu and Cal27) confirmed the low expression of SP140 in HNSCC cell lines and high expression in macrophages (online supplemental figure 2), consistent with SP140’s known role in immune cells. Multiplex immunofluorescence (IF) staining for SP140 showed high expression of SP140 localized in TAMs (CD68+) in HNSCC samples (figure 1D). Interestingly, levels of SP140 were higher in patients’ TAMs with favorable 5-year survival (figure 1E).

SP140 expression is positively associated with high tumor mutation burden, favorable survival outcome, and immune hot tumors

The SP140 expression level positively correlated with mutation load across cancers and specifically in HNSCC (p=9.604E-8 for pan-cancer and p=0.0170 for HNSCC after FDR correction, n=9766 for pan-cancer and n=497 for HNSCC) (figure 1F,G). High levels of SP140 were associated with an increase in lymphocyte-specific protein tyrosine kinase (LCK) protein expression in the TCGA HNSCC dataset (figure 2A). Higher expression of LCK is associated with expression of the immune response transcription module, infiltration of the tumor tissue by activated immune cells, and a favorable response to the immune blockade.31 High levels of SP140 were associated with high protein expression of proapoptotic proteins CASP7 (p<0.001) and CDKN1B (p<0.01) and low levels of PXN protein (p<0.01) in HNSCC (figure 2B–D). The high expression of PXN is associated with poor prognosis in HNSCC.32 High levels of SP140 were associated with better disease-specific survival (log-rank test p=0.01039) and overall survival (log-rank test, p=0.0205) in HNSCC (figure 2E and online supplemental figure 3).

Figure 2

High levels of SP140 in HNSCC tumors are associated with higher expression of proapoptotic markers and favorable survival in HNSCC. (A–D) HNSCC TCGA tumors with available protein expression data (n=212) were divided to two groups of SP140 high and SP140 low based on the median of SP140 expression. The protein levels of LCK, CASP7, CDKN1B, and PXN were compared between the two groups by Student’s t-test, and the p value was calculated after Benjamini-Hochberg correction. (E) Disease-specific survival of patients with high levels of SP140 (n=261) and low levels of SP140 (n=261) in TCGA HNSCC. The p value was calculated based on the log-rank test. **P<0.01, ***P<0.001 after correction for multiple comparisons. HNSCC, head and neck squamous cell carcinoma; RPPA, reverse phase protein array; SP140, speckled protein 140; TCGA, The Cancer Genome Atlas.

We next analyzed the correlation of SP140 expression levels in relation to tumor immune subtypes across 33 cancers in TCGA. This algorithm was developed based on immunogenomic analysis of more than 10 000 tumors and identified an immune response pattern impacting prognosis and classification of the TME to six immune subtypes, including wound healing, IFN-γ-dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGFβ-dominant.23 High SP140 expression was observed in the IFN-γ-dominant tumor type (p<0.0001 after Bonferroni correction) (figure 3A). Across the TCGA pan-cancer cohort, higher levels of SP140 were associated with significantly higher inflammatory macrophage infiltrates (M1) and CD8 T-cells infiltration in most TCGA cancer cohorts (figure 3B,C).

Figure 3

High levels of SP140 in the pan-cancer TCGA tumors are associated with IFN-γ-dominant tumors and high infiltration of proinflammatory TAMs (M1) and CD8 T cells. (A) Association of SP140 gene expression with immune infiltrate subtypes in the TCGA pan-cancer cohort (n=12 839). (B,C) TCGA pan-cancer gene expression data were used to dichotomize SP140 expression based on the median. M1 macrophage and T-cell infiltrations were estimated for each cancer type in SP140 high and SP140 low tumors. Data are presented as scatter plots. The statistical significance is annotated by the number of stars. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001. IFN-γ, interferon gamma; M1, proinflammatory phenotype; SP140, speckled protein 140; TAM, tumor-associated macrophage; TCGA, The Cancer Genome Atlas.

In the HNSCC cohort of TCGA samples, SP140 expression was positively associated with higher inflammatory macrophage infiltrates (M1) and CD8 T-cell infiltration (CD8 T cells: Pearson’s r=0.4522, p value <2.2 E-16; M1 macrophages, Pearson’s r=0.4623, p value <2.2 e-16) (figure 4A,B). We next investigated the expression levels of SP140 associated genes in HNSCC. High SP140 expression in HNSCC was associated with positive enrichment of IFN-γ, inflammatory response pathways, and T-cell activation based on GSEA33 (figure 4C,D, and online supplemental table 2). Expression of SP140 was positively associated with high infiltration of M1 TAMs and CD8 T cells in HNSCC in both HPV+ and HPV HNSCC (online supplemental figure 4).

Figure 4

(A,B) Correlation of SP140 expression with M1 macrophage and CD8 T-cell infiltration levels HNSCC (n=497). (C,D) Genes whose expression significantly correlated with SP140 expression were analyzed by GSEA in HNSCC samples (n = 522). The plot shows the top positively enriched pathways, hallmark inflammatory response, IFN-γ response, and T-cell activation. The Y-axis is the GSEA enrichment score. The X-axis is a list of genes ranked by differential expression correlation with SP140, with black bars representing genes in each gene set. GSEA, gene set enrichment analysis; HNSCC, head and neck squamous cell carcinoma; IFN-γ, interferon gamma; M1, proinflammatory phenotype; SP140, speckled protein 140.

SP140 inhibits STAT1 signaling in macrophages

As SP140 was shown to be IFN-γ inducible and STAT1 plays a role in the activation of IFN-γ signaling and TAMs, we designed mechanistic studies to study the dynamics of SP140-STAT1 interactions in macrophages. We knocked down SP140 using siRNAs in human THP1-derived unpolarized macrophages. We first analyzed the expression of SP140 mRNA and confirmed decreased expression in SP140 siRNA transfected cells (figure 5A). Interestingly, suppression of SP140 increased the expression of STAT1, IL-1Ra (IL-1 receptor antagonist), arginase 1 (ARG1), and IL-6 (figure 5B–D and online supplemental figure 5). However, the knockdown of SP140 did not affect STAT5a levels (online supplemental figure 5). Analysis of ATAC-seq data including macrophages treated with LPS (0 hour), macrophages treated with LPS (4 hours), SP140 KO macrophages (0 hour), and SP140 KO macrophages (4 hours) showed that SP140 acts as an epigenetic repressor at the promotor of STAT1 and controls its expression at the transcriptional level. Analysis of ChIP-seq data from baseline and IFN-γ stimulated macrophages suggests that this repression is correlated with an increase in repressive H3K27me3 modifications in the region of STAT1. (figure 5E). ChIP-qPCR analysis after silencing SP140 confirmed the interaction of SP140 with the STAT1 promotor (online supplemental figure 6). Protein analysis showed that STAT1 and pSTAT1Ty701 are significantly upregulated in SP140 knockdown cells, indicating inhibitory effects of SP140 on STAT1 activation (figure 5F). CRISPR/dCAS9 activation of SP140 confirmed the decrease in pSTAT1 Tyr701 in macrophages after overexpression of SP140 (figure 5G). Interestingly, SP140 knockdown in macrophages was associated with a significant increase in repressive H3K27me3 modifications upstream of the TSS of STAT1. These data indicate the repressor effects of SP140 on STAT1 activation. Mechanistic studies showed that an increase in STAT1 overexpression in macrophages contributes to the TAM immune suppressive activity in the TME and TAM-mediated T-cell deletion in the TME.34

Figure 5

SP140 negatively regulates STAT1 transcription and phosphorylation in macrophages. (A–D) SP140 was downregulated using siRNA in naïve (undifferentiated) macrophages, and cells were collected for RNA expression analysis by qPCR for SP140, STAT1, IL-1Ra and arginase-1, after 24 hours. 18S was used as an endogenous normalizer. The delta–delta Ct method was used to identify the relative gene expression. Data are presented as mean and SD of fold change compared with the control. (E) ATAC-seq read coverage (average count over an arbitrary sliding window) of SP140 KO and non-KO macrophages treated with LPS for 4 hours (along with their respective 0-hour controls) and average H3K27me3 ChIP-seq coverage for IFN-γ stimulated and unstimulated macrophages in the neighborhood of the STAT1 gene. (F) SP140 was downregulated using siRNA in naïve (undifferentiated) macrophages, and cells were collected for protein isolation and western blot 48 hours post transfection. (I) Effect of SP140 overexpressed on Levels of pSTAT1Tyr701 was quantified by flow cytometry in controls (shaded line indicates control macrophages; blue line indicates scramble control (CRISPR/dCAS9 control), and red line indicates test sample (SP140 overexpression, CRISPR/dCAS9)) after 48 hours. *P<0.05, **P<0.01, ***P<0.001. ATAC-seq, assay of transposase accessible chromatin sequencing; CHiP, chromatin immunoprecipitation; IFN-γ, interferon gamma; IL, interleukin; qPCR, quantitative PCR; siRNA, small interfering RNA; SP140, speckled protein 140.

SP140 overexpression induces production of inflammatory cytokines and chemokines and reprograms macrophages to induce TAM-mediated tumor cytotoxicity

Overexpression of SP140 with a CRISPR/dCAS9 activation approach led to an increase in the percentage of M1 (CD80+CD86+) macrophages (figure 6A,B), indicating an enhancement in antigen presentation capacity of macrophages. Overexpression of SP140 induced increases in IFN-γ, CXCL10, and IL-12 in unstimulated and stimulated macrophages (figure 6C–E). Consistently, high levels of SP140 expression were associated with high levels of IFN-γ, IL-12, and CXCL10 in TCGA HNSCC tumors (p<0.0001 for all comparisons) (figure 6F–I). Interestingly, high levels of CXCL10 are reported to be required for an antitumor immune response following immune-checkpoint blockade.35 An increase in CRISPR/dCAS9 activation of SP140 in TAMs isolated from a syngeneic mouse model of HNSCC led to an increase in direct cytotoxicity of TAMs on the HNSCC SCC7 cell line (figure 6J).

Figure 6

Overexpression of SP140 induces the production of inflammatory cytokines and chemokines and reprogram macrophages to induce TAM-mediated tumor cytotoxicity. (A,B) Control CRISPR/dCAS9 (control) or SP140 CRISPR/dCAS9 (SP140 overexpression) was transfected in naïve macrophages, and the expression of CD80 and CD86 was quantified by flow cytometry. (C–E) SP140 CRISPR/dCAS9 (SP140 OE) or scramble control (control CRISPR/dCAS9) were introduced to the naïve macrophages, and a multiplex fluorescence BioLegend assay was used to identify the levels of CXCL10, IFN-γ, and IL-12 (P70) in the supernatant after 48 hours. LPS (100 nM) was administered after 24 hours. (F–I) Correlation of SP140 expression and IFNG, CXCL10, IL-12B, and IL-12A levels in over 500 HNSCCs in the TCGA dataset. (J) Control CRISPR/dCAS9 or SP140 CRISPR/dCAS9 was introduced in TAMs isolated from syngeneic HNSCC tumor and cocultured with murine SCC7 cells. Cell viability after 48 hours of coculture was quantified with a quantitative viability assay kit. Data are presented based on the fold change of non-treated control. *P<0.05, ***P<0.001, ****P<0.0001. HNSCC, head and neck squamous cell carcinoma; IFN-γ, interferon gamma; IL, interleukin; SP140, speckled protein 140; TAM, tumor-associated macrophage; TCGA, The Cancer Genome Atlas.

Correlation of SP140 with immune cell infiltrates and response to anti-PD-1 immune-checkpoint blockade

We next evaluated the expression of SP140 in patients who were treated with anti-PD-1 immune-checkpoint inhibitors in HNSCC and non-small cell lung cancer (n=21). Samples were collected for gene expression analysis before treatment. A significantly higher level of SP140 was associated with response to anti-PD-1 immunotherapy, defined as stable disease or remission after anti-PD-1 treatment over 6 months (p<0.05) (figure 7A). Receiver operating curve analysis showed better performance of SP140 compared with PD-L1 in diagnostic accuracy (area under the curve for SP140: 0.77±0.10 vs area under the curve for PD-L1: 0.55±0.13 (p>0.05) (figure 7B,C). Examining the dataset of metastatic melanoma patients who underwent CTLA-4 blockade immunotherapy,19 we found significantly more favorable survival in patients who had high SP140 levels (n=12) than in patients with low levels of SP140 (n=28) (log-rank test, p=4.98E-3) (figure 7D). High levels of SP140 were associated with a significantly higher immune score (p<0.001), higher CD8+ T cells (p<0.001), higher M1 macrophages, and higher gamma delta T cells (p<0.05) (figure 7E–H). Similar results were found in metastatic melanomas which were treated with anti-PD-1 immune blockade (online supplemental figure 7). Examining the datasets of patients with recurrent glioblastoma who underwent neoadjuvant anti-PD-1 therapy36 and patients with thymic carcinoma who underwent anti-PD-1 therapy,37 we found that higher levels of SP140 were associated with a favorable response to the neoadjuvant anti-PD-1 therapy in recurrent glioblastoma and a favorable response to anti-PD-1 in thymic carcinoma (online supplemental figure 8). These data indicate the vulnerability of SP140 high tumors to immunotherapy.

Figure 7

High expression of SP140 in tumors is associated with a favorable response to immunotherapy. (A) RNA was isolated from pretreatment specimens of HNSCC and lung cancer responders and non-responders to anti-PD-1 immunotherapy. RNA was isolated and levels of SP140 were quantified by RT-qPCR (n=21). 18S was used for normalization. Student’s t-test was used for statistical analysis. (B,C) Receiver operating curves of SP140 and PD-L1 for discriminating responders or non-responder cases. Y-axis represents sensitivity (%) and X-axis represents 100% specificity (%). (D) Patients with metastatic melanoma were dichotomized to high expression and low expression groups based on expression of SP140 (the highest quartile vs the rest). Overall survival of patients with high levels of SP140 (n=12) and low levels of SP140 (n=28) was graphed and analyzed using the Kaplan-Meier estimate. (E,F) Tumors with high expression of SP140 (n=12) versus tumors with low levels of SP140 (n=28) showed higher infiltration of M1 macrophages, CD8 T cells, gamma delta T cells, and overall immune score. The Wilcoxon test was used for statistical analysis, and the p value was corrected for multiple comparisons. **P<0.01, ***P<0.001, ****P<0.0001. HNSCC, head and neck squamous cell carcinoma; M1, proinflammatory phenotype; PD-1, programmed cell death protein 1; PD-L1, programmed death-ligand 1; RT-qPCR, reverse transcription–quantitative PCR; SP140, speckled protein 140.

Discussion

In this study, our results showed that tumors with high levels of SP140 have more infiltrating M1 TAMs and an IFN-γ signature in the TME. Considering the outcome of immunotherapy using anti-PD-1 agents, patients who had a high level of SP140 in the tumor showed a favorable response to immunotherapy, defined as stable disease or remission after anti-PD-1 treatment over 6 months. SP140 mRNA conferred a higher diagnostic accuracy compared with PD-L1 for predicting response to immunotherapy in our cohort. Consistently, a higher level of SP140 in HNSCC, but not PD-L1, was associated with better survival outcomes in patients with HNSCC in the TCGA data. High SP140 level in tumors was associated with high CD8+ T cells in the tumor in both HPV+ and HPV tumors. These data indicate the possible utility of using SP140 expression for patient selection in immunotherapy.

Recently, ChIP-seq data showed that the majority of endogenous SP140 in primary human macrophages is associated with promoters.38 It has been shown that SP140 preferentially occupied heterochromatic regions with low chromatin accessibility and low transcriptional activity, marked by high levels of histone H3 lysine 27 trimethylation (H3K27me3) and low histone H3 lysine 4 trimethylation (H3K4me3), in human macrophages and plays an important role in maintaining the macrophage transcriptional landscape.26 38 Interestingly, our bioinformatic analysis showed an association of SP140 with STAT1 promoter in macrophages. Given that the siRNA-mediated knockdown of SP140 in human macrophages resulted in an increased level of STAT1, IL-6, IL-1R, and ARG1. In HNSCC cells, high levels of SP140 were associated with high protein expression of proapoptotic and tumor suppressor proteins CASP7, CDKN1B, and BIM and low levels of PXN protein. A recent pan-cancer analysis showed that abnormal PXN expression is associated with a poor prognosis.32

In this study, we found an association between high levels of SP140 and high non-synonymous mutation burden. As previously described, tumors with a high mutation burden express a large number of abnormal proteins, neo-antigens, that could elicit T-cell immune response.39 In two independent immunotherapy cohorts of HNSCC/lung cancer and metastatic melanoma, tumors with a high mutation burden showed an improved clinical response to immune-checkpoint blockade and survival.40 Given the high accuracy of SP140 to predict response to anti-PD-1 observed in this study, SP140 levels could be potentially used for the prediction of mutation burden and response to immunotherapy. Furthermore, using the TME immune profile classification developed based on immunogenomic analysis of more than 10 000 tumors,23 we found an association of a high level of SP140 with IFN-γ-dominant tumor immune signature. These data are in accordance with the higher presence of M1 macrophages and CD8+ T cells in SP140 high tumors and their better response to immunotherapy.

In the present study, we showed that SP140 suppresses STAT1 at its transcriptional level and suppresses its phosphorylation at post-transcriptional levels, indicating that inhibition of SP140 could activate STAT1. The role of STAT1 in tumor progression and tumor immunity is not straightforward and is context dependent. Although the tumor-suppressive functions of STAT1 have already been demonstrated in several studies,41–43 in a mechanistic study performed by the Gabrilovich group,34 it was determined that STAT1, but not STAT3 or STAT6, is responsible for TAM immunosuppressive activity. These results are in accordance with the antitumor activity of macrophages high in SP140 and low in STAT1 as observed in this study. Interestingly, in a recent study, it has been shown that persistent activation of STAT1 is associated with the in vivo acquisition of transcriptomic features associated with relapse after radiation and anti-CTLA4 therapy.44 They also observed that prolonged STAT1 activation on interferon treatment is associated with resistance and clinical progression after anti-PD-1 therapy.44 Our observation regarding the inhibitory role of SP140 on STAT1 might be an explanation for the better response to immunotherapy in SP140 high tumors and inhibitory effects of SP140 on STAT1 activation.

As previously reported, SP140 expression is mostly associated with unpolarized (naive macrophages) and classically activated macrophages compared with the alternatively activated macrophages. Mechanistically, overexpression of SP140 increased CD80, the T-cell costimulatory molecule on macrophages and increased IFN-γ production by macrophages and their cytotoxicity of macrophages toward HNSCC cancer cell lines. Future preclinical and clinical studies are needed to investigate the potential use of SP140 activation as a therapeutic approach for cancer immunotherapy.

The results of our study showed that high levels of SP140 were associated with the presence of a higher inflammatory response, higher IFN-γ signature in the tumors, favorable survival in HNSCC, and favorable response to immunotherapy in HNSCC/lung squamous cancer, metastatic melanoma, thymic carcinoma, and recurrent glioblastoma who received neoadjuvant anti-PD-1 therapy. SP140 was expressed in higher abundance in undifferentiated (naïve) macrophages and classically activated macrophages. SP140 plays an inhibitory role in STAT1 activation. Overexpression of SP140 induced inflammatory and antitumor macrophage activity. SP140 could potentially serve as both a disease-associated biomarker to a priori identify immunotherapy responders and could be therapeutically upregulated as an adjunct to other immunotherapy modalities to elicit enhanced antitumor responses.

Data availability statement

Data are available in a public, open access repository. Data are available in a public, open-access repository from the following sources: https://gdc.cancer.gov/node/977 ; https://www.cbioportal.org/datasets; The Gene Expression Omnibus under accession number GSE121810, GSE181815, GSE89177, GSM1625982, GSM1625981.

Ethics statements

Patient consent for publication

Ethics approval

This study involves human subjects and was approved by Columbia University Medical Center's institutional review board (IRB) (IRB number AAAS0169). We used a retrospective cohort of samples. The IRB waived the informed consent requirement. The samples were obtained based on IRB-approved protocols.

Acknowledgments

We thank the Herbert Irving Comprehensive Cancer Center Database Shared Resource and the Herbert Irving Comprehensive Cancer Center Molecular Pathology Core.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter @3kkranthi, @FatemehHeravi

  • Contributors KKKT and FM-H planned and performed the majority of experiments and associated data analysis and wrote the paper; JM-B, MC, AT, and JN provided assistance in data analysis. SD, HN, RC, AM, JN, and RC provided technical assistance and data interpretation. AY and Y-SLC provided assistance with the characterization and interpretation of patient data and tissue specimens. FM-H provided reagents and resources for the project, supervised the whole project and performed data interpretation. All authors contributed to the editing of the manuscript and critical review. FM-H acts as the guarantor.

  • Funding This work was supported by NIH/NIDCR (DE031112-01 and DE031112-01), the American Association of Cancer Research (AACR 20-60-51), and the Mark Foundation Grant, NIH/NCAT (KL2TR001874 and UL1TR001873) and NIH/NCI (P30CA013696).

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.