# The immune contexture of primary central nervous system diffuse large B cell lymphoma associates with patient survival and specific cell signaling

Melissa Alame, Emmanuel Cornillot, Valère Cacheux, Valérie Rigau, Valérie Costes-Martineau, Vanessa Lacheretz-Szablewski, Jacques Colinge

► To cite this version:

Melissa Alame, Emmanuel Cornillot, Valère Cacheux, Valérie Rigau, Valérie Costes-Martineau, et al. The immune contexture of primary central nervous system diffuse large B cell lymphoma as-sociates with patient survival and specific cell signaling. Theranostics, 2021, 11 (8), pp.3565-3579. 10.7150/thno.54343. hal-03672867

HAL Id: hal-03672867

https://hal.science/hal-03672867

Submitted on 19 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Research Paper

# The immune contexture of primary central nervous system diffuse large B cell lymphoma associates with patient survival and specific cell signaling

Melissa Alame<sup>1,2,3,4</sup>, Emmanuel Cornillot<sup>1,3,4*</sup>, Valère Cacheux<sup>2,5*</sup>, Valérie Rigau<sup>5,6</sup>, Valérie Costes-Martineau<sup>5,6</sup>, Vanessa Lacheretz-Szablewski<sup>5,6✉</sup>, Jacques Colinge<sup>1,4,5✉</sup>

1. Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, Parc Euromédecine, 208 rue des Apothicaires, 34298 Montpellier, France.

2. Biological Hematology Department, Montpellier University Hospital, Saint Eloi Hospital, 34275 Montpellier, France.

3. Université de Montpellier, Faculté de Pharmacie, 15 avenue Charles Flahault, 34093 Montpellier, France.

4. Institut régional du Cancer Montpellier (ICM), Parc Euromédecine, 208 rue des Apothicaires, 34298 Montpellier, France.

5. Université de Montpellier, Faculté de Médecine, 2 rue école de Médecine, 34060 Montpellier, France.

6. Département de pathologie et d'oncologie biologique, Montpellier University Hospital, Gui De Chauliac Hospital, 34000 Montpellier, France.

Received: 2020.10.11; Accepted: 2020.11.19; Published: 2021.01.22

## Abstract

**Rationale:** Primary central nervous system diffuse large B-cell lymphoma (PCNSL) is a rare and aggressive entity that resides in an immune-privileged site. The tumor microenvironment (TME) and the disruption of the immune surveillance influence lymphoma pathogenesis and immunotherapy resistance. Despite growing knowledge on heterogeneous therapeutic responses, no comprehensive description of the PCNSL TME is available. We hence investigated the immune subtypes of PCNSL and their association with molecular signaling and survival.

**Methods:** Analysis of PCNSL transcriptomes (sequencing, n = 20; microarrays, n = 34). Integrated correlation analysis and signaling pathway topology enabled us to infer intercellular interactions. Immunohistopathology and digital imaging were used to validate bioinformatic results.

**Results:** Transcriptomics revealed three immune subtypes: immune-rich, poor, and intermediate. The immune-rich subtype was associated to better survival and characterized by hyper-activation of STAT3 signaling and inflammatory signaling, e.g., IFNγ and TNF-α, resembling the hot subtype described in primary testicular lymphoma and solid cancer. WNT/β-catenin, HIPPO, and NOTCH signaling were hyper-activated in the immune-poor subtype. HLA down-modulation was clearly associated with a low or intermediate immune infiltration and the absence of T-cell activation. Moreover, HLA class I down-regulation was also correlated with worse survival with implications on immune-intermediate PCNSL that frequently feature reduced HLA expression. A ligand-receptor intercellular network revealed high expression of two immune checkpoints, i.e., CTLA-4/CD86 and TIM-3/LAGLS9. TIM-3 and galectin-9 proteins were clearly upregulated in PCNSL.

**Conclusion:** Altogether, our study reveals that patient stratification according to immune subtypes, HLA status, and immune checkpoint molecule quantification should be considered prior to immune checkpoint inhibitor therapy. Moreover, TIM-3 protein should be considered an axis for future therapeutic development.

Key words: PCNSL, immune contexture, bioinformatics, immunotherapy, cellular interactions

## Introduction

Primary central nervous system diffuse large B-cell lymphoma (PCNSL) is a rare and aggressive extra-nodal non-Hodgkin's lymphoma (NHL) confined to the brain, spinal cord, leptomeninges, or

eyes. Newly diagnosed PCNSLs account for 4% of brain tumors and for 4-6% of extra-nodal lymphomas in immunocompetent patients. Incidence is still increasing in elderly people [1]. The emergence of

<footnote>*Equal contributions to this work.</footnote>

<footnote>✉ Corresponding authors: E-mail: jacques.colinge@inserm.fr; v-szablewski@chu-montpellier.fr.</footnote>

<footnote>© The author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.</footnote>
high-dose methotrexate (HD-MTX) and rituximab-based regimens in PCNSL therapy has drastically improved patient survival. However, only 50% of patients respond and 10-15% demonstrate a primary refractory disease [2], highlighting an unmet need for alternative therapeutic options.

The contribution of the tumor microenvironment (TME) to tumor aggressiveness, progression, and therapy resistance has been recognized in most tumors. By targeting the immune component, immunotherapies have revolutionized the treatment of cancers. In particular, antibodies directed against immune checkpoints (ICs) or ligands thereof, e.g., PD-1/PD-L1 or CTLA-4, chimeric antigen receptor T-cells (CAR-T), and T-cell engager antibodies have demonstrated clinical benefit in B-cell malignancies [3-5]. Most notably, anti-PD-1 antibodies have been approved for therapy of relapsed/refractory classical Hodgkin's lymphoma (HL) harboring increased PD-1+ tumor-infiltrating lymphocytes (TILs), high PD-L1 expression, and 9p21.1 (PD-L1 gene locus) copy number alteration [6]. In this context, the exclusive brain localization of PCNSL has aroused interest for thorough TME studies.

Significant efforts have been devoted to the identification of common genetic alterations and activating oncogenic signaling in PCNSL [7-14]. These alterations mainly involve nuclear factor-kappa B (NF-κB), B cell receptor (BCR), Toll-like receptors (TLR), Mitogen-activated protein kinase (MAPK) signaling, the DNA damage response, apoptosis, and cell cycle control. Moreover, genomic studies have suggested TME dysfunctions, such as immune evasion mechanisms, e.g., HLA (6p21) or B2M (15q21.2) copy loss, PD-L1 (9p24.1) copy number gain, and immune communication impairments, e.g., IL17REL (22q13.33) copy gain or IFNGR1 (6p23.3-q24) copy loss. High, intermediate, and low tumor mutational burden (TMB) were respectively identified in 19%, 71.5%, and 9.5% of patients with PCNSL in a 42-patient cohort [12], thus suggesting that IC inhibitor (ICI) therapy could be efficient against PCNSL [15]. However, the correlation between high PD-L1 protein expression and TME was negligible in this study.

Multiple studies have started to unravel the main cues of the PCNSL TME, although a comprehensive description is still lacking. The presence of PD-1+ TILs and the expression of PD-L1 in either microglial cells/macrophages (tumor-associated macrophages (TAMs)) or malignant B cells were found correlated with patient outcome [12,13,16]. Some potential biomarkers for ICI therapy have been identified in PCNSL, including intermediate to high TMB, 9p21.1 copy number

alteration, or PD-L1 expression [12,13,16]. Nevertheless, PCNSL has so far been mainly studied as an undifferentiated entity [17] despite a manifest heterogeneity in therapeutic responses. The existence of distinct molecular subtypes has not been addressed

We thus decided that a comprehensive study, linking a complete TME description with oncogenic and immune signaling pathways, could help our understanding of therapy resistance and potentially uncover new therapeutic opportunities. Accordingly, we functionally characterized the immune subtypes of 54 PCNSL patient samples by combining transcriptomic data analysis with histopathology and digital imaging. We showed there are three immune subtypes of PCNSL: immune-rich, poor and intermediate. We also examined the immune evasion mechanisms and the main molecular pathways associated with these different immune subtypes, highlighting new potential therapeutic opportunities, including anti-TIM-3, and the overall clinical relevance of PCNSL immune subtype classification.

## Methods

### Patients and cohorts

