Next Article in Journal
Neuroinflammation in Neurodegenerative Disorders: Current Knowledge and Therapeutic Implications
Next Article in Special Issue
Car T Cells in Solid Tumors: Overcoming Obstacles
Previous Article in Journal
Polydimethylsiloxane Organic–Inorganic Composite Drug Reservoir with Gliclazide
Previous Article in Special Issue
Targeting the Tumor Microenvironment in Breast Cancer: Prognostic and Predictive Significance and Therapeutic Opportunities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Transcriptomic Expression of MSR1 in Solid Tumors Is Associated with Responses to Anti-PD1 and Anti-CTLA4 Therapies

by
Adrián Sanvicente
1,2,
Cristina Díaz-Tejeiro
1,
Cristina Nieto-Jiménez
1,
Lucia Paniagua-Herranz
1,
Igor López Cade
1,3,
Győrffy Balázs
4,5,
Víctor Moreno
6,
Pedro Pérez-Segura
7,
Emiliano Calvo
6,8 and
Alberto Ocaña
1,7,9,*
1
Experimental Therapeutics in Cancer Unit, Medical Oncology Department, Hospital Clínico San Carlos (HCSC), Instituto de Investigación Sanitaria (IdISSC), 28040 Madrid, Spain
2
Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
3
Molecular Oncology Laboratory, Instituto de Investigación Sanitaria San Carlos (IdISSC), 28040 Madrid, Spain
4
Department of Bioinformatics, Semmelweis University, Tűzoltó u. 7-9, H-1094 Budapest, Hungary
5
Research Centre for Natural Sciences, Hungarian Research Network, Magyar Tudosok Korutja 2, H-1117 Budapest, Hungary
6
START Madrid-Fundación Jiménez Díaz (FJD) Early Phase Program, Fundación Jiménez Díaz Hospital, 28040 Madrid, Spain
7
Medical Oncology Department, Hospital Clínico San Carlos (HCSC), Instituto de Investigación Sanitaria San Carlos (IdISSC), 28040 Madrid, Spain
8
START Madrid-HM Centro Integral Oncológico Clara Campal (CIOCC), Early Phase Program, HM Sanchinarro University Hospital, 28050 Madrid, Spain
9
Centro de Investigación Biomédica en Red en Oncología (CIBERONC), 28029 Madrid, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(7), 3987; https://doi.org/10.3390/ijms25073987
Submission received: 5 March 2024 / Revised: 28 March 2024 / Accepted: 29 March 2024 / Published: 3 April 2024

Abstract

:
Immuno-oncology has gained momentum with the approval of antibodies with clinical activities in different indications. Unfortunately, for anti-PD (L)1 agents in monotherapy, only half of the treated population achieves a clinical response. For other agents, such as anti-CTLA4 antibodies, no biomarkers exist, and tolerability can limit administration. In this study, using publicly available genomic datasets, we evaluated the expression of the macrophage scavenger receptor-A (SR-A) (MSR1) and its association with a response to check-point inhibitors (CPI). MSR1 was associated with the presence of macrophages, dendritic cells (DCs) and neutrophils in most of the studied indications. The presence of MSR1 was associated with macrophages with a pro-tumoral phenotype and correlated with TIM3 expression. MSR1 predicted favorable overall survival in patients treated with anti-PD1 (HR: 0.56, FDR: 1%, p = 2.6 × 10−5), anti PD-L1 (HR: 0.66, FDR: 20%, p = 0.00098) and anti-CTLA4 (HR: 0.37, FDR: 1%, p = 4.8 × 10−5). When specifically studying skin cutaneous melanoma (SKCM), we observed similar effects for anti-PD1 (HR: 0.65, FDR: 50%, p = 0.0072) and anti-CTLA4 (HR: 0.35, FDR: 1%, p = 4.1 × 10−5). In a different dataset of SKCM patients, the expression of MSR1 predicted a clinical response to anti-CTLA4 (AUC: 0.61, p = 2.9 × 10−2). Here, we describe the expression of MSR1 in some solid tumors and its association with innate cells and M2 phenotype macrophages. Of note, the presence of MSR1 predicted a response to CPI and, particularly, anti-CTLA4 therapies in different cohorts of patients. Future studies should prospectively explore the association of MSR1 expression and the response to anti-CTLA4 strategies in solid tumors.

1. Introduction

Immuno-oncology has experienced an increase in relevance with the approval of antibodies that have demonstrated clinical activity in different indications [1]. Among these compounds, we can include agents that target programmed death-ligand 1 (PD (L)1) or cytotoxic T-lymphocyte-associated protein 4 (CTLA4) molecules [2]. Very recently, an anti-lymphocyte-activation gene 3 (LAG3) antibody has also been approved for the treatment of patients with metastatic melanoma in the first-line setting [3]. These immunomodulatory proteins induce a suppressor action on effector T cells, therefore repressing a proper immune response [4]. In this context, patients who have greater responses to check-point inhibitors (CPI), mainly anti PD (L)1 agents, are those that have high expression levels of the PD1 ligand (PD L1) within their tumors [1]. Unfortunately, even using this selection criteria, a high number of patients do not respond to these therapies or they experience early progressions within the first months of treatment [5,6]. It is considered that even in a best-case scenario, the responses to CPIs alone will not benefit more than half of the treated population [5,6].
In parallel with this, no biomarker currently exists to predict the efficacy of anti-CTLA4 therapies. Therefore, the identification of molecular markers that can predict responses to these therapies is mandatory to better select patients that would benefit from these agents while avoiding undesirable side effects. Of note, it is very well known that this family of agents has poorer tolerability than anti-PD (L)1 compounds.
Some tumors display an immunosuppressive environment where the action of immune effector cells is repressed through the presence of different CPIs [2]. In addition to the negative effect of tumoral cells on the activation of effector cells, some immune cells such as myeloid-derived cells or Tregs also play a role promoting an immunosuppressive environment [7,8]. In line with this, some cells that belong to innate immunity have demonstrated the potential to play a negative role, as confirmed when agents that act on some of their suppressor receptors augment the efficacy of anti-PD (L) 1 agents [9]. An additional example is the presence of M2 or pro-tumoral macrophages that prevent an effector immune response and have no antigen-presenting capacity [10,11,12].
MSR1, also known as macrophage scavenger receptor-A (SR-A) or cluster of differentiation 204 (CD204), is mainly associated with the uptake and degradation of acetylated low-density lipoprotein (acetyl-LDL) but not non-modified low-density lipoprotein (LDL) [13,14,15]. This produces an increase in intracellular cholesterol deposition, leading to effects similar to those observed in familial hypercholesterolemia [16,17]. Beyond this role, MSR1 has been implicated in different human pathologies and very recently in cancer progression [18,19]. The expression of MSR1 has been linked to significantly poor prognoses and the increased severity of multiple forms of cancer [19,20].
In this article, we evaluated the expression of MSR1 in cancer, observing an association with the presence of macrophages, DCs and neutrophils, in addition to a clear prediction of the response to anti-PD (L) 1 and CTLA4 inhibitors.

