Next Article in Journal
Chronological and Biological Aging in Amyotrophic Lateral Sclerosis and the Potential of Senolytic Therapies
Previous Article in Journal
Possibility of Cell Block Specimens from Overnight-Stored Bile for Next-Generation Sequencing of Cholangiocarcinoma
Previous Article in Special Issue
Early Pharmacodynamic Changes Measured Using RNA Sequencing of Peripheral Blood from Patients in a Phase I Study with Mitazalimab, a Potent CD40 Agonistic Monoclonal Antibody
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes

1
Department of Immunotechnology, Lund University, 223 81 Lund, Sweden
2
CREATE Health Cancer Center, Lund University, 223 81 Lund, Sweden
3
Division of Clinical Genetics, Department of Laboratory Medicine, Lund University, 221 84 Lund, Sweden
4
Division of Oncology, Department of Clinical Sciences Lund, Lund University, 221 84 Lund, Sweden
5
Department of Translational Medicine, Lund University, 205 02 Malmö, Sweden
6
Department of Urology, Skåne University Hospital, 205 02 Malmö, Sweden
7
Alligator Bioscience AB, Medicon Village, 223 63 Lund, Sweden
8
Department of Clinical Genetics, University and Regional Laboratories Region Skåne, 221 85 Lund, Sweden
*
Author to whom correspondence should be addressed.
Cells 2024, 13(11), 926; https://doi.org/10.3390/cells13110926
Submission received: 11 April 2024 / Revised: 16 May 2024 / Accepted: 20 May 2024 / Published: 28 May 2024

Abstract

:
Bladder cancer is a heterogenous disease, and molecular subtyping is a promising method to capture this variability. Currently, the immune compartment in relation to subtypes is poorly characterized. Here, we analyzed the immune compartment in bladder tumors and normal bladder urothelium with a focus on T cell subpopulations using flow cytometry and RNA sequencing. The results were investigated in relation to tumor invasiveness (NMIBC/MIBC) and molecular subtypes according to the Lund Taxonomy system. Whereas the NMIBC/MIBC differed in the overall immune infiltration only, the molecular subtypes differed both in terms of immune infiltration and immune compartment compositions. The Basal/Squamous (Ba/Sq) and genomically unstable (GU) tumors displayed increased immune infiltration compared to urothelial-like (Uro) tumors. Additionally, the GU tumors had a higher proportion of regulatory T cells within the immune compartment compared to Uro tumors. Furthermore, sequencing showed higher levels of exhaustion in CD8+ T cells from GU tumors compared to both Uro tumors and the control. Although no such difference was detected at the transcriptomic level in Uro tumors compared to the controls, CD8+ T cells in Uro tumors showed higher expression of several exhaustion markers at the protein level. Taken together, our findings indicate that depending on the molecular subtype, different immunotherapeutic interventions might be warranted.

1. Introduction

In 2020, 573,000 patients worldwide were diagnosed with bladder cancer, and 213,000 succumbed to the disease, thus making bladder cancer the second most prevalent urogenital malignant disease globally [1]. Bladder cancer is traditionally divided into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC). Around 25% of bladder cancer patients have muscle-invasive disease at diagnosis, and despite drastic treatment, including neoadjuvant chemotherapy plus radical cystectomy, and sometimes also adjuvant immunotherapy, the 5-year survival rate is only around 50% [2,3]. Patients diagnosed with NMIBC have a better 5-year survival rate, but approximately every second patient suffers from at least one local recurrence, and, in the long term, there is a risk of progression to MIBC [4]. The high recurrence rate in NMIBC urges frequent monitoring and multiple treatment occasions, making bladder cancer one of the most expensive cancers per patient for society [5]. Treatment of bladder cancer has largely remained the same during the last decades, and survival has not improved to the same extent as in many other cancer forms [2].
One reason for the lack of progress is the extensive heterogeneity of bladder cancer and the inability of the traditional classification to capture this, leading to low-resolution patient stratification. Molecular classification systems based on gene expression profiling have been developed to account for disease heterogeneity and thereby better stratify patients. The Lund Taxonomy molecular classification system is solely based on cancer cells, and it is the only classification system that is applicable to both NMIBC and MIBC [6]. The molecular subtypes include, e.g., the urothelial-like (Uro), Basal/Squamous (Ba/Sq), and genomically unstable (GU) subtypes, and these have been shown to differentially respond to therapy [7]. For example, the response rate to immunotherapy has been shown to be around 50% for the GU subtype, whereas the Uro and Ba/Sq subtypes had an overall response rate of around 20% [8].
It is additionally well-established that the immune microenvironment has a major impact on bladder cancer. For example, CD8+ T cells in the tumor have been shown to correlate with overall survival in MIBC [9]. In contrast, regulatory T cells (Tregs) are generally known as tumor-promoting cells in cancers, including bladder cancer [10]. However, there are also studies showing that infiltration of Tregs correlates with improved survival of bladder cancer patients [11,12]. The infiltration of specific immune cells also impacts the response to immunotherapy, and several studies have, for example, demonstrated that the presence of CD8+ T cells in tumors correlate with the response [8,13]. Despite some efforts, the immune microenvironment of the different molecular subtypes is still poorly understood. So far, investigations have mainly relied on bulk RNA sequencing data and tissue staining using a limited set of markers in parallel [8,14,15]. Detailed studies where distinct immune populations are precisely defined in relation to specific cancer cell molecular subtypes are lacking.
In this study, we investigated the immune compartment in bladder tumors stratified according to invasiveness and molecular subtypes using the Lund Taxonomy system [7]. We assessed the infiltration of CD45+ immune cells, as well as the relative abundance of helper T cells (Th), Tregs, and cytotoxic CD8+ T cells. Additionally, we assessed transcriptomic differences of CD8+ T cells in GU and Uro tumors. We show that the molecular subtypes have different levels of immune infiltration as well as different cellular composition of the immune compartment. These findings indicate that bladder cancer molecular subtypes shape characteristic immune microenvironments. Dissecting and understanding these different environments could aid patient stratification and treatment choice, especially regarding immunotherapies.