Surgically resected tumors from immunocompetent patients with PCNSL were retrospectively retrieved from the Department of Pathology, Montpellier University Hospital (France). Snap-frozen tissues were obtained during surgery from 20 patients. The newly diagnosed patients with PCNSL provided written informed consent for tissue collection and subsequent research purposes. Patients with prior or concurrent low grade B-cell lymphoma and central nervous system (CNS) metastasis of diffuse large B-cell lymphoma (DLBCL) were excluded. This project was approved by the research ethics boards of our institution (Centre des Ressources Biologiques, CRB, Montpellier) according to the Declaration of Helsinki (AC-2010-1200 and AC-2013-2033). Most of our patients (18/20) received Rituximab and Methotrexate based therapy. The other two patients received corticotherapy and progressed rapidly.

### Transcriptomic data

RNA was extracted from the 20 fresh tissue samples and quality-controlled. DNA libraries were prepared with the NEBNext Ultra II RNA-Seq kit and sequenced on a NextSeq 500 (Illumina) system using 75bp single reads. With our pipeline, Fastq files were aligned against the human genome (Ensembl GRCh38) (STAR using default parameters and 2-pass mode, read count extraction with HTSeq-count). Data are accessible from GEO (GSE155398).
To increase the size of our PCNSL cohort, we retrieved 34 additional PCNSL processed microarrays from Gene Expression Omnibus (GEO) (GSE34771). The authors of this dataset employed MAS 5.0 (Affymetrix Inc.) for background adjustment and log-transformation, and they performed quantile normalization [18]. The subsequent analysis and clustering of normalized RNA-sequencing and microarray data was performed using the same algorithms although separately to avoid impossible to correct batch effects between different technologies.

As PCNSL tumor cells are of the same origin than DLBCL tumor cells, 48 DLBCL transcriptomes (mRNA-sequencing) were also retrieved from The Cancer Genome Atlas (TCGA) and used as a point of comparison for PCNSL TME composition.

## Differential gene expression analyses

Differential gene expression was analyzed with edgeR [19], imposing p-value < 0.01, FDR < 0.01, minimum fold-change of 2, and a minimum average of 20 (normalized) read counts for all the samples. Heatmaps were generated with ComplexHeatmap [20]. Dendrograms for the clustering of samples (columns) and genes (rows) were constructed with Ward's method based on Euclidean distance. For color assignment, a cutoff was applied to the 2.5% highest and lowest values.

## Pathway enrichment analysis

In order to perform whole-transcriptomic and HLA status-related analysis, we implemented hypergeometric testing on differentially expressed genes to search for enriched GO terms and Reactome pathways respectively. This was followed by Benjamini-Hochberg multiple hypothesis correction with an imposed FDR < 0.05, and a minimum of 5 query genes in a Reactome pathway or 3 query genes in a GO term.

Gene signatures of described signaling pathways linked to either immune response, e.g. gamma interferon (IFNγ) and Tumor Necrosis Factor (TNFα), or oncogenesis, e.g. Protein P53 (P53) and MAPK, were retrieved from the hallmark gene sets in the Molecular Signatures Database (MSigDB v7.1) [21]. Gene Set Enrichment Analysis (GSEA) was performed with the Fast GSEA Bioconductor library [21,22]. Gene signatures listed in this study were only considered as significantly enriched (adjusted p-value < 0.05) when scored together with all the hallmark signatures.

## Ligand-receptor interactions and pathways

We recently defined a ligand-receptor (L-R) database [23] and a method to infer L-R interactions from bulk RNA expression [24]. This database contains L-R pairs, which have been validated in vitro

and/or in vivo. The methodology originally developed for mRNA-sequencing data was adapted to microarray data by using the maximum Spearman correlation between the different probes of a given ligand and a given receptor. L-R interactions were predicted separately in the two different datasets (sequencing and microarray). In total, 673 and 393 L-R pairs were determined from microarray and RNA-sequencing data respectively. Intersection identified 165 unique high-confident L-R pairs. Brain tissue interactions were discarded, considered as non-specific interactions, and in all 128 L-R pairs were conserved for further analysis. Subsequently, for these a score named L-R score was computed as previously described [24].

## Histological and immunohistochemical analysis

All cases were reviewed by expert pathologists (VLS, VR, and VCM). The diagnosis of PCNSL was made on Hematoxylin-Eosin (HE) tissue staining and was based on the WHO 2016 classification of hematopoietic and lymphoid tissue [17]. For immunohistochemical examination, 3µm-thick tissue sections from formalin-fixed paraffin-embedded (FFPE) blocks were subjected to antigen retrieval and immunostained on a Ventana Benchmark XT autostainer (Ventana Tucson, AZ, USA). The antibodies were used after appropriate antigen retrieval according to the manufacturer's instructions (Table S1). Association with Epstein-Barr virus (EBV) was examined by in situ hybridization (ISH) using Epstein-Barr encoding region (EBER). MUM1, MYC, and P53 protein expression were considered as positive if nuclear staining was observed in at least 30%, 40%, and 10% of neoplastic cells respectively. Besides PCNSL samples, we could access one postmortem normal brain sample for comparison. This sample was collected for autopsy purposes and processed identically by pathologists.

## Digital imaging

DAB-positive stained cells (galectin-9 and TIM-3 staining) were automatically counted using the open-source software QuPath [25]. Given both TIM-3 and galectin-9 stained tumor cells, macrophages, and endothelial cells - cell types with identifiable cell morphologies - their quantification was evaluated by counting DAB-positive pixels to improve accuracy. Intensity thresholds for pixel detection and classification were manually set for each staining type and performed identically for the all samples. For further analyses, pixel densities were estimated as the percentage of positive pixels per mm² of surface area [25]. All steps were performed under the supervision
of an expert pathologist (VLS). Necrosis, tissue folds, and entrapped normal structures were carefully removed.

## Interphase fluorescence in situ hybridization

Interphase fluorescence in situ hybridization (FISH) was performed on 3-µm thick tissue sections using split signal FISH DNA probes for BCL2/18q21 (probe Y5407; DAKO A/S), BCL6/3q27 (probe Y5408; DAKO A/S), MYC/8q24 (probe Y5410; DAKO A/S), and PDL1/9p24.1 (PDL1, CD274 Break Apart Probe; Empire Genomics). Digital images were captured with a Metafer Slide Scanning Platform using a Leica Axioplan fluorescence microscope (Zeiss Axio Imager M1) equipped with a charge-coupled device (CCD) camera coupled to and driven by ISIS software (MetaSystem, FISH Imaging System, Germany). A total of 100 nuclei were evaluated independently by three specialists (VLS, MA, and VC). Cases were considered positive when more than 15% of the cells exhibited abnormalities in the tissue sections.

## Statistical and survival analysis

Statistical analysis was performed using the R library survminer package. Estimation of overall survival (OS) and relapse-free survival (RFS) were generated using Kaplan-Meier method. Time-to-event distributions were compared by means of a log-rank test. Cox regression model and multivariate analysis could not be performed due to the small size of the cohort. The dependence between clinical variables and the immune subtypes was assessed by the χ² test.

# Results

## Bulk transcriptomic analysis highlights immune cell signatures in PCNSL

We conducted bulk transcriptomic analysis of our PCNSL cohort (n = 20) by mRNA-sequencing. Unsupervised clustering of the complete transcriptomes revealed a group of genes almost exclusively involved in immune activation (Figure S1, Table S2). The expression of these genes separated tumors into three distinct groups with high, intermediate, and low immune gene expression levels (Figure 1A). This initial observation indicated the need to examine the immune infiltrate of each tumor. In order to investigate this, we increased the size of our cohort by retrieving 34 additional PCNSL microarrays from GEO [18]. In addition, we obtained 48 DLBCL transcriptomes (mRNA-sequencing) from The Cancer Genome Atlas (TCGA). DLBCL provided a comparison of PCNSL with its nodal counterpart.

We investigated 13 previously reported immune cell signatures [26] to characterize both the PCNSL and DLBCL immune infiltrates. To cover the stromal

part of the TME we included signatures for cancer-associated fibroblasts (CAFs) and endothelial cells (ECs) [27]. Signature scores were computed as gene z-score averages. Samples clustering were performed using T cells because of their well-known effector role in tumor rejection.

Based on the number of T cells, both mRNA-sequencing and microarray data clustered into three distinct groups (Figure 1B-C). Considering both data types, cluster 1 (10 patients with PCNSL and one patient with DLBCL) was devoid of activated CD4+ and CD8+ T cells, but contained heterogeneous amounts of Th17 cells and macrophages. Cluster 2 (20 patients with PCNSL and 25 patients with DLBCL) was enriched in a higher number of lymphoid cells, e.g., activated CD4+ and CD8+ T cells, regulatory T cells (Tregs), and myeloid cells such as TAMs, myeloid-derived suppressive cells (MDSCs), or activated dendritic cells (DCs). The last cluster, cluster 3 (24 patients with PCNSL and 22 patients with DLBCL), assumed no specific pattern with a variable immune cell presence and at a much lower level than cluster 2.