2. Results

2.1. MSR1 Is Highly Expressed in Solid Tumors

We first mapped the expression of MSR1 on all available tumor types using publicly available datasets, as described in the material and methods section and in Supplemental Table S1. MSR1 was highly expressed, with a statistically significant difference in some tumors, including glioblastoma (GBM), kidney renal cell carcinoma (KIRC), breast cancer (BRCA), pancreatic adenocarcinoma (PAAD), ovarian carcinoma (OV), stomach adenocarcinoma (STAD), skin cutaneous melanoma (SKCM) and esophagus carcinoma (ESCA), among others (Figure 1A,B). On the other hand, we observed a higher expression in normal tissue compared with transformed tissue in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) (Figure 1A). The tumors with the highest fold changes compared with normal tissues included GBM (≥100), followed by low-grade glioma (LGG) and PAAD, which had ≥10-fold changes. ESCA, STAD, SKCM and OV displayed ≥5-fold changes, and KIRCC, diffuse large B cell lymphoma (DLBCL), thymoma (THYM), head and neck squamous cell carcinoma (HNSCC) and BRCA had ≥2-fold changes (Figure 1B). We did not observe relevant differences between the breast cancer subtypes or ovarian cancer subgroups (Supplementary Figure S1). For further analysis, we set exclusion criteria based on a statistically significant difference between the tumor and normal tissue and a fold change expression equal to or greater than two. With these criteria, we selected the following tumor types for further study: GBM, PAAD, ESCA, STAD, SKCM, OV, KIRC, THYM, HNSC and BRCA. We excluded DLBC as we intended to focus only on solid tumors.

2.2. MSR1 Expression Is Correlated with The Presence of Innate Immune Cells

In a second step, we aimed to explore the expression of MSR1 in relation to the presence of several immunologic cells. We evaluated the following ten selected tumor types: GBM, PAAD, ESCA, STAD, SKCM, OV, KIRC, THYM, HNSC and BRCA, following the criteria described before. As can be seen in Figure 2, a clear and very strong association was observed between the expression levels of MSR1 and the presence of macrophages in almost every tumor type studied, with SKCM being the one with a weaker but still positive correlation (GBM: Rho = 0.645, p = 1.82× 1017; PAAD: Rho = 0.640, p = 4.28 × 10−21; ESCA: Rho = 0.617, p = 3.14 × 10−20; STAD: Rho = 0.675, p = 1.13 × 10−51; SKCM: Rho = 0.289, p = 2.94 × 10−10; OV: Rho = 0.81, p = 3.30 × 10−59; KIRC: Rho = 0.536, p = 1.30 × 10−35; THYM: Rho = 0.493, p = 2.14 × 10−8; HNSC: Rho = 0.513, p = 2.40 × 10−34; and BRCA: Rho = 0.579, p = 3.68 × 10−90). A similar association was identified for the DCs in all tumor types, except for THYM (GBM: Rho = 0.519, p = 7.89 × 10−11; PAAD: Rho = 0.791, p = 6.86 × 10−38; ESCA: Rho = 0.578, p = 1.93 × 10−17; STAD: Rho = 0.708, p = 5.71 × 10−59; SKCM: Rho = 0.369, p = 3.50 × 10−16; OV: Rho = 0.502, p = 2.47 × 10−17; KIRC: Rho = 0.610, p = 2.44 × 10−48; THYM: Rho = −0.320, p = 4.81 × 10−4; HNSC: Rho = 0.360, p = 1.53 × 10−16; and BRCA: Rho = 0.282, p = 1.21 × 10−19) and in the neutrophils (GBM: Rho = 0.609, p = 2.80 × 10−15; PAAD: Rho = 0.694, p = 6.39 × 10−26; ESCA: Rho = 0.455, p = 1.39 × 10−10; STAD: Rho = 0.477, p = 6.30 × 10−23; SKCM: Rho = 0.646, p = 3.08 × 10−55; OV: Rho = 0.557, p = 1.17 × 10−21; KIRC: Rho = 0.738, p = 2.22 × 10−80; THYM: Rho = −0.081, p = 3.90 × 10−1; HNSC: Rho = 0.581, p = 9.22 × 10−46; and BRCA: Rho = 0.525, p = 1.75 × 10−71) (Figure 2A,B).
A strong association with Tregs was only observed for PAAD (Rho = 0.613), and mild associations were observed for STAD (Rho = 0.443), HNSC (Rho = 0.527) and BRCA (Rho = 0.425). There were also associations with CD8 T cells in some tumor types, such as PAAD (Rho = 0.604) and STAD (Rho = 0.669) (Figure 2A). Weaker correlations were identified for the other immune cells evaluated, including B cells and T CD4.
Of note in all these settings, low purity scores were detected, and this was what demonstrated that MSR1 was mainly present in immune cells and not in tumoral cells (Supplementary Figure S2).

2.3. Correlation of MSR1 Expression Level with Macrophage Subtypes and Immune-Suppressive Molecules

