Abstract
Background/Aim: The most frequently altered epigenetic modifier in head and neck squamous carcinoma (HNSC) is the histone methyltransferase KMT2D. KMT2D catalyzes methylation of histone H3K4 resulting in open chromatin and the activation of target genes. Tumor-associated macrophages (TAMs) promote cancer growth by causing T lymphocyte exhaustion. C-C motif chemokine ligand 2 (CCL2) is a potent TAM chemotactic factor. In HNSC, TAMs have been associated with unfavorable patient outcomes and metastasis. The aim of this study was to determine the role of KMT2D in HNSC using genetically engineered in vivo models. Materials and Methods: KMT2D protein expression was correlated with lymph node metastasis in human HNSC using immunohistochemistry. Genetically engineered KMT2D and CCL2 knockout models of HNSC were created in vivo. HNSC was characterized using qRT-PCR, histopathology, and immunohistochemistry/immunofluorescence microscopy. We also analyzed the effects of KMT2D expression on the proliferation and migration of human HNSC lines. The regulation of the CCL2 gene by KMT2D was characterized using chromatin immunoprecipitation-sequencing assay of transposase accessible chromatin-sequencing, and chromatin conformation capture-sequencing. Results: Human HNSC cases with high KMT2D expression exhibited significantly increased lymph node metastasis. Reduced KMT2D expression in our genetically engineered model correlated with reduced lymph node metastasis, longer latency, and slow tumor growth. CCL2 expression was decreased in KMT2D deficient HNSC, which correlated with a reduced TAM gene expression signature. Genomic experiments demonstrated that KMT2D directly targeted the CCL2 gene. A new genetically engineered in vivo model of CCL2-null HNSC was created, recapitulating the KMT2D deficient phenotype and showing a decreased T lymphocyte exhaustion signature. Conclusion: KMT2D regulates CCL2-mediated immune response and metastasis in HNSC.
- T lymphocyte
- C-C motif chemokine ligand 2
- lymphocyte activation gene 3
- T cell immunoglobulin mucin receptor 3
- genomics
Head and neck squamous carcinoma (HNSC) is the one of the most frequent cancers worldwide, with 890,000 cases diagnosed and 450,000 deaths each year (1). The risk factors for HNSC are tobacco and alcohol consumption. A subset of HNSC, primarily of the oropharynx, has been associated with human papillomavirus infection (2). Most patients with HNSC present with advanced tumors and lymph node metastasis. Recurrence after treatment is a frequent clinical problem, and the five-year survival rate for HNSC patients is low. Treatment for advanced tumors consists of concurrent chemoradiation followed by surgery. Epidermal growth factor receptor and immune checkpoint antibodies have been approved for treatment of recurrent HNSC (3).
Genomic studies have revealed that HNSC is a heterogeneous disease at both the molecular and cellular levels (4, 5). The most frequently mutated epigenetic modifier in HNSC is the lysine methyltransferase KMT2D (6-9). KMT2D protein catalyzes methylation of histone H3K4 resulting in open chromatin and activation of target genes (10-12). Fourteen percent of HNSC cases exhibit gain of KMT2D copy number (4, 13-15), and high KMT2D expression predicted tumor recurrence in HNSC (16). Despite the frequency of altered KMT2D expression in HNSC, the mechanism by which this enzyme regulates progression of these cancers is unclear.
Tumor-associated macrophages (TAMs) promote tumor growth by fostering immunosuppression via T lymphocyte exhaustion (17). Unfavorable patient outcomes and metastasis have been associated with increased TAM (18). Macrophages have been phenotypically classified by gene expression signatures which are regulated by different microenvironmental signals (19). Anti-tumor macrophage polarization is stimulated by pro-inflammatory cytokines, while immunosuppressive pro-tumorigenic TAMs are promoted by anti-inflammatory chemokines (20, 21). High levels of TAMs predicted poor immune response (22). Increased TAMs are frequently observed in HNSC and create an immunosuppressive pro-tumorigenic environment (17). C-C motif chemokine ligand 2 (CCL2) is a potent macrophage chemotactic factor (23). Tumors expressing low CCL2 levels exhibit fewer TAMs and improved immune response (24).
Our present study demonstrated that human HNSC cases with high KMT2D expression exhibit significantly increased lymph node metastasis. We created a novel genetically engineered in vivo model of conditional KMT2D inactivation in HNSC. CCL2 expression was decreased in KMT2D deficient HNSC, which correlated with reduced TAM signature and lymph node metastasis. A new genetically engineered in vivo model of CCL2 null HNSC recapitulated the KMT2D deficient phenotype and decreased T lymphocyte exhaustion signature. We concluded that KMT2D regulates immune response and lymph node metastasis in HNSC.
Materials and Methods
Mouse procedures. Mouse procedures were approved by the institutional animal care committee. To characterize the effects of KMT2D deficiency on lymph node metastasis, we created a genetically engineered preclinical model with reduced KMT2D expression in stratified epithelia. Crossing B6N.Cg-Tg(KRT14-cre)1Amc/J with Kmt2dtm1.1Kaig/J mice (The Jackson Laboratories, Bar Harbor, ME, USA) resulted in Cre mediated deletion of KMT2D exons 16-19 in the stratified epithelia of offspring, producing loss of heterozygosity in mucosal epithelium. B6.129S4-Ccl2tm1Rol/J null mutant mice were purchased from The Jackson Laboratories. HNSC induction and tissue processing were performed according to our published protocol (25).
Cell culture and transduction. Human HNSC cell lines SCC9, SCC15, SCC25, SCC71, SCC90, and SCC104 were cultured according to our published protocol (25). Cultures tested negative for mycoplasma contamination using a commercially available kit (MilliporeSigma, St. Louis, MO, USA). SCC15, SCC25, and SCC9 cells were incubated with KMT2D CRISPR-Cas9 or control lentiviruses for 6 h in serum free medium according to manufacturer’s protocol (Horizon Discovery, Lafayette, CO, USA) at 37°C. Puromycin resistant clones were selected and expanded for analysis using western blot, proliferation, and migration experiments. In proliferation experiments, 3×104 cells from control and KMT2D CRISPR-Cas9 clones were plated in triplicate and counted daily for four days. In migration experiments, 2×105 from control and KMT2D CRISPR-Cas9 clones were plated in triplicate in transwell migration chambers for 16 h. Cells which migrated through the membranes were fixed in 70% ethanol, stained with hematoxylin, photographed under light microscopy, and quantitated using digital image analysis (Akoya Biosciences, Marlborough, MA, USA).
qRT-PCR. RNA extraction and reverse transcription were performed according to manufacturer’s instructions (Invitrogen, Carlsbad, CA, USA). cDNA was amplified using mouse primers: KMT2D 5′-AGAACAATGGACGGCCTGAG-3′ and 5′-TCAAGGTATGGG GCCGTTTA-3′, ITGAM 5′-AGCTTGGCTTTTTCAAGCGG-3′ and 5′-CTTACCCACAAGGGTCCCAC-3′, CD33 5′-GGAAA GGACCATCCAGCTCA-3′ and 5′-CAGATGCAGTTGGAGAT AGGCA-3′, CD14 5′-GAAGCAGATCTGGGGCAGTT-3′ and 5′-ACGTTGCGGAGGTTCAAGAT-3′, LY6C 5′-ACCCTTCTCTGA GGATGGACA-3′ and 5′-GATCCCTGATTGGCACACCA-3′, CD80 5′-AGTTTCTCTTTTTCAGGTTGTGAA-3′ and 5′-ACATGATGG GGAAAGCCAGG-3′, CCL2 5′-GACCCCAAGAAGGAATG GGT-3′ and 5′-ACCTTAGGGCAGATGCAGTT-3′, LAG3 5′-AATCCTTCGGGTTACCTGGC-3′ and 5′-CGTACACTGTCGCT CCAAGA-3′, TIM3 5′-GCAGTCAGGATTCGCTCTGA-3′ and 5′-TGGGAGATCCGATCCCTGAA-3′, β-actin 5′-AAAAGCCACCC CCACTCCTAAG-3′ and 5′-TCAAGTCAGTGTACAGG CCAGC-3′. qRT-PCR was performed according to our published protocol (25).
Western blot. Protein extracts from human HNSC lines were separated in 4-20% gradient SDS-PAGE, blotted with Kmt2d or β-actin antibodies, and analyzed according to our published protocol (25).
Histopathology, immunofluorescence, and immunohistochemistry. Human HNSC tissue microarrays were purchased from Tissue Array (Rockville, MD, USA). Formalin fixed mouse and human tumor tissues were processed and analyzed by histopathology, immunofluorescence, and immunohistochemistry according to our published protocol (25).
Cell death analysis. Fixed tumor sections were processed and analyzed for programmed cell death according to manufacturer’s protocol (Roche Applied Sciences, Indianapolis, IN, USA).
Fluorescence activated cell sorting. HNSC were sorted using phycoerythrin conjugated CD14 antibody and analyzed according to our published protocol (25).
Chromatin immunoprecipitation-sequencing. SCC15 cells were formaldehyde cross-linked using 1% formaldehyde. Cells were washed with phosphate buffered saline and resuspended in lysis buffer. The nuclei were collected, resuspended in the immunoprecipitation buffer and sonicated to obtain fragments of 200-300 bp. Chromatin concentrations were determined using spectrophotometry. Immunoprecipitation was performed with 10 μl KMT2D antibody and incubated overnight at 4°C. Protein G beads were added and incubated at 4°C for 2 h. Beads were washed three times and cross-links were reversed overnight at 65°C. Libraries for ChIP-seq were prepared and AmpureXP beads (Beckman Coulter, Indianapolis, IN, USA) were used for size selection. Libraries were quantified and sequenced using Illumina NovaSeq 6000 (San Diego, CA, USA). Reads were de-multiplexed and analyzed for quality control using bcl2fastq and FastQC. FASTQ reads were mapped to the reference genome using Bowtie2. Following conversion using SAMtools, BAM files were indexed and PCR duplicates were excluded using the rmdup function. Occupancy maps were generated using DeepTools and the matrix was used to generate profile plots. BAM files were converted to BigWig format. ComputeMatrix was used to quantify occupancy of reads followed by plotting using plotProfile and plotHeatmap.
Assay of transposase accessible chromatin-sequencing. SCC15 cells were lysed and nuclei were incubated with 12.5 μl of 2X transposase buffer, 2 μl of Tn5 transposase (Illumina, San Diego, CA, USA), and 10.5 μl distilled water at 37°C for 1 h followed by purification. Ten cycles of PCR were performed using Illumina dual index primers followed by purification. Following size selection, library quality and fragment size distribution were analyzed using TapeStation (Agilent, Santa Clara, CA, USA). Libraries were quantified and sequencing was performed on the Illumina NovaSeq 6000. After trimming of adapter sequences, reads were mapped using Bowtie2. Picard was used to remove PCR duplicated and low-quality reads, and MACS2 was used for peak calling.
Chromatin conformation capture-sequencing. One million SCC15 cells were formaldehyde fixed followed by lysis and DpnII digestion (New England BioLabs, Ipswich, MA, USA). PCR amplification was performed, and libraries were sequenced on a NovaSeq 6000 instrument (Illumina). Data were processed using Juicer software and Juicer tools was used to create hic-files. Juicebox was used to view count maps with K-R normalization. For Hi-C map generation, read-pairs of quality ≥30 were used.
Antibodies. Antibodies against p16 (AB241543), H3K9me3 (AB8898) and PCNA (AB92552) were purchased from Abcam (Waltam, MA, USA); The antibody against β-actin (4970) was purchased from Cell Signaling (Danvers, MA, USA); The antibody against KMT2D (HPA035977) was obtained from MilliporeSigma; The antibody against CCL2 (NBP1-07035) was procured from Novus Biologicals (Centennial, CO, USA).
Statistical analysis. Student t-test was used for two groups of parametric data. Fisher exact test was used for categorical data. p-Values <0.05 were considered statistically significant. GraphPad Prism (GraphPad, Boston, MA, USA) was used for the analyses.
Results
Given that KMT2D gain of copy number was detected in HNSC genomic data, we analyzed KMT2D protein expression in human HNSC lines and pathology specimens. KMT2D protein expression was detected in an immortalized mucosal cell line (HOK16B) and at 2-3 fold higher levels in 6 HNSC lines using western blot (Figure 1A). Thirty-one percent (20/64) of human HNSC pathology specimens with low KMT2D expression demonstrated lymph node metastasis (Figure 1B and D). In contrast, 56% (18/32) of human HNSC pathology specimens with high KMT2D expression exhibited lymph node metastasis (Figure 1C and E; p<0.03; Fisher exact test). KMT2D expression did not significantly correlate with patient sex, age, or tumor grade. Our results demonstrate that high KMT2D expression correlates with lymph node metastasis in human HNSC.
KMT2D expression correlates with lymph node metastasis in HNSC. (A) Western blot showing KMT2D protein expression in human HNSC cell lines. β-actin expression was used as gel loading control. Molecular mass is shown in kD. Immunohistochemistry showing KMT2D protein expression in human archival HNSC pathology specimens. Representative HNSC cases with low (B,D) and high (C,E) KMT2D protein expression are shown. The fraction of cases with lymph node metastasis in KMT2Dlow and KMT2Dhigh HNSC is shown. p-Value is shown.
To characterize the effects of reduced KMT2D expression on lymph node metastasis, we created a genetically engineered preclinical model with reduced KMT2D expression in stratified epithelia. Crossing B6N.Cg-Tg(KRT14-cre)1Amc/J with Kmt2dtm1.1Kaig/J mice resulted in Cre mediated deletion of KMT2D exons 16-19 in the stratified epithelia of offspring, resulting in loss of heterozygosity in mucosal epithelium. qRT-PCR results showing reduced KMT2D mRNA expression in K14Cre;KMT2D+/f mucosa is shown in Figure 2A. The homozygous null mutant mouse was not viable. We used our published carcinogenesis model to generate 20 K14Cre;KMT2D+/+ and 20 K14Cre;KMT2D+/f HNSC (25). Primary tumor latency was increased in K14Cre;KMT2D+/f compared to K14Cre;KMT2D+/+ HNSC (mean 158 vs. 131 days; p<0.04; Figure 2B; t-test). K14Cre;KMT2D+/f primary tumors exhibited slower growth than K14Cre;KMT2D+/+ HNSC (mean 185 mm3 vs. 240 mm3; p<0.05; Figure 2C). The histopathology of primary tumors from both genotypes are shown in Figure 2D and F. Lymph node metastases were significantly reduced in K14Cre;KMT2D+/f cancers compared to K14Cre;KMT2D+/+ HNSC (16% vs. 62% positive nodes; p<0.0001; Figure 2E, G and H). Our results indicate that reduced KMT2D expression correlated with reduced tumor growth and lymph node metastasis in HNSC.
KMT2D deficient HNSC exhibit reduced lymph node metastasis. (A) qRT-PCR showing KMT2D mRNA expression in K14Cre;KMT2D+/+ and K14Cre;KMT2D+/f mucosa. (B) Increased tumor latency in K14Cre;KMT2D+/f compared to K14Cre;KMT2D+/+ HNSC. Percent tumor free mice is shown on the y axis and days latency is shown on the x axis. (C) Reduced tumor growth in K14Cre;KMT2D+/f compared to K14Cre;KMT2D+/+ HNSC. Tumor volume is shown on the y axis and weeks tumor growth is shown on the x axis. Error bars indicate SEM. H&E staining showing the histopathology of K14Cre;KMT2D+/+ primary (D) and metastatic (E) tumors. H&E staining showing the histopathology of K14Cre;KMT2D+/f primary (F) and metastatic (G) tumors. Scale bar=10 μm. Tumor cells are annotated by T, stromal cells are indicated by S, lymphocytes are shown by L, and differentiated tumor cells are shown by diff. (H) Quantitation of lymph node metastases in K14Cre;KMT2D+/+ and K14Cre;KMT2D+/f HNSC. p-Values are shown.
To determine the mechanism for longer latency and reduced growth of KMT2D+/f HNSC, we examined expression of the senescence markers p16 and H3K9me3 in K14Cre;KMT2D+/+ and K14Cre;KMT2D+/f tumor tissue sections. The p16 positive cell fraction was low in both groups (0.1% vs. 0.2%; Figure 3A and B). The H3K9me3 positive cell fraction also was low in both groups (0.3% vs. 0.1%; Figure 3C and D). We then determined the apoptotic cell fraction in K14Cre;KTM2D+/+ and K14Cre;KMT2D+/f HNSC using TUNEL analysis. The apoptotic cell fraction was low in both groups (0.2% vs. 0.3%; Figure 3E and F). Statistical analysis using t-test indicated that there were no statistically significant differences in the number of senescent or apoptotic cells between K14Cre;KMT2D+/+ and K14Cre;KMT2D+/f HNSC. Finally we examined the proliferating cell fraction in K14Cre;KTM2D+/+ and K14Cre;KMT2D+/f HNSC using PCNA immunohistochemistry. The proliferating cell fraction was significantly reduced in K14Cre;KMT2D+/f compared to K14Cre;KMT2D+/+ HNSC (27% vs. 64%; p<0.01, t-test; Figure 3G-I). These results indicate that reduced cellular proliferation correlated with decreased metastasis in K14Cre;KMT2D+/f HNSC.
Reduced proliferating cells in K14Cre;KMT2D+/f HNSC. Immunohistochemistry showing expression of the senescence marker p16 in K14Cre;KMT2D+/+ (A) and K14Cre;KMT2D+/f (B) HNSC. Nuclei were counterstained with hematoxylin. Scale bar=10 μm. Immunofluorescence microscopy showing the expression of the senescence marker H3K9me3 in K14Cre;KMT2D+/+ (C) and K14Cre;KMT2D+/f (D) HNSC. Nuclei were counterstained with DAPI. Programmed cell death was determined using TUNEL analysis in K14Cre;KMT2D+/+ (E) and K14Cre;KMT2D+/f (F) HNSC. Proliferation was determined by PCNA expression in K14Cre;KMT2D+/+ (G) and K14Cre;KMT2D+/f (H) HNSC using immunohistochemistry. (I) Quantitation of PCNA+ cells in K14Cre;KMT2D+/+ and K14Cre;KMT2D+/f HNSC. Error bars indicate SEM.
To determine if reduced KMT2D expression directly regulated proliferation and migration of HNSC cells, we transduced SCC15, SCC25, and SCC9 cells with control and KMT2D CRISPR-Cas9 lentiviruses. Reduced KMT2D protein expression in puromycin resistant CRISPR-Cas9 clones is shown in Figure 4A. Control and KMT2D CRISPR-Cas9 clones were then used in proliferation and migration experiments in vitro. Cultured SCC15, SCC25, and SCC9 puromycin resistant control and KMT2D CRISPR-Cas9 clones are shown in Figure 4B, C, F, G, J, and K. Migrated SCC15, SCC25, and SCC9 puromycin resistant control and KMT2D CRISPR-Cas9 clones are shown in Figure 4D, E, H, I, L, and M. Reduced KMT2D expression did not significantly affect proliferation or migration of SCC15 cells (Figure 4N and O; t-test).
Reduced KMT2D expression does not inhibit HNSC cell line proliferation or migration. (A) Western blot indicating KMT2D expression in human HNSC lines SCC15, SCC25, and SCC9 in puromycin resistant (puro) control and CRISPR-Cas9 (KO) clones. β-actin expression was used as the gel loading control. Phase contrast microscopy images showing cultured SCC15, SCC25, and SCC9 puromycin resistant (puro) control (B, F, J) and CRISPR-Cas9 (KO) clones (C, G, K). Hematoxylin staining and brightfield microscopy images showing migrating SCC15, SCC25, and SCC9 puromycin resistant (puro) control (D, H, L) and CRISPR-Cas9 (KO) clones (E, I, M). (N) Proliferation time course experiments of SCC15, SCC25, and SCC9 puromycin resistant (puro) control and CRISPR-Cas9 (KO) clones. (O) Quantitation of migrating SCC15, SCC25, and SCC9 puromycin resistant (puro) control and CRISPR-Cas9 (KO) clones. Error bars indicate SEM.
Given that reduced KMT2D expression did not have cell autonomous effects on HNSC proliferation or migration, we investigated the tumor microenvironment in our genetically engineered model of KMT2D deficient HNSC. TAM gene expression signature ITGAM (−2.1 fold), CD33 (−2.8 fold), CD14 (−4 fold), LY6C (−4.3 fold), and CD80 (−4.6 fold) was reduced in K14Cre;KMT2D+/f compared to K14Cre; KMT2D+/+ HNSC (p<0.002; t-test; Figure 5A). We sorted CD14+ TAMs from K14Cre;KMT2D+/+ and K14Cre; KMT2D+/f HNSC using flow cytometry. The CD14+ TAM fraction was significantly reduced in K14Cre;KMT2D+/f HNSC (0.07% vs. 0.2%; p<0.04; t-test; Figure 5B and C) compared to K14Cre;KMT2D+/+ cancers. CCL2, secreted by cancer cells, is a major chemotactic factor for TAMs. CCL2 mRNA (−4.6 fold) and protein expression was significantly reduced in K14Cre;KMT2D+/f compared to K14Cre; KMT2D+/+ HNSC (p<0.002; t-test; Figure 5A, D and E). To determine if KMT2D directly targeted the CCL2 gene, we performed ChIP-seq, ATAC-seq, and HiC genomic studies. Given the limited number of mouse HNSC cells available for these experiments, we used the human HNSC line SCC15 for genomic studies. The human CCL2 gene is located at chr17:34255274. The CCL2 promoter region is located in a micro-compartment domain in SCC15 cells as shown by HiC chromatin conformation capture (Figure 5F), which has been proposed to regulate gene expression (26). KMT2D bound directly to the CCL2 promoter as shown using ChIP-seq which overlapped with open chromatin as shown using ATAC-seq (Figure 5F). These results indicate that KMT2D directly targets the CCL2 gene; KMT2D deficiency in HNSC resulted in decreased CCL2 expression, TAM signature, and cell fraction.
Reduced TAM gene expression signature in K14Cre;KMT2D+/f HNSC. (A) qRT-PCR results indicating relative ITGAM, CD33, CD14, LY6C, CD80, and CCL2 expression in KMT2D+/+ and KMT2D+/f HNSC. Error bars indicate SEM. Fluorescence activated cell sorting of CD14+ TAM from K14Cre;KMT2D+/+ (B) and K14Cre;KMT2D+/f (C) HNSC. CD14 log and forward scatter linear scales are shown. Immunohistochemistry showing CCL2 expression in K14Cre;KMT2D+/+ (D) and K14Cre;KMT2D+/f (E) HNSC. Nuclei were counterstained with hematoxylin. Scale bar=10 μm. (F) The CCL2 gene is a direct target of KMT2D. The CCL2 gene is located in a micro-compartment domain as shown by Hi-C analysis. ChIP-seq results showing KMT2D binding to the CCL2 gene which correlated with open chromatin conformation as shown by ATAC-seq. Exon-intron structure of the CCL2 gene is shown.
To determine if reduced CCL2 expression could recapitulate the phenotype of K14Cre;KMT2D+/f cancers, we induced HNSC in 20 CCL2−/− and 20 CCL2+/+ mice. qRT-PCR results indicated reduced CCL2 mRNA expression in CCL2−/− mucosa (Figure 6A). There were no significant differences in primary tumor latency in CCL2−/− compared to CCL2+/+ HNSC (mean 128 vs. 122 days; t-test; Figure 6B). CCL2−/− primary tumors exhibited slower growth than CCL2+/+ HNSC (mean 43 mm3 vs. 225 mm3; t-test; p<0.0006; Figure 6C). The histopathology of primary tumors from both genotypes are shown in Figure 6D and F. Lymph node metastases were significantly reduced in CCL2−/− cancers compared to CCL2+/+ HNSC (19% vs. 56% positive nodes; t-test; p<0.004; Figure 6E, G and H). Our results indicate that reduced CCL2 expression inhibited tumor growth and lymph node metastasis in HNSC.
Reduced growth and metastasis of CCL2−/− HNSC. (A) qRT-PCR results showing CCL2 mRNA expression in CCL2+/+ and CCL2−/− mucosa. (B) Tumor latency in CCL2−/− compared to CCL2+/+ HNSC. Percent tumor free mice is shown on the y axis and days latency is shown on the x axis. (C) Reduced tumor growth in CCL2−/− compared to CCL2+/+ HNSC. Tumor volume is shown on the y axis and weeks tumor growth is shown on the x axis. Error bars indicate SEM. H&E staining showing the histopathology of CCL2+/+ primary (D) and metastatic (E) HNSC. H&E staining showing the histopathology of CCL2−/− primary (F) and metastatic (G) HNSC. Scale bar = 10 μm. (H) Quantitation of lymph node metastases in CCL2+/+ and CCL2−/− HNSC. p-Values are shown.
To determine the mechanism for the reduced growth and metastasis of CCL2−/− HNSC, we examined expression of the senescence markers p16 and H3K9me3 in CCL2+/+ and CCL2−/− tumor tissue sections. The p16 positive cell fraction was low in both groups (0.2% vs. 0.2%; Figure 7A and B). The H3K9me3 positive cell fraction also was low in both groups (0.2% vs. 0.3%; Figure 7C and D). We then determined the apoptotic cell fraction in CCL2+/+ and CCL2−/− HNSC using TUNEL analysis. The apoptotic cell fraction was low in both groups (0.3% vs. 0.1%; Figure 7E and F). There were no statistically significant differences in the number of senescent or apoptotic cells between CCL2+/+ and CCL2−/− HNSC. Finally, we examined the proliferating cell fraction in CCL2+/+ and CCL2+/− HNSC using PCNA immunohistochemistry. The proliferating cell fraction was significantly reduced in CCL2−/− compared to CCL2+/+ HNSC (18% vs. 77%; p<0.003; t-test; Figure 7G-I). CCL2−/− HNSC expressed reduced T lymphocyte exhaustion signature (LAG3:−2.6-fold, TIM3: −3.4-fold; p<0.009; t-test; Figure 7J), which correlated with reduced TAM signature. These results indicate that reduced CCL2 expression was associated reduced T lymphocyte exhaustion (improved immune response), decreased proliferating tumor cells, and decreased metastatic HNSC.
Reduced cellular proliferation in CCL2−/− HNSC. Immunofluorescence microscopy showing the expression of the senescence marker p16 in CCL2+/+ (A) and CCL2−/− (B) HNSC. Nuclei were counterstained with DAPI. Scale bar=10 μm. Immunofluorescence microscopy showing the expression of the senescence marker H3K9me3 in CCL2+/+ (C) and CCL2−/− (D) HNSC. Programmed cell death was determined using TUNEL analysis in CCL2+/+ (E) and CCL2−/− (F) HNSC. Proliferating cell fraction was determined by PCNA expression CCL2+/+ (G) and CCL2−/− (H) HNSC using immunohistochemistry. Nuclei were counterstained with hematoxylin. (I) Quantitation of PCNA+ cells in CCL2+/+ and CCL2−/− HNSC. (J) Reduced cytotoxic T cell exhaustion gene expression signature in CCL2−/− HNSC. qRT-PCR showing the relative LAG3 and TIM3 mRNA expression in CCL2−/− compared to CCL2+/+ HNSC. Error bars indicate SEM.
Discussion
Our results indicated that KMT2D directly targeted the CCL2 gene in HNSC. CCL2 is a potent chemotactic factor for TAMs, and KMT2D deficiency correlated with reduced TAM gene expression signature in our HNSC model. TAMs have been associated with poor clinical outcome in HNSC (27). CCL2 has been associated with TAM polarization in cancers (28). HNSC with high CCL2 expression correlated with decreased overall patient survival (29). HNSC radiotherapy has been associated with increased CCL2 secretion and induction of an immunosuppressive tumor microenvironment (30). CCL2 was associated with cisplatin sensitivity in head and neck cancer cells (31). CCL2 secretion by pancreatic cancer cells recruited TAMs, which promoted aggressive tumor phenotype (32). CCL2 expression was increased in fibroblasts exposed to tobacco carcinogens (33). Furthermore, it has been shown that uptake of HNSC extracellular vesicles by monocytes resulted in increased CCL2 expression (34).
The immunosuppressive TAM phenotype inhibits cytotoxic T lymphocytes (35). In vivo TAM depletion was shown to reduce cytotoxic T lymphocyte exhaustion gene expression signatures in mouse models of melanoma and mammary tumors (36). Our study demonstrated that CCL2 null HNSC exhibited reduced cytotoxic T cell exhaustion gene expression signature. Exhausted cytotoxic T lymphocytes express plasma membrane bound inhibitory receptors including lymphocyte activation gene 3 (LAG3) and T cell immunoglobulin mucin receptor 3 (TIM3; 37) as shown in our current study. TAMs also produce IL10 and TGFβ which convert CD4+ T lymphocytes to regulatory T cells, which suppress cytotoxic T lymphocyte function (38).
Our study determined that CCL2 null mutation recapitulated reduced proliferation, tumor growth, and metastasis identified in KMT2D deficient HNSC. Reduced KMT2D expression inhibited proliferation and tumor formation by gastric cancer cells (39). Low KMT2D expression correlated with fibrosis in gastric cancer (40). KMT2D was mutated in sinonasal squamous cell carcinomas derived from inverted papillomas (41). KMT2D somatic and germline polymorphisms may also be associated with prostate cancer (42). A small molecule inhibitor of CCR2, the receptor for CCL2, correlated with increased pancreatic patient survival (43). Our future studies will determine how KMT2D regulates TAM polarization and HNSC phenotype using single cell genomics.
Conclusion
KMT2D regulates CCL2 mediated immune response and metastasis in HNSC.
Acknowledgements
The Authors thank Drs. Ke Ma, Zarema Arbieva, Alvaro Hernandez, and Mark Maienschein-Cline (University of Illinois Research Resources Center and Carver Biotechnology Center) for assistance with confocal microscopy, genomics, and bioinformatics.
Footnotes
Authors’ Contributions
J.W., C.C., and A.M.L. performed experiments, analyzed data, and wrote the manuscript. D.L.C. conceived the study, analyzed data, and edited the manuscript. Part of this study was the master’s thesis project of one of the authors (C.C.) and can be found at: indigo.uic.edu.
Conflicts of Interest
The Authors declare no conflicts of interest in relation to this study.
Funding
This study was supported by the University of Illinois Cancer Center.
- Received May 20, 2024.
- Revision received June 20, 2024.
- Accepted June 24, 2024.
- Copyright © 2024 The Author(s). Published by the International Institute of Anticancer Research.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY-NC-ND) 4.0 international license (https://creativecommons.org/licenses/by-nc-nd/4.0).