Activated CD8+ T cell, macrophage, and Treg scores were higher in cluster 2 than in clusters 1 and 3 merged (Figure 1D). Moreover, the PCNSL classification obtained by TME cell gene signatures (Figure 1B-C) was consistent with the groups defined in Figure 1A. Hence, cluster 1 was termed immune poor, cluster 2 immune rich, and cluster 3 immune intermediate. These terms defined three immune subtypes of PCNSL. Fifty-two percent (25/48) of DLBCL tumors harbored an immune-rich phenotype versus only 37% (20/54) of PCNSL tumors. In contrast, 19% (10/54) of PCNSL tumors were devoid of CD4+ and CD8+ T cells (immune poor), whereas only one DLBCL sample (2%, 1/48) fell in this group. High proportions of both PCNSL (44%, 24/54) and DLBCL (46%, 22/48) tumors adopted an intermediate-immune phenotype. Compared to normal brain tissues, immune-poor PCNSL immune infiltrate harbored comparable amounts of cells with slight compositional variations. Immune-rich and intermediate PCNSLs included much more activated T cells as expected (Figure S2).

PCNSL tumors lacking an immune-rich infiltrate relapsed significantly earlier (Figure 1E). We found a similar trend with overall survival (OS); PCNSL patient tumors lacking an immune-rich infiltrate had a generally decreased OS (Figure 1E). We did not find any association between the immune pattern and the clinical, cytogenetic, or immunohistological variables, with exception of patient outcome (Table 1). Note that we compared the immune-rich subtype to the other two subtypes according to OS and relapse-free
survival (RFS) for our cohort only since the available clinical data for the PCNSL microarray cohort [18] were too limited (Table 2).

HLA down-regulation correlates with the immune microenvironment in PCNSL

Presentation of neoantigens via HLA molecules on the surface of malignant B-cells should induce an antitumor immune response. However, B lymphoma

cells can evade this response through various mechanisms, including loss or aberrant expression of HLA molecules [28]. Loss of either HLA class I or II expression has been demonstrated in PCNSL [28–31]. Nevertheless, these results remained incomplete regarding the link between the immune infiltrate and loss of HLA expression in PCNSL [28,30,32]. We thus investigated the global gene expression patterns related to HLA class I and II expression loss and

<figure><img src="image_3.png" /><figcaption>**Figure 1.** PCNSL with an immune-rich TME defines a patient subgroup with a better outcome. **A.** Three PCNSL clusters exist based on immune gene expression, denoted *high*, *intermediate*, and *low*. **B.** TME cell gene signatures across PCSNL (n = 20) and DLBCL (n = 48) tumor sample transcriptomes reveal three groups of tumors: the *immune-rich*, *immune-poor* and *immune-intermediate* subtypes (clusters 2, 1, 3 respectively). PCNSL were classified in agreement with the gene clusters of panel A. **C.** Applying identical TME cell gene signatures to PCNSL microarray data (n = 34) reveals four clusters. Merging the two small rightmost clusters yielded three groups of tumors (denoted 1-3) with immunological subtypes comparable to those of panel B. **D.** Quantification of activated CD8+ T cells, regulatory T cells, macrophages, and CAFs in the immune-rich tumors versus the other subtype tumors (Wilcoxon two-sided tests, n = 20 = 8+12). **E.** The immune-rich subtype features a more favorable outcome (Kaplan-Meier curve, log-rank test, n = 20 = 8+12). **Activ CD4:** activated CD4+ T cells, **Activ CD8:** activated CD8+ T cells, **Tregs:** regulatory T cells, **Tfh:** T follicular helper cells, **Th1:** type I helper cells, **Th2:** type 2 helper cells, **Th17:** type 17 helper cells, **Activ DC:** activated dendritic cells, **CD56 bright:** CD56 bright NK cells, **CD56 dim:** CD56 dim NK cells, **MDSC:** Myeloid-derived dendritic cells, **Mono:** monocytes, **TAM:** tumor-associated macrophage, **EC:** endothelial cells, **CAF:** cancer-associated fibroblast.</figcaption></figure>
potential associations with the immune subtypes identified in this study. Sixteen *HLA* genes were retrieved from the literature [33] out of which 15 were found in both PCNSL cohorts in this study, i.e., in mRNA-sequencing and microarrays. For each gene, we defined low and high expression according to its median, which was performed in mRNA-sequencing and microarray data independently. Results from both cohorts were pooled and hierarchical clustering of the total 54 samples identified four well-defined clusters (Figure 2A). Seventy-eight percent (42/54) of PCNSL tumor samples featured low HLA class I or II molecule gene expression (clusters 2, 3 and 4, Figure 2A-B), and 31% (17/54) exhibited low expressions of both (cluster 4). Twenty-two percent (12/54) featured only low HLA class I molecule gene expression (cluster 2), and lastly, 22% (12/54) had only low HLA class II molecule gene expression (cluster 3, Figure 2A-B). Out of the 22% (13/54) of tumors featuring HLA expression above median expression (cluster 1), 77% (10/13) were highly infiltrated by immune cells (immune-rich subtype), whereas none of the PNCSL tumor samples with an immune-rich TME harbored HLA down-regulation (Fisher exact one-sided test, p-value = 0.001). The HLA class I and II down-regulated cluster (cluster 4) was composed of 41% (7/17) immune-poor and 59% (10/17) immune-intermediate PCNSL tumors (Figure 2B). Sequencing data were submitted to differential gene expression analysis between clusters 4 (HLA class I and II below median, n = 8) and 1 (HLA class I and II above median, n = 6) and revealed pathways related to immune activation (Figure 2C-D, Table S3). Moreover, patients with complete reduction of HLA expression relapsed earlier than patients who maintained HLA expression (Figure 2E), confirming a link between HLA expression maintenance and immune-rich TME in PCNSL. A trend between a complete reduction of HLA gene expression and overall survival was also observed (Figure S3A).

We next investigated the expression of the different classes of HLA separately. HLA class I gene down-regulation was more frequent in the immune-poor than the immune-rich PCNSL tumor samples (Fisher exact two-sided test, p-value = 9×10⁻⁴). Differential gene expression analysis on sequencing data between tumors with a low (cluster 2 and 4, n = 11) versus higher than median (cluster 3, n = 6) HLA class I gene expression revealed pathways exclusively related to immune activation, e.g., costimulatory molecules, cytokines, T cell activation, and cytotoxicity (Figure 2F-G, Table S4). Moreover, survival analysis showed that HLA class I gene expression diminution was associated with earlier relapse in patients with PCNSL (Figure 2H) and

tended to be associated with OS despite limited significance due to the small cohort size (Figure S3B).