As MSR1 has been described as a classical receptor expressed in macrophages, we explored whether its presence was particularly associated with a specific macrophage subtype. As can be seen in Figure 3A, a clear correlation was observed with the expression of MSR1 and the presence of M2 macrophages in the following particular tumor types [10,11,12]: GBM (Rho = 0.427), PAAD (Rho = 0.801), ESCA (Rho = 0.759), STAD (Rho = 0.827), SKCM (Rho = 0.684), OV (Rho = 0.770), KIRC (Rho = 0.746), THYM (Rho = 0.451), HNSC (Rho = 0.781) and BRCA (Rho = 0.717). Tumor-associated macrophages are typically classified as M2s, and they usually do not present antigens to T cells promoting the presence of immune-suppressive molecules. With the intent of characterizing the immune microenvironment associated with the presence of MSR1, we explored the association of MSR1 with different molecules related to immune activation or suppression. As can be seen in Figure 3B, a strong association was observed for the presence of TIM3 with GBM (Rho = 0.727), PAAD (Rho = 0.906), ESCA (Rho = 0.824), STAD (Rho = 0.879), SKCM (Rho = 0.813), OV (Rho = 0.893), THYM (Rho = 0.643), HNSC (Rho = 0.832) and BRCA (Rho = 0.857). Of note, there were also positive correlations for 4-1BB with GBM (Rho = 0.424), PAAD (Rho = 0.643), KIRC (Rho = 0.529), THYM (Rho = 0.492) and BRCA (Rho= 0.512), and for PDL1, there were positive correlations with PAAD (Rho = 0.607), STAD (Rho = 0.484), SKCM (Rho = 0.544), OV (Rho = 0.443), KIRC (Rho = 0.406) and BRCA (Rho = 0.489). Other immune activators, expressed in B cells and DC such as OX40 and CD40, did not correlate with the presence of MSR1 (Figure 3B).
Finally, in the available tumor types including PAAD, SKCM and BRCA, we evaluated the co-expression of MSR1 and co-stimulatory molecules for the same cell types in macrophages and DCs using single cell analysis (Figure 3C). Of note, no co-expression was identified in those cells, demonstrating that TIM3, 4-1BB and PDL1 were not expressed in the same cell lines. These data globally describe the presence of MSR1 in pro-tumoral or M2 macrophages and its association with co-inhibitory molecules such as TIM3, 4-1BB and PDL1.

2.4. MSR1 Predicts Favorable Outcomes in Patients Treated with Check-Point Inhibitors

Next, we assessed the role of MSR1 in the response to CPIs. To do so, we used publicly available data from patients treated with anti-PD(L) 1 and anti-CTLA4 antibodies. As can be seen in Figure 4A, the patients treated with CPIs, including anti-PD(L)1 or CTLA4, had favorable survival rates, as follows: anti-PD1 (HR: 0.56, FDR: 1%, p = 2.6 × 10−5), anti-PD-L1 (HR: 0.66, FDR:20%, p = 0.00098) and anti-CTLA4 (HR: 0.37, FDR:1%, p = 4.8 × 10−5). Of note, 55% of the samples included in this analysis belonged to the same tumor types evaluated earlier. We then aimed to confirm the data in a specific tumor type, where a higher number of patients were available. As can be seen in Figure 4B, the expression of MSR1 was associated with favorable outcomes in the SKCM patients treated with CPIs, as follows: anti-PD1 (HR: 0.65, FDR: 50%, p = 0.0072) and anti-CTLA4 (HR: 0.35, FDR: 1%, p = 4.1 × 10−5). Finally, we compared MSR1 expression with PD-L1 expression as a potential biomarker, using the same population. As shown in Figure 4C, we could observe that PD-L1 was associated with the response to CPIs with a similar accuracy as MSR1 was for anti-PD (L)1 treatment, and much poorer response was associated with anti-CTLA4, as follows: anti-PD1 (HR: 0.54, FDR: 1%, p = 6.3 × 10−5), anti-PD-L1 (HR: 0.54, FDR: 1%, p = 2.8 × 10−7) and anti-CTLA4 (HR: 0.51, FDR: 50%, p = 0.0078).

2.5. MSR1 Expression Predicts Clinical Response in a Different Dataset

To confirm the previous findings, we used a different dataset to evaluate whether patients treated with CPIs experienced a better clinical response based on their expression levels of MSR1. As can be seen in Figure 5A, high expression levels of MSR1 predicted the response to CPI: for the response to any CPI we found: AUC: 0.552, p = 2.5 × 10−3; for the response to anti-PD1, we found AUC: 0.581, p = 1.3 × 10−3; and for anti-CTLA4, we found AUC: 0.612, p = 2 × 10−2. In a specific cohort of SKCM patients, we observed results that trended in the same direction, as follows: for any CPI therapy, AUC: 0.56, p = 2.6 × 10−2; for anti-PD1, AUC: 0.556, p = 4.9 × 10−2; and for anti-CTLA4, AUC: 0.61, p = 2.9 × 10−2 (Figure 5B).

3. Discussion

