Background kalirin RhoGEF kinase (KALRN) is mutated in a wide range of cancers. Nevertheless, the association between KALRN mutations and the pathogenesis of cancer remains unexplored. Identification of biomarkers for cancer immunotherapy response is crucial because immunotherapies only show beneficial effects in a subset of patients with cancer.
Methods We explored the correlation between KALRN mutations and antitumor immunity in 10 cancer cohorts from The Cancer Genome Atlas program by the bioinformatics approach. Moreover, we verified the findings from the bioinformatics analysis with in vitro and in vivo experiments. We explored the correlation between KALRN mutations and immunotherapy response in five cancer cohorts receiving immune checkpoint blockade therapy.
Results Antitumor immune signatures were more enriched in KALRN-mutated than KALRN-wildtype cancers. Moreover, KALRN mutations displayed significant correlations with increased tumor mutation burden and the microsatellite instability or DNA damage repair deficiency genomic properties, which may explain the high antitumor immunity in KALRN-mutated cancers. Also, programmed cell death 1 ligand (PD-L1) expression was markedly upregulated in KALRN-mutated versus KALRN-wildtype cancers. The increased antitumor immune signatures and PD-L1 expression in KALRN-mutated cancers may favor the response to immune checkpoint blockade therapy in this cancer subtype, as evidenced in five cancer cohorts receiving antiprogrammed cell death protein 1 (PD-1)/PD-L1/cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) immunotherapy. Furthermore, the significant association between KALRN mutations and increased antitumor immunity was associated with the fact that KALRN mutations compromised the function of KALRN in targeting Rho GTPases for the regulation of DNA damage repair pathways. In vitro and in vivo experiments validated the association of KALRN deficiency with antitumor immunity and the response to immune checkpoint inhibitors.
Conclusions The KALRN mutation is a useful biomarker for predicting the response to immunotherapy in patients with cancer.
- molecular immunogene
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
Recently, immunotherapy has achieved success in treating various cancers.1–4 In particular, immune checkpoint blockade (ICB) methods, such as inhibition of cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), programmed cell death protein 1 (PD-1), and programmed cell death 1 ligand (PD-L1), have been used as therapy for many advanced malignancies.5 Nevertheless, many patients with cancer failed to respond to immunotherapy. Numerous studies have identified specific genetic or genomic features associated with cancer immunotherapy response, such as PD-L1 expression6 and tumor mutation burden (TMB).7 Mutations of certain genes, such as TP53 8 and ARID1A,9 have been associated with cancer immunotherapy response.
Kalirin RhoGEF kinase (KALRN) encodes a protein that activates specific Rho GTPase family members to regulate neurons and the actin cytoskeleton.10 This gene is mutated in a wide range of cancers,11 12 for example, melanoma, lung cancer, uterine corpus endometrial carcinoma (UCEC), glioblastoma multiforme (GBM) and colorectal cancer (COAD). Nevertheless, very few reports of the association between KALRN mutations and carcinogenesis have been published.13 14 In this study, we investigated the associations between KALRN mutations and antitumor immune signatures (NK cells, CD8 +T cells and immune cytolytic activity (ICA)) in 10 individual cancer types (CESC, COAD, ESCA, GBM, PRAD, READ, SARC, SKCM, STAD, and UCEC) and in them all together (the pan cancer) by analyzing The Cancer Genome Atlas (TCGA, https://cancergenome.nih.gov) cohorts (online supplemental table S1). For these cohorts, we also explored the associations between KALRN mutations and cancer immunotherapy response-related biomarkers, such as PD-L1 expression and TMB. Moreover, based on five cancer cohorts receiving immunotherapy (online supplemental table S1), we investigated the association between KALRN mutations and cancer immunotherapy response. We validated our findings from the bioinformatics analysis by performing in vitro experiments in four tumor cell lines (MGC803, SJSA1, SW620 and B16F10), which were associated with four of the 10 cancer types, and in vivo experiments with mouse tumor models. Our study demonstrates that KALRN mutations can promote antitumor immunity and cancer immunotherapy response and are a useful biomarker for stratifying patients with cancer responsive to immunotherapy.
Evaluation and comparisons of immune signature enrichment scores
We calculated the enrichment score for an immune signature in a tumor sample by single-sample gene-set enrichment analysis (ssGSEA)15 of the marker genes of the immune signature. Three immune signatures were analyzed, namely NK cells (marker genes KLRC1 and KLRF1), CD8 +T cells (CD8A) and ICA (PRF1 and GZMA).11 12 16 We compared the immune signature enrichment scores between two classes of samples by using the Mann-Whitney U test (one tailed). In addition, we compared the ratios of immune-stimulatory to immune-inhibitory signatures (CD8+/CD4 +regulatory T cells, M1/M2 macrophages, and proinflammatory/anti-inflammatory cytokines) between two classes of samples by using the Mann-Whitney U test (one tailed). The ratios were the log2-transformed mean expression levels of the marker genes of immune-stimulatory signatures over those of immune-inhibitory signatures. We carried out these analyzes with R programming software (ssGSEA scores were calculated by using the R package ‘GSVA’.15)
We identified the protein–protein interaction network of KALRN with the webtool STRING (V.11.0, https://string-db.org/) by using the default settings.
Logistic regression analysis
We used logistic regression with two predictors (KALRN mutation and Rho GTPase gene mutation) to predict the immune signature scores (high (upper third) vs low (bottom third) scores) and immunotherapy response (responsive vs non-responsive), respectively. Both predictors were binary variables: 1 (mutation) or 0 (no mutation). In performing the logistic regression analyzes, we used the R function ‘glm’ to fit the binary model and the function ‘lm.beta’ in the R package ‘QuantPsyc’ to calculate the standardized regression coefficients (β values).
We compared the expression levels of genes and proteins between two classes of samples by using Student’s t test (two tailed). We performed the analyzes with log2-transformed values of the normalized RNA-Seq gene expression values and the downloaded normalized protein expression data. We used Fisher’s exact test to examine the association between two categorical variables. All these statistical analyzes were performed in the R programming environment (V.4.0.2).
We performed the univariate survival analyzes by using the function ‘survfit’ in the R package ‘survival’, used Kaplan-Meier survival curves to display the survival time differences, and applied the log-rank test to evaluate the significance of the survival time differences.
In vitro experiments
Antibodies, reagents and cell lines
All antibodies were used at a dilution of 1:1000 unless otherwise specified. Anti-PD-L1 (ab213524) and anti-CD8α were purchased from Abcam (Burlingame, CA). GAPDH (sc32233) was purchased from Santa Cruz Biotechnology, California, USA. Anti-MSH2 (15 520–1-AP) was purchased from Proteintech Group (Rosemont, Illinois, USA). Anti-KALRN (BS3514) was purchased from Bioworld Technology (St. Louis, MO). PE anti-mouse TNF-α Antibody (catalog no. 506305), APC anti-mouse IFN-γ Antibody (catalog no. 505809), APC anti-mouse CD366 (Tim-3) Antibody (catalog no. 134007), APC anti-mouse CD279 (PD-1) Antibody (catalog no. 135209), APC anti-mouse CD223 (LAG-3) Antibody (catalog no. 125209) were purchased from BioLegend (San Diego, CA). Dulbecco’s modified Eagle’s medium (12 100–061), Alpha minimum essential medium (MEM) (11900024) and horse serum (16050122) were purchased from Thermo Fisher Scientific (Waltham, Massachusetts, USA). Inositol (I7508), 2-mercaptoethanol (M3148) and folic acid (F8758) were purchased from Sigma-Aldrich (St. Louis, Missouri, USA). The human cancer cell lines MGC803 (gastric cancer), SJSA1 (osteosarcoma) and SW620 (colon cancer) and the murine cell line B16F10 (melanoma) were from the American Type Culture Collection. They were cultured in a complete medium (Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum) in a humidified incubator at 37°C and 5% CO2. NK92 cells (KeyGEN BioTECH, Nanjing, China) were cultured in Alpha MEM with 2 mM L-glutamine, 1.5 g/L sodium bicarbonate, 0.2 mM inositol, 0.1 mM 2-mercaptoethanol, 0.02 mM folic acid, 100–200 U/mL recombinant human IL-2 (PeproTech, Rocky Hill, New Jersey, USA), and a final concentration of 12.5% horse serum and 12.5% fetal bovine serum.
Knockdown of KALRN with small interfering RNA (siRNA)
Cancer cells were transfected with KALRN siRNA or control siRNA by using Effectene Transfection Reagent (Qiagen, Hilden, Germany, B00118) according to the manufacturer’s instructions. The medium was replaced after 24 hours incubation with fresh medium, and the cells were maintained for a further 24 hours. The transfection efficiency was detected with quantitative PCR or Western blotting. KALRN siRNA and control siRNA were synthesized by GenePharma (Shanghai, China). Their sequences were as follows: KALRN siRNA: 1, 5'-GCUUCCACUGAAGUACCUATT-3' (sense) and 5'-UAGGUACUUCAGUGGAAGCTT-3' (antisense); 2, 5'-GCAAUCGCCCAUUGAGUAUTT-3' (sense) and 5'-AUACUCAAUGGGCGAUUGCTT-3' (antisense); 3, 5'-GCUUCGACCUUGGACACUUTT-3' (sense) and 5'-AAGUGUCCAAGGUCGAAGCTT-3' (antisense); control siRNA: 5'-UUCUCCGAACGUGUCACGUTT-3' (sense) and 5'-ACGUGACACGUUCGGAGAATT-3' (antisense).
Lentivirus generation and infection
Heteroduplexes (small hairpin RNA or shRNA) containing different 19 bp RNA segments to target different locations of the KALRN gene were selected. The heteroduplexes, supplied as 58-nucleotide oligomers, were annealed, and the downstream of the U6 promoter was inserted into the pLKO.1 plasmid to generate pLKO.1/ShKALRN. Recombinant and control lentiviruses were produced by transiently transfecting pLKO.1/vector and pLKO.1/ShKALRN, respectively. The lentiviruses were transfected into 293 T cells, along with the packaging plasmid psPAX2 and envelope plasmid pMD2.G, by using the calcium phosphate precipitation method. After 48 hours, lentiviral particles were collected and concentrated from the supernatant by ultracentrifugation. Effective lentiviral shRNA was screened by infecting these viruses with B16F10 cells, and their inhibitory effect on KALRN expression was analyzed by Western blotting. The lentivirus containing the ShKALRN RNA target sequences and a control virus were used for the animal study. The coding strand sequence of the shRNA-encoding oligonucleotides was 5'-CCGGGTGGAGTTAATGTGCCTTGTTCTCGAGAACAAGGCACATTAACTCCACTTTTTG-3' for KALRN.
Cell extracts were generated by using lysis buffer supplemented with protease inhibitor cocktail immediately before use. Proteins were denatured by addition of 6 volumes of SDS sample buffer and boiling at 100°C for 10 min and were then separated by SDS-PAGE. After electrophoresis, proteins were electrotransferred onto nitrocellulose membranes. The membranes were blocked in 5% nonfat dry milk buffer at 25°C for 1 hour and then incubated overnight at 4°C with specific antibodies. Binding of the primary antibody was detected by using peroxidase-conjugated secondary antibodies and a SuperSignal West Dura chemiluminescent substrate kit (Thermo Scientific, 34075) according to the manufacturer’s instructions. The Western blotting result was semiquantified by using ImageJ software to measure the intensities of the bands. The band densities were normalized to the background, and the relative optical density ratios were calculated relative to the housekeeping gene GAPDH.
Total RNA was extracted from cells and reversely transcribed to cDNA. Quantitative PCR was performed with the BioRad CFX96 Touch Real-Time PCR Detection System (BioRad, Richmond, California, USA) using iQ SYBR1 Green Supermix (BioRad, Richmond, California, USA). The threshold cycle numbers were obtained by using BioRad CFX manager software. The program for amplification was one cycle of 95°C for 2 min, followed by 40 cycles of 95°C for 10 s, 60°C for 30 s, and 72°C for 30 s. The relative amount of each gene was normalized to the amount of GAPDH. The primer sequences were as follows: hKALRN: 5'-GATTGTCATCTTCAGTGA-3' (forward) and 5'-AGGACCAAGTAATTCATC-3' (reverse); hGAPDH: 5'-TGTGGGCATCAATGGATTTGG-3' (forward) and 5'-ACACCATGTATTCCGG GTCAAT-3' (reverse).
Co-culture of tumor cells with NK92 cells
A transwell chamber (Corning, New York, USA) was inserted into a six-well plate to construct a co-culture system. Tumor cells (MGC803, SJSA1, SW620 and B16F10) were seeded on the six-well plate at a density of 5×104 cells/well, and NK92 cells were seeded on the membrane (polyethylene terephthalate, pore size of 0.4 µm) of the transwell chamber at a density of 5×104 cells/chamber. Tumor cells and NK92 cells were co-cultured in a humidified incubator at 37°C and 5% CO2 atmosphere for 48 hours.
EdU proliferation assay
After co-culture of the tumor cells with NK92 cells for 48 hours, the proliferation capacity of the NK92 cells was measured by an EdU (5-ethynyl-2'-deoxyuridine; Invitrogen, California, USA) proliferation assay. NK92 cells were plated in 96-well plates with a density of 2×103 cells/well for 24 hours. Before fixation, permeabilization and EdU staining, the cells were incubated with 10 µM EdU at 37°C for 24 hours. The cell nuclei were stained with 4',6-diamidino-2-phenylindole (DAPI) at a concentration of 1 µg/mL for 20 min. The proportion of NK92 cells incorporating EdU was determined with fluorescence microscopy. Each assay was performed in triplicate wells.
In vivo experiments
In vivo mouse models
B16F10 tumor cells were transduced with ShCon (scramble) or ShKALRN lentivirus and selected by puromycin for 7 days. The stably transfected B16F10 tumor cells (2×106) were subcutaneously injected into the right flank of recipient mice after shaving the injection site. After 5 days, when the tumor volume was approximately 4–5 mm3, half of the ShCon and ShKALRN mice were treated with 100 U/L PD1/PDL1 inhibitor BMS-1 (concentration 500 mg/mL; i.p.) (MCE Cat. No. HY-19991) every 3 days. The tumors were isolated from mice after 15 days. Tumor volumes did not exceed the maximum allowable size according to the LJI IACUC animal experimental protocol. The tumor volume was measured every 2 days after the tumor appeared on the fifth day and was calculated according to the formula: V=1/2 × width2×length.
Isolation of tumor-infiltrating lymphocytes (TILs)
From subcutaneously established B16F10 melanoma-bearing mice, tumors were excised and chopped with tweezers and scissors and were then digested with 2 mg/mL collagenase (type IV, Sigma, UK) for 45 min. Following digestion, tumors were filtered through 70 µM cell strainers. The cell suspension was washed twice in culture medium by centrifugation at 1500 rpm and 4°C for 10 min. After the washing, the cells were resuspended in 6 mL of PBS and were layered over 3 mL of 30%–100% gradient percoll (Beijing Solarbio Science & Technology, Beijing, China); this was followed by centrifugation at 2600 rpm for 30 min at 25°C. The enriched TILs were obtained at the interface as a thin buffy layer, were washed with PBS twice, and finally were re-suspended in FACS staining buffer for further staining procedures.
TILs were stained with CD8 and PD-1 or LAG3 (Biolegend, San Diego, California, USA) and were analyzed by flow cytometry. TILs were restimulated with cell stimulation cocktail (eBioscience, San Diego, California, USA), and expression of IFN-γ and TNF-α (Biolegend) was analyzed by flow cytometry. Staining for cell surface markers was performed by incubating cells with antibody (1:100 dilution) in FACS buffer (0.1% BSA in PBS) for 30 min at 4°C. For intracellular cytokine (IFN-γ and TNF-α) staining, surface markers were stained before fixation/permeabilization (Intracellular Fixation & Permeabilization Buffer Set, ThermoFisher).
Immunofluorescence of CD8, PD-L1 and CD49b
Paraffin-embedded mice tumor tissue (subcutaneously established B16F10 melanoma) sections (3 µm thick) were subjected to immunofluorescence with CD8 (Abcam), PD-L1 (Abcam), or CD49b (Abcam) primary antibodies. Before immunostaining, tumor tissue sections were deparaffinized with xylene, rehydrated and unmasked in sodium citrate buffer (10 mM, pH 6.0), and treated with a glycine solution (2 mg/mL) to quench autofluorescence. After blockade of non-specific antibody binding sites with a solution of 1% BSA in PBS for 1 hour at 25°C, tissue slides were incubated overnight at 4°C with anti-CD8, anti-PD-L1, or CD49b rabbit primary antibodies (1:200 dilution) in blocking solution. Irrelevant isotype-matched rabbit IgG (Sigma–Aldrich) was applied in parallel to serial sections as negative controls. The following day, the immunoreactions were revealed by using Alexa Fluor-488-conjugated goat anti-rabbit IgG (for CD8 and CD49b) or Alexa Fluor-594-conjugated goat anti-rabbit IgG (for PD-L1) (Invitrogen, San Diego, California, USA) secondary antibodies (1:200 dilution) in blocking solution. The immunolabeled slides were examined with a fluorescence microscope after nuclear counterstaining with DAPI. Green, red and blue channel fluorescence images were acquired with a Leica DFC310 FX 1.4-megapixel digital color camera equipped with LAS V.3.8 software (Leica Microsystems, Wetzlar, Germany). Overlay images were reconstructed by using the free-share ImageJ software.
KALRN mutations correlate with elevated anti-tumor immunity
We found that the immune signatures (NK cells, CD8 +T cells and ICA) consistently displayed higher enrichment levels in KALRN-mutated than KALRN-wildtype cancers in the pan cancer analysis (Mann-Whitney U test, one tailed, p<0.001) (figure 1A). Moreover, the NK cell, CD8 +T cell, and ICA enrichment scores were significantly higher in KALRN-mutated than KALRN-wildtype cancers in four, five, and five individual cancer types, respectively (Mann-Whitney U test, one tailed, p<0.05) (figure 1A). We compared the ratios of immune-stimulatory to immune-inhibitory signatures in the pan cancer and in individual cancer types. We found that the ratios (CD8+/CD4 +regulatory T cells, M1/M2 macrophages, and proinflammatory/anti-inflammatory cytokines) were significantly higher in KALRN-mutated than KALRN-wildtype cancers in the pan cancer and in multiple individual cancer types (Mann-Whitney U test, p<0.05) (figure 1B). The EdU proliferation assay showed that NK cells co-cultured with KALRN-knockdown tumor cells had a significantly stronger proliferation capacity than NK cells co-cultured with KALRN-wildtype tumor cells (figure 1C). Taken together, these data suggest that KALRN mutations enhance antitumor immunity in cancer.
KALRN mutations correlate with increased TMB and DNA damage repair deficiency
The TMB has been recognized as a predictive biomarker for cancer immunotherapy response.7 We found that KALRN-mutated cancers exhibited a significantly higher TMB (defined as the total count of gene somatic mutations) than KALRN-wildtype cancers in the pan cancer and in all 10 individual cancer types (Mann-Whitney U test, one tailed, p<0.01) (figure 2A). In addition, KALRN-mutated cancers displayed more predicted neoantigens16 than KALRN-wildtype cancers in the pan cancer and in five individual cancer types (Mann-Whitney U test, p<0.01) (figure 2A).
Among the five cancer types (COAD, ESCA, READ, STAD and UCEC) with a prevalent MSI subtype, MSI-high (MSI-H) cancers had significantly higher mutation rates of KALRN than MSI-low (MSI-L) and microsatellite stable (MSS) cancers in the pan cancer and in four individual cancer types (Fisher’s exact test, p<0.01) (figure 2B). Moreover, we found that KALRN mutations co-occurred with mutations of DNA mismatch repair pathway genes (PMS1, PMS2, MLH1, MSH2, MSH3 and MSH6) and DNA replication and repair genes (POLD1 and POLE) in the pan cancer and in multiple individual cancer types (Fisher’s exact test, p<0.05) (figure 2C). By analyzing the TCGA protein expression profiling datasets with MSH2 and MSH6 expression values available, we found that MSH2 and MSH6 were downregulated in KALRN-mutated vs KALRN-wildtype cancers in UCSC and that MSH6 was downregulated in KALRN-mutated cancers in COAD (Student’s t test, p<0.05). Furthermore, in vitro experiments showed that MSH2 expression was reduced in KALRN-knockdown relative to KALRN-wildtype tumor cells (MGC803, SJSA1, SW620 and B16F10) (figure 2D). Overall, these results suggest a significant association between KALRN mutations and deficient DNA damage repair, which could explain the significantly increased TMB in KALRN-mutated cancers.
Associations of the KALRN downstream target Rho GTPases with TMB, DNA damage repair and tumor immunity
KALRN is a member of the Rho-GEF family that activates Rho GTPases by promoting GDP/GTP exchange.17 As a result, KALRN inactivation may compromise the activation of Rho GTPases. Indeed, KALRN expression was significantly downregulated in KALRN-mutated vs KALRN-wildtype cancers in the pan cancer analysis (Student’s t test, two tailed, p=1.11 × 10-5), suggesting that KALRN mutations are inactivating mutations that impair KALRN activation. Network analysis with STRING18 showed that KALRN interacted with multiple members of the Rho GTPase family, including RAC1, RAC2, RAC3, RHOA, RHOB, RHOC, RHOD, RHOG and CDC42 (figure 3A). Moreover, we found that KALRN mutations likely co-occurred with mutations of these Rho GTPase genes in the pan cancer and in eight individual cancer types (Fisher’s exact test, p<0.05) (figure 3B). Previous studies demonstrated that Rho GTPases might play a role in regulating the DNA damage response.19 As expected, Rho GTPase gene-mutated cancers had a significantly higher TMB than Rho GTPase gene-wildtype cancers in the pan cancer and in eight individual cancer types (Mann-Whitney U test, one tailed, p<0.01) (figure 3C). Rho GTPase gene-mutated cancers exhibited more predicted neoantigens16 than Rho GTPase gene-wildtype cancers in the pan cancer and in five individual cancer types (Mann-Whitney U test, p<0.05) (figure 3C). The mutations of Rho GTPase genes co-occurred with the mutations of DNA damage repair pathway genes (PMS1, PMS2, MLH1, MSH2, MSH3, MSH6, POLD1 and POLE) in the pan cancer and in multiple individual cancer types (Fisher’s exact test, p<0.05) (figure 3D). Moreover, the mutations of Rho GTPase genes more frequently occurred in MSI-H than in MSI-L/MSS cancers in the pan cancer and in two individual cancer types (STAD and COAD; Fisher’s exact test, p<0.05) (figure 3E). Collectively, these results confirmed the critical role of Rho GTPases in mediating DNA damage repair pathways.
Because Rho GTPase deficiency may cause DNA damage repair deficiency that, in turn, would result in increased TMB and neoantigens, we expected that Rho GTPase gene-mutated cancers would display enhanced antitumor immune signatures. As expected, the three immune signatures (NK cells, CD8 +T cells and ICA) were consistently stronger in Rho GTPase gene-mutated than Rho GTPase gene-wildtype cancers in the pan cancer (Mann-Whitney U test, one tailed, p<0.01), as well as in diverse individual cancer types (Mann-Whitney U test, one tailed, p<0.05) (figure 3F). Moreover, Rho GTPase gene-mutated cancers had significantly higher ratios of immune-stimulatory signatures to immune-inhibitory signatures (CD8+/CD4 +regulatory T cells and M1/M2 macrophages) in the pan cancer and in diverse individual cancer types (Mann-Whitney U test, one tailed, p≤0.05) (figure 3F).
Together, these data suggest that the significant correlations of KALRN mutations with increased TMB/neoantigens, deficient DNA damage repair pathways, and enhanced antitumor immune infiltrates could be consequences of the deregulation of the Rho GTPases targeted by KALRN.
KALRN mutations correlate with elevated PD-L1 expression and favorable immunotherapy response
We found that PD-L1 showed significantly higher expression levels in KALRN-mutated than KALRN-wildtype cancers in the pan cancer and in four individual cancer types (Student’s t test, p<0.05) (figure 4A). In vitro experiments showed that PD-L1 expression was significantly upregulated in KALRN-knockdown versus KALRN-wildtype tumor cells (MGC803, SJSA1, SW620 and B16F10) (figure 4B). These data suggest that PD-1/PD-L1-directed immunotherapy is likely to be more effective against KALRN-mutated cancers because PD-L1 expression is a biomarker for an active response to anti-PD-1/PD-L1 immunotherapy.6 To prove this hypothesis, we investigated the association between KALRN mutations and ICB therapy response in five cancer cohorts receiving anti-PD-1/PD-L1/CTLA-4 immunotherapy, namely the Allen cohort1 (melanoma), Hellmann cohort20 (lung cancer), Rizvi cohort21 (lung cancer), Riaz cohort22 (melanoma) and Hugo cohort23 (melanoma). We found that KALRN-mutated cancers had significantly higher response rates than KALRN-wildtype cancers in these cohorts (37.04% vs 10.96% in the Allen cohort, 45% vs 11.76% in the Hugo cohort, 80% vs 27.5% in the Riaz cohort, 100% vs 40.74% in the Rizvi cohort, and 70% vs 29.31% in the Hellmann cohort; Fisher’s exact test, p<0.1) (figure 4C). Because of the appreciably higher response rates, KALRN-mutated cancers displayed a more favorable overall or progress-free survival tendency than KALRN-mutated cancers (log-rank test, p=0.262, 0.067, 0.055, 0.176 and 0.047 for the Allen, Hugo, Riaz, Rizvi and Hellmann cohorts, respectively) (figure 4D). However, the correlation between KALRN mutations and overall survival was much less significant in TCGA melanoma or lung cancer cohorts not receiving immunotherapy (figure 4E). These results indicate that the more favorable ICB therapy response in KALRN-mutated cancer cohorts is attributed to the high antitumor immune infiltrates and PD-L1 expression in this subtype, suggesting that KALRN mutation is a predictive biomarker for an active response to ICB therapy.
Validation by in vivo experiments
To verify whether KALRN deficiency can enhance anti-tumor immune signatures and sensitization to ICB therapy, we performed in vivo experiments with mouse tumor models. We found that KALRN loss markedly increased B16F10-tumor volume and progression (figure 5A). KALRN-depleted tumors displayed significant increases in CD8 +T cell infiltration (figure 5B) and PD-L1 expression (figure 5C); this is consistent with the previous results. Moreover, in KALRN-depleted tumors, the expression levels of IFN-γ and TNF-α were significantly reduced (figure 5D), indicating that the activity of tumor-infiltrating CD8 +lymphocytes is impaired in KALRN-depleted tumors. The expression levels of PD-1 and LAG3 were also markedly enhanced in KALRN-deficient tumors (figure 5E,F), suggesting that KALRN deficiency promotes CD8 +T cell exhaustion. Likewise, KALRN-depleted tumors displayed significant increases in NK cell infiltration (figure 5G) but reduced activity of tumor-infiltrating NK lymphocytes (figure 5H). Moreover, the volume of the KALRN-depleted tumors markedly decreased relative to that of the KALRN-wildtype tumors after use of BMS-1, a small-molecule inhibitor of PD-1/PD-L1 (figure 5I); this result indicates higher sensitivity of the KALRN-depleted tumor to the PD-1/PD-L1 inhibitor. This is consistent with the previous finding that KALRN mutations correlated positively with the response to ICB therapy. Furthermore, although the activities of tumor-infiltrating CD8 +and NK lymphocytes were lower in KALRN-depleted than KALRN-wildtype tumors, they were significantly higher in KALRN-depleted tumors after treatment with BMS-1 (figure 5J,K). These data suggest that the PD-1/PD-L1 inhibitor inhibits the exhaustion of tumor-infiltrating CD8+ and NK lymphocytes in KALRN-depleted tumors, thereby promoting the killing of tumor cells by the immune system.
For the first time, we report the significant role of KALRN mutations in promoting anti-tumor immunity and immunotherapeutic vulnerabilities in diverse cancer types. KALRN mutations can enhance antitumor immunity because they increase the tumor mutation load to create more neoantigens as a result of DNA damage repair deficiency. KALRN-mutated cancers express more PD-L1, in addition to the increased anti-tumor immune infiltrates. Because both tumor-infiltrating lymphocyte density and PD-L1 expression are important factors in determining ICB therapy responsiveness,6 24 it would be reasonable for KALRN-mutated cancers to exhibit increased susceptibility to ICB therapy. This inference was evidenced in five cancer cohorts receiving ICB therapy (figure 4C,D).
We investigated the underlying mechanism of KALRN mutations in enhancing antitumor immunity and immunotherapy response (figure 6). The significant association between KALRN mutations and heightened antitumor immunity may be attributed to the fact that KALRN inactivation compromises the function of the target Rho GTPases in regulating DNA damage repair pathways. Previous studies have revealed the function of Rho GTPases in modulating DNA damage response.24–26 In this study, we demonstrated a significant correlation between Rho GTPases and DNA damage repair pathways (figure 3C–E). Rho GTPases are involved in a wide variety of cellular functions and oncogenic processes, including the cell cycle, cell death, cell migration, cell polarity and cell adhesion.27 Notably, previous studies have shown that inhibition of Rho GTPases can promote an antitumor immune response by increasing the ICA of cytotoxic CD8 +T lymphocytes.27 This finding is in line with our results showing that the antitumor immune signatures grew stronger in cancers harboring mutations of Rho GTPase genes (figure 3F). Interestingly, we found that the cancers with mutations of Rho GTPase genes overexpressed PD-L1 in the pan cancer and multiple individual cancer types (Mann-Whitney U test, one tailed, p<0.05). This suggests that mutations of Rho GTPase genes may favor the ICB therapy response. Indeed, we found that cancers harboring mutations of Rho GTPase genes showed higher ICB therapy response rates than those without such mutations in the five cancer cohorts receiving ICB therapy (42.86% vs 22.78% in the Allen cohort, 70% vs 48.15% in the Hugo cohort, 83.33% vs 25.64% in the Riaz cohort, 100% vs 45% in the Rizvi cohort, and 50% vs 28.26% in the Hellmann cohort). Nevertheless, logistic regression analyzes showed that the KALRN mutation was a significant positive predictor for the immunotherapy response in three of the five cancer cohorts after correction for the Rho GTPase gene mutation (p<0.05) (online supplemental figure S1A). In addition, the KALRN mutation was a significant positive predictor for the immune signature scores in TCGA pan cancer after correcting for the Rho GTPase gene mutation (p<0.001) (online supplemental figure S1B). These results indicate that KALRN mutations may promote antitumor immunity and immunotherapy response via other approaches besides the upstream regulation of Rho GTPases. As expected, the Rho GTPase gene mutation was also a positive predictor for the immunotherapy response and immune signature scores, whereas its β coefficient was smaller than that of the KALRN mutation in almost all cases (online supplemental figure S1). In addition, we found that the association between the Rho GTPase gene mutation and survival was less significant in the KALRN-wildtype cancer subset than in all cancers in three of the cancer cohorts receiving ICB therapy (online supplemental figure S2). Overall, these results indicate that Rho GTPase mutations may not correlate independently with antitumor immunity and immunotherapy response but rather in a KALRN-dependent manner.
A limitation of this study is that the sample sizes were not sufficiently large in the cancer cohorts receiving immunotherapy. To overcome this limitation, we need to validate our findings by using more significant numbers of cancer samples with the immunotherapy response data available in the future.
The KALRN mutation correlates with enhanced antitumor immunity and immunotherapy response and, thus, is a predictive biomarker for the response to immunotherapy.
We thank Ms. Yin He for help in performing the experiments.
ML and YM contributed equally.
Contributors ML performed experiments and data analyses, and edited the manuscript. YM performed experiments and data analyses, and edited the manuscript. YZ performed experiments. QL performed experiments. CC performed experiments. LQ performed experiments. XW conceived this study, designed the methods and wrote the manuscript. All the authors read and approved the final manuscript.
Funding This work was supported by the China Pharmaceutical University (grant number 3150120001 to XW) and Jiangsu Province Graduate Research and Innovation Program Project (grant number KYCX20_0658 to ML).
Competing interests None declared.
Patient consent for publication Not required.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available in a public, open access repository. Data are available on reasonable request. All data relevant to the study are included in the article or uploaded as online supplemental information. The datasets for the 10 TCGA cancer cohorts can be downloaded from the genomic data commons data portal (https://portal.gdc.cancer.gov/). The 10 cancer cohorts included cervical squamous-cell carcinoma and endocervical adeno-carcinoma (CESC), colon adeno-carcinom (COAD), Esophageal Carcinoma (ESCA), Glioblastoma Multiforme (GBM), Prostate Adenocarcinoma (PRAD), Rectum Adenocarcinoma (READ), Sarcoma (SARC), Skincutaneous Melanoma (SKCM), Stomach Adenocarcinoma (STAD), and Uterine Corpus Endometrial Carcinoma (UCEC). For each cancer type, we obtained its gene somatic mutations, RNA-Seq gene expression profiles, protein expression profiles, and clinical data. For pan cancer, we downloaded the related datasets from UCSC (https://xenabrowser.net/datapages/). We obtained the data for predicted neoantigens in the TCGA cancer cohorts from the publication , in which a total of 6 cancer types (CESC, PRAD, GBM, STAD, SKCM, UCEC) had predicted neoantigen data available. For the four cancer cohorts receiving anti-PD-1/PD-L1/CTLA-4 immunotherapy (Allen cohort, Hugo cohort, Riaz cohort, Rizvi cohort, and Hellmann cohort), we obtained their genomics and clinical data from the associated publications. A summary of these datasets is shown in online supplemental table S1.
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.
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.