2. Materials and Methods

2.1. Sample Preparation and Classification

Tumor tissue samples were cold-cup biopsies obtained from the exophytic part of the tumor during transurethral resection (TURB) of treatment-naïve patients with primary bladder cancer. Tumors were classified according to the Lund Taxonomy based on IHC and/or bulk RNA seq data as part of the UROSCANSEQ prospective sequencing initiative. Non-malignant control bladder tissues were obtained during open partial cystectomy for benign conditions or non-urothelial tumors (patient characteristics presented in Table 1). Immediately after surgery, biopsies were placed in tissue storage solution (Miltenyi Biotec, Somerville, MA, USA) and transported to the laboratory on ice. Subsequently, tissue samples were cut into small pieces in RPMI 1640 medium (Thermo Fisher Scientific, Bremen, Germany) supplemented with 0.1 mg/mL gentamycin (Sigma Aldrich, St Louis, MO, USA) and prepared into single-cell suspensions through enzymatic digestion at 37 °C for 20 min using Collagenase IV (Sigma Aldrich) (2 mg/mL) and DNase I (Sigma Aldrich) (200 Kunitz units/mL), followed by filtering through 70 µm cell strainer (BD Biosciences, San Jose, CA, USA). The study was approved by the Research Ethics Board of Lund University (Dnr 2017/269 and 2018/963) and Stockholm University (Dnr 2020-05559 and 2022-03081-02).

2.2. Flow Cytometry, Sorting, RNA Extraction, and Sequencing

Cells were stained with Fixable Viability Stain (FVS) 620 (BD Biosciences) according to manufacturer’s protocol to enable viability evaluation. Cells were then washed, blocked using mouse IgG (Jackson ImmunoResearch, West Grove, PA), and finally stained with an antibody panel for 20 min at +4 °C in Brilliant Stain Buffer (BD Biosciences) (Table A1). Stained cells were acquired on a BD FACS Aria II (BD Biosciences) for analysis and sorting of CD8+ T cells, Th cells, and Tregs. For the phenotypic analysis, doublets were excluded, and viable T cells (FVS CD45+CD3+) were gated. CD8+ T cells, Tregs, and Th cells were subsequently identified as CD8+CD4, CD4+CD25+CD127−/dim, and CD4+CD25, respectively (Figure 1A). For sorting, cells expressing lineage (Lin) markers (CD14, CD16, CD19, CD20, CD56, and CD66b) were excluded before gating CD8+ T cells, as described above. T cell populations were sorted from tumors as well as from non-malignant control tissues from patients undergoing bladder surgery for reasons other than bladder cancer (purity > 95%). Cells were sorted directly into extraction buffer from the Arcturus PicoPure RNA Isolation Kit (Thermo Fisher Scientific). RNA sequencing was performed as previously described [16].

2.3. Bulk RNA Seq Analysis of Sorted CD8+ T Cells

RNA sequencing reads were aligned to the human reference genome hg19 using STAR 2.5.0a [17]. Reads within genes were counted using featureCounts (version 1.6.3) [18]. In brief, a gene was considered detected at read counts >5, and samples with <7000 detected genes were excluded from the analysis. Genes with >1 counts per million reads mapped (CPM) in 2 or more samples were kept for further analysis using EdgeR (version 3.38.2) [19]. Data normalization and differential gene expression analysis were carried out using DESeq2 (version 1.36.0) with the apeglm shrinkage estimator [20,21]. Differentially expressed genes, defined as adjusted p-values (padj) < 0.05 and absolute log2 fold change >1, were extracted. Gene set enrichment analysis (GSEA) is a method to investigate the relative enrichment of a set of genes between groups. Here, we used GSEA to compare T cell exhaustion between the molecular subtypes. The fgsea R package (version 1.22.0) with genes ranked according to the log2 fold changes obtained from the differential expression analysis was used to perform GSEA [22]. Enrichment analysis was performed with the differentially expressed genes (padj < 0.05 and absolute log2 fold change >1) using the enrichR R package (version 3.1) together with the Reactome pathways database [23]. For GSEA/enrichment analysis, results with p-value/padj < 0.05 were considered significant.

2.4. Statistical Analysis of Flow Cytometry Data

For the statistical analysis of the flow cytometry data, non-parametric Kruskal–Wallis analysis and the two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test or the Mann–Whitney U test were used as indicated. p-value/padj < 0.05 was considered significant.

3. Results

3.1. Invasiveness of Bladder Cancer Is Associated with Higher Immune Infiltration but Poorly Explains Tumor Immune Microenvironment Heterogeneity

We investigated the overall immune infiltration, as well as the proportion of T cells, CD8+ T cells, Th, and Tregs out of CD45+ cells in tumor biopsies from patients with bladder cancer stratified according to NMIBC/MIBC, as well as in control bladder tissue from patients without bladder cancer (gating strategy in Figure 1A). As expected, immune infiltration was significantly higher in both NMIBC and MIBC compared to control tissue (Figure 1B). In line with previous findings, muscle-invasive tumors were associated with higher immune infiltration (Figure 1B) [15,24]. However, no difference in the immune cell composition was detected between NMIBC and MIBC, as the proportion of CD3+, CD8+, CD4+, Th, or Tregs out of all CD45+ cells was not significantly different (Figure 1C–G). Stratifying according to tumor stage yielded similar results (Figure A1). At the same time, a large variation was seen among the NMIBC and MIBC samples, respectively. Furthermore, we did not observe any differences in immune composition in relation to bacteriuria status prior to surgery or in relation to sex (Figure A1 and Figure A2). Taken together, invasiveness was associated with overall immune infiltration in bladder cancer, but the NMIBC and MIBC dichotomy could not explain the heterogeneity within the tumor immune compartment.
Figure 1. Gating strategy for analysis of Tregs, Th cells, and CD8+ T cells (A). Immune cell infiltration measured as CD45+ out of all viable cells (B). Frequency of T cells and T cell subsets out of all CD45+ cells (CG). Treg = regulatory T cells, Th = T helper, NMIBC = non-muscle-invasive bladder cancer, MIBC = muscle-invasive bladder cancer. Statistical significance was determined by the Kruskal–Wallis test followed by the two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Figure 1. Gating strategy for analysis of Tregs, Th cells, and CD8+ T cells (A). Immune cell infiltration measured as CD45+ out of all viable cells (B). Frequency of T cells and T cell subsets out of all CD45+ cells (CG). Treg = regulatory T cells, Th = T helper, NMIBC = non-muscle-invasive bladder cancer, MIBC = muscle-invasive bladder cancer. Statistical significance was determined by the Kruskal–Wallis test followed by the two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Cells 13 00926 g001