In the current study, we mapped the presence of MSR1 in solid tumors. The tumor types with the highest MSR1 expression levels included GBM, KIRC, BRCA, PAAD, OV, STAD, SKCM and ESCA, among others. The tumors with the highest fold-change differences compared with normal tissue included GBM, LGG and PAAD. Of note, the expression level of MSR1 in non-transformed lung tissue was higher than that in the tumors.
When evaluating cells related to the expression of MSR1, we observed that a negative correlation with tumor purity was observed, suggesting that the protein was mainly present in the tumor’s microenvironment. Indeed, positive associations were identified with macrophages, DCs and neutrophils. No associations with adaptive cells, including T or B cells, were observed. These data indicated that MSR1 mRNA was mainly present in these cells in the tumors’ microenvironments, except for lung cancer, where it existed at high levels in immune cells within the normal tissue. These findings are in line with the literature regarding several normal and physiological conditions beyond cancer, and they reinforce the concept and role of antigen-presenting cells (APCs) in non-transformed lung [21].
We next observed an association between MSR1 expression and M2 or pro-tumoral macrophages [22,23]. This was in line with previous publications that reported the presence of MSR1 in tumor-associated macrophages (TAM) and correlated that presence with detrimental prognoses for several tumor types [22,23,24,25,26,27,28,29]. Other studies have described their expression with high-grade metastases and tumor aggressiveness [19].
Positive correlations with TIM3 and PD-L1 were identified in the tumor types studied, signifying that MSR1 is abundant in an immune-suppressive environment. To our knowledge, only one study on glioma has suggested the co-expression of MSR1 and TIM3 [28]. A single-cell analysis showed that TIM3 was not present in the same cell populations with higher levels of MSR1 transcripts, with mainly macrophages or dendritic cells. Studies describing an association between PDL1 expression and MSR1 have been published, confirming our findings [30]. An interesting observation was the presence of 4-1BB, which could suggest the existence of already activated and exhausted T cells.
Our study discovered the association of MSR1 expression with response to CPIs, particularly anti-CTLA4 antibodies. Indeed, the presence of MSR1 predicted outcome with more accuracy than PD-L1, particularly for anti-CTLA4 therapies, opening the possibility for this biomarker to be specifically evaluated in patients treated with anti-CTLA4 therapies.
Finally, we used melanoma as a model tumor type to confirm the capacity to predict the response to CPI and, particularly, anti-CTLA4 antibodies. In this context, given the fact that there are no available biomarkers for predicting the response to anti-CTLA4 therapies, the identification of markers that could help to select patients in this clinical setting is a priority, particularly given the toxicity observed with anti-CTLA4 agents.
We acknowledge that our study has limitations. This was an analysis performed using published available datasets in which bulk and single-cell information was extracted. In this context, the confirmation of these data in human samples using immunohistochemistry techniques would undoubtedly help to confirm these results. In addition, we recognize that for the outcome analysis, we could not evaluate some of the tumor types that we used for the correlational study. However, the analysis performed included sixty percent of the populations with those tumor types, and we also used melanoma as a tumor type to confirm the data identified in the whole population.
In conclusion, we described as a potential novel biomarker associated with the response to CPIs, particularly anti-CTLA4, that is mainly expressed in antigen-presenting cells (M2 macrophages) and that is present in an immunosuppressive microenvironment. Further studies should be performed to confirm MSR1’s prediction capacity in patients.

4. Materials and Methods

4.1. Data Collection and Processing