<table><caption>**Table 1.** Comparison of clinical data in the three immune groups of PCNSL (RNA-seq data)</caption><thead><tr><th></th><th>Patients</th><th>Immune<br>poor (n=4)</th><th>Intermediate<br>(n=8)</th><th>Immune<br>rich (n=8)</th><th>χ² test<br>(p-value)</th></tr></thead><tbody><tr><td><strong>Patients</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean age, years (range)</td><td>64<br>(43-81)</td><td>74 (62-81)</td><td>63 (43-80)</td><td>60 (43-70)</td><td>0.06#</td></tr><tr><td>Male, n (%)</td><td>10 (50)</td><td>2 (50)</td><td>4 (50)</td><td>4 (50)</td><td rowspan="2">1</td></tr><tr><td>Female, n (%)</td><td>10 (50)</td><td>2 (50)</td><td>4 (50)</td><td>4 (50)</td></tr><tr><td><strong>MSKCC prognostic class</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Class 1, n (%)</td><td>4 (20)</td><td>0</td><td>3 (37)</td><td>1 (13)</td><td rowspan="3">0.08</td></tr><tr><td>Class 2, n (%)</td><td>7 (35)</td><td>2 (50)</td><td>0</td><td>5 (63)</td></tr><tr><td>Class 3, n (%)</td><td>9 (40)</td><td>2 (50)</td><td>5 (63)</td><td>2 (25)</td></tr><tr><td><strong>Outcome</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CR at last follow up, n<br>(%)</td><td>9 (45)</td><td>2 (50)</td><td>1 (12.5)</td><td>6 (75)</td><td>0.02*</td></tr><tr><td>AWD at last follow up,<br>n (%)</td><td>2 (10)</td><td>0 (0)</td><td>1 (13.5)</td><td>1 (12.5)</td><td></td></tr><tr><td>DOD, n (%)</td><td>9 (45)</td><td>2 (50)</td><td>6 (75)</td><td>1 (13.5)</td><td></td></tr><tr><td>Relapse or progression,<br>n (%)</td><td>12 (60)</td><td>2 (50)</td><td>7 (88)</td><td>3 (38)</td><td>0.11</td></tr><tr><td><strong>Tumor cell phenotype</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Non-GC phenotype-<br>positive, n (%)</td><td>15 (75)</td><td>4 (100)</td><td>5 (63)</td><td>6 (75)</td><td>0.56</td></tr><tr><td>GC phenotype-<br>positive, n (%)</td><td>4 (20)</td><td>0 (0)</td><td>2 (25)</td><td>2 (25)</td><td></td></tr><tr><td>unclassified</td><td>1 (5)</td><td>0 (0)</td><td>1 (12)</td><td>0 (0)</td><td></td></tr><tr><td>cMYC-positive, n (%)</td><td>9 (45)</td><td>2 (50)</td><td>3 (38)</td><td>4 (50)</td><td>0.86</td></tr><tr><td>BCL2-positive, n (%)</td><td>15 (75)</td><td>2 (50)</td><td>7 (88)</td><td>6 (75)</td><td>0.38</td></tr><tr><td>cMYC/BCL2-<br>positive, n (%)</td><td>8 (40)</td><td>2 (50)</td><td>3 (38)</td><td>3 (38)</td><td>0.9</td></tr><tr><td>BCL6-positive, n (%)</td><td>17 (85)</td><td>4 (100)</td><td>5 (63)</td><td>8 (100)</td><td>0.07</td></tr><tr><td>CD10-positive, n (%)</td><td>4 (20)</td><td>0 (0)</td><td>2 (25)</td><td>2 (25)</td><td>0.53</td></tr><tr><td>EBER-positive, n (%)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1</td></tr><tr><td>P53-positive, n (%)</td><td>5 (25)</td><td>2 (50)</td><td>0 (0)</td><td>3 (38)</td><td>0.09</td></tr><tr><td>PDL1-positive, n (%)</td><td>2 (10)</td><td>0 (0)</td><td>0 (0)</td><td>2 (25)</td><td>0.28</td></tr><tr><td><strong>Cytogenetic</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BCL2-break positive, n<br>(%)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1</td></tr><tr><td>BCL6-break positive, n<br>(%)</td><td>5 (25)</td><td>0 (0)</td><td>2 (25)</td><td>3 (38)</td><td>0.36</td></tr><tr><td>cMYC-break positive, n<br>(%)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1</td></tr><tr><td>PDL1-break positive, n<br>(%)</td><td>1 (0.05)</td><td>0 (0)</td><td>0 (0)</td><td>1 (13)</td><td>NA</td></tr></tbody></table>

MSKCC: Memorial Sloan-Kettering Cancer Center; CR: complete remission; AWD: alive with disease; DOD: dead of disease; # Comparison of age of patients with immune-poor tumors to age of patients with immune-rich tumors (Wilcoxon two-sided test).

<table><caption>**Table 2.** Comparison of the available clinical data in the three immune groups of PCNSL (Microarray data)</caption><thead><tr><th></th><th>Patients</th><th>Immune<br>poor (n=6)</th><th>Intermediate<br>(n=16)</th><th>Immune<br>rich (n=12)</th><th></th></tr></thead><tbody><tr><td><strong>Patients</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean age, years (range)</td><td>64 (44-76)</td><td>64 (47-74)</td><td>63.4 (50-72)</td><td>65 (44-76)</td><td>0.56#</td></tr><tr><td>Male, n (%)</td><td>19 (56)</td><td>4 (66)</td><td>8 (50)</td><td>7 (58)</td><td>0.76</td></tr><tr><td>Female, n (%)</td><td>15 (44)</td><td>2 (33)</td><td>8 (50)</td><td>5 (52)</td><td></td></tr><tr><td><strong>MSKCC prognostic class</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Class 1, n (%)</td><td>3 (9)</td><td>1 (20)</td><td>1 (6)</td><td>1 (8)</td><td>0.78</td></tr><tr><td>Class 2, n (%)</td><td>13 (38)</td><td>3 (60)</td><td>5 (31)</td><td>5 (41)</td><td></td></tr><tr><td>Class 3, n (%)</td><td>9 (53)</td><td>2 (40)</td><td>10 (62)</td><td>6 (50)</td><td></td></tr></tbody></table>

MSKCC: Memorial Sloan-Kettering Cancer Center; # Comparison of age of patients with immune-poor tumors to age of patients with immune-rich tumors (Wilcoxon two-sided test).
<figure><img src="image_4.png" /><figcaption>**Figure 2.** HLA expression is related to an immune-rich TME in PCNSL. A. HLA class I and II gene expression in PCNSL. High/low HLA gene expression defined the HLA status: HLA normal, HLA I, II, and I & II down (see heatmap legend). Expression was assessed independently in each cohort; cohort 1 is mRNA-sequencing and cohort 2 is microarray data. B. Samples from the distinct immune subtypes were counted in each HLA gene expression cluster of panel A. C. Differentially expressed genes according to distinct HLA status (HLA normal in pink versus HLA I, II, and I & II down in green): 103 significantly deregulated genes were selected that perfectly segregated the samples according to their HLA status (FDR < 0.01, log2-FC > 4 (absolute value), average read counts > 20). D. Main Gene Ontology Biological Process (GOBP) terms found significantly enriched in genes in panel C (hypergeometric test, FDR < 0.05, at least 3 deregulated genes in each GO term). E. HLA-down status (I, II, I & II) associates with earlier relapse in PCNSL (Kaplan-Meier curves, log-rank test, n = 14, normal = above median, down = below median). F. Differentially expressed genes between PCNSL with normal HLA class I gene expression (orange) and HLA I down (dark blue): 62 significantly deregulated genes were selected which perfectly segregated the two groups (FDR < 0.01, log2-FC > 4 (absolute value), average read counts > 20). G. Main GOBP terms found significantly enriched in genes in (f) (hypergeometric test, FDR < 0.05, at least 3 deregulated genes in each GO term). H. HLA I-down status associates with earlier relapse in PCNSL (Kaplan-Meier curves, log-rank test, n = 17, normal = above median, down = below median).</figcaption></figure>

Dependence between HLA class II down-regulation and the immune subtypes was significant (Fisher exact two-sided test, p-value = 0.05). Small homozygous deletions are known to affect the *HLA-DR* gene in PCNSL [31] and *HLA-DRA* expression was reported as a prognostic factor in DLBCL [34]. We hence compared *HLA-DRA* expression between the three immune subtypes and found it more highly

expressed in the immune-rich group than in the immune-poor PCNSL (Fisher exact one-sided test, p-value = 4.4·10⁻⁴) or the immune-intermediate (Fisher exact one-sided test, p-value = 0.01) tumors. In particular, *HLA-DRA* expression was correlated with activated CD8+ T cells (spearman correlation r = 0.71, p-value = 4.10⁻⁴) and Th1 cells (spearman correlation r = 0.53, p-value = 0.0027). Differentially expressed
genes between *HLA-DRA* normal versus loss of gene expression samples identified genes related to the innate immune response, cytokines, antigen presentation, and extracellular matrix (ECM) organization (Figure S4A-B, Table S5). Contrary to patients with DLBCL, patients with PCNSL showed no association between *HLA-DRA* gene expression and outcome (Figure S4C-D).

Immune subtypes of PCNSL and related
signaling pathways

Next, we explored the potential relations
between the PCNSL immune subtypes identified in
our present study and the deregulated signaling
pathways discussed in the literature on PCNSL
[7-9,11-14,35,36]. The NF-κB, Signal transducer and
activator of transcription 3 (STAT3), IFNγ,

Phosphoinositide 3-kinase/Protein kinase B (PI3K/AKT), Kirsten ras sarcoma viral oncogene (KRAS), Vascular endothelial growth factor (VEGF), P53, MAPK, Salvador-Warts-Hippo (HIPPO), Interleukin-10 (IL-10), TNF-α, WNT/β-catenin, Transforming growth factor beta (TFG-β), and NOTCH pathways were investigated. We retrieved their gene signatures from the Molecular Signatures DataBase (MSigDB) [21]. Average z-scores of signature genes were computed as previously described [37]. Our investigation aimed at assessing whether the cell type composition may influence the activation of certain signaling pathways.