3.2. Molecular Subtypes of Bladder Cancer Harbor Distinct Immune Microenvironments

Response to immunotherapy varies between molecular subtypes, and to better understand the immune environment underpinning these differences, we stratified tumor samples according to the Lund Taxonomy, which depicts the cancer cell phenotype independent of infiltration, and we analyzed the immune infiltration and abundance of T cell subsets. Immune infiltration was significantly higher in the Ba/Sq and GU subtypes compared to the Uro subtype (Figure 2A), and the same trend was seen among MIBC samples only (padj = 0.0608 GU vs. Uro) (Figure 2B). T cell infiltration (CD3+ cells out of viable cells) was also higher in the Ba/Sq and GU subtypes compared to the Uro subtype (Figure A3). In contrast, no statistically significant difference in the proportions of T cells or CD8+ T cells within the immune compartment were detected between the molecular subtypes (Figure 2C,D). Interestingly, the proportion of CD4+ T cells within the immune compartment, as well as the CD4/CD8 ratio, were higher in the GU and Uro subtypes compared to the Ba/Sq subtype (Figure 2E,F). Moreover, upon further analysis of the CD4+ T cells, it was shown that Tregs, but not Th cells, were significantly enriched within the immune compartment of the GU subtype compared to the Uro subtype (Figure 2G,H). Similarly, the proportion of Tregs out of all cells was higher in the GU subtype compared to the Uro subtype (Figure A3). Taken together, although the immune infiltration and the proportions of immune cells varied within each molecular subtype, statistically significant differences were detected.

3.3. Exhaustion Gene Profile Enriched in CD8+ T Cells from GU Tumors

To further investigate CD8+ T cells from the different molecular subtypes, we sorted and RNA sequenced CD8+ T cells from bladder tumor biopsies and control bladder tissues. After filtering out low-quality samples, adequate sample numbers were obtained for CD8+ T cells from control tissue (n = 4), GU (n = 4) and Uro tumors (n = 4). Principal component analysis revealed that the CD8+ T cell populations grouped according to control/cancer and, to some degree, according to molecular subtype (Figure 3A). Differential gene expression analysis resulted in 35 genes expressed at significantly different levels between GU and Uro tumors, several of which were related to T cell exhaustion, such as LAYN, CXCL13, and HAVCR2 (Figure 3B). To further investigate T cell exhaustion, GSEA using an exhaustion profile was performed [25]. The results revealed that CD8+ T cells from GU tumors were significantly enriched for the exhaustion profile when compared to CD8+ T cells from Uro tumors and control tissue, whereas no difference was observed between CD8+ T cells between Uro tumors and control tissue (Figure 3C). Furthermore, several genes were found to be differentially expressed by CD8+ T cells from GU and Uro tumors compared to the control tissue, and pathway analysis revealed that these genes were mostly related to cell cycle and immune signaling (Figure 3D,E and Figure A4). Interestingly, CD8+ T cells in GU tumors showed higher expression of the gene VDR, encoding the vitamin D receptor, previously associated with reduced exhaustion [26], compared to both Uro tumors and controls (Figure 3D). For the Uro subtype, we further investigated the presence of exhaustion markers on the cell surface. Contrary to the findings on the transcriptomic level, we observed a significant increase in CTLA-4, LAG-3, PD-1, and TIM-3-positive CD8+ T cells in the Uro tumors compared to the controls. Sample numbers for the other subtypes were not adequate to perform statistical analysis, as the assessment was performed on fresh samples at a timepoint when the subtype was unknown (Figure 3F).

4. Discussion

