Abstract
The development of brain metastases hallmarks disease progression in 20-40% of melanoma patients and is a serious obstacle to therapy. Understanding the processes involved in the development and maintenance of melanoma brain metastases (MBM) is critical for the discovery of novel therapeutic strategies. Here, we generated transcriptome and methylome profiles of MBM showing high or low abundance of infiltrated Iba1high tumor-associated microglia and macrophages (TAMs). Our survey identified potential prognostic markers of favorable disease course and response to immune checkpoint inhibitor (ICi) therapy, among them APBB1IP and the interferon-responsive gene ITGB7. In MBM with high ITGB7/APBB1IP levels, the accumulation of TAMs correlated significantly with the immune score. Signature-based deconvolution of MBM via single sample GSEA revealed enrichment of interferon-response and immune signatures and revealed inflammation, stress and MET receptor signaling. MET receptor phosphorylation/activation maybe elicited by inflammatory processes in brain metastatic melanoma cells via stroma cell-released HGF. We observed phospho-METY1234/1235 in a subset of MBM and observed marked response of brain metastasis-derived cell lines (BMCs) that lacked druggable BRAF mutations or developed resistance to BRAF inhibitors (BRAFi) in vivo to MET inhibitors PHA-665752 and ARQ197 (tivantinib). In summary, the activation of MET receptor in brain colonizing melanoma cells by stromal cell-released HGF may promote tumor cells self-maintenance and expansion might counteract ICi therapy. Therefore, therapeutic targeting of MET possibly serves as promising strategy to control intracranial progressive disease and improve patient survival.
Introduction
The interaction of brain colonizing tumor cells with the tumor microenvironment (TME), mainly comprising innate and adaptive immune cells, microglia, astrocytes, neurons and oligodendrocytes crucially determines the developmental stages of brain metastases (BM). Brain metastases are observed in 20 – 40% of melanoma patients during the course of disease and micrometastases are evident in more than 75% of autopsied brains1. Hence, only a subset of melanoma cells that entered the brain develop symptomatic and detectable BM during the lifetime of melanoma patients. Unlike peripheral metastases, the emergence of BM depends on a plethora of environmental cues such as the spatiotemporal availability of factors that are provided by cells of the TME, supporting or repressing tumor cell growth2. Moreover, single-cell RNA sequencing (scRNAseq) studies have confirmed regional heterogeneity of astrocytes3, oligodendrocytes4 and microglia5,6 in healthy human brains. Particularly, astrocytes and microglia adopt a reactive cell state7,8 that accompanies secretion of pro- and anti-inflammatory factors under pathological conditions9, 10. It is therefore possible that subfractions of astrocytes and microglia communicate and react with tumor cells in different ways. Probably, neuroinflammation precedes colonization of the brain by tumor cells. However, tumor cells invading the brain amplify inflammatory processes mediated by astrocytes and infiltrating tumor-associated microglia and macrophages (TAMs)11,12. Recently, signaling mediated by hepatocyte growth factor (HGF) and the related receptor MET (c-MET, HGFR) was identified as the trigger of reactive microglia13. HGF is thereby secreted by microglia in the context of trauma but also under normal conditions and seems to play a special role in the growth and self-renewal of neural stem cells in the subventricular zone (SVZ) of rat brains14,15. Therefore, metastatic melanoma cells expressing MET might scavenge HGF from the brain for activation of processes downstream of MET mediating survival and proliferation.
Here, we used transcriptome and methylome profiling to unravel the epigenetic and transcriptomic landscapes of MBM that featured infiltration of TAMs with emphasis on the potential role of microglia in the activation of the HGF/MET receptor signaling pathway. The MET receptor inhibitors PHA-665752 and tivantinib (ARQ197) effectively blocked the growth of brain metastases derived cells (BMCs). Hence, targeting MET receptor signaling might serve as a potent therapeutic target for brain metastases lacking druggable BRAFV600 mutations.
Results
A microglia-specific gene cluster discriminates MBM
Microglia are a unique population of antigen-presenting cells in the central nervous system (CNS) that are capable of clearing the brain of microbes, dead cells and protein aggregates16. Besides, microglia play a crucial role during injury repair and display an exceptional role in immune surveillance and tumor clearance17,18. Although the role of tumor-associated microglia and macrophages (TAMs) in primary brain tumors such as glioblastoma19–22 has been intensively studied, their role in the progression of brain metastases remains poorly understood.
We performed immunohistochemistry (IHC) of our MBM cohort (Supplementary table 1 and23) to determine the levels of activated, Iba1high TAMs. Although Iba1/AIF1 serves as a well-established marker, reactive microglia cannot be distinguished from brain infiltrated macrophages24. Initial studies of MBM revealed that Iba1/AIF1 levels classified tumors into highly and lowly TAM infiltrated (Figure 1a, Supplementary figure 1a). Moreover, we observed overlapping patterns of infiltration of Iba1high TAMs and CD3+ T cells (Figure 1b). As CD3 only provided information about levels of T cell infiltration, we used the ESTIMATE algorithm25 to gain insight into the overall degree of immune cell infiltration of MBM. In line with our previous observation, tumors with intensive TAM and T cell infiltration exhibited a high immune score (Pts 3, 4, 10, 12) whereas MBM with low levels of Iba1high/ CD3+ cell infiltration (Pts 1, 2) or low expression of Iba1/AIF1 (Supplementary figures 1b, c) showed low immune scores (Figure 1c). As expected, brain metastases derived cell lines (BMCs) with absence of immune cells featured lowest scores (Supplementary figures 1b, c). The brain has long been considered a sanctuary where tumor cells can grow undisturbed and protected from attack by immune cells. Therefore, we next investigated expression levels of Iba1/AIF1 in brain (MBM) and extracranial metastases (EM). We observed AIF1 expression in both as well as a high correlation with immune score (Figure 1d), suggesting a relationship between levels of infiltration of TAMs and immune cells not only in the brain. As high levels of immune cell/T cell infiltration are generally associated with good prognosis26, we determined the probability of survival related to Iba1/AIF1 expression of patient’s with (study EGAS00001003672) and without (TCGA-SKCM) MBM. We observed beneficial effects of high Iba1/AIF1 levels in the TCGA cohort (HR=0.46 (0.35 – 0.62), logrank p=1.3e-07) (Figures 1e, f), however, Iba1/AIF1 level had no beneficial effects on the survival of MBM patients. Since no data on TAM-infiltrated MBM are available, we performed comparative methylome and transcriptome profiling of Iba1high (n=5-10) and Iba1low (n=2-6) tumors and identified a set of 417 differentially methylated genomic regions (DMRs) that corresponded to 294 MBM expressed genes (Figure 1g) a core set of markers (n=31) sufficient to split tumors (Figure 1h; Supplementary tables 2, 3). Among them, we identified the integrin family member and gut-homing receptor ITGB7 -which we described in our previous study as a distinguishing mark between BRAF and NRAS mutant MBM23 - and APBB1IP (amyloid b precursor protein-binding family b member 1 interacting protein). Both are associated with better prognosis in patients with colorectal cancer27,28 and clustered with known TAM-associated genes such as P2RY12 and AIF1 (Figure 1h). Remarkably, all clustered tumors were associated with a high immune score. A correlation analysis of clustered genes revealed a high degree of correlation among each other (Figure 1i) and association with hepatocyte growth factor (HGF) that was recently connected with microglia activation13. However, only some of the identified markers within the gene cluster were specifically expressed in microglia but not in brain infiltrating macrophages or other brain cells such as APBB1IP (Figure 1j). The latter gene which has been identified as a conserved microglial gene29 and binding partner of amyloid precursor protein (APP), Tau, 14-3-3g, and glycogen synthase kinase 3 b (GSK3 b) was associated with actin dynamics and retinoic acid signaling30,31. Expression of APBB1IP was significantly (MBM: R=0.86, p<2.2e-16) correlated with immune score (Figure 1k) and survival of melanoma patients (Supplementary figures 1d-e). Moreover, our survey identified a differentially methylated side (Supplementary table 4) within the promoter of PD-L2 (PDCD1LG2) that may predict progression-free survival in melanoma patients receiving anti-PD-1 immunotherapy32. PD-L2 expression was associated with favorable survival (p=0.020) of patients with MBM (Supplementary figure 1f). We found additional genes among our cluster that were expressed in TAMs and significantly associated with immune score (Supplementary figures 1g-n).
Expression of ITGB7 serves as indicator of immune cell infiltration
Recent studies have shown that ITGB7 plays a critical role in the recruitment of T cells to the intestine and that downregulation of ITGB7 is important in protecting intestinal tumors from attack by activated T cells27,33. Hence, we sought to investigate ITGB7 in more detail. Mining of publicly available immune cell data (studies GSE14677134, DICE database35) revealed expression of ITGB7 across different immune cell stages including naïve and memory subsets of T cells, B cells and NK cells (Figure 2a and Supplementary figure 2a). We found that ITGB7 was rather expressed in MBM with infiltration of immune cells and particularly within immune cell dense areas (Supplementary figure 2b). Co-staining revealed accumulation of CD3+ T cells as well as of Iba1high TAMs (Figure 2b). Ranking of MBM regarding levels of ITGB7 expression showed co-occurrence in the expression of CD4, CD274, Sushi Domain Containing 3 (SUSD3) and ITGB7 level (Figure 2c) and validated a possible, previously observed23 correlation of ITGB7 and SUSD3. Moreover ITGB7, SUSD3 and APBB1IP showed expression across different immune cell types except for monocytes and NK cells (Supplementary figures 2c-f). Global (850k) methylome profiling uncovered four epigenetic regulation sites of ITGB7 (Supplementary table 4) with two sites that were associated with expression levels and immune score (Figures 2d, left and center panel, Supplementary figure 3a), located in a proximal enhancer-like region (probe cg26689077) or near by the promotor of ITGB7 (probe cg01033299). The latter site was also identified in the TCGA-SKCM cohort. The sites did not correlate with the BRAF mutation status of MBM (Figures 2d, right panel) in contrast to additional two sides that were found within intergenic regions including an CpG island located between exons 4 and 5 (probes cg11510999 and cg18320160; Supplementary figures 3b-e). Hence, methylome profiling of MBM identified two DMRs within the ITGB7 gene that might serve as indicators of the degree of immune cell infiltration.
A recent study demonstrated that MBM feature a lower T cell content than matched extracranial metastases, however response rates to ICi of both were comparable36. Assuming that ITGB7 expression might be crucial for T cell recruitment, we ascertained ITGB7 levels in MBM (n=79) and EM (n=59; study EGAS00001003672). We observed that ITGB7 was expressed in both metastatic subtypes and was significantly correlated (MBM: R=0.51, p=1.8e-06; EM: R=0.69, p=1.1e-09) with the tumoŕs immune scores (Figure 2e). As we suggest that ITGB7 expression might indicate the degree of immune cell infiltration and possibly serve as indicator of response to ICi, we next performed correlation analysis of ITGB7 and known markers of T cells and B cells. We observed a high concordance with immune cell-related but not tumor cell-related genes (NGFR, MITF, MLANA, SLC45A2) and correlation with expression of PDCD1LG2 (PD-L2) and SUSD3, irrespective of the side of metastasis (Figures 2f, g). In line with previous observations, ITGB7 was expressed in primary and metastatic tumors (TCGA-SKCM) and like SUSD3 was associated with favored survival (Figure 2i). In summary, our survey identified a set of markers that are potentially associated with the level of TAM/immune cell infiltration, particularly ITGB7 might serve as a marker for a favorable course of the disease.
A signature-based deconvolution revealed MET receptor signaling in microglia-enriched MBM
Our previous survey identified a set of markers that potentially characterize a molecular subset of MBM, likely showing a favorable course and response to ICi therapy37,38. To further characterize this set of tumors, we performed single-sample Gene set Enrichment-Analysis (ssGSEA) using defined immune-related and gene signatures specifying signaling processes such as MET receptor or STAT3 signaling among others that are reported to be involved in immune-response mechanisms (Supplementary table 5). We observed that MBM featuring a high immune score showed activation of MET and STAT3 signaling, increased tumor inflammation, stress and senescence (SenMayo39) (Figure 3a). Moreover, deconvolution revealed the presence of reactive microglia, astrocytes and immune cell subsets, among them stem cell-like CD8+ T cells (TCF7)40 in tumors, absent in BMCs. CD8+ (TCF7) T cells are necessary for long-term maintenance of T cell responses and predicted positive clinical outcome41,42. Signatures clearly discriminated MBM and BMCs and reinforced the differences of Iba1high (Pts 3, 4) and Iba1low/neg (Pts 1, 2) tumors. We therefore suggest that the activation of MET- or STAT3-mediated signaling processes or those related to stress/senescence or inflammation strongly depend on the composition of the tumor microenvironment, likely determining the response to therapeutic interventions. Although infiltration of TAMs is not evident in all MBM, microglia infiltration seems to be an early occurring process observed ∼21d after intracranial injection of BMCs into brains of immune compromised Crl:CD1-Foxn1nu mice23 (Figure 3b). Moreover, we observed activation of Stat3 signaling in tumor adjacent cells (Figure 3b), suggesting that brain microenvironmental cells are activated after a short time of tumor-stroma interaction and establish an inflammatory environment. We performed ssGSEA and applied the above mentioned signatures and observed a comparable pattern of enrichment in a more comprehensive and independent set of MBM (study EGAS0000100367243, n=79 MBM) (Supplementary figure 4a).
HGF or scatter factor (SF) is the only identified ligand of MET, plays a pivotal role during neural development, regulating growth and survival of neurons15,44 and likely serves as inducer of reactive microglia by an autocrine loop in response to trauma or neurodegenerative disorders13. Therefore, MET expressing, brain colonizing melanoma cells may benefit and take advantage of the HGF-controlled systems naturally occurring in the brain. We observed HGF expression among tumors of different data sets comprising MBM, EM and primary tumors (studies EGAS00001005976; TCGA-SKCM; EGAS00001003672) with no significant difference of HGF levels in tumor subsets (Figures 3d, e). Investigation of immune cell and brain cell data (DICE database35 and study GSE7372145) revealed high expression of HGF in monocytes and astrocytes (Figures 3f, g), suggesting a potential role of different stroma cell populations for activating HGF/MET signaling in brain-infiltrating tumor cells. Assuming that the degree of microglia infiltration determines signaling processes in MBM cells, we explored expression levels of MET signaling-associated genes in tumors with high and low level of infiltrated TAMs and found levels of HGF, PIK3CG, PTK2B, STAT3 and MAP4K1 significantly correlated with microglia score (Supplementary figure 4c-e) that was defined as average expression level or β-value of microglia markers APBB1IP, SYK, HCK and P2RY12 (Supplementary table 6).
HGF might be released by immune cells as well as homeostatic and reactive microglia or astrocytes. We surveyed the Seattle Alzheimeŕs Disease Brain Atlas which is implemented in the Allen brain atlas database (https://portal.brain-map.org/) and observed that dementia fostered expansion of microglia with increased expression of HGF (Figure 3h, center and right panels). Reactive microglia and immune cell released HGF might hence be responsible for activation of growth factor/survival signaling in adjacent tumor cells. In line with previous studies, we observed a significant correlation of HGF expression with immune score in brain (BM, R=0.49, p=5.3e-06) and extracranial metastases (EM, R=0.41, p=1.5e-03), (Figure 3i).
Expression and activation of MET receptor classifies a molecular subset of MBM
Understanding the molecular mechanisms that establish cellular dependencies and thus control the development and maintenance of brain metastases is critical for their therapeutic manipulation. Recently, we identified that the expression of Ecad and NGFR sufficiently discriminated molecular subsets of MBM23. These subsets likely distinctly interact with microenvironmental cells and respond to therapeutics (Figure 4a). To identify potential druggable targets, we surveyed the pan-MBM, NGFR and Ecad-specific gene sets for cell surface receptors that may serve as crucial key factors that control tumor cell maintenance and expansion and identified 24 receptors that distinguished Ecad+ and NGFR+ tumors (Figure 4b). Particularly ADIPOR1 (adiponectin receptor 1, p=1.9e-02), SIRPA (signal regulatory protein alpha, p=1.1e-05) and PLXNC1 (plexin C1) showed significantly increased expression in Ecad+ MBM and EM but comparable levels among MBM and EM (Figure 4c). In addition, Ecad+ MBM featured increased levels of MET receptor in (p=1.4e-04). MET was significantly (p=2.7e-05) higher expressed in MBM than EM (Figure 4d, left and center panels). The MET tyrosine kinase receptor pathway serves as a potent survival and maintenance factor for MBM and might be a promising therapeutic target46. MET expression was associated with increased cell cycle progression and proliferation (Figure 4d, right panel) and defined yet another subset of MBM (Figure 4e). Next we assessed whether expressed MET indeed participated in active signaling processes. Phosphorylation of MET at tyrosine residues 1234/1235 (pMETY1234/1235) is critical for kinase activation and initiation of downstream processes and was evident in nearly all METhigh MBM investigated, independent of the BRAF mutation status (Figures 4f, g). MET receptor alterations are evident in 9 % of all SKCM melanoma cases, including amplification as observed in 1.13 – 17.2% or 11% of melanoma (TCGA-SKCM, study by Ramani et al.47). However, targeted DNA sequencing (TargetSeq) and fluorescence in-situ hybridization (FISH) revealed absence of MET activating mutations and a tendential MET amplification in only one case (Pat 5, Supplementary figure 5a, b). However, all but one tumor (Pat 14) showed high polysomy. As we assume that environmental cells foster activation of MET receptor signaling in a subset of tumor cells, we performed co-IHC for pMETY1234/1235 and Iba1. We observed pMETY1234/1235 positive tumor cells in close proximity to Iba1high TAMs (Figure 4h), though MET receptor was not activated in Iba1high microglia that resided in adjacent normal tissue (Supplementary figure 5c, upper panel). However, MET receptor activation was also evident in scattered tumor cells in the absence of adjacent Iba1high TAMs (Supplementary figure 5c, lower panel) suggesting paracrine mechanisms or additional sources of HGF such as immune cells or astrocytes. Considering that HGF levels, like those of other growth factors provided by stromal cells, might depend on spatial factors, we examined the Allan Brain Atlas database and found that HGF is comparably expressed in different brain sections (frontal lobe (FL), parietal lobe (PL), temporal lobe (TL), occipital lobe (OL)) but is lowly abundant in the brainstem (pons) (supplemental Figure 5d). The spatially dependent expression of growth factors in the brain may therefore determine the dependencies of the tumor cells.
Interferon signaling determines response of MBM to immune checkpoint inhibitor therapy
Interferon-gamma signaling has been identified as an important mechanism for upregulation of PD-L1 on melanoma cells and escape from immune recognition. On the other hand, recent studies uncovered that high interferon-gamma-related gene expression signature scores (IFN-γ score) were associated with low risk of melanoma relapse from neoadjuvant ipilimumab plus nivolumab therapy48,49.
In our recent study, we observed significant enrichment of interferon and inflammatory response (“Hallmark”, MsigDB50) signatures in MBM with high level of tumor infiltrating lymphocytes (TILhigh)23 that have been attributed with favored survival in a pre-clinical melanoma model49. We found overlapping expression of ITGB7, SUSD3 and HGF and Hallmark interferon-response genes, separating MBM of our cohort and MBM of study EGAS00001003672 (Figure 4i, Supplementary figure 6a). Expression of ITGB7 significantly correlated with levels of interferon regulatory factor 1 (IRF1) and IRF8 in MBM (BM) and extracranial metastases (EM) of study EGAS00001003672 (Supplementary figure 6b-d). Moreover, we observed high correlation of levels of HGF, IRF1 and IRF8 in MBM (Supplementary figure 6e). As microglia serve as a source of soluble receptor ligands such as Hgf, we next surveyed data of interferon-gamma treated (1 U/mL IFNγ, 24h) murine microglia cells (BV2, GSE132739). Indeed, we found significant upregulation of Itgb7 (p = 2.9e-03) and Hgf (p=4.4e-02) but downregulation of Susd3 (p=4.0e-02) in BV2 cells (Supplementary figure 6f). For control, we investigated levels of known interferon-responsive genes that were significantly increased 24h after interferon treatment, Mx1 (p=1.2e-02), PD-L1/Cd274 (p=3.6e-02), Irf1 (p=3.1e-02) Cxcl9 (p=4.0e-03) and Aif1 (p=3.9e-04) (Supplementary figure 6g). In order to classify MBM of our study into anti-PD-L1 responsible and non-responsible and for linking ITGB7, SUSD3 and HGF with therapy response, we performed ssGSEA and applied interferon responsive and additional immune response gene signatures (of study GSE18634451). Our survey validated that ITGB7, SUSD3 and HGF were highly expressed in MBM that featured enrichment of interferon responsive genes/signatures (Pts. 3-6, 12; Supplementary figure 6h). Hence, we suggest that ITGB7, SUSD3 and HGF like PD-L1 are among the interferon-regulated genes triggered by immune cell-released interferon-gamma and may be involved in immune response mechanisms of MBM.
The targeting of MET receptor serves as a promising strategy to control MBM growth
Although a subset of MBM exhibit immune cell subset enrichment and interferon response signatures and respond to ICi therapy, MET-expressing brain metastatic melanoma cells may benefit from HGF released by stromal cells to drive progression. Hence, activation of MET signaling may depend on the degree of tumor-stroma interaction, possibly counteracting the beneficial impact of immune checkpoint inhibition (ICi). Resistance-mediating processes include the phosphorylation of ribosomal protein S6 (pS6), which is downstream of MET and mTOR signaling52 and was observed in progressive BRAFi-resistant melanomas53,54.
We assessed pS6 phosphorylation of serine residues 235/236 and found co-occurrence of activated MET receptor and of pS6235/236 in MITF positive tumors (Figure 5a and Supplementary figures 7a, b). Moreover, pS6 phosphorylation was evident in a BRAFwt (T2002) and mutated (V600E, BMC53) cell lines probably suggesting a general activation of pS6 signaling irrespective of the presence of mutated BRAF (Figure 5b). As MET signaling might serve as mediator of a resistance-mediating program, we assessed the efficacies of the ATP-competitive inhibitor PHA-665752 and the non-ATP-competitive, clinical phase I/II MET receptor inhibitor (METi) tivantinib (ARQ197) in BMCs that showed variable levels of MET expression (Figure 5c and Supplementary figure 7c). ARQ197 failed to improve the outcome and overall survival of patients with hepatocellular carcinoma55 but may potentially be effective in melanoma patients. The initial testing revealed a general response of BMCs (BMC1-M1, BMC53), T2002 cells and conventional cell lines (A375, A2058, MeWo) to both inhibitors irrespective of the BRAF mutation status (Figures 5d, e). As we observed a mutually exclusive rather than co-expression of MET receptor and NGFR (nerve growth factor receptor), we tested whether the manipulation of NGFR levels might affect the response to PHA-665752 (PHA). We observed that overexpression of NGFR in A375 cells (A375NGFR) sensitized to METi compared to RFP expressing control (A375RFP) or MeWo cells (Figure 5f). Next, we asked whether METi targeting may serve as alternative therapeutic strategy for BRAFi resistant (BMC4) or cells with non-BRAFV600 mutations (BMC2) showing only moderate or no response to dabrafenib (Figure 5g) as indicated by IC50 values (BMC4, IC50 = 226.4 nM and BMC2, IC50 = 3029.5 nM). The broad range (1nM – 10µM) testing of tivantinib in BMCs, T2002 and conventional melanoma cells and PHA in BMCs revealed that all cell lines responded to both METi, irrespective of the mutation status. However, we observed that the non-(brain) metastatic cell lines A375, A2058, and T2002 were more sensitive to treatment with tivantinib than BMCs (Figures h-j, Supplementary figure 7d). The median IC50 value of BMCs was ∼600 nM (range: 406.5 – 800.1). Cell lines lacking BRAF and NRAS mutations (MeWo, T2002) showed highest responses to tivantinib (Figure 5k).
In summary, brain metastatic as conventional melanoma cell lines responded to METi, suggesting that targeting of MET signaling might be a promising tool for the treatment of non-BRAFV600 and BRAFV600 mutated MBM that acquired resistance to BRAFi or for combinatorial of METi and ICi in NRAS mutated tumors.
Discussion
The spatiotemporal development of primary and secondary brain tumors is strongly determined by the crosstalk of tumor and brain micronenvironmental cells, particularly macrophages, astrocytes and microglia56 and the consequential activation of inflammatory processes57. Although the neuro-inflammatory processes that are activated alongside development and progression of primary brain tumors such as glioblastoma have been intensively studied, the mechanisms that accompany emergence of brain metastases on the other side are not well investigated.
Here, we used combined transcriptome and methylome profiling to unravel the molecular features of MBM of different progression stages showing high and low level of tumor-associated macrophages/microglia (TAMs) infiltration, irrespective of the phenotype (Ecad, NGFR). Generally, TAMs foster development and progression of primary brain tumors58,59, however their functional role in MBM may be different. We observed that MBM containing a high proportion of TAMs were associated with a high immune score and infiltration of CD3+ T cells. The profiling of Iba1/AIF1high tumors revealed a cluster of genes, among them ITGB7, APBB1IP as SUSD3 and PD-L2, that were widely expressed among immune cell subtypes and previously associated with increased immune T cell infiltration23,27,28,33 and favored outcome. Previous mouse studies demonstrated a pivotal role of ITGB7 for intestinal T cell recruitment and correlated low levels of Itgb7 with colorectal cancer progression and maintenance of intestinal stem cells via Ecad-mediated interaction27,33. However, we observed strong protein expression of ITGB7 in immune cells adjacent to tumor cells of Ecad+ and NGFR+ tumors, suggesting a broader function of ITGB7 in different subtypes of metastases and cancers. Enhanced expression of ITGB7 might be a prerequisite for immune cell invasion. Therefore, epigenetic marks that correlate with the expression of ITGB7 and other genes mentioned above may be of prognostic importance, and the expression of these markers could determine the pathways of intracranial progression.
As previously described for the Ecad+ and NGFR+ subtypes of MBM, whether tumors are enriched or depleted in TAMs and immune cell subsets is critical and may determine response to therapeutic intervention. The subsequent ssGSEA-based characterization of MBM of studies performed by us and others revealed molecular programs fostering or accompanying the TAM+/TIL+ tumor subtype. TAM+/TIL+ tumors featured activation of MET and STAT3 signaling, increased stress response, tumor inflammation, senescence and activation of microglia and astrocytes but also activated interferon signaling. STAT3 activation in tumor-adjacent astrocytes in response to brain damage or tumor cells is well-investigated process60,61 and was rapidly induced in response to brain infiltrating BMCs. Hence, enrichment of STAT3 signature genes was likely attributed to tumor-adjacent astrocytes and infiltrated immune cells62.
The HGF/MET receptor signaling plays a pivotal role during brain development and neuro-regeneration, homeostasis of microglia and neurons14,44 but is also involved in microglia activation in response to trauma13,44,63. Brain infiltrating melanoma cells hence may engage the HGF/MET signaling of brain cells and utilize for regulation of survival and proliferation. We observed expression of MET receptor in the subset of E-cadherin (Ecad) expressing tumors23, suggesting that Ecad+ but not NGFR+ cells may depend on HGF/MET signaling. As we observed that HGF is expressed by immune cell subsets and homeostatic or reactive astrocytes and microglia, we investigated the level of activated/phosphorylated MET receptor in TAM-adjacent tumor cells. We found that tumor cells but not Iba1high TAMs that resided in tumor cell-free adjacent stroma showed activation of MET, however MET was also activated in the absence of adjacent microglia in some tumor cells, suggesting a paracrine effect of HGF. In line with our previous study23, we observed enrichment of interferon-response signatures in the subset of TILhigh/immune score (IS)high tumors. We observed enrichment of interferon-response genes in MBM with high levels of ITGB7 expression and observed significant response of Itgb7 and Hgf among known interferon-inducible genes such as Cd27464 and Mx165 in interferon-gamma treated BV2 murine microglia cells (unpublished study GSE132739). Hence, T cell-provided interferon-gamma might not only induce expression of Cd274/PD-L1 but may also activate expression of HGF and ITGB7. Therefore, autocrine MET receptor signaling might be triggered in response to immune cell released interferon-gamma and/or paracrine activation of MET signaling may occur via (INFG-activated) reactive glia-released HGF.
Our study bridges the gap between the immune cell phenotype of MBM and the activation of potentially therapeutic counteracting signaling pathways. The infiltration of TAMs and immune cells thus represents a double-sided sword and on the one hand is associated with an effective response to immune checkpoint inhibitors, but on the other hand can support the growth of MET expressing tumor cells via secreted factors such as HGF.
Therefore, we finally assessed the potential role of small molecule inhibitors of MET receptor (METi) for targeting of MBM that lack druggable BRAFV600 mutations or developed refractory disease. To this end, we took advantage of well-characterized BMCs serving as in vitro model systems. We observed that the ATP-competitive inhibitor PHA-665752 and the non-ATP-competitive, clinical phase II inhibitor ARQ197 (tivantinib) elicited response in BMCs and conventional melanoma cell lines irrespective of the BRAF/NRAS mutation status. However, although being effective at doses of 100 – 200 nM in MeWo and A375 cells, ARQ197 showed a median IC50 value of ∼1 µM in BMCs, suggesting a general difference among brain metastatic and long-term maintained conventional cell lines established from either non metastatic (A375) or locally metastatic (MeWo) cells.
In summary, we have shown that MET receptor signaling is active in a subset of MBM, conferring a survival/growth benefit independent of BRAF/NRAS mutation status. MET activation may occur in response to HGF released by TAM/immune cells and could counteract therapeutic interventions. Furthermore, we suggest interferon-induced expression of HGF in tumor cells triggered by interferon-gamma provided by stromal cells mediates autocrine activation of MET-signaling tumor cells (Figure 6). In addition, we demonstrated that methylome profiling of MBM has high potential to identify gene regulatory sites that may predict favorable progression of intracranial disease. In the present study, we identified epigenetic regulatory sites in a group of genes comprising ITGB7, APBB1IP, SUSD3 and PD-L2 (PDCD1LG2).
Limitations of the study
The present study is not without limitations. Although we suggested that HGF/MET signaling is activated in tumor cells in close proximity to infiltrated microglia, we have not provided evidence that growth of established brain tumors and phosphorylation of MET decreases in response to METi. Moreover, whether METi are capable of passing the blood-brain barrier and not affect normal homeostatic processes e.g. those crucial for neuron survival needs to be investigated.
Methods
Patient cohorts
All procedures performed in this study were in accordance with the ethical standards of the respective institutional research committees and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All patients gave written informed consent for the collection and scientific use of tumor material which was collected at the Biobank of the Charité – Comprehensive Cancer Center (CCCC). The study was approved by the Ethics Committee of the Charité (EA1/152/10; EA1/107/17; EA4/028/18 and EA1/107/17 and EA1/075/19) and Universitätsmedizin Greifswald (BB 001/23).
Cultivation of MBM-derived and conventional melanoma cell lines
All cell lines were cultured as previously reported23. Briefly, all brain metastases-derived cell lines (BMCs) and conventional melanoma cell lines were kept at 37°C/ 5% CO2 and 95% humidity in cell culture medium (DMEM, 4.5 g/L glucose, stabilized glutamine/GlutaMax, pyruvate, Gibco/ThermoFisher) supplemented with 10% fetal bovine (FBS, Gibco) serum and 1% penicillin/streptomycin (P/S) (Gibco/ThermoFisher) and routinely passaged. BMCs were established from intraoperative tumors as previously reported23.
Live cell imaging-based drug sensitivity assays
Drug treatments were performed 24 h after seeding of 2,500-5,000 cells/96-well in 100 µl medium. The response of BMCs and conventional melanoma cell lines to dabrafenib, PHA-665752 or ARQ197 (all purchased from Selleckchem) in a range of 1nM-10µM of eight technical replicates was determined by live cell imaging. Images were taken every three hours using a 10x objective and the general label-free mode, two pictures of eight technical replicates per condition were taken. Drug response was assessed by changes in the cellular density over time. The cell density was determined by a confluence mask tool as part of the IncucyteS3 software. IC50 values were calculated by curve-fitting (https://search.r-project.org/CRAN/refmans/REAT/html/curvefit.html) based on confluence measurements at day 3.
In vivo experiments
All animal experiments were performed in accordance with the German Animal Protection Law under the permission number G0130/20 obtained via the Berlin Ministry of Health and Social Affairs (LaGeSo). ARRIVE 2.0 Guidelines were strictly followed and performed as previously reported23. Briefly, 2.5×104 BMC1-M4 and BMC2 cells were stereotactically inoculated into brains of female Crl:CD1-Foxn1nu nude mice (8-9 weeks of age, 24-26g, Charles River Laboratories) were with using a 1µl Hamilton syringe and a stereotactic frame as described previously66. Tumor growth was tracked by MRI and animals were sacrificed by perfusion with 4% PFA in deep anesthesia after tumors reached a volume of 20 mm³. Following, whole brains were removed, dehydrated, paraffin embedded and sections of 2 µm were used for downstream analyses.
RNA isolation and sequencing
Isolation of total RNA from snap frozen tumors and RNA sequencing was performed as previously reported23. Briefly, 100 ng of total RNA was used for library preparation with TruSeq Stranded total RNA Sample Preparation-Kit and Ribo-Zero Gold Kit (Illumina). Paired-end (2×100 bp) sequencing of RNA libraries with integrity numbers (RIN) ≥7 was performed on NovaSeq6000 platform at Cegat GmbH, Tuebingen (Germany). Following demultiplexing of sequenced reads and adapter trimming67, FASTQ files were obtained. Raw counts of protein-coding genes were normalized using the DESeq2 (https://bioconductor.org/packages/release/bioc/html/DESeq2.html) package68. Differential expression of genes between groups was determined after fitting models of negative binomial distributions to the raw counts. Raw p-values were FDR (false discovery rate)-adjusted for multiple testing and a value below 0.05 for the adjusted p-values were used to determine significant differentially expressed genes.
Gene-set enrichment GSEA/Single-sample GSEA/Scores
GSEA was performed using the most current BROAD javaGSEA standalone version (http://www.broadinstitute.org/gsea/downloads.jsp) and gene signatures of the molecular signature database MsigDB50,69, v7.4. In addition, we performed GSVA/ssGSEA using R packages GSVA70, GSRI, GSVAdata and org.Hs.eg.db and a customized collection of gene signatures including the signatures provided by Biermann et al.71 and own signatures as defined by selected Ecad23, NGFR23, microglia or TME core genes (this study). All gene signatures are shown in Supplementary table 5. Microglia scores were defined as the mean β-value of probes cg24400465 (APBB1IP), cg05128364 (SYK), cg21704050 (P2RY12) and cg03498995 (HCK) or expression levels (log2 FPKM) of these markers. The proliferation index in Figure 4d was defined as mean expression level of the cell cycle regulators PCNA, MKI67, CCNB1 and CCNB2.
Fluorescence in situ hybridization (FISH)
FISH analysis was performed on 4 µm sections of FFPE blocks. Slides were deparaffinized, dehydrated and incubated in pre-treatment solution (Dako, Denmark) for 10 min at 95–99°C. Samples were treated with pepsin solution for 6 min at 37°C. For hybridization, a Vysis MET SpectrumRed/ Vysis CEP 7 (D7Z1) SpectrumGreen Probe (Abbott, Chicago, USA) was used. Incubation took place overnight at 37°C, followed by counterstaining with 4,6-diamidino-2-phenylindole (DAPI). For each case, signals were counted in 50 non-overlapping tumor cells using a fluorescence microscope (BX63 Automated Fluorescence Microscope, Olympus Corporation, Tokyo, Japan). Computer-based documentation and image analysis was performed with the SoloWeb imaging system (BioView Ltd, Israel). MET high-level amplification (MET FISH+) was defined as (a) MET/CEN7 ratio ≥2.0, (b) average MET copy number/cell ≥ 6 or (c) ≥10% of tumor cells with ≥15 MET copies/cell as described in Schildhaus et al72.
Quantitative real-time RT-PCR
RNA isolation from frozen cell pellets was performed with the RNeasy Mini Kit (Qiagen, Germany) and, following the manufacturers protocol as previously reported23.qRT-PCR was carried out on a Step one plus PCR cycler (Applied Biosystems, Germany) for 30–40 cycles. Primers were designed for 55–60°C annealing temperatures. Relative expression levels were calculated with the ΔΔCT method73, normalized to β-actin. Primer sequences are shown in Supplementary table 7.
Immunohistochemistry (IHC)/Immunofluorescence (IF)
Automated immunohistochemical staining was performed on formalin-fixed, paraffin-embedded (FFPE) tissue sections using the BenchMark Ultra (Ventana) autostainer. The following primary antibodies were used: CD3 (anti-CD3ε, Agilent, catalog number: #A045201-2, rabbit, dilution: 1:100), pMET (phospho-MET, Tyr1234/1235, Cell signaling, catalog number: #3077, rabbit, dilution: 1:100), pS6 (phospho-S6 ribosomal protein Ser235/236, Cell signaling, catalog number: #2211, rabbit, dilution: 1:100), IBA1 (IBA1/AIF-1, ionized calcium-binding adaptor molecule 1, Cell signaling, catalog number: #17198, rabbit, dilution: 1:100), ITGB7 (Integrin beta 7, Thermo Fisher, catalog number: #BS-1051R, rabbit, dilution: 1:100) and pSTAT3 (phospho-STAT3, Tyr705, Cell signaling, catalog number: #9145, rabbit, dilution: 1:100) and MITF (clones C5L+LD5, Zytomed, catalog number: Z2161MP, mouse, dilution: 1:100). Primary antibodies were applied and developed using the iVIEW DAB Detection Kit (Ventana Medical Systems) or the ultraView Universal Alkaline Phosphatase Red Detection Kit (Ventana Medical Systems). All slides were counterstained with hematoxylin for 8 minutes. IF of mouse brain sections was performed with IBA1 (IBA1/AIF-1, ionized calcium-binding adaptor molecule 1, Cell signaling, catalog number: #17198, rabbit, dilution: 1:100), KBA.62, NovusBiologicals, catalog number: NBP2-45285, mAb mouse, 1:100; GFAP-AlexaFluor594, BioLegend, catalog number: 644708, mAb mouse.
Data availability
Whole transcriptome and methylome data were deposited in the European Genome-Phenome Archive (EGA), under accession numbers EGAS00001005975, EGAS00001005976 (https://ega-archive.org/studies/). The data are available under controlled access. Supplementary tables of our recent study have been deposited at Zenodo (https://doi.org/10.5281/zenodo.10006881). Supplementary tables of our previous study containing a full list of patient’s characteristics (Supplementary table 1) have been deposited at Zenodo (https://zenodo.org/record/7013097 and https://doi.org/10.5281/zenodo.7249214).
Author contributions
TR performed data analyses, prepared figures and wrote the manuscript; ES collected tumors, established BMCs and performed experiments; KP performed histological analyses; MEW, SN and HWSS performed craniotomy and provided MBM; AV performed experiments/ drug response assays; HR provided expertise on microglia markers; AL performed MET-FISH analyses; SKR provided funding; KJ performed bioinformatics data analysis and expertise on data analysis; JR wrote the manuscript, provided resources and funding. All authors have read and agreed to the published version of the manuscript.
Competing interests
The authors declare no conflicting interests.
Funding
JR is an alumnus of the BIH-Charité Clinical Scientist Program funded by the Charité – Universitätsmedizin Berlin and the Berlin Institute of Health. We thank the German Cancer Consortium (DKTK), Partner site Berlin for technical support. AV received funding by the ÖAW (DOC Fellowship: DOC/26523).
Supplementary figure legends
Supplementary figure 1: Iba1/AIF1 expression separates MBM a.) Immunofluorescence (IF) for Iba1 (red) and GFAP (labeling of reactive and normal astrocytes) of a NRASQ61R mutated MBM (Pat 15) that progressed on treatment with immune checkpoint inhibitors (ICi; ipilimumab) showed strong infiltration of tumor-associated microglia/macrophages (TAMs) and the presence of aggregates of microglia and reactive astrocytes. MBM (without adjacent stromal cells) of Pat 3 showed high infiltration of TAMs. DAPI served as nuclear dye. b.) Box plots representing the levels (FPKM, log2) of Iba1/AIF1 in MBM and BMC of study EGAS00001005976 and MBM and EM of study EGAS00001003672. c.) Per sample representation of expression levels (FPKM, log2) of Iba1/AIF1 in MBM, brain metastases-derived cell lines (BMCs) and brain controls (BC) of studies mentioned in (b). d-f.) Survival analyses of MBM patients (n=67) of study EGAS00001003672 or of TCGA-SKCM study (n=459), featuring high or low level of APBB1IP (d, e) or PD-L2 (PDCD1LG2) expression. Analysis revealed a significant (logrank p=5.5e-07) favorable disease course (HR=0.50, Cox-regression analysis) of APBB1IPhigh melanoma (e) and favorable outcome associated with high levels of PD-L2 expression in MBM (f). Survival of MBM patients was not significantly affected by APBB1IP levels (d). g.-m.) Correlation and cell type-specificity of microglia markers SYK, HCK, P2RY12 and AIF1. n.) Dot plot shows significant (R=0.85, p=1.3e-04) correlation of ITGB7 expression of MBM (n=16, study EGAS00001005976) and immune score. Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (b).
Supplementary figure 2: The indicators of favorable disease course ITGB7, SUSD3 and APBB1IP are broadly expressed among immune cell subsets. a.) Box plots showing expression levels of ITGB7 among T cell subsets (CD4, T helper cells), CD8 (cytotoxic T cells), NK cells (natural killer cells), cDC (conventional dendritic cells), pDC (plasmacytoid dendritic cells) lymphoid follicle-residing B cells (FollicularB cell subsets), innate lymphoid cells, type 3 (ILC3), mast cells, macrophages (Macro), monocyte subsets (Mono), myofibrils (non-immune related cells), plasma cells (PlasmaB) and tumor-associated macrophages (TAMs) as provided by study GSE146771. b.) IHC of MBM of indicated patients for ITGB7. ITGB7 expression is evident in lymphocyte-enriched areas. c.-d.) Boxplot showing expression levels of SUSD3 in immune cell subsets of (a), and shows levels obtained from DICE (Database of Immune Cell Expression, Expression quantitative trait loci (eQTLs) and Epigenomics), suggesting broad but variable expression among immune cell subsets. e.) Comparative illustration of levels of ITGB7 and SUSD3 as requested from DICE. f.) Box plot showing expression levels of APBB1IP among immune cell subsets of the aforementioned GEO study, suggesting a broad expression among immune cell types. Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (a, c, d, f).
Supplementary figure 3: Methylome profiling uncovered epigenetic regulatory sites in the ITGB7 gene. a.) Schematic representation of the ITGB7 gene, showing exons, intronic regions and sites of epigenetic marks as depicted by indicated probes. ITGB7 expression levels are associated with methylation at sites covered by probes cg26689077 and cg01033299 located within a proximal enhancer-like signature and intronic region close to 5’-UTR. Additional regions as covered by probes cg18320160 and cg11510999 are associated with the BRAF mutation status of tumors. b.-c.) Dot plots showing no significant correlation of methylation status (indicated by β-values) at sites covered by probes cg18320160 and cg11510999 and immune score. d.-e.) Box pots indicating a significant association of β-values, determined by aforementioned probes (cg11510999, p = 3e-06; cg18320160, p= 0.003) and BRAF status (BRAFV600 vs. wt/NRASQ61) of MBM. Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (d, e).
Supplementary figure 4: Single sample GSEA revealed classification of immune molecular subtypes of MBM. a.) Single-sample GSEA (ssGSEA)-based deconvolution of MBM (n=79) of study EGAS00001003672 using customized gene signatures indicating “Signaling” processes, cellular subsets and stages of microglia and astrocyte and immune cell subsets. ssGSEA demonstrated distinct separation of MBM with high, median or low immunescore regarding expression levels of signature genes, BMCs served as controls. ssGSEA uncovered differentially activated pathways and processes such as MET and STAT3 and interferon signaling, senescence (SenMayo), stress response and tumor inflammation in tumors enriched for reactive microglia and astrocytes and innate and acquired immune cells subsets. b.) Boxplot showing expression levels of HGF in indicated immune cell subsets obtained from DICE. c.-e.) Expression of MET receptor pathway genes in MBM of study EGAS00001003672 showing high or low enrichment of microglia, as determined by levels of microglia-specific genes (SYK, HCK, AIF1/Iba1= microglia score suggests a significant correlation of microglia infiltration and activation of MET receptor signaling. HGF, hepatocyte growth factor (p=1.8e-09); PIK3CG, Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Gamma (p<2.2e-16); PTK2B, Protein Tyrosine Kinase 2 Beta (p<2.2e-16); STAT3, Signal Transducer And Activator Of Transcription 3 (p=5.3e-11); MAP4K1; Mitogen-Activated Protein Kinase Kinase Kinase Kinase 1 (p=3.4e-12). Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (b-e).
Supplementary figure 5: MET-FISH analysis revealed absence of MET receptor amplifications in MBM. a.) Hematoxylin and eosin (H&E) staining shows tumor cell histology (upper row). Fluorescence in-situ hybridization with a MET-specific probe (red) revealed no specific amplifications of the MET gene in MBM (n=7) as compared with centromere control (green) and irrespective of the BRAF/NRAS mutation status. DAPI served as nuclear dye, bars indicate 50 µm. b.) Quantitative representation of FISH analysis, indicating the number of MET copies per nucleus and ratio of MET and CEP7 (Centromer 7). c.) IHC for Iba1 (red) and pMETY1234/1235 (brown) revealed absence of activated MET in Iba1high microglia residing in adjacent tissue (upper panel) and MET activation in tumor cells without neighboring TAMs. d.) HGF expression in brain cells residing within the different lobes (FL, frontal; PL, parietal; TL, temporal; OL, occipital lobe) and pons as retrieved from the Allan Brain Atlas (https://portal.brain-map.org/). Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (d).
Supplementary figure 6: Expression of interferon-related genes is enriched in MBM of ITGB7high/IScorehigh phenotype. a.) Heat map indicating expression levels and subset-association of interferon-inducible genes, mediators of interferon signaling and relevant immune cell-expressed markers such as CD3E, CD4, CD8A in MBM (n=79) of study EGAS00001003672. Molecular subsets, category of genes and strength of expression are color coded. b.-d.) Dot plots indicating the significant correlation of ITGB7, IRF1 (Interferon Regulatory Factor 1), IRF8 and IFNG in MBM and EM of the aforementioned study. e.) Dot plots indicating the significant correlation of expression of HGF, IRF1 and IRF8 in MBM and EM of the aforementioned study. f.-g.) Investigation of expression data of murine BV2 microglia cells (study GSE132739) following interferon (1 U/mL IFNγ, 24h) or control treatment revealed interferon-responsible genes. Interferon treatment significantly increased levels of Itgb7 (p=2.9e-03), Hgf (p=4.4e-02), Mx1 (p=1.2e-02), Cd274 (p=3.6e-02), Irf1 (p=3.1e-02), Cxcl9 (p=4.0e-03) and Aif1 (p=4.0e-03). However, Susd3 was significantly downregulated upon interferon treatment (p=4.0e-02). h.) ssGSEA analysis of MBM with defined signatures showing enrichment of interferon-related signaling among other indicated processes. Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (f, g).
Supplementary figure 7: The mTOR/pS6 signaling is activated in MBM. a.) IHC of MBM (Pat 14) for levels of activated/phosphorylated MET receptor (pMETY1234/1235) and mTOR/pS6 (pS6S235/236) signaling revealed co-occurrence of both. b.) Co-occurrence of pS6S235/236 and MITF. Bars indicate 50 µm. c.) Confocal microscopy imaging of BMC1-M1 and BMC53 cells for levels of NGFR and MET showing a mutually exclusive expression pattern or low level of MET in NGFR+ cells. d.) Live cell imaging tracked dose-response of BMC1-M1 cells to increasing doses of ARQ197.
Acknowledgments
We gratefully thank Cathrin Müller for excellent technical assistance.