We observed distinct activations of the signaling
pathways depending on the immune subtypes of
PCNSL (Figure 3A). The immune-rich tumors
logically featured strong activation of signaling

<figure><img src="image_5.png" /><figcaption>**Figure 3.** The immune landscapes of PCNSL unravel distinct oncogenic signaling. A. Signature scores (average z-scores of signature genes) of different oncogenic signaling pathways according to the immune subtypes defined in Figure 1bc. B. Gene Set Enrichment Analysis performed between the immune-rich and other subtypes (intermediate & poor). We report the normalized enrichment score (NES) and indicate statistical significance (FDR < 0.05) by a pink dot (otherwise blue). C. Quantification of NOTCH, WNT, and HIPPO signature scores in the different immune subtypes (Wilcoxon one-sided tests, p-value < 0.05, n = 4+8+8 = 20, one immune-poor subtype outlier was removed (significant according to Grubbs and Dixon tests). D. Correlation between the main cell-type scores (EC, CAF, Act CD8, Tregs, Act DC, TAM, MDSC) and oncogenic signaling z-scores. Highly significant correlations are indicated only (spearman rank correlation coefficient, r > 0.6, p-value < 0.01). E. Correlation between STAT3 and IFNγ signature scores and CD274 gene expression in the PCNSL cohort (spearman rank correlation coefficient, r > 0.6, p-value < 0.0001). Note that correlation is observed with both mRNA-sequencing (RNA-seq) and microarray (MA) data. F. High IFNγ and STAT3 signaling z-scores were associated with lower relapse-free survival (RFS) in PCNSL (Kaplan-Meier curves, log-rank test, n = 20, high = above median, low = below median).</figcaption></figure>
pathways known to be involved in the immune response, e.g., STAT3, IFNγ, IL-10, TNF-α, and NF-κB, as well as in stromal signaling, i.e., TGF-β. Activation of KRAS and P53 signaling appeared also favored in immune-rich PCNSL samples, although in a less pronounced manner and not exclusively. In contrast, WNT/β-catenin, NOTCH, and HIPPO pathways were more activated in the immune-poor PCNSL samples. GSEA was performed to confirm the differences observed between the immune-rich and the other PCNSL subtypes (Figure 3B). Signature scores for NOTCH, WNT/β-catenin, and HIPPO for each immune subtype are featured in Figure 3C.

We then assessed correlations between the main cell types we used to profile PCNSL TME (Figure 1B-C), i.e., activated CD8+ T cells, activated CD4+ T cells, Tregs, macrophages, DC, EC, and CAF, and the different signaling pathways (Figure 3D). The TGF-β z-score was correlated with CAF abundance as expected, while the four immune cell types (activated CD8+ T cells, Tregs, TAMs, and DCs) were highly correlated with IFNγ, TNF-α, IL-10, and STAT3 z-scores. Angiogenesis, as well as WNT/β-catenin, NOTCH, and HIPPO pathways correlated with the EC signature score.

STAT3 and IFNγ signaling pathways are generally considered coupled to PD-L1 up-regulation in solid tumors and DLBCL [38]. Although a consensual good prognostic value for IFNγ activity has been acknowledged, hyper-activation of STAT3 signaling has been associated with varying survival outcomes [38-40]. Given a link between PD-L1-expressing TAMs and STAT3 signaling has been reported in PCNSL [41], we exploited a PCNSL retrospective cohort of 57 patients [16] that we previously described in the assessment of an association between PD-L1-expressing TAMs and survival. We defined two histological scores for this purpose: TAM density, evaluated by CD68-positive staining, and the percentage of PD-L1-positive TAMs. High PD-L1 protein expression in TAMs was defined as the simultaneous detection of high TAM density (grades 2-3) and high PD-L1 expression in these TAMs (>50% of CD68-positive staining). High PD-L1-expressing TAMs were found to be strongly associated with better survival and a lower risk of relapse in patients with PCNSL (Figure S5A). At the gene expression level, the PD-L1 transcript (CD274) was highly correlated with STAT3 and IFNγ signature scores in RNA-sequencing and microarray data (Figure 3E). Segregating patients according to STAT3 and IFNγ median expressions resulted in an association of high STAT3 and IFNγ expressions with lower relapse risk (Figure 3F). We also observed a trend with OS but not a significant association (Figure S5B-C).

## Mapping PCNSL intercellular interactions

Intercellular interactions, particularly ICs, within the TME are known to contribute to tumor progression and therapy resistance. In this context, specific ligand-receptor (L-R) interactions, e.g., PD-1/PD-L1, have been extensively studied in solid tumors and PCNSL [12, 16, 41, 42]. We recently proposed an algorithm to infer L-R interactions from bulk transcriptomics [24]. Here, this algorithm was applied to mRNA-sequencing and microarray data separately (Figure 4A, Tables S6-7) to identify candidate L-R pairs that arise in the PCNSL TME. This resulted in the identification of 165 confident L-R pairs, from which brain cell-related interactions were discarded as potential background noise, yielding a total of 128 PCNSL-specific confident L-R pairs (Table S8). Since our algorithm associates each L-R pair to signaling pathways in order to test receptor downstream activity, we summarized recurrent pathways in Figure 4B. These pathways largely reflect the oncogenic signaling pathways already investigated above in Figure 3. We then searched for immune infiltration-related L-R pairs by computing a score (the L-R score), indicative of the L-R pair co-expression level in a given sample [24] and correlated that score with immune infiltration (sum of all the immune cell signatures) as reported in Figure 1. We found 46 correlated L-R pairs (Spearman correlation, r>0.5, adjusted p-value < 0.05, Figure 4C). Among these, we denoted several ICs, e.g. CD86/CTLA4, LGALS9/HAVCR2, LILRB2 and its ligands, as well as inflammatory pairs, e.g., B2M/CD247, and CCR5 and its ligands. Other pairs, e.g., ANGPT1/TEK, DLL1/NOTCH1, THBS1/ITGB1, were linked to angiogenesis according to the pathways used by our algorithm. Some L-R pairs, e.g., SELPLG/ITGB2, PDGFB/LRP1, C3/C3AR1, and LGALS9/HAVCR2, were expressed in PCNSL devoid of lymphoid cells (Figure 4C).

Among the 128 PCNSL-specific confident L-R pairs, we found 26 pairs significantly associated with RFS. These included eight pairs significantly associated with both OS and RFS (Figure 4A and Table S8). LGALS9/HAVCR2 L-R scores were higher in the immune-rich PCNSL samples (Figure 4D) and were significantly associated with RFS (Figure 4E). T-cell immunoglobulin mucin receptor 3 (TIM-3, HAVCR2 gene) is an IC receptor that mainly plays a role in T-cell exhaustion. It suppresses T cell responses upon its binding to galectin-9 (LGALS9 gene). However, TIM-3 has demonstrated several behaviors depending on context [43] and hence its role in the brain TME of PCNSL demands
investigation. Protein expression of TIM-3 and its ligand, galectin-9, was quantified by digital imaging in 32 PCNSL patient tumor samples from a retrospective cohort we published recently [16] and in one postmortem normal brain sample (Figure 5A-B). In the normal brain tissue, we observed galectin-9 protein expression in glial cell ramification, ECs, and rare macrophages, while TIM-3 was rarely expressed (only by a very few ECs and macrophages) (Figure 5A). In the tumors, we found up-regulation of both TIM-3 and galectin-9 expression (Figure 5B-C). TIM-3 and galectin-9 protein expression were also correlated in tumor samples (Figure 5D). TIM-3 was mainly expressed by tumor cells, TAMs, and small lymphocytes, whereas galectin-9 was mainly expressed by TAMs, ECs, glial cells, and gemistocytes. Notably, galectin-9 was strongly expressed in the tumoral area and glia, characterized areas of brain inflammation. The characterization of the cell type morphology was assessed by expert pathologists (VLS, VCM). Given the clinical relevance of TIM-3/galectin-9, we looked at its role within the TME of PCNSL. TIM-3 (*HAVCR2*) and galectin-9 (*LGALS9*) gene expression were highly correlated with most IC ligands and/or receptors, e.g., HVEM, HVEML, LAG3, PD-L1, IDO1, and CD86 (Figure 5E). We also found that *LGALS9* and *HAVCR2* gene expressions were highly correlated with gene signatures of cell types that were enriched in the immune-rich subtype of PCNSL, i.e., TAMs, Th1, MDSC, and TFh (Figure 5F). Finally, *LGALS9* and *HAVCR2* expression was highly correlated with STAT3 and IFNγ z-scores, and correlated well with TNF-α, MAPK, and IL-10 signaling. *HAVCR2* gene expression was also correlated with TGF-β and P53 (Figure 5G).

Discussion

Bulk transcriptomic analysis of PCNSL patient samples allowed us to recognize three well-defined subtypes related to the immune infiltrate abundance and composition based on T cell presence. Other cell types [44–46] have recently been shown to play an important role in TME or therapeutic resistance, but their importance for PCNSL has to be further investigated. We defined the T cell-based subtypes as immune-rich, intermediate, and poor subtypes. These subtypes were independent of clinical parameters except for the patient outcome. Moreover, comparing PCNSL and DLBCL, we found immune-poor tumors enriched in PCNSL, whereas DLBCL tumors were relatively more abundant in the immune-rich subtype. This indicates that the particular brain microenvironment could influence the immune response in PCNSL physiopathology.

We found specific associations when correlating the three PCNSL immune subtypes with signaling pathways commonly deregulated in this entity [7-9,11-14,35,36]. Indeed, inflammatory pathways, e.g., IFNγ or NF-κB, as well as anti-inflammatory pathways, e.g., STAT3 or TNFα, were found active in immune-rich PCNSL. In addition, IL-10, TGF-β signaling, and angiogenesis mediated by VEGF are all STAT3 activators [47] and were also found up-regulated in the immune-rich PCNSL. Immune-poor PCNSL harbored activated WNT/β-catenin, HIPPO, and NOTCH signaling. All the deregulated pathways in the immune-poor subtype were correlated with ECs and are known to play a role in angiogenesis [48,49]. No clear association was found with the immune-intermediate subtype, which is not surprising given its heterogeneous immune cell composition.

Pathways related to ICs are of particular interest to tumor biology and the clinics [50]. The PD-1/PD-L1 interaction has been described in PCNSL [16,41] and PD-L1 was found to be expressed by TAMs [16] in this tumor. We show here that high PD-L1-expressing TAMs were associated with better survival and lower relapse. Interestingly, we also showed that two important PD-1/PD-L1-regulated mechanisms (STAT3 and IFNγ signaling) were active in the immune-rich subtype and were associated with improved survival. This highlights the prognostic value of PD-L1-expressing TAMs in PCNSL. This result is in agreement with PD-L1 expression in TAMs and its correlation with STAT3 in large DLBCL cohorts, where it was also associated with improved outcome [38]. This result could help identifying an immunologically distinct patient subgroup amenable to ICI therapy.

Another IC interaction, *LGALS9* (galectin-9)/*HAVCR2* (TIM-3), was identified by a new algorithm that we recently described [24]. High *LGALS9/HAVCR2* co-expression scores (the so-called L-R scores) were found in every immune subtype of PCNSL, although they were higher in the immune-rich tumors. Concomitant up-regulation of TIM-3 and galectin-9 proteins occurred in PCNSL tumors compared to normal brain tissue, an observation already reported for other tumors in the brain [43,51,52]. In contrast to PD-1, we show that TIM-3 was mostly expressed by both tumor cells and microglial/TAMs in PCNSL whereas it is mainly expressed by CD8+ T cells, and microglial cells in brain tumors, such as glioma [43]. This result suggests different intercellular communication and immune surveillance escape mechanisms in PCNSL compared to solid brain tumors. *HAVCR2* gene expression was associated to the presence of other IC molecules,
inflammatory pathways, e.g., STAT3 or IFNγ, and the abundance of immune cells found in the immune-rich subtype of PCNSL. Previous reports have suggested a regulatory role of TIM-3 on the expression of other IC molecules, such as PD-L1 in glioma [43], thus implying that anti-TIM-3 agents could be used to strengthen anti-PD-1 agents.

The size of our cohort did not allow to assess an association with patient survival. However, high TIM-3 expression in either T cells [53] or tumor cells [54] was reported in DLBCL and associated with a poor outcome in DLBCL patients. Thus, it would be of

clinical relevance to validate the predictive value of TIM-3 expression in the different cellular components on survival in a larger PCNSL cohort.

The immune-poor subtype of PCNSL featured low HLA class I and II molecule gene expression, whereas immune-rich PCNSL associated with higher HLA expression. The immune-intermediate subtype predominantly harbored low expression of HLA class I or both class I and II genes. Low HLA class I gene expression was associated with low immune infiltration, an absence of T-cell activation, and earlier relapse, in all corroborating that low HLA class I gene

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Ligand-receptor interactions within the PCNSL microenvironment. A. Ligand-receptor (L-R) pair selection strategy: our algorithm selected 673 and 393 confident L-R pairs from the mRNA-sequencing (n = 20) and microarray (n = 34) datasets respectively. A total of 128 confident and PCNSL-specific L-R pairs were selected from both datasets. Out of these 128 L-R pairs, survival analysis found 26 L-R pairs significantly associated to RFS and 8 to overall survival (OS). B. Functional categories associated to the L-R pairs selected. C. Immune infiltrate-associated L-R pairs in our PCNSL cohort (n = 20). D. LGALS9/HAVCR2 L-R scores are higher in the immune-rich subtype of PCNSL (Wilcoxon one-sided tests, p-value < 0.01, n = 20 = 8+12 and n = 2+22 for mRNA-sequencing and microarrays data respectively). E. LGALS9/HAVCR2 above median L-R scores are associated with better RFS (Kaplan-Meier curve, log-rank test, p-value < 0.05).</figcaption></figure>
expression results in poor prognosis in solid cancers [55,56]. This observation is also in agreement with the bad prognosis linked to low immune infiltration we observed in PCNSL. In our study, *HLA DRA* gene down-regulation was associated with antigen presentation, innate and adaptive immune responses, activated CD8+ T cells and Th1 cells. Although *HLA-DRA* gene expression was recognized as an independent adverse prognosis factor for RFS in DLBCL [57], we did not find any association with patient outcome in PCNSL, a result that may be linked to the poorer overall outcome of patients with PCNSL compared to DLBCL [58].

Altogether, our results highlight the similarities between the immune-rich, intermediate, and poor subtypes of PCNSL and the hot, intermediate, and cold subtypes observed in primary testicular lymphoma (PTL) [53] and in solid tumors [59]. Despite the obvious differences in terms of spatial architecture between PCNSL and solid tumors, analogies could be drawn between PCNSL immune-intermediate tumors displaying reduced HLA gene expression and solid and altered tumors [60,61]. Herein, knowledge on these PCNSL immune subtypes could help in stratifying patients prior to treatment selection. Most of our patients (18/20) were

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Galectin-9/TIM-3 crosstalk is up-regulated and linked to immune activities in PCNSL. A. TIM-3 protein expression is up-regulated in PCNSL compared to normal brain tissue. B. Galectin-9 (Gal-9) protein expression is up-regulated in PCNSL tissue compared to normal brain tissue. C. TIM-3 and galectin-9 expression are up-regulated in PCNSL tumors in comparison to normal brain tissue (student t-test). D. TIM-3 and galectin-9 expression are correlated in our PCNSL cohort (spearman rank correlation, n = 20). E. Correlation between either HAVCR2 (TIM-3) or LGALS9 (galectin-9) gene expression and other immune checkpoints (spearman rank correlation, adjusted p-value < 0.01). F. Correlation between either HAVCR2 (TIM-3) or LGALS9 (galectin-9) expression and the cell-type z-scores (spearman rank correlation, adjusted p-value < 0.05). G. Correlation between either HAVCR2 (TIM-3) or LGALS9 (galectin-9) gene expression and PCNSL related signaling pathway z-scores (spearman rank correlation, adjusted p-value < 0.05).</figcaption></figure>
treated with immunochemotherapy containing rituximab. Despite the modest size of our cohort due to the rarity of PCNSL, we did however show longer RFS for patients harboring an immune-rich tumor and preserved HLA class I and II gene expression. This result emphasizes the informative value of PCNSL immune patterns. One limitation of our study remains the lack of a validation cohort. Nonetheless, a T cell-inflamed signature has been found associated to favorable patient outcome in a PTL cohort treated with a rituximab-containing immunochemotherapy [53].

Several PCNSL clinical trials with ICIs are ongoing. Despite a high response rate in HL, melanoma, and lung cancers [62,63], response to ICIs remains generally heterogeneous [64]. The status of HLA machinery [65], as well as the presence of active signaling pathways (such as STAT3), are considered as immune evasion mechanisms that influence responses to ICIs [66]. Therefore, patients with PCNSL harboring the immune-rich subtype could be further stratified by assessing potential hyper-activation of STAT3 signaling. In such a case, a combined regimen including an ICI and a STAT3 inhibitor could be envisioned [67]. Patients with immune-poor PCNSL are obviously less likely to benefit from ICI. Tumors featuring an intermediate-immune subtype and down-regulation of either HLA class I or II might not respond optimally to ICI monotherapy. Nevertheless, novel immunotherapies that have emerged in lymphoma, such as targeting CAR-T cells, cancer vaccines, or bispecific antibodies [3-5] could help in restoring an immune response in such intermediary immune subtypes.

NF-κB signaling was also found hyper-activated in immune-rich PCNSL. Previous studies have shown that *MYD88* gene mutations are highly prevalent and support lymphoma growth through NF-κB signaling in PCNSL [10-13]. It would hence be reasonable to co-assess the genomic and immune phenotypic statuses of PCNSL to propose combined ibrutinib/ICI therapies to immune rich/*MYD88*-mutated PCNSL patients (ongoing trial NCT03770416).

Finally, reported epigenetic profiles have revealed different methylation profiles among patients with PCNSL [68]. In particular, the *PTPN6* gene promoter region was found highly methylated in 48.5% of PCNSL tumors, leading to STAT3 hyper-activation [69]. Given that epigenetic regulators can cross the blood-brain barrier (BBB), such as DNA methyltransferase (DNMT) inhibitors, one might hypothesize that DNMT inhibitors could restore a normal STAT3 expression in immune-infiltrated PCNSL tissues. For immune-poor PCNSL, another option would be to re-establish the HLA class II

expression in lymphoma cells by the use of histone deacetylase (HDAC) inhibitors [70] that can cross the BBB. These speculations further highlight the profound consequences of PCNSL immune subtype stratification, for instance in the addressing of immune-intermediate tumors.

In conclusion, we characterized the immune landscape of human PCNSL by combining bulk transcriptomic analysis, histopathology, and digital imaging. The immune-rich subtype is associated with HLA expression preservation, activation of specific signaling, e.g., IFNγ, STAT3, or NF-κB, and expression of inhibitory ICs, such as PD-1/PD-L1 and TIM-3/Galectin-9. The immune-poor subtype of PCNSL is characterized by active WNT/β-catenin, NOTCH, and HIPPO signaling, limited presence of active T-cells, and down-regulation of HLA expression. Several immune evasion mechanisms and new potential therapeutic opportunities, including anti-TIM-3, highlight the clinical relevance of PCNSL immune subtype classification.

## Supplementary Material

Supplementary figures and tables. http://www.thno.org/v11p3565s1.pdf

## Acknowledgements

We thank the Centre des Ressources Biologiques of Montpellier University Hospital for access to tissue samples. We also thank Dr Simon Cabello-Aguilar for help preparing this manuscript. MA was supported by a GIRCI SOOM API-K 2016-811-DRC-AC grant, and JC was supported by Fondation ARC PJA 20141201975 and Labex EpiGenMed ANR 10-LABX-0012 grants.

## Competing Interests

The authors have declared that no competing interest exists.

## References

1. Villano JL, Koshy M, Shaikh H, Dolecek TA, McCarthy BJ. Age, gender, and racial differences in incidence and survival in primary CNS lymphoma. Br J Cancer. 2011; 105: 1414-8.

2. Jahnke K, Thiel E, Martus P, et al. Relapse of primary central nervous system lymphoma: Clinical features, outcome and prognostic factors. J Neurooncol. 2006; 80: 159-65.

3. Xu-Monette ZY, Zhou J, Young KH. PD-1 expression and clinical PD-1 blockade in B-cell lymphomas. Blood. 2018; 131: 68-83.

4. June CH, O'Connor RS, Kawalekar OU, Ghassemi S, Milone MC. CAR T cell immunotherapy for human cancer. Science. 2018; 359: 1361-5.

5. Velasquez MP, Bonifant CL, Gottschalk S. Redirecting T cells to hematological malignancies with bispecific antibodies. Blood. 2018; 131: 30-8.

6. Mottok A, Steidl C. Biology of classical Hodgkin lymphoma: implications for prognosis and novel therapies. Blood. 2018; 131: 1654-65.

7. Braggio E, Van Wier S, Ojha J, et al. Genome-wide analysis uncovers novel recurrent alterations in primary central nervous system lymphomas. Clin Cancer Res. 2015; 21: 3986-94.

8. Booman M, Szuhai K, Rosenwald A, et al. Genomic alterations and gene expression in primary diffuse large B-cell lymphomas of immune-privileged sites: The importance of apoptosis and immunomodulatory pathways. J Pathol. 2008; 216: 209-17.
9. Fukumura K, Kawazu M, Kojima S, et al. Genomic characterization of primary central nervous system lymphoma. Acta Neuropathol. 2016; 131: 865-75.

20. Zhou Y, Liu W, Xu Z, et al. Analysis of Genomic Alteration in Primary Central Nervous System Lymphoma and the Expression of Some Related Genes. Neoplasia. 2018; 20: 1059-69.

11. Zorofchian S, El-Achi H, Yan Y, Esquenazi Y, Ballester LY. Characterization of genomic alterations in primary central nervous system lymphomas. J Neurooncol. 2018; 140: 509-17.

12. Ou A, Sumrall A, Phphanich S, et al. Primary CNS lymphoma commonly expresses immune response biomarkers. Neurooncol Adv. 2020; 2: vdaad018.

13. Chapuy B, Roemer MGM, Stewart C, et al. Targetable genetic features of primary testicular and primary central nervous system lymphomas. Blood. 2016; 127: 869-81.

14. Sung CO, Kim SC, Karnan S, et al. Genomic profiling combined with gene expression profiling in primary central nervous system lymphoma. Blood. 2011; 117: 1291-300.

15. Samstein RM, Lee CH, Shoushtari AN, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019; 51:202-6.

16. Alame M, Pirel M, Costes-Martineau V, et al. Characterisation of tumour microenvironment and immune checkpoints in primary central nervous system diffuse large B cell lymphomas. Virchows Archiv. 2020; 476: 891-902.

17. Swerdlow SH, Campo E, Harris NL, Jaffe ES, Pileri SA, Stein H TJ. WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues. revised 4th edition. Swerdlow SH, Campo E, Harris NL, Jaffe ES, Pileri SA, Stein H TJ, Ed. 2017.

18. Kawaguchi A, Iwadate Y, Komohara Y, et al. Gene expression signature-based prognostic risk score in patients with primary central nervous system lymphoma. Clin Cancer Res. 2012; 18: 5672-81.

19. McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012; 40: 4288-97.

20. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016; 32: 2847-9.

21. Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102: 15545-50.

22. Sergushichev AA. An algorithm for fast pranked gene set enrichment analysis using cumulative statistic calculation. bioRxiv. 2016; 060012.

23. Cabello-Aguilar S, Alame M, Kon-Sun-Tack F, Fau C, Lacroix M, Colinge J. SingleCellSignalR: inference of intercellular networks from single-cell transcriptomics. Nucleic Acids Res. 2020.

24. Alame M, Cornillot E, Cacheux V, et al. The molecular landscape and microenvironment of salivary duct carcinoma reveal new therapeutic opportunities. Theranostics. 2020; 10: 4383-94.

25. Bankhead P, Loughrey MB, Fernández JA, et al. QuPath: Open source software for digital pathology image analysis. Sci Rep. 2017; 7.

26. Charoentong P, Finotello F, Angelova M, et al. Pan-cancer Immunogenic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 2017; 18: 248-62.

27. Becht E, Giraldo NA, Lacroix L, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016; 17: 218.

28. Riemersma SA, Oudejans JJ, Vonk MJ, et al. High numbers of tumour-infiltrating activated cytotoxic T lymphocytes, and frequent loss of HLA class I and II expression, are features of aggressive B cell lymphomas of the brain and testis. Journal Pathol. 2005; 206: 328-36.

29. Nijland M, Veenstra RN, Visser L, et al. HLA dependent immune escape mechanisms in B-cell lymphomas: Implications for immune checkpoint inhibitor therapy? Oncoimmunology. 2017; 6: e1295202.

30. Chang C, Lin CH, Cheng AL, Medeiros LJ, Chang KC. Primary central nervous system diffuse large B-cell lymphoma has poorer immune cell infiltration and prognosis than its peripheral counterpart. Histopathology. 2015; 67: 625-35.

31. Booman M, Douwes J, Glas AM, et al. Mechanisms and effects of loss of human leukocyte antigen class II expression in immune-privileged site-associated B-cell lymphoma. Clin Cancer Res. 2006; 12: 2698-705.

32. Muris JJF, Meijer CJLM, Cillessen SAGM, et al. Prognostic significance of activated cytotoxic T-lymphocytes in primary nodal diffuse large B-cell lymphomas. Leukemia. 2004; 18: 589-96.

33. Ennishi D, Takata K, Béguelin W, et al. Molecular and genetic characterization of MHC deficiency identifies ezh2 as therapeutic target for enhancing immune recognition. Cancer Discov. 2019; 9: 546-63.

34. Rimsza LM, Roberts RA, Miller TP, et al. Loss of MHC class II gene and protein expression in diffuse large B-cell lymphoma is related to decreased tumor immunosurveillance and poor patient survival regardless of other prognostic factors: A follow-up study from the Leukemia and Lymphoma Molecular Profiling Project. Blood. 2004; 103: 4251-8.

35. Montesinos-Rongen M, Brunn A, Bentink S, et al. Gene expression profiling suggests primary central nervous system lymphomas to be derived from a late germin center B cell. Leukemia. 2008; 22: 400-5.

36. Rubenstein JL, Fridlyand J, Shen A, et al. Gene expression and angiotropism in primary CNS lymphoma. Blood. 2006; 107: 3716-23.

37. Bolen CR, McCord R, Huet S, et al. Mutation load and an effector T-cell gene signature may distinguish immunologically distinct and clinically relevant lymphoma subsets. Blood Adv. 2017; 1: 1884-90.

38. McCord R, Bolen CR, Koeppen H, et al. PD-L1 and tumor-associated macrophages in *de novo* DLBCL. Blood Adv. 2019; 3: 531-40.

39. MaHa MA, Matta A, Kaur J, et al. Prognostic significance of nuclear pSTAT3 in oral cancer. Head Neck. 2011; 33: 482-9.

40. Kusaba T, Nakayama T, Yamazumi K et al. Activation of STAT3 is a marker of poor prognosis in human colorectal cancer. Oncol Rep. 2006; 15: 1445-51.

41. Miyasato Y, Takashima Y, Takeya H, et al. The expression of PD-1 ligands and IDO1 by macrophage/microglia in primary central nervous system lymphoma. J Clin Exp Hematopathol. 2018; 58: 95-101.

42. Sugita Y, Furuta T, Ohshima K, et al. The perivascular microenvironment in Epstein-Barr virus positive primary central nervous system lymphoma: The role of programmed cell death 1 and programmed cell death ligand 1. Neuropathology. 2018; 38: 125-34.

43. Kim H-S, Chang CY, Yoon HJ, et al. Glial TIM-3 modulates immune responses in the brain tumor microenvironment. Cancer Res. 2020; 80: 1833-45.

44. Ruffell B, Coussens LM. Macrophages and therapeutic resistance in cancer. Cancer Cell. 2015; 27: 462-72.

45. Guillerey C, Huntington ND, Smyth MJ. Targeting natural killer cells in cancer immunotherapy. Nat Immunol. 2016; 17: 1025-36.

46. Petitprez F, de Reyniès A, Keung EZ, et al. B cells are associated with survival and immunotherapy response in sarcoma. Nature. 2020; 577: 556-60.

47. Yu H, Kortylewski M, Pallord D. Crosstalk between cancer and immune cells: role of STAT3 in the tumour microenvironment. Nat Rev Immunol. 2007; 7: 41-51.

48. Saltarella I, Frassanito MA, Lamanuzzi A, et al. Homotypic and Heterotypic Activation of the Notch Pathway in Multiple Myeloma-Enhanced Angiogenesis: A Novel Therapeutic Target? Neoplasia. 2019; 21: 93-105.

49. Wang X, Freire Valls A, Schermann G, et al. YAP/TAZ Orchestrate VEGF Signaling during Developmental Angiogenesis. Dev Cell. 2017; 42: 462-78.

50. Wester HJ, Keller U, Schottelius M, et al. Disclosing the CXCR4 expression in lymphoproliferative diseases by targeted molecular imaging. Theranostics. 2015; 5: 618-30.

51. Li G, Wang Z, Zhang C, et al. Molecular and clinical characterization of TIM-3 in glioma through 1,024 samples. Oncoimmunology. 2017; 6: e1328339.

52. Yuan F, Ming H, Wang Y, et al. Molecular and clinical characterization of Galectin-9 in glioma through 1,027 samples. J Cell Physiol. 2020; 235: 4326-34.

53. Leivonen SK, Pollari M, Brück O, et al. T-cell inflamed tumor microenvironment predicts favorable prognosis in primary testicular lymphoma. Haematologica. 2019; 104: 338-46.

54. Chen BJ, Dashmaamoorthy R, Galera P, et al. The immune checkpoint molecules PD-1, PD-L1, TIM-3 and LAG-3 in diffuse large B-cell lymphoma. Oncotarget. 2019; 10: 2030-40.

55. Mizukami Y, Kono K, Maruyama T, et al. Downregulation of HLA Class I molecules in the tumour is associated with a poor prognosis in patients with oesophageal squamous cell carcinoma. Br J Cancer. 2008; 99: 1462-7.

56. Pedersen MH, Hood BL, Beck HC, Conrads TP, Ditzel HJ, Leth-Larsen R. Downregulation of antigen presentation-associated pathway proteins is linked to poor outcome in triple-negative breast cancer patient tumors. OncolImmunology. 2017; 6: e1305531.

57. Tada K, Maeshima AM, Hiraoka N, et al. Prognostic significance of HLA class I and II expression in patients with diffuse large B cell lymphoma treated with standard chemoimmunotherapy. Cancer Immunol Immunother. 2016; 65: 1213-22.

58. Rubenstein JL, Gupta NK, Mannis GN, LaMarre AK, Treseler P. How I treat CNS lymphomas. Blood. 2013; 122: 2318-30.

59. Cristescu R, Mogg R, Ayers M, et al. Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy. Science. 2018; 362: eaa3593.

60. Galon J, Bruni D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat Rev Drug Discov. 2019; 18: 197-218.

61. Hege PS, Chen DS. Top 10 Challenges in Cancer Immunotherapy. Immunity. 2020; 52: 17-35.

62. Ansell SM, Lesokhin AM, Borrello I, et al. PD-1 blockade with nivolumab in relapsed or refractory Hodgkin's lymphoma. New Engl J Med. 2015; 372: 311-9.

63. Brahmer JR, Tykodi SS, Chow LQM, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. New Engl J Med. 2012; 366: 2455-65.

64. Nishino M, Ramaiya NH, Hatabu H, Hodi FS. Monitoring immune-checkpoint blockade: response evaluation and biomarker development. Nat Rev Clin Oncol. 2017; 14: 655-68.

65. Gettinger S, Choi J, Hastings K, et al. Impaired HLA class I antigen processing and presentation as a mechanism of acquired resistance to immune checkpoint inhibitors in lung cancer. Cancer Discov. 2017; 7: 1420-35.

66. Johnson DB, Estrada M V., Salgado R, et al. Melanoma-specific MHC-II expression represents a tumour-autonomous phenotype and predicts response to anti-PD-1/PD-L1 therapy. Nat Comm. 2016; 7: 10582.

67. Johnson DE, O'Keefe RA, Grandis JR. Targeting the IL-6/JAK/STAT3 signalling axis in cancer. Nature Rev Clin Oncol. 2018; 15: 234-48.

68. Nakamura T, Yamashita S, Fukumura K, et al. Genome-wide DNA methylation profiling identifies primary central nervous system lymphoma as a distinct entity different from systemic diffuse large B-cell lymphoma. Acta Neuropathol. 2017; 133: 321-4.
69. Liu J, Yaming W, Sun X, et al. Promoter methylation attenuates SHP1 expression and function in patients with primary central nervous system lymphoma. Oncol Rep. 2017; 37: 887-94.

70. Younes A, Berdeja JG, Patel MR, et al. Safety, tolerability, and preliminary activity of CUDC-907, a first-in-class, oral, dual inhibitor of HDAC and PI3K, in patients with relapsed or refractory lymphoma or multiple myeloma: an open-label, dose-escalation, phase 1 trial. Lancet Oncol. 2016; 17: 622-31.