In this study, we investigated the immune microenvironment in bladder tumors using flow cytometry and RNA sequencing. We show that the level of immune infiltration differs between NMIBC and MIBC as well as between molecular subtypes. Importantly, the immune microenvironment composition was significantly different between the molecular subtypes, with lower CD4+ T cell frequencies in Ba/Sq and GU tumors having a higher proportion of Tregs than Uro tumors. These differences could impact the functional response to cancer therapies, especially immunotherapy, and they should be considered when designing future clinical trials.
As bladder cancer has traditionally been divided according to the invasiveness of the disease, and because this stratification largely decides whether organ-sparing or radical treatment is suitable, we initially stratified the samples accordingly. We found that immune infiltration varied greatly between patients and that muscle-invasive disease was associated with higher overall immune infiltration, as demonstrated previously using DNA methylation analysis and tissue microarrays (TMAs) [15,24,27]. However, the observed variation in the proportion of T cell subsets (out of CD45+ cells) could not be explained by the MIBC/NMIBC dichotomy. This suggests that while the amount of immune infiltration increased with increasing tumor stage, there are changes in the composition of the immune compartment that relate to underlying biology rather than invasion depth. Furthermore, incidence rate and survival in bladder cancer have been shown to differ between the sexes [28]. For example, increased mortality has been observed for women within the first two years after diagnosis [29,30]. Moreover, sex differences have been implicated in anti-tumor immune responses in lung cancer [31]. In our study, we did not observe any significant differences in immune infiltration or composition in relation to sex, thus suggesting that the reported differences in bladder cancer may be related to infiltrating myeloid cells [32] or have other or additional underlying causes. It should be noted, though, that the incidence rate for bladder cancer is higher among males; hence, there are more samples from males in the study.
Bladder cancer molecular subtyping is a stratification system with higher resolution, and we show that the median immune infiltration was highest in the Ba/Sq subtype, followed by the luminal subtypes GU and Uro, which is in line with previous findings [15,33,34]. While MIBC comprises tumors of the Uro, Ba/Sq, and GU subtypes, the Uro subtype is by far the most common subtype among NMIBCs. By restricting the analysis to only MIBC cases, we demonstrate that the increased immune infiltration in Ba/Sq and GU tumors compared to Uro tumors still persists, indicating that the difference in NMIBC/MIBC prevalence does not explain the observed difference. We also observed an increase in infiltration of T cells in Ba/Sq and GU tumors compared to Uro tumors and this could potentially also be true for other immune populations. The agreement with previous findings on immune infiltration based on larger patient cohorts and using bulk RNA seq [7] or tissue staining [15] supports that our sample collection accurately represents the different molecular subtypes.
By further disentangling the immune compartments, we showed that the proportion of CD4+ T cells out of all CD45+ cells is higher in GU and Uro tumors compared to Ba/Sq tumors and controls. A comparison of the CD4+/CD8+ T cell ratio between the subtypes showed differences between subtypes, with Ba/Sq having the lowest ratio. We furthermore showed that Tregs are more abundant (proportion out of all cells) in GU compared to Uro tumors, which is in line with previous findings based on staining of tissue microarrays [15]. Importantly, our approach allowed us to also investigate the relative abundance of Tregs within the immune compartment. We demonstrated that the proportion of Tregs out of CD45+ cells is higher in GU tumors compared to Uro tumors. Although it has previously been suggested that Tregs in bladder cancer, unlike many other solid cancers [35], can be associated with better prognosis [11], the immunosuppressive role of Tregs is well-established, and the correlation could be due to Treg infiltration being associated with the infiltration of cells that contribute to the anti-tumor immune response. Despite the association, the consensus is that Tregs should be inhibited/eliminated in order to promote anti-tumor immune responses. This is, in fact, a strategy currently being pursued in clinical trials [36]. Targeting Tregs might constitute an additional therapeutic axis to further improve the anti-tumor immune response in bladder cancer, and, given the relatively high Treg abundance, GU tumors could be particularly responsive.
For the CD8+ T cells, we observed no significant difference in frequency in CD45+ cells among the different subtypes, but CD8+ T cells from GU, Uro, and control non-malignant bladder tissue (from patients not suffering from bladder cancer) were further compared using sequencing. Enrichment analysis revealed an upregulation of cell-cycle-related pathways in CD8+ T cells from both GU and Uro tumors compared to the control, suggesting local proliferation of T cells as a potential source of tumor-infiltrating T cells. Interestingly, CD8+ T cells in GU tumors also showed a higher expression of the vitamin-D-receptor-encoding gene, and triggering this receptor with the active form of vitamin D was recently shown to inhibit the transcription of several exhaustion markers in cytotoxic T cells [26]. In the study, treatment of patients with non-small-cell lung cancer with Rocaltrol (active vitamin D3) decreased the expression of PD-1, TIM-3, and TIGIT and resulted in cytokine production associated with anti-tumor immunity. In addition to VDR, several exhaustion-related genes were upregulated in GU tumors compared to both Uro tumors and control tissues, whereas no significant difference was observed between Uro tumors and control tissues. In contrast, the exhaustion markers analyzed at the protein level using flow cytometry were indeed shown to be higher in Uro tumors compared to the control, suggesting that some degree of exhaustion of CD8+ T cells also occurs in the Uro subtype. This observed discrepancy in protein- versus transcriptomic-level expression could perhaps be explained by the large variation in expression of the investigated exhaustion markers within the Uro cohort and that the limited number of samples sequenced expressed low levels of these. Only one GU tumor was available for flow cytometry analysis, and no comparisons could thus be made of the protein expression of exhaustion markers. In contrast to our findings, Oh et al. observed no significant difference in the frequency of exhausted CD8+ T cells (their CD8ENTPD1 cluster formed upon analysis of single-cell RNA seq data) in bladder tumors compared to control bladder tissues [37]. However, their control was bladder tissue adjacent to the tumor, which, unlike the control tissue used in the present study, could be affected by the tumor, and this might explain why we observed a significant increase in CD8+ T cell exhaustion in bladder tumors compared to control non-malignant tissue. Taken together, our results suggest that CD8+ T cell exhaustion can be a phenomenon in both Uro and GU subtypes.
Given the differences in immune infiltration and composition between subtypes, different immunomodulatory interventions might be warranted for specific subtypes. Of note, there are now ongoing efforts to investigate the utility of stratifying patients according to molecular subtypes clinically (GUSTO trial: NIHR128103). More work in this area could reveal subtype-specific immune modulatory strategies, which could improve patient stratification and treatment selection.
The current study is not without limitations. No systematic exclusion of patient samples was made based on patient comorbidities, and any effects thereof cannot be ruled out. Also, it should be noted that the sample cohort comprised a higher proportion of MIBC compared to the incidence rate in the population. The reason for this is that the size of some tumors, predominantly NMIBC, was inadequate to allow for a detailed analysis of T cell subsets. Furthermore, the limited number of samples from subtypes other than the Uro subtype decreased the statistical precision in the comparisons between subtypes. Similarly, the number of control bladder tissue samples was limited, as these are obtained during the rare event of performing bladder surgery for reasons other than bladder cancer. This control tissue was selected to ensure that the control samples were truly not affected by bladder cancer, thus representing controls more similar to healthy bladder tissue compared to, e.g., tissue adjacent to the tumor. Thus, although a larger cohort is needed to validate our findings, the data related to the rare control material are considered valuable, and the results regarding the rarer molecular subtypes offer new and potentially useful information for these patient groups.