Processed TCGA (The Cancer Genome Atlas; https://www.cancer.gov/ccg/research/genome-sequencing/tcga; last accessed 23 January 2023) PanCancer data were downloaded using different web tools. This dataset contains whole exome sequencing and RNA-Seq information from patients’ tumors and their matched normal tissues. The publicly available web tools that were used to evaluate the expression of the selected genes and the mutational information included Gepia2 [31] (http://gepia2.cancer-pku.cn/#index; last accessed: 23 January 2023), Gent2 [32] (http://gent2.appex.kr/gent2/; last accessed 23 January 2023) and CBioportal [33] (https://cbioportal.org, last accessed 24 January 2023). Tumor types where a fold change greater than 2 between the expression of MSR1 in the tumor vs. normal tissues were selected.

4.2. Immune Cell Infiltration and Gene Expression Correlation

The Tumor Immune Estimation Resource (TIMER) platform [23] was used to investigate the association between gene expression and tumor purity and the association between the presence of tumor immune infiltrates, including CD4+ T cells, CD8+ T cells, Tregs, macrophages, dendritic cells, neutrophils and B cells. TIMER applies a deconvolution method (previously described in [34,35] (http:/timer.cistrome.org/, accessed 25 January 2023)) to infer the abundance of tumor-infiltrating immune cells from gene expression profiles. It contains 10,897 samples from diverse cancer types from the TCGA.
Single-cell data were obtained from the MNP-VERSE Seurat package (downloaded from https://gustaveroussy.github.io/FG-Lab/ accessed 17 February 2023) [36]. The data were preprocessed by filtering and normalizing using the Seurat R package. Cell clustering, dimensionality reduction and visualization were performed using R software (version 4.1.0).

4.3. Gene Correlations

The Spearman correlation coefficients between the expressions of every pair of genes were used for the correlation analysis between them. Data from TCGA [37] were included in the analysis.

4.4. Outcome and Prognosis Analysis

For the evaluation of the activities of the CPIs, we used datasets that included patients treated with these agents. The datasets were identified in the gene expression omnibus (GEO) using the keywords “gene expression”, “PD1”, “CTLA4” and “immunotherapy”, as well as the names of the available immunotherapy agents, as described elsewhere by members of this research group [38]. In this cohort, we evaluated the association of MSR1 expression with overall survival (OS). Patients were separated into two cohorts according to the best cut-off values and analyzed based on the administered therapy, and the anti-PD1 treatment population included 797 patients while the anti-CTLA4 cohort included 131 patients. The Kaplan–Meier (KM) plots were presented with hazard ratios (HRs), 95% confidence intervals (CIs) and log-rank p-values (p). Genes with an HR of <1 and a p-value of <0.05 were considered predictors of favorable outcomes, while genes with an HR of >1 and a p-value of <0.05 were considered predictors of detrimental outcomes. This was performed with the KM Plotter Online tool [39,40] (https:/kmplot.com/analysis/, accessed 25 January 2023). As described in the KM Plotter to avoid missing correlations due to the use of a specific cut-off, all available cut-off values between the lower and upper quartiles of expression were used for the selected genes (as binary variables), and false discovery rate (FDR) using the Benjamini–Hochberg method was computed to correct for the multiple hypothesis testing. The cut-off value with the highest significance (the lowest FDR value) was determined. The ROC plotter online tool [41] was used to correlate gene expression and response with immunotherapy (anti-PD1 or anti-CTLA4) in a cohort of different solid tumors that included metastatic and primary tumors. The area under the curve (AUC) was computed to evaluate the clinical utility (the sensitivity and specificity) of MSR1 as a biomarker. An AUC of 0.5 corresponded to no classification power at all, while an AUC value of 1 denoted a perfect biomarker.

4.5. Datasets Used

Information describing all datasets used in this analysis is provided in Supplementary Table S1.

4.6. Graphical Design

The histograms, bar charts and heatmaps were plotted using GraphPad Prism software (Prism 9.5.1—7 February 2023) (GraphPad Software, San Diego, CA, USA).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25073987/s1.

Author Contributions

Main contribution to the conception of the work, A.S.; design of the work, A.O.; analysis and interpretation of the data, A.S., C.D.-T., C.N.-J., L.P.-H. and G.B.; substantial revision of the manuscript, I.L.C., V.M., P.P.-S., E.C. and A.O. All authors have read and agreed to the published version of the manuscript.

Funding

The results were supported by the Instituto de Salud Carlos III (PI19/00808), ACEPAIN, Diputación de Albacete, CIBERONC and the CRIS Cancer Foundation (to A.O.). The work carried out in the authors’ laboratories received support from the European Community through the regional development funding program (FEDER).

Institutional Review Board Statement

Approved by research board but no ethical approval needed as not applicable” for studies not involving humans or animals.

Informed Consent Statement

Patient consent was waived due to the nature of the original source of the data. This study do not obtain patients data directly but from publicly available datasets previously developed.

Data Availability Statement

All data generated or analyzed during this study are included in this published article [and its Supplementary Information Files].

Acknowledgments

The author would like to thank his wife Raquel Paéz and his children Mateo, Pablo, Rodrigo and Inés Ocana Páez.

Conflicts of Interest

V.M. declares consulting fees from Roche, Bayer, BMS, Janssen, Syneos, Affimed and Astra Zeneca. The principal investigator received institutional funding from AbbVie, AceaBio, Adaptimmune, ADC Therapeutics, Aduro, Agenus, Amcure, Amgen, Astellas, AstraZeneca, Bayer, Beigene BioInvent International AB, BMS, Boehringer, Boston, Celgene, Daichii Sankyo, DEBIOPHARM, Eisai, e-Terapeutics, Exelisis, Forma Therapeutics, Genmab, GSK, Harpoon, Hutchison, Immutep, Incyte, Inovio, Iovance, Janssen, Kyowa Kirin, Lilly, Loxo, MedSir, Menarini, Merck, Merus, Millennium, MSD, Nanobiotix, Nektar, Novartis, Odonate Therapeutics, Pfizer, PharmaMar, Principia, PsiOxus, Puma, Regeneron, Rigontec, Roche, Sanofi, Sierra Oncology, Synthon, Taiho, Takeda, Tesaro, Transgene, Turning Point Therapeutics and Upshersmith. E.C. declares advisory board participation (financial interest) from Adcendo, Alkermes, Amunix, Anaveon, Amcure, Astra Zeneca, BMS, Janssen, MonTa, MSD, Nanobiotix, Nouscom, Novartis, OncoDNA, PharmaMar, Roche/Genentech, Sanofi, Servier, SyneosHealth, TargImmune and T-knife, and they received research grants from Achilles and BeiGene, with steering committee (financial interest) and IDMC participation from BeiGene (IDMC steering committee), the EORTC IDMC chair (non-financial interest), MedSIR (steering committee), and Novartis (steering committee), as well as honoraria for scientific board participation from Adcendo, Chugai Pharmaceuticals, and PsiOxus Therapeutics, along with their employment at the HM Hospitals Group and the START Program of Early Phase Clinical Drug Development in Oncology, the START corporation, Oncoart Associated, and International Cancer Consultants. Further, E.C. is the founder and president of the not-for-profit Investigational Therapeutics in Oncological Sciences (INTHEOS) Foundation. A.O. has been a consultant for Servier, WWCT, viewpoint, Daichii-Sankyo, CancerAppy and Entrechem, and they are a former employee of Symphogen and Servier and are currently a consultant at NMS. The authors have no conflicts of interest to declare in relation to this manuscript.

References

  1. Kraehenbuehl, L.; Weng, C.H.; Eghbali, S.; Wolchok, J.D.; Merghoub, T. Enhancing immunotherapy in cancer by targeting emerging immunomodulatory pathways. Nat. Rev. Clin. Oncol. 2021, 19, 37–50. [Google Scholar] [CrossRef] [PubMed]
  2. Sharma, P.; Allison, J.P. Immune checkpoint therapy: Forging ahead. Sci. Transl. Med. 2022, 14, eadf2947. [Google Scholar] [CrossRef] [PubMed]
  3. Tawbi, H.A.; Schadendorf, D.; Lipson, E.J.; Ascierto, P.A.; Matamala, L.; Gutiérrez, E.C.; Rutkowski, P.; Gogas, H.J.; Lao, C.D.; De Menezes, J.J.; et al. Relatlimab and Nivolumab versus Nivolumab in Untreated Advanced Melanoma. N. Engl. J. Med. 2022, 386, 24–34. [Google Scholar] [CrossRef] [PubMed]
  4. Ribas, A.; Wolchok, J.D. Cancer Immunotherapy Using Checkpoint Blockade. Science 2018, 359, 1350–1355. [Google Scholar] [CrossRef] [PubMed]
  5. Cho, B.C.; Abreu, D.R.; Hussein, M.; Cobo, M.; Patel, A.J.; Secen, N.; Lee, K.H.; Massuti, B.; Hiret, S.; Yang, J.C.H.; et al. Tiragolumab plus atezolizumab versus placebo plus atezolizumab as a first-line treatment for PD-L1-selected non-small-cell lung cancer (CITYSCAPE): Primary and follow-up analyses of a randomised, double-blind, phase 2 study. Lancet Oncol. 2022, 23, 781–792. [Google Scholar] [CrossRef] [PubMed]
  6. Herbst, R.S.; Majem, M.; Barlesi, F.; Carcereny, E.; Chu, Q.; Monnet, I.; Sanchez-Hernandez, A.; Dakhil, S.; Camidge, D.R.; Winzer, L.; et al. COAST: An Open-Label, Phase II, Multidrug Platform Study of Durvalumab Alone or in Combination with Oleclumab or Monalizumab in Patients with Unresectable, Stage III Non–Small-Cell Lung Cancer. J. Clin. Oncol. 2022, 40, 3383–3393. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, L.; Sfakianos, J.P.; Beaumont, K.G.; Akturk, G.; Horowitz, A.; Sebra, R.P.; Farkas, A.M.; Gnjatic, S.; Hake, A.; Izadmehr, S.; et al. Myeloid cell–associated resistance to PD-1/PD-L1 blockade in urothelial cancer revealed through bulk and single-cell RNA sequencing. Clin. Cancer Res. 2021, 27, 4287–4300. [Google Scholar] [CrossRef] [PubMed]
  8. McRitchie, B.R.; Akkaya, B. Exhaust the exhausters: Targeting regulatory T cells in the tumor microenvironment. Front. Immunol. 2022, 13, 940052. [Google Scholar] [CrossRef] [PubMed]
  9. Siu, L.L.; Wang, D.; Hilton, J.; Geva, R.; Rasco, D.; Perets, R.; Abraham, A.K.; Wilson, D.C.; Markensohn, J.F.; Lunceford, J.; et al. First-in-Class Anti-immunoglobulin–like Transcript 4 Myeloid-Specific Antibody MK-4830 Abrogates a PD-1 Resistance Mechanism in Patients with Advanced Solid Tumors. Clin. Cancer Res. 2021, 28, 57–70. [Google Scholar] [CrossRef]
  10. Yunna, C.; Mengru, H.; Lei, W.; Weidong, C. Macrophage M1/M2 polarization. Eur. J. Pharmacol. 2020, 877, 173090. [Google Scholar] [CrossRef]
  11. Mills, C.D. M1 and M2 Macrophages: Oracles of Health and Disease. Crit. Rev. Immunol. 2012, 32, 463–488. [Google Scholar] [CrossRef] [PubMed]
  12. Chávez-Galán, L.; Olleros, M.L.; Vesin, D.; Garcia, I. Much More than M1 and M2 Macrophages, There are also CD169+ and TCR+ Macrophages. Front. Immunol. 2015, 6, 263. [Google Scholar] [CrossRef] [PubMed]
  13. Matsumoto, A.; Naito, M.; Itakura, H.; Ikemoto, S.; Asaoka, H.; Hayakawa, I.; Kanamori, H.; Aburatani, H.; Takaku, F.; Suzuki, H. Human macrophage scavenger receptors: Primary structure, expression, and localization in atherosclerotic lesions. Proc. Natl. Acad. Sci. USA 1990, 87, 9133–9137. [Google Scholar] [CrossRef] [PubMed]
  14. Goldstein, J.L.; Ho, Y.K.; Basu, S.K.; Brown, M.S. Binding site on macrophages that mediates uptake and degradation of acetylated low density lipoprotein, producing massive cholesterol deposition. Proc. Natl. Acad. Sci. USA 1979, 76, 333–337. [Google Scholar] [CrossRef] [PubMed]
  15. Kunjathoor, V.V.; Febbraio, M.; Podrez, E.A.; Moore, K.J.; Andersson, L.; Koehn, S.; Rhee, J.S.; Silverstein, R.; Hoff, H.F.; Freeman, M.W. Scavenger receptors class A-I/II and CD36 are the principal receptors responsible for the uptake of modified low density lipoprotein leading to lipid loading in macrophages. J. Biol. Chem. 2002, 277, 49982–49988. [Google Scholar] [CrossRef] [PubMed]
  16. Brown, M.S.; Goldstein, J.L.; Krieger, M.; Ho, Y.K.; Anderson, R.G. Reversible accumulation of cholesteryl esters in macrophages incubated with acetylated lipoproteins. J. Cell Biol. 1979, 82, 597–613. [Google Scholar] [CrossRef]
  17. Zani, I.A.; Stephen, S.L.; Mughal, N.A.; Russell, D.; Homer-Vanniasinkam, S.; Wheatcroft, S.B.; Ponnambalam, S. Scavenger Receptor Structure and Function in Health and Disease. Cells 2015, 4, 178–201. [Google Scholar] [CrossRef]
  18. Canton, J.; Neculai, D.; Grinstein, S. Scavenger receptors in homeostasis and immunity. Nat. Rev. Immunol. 2013, 13, 621–634. [Google Scholar] [CrossRef] [PubMed]
  19. Gudgeon, J.; Marín-Rubio, J.L.; Trost, M. The role of macrophage scavenger receptor 1 (MSR1) in inflammatory disorders and cancer. Front. Immunol. 2022, 13, 1012002. [Google Scholar] [CrossRef]
  20. Ji, Q.; Huang, K.; Jiang, Y.; Lei, K.; Tu, Z.; Luo, H.; Zhu, X. Comprehensive analysis of the prognostic and role in immune cell infiltration of MSR1 expression in lower-grade gliomas. Cancer Med. 2022, 11, 2020–2035. [Google Scholar] [CrossRef]
  21. Zagorulya, M.; Duong, E.; Spranger, S. Impact of anatomic site on antigen-presenting cells in cancer. J. Immunother. Cancer 2020, 8, e001204. [Google Scholar] [CrossRef]
  22. Hirayama, S.; Ishii, G.; Nagai, K.; Ono, S.; Kojima, M.; Yamauchi, C.; Aokage, K.; Hishida, T.; Yoshida, J.; Suzuki, K.; et al. Prognostic Impact of CD204-Positive Macrophages in Lung Squamous Cell Carcinoma: Possible Contribution of Cd204-Positive Macrophages to the Tumor-Promoting Microenvironment. J. Thorac. Oncol. 2012, 7, 1790–1797. [Google Scholar] [CrossRef]
  23. Ohtaki, Y.; Ishii, G.; Nagai, K.; Ashimine, S.; Kuwata, T.; Hishida, T.; Nishimura, M.; Yoshida, J.; Takeyoshi, I.; Ochiai, A. Stromal macrophage expressing CD204 is associated with tumor aggressiveness in lung adenocarcinoma. J. Thorac. Oncol. 2010, 5, 1507–1515. [Google Scholar] [CrossRef]
  24. Li, Z.; Maeda, D.; Yoshida, M.; Umakoshi, M.; Nanjo, H.; Shiraishi, K.; Saito, M.; Kohno, T.; Konno, H.; Saito, H.; et al. The intratumoral distribution influences the prognostic impact of CD68- and CD204-positive macrophages in non-small cell lung cancer. Lung Cancer 2018, 123, 127–135. [Google Scholar] [CrossRef]
  25. Sun, Y.; Xu, S. Tumor-Associated CD204-Positive Macrophage Is a Prognostic Marker in Clinical Stage I Lung Adenocarcinoma. BioMed Res. Int. 2018, 2018, 8459193. [Google Scholar] [CrossRef]
  26. Kawachi, A.; Yoshida, H.; Kitano, S.; Ino, Y.; Kato, T.; Hiraoka, N. Tumor-associated CD204+ M2 macrophages are unfavorable prognostic indicators in uterine cervical adenocarcinoma. Cancer Sci. 2018, 109, 863–870. [Google Scholar] [CrossRef]
  27. Miyasato, Y.; Shiota, T.; Ohnishi, K.; Pan, C.; Yano, H.; Horlad, H.; Yamamoto, Y.; Yamamoto-Ibusuki, M.; Iwase, H.; Takeya, M.; et al. High density of CD204-positive macrophages predicts worse clinical prognosis in patients with breast cancer. Cancer Sci. 2017, 108, 1693–1700. [Google Scholar] [CrossRef]
  28. Yuan, Y.; Zhao, Q.; Zhao, S.; Zhang, P.; Zhao, H.; Li, Z.; Du, Y.; Tian, X.; Lu, J. Characterization of transcriptome profile and clinical features of a novel immunotherapy target CD204 in diffuse glioma. Cancer Med. 2019, 8, 3811–3821. [Google Scholar] [CrossRef]
  29. Ikarashi, D.; Kitano, S.; Tsuyukubo, T.; Takenouchi, K.; Nakayama, T.; Onagi, H.; Sakaguchi, A.; Yamashita, M.; Mizugaki, H.; Maekawa, S.; et al. Pretreatment tumour immune microenvironment predicts clinical response and prognosis of muscle-invasive bladder cancer in the neoadjuvant chemotherapy setting. Br. J. Cancer 2022, 126, 606–614. [Google Scholar] [CrossRef]
  30. La Fleur, L.; Boura, V.F.; Alexeyenko, A.; Berglund, A.; Pontén, V.; Mattsson, J.S.; Djureinovic, D.; Persson, J.; Brunnström, H.; Isaksson, J.; et al. Expression of scavenger receptor MARCO defines a targetable tumor-associated macrophage subset in non-small cell lung cancer. Int. J. Cancer 2018, 143, 1741–1752. [Google Scholar] [CrossRef]
  31. Tang, Z.; Kang, B.; Li, C.; Chen, T.; Zhang, Z. GEPIA2: An enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019, 47, W556–W560. [Google Scholar] [CrossRef] [PubMed]
  32. Park, S.J.; Yoon, B.H.; Kim, S.K.; Kim, S.Y. GENT2: An updated gene expression database for normal and tumor tissues. BMC Med. Genom. 2019, 12 (Suppl. 5), 101. [Google Scholar] [CrossRef] [PubMed]
  33. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef]
  34. Li, T.; Fu, J.; Zeng, Z.; Cohen, D.; Li, J.; Chen, Q.; Li, B.; Liu, X.S. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020, 48, W509–W514. [Google Scholar] [CrossRef] [PubMed]
  35. Li, T.; Fan, J.; Wang, B.; Traugh, N.; Chen, Q.; Liu, J.S.; Li, B.; Liu, X.S. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017, 77, e108–e110. [Google Scholar] [CrossRef] [PubMed]
  36. Mulder, K.; Patel, A.A.; Kong, W.T.; Piot, C.; Halitzki, E.; Dunsmore, G.; Khalilnezhad, S.; Irac, S.E.; Dubuisson, A.; Chevrier, M.; et al. Cross-tissue single-cell landscape of human monocytes and macrophages in health and disease. Immunity 2021, 54, 1883–1900.e5. [Google Scholar] [CrossRef]
  37. The Cancer Genome Atlas Research Network; Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.M.; Ozenberger, B.A.; Ellrott, K.; Shmulevich, I.; Sander, C.; Stuart, J.M. The cancer genome atlas pan-cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef] [PubMed]
  38. Kovács, S.A.; Győrffy, B. Transcriptomic datasets of cancer patients treated with immune-checkpoint inhibitors: A systematic review. J. Transl. Med. 2022, 20, 249. [Google Scholar] [CrossRef] [PubMed]
  39. Lánczky, A.; Győrffy, B. Web-Based Survival Analysis Tool Tailored for Medical Research (KMplot): Development and Implementation. J. Med. Internet Res. 2021, 23, e27633. [Google Scholar] [CrossRef]
  40. Győrffy, B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput. Struct. Biotechnol. J. 2021, 19, 4101–4109. [Google Scholar] [CrossRef]
  41. Fekete, J.T.; Győrffy, B. ROCplot.org: Validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3,104 breast cancer patients. Int. J. Cancer 2019, 145, 3140–3151. [Google Scholar] [CrossRef] [PubMed]