5. Conclusions

Molecular subtypes in bladder cancer differ in relation to immune infiltration and composition, suggesting that distinct biology and, thus, putative targets exist in the respective subtypes. In particular, targeting Tregs could provide an additional therapeutic avenue in the GU subtype, further boosting the effectiveness of immunotherapy. These results support the further implementation of molecular subtypes in clinical trials for established as well as explorative therapies to tailor treatments and achieve improved survival.

Author Contributions

Conceptualization, K.L., F.L., C.A.K.B. and T.F.; methodology, V.S., K.L., F.L., K.F.A., H.L., P.E. (Peter Ellmark), K.H. and M.W.; investigation, V.S., K.F.A. and H.L.; formal analysis, V.S., S.R., K.F.A., H.L., K.L., F.L., G.S. and P.E. (Pontus Eriksson); data curation, V.S., K.L., F.L., M.W. and P.E. (Pontus Eriksson); writing—original draft preparation, V.S., mainly supported by K.L. and S.R.; writing—review and editing, all authors; visualization, V.S. and S.R.; supervision, K.L. and F.L.; project administration, K.L.; funding acquisition, K.L., C.A.K.B., F.L. and T.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Cancera Foundation, Gunnar Nilssons Cancerstiftelse, the Mats Paulsson Foundation, the Stefan Paulsson Foundation, the Per-Erik och Ulla Schybergs stiftelse/Royal Physiographic Society of Lund, the ISREC Foundation (through a joint grant to the Swiss Cancer Center, Lausanne and the CREATE Health Cancer Center, Lund, from the Biltema foundation), the Faculty of Engineering at Lund University (LTH), the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 847583 and the Swedish Cancer Society (grant number CAN 2023/2807), the Swedish Research Council (2021-00859), the Lund Medical Faculty (ALF), Skåne University Hospital Research Funds, the Cancer Research Fund at Malmö General Hospital, the Hjelm Family Foundation for Medical Research, the Gösta Jönsson Research Foundation, the Foundation of Urological Research (Ove and Carin Carlsson bladder cancer donation), and the Hillevi Fries Research Foundation.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Research Ethics Board of Lund University (Dnr 2017/269 4 April 2017 and 2018/963 4 December 2018) and Stockholm University (Dnr 2020-05559 11 November 2020 and 2022-03081-02 6 August 2022).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

We would like to thank Marianne Rissler from the Division of Clinical Genetics, Department of Laboratory Medicine of Lund University for her help with RNA extraction and sequencing procedures.

Conflicts of Interest

Authors Karin Hägerbrand and Peter Ellmark were employed by the company Alligator Bioscience AB, Medicon Village. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Table A1. Antibodies used for flow cytometry analysis and sorting of T cell populations from tumor and control bladder tissues.
Table A1. Antibodies used for flow cytometry analysis and sorting of T cell populations from tumor and control bladder tissues.
Staining AgentCloneSupplier
BD Horizon™ Fixable Viability Stain 620-BD biosciences
FITC Mouse Anti-human CD45HI30BD biosciences
PerCP-Cy5.5 Mouse Anti-human CD19HIB19Biolegend (San Diego, CA, USA)
PerCP-Cy5.5 Mouse Anti-human CD202H7BD biosciences
PerCP-Cy5.5 Mouse Anti-Human CD56B159BD biosciences
PerCP-Cy5.5 Mouse Anti-Human CD66bG10F5BD biosciences
PerCP-Cy5.5 Mouse Anti-Human CD14MφP9BD biosciences
PerCP-Cy5.5 Mouse Anti-Human CD163G8BD biosciences
BV785 Mouse Anti-Human CD3SK7Biolegend
PE-Cy7 Mouse Anti-Human CD8RPA-T8BD biosciences
BV605 Mouse Anti-Human CD4RPA-T4BD biosciences
BV421 Mouse Anti-Human CD127HIL-7R-M21BD biosciences
APC Mouse Anti-Human CD25M-A251BD biosciences
PE Mouse Anti-Human CTLA4L3D10Biolegend
PerCp.Cy 5.5 Mouse Anti-Human TIGITA15153GBiolegend
PE-Cy7 Mouse Anti-Human PD-1EH12.2H7Biolegend
BV510 Mouse Anti-Human CD4RPA-T4Biolegend
BV711 Mouse Anti-Human TIM37D3BD biosciences
BV785 Mouse Anti-Human CD3SK7Biolegend
R718/APC-R700 Mouse Anti-Human LAG-3T47-530BD biosciences
APC-Cy7 Mouse Anti-Human CD8SK1BD biosciences
AF700 Mouse Anti-Human CD4RPA-T4BD biosciences
BV785 Mouse Anti-Human CD8RPA-T8Biolegend
BV510 Mouse Anti-Human CD3SK7Biolegend
Figure A1. Infiltration of CD45+ immune cells out of all viable cells and frequency of T cell populations out of CD45+ immune cells in tumor samples stratified according tumor stage (AD) and bacteriuria (EH).
Figure A1. Infiltration of CD45+ immune cells out of all viable cells and frequency of T cell populations out of CD45+ immune cells in tumor samples stratified according tumor stage (AD) and bacteriuria (EH).
Cells 13 00926 g0a1
Figure A2. Infiltration of CD45+ immune cells out of all viable cells and frequency of T cell populations out of CD45+ immune cells in tumor samples according to sex (AD,F,G). Ratio of CD4+/CD8+ T cells (E) Statistical significance was determined by Kruskal–Wallis test followed by two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Figure A2. Infiltration of CD45+ immune cells out of all viable cells and frequency of T cell populations out of CD45+ immune cells in tumor samples according to sex (AD,F,G). Ratio of CD4+/CD8+ T cells (E) Statistical significance was determined by Kruskal–Wallis test followed by two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Cells 13 00926 g0a2
Figure A3. T cell populations in tumor samples classified according to the Lund Taxonomy molecular subtypes. Percentage of T cells (A), CD8+ T cells (B), CD4+ T cells (C), Tregs (D), and Th (E) out of all viable cells. Statistical significance was determined by Kruskal–Wallis test followed by two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Figure A3. T cell populations in tumor samples classified according to the Lund Taxonomy molecular subtypes. Percentage of T cells (A), CD8+ T cells (B), CD4+ T cells (C), Tregs (D), and Th (E) out of all viable cells. Statistical significance was determined by Kruskal–Wallis test followed by two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Cells 13 00926 g0a3
Figure A4. Pathway analysis of differentially expressed genes between CD8+ T cells from GU subtype and the control (A) and Uro subtype and the controls (B). padj = 0.05 marked in red.
Figure A4. Pathway analysis of differentially expressed genes between CD8+ T cells from GU subtype and the control (A) and Uro subtype and the controls (B). padj = 0.05 marked in red.
Cells 13 00926 g0a4

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Patel, V.G.; Oh, W.K.; Galsky, M.D. Treatment of muscle-invasive and advanced bladder cancer in 2020. CA A Cancer J. Clin. 2020, 70, 404–423. [Google Scholar] [CrossRef] [PubMed]
  3. Witjes, J.A.; Bruins, H.M.; Cathomas, R.; Compérat, E.M.; Cowan, N.C.; Gakis, G.; Hernández, V.; Linares Espinós, E.; Lorch, A.; Neuzillet, Y.; et al. European Association of Urology Guidelines on Muscle-invasive and Metastatic Bladder Cancer: Summary of the 2020 Guidelines. Eur. Urol. 2021, 79, 82–104. [Google Scholar] [CrossRef] [PubMed]
  4. van der Heijden, A.G.; Witjes, J.A. Recurrence, Progression, and Follow-Up in Non–Muscle-Invasive Bladder Cancer. Eur. Urol. Suppl. 2009, 8, 556–562. [Google Scholar] [CrossRef]
  5. Yeung, C.; Dinh, T.; Lee, J. The health economics of bladder cancer: An updated review of the published literature. Pharmacoeconomics 2014, 32, 1093–1104. [Google Scholar] [CrossRef] [PubMed]
  6. Höglund, M.; Bernardo, C.; Sjödahl, G.; Eriksson, P.; Axelson, H.; Liedberg, F. The Lund taxonomy for bladder cancer classification—From gene expression clustering to cancer cell molecular phenotypes, and back again. J. Pathol. 2023, 259, 369–375. [Google Scholar] [CrossRef] [PubMed]
  7. Sjödahl, G.; Lauss, M.; Lövgren, K.; Chebil, G.; Gudjonsson, S.; Veerla, S.; Patschan, O.; Aine, M.; Fernö, M.; Ringnér, M.; et al. A Molecular Taxonomy for Urothelial Carcinoma. Clin. Cancer Res. 2012, 18, 3377–3386. [Google Scholar] [CrossRef] [PubMed]
  8. Mariathasan, S.; Turley, S.J.; Nickles, D.; Castiglioni, A.; Yuen, K.; Wang, Y.; Kadel, E.E., III; Koeppen, H.; Astarita, J.L.; Cubas, R.; et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature 2018, 554, 544–548. [Google Scholar] [CrossRef] [PubMed]
  9. Sharma, P.; Shen, Y.; Wen, S.; Yamada, S.; Jungbluth, A.A.; Gnjatic, S.; Bajorin, D.F.; Reuter, V.E.; Herr, H.; Old, L.J. CD8 tumor-infiltrating lymphocytes are predictive of survival in muscle-invasive urothelial carcinoma. Proc. Natl. Acad. Sci. USA 2007, 104, 3967–3972. [Google Scholar] [CrossRef]
  10. Joseph, M.; Enting, D. Immune responses in bladder cancer-role of immune cell populations, prognostic factors and therapeutic implications. Front. Oncol. 2019, 9, 1270. [Google Scholar] [CrossRef]
  11. Winerdal, M.E.; Marits, P.; Winerdal, M.; Hasan, M.; Rosenblatt, R.; Tolf, A.; Selling, K.; Sherif, A.; Winqvist, O. FOXP3 and survival in urinary bladder cancer. BJU Int. 2011, 108, 1672–1678. [Google Scholar] [CrossRef] [PubMed]
  12. Winerdal, M.E.; Krantz, D.; Hartana, C.A.; Zirakzadeh, A.A.; Linton, L.; Bergman, E.A.; Rosenblatt, R.; Vasko, J.; Alamdari, F.; Hansson, J.; et al. Urinary Bladder Cancer Tregs Suppress MMP2 and Potentially Regulate Invasiveness. Cancer Immunol. Res. 2018, 6, 528–538. [Google Scholar] [CrossRef]
  13. Powles, T.; Kockx, M.; Rodriguez-Vida, A.; Duran, I.; Crabb, S.J.; Van Der Heijden, M.S.; Szabados, B.; Pous, A.F.; Gravis, G.; Herranz, U.A.; et al. Clinical efficacy and biomarker analysis of neoadjuvant atezolizumab in operable urothelial carcinoma in the ABACUS trial. Nat. Med. 2019, 25, 1706–1714. [Google Scholar] [CrossRef] [PubMed]
  14. Li, H.; Zhang, Q.; Shuman, L.; Kaag, M.; Raman, J.D.; Merrill, S.; DeGraff, D.J.; Warrick, J.I.; Chen, G. Evaluation of PD-L1 and other immune markers in bladder urothelial carcinoma stratified by histologic variants and molecular subtypes. Sci. Rep. 2020, 10, 1439. [Google Scholar] [CrossRef] [PubMed]
  15. Sjödahl, G.; Lövgren, K.; Lauss, M.; Chebil, G.; Patschan, O.; Gudjonsson, S.; Månsson, W.; Fernö, M.; Leandersson, K.; Lindgren, D.; et al. Infiltration of CD3⁺ and CD68⁺ cells in bladder cancer is subtype specific and affects the outcome of patients with muscle-invasive tumors. Urol. Oncol. 2014, 32, 791–797. [Google Scholar] [CrossRef] [PubMed]
  16. Abolhalaj, M.; Sincic, V.; Lilljebjörn, H.; Sandén, C.; Aab, A.; Hägerbrand, K.; Ellmark, P.; Borrebaeck, C.A.K.; Fioretos, T.; Lundberg, K. Transcriptional profiling demonstrates altered characteristics of CD8(+) cytotoxic T-cells and regulatory T-cells in TP53-mutated acute myeloid leukemia. Cancer Med. 2022, 11, 3023–3032. [Google Scholar] [CrossRef] [PubMed]
  17. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  18. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [PubMed]
  19. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
  20. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  21. Zhu, A.; Ibrahim, J.G.; Love, M.I. Heavy-tailed prior distributions for sequence count data: Removing the noise and preserving large differences. Bioinformatics 2019, 35, 2084–2092. [Google Scholar] [CrossRef] [PubMed]
  22. Sergushichev, A.A. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation. bioRxiv 2016, 060012. [Google Scholar]
  23. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [PubMed]
  24. Lipponen, P.K.; Eskelinen, M.J.; Jauhiainen, K.; Harju, E.; Terho, R. Tumour infiltrating lymphocytes as an independent prognostic factor in transitional cell bladder cancer. Eur. J. Cancer 1993, 29, 69–75. [Google Scholar] [CrossRef] [PubMed]
  25. Tirosh, I.; Izar, B.; Prakadan, S.M.; Wadsworth, M.H.; Treacy, D.; Trombetta, J.J.; Rotem, A.; Rodman, C.; Lian, C.; Murphy, G.; et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 2016, 352, 189–196. [Google Scholar] [CrossRef] [PubMed]
  26. Li, P.; Zhu, X.; Cao, G.; Wu, R.; Li, K.; Yuan, W.; Chen, B.; Sun, G.; Xia, X.; Zhang, H.; et al. 1α,25(OH)(2)D(3) reverses exhaustion and enhances antitumor immunity of human cytotoxic T cells. J. Immunother. Cancer 2022, 10, e003477. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, J.-Q.; Salas, L.A.; Wiencke, J.K.; Koestler, D.C.; Molinaro, A.M.; Andrew, A.S.; Seigne, J.D.; Karagas, M.R.; Kelsey, K.T.; Christensen, B.C. Matched analysis of detailed peripheral blood and tumor immune microenvironment profiles in bladder cancer. Epigenomics 2024, 16, 41–56. [Google Scholar] [CrossRef]
  28. Toren, P.; Wilkins, A.; Patel, K.; Burley, A.; Gris, T.; Kockelbergh, R.; Lodhi, T.; Choudhury, A.; Bryan, R.T. The sex gap in bladder cancer survival—A missing link in bladder cancer care? Nat. Rev. Urol. 2024, 21, 181–192. [Google Scholar] [CrossRef]
  29. Andreassen, B.K.; Grimsrud, T.K.; Haug, E.S. Bladder cancer survival: Women better off in the long run. Eur. J. Cancer 2018, 95, 52–58. [Google Scholar] [CrossRef]
  30. Richters, A.; Dickman, P.W.; Witjes, J.A.; Boormans, J.L.; Kiemeney, L.A.L.M.; Aben, K.K.H. Bladder cancer survival: Women only fare worse in the first two years after diagnosis. Urol. Oncol. Semin. Orig. Investig. 2019, 37, 853–861. [Google Scholar] [CrossRef]
  31. Conforti, F.; Pala, L.; Pagan, E.; Bagnardi, V.; De Pas, T.; Queirolo, P.; Pennacchioli, E.; Catania, C.; Cocorocchio, E.; Ferrucci, P.F.; et al. Sex-Based Dimorphism of Anticancer Immune Response and Molecular Mechanisms of Immune Evasion. Clin. Cancer Res. 2021, 27, 4311–4324. [Google Scholar] [CrossRef] [PubMed]
  32. Chenard, S.; Jackson, C.; Vidotto, T.; Chen, L.; Hardy, C.; Jamaspishvilli, T.; Berman, D.; Siemens, D.R.; Koti, M. Sexual Dimorphism in Outcomes of Non-muscle-invasive Bladder Cancer: A Role of CD163+ Macrophages, B cells, and PD-L1 Immune Checkpoint. Eur. Urol. Open Sci. 2021, 29, 50–58. [Google Scholar] [CrossRef] [PubMed]
  33. Robertson, A.G.; Kim, J.; Al-Ahmadie, H.; Bellmunt, J.; Guo, G.; Cherniack, A.D.; Hinoue, T.; Laird, P.W.; Hoadley, K.A.; Akbani, R.; et al. Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer. Cell 2017, 171, 540–556.e25. [Google Scholar] [CrossRef] [PubMed]
  34. Benítez, R.; Yu, K.; Sirota, M.; Malats, N.; Pineda, S. Characterization of the tumor-infiltrating immune repertoire in muscle invasive bladder cancer. Front. Immunol. 2023, 14, 986598. [Google Scholar] [CrossRef] [PubMed]
  35. Saleh, R.; Elkord, E. FoxP3+ T regulatory cells in cancer: Prognostic biomarkers and therapeutic targets. Cancer Lett. 2020, 490, 174–185. [Google Scholar] [CrossRef] [PubMed]
  36. Li, Q.; Lu, J.; Li, J.; Zhang, B.; Wu, Y.; Ying, T. Antibody-based cancer immunotherapy by targeting regulatory T cells. Front. Oncol. 2023, 13, 1157345. [Google Scholar] [CrossRef]
  37. Oh, D.Y.; Kwek, S.S.; Raju, S.S.; Li, T.; McCarthy, E.; Chow, E.; Aran, D.; Ilano, A.; Pai, C.-C.S.; Rancan, C.; et al. Intratumoral CD4+ T Cells Mediate Anti-tumor Cytotoxicity in Human Bladder Cancer. Cell 2020, 181, 1612–1625.e13. [Google Scholar] [CrossRef]