Figure 1. MSR1 expression profiles across all tumor samples and paired normal tissues. (A) MSR1 expression as transcript per million (TPM) in several tumor types represented using a heatmap and bar graph using the GEPIA2 database. Statistically significant differences between the normal tissue and tumoral ones are marked by an asterisk. (B) Heatmap showing the fold changes (FCs) between the expression of MSR1 in the tumoral and normal tissues for the previous tumor types. The inclusion criteria for the selection of cancers for further analysis were set as a FC of >2 and a statistical difference between the expression in the tumor and normal tissues.
Figure 1. MSR1 expression profiles across all tumor samples and paired normal tissues. (A) MSR1 expression as transcript per million (TPM) in several tumor types represented using a heatmap and bar graph using the GEPIA2 database. Statistically significant differences between the normal tissue and tumoral ones are marked by an asterisk. (B) Heatmap showing the fold changes (FCs) between the expression of MSR1 in the tumoral and normal tissues for the previous tumor types. The inclusion criteria for the selection of cancers for further analysis were set as a FC of >2 and a statistical difference between the expression in the tumor and normal tissues.
Ijms 25 03987 g001
Figure 2. Association of MSR1 expression levels with immune infiltrates in the ten selected cancer types. (A) Correlation analysis between MSR1 expression and relevant immune population infiltrates in different tumor types using TIMER2.0. A Spearman’s correlation was used with a purity adjustment. (B) Dot plot details of the correlations displayed in (A): macrophages, DCs and neutrophils. The MSR1 expression level is represented on the Y axis as log2 (transcripts per million, TPM) and the population infiltration level is shown on the X axis.
Figure 2. Association of MSR1 expression levels with immune infiltrates in the ten selected cancer types. (A) Correlation analysis between MSR1 expression and relevant immune population infiltrates in different tumor types using TIMER2.0. A Spearman’s correlation was used with a purity adjustment. (B) Dot plot details of the correlations displayed in (A): macrophages, DCs and neutrophils. The MSR1 expression level is represented on the Y axis as log2 (transcripts per million, TPM) and the population infiltration level is shown on the X axis.
Ijms 25 03987 g002
Figure 3. Association of MSR1 expression levels with macrophage subtypes (M0, M1 and M2 classifications) and co-stimulatory immune checkpoints for the ten selected cancer types. (A) Correlation analysis between MSR1 expression and macrophage subtypes (M2, M1 and M0). (B) Correlation between the expression of MSR1 and immunomodulatory proteins. (C) Single cell correlation analysis between MSR1 and the co-inhibitory genes described in (B) for those tumors with available data.
Figure 3. Association of MSR1 expression levels with macrophage subtypes (M0, M1 and M2 classifications) and co-stimulatory immune checkpoints for the ten selected cancer types. (A) Correlation analysis between MSR1 expression and macrophage subtypes (M2, M1 and M0). (B) Correlation between the expression of MSR1 and immunomodulatory proteins. (C) Single cell correlation analysis between MSR1 and the co-inhibitory genes described in (B) for those tumors with available data.
Ijms 25 03987 g003
Figure 4. Association of MSR1 expression levels with overall survival (OS) in patients treated with immune check-point inhibitors (ICI), including anti-PD1, anti-PDL1 and anti-CTLA4, respectively. (A) Kaplan–Meier survival plots comparing the high vs. low expression levels of MSR1, including the patients with all available tumors. (B) Kaplan–Meier survival plots comparing the high vs. low expression levels of MSR1 for the SKCM patients. (C) Kaplan–Meier survival plots comparing the high vs. low expression levels of CD274 with OS, including the patients with all available tumors.
Figure 4. Association of MSR1 expression levels with overall survival (OS) in patients treated with immune check-point inhibitors (ICI), including anti-PD1, anti-PDL1 and anti-CTLA4, respectively. (A) Kaplan–Meier survival plots comparing the high vs. low expression levels of MSR1, including the patients with all available tumors. (B) Kaplan–Meier survival plots comparing the high vs. low expression levels of MSR1 for the SKCM patients. (C) Kaplan–Meier survival plots comparing the high vs. low expression levels of CD274 with OS, including the patients with all available tumors.
Ijms 25 03987 g004
Figure 5. Association of MSR1 expression levels with responses in patients treated with ICIs including anti-PD1, anti-PDL1 and anti-CTLA4 (measured using area under the curve (AUC) and p values). The data were obtained using ROCPLOT. (A) For all available cancers. (B) For the SKCM patients.
Figure 5. Association of MSR1 expression levels with responses in patients treated with ICIs including anti-PD1, anti-PDL1 and anti-CTLA4 (measured using area under the curve (AUC) and p values). The data were obtained using ROCPLOT. (A) For all available cancers. (B) For the SKCM patients.
Ijms 25 03987 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sanvicente, A.; Díaz-Tejeiro, C.; Nieto-Jiménez, C.; Paniagua-Herranz, L.; López Cade, I.; Balázs, G.; Moreno, V.; Pérez-Segura, P.; Calvo, E.; Ocaña, A. In Silico Transcriptomic Expression of MSR1 in Solid Tumors Is Associated with Responses to Anti-PD1 and Anti-CTLA4 Therapies. Int. J. Mol. Sci. 2024, 25, 3987. https://doi.org/10.3390/ijms25073987