Figure 2. Frequency of immune cells and T cell populations in tumor samples classified according to the Lund Taxonomy molecular subtypes. Infiltration of CD45+ immune cells among all viable cells (A). Immune infiltration in MIBC samples only (B). Percentages of T cells, CD4+ T cells, and CD8+ T cells out of all CD45+ cells (CE). Ratio of CD4+/CD8+ T cells (F). Percentages of Tregs and Th cells out of CD45+ cells (G,H). Statistical significance was determined through the Kruskal–Wallis test followed by the two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Figure 2. Frequency of immune cells and T cell populations in tumor samples classified according to the Lund Taxonomy molecular subtypes. Infiltration of CD45+ immune cells among all viable cells (A). Immune infiltration in MIBC samples only (B). Percentages of T cells, CD4+ T cells, and CD8+ T cells out of all CD45+ cells (CE). Ratio of CD4+/CD8+ T cells (F). Percentages of Tregs and Th cells out of CD45+ cells (G,H). Statistical significance was determined through the Kruskal–Wallis test followed by the two-stage linear step-up procedure of Benjamini, Krieger, and Yekutieli post-test.
Cells 13 00926 g002
Figure 3. Transcriptional analysis of CD8+ T cells from non-malignant bladder tissue and from GU and Uro tumor biopsies. Principal component analysis (A). Heatmap of differentially expressed genes between CD8+ T cells from GU and Uro subtypes (B). GSEA using an exhaustion profile (Tirosh et al.’s 28-gene exhaustion profile) [25] (C). Volcano plot with the differentially expressed genes (GU vs. control) marked in red, with the top 40 annotated (D). Volcano plot with the differentially expressed genes (Uro vs. control) marked in red, with the top 40 annotated (E). Percentage of CD8+ cells expressing the exhaustion markers CTLA-4, LAG3, PD-1, TIGIT, or TIM-3 measured by flow cytometry (F). Statistical significance assessed by Mann–Whitney U test.
Figure 3. Transcriptional analysis of CD8+ T cells from non-malignant bladder tissue and from GU and Uro tumor biopsies. Principal component analysis (A). Heatmap of differentially expressed genes between CD8+ T cells from GU and Uro subtypes (B). GSEA using an exhaustion profile (Tirosh et al.’s 28-gene exhaustion profile) [25] (C). Volcano plot with the differentially expressed genes (GU vs. control) marked in red, with the top 40 annotated (D). Volcano plot with the differentially expressed genes (Uro vs. control) marked in red, with the top 40 annotated (E). Percentage of CD8+ cells expressing the exhaustion markers CTLA-4, LAG3, PD-1, TIGIT, or TIM-3 measured by flow cytometry (F). Statistical significance assessed by Mann–Whitney U test.
Cells 13 00926 g003
Table 1. Cohort characteristics. Variables depicted as absolute numbers and percentage of all tumor samples. Spread described as interquartile range (IQR). Samples were collected between the years 2018 and 2022.
Table 1. Cohort characteristics. Variables depicted as absolute numbers and percentage of all tumor samples. Spread described as interquartile range (IQR). Samples were collected between the years 2018 and 2022.
Patient Samplesn%
Tumor4886
Control814
Healthy controls
Age, median (IQR)55 (46–59.5)
Sexn%
Male450
Female450
Tumor samples
Age, median (IQR)76 (68.5–79.5)
Sexn%
Male3981
Female919
Stage (WHO 1999)n%
Ta1735
T11123
T2817
T3919
T436
Graden%
G112
G21531
G33267
Clinical nodal stagen%
N04185
N136
N224
N324
Metastasisn%
MX1735
M02960
M124
Molecular subtypen%
Basal/Squamous36
GU817
Uro3777
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

Sincic, V.; Arlenhold, K.F.; Richtmann, S.; Lilljebjörn, H.; Eriksson, P.; Sjödahl, G.; Wokander, M.; Hägerbrand, K.; Ellmark, P.; Fioretos, T.; et al. Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes. Cells 2024, 13, 926. https://doi.org/10.3390/cells13110926

AMA Style

Sincic V, Arlenhold KF, Richtmann S, Lilljebjörn H, Eriksson P, Sjödahl G, Wokander M, Hägerbrand K, Ellmark P, Fioretos T, et al. Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes. Cells. 2024; 13(11):926. https://doi.org/10.3390/cells13110926

Chicago/Turabian Style

Sincic, Viktor, Ken F. Arlenhold, Sarah Richtmann, Henrik Lilljebjörn, Pontus Eriksson, Gottfrid Sjödahl, Mats Wokander, Karin Hägerbrand, Peter Ellmark, Thoas Fioretos, and et al. 2024. "Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes" Cells 13, no. 11: 926. https://doi.org/10.3390/cells13110926

APA Style

Sincic, V., Arlenhold, K. F., Richtmann, S., Lilljebjörn, H., Eriksson, P., Sjödahl, G., Wokander, M., Hägerbrand, K., Ellmark, P., Fioretos, T., Borrebaeck, C. A. K., Liedberg, F., & Lundberg, K. (2024). Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes. Cells, 13(11), 926. https://doi.org/10.3390/cells13110926

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