AMA Style

Sanvicente A, Díaz-Tejeiro C, Nieto-Jiménez C, Paniagua-Herranz L, López Cade I, Balázs G, Moreno V, Pérez-Segura P, Calvo E, Ocaña A. In Silico Transcriptomic Expression of MSR1 in Solid Tumors Is Associated with Responses to Anti-PD1 and Anti-CTLA4 Therapies. International Journal of Molecular Sciences. 2024; 25(7):3987. https://doi.org/10.3390/ijms25073987

Chicago/Turabian Style

Sanvicente, Adrián, Cristina Díaz-Tejeiro, Cristina Nieto-Jiménez, Lucia Paniagua-Herranz, Igor López Cade, Győrffy Balázs, Víctor Moreno, Pedro Pérez-Segura, Emiliano Calvo, and Alberto Ocaña. 2024. "In Silico Transcriptomic Expression of MSR1 in Solid Tumors Is Associated with Responses to Anti-PD1 and Anti-CTLA4 Therapies" International Journal of Molecular Sciences 25, no. 7: 3987. https://doi.org/10.3390/ijms25073987

APA Style

Sanvicente, A., Díaz-Tejeiro, C., Nieto-Jiménez, C., Paniagua-Herranz, L., López Cade, I., Balázs, G., Moreno, V., Pérez-Segura, P., Calvo, E., & Ocaña, A. (2024). In Silico Transcriptomic Expression of MSR1 in Solid Tumors Is Associated with Responses to Anti-PD1 and Anti-CTLA4 Therapies. International Journal of Molecular Sciences, 25(7), 3987. https://doi.org/10.3390/ijms25073987

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop