Next Article in Journal
Review of Treatment Options for the Management of Advanced Stage Hodgkin Lymphoma
Next Article in Special Issue
Clearing up Clear Cell: Clarifying the Immuno-Oncology Treatment Landscape for Metastatic Clear Cell RCC
Previous Article in Journal
Development and Implementation of Proton Therapy for Hodgkin Lymphoma: Challenges and Perspectives
Previous Article in Special Issue
Immune Checkpoints Inhibitors and Chemotherapy as First-Line Treatment for Metastatic Urothelial Carcinoma: A Network Meta-Analysis of Randomized Phase III Clinical Trials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma

by
Nicholas H. Chakiryan
1,*,
Gregory J. Kimmel
2,
Youngchul Kim
3,
Joseph O. Johnson
4,
Noel Clark
5,
Ali Hajiran
1,
Andrew Chang
1,
Ahmet M. Aydin
1,
Logan Zemp
1,
Esther Katende
1,
Jad Chahoud
1,
Meghan C. Ferrall-Fairbanks
2,
Philippe E. Spiess
1,
Natasha Francis
1,
Michelle Fournier
1,
Jasreman Dhillon
6,
Jong Y. Park
7,
Liang Wang
8,
James J. Mulé
9,
Philipp M. Altrock
2 and
Brandon J. Manley
1
add Show full author list remove Hide full author list
1
Department of Genitourinary Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
2
Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
3
Department of Biostatistics and Bioinformatics, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
4
Analytic Microcopy Shared Resource, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
5
Tissue Core Shared Resource, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
6
Department of Pathology, H. Lee Moffitt Cancer Center, Tampa, FL 33612, USA
7
Department of Cancer Epidemiology, H. Lee Moffitt Cancer Center, Tampa, FL 33612, USA
8
Department of Tumor Biology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
9
Immunology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(15), 3743; https://doi.org/10.3390/cancers13153743
Submission received: 24 May 2021 / Revised: 12 July 2021 / Accepted: 22 July 2021 / Published: 26 July 2021
(This article belongs to the Special Issue Advances in Immunotherapy for Genitourinary Malignancies)

Abstract

:

Simple Summary

Cancer-associated fibroblasts (CAFs) are highly prevalent cells in the clear cell renal cell carcinoma (ccRCC) tumor immune microenvironment. CAFs are thought to potentiate tumor proliferation primarily through paracrine interactions, as evidenced by laboratory-based studies. We sought to corroborate these findings using surgically removed tissue samples from 96 patients with metastatic ccRCC and associate geospatial relationships between CAFs and rapidly proliferating tumor cells with survival outcomes. We found that CAFs exhibited more geospatial clustering with proliferating tumor cells than with dying tumor cells, and patients whose samples exhibited higher tumor cell proliferation had worse overall survival and were more likely to be resistant to systemic tyrosine-kinase-inhibiting targeted therapies. Immunotherapy resistance was not associated with the geospatial metrics measured in this analysis. Overall, these findings suggest that close proximity to CAFs potentiates tumor cell proliferation, worsening survival and conferring resistance to targeted therapies.

Abstract

Cancer-associated fibroblasts (CAF) are highly prevalent cells in the tumor microenvironment in clear cell renal cell carcinoma (ccRCC). CAFs exhibit a pro-tumor effect in vitro and have been implicated in tumor cell proliferation, metastasis, and treatment resistance. Our objective is to analyze the geospatial distribution of CAFs with proliferating and apoptotic tumor cells in the ccRCC tumor microenvironment and determine associations with survival and systemic treatment. Pre-treatment primary tumor samples were collected from 96 patients with metastatic ccRCC. Three adjacent slices were obtained from 2 tumor-core regions of interest (ROI) per patient, and immunohistochemistry (IHC) staining was performed for αSMA, Ki-67, and caspase-3 to detect CAFs, proliferating cells, and apoptotic cells, respectively. H-scores and cellular density were generated for each marker. ROIs were aligned, and spatial point patterns were generated, which were then used to perform spatial analyses using a normalized Ripley’s K function at a radius of 25 μm (nK(25)). The survival analyses used an optimal cut-point method, maximizing the log-rank statistic, to stratify the IHC-derived metrics into high and low groups. Multivariable Cox regression analyses were performed accounting for age and International Metastatic RCC Database Consortium (IMDC) risk category. Survival outcomes included overall survival (OS) from the date of diagnosis, OS from the date of immunotherapy initiation (OS-IT), and OS from the date of targeted therapy initiation (OS-TT). Therapy resistance was defined as progression-free survival (PFS) <6 months, and therapy response was defined as PFS >9 months. CAFs exhibited higher cellular clustering with Ki-67+ cells than with caspase-3+ cells (nK(25): Ki-67 1.19; caspase-3 1.05; p = 0.04). The median nearest neighbor (NN) distance from CAFs to Ki-67+ cells was shorter compared to caspase-3+ cells (15 μm vs. 37 μm, respectively; p < 0.001). Multivariable Cox regression analyses demonstrated that both high Ki-67+ density and H-score were associated with worse OS, OS-IT, and OS-TT. Regarding αSMA+CAFs, only a high H-score was associated with worse OS, OS-IT, and OS-TT. For caspase-3+, high H-score and density were associated with worse OS and OS-TT. Patients whose tumors were resistant to targeted therapy (TT) had higher Ki-67 density and H-scores than those who had TT responses. Overall, this ex vivo geospatial analysis of CAF distribution suggests that close proximity clustering of tumor cells and CAFs potentiates tumor cell proliferation, resulting in worse OS and resistance to TT in metastatic ccRCC.

1. Introduction

In normal tissue, fibroblast activation and the subsequent release of cytokines, angiogenic mediators, and growth factors are physiologic responses to tissue injury or stress [1,2]. Cancer-associated fibroblasts (CAF) are fibroblasts that have been permanently activated by adjacent tumor cells, which repurpose physiologic fibroblast activity into a pro-tumor survival advantage [1,2,3]. CAFs are highly prevalent cells in the tumor microenvironment in clear cell renal cell carcinoma (ccRCC) and have been implicated in facilitating tumor cell proliferation, angiogenesis, metastasis, and therapy resistance [4,5,6]. Several well-described pro-tumor properties of CAFs are mediated through hypoxia-inducible-factor-1 (HIF1), a pathway that drives oncogenesis in kidney cancer and is upregulated in the majority of ccRCC tumors via alteration of the von-Hippel Lindau (VHL) gene [7]. Additionally, molecular receptors mediating the HIF1 pathway are the primary targets for tyrosine-kinase–inhibiting targeted therapies (TT), which are frequently used in metastatic ccRCC [8]. Previous work has demonstrated that increased CAF density in ccRCC tumors is associated with worse overall survival (OS) [4,6].
The induction of tumor cell proliferation has been postulated as a major mechanism of CAF-mediated pro-tumor activity [1]. Immunohistochemical (IHC) staining with Ki-67, a nuclear protein that is present during active phases of the cell cycle and absent from resting cells, has been shown to be an excellent marker for identifying rapidly proliferating tumor cells [9,10]. This staining is clinically relevant for grading and prognosis in several primary cancer sites, including breast and colorectal cancers [11,12,13]. Currently, there is no role for Ki-67 staining in ccRCC guideline-based clinical practice, though a handful of studies have associated high Ki-67 staining with more advanced disease stage, worse OS, and worse cancer-specific survival (CSS) [14,15,16,17].
Despite the well-described in vitro relationship between CAFs and increased tumor cell proliferation, no prior study has investigated the interplay between CAFs and rapidly proliferating tumor cells in ccRCC tumor samples. The discovery of distinct infiltration patterns and spatial relationships between CAFs and rapidly proliferating tumor cells could support previous in vitro study findings with ex vivo evidence and inform ongoing research into microenvironment-modulating antineoplastic therapies targeting fibroblasts [18]. Our primary objective was to investigate cellular distribution patterns and spatial relationships between CAFs and proliferating and apoptotic tumor cells in primary tumor samples from patients with metastatic ccRCC. Additionally, to evaluate the associations between these measures and OS. Our secondary objective was to determine associations between treatment outcomes and targeted therapy and immunotherapy, as well as define spatial relationships between CAFs, proliferating tumors cells, and apoptotic tumor cells.

2. Materials and Methods

2.1. Patient and Sample Selection

Samples were included from patients with primary ccRCC who had metastatic ccRCC at the time of sample collection and whose tumor specimens were available in formalin-fixed paraffin-embedded blocks. Included patients had received either TT, immunotherapy (IT), or combination TT/IT as systemic treatment from October 2004 to September 2020. Written informed consent was obtained from all tissue donors. All tumor and normal tissue samples were obtained through protocols approved by the institutional review board (H. Lee Moffitt Cancer Center and Research Institute’s Total Cancer Care protocol MCC# 14690; Advarra IRB Pro00014441). The general workflow of our methods is shown in Figure 1.

2.2. Immunohistochemical Specimen Preparation

Three immediately adjacent slides (3 µm thickness) were prepared from each tissue block. IHC staining for alpha-smooth muscle actin ((αSMA a marker for activated fibroblasts), Ki-67 (a marker for proliferating cells), and caspase-3 (a marker for cells undergoing apoptosis) were performed; 1 stain was used on each slide. Slides were stained using a Ventana Discovery XT automated system (Ventana Medical Systems, Tucson, AZ, USA) as per manufacturer’s protocol, with proprietary reagents. Briefly, slides were deparaffinized on the automated system with the EZ Prep solution (Ventana). A heat-induced antigen retrieval method was used in Cell Conditioning (Ventana). The rabbit primary antibody that reacts to Ki-67 (#790-4286 (Ventana)) was used at a prediluted strength and incubated for 16 min. The rabbit primary antibody that reacts to Cleaved Caspase 3 (#9661 (Cell Signaling, Danvers, MA, USA)) was used at a 1:4000 concentration in Dako antibody diluent (Carpenteria, CA, USA) and incubated for 60 min. The rabbit primary antibody that reacts to αSMA (#ab32575 (Abcam, Cambridge, MA, USA)) was used at a 1:250 concentration in Dako antibody diluent (Carpenteria, CA, USA) and incubated for 32 min. For all stains, the Ventana OmniMap Anti-Rabbit Secondary Antibody was used for 16 min. The detection system used was the Ventana ChromoMap kit, and slides were then counterstained with hematoxylin. Slides were then dehydrated and cover-slipped as per normal laboratory protocol.

2.3. Quantitative Digital Image Analysis

Slides were digitized with a Leica Aperio AT2 slide scanner (Vista, CA, USA) using a 20 X/0.75 NA objective lens. An experienced genitourinary pathologist (JD) used the annotation pen tool in the Aperio Imagescope software to define the tumor-core zone in each hematoxylin and eosin image. High-resolution Aperio SVS images representing the 3 IHC stains were imported into Visiopharm (Hoersholm, Denmark) for quantitative digital image analysis. We leveraged a tissue alignment algorithm available in Visiopharm to optimally align 3 adjacent slides of each patient sample set. These slides were visually inspected to ensure appropriate alignment. Two equally sized regions of interest (ROI) were selected from the tumor-core zone from each image; tumor cells were evenly distributed throughout each ROI, and cellular appearance was consistent with that seen in the remainder of the slide. The ROI sizes were standardized at 3426 × 1379 pixels, at a pixel resolution of 0.502 µm/pixel. Thresholds for staining positivity were set by an experienced digital pathology image analyst (JJ) and confirmed by a study pathologist (JD). These thresholds were used with Visiopharm’s cell detection algorithms to identify and categorize cells into negative, weak, moderate, and strong bins on the basis of staining intensity. This intensity-based distribution of cells is comparable to the qualitative method used by pathologists, in which staining intensity is categorized as 0, 1+, 2+, and 3+ [19]. Percent positivity and H-Scores were calculated for each ROI using this intensity data. H-score is an approach that globally quantifies intensity and percent positivity throughout the entire ROI into 1 score, according to the following formula: H-score = (1 × (%cells weak)) + (2 × (%cells moderate)) + (3 × (%cells strong)) [20]. H-scores range from 0 to 300, with 0 representing no cell staining for the marker of interest and 300 representing every cell staining with the highest intensity [20]. Additionally, Cartesian coordinates for the (x,y) location of each cell’s central mass, with each cell’s associated marker status, were abstracted from the digital image using Visiopharm.

2.4. Cellular Distribution and Spatial Analysis

For each ROI, the per-cell Cartesian coordinate and marker positivity data were converted into spatial point patterns. Cell density was calculated as the number of positive cells per mm2. ROIs containing ≥10 cells positive for a relevant marker were considered eligible for spatial analysis. As there is no previously validated standard for this cutoff, the ≥10 cell cutoff was agreed upon through the consensus of the authors. Cells classified as exhibiting strong staining intensity, as defined above, were considered positive for the purpose of spatial analysis. Cellular clustering was quantified using Ripley’s K function, a methodology for quantifying spatial heterogeneity that is most commonly used in ecology and economics, with isotropic edge correction, and the following normalization was applied: nK(r) = K(r)/πr2, as described previously [21,22]. As such, the expected value of nK(r) for complete spatial randomness is 1.0, assuming a homogenous Poisson process [21]. Values of nK(r) > 1.0 represent cellular clustering and values < 1.0 represent cellular dispersion. The range of possible values for nK(r) is 0 to infinity. The nK(r) value is an observed over expected ratio (i.e., αSMA/Ki-67 nK(25 µm) = 1.30 can be interpreted as: “There were 30% more Ki-67+ cells within a 25 µm radius of each αSMA cell than would be expected if the cells were randomly distributed.”).
Initially, two search-circle radii were utilized for the spatial analysis. To reflect cellular clustering at a localized distance, nK(r) at a radius of 25 µm was used in this analysis and will henceforth be referred to as nK(25). The search-circle radius value of 25 µm was selected, as it represents approximately double that of a typical ccRCC tumor cell radius and, as such, should represent the area in the immediate vicinity of the cell. To reflect a more global view of cellular clustering, a radius of 125 µm was utilized, which will be referenced as nK(125). Subsequently, we determined that Spearman’s correlations demonstrated a very strong correlation between nK(25) and nK(125) (r > 0.90). Given this redundancy, we utilized nK(25), not nK(125), for the remainder of the analysis.
As a second independent measure of spatial analysis, linear nearest neighbor distances were determined from each cell to its nearest neighbor cell among each of the 3 cell types. Scaled two-dimensional kernel density plots were generated to analyze the distribution of nearest neighbor distances from each αSMA+ cell to its nearest neighbor caspase-3+ and Ki-67+ cell. These nearest neighbor density plots were generated with the cohort stratified by ccRCC versus normal kidney, treatment response versus resistance to TT, and treatment response versus resistance to IT.

2.5. Variable Definitions

Age was defined as life years at the time of ccRCC diagnosis. Tumor grade was defined per the histologic classification criteria proposed by the International Society of Urologic Pathologists and implemented by the World Health Organization (ISUP/WHO) [23]. International Metastatic RCC Database (IMDC) scores were determined for each patient and categorized into good, intermediate, and poor risk groups, as previously described [24]. OS was determined from the date of metastatic ccRCC diagnosis to the date of death or censoring at the last follow-up. Treatments classified as IT included both immune checkpoint inhibitors and high-dose interleukin-2 (IL-2). Treatments classified as TT included small-molecule tyrosine-kinase-inhibiting therapies approved for first-line treatment in ccRCC [17]. OS from the date of the first receipt of IT (OS-IT) and OS from the date of first receipt of TT (OS-TT) were also determined. Response to therapy was defined as clinical progression-free survival greater than 9 months from treatment initiation, and resistance was defined as progression <6 months from treatment initiation. Patients who progressed between 6 and 9 months after initiating therapy were not classified as responsive or resistant. Statistical significance was defined as a two-tailed alpha-risk of 0.05 or less.

2.6. Statistical Analysis

In addition to demographic, clinical, and pathologic characteristics, patient-level data included marker densities and intensities, H-scores, and uni- and bi-variate spatial distribution metrics for each of the 3 included IHC markers. As each patient had 2 ROIs analyzed, the IHC-derived metrics were averaged such that each patient had 1 value for each metric. Spearman’s correlation coefficients were determined between each pairwise combination to assess the interactions between the IHC-derived metrics.
As standardized cutoffs do not exist for the IHC-derived metrics utilized in this analysis, optimal cut-points were determined for each metric, maximizing the log-rank test statistics for OS, as previously described [25]. The survival analysis used a multivariable Cox proportional hazards regression for each metric, utilizing age and IMDC risk category as covariates. False discovery rate- (FDR-)adjusted p values were calculated for multiple comparison correction. This analysis was repeated for the OS-IT and OS-TT endpoints.
Patients were then stratified by response or non-response to IT and TT, and values of each IHC-derived metric were compared between responders and non-responders using a two-sample t-test.
The distribution of nearest neighbor distances from αSMA+ cells to their nearest caspase-3+ cell was compared to that of the distance to the nearest Ki-67- cell using a Wilcoxon test. This process was repeated after stratifying the cohort by IT response and non-response and TT response and non-response.

2.7. Statistical Software

All statistical and spatial analyses were performed using R version 4.0.2 (The R Foundation for Statistical Computing (Vienna, Austria)). The maxstat package was used to determine optimal cut-points, survival and survminer packages used for the survival analysis, and the spatstat, MASS, and FNN packages for spatial analysis.

3. Results

3.1. Study Population

The study population included 96 patients (median age, 60 years (IQR 33–87); male, 68 (71%); white race, 90 (94%)) (Table 1). Fifty-one patients (53%) underwent first-line IT (18 responders, 25 resistant, and 8 indeterminate); 42 patients (44%) had first-line TT (27 responders, 9 resistant, and 6 indeterminate); and 3 patients (3%) had first-line combination IT/TT. Median follow-up was 54 months (IQR, 41–97) for the 37 patients (38%) who were alive at last contact.

3.2. CAFs Are Highly Prevalent in the ccRCC Tumor Microenvironment

αSMA+ CAFs had the highest median cell density and H-score (density, 1904 cells/mm2 (IQR 1266–3382); H-score 35.5 (IQR 24.0–51.2)), as compared with respective values for Ki-67+ and caspase-3+ tumor cells (Ki-67+ density, 1087 cells/mm2 (IQR, 433–1909); Ki-67 H-score 7.0 (IQR 3.1–15.5); caspase-3+ density, 203 cells/mm2 (IQR 64–547); caspase-3+ H-score, 0.8 (IQR, 0.3–2.2)) (Wilcoxon p  <  0.001 for all comparisons within density and H-score groups) (Figure 2A).

3.3. CAFs Are Spatially Clustered with Proliferating Tumor Cells

Using the normalized Ripley’s K function, we found that αSMA+ CAFs were more spatially clustered with Ki-67+ tumor cells than they were with caspase-3+ tumor cells (nK(25), 1.15 vs. 1.09, respectively; p = 0.045) (Figure 2B). Caspase-3+ tumor cells were more spatially clustered with Ki-67+ tumor cells than they were with αSMA+ CAFs (nK(25), 1.16 vs. 1.09, respectively; p = 0.035). The amount of cellular clustering demonstrated between αSMA+ CAFs and Ki-67+ cells was similar to that between Ki-67+ tumor cells and caspase-3+ tumor cells (p = 0.434).
The median nearest-neighbor (NN) distance from αSMA+ CAFs to the nearest Ki-67+ cell was shorter than the distance to the nearest caspase-3+ cell (14.4 vs. 34.9 µm; p < 0.001) (Figure 2C). The median NN distance from Ki-67+ cells to the nearest caspase-3+ cell was shorter than the distance from αSMA+ CAFs to the nearest caspase-3+ cell (31.3 vs. 34.9 µm; p < 0.001).
Overall, the density and spatial analyses suggest a tumor architecture in which αSMA+ CAFs and Ki-67+ tumor cells demonstrate the highest clustering and shortest NN distance, αSMA+ CAFs and caspase-3+ tumor cells have the least amount of clustering and longest NN distance, and cell densities are the highest for αSMA+ CAFs and lowest for caspase-3+ tumor cells. Figure 3 shows an example of the hypothetical median tumor architecture of these 3 markers, as suggested by the distribution and spatial analysis metrics for the entire cohort.

3.4. CAF Density and αSMA H-Score Are Positively Correlated with Proliferating and Apoptotic Tumor Cells

Spearman’s correlations demonstrated that the αSMA H-score had a moderate positive correlation (r = 0.4–0.6) with αSMA density, Ki-67 H-score, and Caspase H-score and a weak positive correlation (r = 0.2–0.4) with caspase-3 density and Ki-67 density (Figure 2D). αSMA density had a moderate positive correlation with Ki-67 density, caspase-3 density, and αSMA H-sore and a weak positive correlation with Ki-67 H-score and caspase-3 H-score. Clustering metrics and density/H-score metrics were consistently found to have either no correlation or a weak negative correlation (r = −0.4–0). Clustering metrics at a local or global scale (nK(25) or nK(125), respectively) were found to have a strong positive correlation with each other within marker-types (r = 0.6–1.0). Complete information regarding pairwise correlations between the IHC-derived metrics can be found in Figure 2D.
None of the IHC-derived metrics were significantly correlated with primary tumor size or grade, confirming that these measures are not simply surrogates for conventionally measures (Supplemental Figure S1).

3.5. High αSMA, Ki-67, and Caspase-3 H-Scores Are Associated with Significantly Worse OS

Multivariable Cox proportional hazards regression analysis showed that patients with a high αSMA H-core had significantly worse OS (HR 1.45 (1.05–2.01); p = 0.02), OS-IT (HR 1.51 (1.05–2.17); p = 0.03), and OS-TT (HR 1.46 (1.06–2.02); p = 0.012 than patients with a low αSMA H-score. Interestingly, stratification by αSMA CAF cell density alone had no statistically significant association with OS, OS-IT, or OS-TT (p = 0.32, 0.17, and 0.24, respectively) compared with patients with low αSMA+ CAF cell density (Figure 4A).
Patients with a high Ki-67 H-score had significantly worse OS (HR 1.79 (1.40–2.29); p < 0.001), OS-IT (HR 1.47 (1.13–1.91); p < 0.001), and OS-TT (HR 1.71 (1.35–2.17); p < 0.001) than patients with a low Ki-67 H-score. Patients with high Ki-67+ cell density had significantly worse OS (HR, 1.30 (1.06–1.61); p = 0.01), and OS-TT (HR, 1.27 (1.04–1.57); p = 0.02), with no significant difference in OS-IT (HR, 1.21 (0.98–1.50); p = 0.08), than patients with low Ki-67+ cell density (Figure 4B).
Patients with a high caspase-3 H-score had significantly worse OS (HR 1.33 (1.01–1.75); p = 0.04) and OS-TT (HR 1.46 (1.14–1.88); p < 0.001), with no significant difference in OS-IT (HR 1.19 (0.89–1.58); p = 0.24), than patients with low caspase-3 H-scores. Patients with high caspase-3+ cell density had significantly worse OS (HR 1.26 (1.02–1.56); p = 0.03), and OS-TT (HR 1.25 (1.02–1.53); p = 0.03), with no significant difference in OS-IT (HR 1.18 (0.95–1.46); p = 0.13) (Figure 4C).
Using Ripley’s K function, we did not find any spatial clustering metrics predictive of OS, OS-IT, or OS-TT. The results of the multivariable Cox regression analysis assessing all IHC-derived H-scores, densities, and spatial clustering metrics are available in Table 2.

3.6. High Ki-67 Density and H-Score Are Associated with Resistance to First-Line TT

To examine specific outcomes to first-line therapy, we stratified patients by response and resistance to TT or IT. Patients who were resistant to TT had higher Ki-67 H-scores and higher Ki-67+ cell densities (Wilcoxon p values = 0.013 and 0.035, respectively; Figure 5). None of the other IHC-derived metrics were associated with TT or IT response or resistance in the first-line setting.

3.7. Distinct Geospatial Distributions of CAFs Are Associated with Treatment Response

Examining our cellular marker distributions, we found ccRCC and normal kidney samples to have dramatically different kernel density distributions with respect to their NN distances from each marked cell type (Figure 5A(i,ii)). ccRCC samples had a shorter NN distance from αSMA+ to Ki-67+ cells, compared with normal kidney samples (13 vs. 31 µm, respectively; p < 0.001), and longer NN distance from αSMA+ to caspase-3+ cells (32 vs. 9 µm, respectively; p < 0.001) (Figure 6A(iii)).
TT-responsive and TT-resistant patient samples had significantly different kernel density distributions with respect to their NN distances from each marked cell type (Figure 5B(ii) and Figure 6B(i)). TT-responsive patient samples had a longer NN distance from αSMA+ to caspase-3+ cells (31 vs. 21 µm, respectively; p < 0.001) (Figure 6B(iii)).
IT-responsive and IT-resistant patient samples had significantly different kernel density distributions with respect to their NN distances from each marked cell type (Figure 6C(i,ii)). IT-responsive patient samples had shorter NN distance from αSMA+ to caspase-3+ cells (24 vs. 35 µm, respectively; p < 0.001) (Figure 6C(iii)).

4. Discussion

Our geospatial analysis of CAF distribution demonstrated that close-proximity clustering of CAFs and tumor cells might potentiate tumor cell proliferation. It also shows that a higher density and staining intensity of CAFs and proliferating tumor cells is associated with poor OS and systemic treatment outcomes. To our knowledge, this is the first ex vivo analysis associating CAFs with the potentiation of tumor cell proliferation in metastatic ccRCC.
We found that CAFs were significantly more clustered with proliferating than apoptotic tumor cells. Overall, these findings, as illustrated in Figure 3, show that close tumor cell proximity to CAFs potentiates proliferation. Additionally, this spatial architecture suggests that tumor cells may be in competition with one another for resources provided by CAFs, with the tumor cells further away from CAFs succumbing to apoptosis; however, this analysis was not designed to directly test this hypothesis.
Though this is the first study to directly measure ex vivo clustering of CAFs with proliferating tumor cells in ccRCC, prior studies have demonstrated this effect in vitro, and several ex vivo analyses have discovered pathways that presumably rely on close proximity [1,2,3,7,18]. Notably, CAFs have been found to undergo a metabolic shift when exposed to adjacent cancer cells. This shift results in a Warburg-like glycolytic metabolism in CAFs supplying lactate for cancer cells to fuel the Krebs cycle, which leads to anabolic growth and tumor cell proliferation [7]. Presumably, this is a paracrine effect, with CAF-adjacent tumor cells competing with one another for lactate. This metabolic shift is impacted via the HIF1 pathway, which is of particular relevance in ccRCC, a disease in which the vast majority of tumors harbor a somatic alteration in the VHL gene. The VHL alteration results in dramatically increased intracellular HIFα, which would be expected to further enhance the aforementioned metabolic shift [26]. This well-described effect of CAFs on the tumor microenvironment lends significant biologic plausibility to the geospatial findings identified in this analysis.
Similarly, hypoxia in the tumor microenvironment has been identified as a significant determinant of the extracellular matrix composition in tumors, resulting in increased HIF1A expression. This subsequently increases the production of growth factors that can potentiate tumor cell proliferation and trigger fibroblast activation and fibrosis [27]. In the >80% of ccRCC tumors that harbor somatic VHL mutations, the resulting derangement of HIFα metabolism would be expected to result in a tumor microenvironment consistent with extreme hypoxic conditions. This resultant state would be in place regardless of the true state of oxygen availability within the tissue [26]. This is another example of significant overlap existing between ccRCC and CAF-potentiated molecular pathways.
Additionally, this analysis of a metastatic ccRCC cohort demonstrated strong associations between increased CAF staining intensity and inferior OS while accounting for age and IMDC risk category. This finding corroborates previous IHC studies of ccRCC patients that identified an association between increasing CAF density and OS, as well as a more advanced stage at diagnosis. However, these cohorts did not include stage IV patients [4,5,6].
While using Ki-67 staining, this analysis also identified worse OS associated with increased density of proliferating tumor cells. Though not yet adopted into clinical practice, increased Ki-67 staining has been shown in several studies to be associated with worse survival in ccRCC patients across a wide range of clinical stages [14,15,16,28]. Cell cycle proliferation (CCP) scoring, an RNA-seq-based gene signature score quantifying the expression level of several genes associated with tumor cell proliferation, has been recently described as a promising prognostic biomarker for predicting poor survival and adverse pathology in patients with clinically localized ccRCC [29,30]. Ki-67 IHC staining is a simpler and less costly method than the CCP score and may yield similar biologic information. Further studies are warranted to determine agreement and concordance between Ki-67 IHC staining and CCP score.
In addition to survival metrics, we identified several promising associations between geospatial distributions and treatment response and resistance. Patients who were resistant to TT had increased Ki-67 density and H-scores, which to our knowledge, has not been previously reported. Additionally, two-dimensional kernel-density plots depicting the NN distances of CAFs to their nearest Ki-67+ and caspase-3+ cellular neighbors yielded density distributions that are clearly disparate between TT- and IT-responding and resistant tumors. The difference between these distributions was primarily driven by the distance from CAFs to their nearest caspase-3+ neighbor, with TT responders having significantly longer CAF-caspase NN distance than patients with TT resistance and IT responders having significantly shorter CAF-caspase NN distance than patients with IT resistance. These novel findings suggest the possibility of distinct geospatial CAF architectures associated with response to TT versus IT.
There are several limitations of our study that deserve mentioning. First, the patient cohort studied was heterogeneous and included a variety of systemic treatment agents and was determined partially on the availability of tumor samples. This limited the applicability of our results with contemporary systematic treatment regimens. Secondly, we used a limited number of cellular markers, which likely oversimplified the underlying biology of the tumor microenvironment. There are limits to the practicality of examining an exhaustive list of cellular markers in a novel study such as this and interpreting their possible clinical impact. However, this simplification did allow us to examine these specific markers in a robust fashion and identify possible spatial interactions deserving of future investigations. Additionally, αSMA staining is not 100% specific for CAFs but is a marker for activated fibroblasts that can also stain pericytes and smooth muscle cells. Additionally, CAFs are a heterogeneous group of cells with several distinct subtypes. This study was not designed to assess subclassifications of CAFs. Additionally, non-malignant cells such as infiltrating immune cells can also stain positive for Ki-67, potentially resulting in false-positive staining in some cases. Lastly, we limited our study to using 2 of the more commonly employed spatial metrics, Ripley’s K and NN distance, while acknowledging that many other spatial metrics have been defined in the literature.

5. Conclusions

This ex vivo geospatial analysis of CAF distribution in ccRCC samples suggests that close-proximity clustering of tumor cells and αSMA+CAFs is associated with tumor cell proliferation. Increased tumor cell proliferation was associated with worse OS and resistance to TT. Patients with high αSMA+CAF density and tumor cell proliferation had significantly worse OS from the time of immunotherapy initiation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13153743/s1, Supplemental Figure S1: A. Scatter plots and best-fit lines correlating tumor size with the immunohistochemistry-derived metrics used in the analysis. B. Correlation coefficients and their associated 95% confidence intervals and p values correlating tumor size with each immunohistochemistry-derived metric.

Author Contributions

Analysis planning and supervision, B.J.M., P.M.A., J.J.M., N.C. and J.O.J.; analysis, N.H.C., Y.K., J.O.J. and J.D.; manuscript draft, N.H.C.; critical manuscript editing, G.J.K., A.H., A.C., A.M.A., L.Z., J.C., M.C.F.-F., P.E.S., N.F., E.K., M.F., J.D., J.Y.P., L.W., J.J.M., P.M.A. and B.J.M.; administrative support, E.K. All authors have read and agreed to the published version of the manuscript.

Funding

The described analysis and the APC were directly funded by the Department of Defense United States Army Medical Research Acquisition Activity, grant number KC180036.

Institutional Review Board Statement

All tumor samples were obtained through protocols approved by the institutional review board (H. Lee Moffitt Cancer Center and Research Institute’s Total Cancer Care protocol MCC# 14690; Advarra IRB Pro00014441).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This study was supported by the Molecular Genomics, Tissue Core, and Biostatistics and Bioinformatics Shared Resources at the H. Lee Moffitt Cancer Center and Research Institute, an NCI designated Comprehensive Cancer Center (P30-CA076292). Editorial assistance was provided by Paul Fletcher and Daley Drucker at the Moffitt Cancer Center’s Office of Scientific Publishing. No compensation was given beyond their regular salaries.

Conflicts of Interest

The corresponding author certifies that all conflicts of interest, including specific financial interests and relationships and affiliations relevant to the subject matter or materials discussed in the manuscript (i.e., employment/affiliation, grants or funding, consultancies, honoraria, stock ownership or options, expert testimony, royalties, or patents filed, received, or pending), are the following: G.J.K., Y.K., J.J., N.C., A.H., A.C., A.M.A., L.Z., E.K., J.C., M.C.F.-F., N.F., M.F., J.D., J.Y.P., L.W., and P.M.A. have no disclosures; B.J.M. is an NCCN Kidney Cancer Panel Member; P.E.S. is an NCCN Bladder and Penile Cancer Panel Member and Vice-Chair; J.M. is an Associate Center Director at the Moffitt Cancer Center, has ownership interest in Fulgent Genetics, Inc., Aleta Biotherapeutics, Inc., Cold Genesys, Inc., Myst Pharma, Inc., and Tailored Therapeutics, Inc., and is a consultant/advisory board member for ONCoPEP, Inc., Cold Genesys, Inc., Morphogenesis, Inc., Mersana Therapeutics, Inc., GammaDelta Therapeutics, Ltd., Myst Pharma, Inc., Tailored Therapeutics, Inc., Verseau Therapeutics, Inc., Iovance Biotherapeutics, Inc., Vault Pharma, Inc., Noble Life Sciences Partners, Fulgent Genetics, Inc., UbiVac, LLC, Vycellix, Inc., and Aleta Biotherapeutics, Inc.

References

  1. Kalluri, R. The biology and function of fibroblasts in cancer. Nat. Rev. Cancer 2016, 16, 582–598. [Google Scholar] [CrossRef]
  2. Tao, L.; Huang, G.; Song, H.; Chen, Y.; Chen, L. Cancer associated fibroblasts: An essential role in the tumor microenvironment. Oncol. Lett. 2017, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  3. Harper, J.; Sainson, R.C. Regulation of the anti-tumour immune response by cancer-associated fibroblasts. Semin. Cancer Biol. 2014, 25, 69–77. [Google Scholar] [CrossRef]
  4. Errarte, P.; Larrinaga, G.; Lopez, J.I. The role of cancer-associated fibroblasts in renal cell carcinoma. An example of tumor modulation through tumor/non-tumor cell interactions. J. Adv. Res. 2020, 21, 103–108. [Google Scholar] [CrossRef] [PubMed]
  5. Bakhtyar, N.; Wong, N.; Kapoor, A.; Cutz, J.C.; Hill, B.; Ghert, M.; Tang, D. Clear cell renal cell carcinoma induces fibroblast-mediated production of stromal periostin. Eur. J. Cancer 2013, 49, 3537–3546. [Google Scholar] [CrossRef] [PubMed]
  6. Lopez, J.I.; Errarte, P.; Erramuzpe, A.; Guarch, R.; Cortes, J.M.; Angulo, J.C.; Pulido, R.; Irazusta, J.; Llarena, R.; Larrinaga, G. Fibroblast activation protein predicts prognosis in clear cell renal cell carcinoma. Hum. Pathol. 2016, 54, 100–105. [Google Scholar] [CrossRef] [PubMed]
  7. Fiaschi, T.; Marini, A.; Giannoni, E.; Taddei, M.L.; Gandellini, P.; De Donatis, A.; Lanciotti, M.; Serni, S.; Cirri, P.; Chiarugi, P. Reciprocal metabolic reprogramming through lactate shuttle coordinately influences tumor-stroma interplay. Cancer Res. 2012, 72, 5130–5140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Rini, B.I.; Small, E.J. Biology and clinical development of vascular endothelial growth factor-targeted therapy in renal cell carcinoma. J. Clin. Oncol. 2005, 23, 1028–1043. [Google Scholar] [CrossRef]
  9. Scholzen, T.; Gerdes, J. The Ki-67 protein: From the known and the unknown. J. Cell Physiol. 2000, 182, 311–322. [Google Scholar] [CrossRef]
  10. Li, L.T.; Jiang, G.; Chen, Q.; Zheng, J.N. Ki67 is a promising molecular target in the diagnosis of cancer (review). Mol. Med. Rep. 2015, 11, 1566–1572. [Google Scholar] [CrossRef] [Green Version]
  11. Nielsen, T.O.; Leung, S.C.Y.; Rimm, D.L.; Dodson, A.; Acs, B.; Badve, S.; Denkert, C.; Ellis, M.J.; Fineberg, S.; Flowers, M.; et al. Assessment of Ki67 in Breast Cancer: Updated Recommendations from the International Ki67 in Breast Cancer Working Group. J. Natl. Cancer Inst. 2020. [Google Scholar] [CrossRef]
  12. Nishishita, R.; Morohashi, S.; Seino, H.; Wu, Y.; Yoshizawa, T.; Haga, T.; Saito, K.; Hakamada, K.; Fukuda, S.; Kijima, H. Expression of cancer-associated fibroblast markers in advanced colorectal cancer. Oncol. Lett. 2018, 15, 6195–6202. [Google Scholar] [CrossRef] [Green Version]
  13. Soliman, N.A.; Yussif, S.M. Ki-67 as a prognostic marker according to breast cancer molecular subtype. Cancer Biol. Med. 2016, 13, 496–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Xie, Y.; Chen, L.; Ma, X.; Li, H.; Gu, L.; Gao, Y.; Fan, Y.; Zhang, Y.; Zhang, X. Prognostic and clinicopathological role of high Ki-67 expression in patients with renal cell carcinoma: A systematic review and meta-analysis. Sci. Rep. 2017, 7, 1–9. [Google Scholar] [CrossRef] [Green Version]
  15. Gayed, B.A.; Youssef, R.F.; Bagrodia, A.; Darwish, O.M.; Kapur, P.; Sagalowsky, A.; Lotan, Y.; Margulis, V. Ki67 is an independent predictor of oncological outcomes in patients with localized clear-cell renal cell carcinoma. BJU Int. 2014, 113, 668–673. [Google Scholar] [CrossRef]
  16. Bui, M.H.; Visapaa, H.; Seligson, D.; Kim, H.; Han, K.R.; Huang, Y.; Horvath, S.; Stanbridge, E.J.; Palotie, A.; Figlin, R.A.; et al. Prognostic value of carbonic anhydrase IX and KI67 as predictors of survival for renal clear cell carcinoma. J. Urol. 2004, 171, 2461–2466. [Google Scholar] [CrossRef] [PubMed]
  17. Motzer, R.J.; Jonasch, E.; Michaelson, M.D.; Nandagopal, L.; Gore, J.L.; George, S.; Alva, A.; Haas, N.; Harrison, M.R.; Plimack, E.R.; et al. NCCN Guidelines Insights: Kidney Cancer, Version 2.2020. J. Natl. Compr. Cancer Netw. 2019, 17, 1278–1285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Chen, X.; Song, E. Turning foes to friends: Targeting cancer-associated fibroblasts. Nat. Rev. Drug Discov. 2019, 18, 99–115. [Google Scholar] [CrossRef] [PubMed]
  19. Rizzardi, A.E.; Johnson, A.T.; Vogel, R.I.; Pambuccian, S.E.; Henriksen, J.; Skubitz, A.P.; Metzger, G.J.; Schmechel, S.C. Quantitative comparison of immunohistochemical staining measured by digital image analysis versus pathologist visual scoring. Diagn. Pathol. 2012, 7, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. McClelland, R.A.; Finlay, P.; Walker, K.J.; Nicholson, D.; Robertson, J.F.; Blamey, R.W.; Nicholson, R.I. Automated quantitation of immunocytochemically localized estrogen receptors in human breast cancer. Cancer Res. 1990, 50, 3545–3550. [Google Scholar] [PubMed]
  21. Marcon, E.; Puech, F. A typology of distance-based measures of spatial concentration. Reg. Sci. Urban. Economics 2017, 62, 56–67. [Google Scholar] [CrossRef]
  22. Chakiryan, N.H.; Kimmel, G.J.; Kim, Y.; Hajiran, A.; Aydin, A.M.; Zemp, L.; Katende, E.; Nguyen, J.; Lopez-Blanco, N.; Chahoud, J.; et al. Spatial clustering of CD68+ tumor associated macrophages with tumor cells is associated with worse overall survival in metastatic clear cell renal cell carcinoma. PLoS ONE 2021, 16, e0245415. [Google Scholar] [CrossRef]
  23. Warren, A.Y.; Harrison, D. WHO/ISUP classification, grading and pathological staging of renal cell carcinoma: Standards and controversies. World J. Urol. 2018, 36, 1913–1926. [Google Scholar] [CrossRef] [Green Version]
  24. Heng, D.Y.; Xie, W.; Regan, M.M.; Harshman, L.C.; Bjarnason, G.A.; Vaishampayan, U.N.; Mackenzie, M.; Wood, L.; Donskov, F.; Tan, M.H.; et al. External validation and comparison with other models of the International Metastatic Renal-Cell Carcinoma Database Consortium prognostic model: A population-based study. lancet Oncol. 2013, 14, 141–148. [Google Scholar] [CrossRef] [Green Version]
  25. Faraggi, D.; Simon, R. A simulation study of cross-validation for selecting an optimal cutpoint in univariate survival analysis. Stat. Med. 1996, 15, 2203–2213. [Google Scholar] [CrossRef]
  26. Gnarra, J.R.; Tory, K.; Weng, Y.; Schmidt, L.; Wei, M.H.; Li, H.; Latif, F.; Liu, S.; Chen, F.; Duh, F.M.; et al. Mutations of the VHL tumour suppressor gene in renal carcinoma. Nat. Genet. 1994, 7, 85–90. [Google Scholar] [CrossRef] [PubMed]
  27. Gilkes, D.M.; Semenza, G.L.; Wirtz, D. Hypoxia and the extracellular matrix: Drivers of tumour metastasis. Nat. Rev. Cancer 2014, 14, 430–439. [Google Scholar] [CrossRef] [Green Version]
  28. Mehdi, M.Z.; Nagi, A.H.; Naseem, N. MCM-2 and Ki-67 as proliferation markers in renal cell carcinoma: A quantitative and semi-quantitative analysis. Int. braz j urol 2016, 42, 1121–1128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Morgan, T.M.; Mehra, R.; Tiemeny, P.; Wolf, J.S.; Wu, S.; Sangale, Z.; Brawer, M.; Stone, S.; Wu, C.L.; Feldman, A.S. A Multigene Signature Based on Cell Cycle Proliferation Improves Prediction of Mortality Within 5 Yr of Radical Nephrectomy for Renal Cell Carcinoma. Eur. Urol. 2018, 73, 763–769. [Google Scholar] [CrossRef] [PubMed]
  30. Tosoian, J.J.; Feldman, A.S.; Abbott, M.R.; Mehra, R.; Tiemeny, P.; Wolf, J.S., Jr.; Stone, S.; Wu, S.; Daignault-Newton, S.; Taylor, J.M.G.; et al. Biopsy Cell Cycle Proliferation Score Predicts Adverse Surgical Pathology in Localized Renal Cell Carcinoma. Eur. Urol. 2020, 78, 657–660. [Google Scholar] [CrossRef]
Figure 1. General project workflow for the analysis.
Figure 1. General project workflow for the analysis.
Cancers 13 03743 g001
Figure 2. Geospatial metrics quantifying cell-cell relationships. (A): Boxplot diagrams depicting median cell density (cells/mm2) and H-scores for each of the immunohistochemistry stains used in the analysis (αSMA, Ki-67, and caspase-3). Wilcoxon p values < 0.001 for all comparisons. (B): Boxplot diagrams depicting bivariate spatial clustering, nK(25) and values for each pairwise combination of cells for patient-level data (n = 96 patients). Values > 1.0 reflect clustering of the paired cells, and values < 1.0 reflect dispersion. Wilcoxon p values displayed in the plot. (C): Boxplots depicting median nearest neighbor distances in µm from the indicated pairwise cell combinations for cellular-level data (n > 500,000 cells). Wilcoxon p values are displayed in the plot. (D): Heatmap diagram depicting Spearman’s correlation coefficients between each immunohistochemistry-derived measurement used in the analysis. Note that the deepest blue values reflect a Spearman’s correlation coefficient of −0.4 and the deepest red1.0.
Figure 2. Geospatial metrics quantifying cell-cell relationships. (A): Boxplot diagrams depicting median cell density (cells/mm2) and H-scores for each of the immunohistochemistry stains used in the analysis (αSMA, Ki-67, and caspase-3). Wilcoxon p values < 0.001 for all comparisons. (B): Boxplot diagrams depicting bivariate spatial clustering, nK(25) and values for each pairwise combination of cells for patient-level data (n = 96 patients). Values > 1.0 reflect clustering of the paired cells, and values < 1.0 reflect dispersion. Wilcoxon p values displayed in the plot. (C): Boxplots depicting median nearest neighbor distances in µm from the indicated pairwise cell combinations for cellular-level data (n > 500,000 cells). Wilcoxon p values are displayed in the plot. (D): Heatmap diagram depicting Spearman’s correlation coefficients between each immunohistochemistry-derived measurement used in the analysis. Note that the deepest blue values reflect a Spearman’s correlation coefficient of −0.4 and the deepest red1.0.
Cancers 13 03743 g002
Figure 3. The median tumor architecture identified on immunohistochemistry, derived from the marked cell density and clustering metrics described in Figure 2.
Figure 3. The median tumor architecture identified on immunohistochemistry, derived from the marked cell density and clustering metrics described in Figure 2.
Cancers 13 03743 g003
Figure 4. Kaplan–Meier curves for overall survival, overall survival for the date of immunotherapy initiation, and overall survival from the date of targeted therapy initiation, with groups stratified by (A): αSMA H-score, (B): Ki-67 H-score, and (C): Caspase H-score. Multivariable Cox regression and log-rank p values were reported within the plots.
Figure 4. Kaplan–Meier curves for overall survival, overall survival for the date of immunotherapy initiation, and overall survival from the date of targeted therapy initiation, with groups stratified by (A): αSMA H-score, (B): Ki-67 H-score, and (C): Caspase H-score. Multivariable Cox regression and log-rank p values were reported within the plots.
Cancers 13 03743 g004
Figure 5. Ki-67 H-scores and densities stratified by response to first-line therapy with (A): immunotherapy and (B): targeted therapy. The NA category represents patients that did not receive that respective therapy as a first-line agent.
Figure 5. Ki-67 H-scores and densities stratified by response to first-line therapy with (A): immunotherapy and (B): targeted therapy. The NA category represents patients that did not receive that respective therapy as a first-line agent.
Cancers 13 03743 g005
Figure 6. Scaled kernel density plots (i and ii’s), with the x-axis representing the distance from each αSMA+ cell to its nearest Ki-67+ neighbor (µm) and the y-axis representing the distance from each αSMA+ cell to its nearest caspase-3+ neighbor (µm). Each point on the plot represents an individual αSMA+ cell, and the location of the point corresponds with cells’ NN distance from the nearest Ki-67+ and caspase-3+ neighbor. Two-dimensional kernel density is graphically depicted, overlying the scatter plot, and has been scaled to a range of 0–1.0 in each plot. Dotted red lines reflect median values for the x and y axes. Solid black lines are anchored at fixed values to facilitate visual interpretation between plots. Boxplots (iii’s) compare median nearest-neighbor (NN) distances (µm) from αSMA+ to NN Ki-67+ and αSMA+ to NN caspase-3+ cells, stratified by groups as indicated. (Ai) All ccRCC samples (n = 96), (Aii) Normal kidney samples (n = 10). (Bi) Targeted therapy responders (n = 27). (Bii) Targeted therapy-resistant (n = 9). (Ci) Immunotherapy responders (n = 18). (Cii)-Immunotherapy-resistant (n = 25). (Aiii, Biii and Ciii) Boxplots representing the nearest neighbor distances for each pairwise cell-cell distance for (Aiii) ccRCC verus normal, (Biii) Targeted therapy response versus non-response, and (Ciii) Immunotherapy response versus non-response.
Figure 6. Scaled kernel density plots (i and ii’s), with the x-axis representing the distance from each αSMA+ cell to its nearest Ki-67+ neighbor (µm) and the y-axis representing the distance from each αSMA+ cell to its nearest caspase-3+ neighbor (µm). Each point on the plot represents an individual αSMA+ cell, and the location of the point corresponds with cells’ NN distance from the nearest Ki-67+ and caspase-3+ neighbor. Two-dimensional kernel density is graphically depicted, overlying the scatter plot, and has been scaled to a range of 0–1.0 in each plot. Dotted red lines reflect median values for the x and y axes. Solid black lines are anchored at fixed values to facilitate visual interpretation between plots. Boxplots (iii’s) compare median nearest-neighbor (NN) distances (µm) from αSMA+ to NN Ki-67+ and αSMA+ to NN caspase-3+ cells, stratified by groups as indicated. (Ai) All ccRCC samples (n = 96), (Aii) Normal kidney samples (n = 10). (Bi) Targeted therapy responders (n = 27). (Bii) Targeted therapy-resistant (n = 9). (Ci) Immunotherapy responders (n = 18). (Cii)-Immunotherapy-resistant (n = 25). (Aiii, Biii and Ciii) Boxplots representing the nearest neighbor distances for each pairwise cell-cell distance for (Aiii) ccRCC verus normal, (Biii) Targeted therapy response versus non-response, and (Ciii) Immunotherapy response versus non-response.
Cancers 13 03743 g006
Table 1. Baseline patient, tumor, and treatment characteristics, n = 96 1.
Table 1. Baseline patient, tumor, and treatment characteristics, n = 96 1.
Characteristic.No. (%)
Median age, range59 (55–67)
Gender
Female28 (29)
Male68 (71)
Race
White90 (94)
Black2 (2.1)
Asian1 (1.0)
Other3 (3.1)
IMDC Risk Category
Favorable20 (21)
Intermediate56 (58)
Poor9 (9)
Indeterminate11 (12)
Primary tumor size (cm)8.2 (6.0, 11.0)
ISUP Grade
211 (11)
360 (62)
425 (26)
Sarcomatoid Variant
No87 (91)
Yes9 (9.4)
Rhabdoid Variant
No88 (92)
Yes8 (8.3)
First Line Therapy
IT51 (53)
TT42 (44)
Both3 (3.1)
Response Category
IT Resistance25 (32)
IT Response18 (23)
IT Indeterminate8 (8.3)
TT Resistance9 (9.4)
TT Response27 (28)
TT Indeterminate6 (6.3)
1 Statistics presented: median (IQR); Abbreviations: IMDC, International Metastatic RCC Database Consortium; ISUP, International Society of Urologic Pathologists; IT, immunotherapy; TT, targeted therapy.
Table 2. Multivariable Cox regressions using the covariates age and International Metastatic RCC Database Consortium risk category for the outcomes of (A.) overall survival, (B.) overall survival from the date of immunotherapy initiation, and (C.) overall survival from the date of targeted therapy initiation.
Table 2. Multivariable Cox regressions using the covariates age and International Metastatic RCC Database Consortium risk category for the outcomes of (A.) overall survival, (B.) overall survival from the date of immunotherapy initiation, and (C.) overall survival from the date of targeted therapy initiation.
A.
MetricCutoffN HighN LowSurvival High (mo.)Survival Low (mo.)Cox HR95% CIp
SMA H-score37.028128463.297.61.451.05–2.010.02
Ki 67 H-score7.088356154.797.61.791.40–2.29<0.001
Caspase H-Score0.834603685.868.61.331.01–1.750.04
SMA Density1892.363484885.163.21.170.86–1.600.32
Ki 67 Density1097.034484864.485.11.301.06–1.610.01
Caspase Density205.198474985.868.61.261.02–1.560.03
SMA-Ki67 nK(25)1.167464869.668.60.980.71–1.340.90
SMA-Caspase nK(25)1.106384485.167.60.920.65–1.300.64
Caspase-Ki67 nK(25)1.163424267.689.70.820.58–1.160.26
B.
MetricCutoffN HighN LowSurvival High (mo.)Survival Low (mo.)Cox HR95% CIp
SMA H-score36.346106032.647.61.511.05–2.170.03
Ki 67 H-score7.06274333.356.51.471.13–1.91<0.001
Caspase H-Score1.137135746.746.71.190.89–1.580.24
SMA Density1911.912353546.735.91.240.91–1.680.17
Ki 67 Density1097.034363440.546.71.210.98–1.500.08
Caspase Density189.285353586.140.51.180.95–1.460.13
SMA-Ki67 nK(25)1.139353340.546.70.940.67–1.310.70
SMA-Caspase nK(25)1.085283247.646.70.830.56–1.260.39
Caspase-Ki67 nK(25)1.131313154.140.50.720.47–1.090.12
C.
MetricCutoffN HighN LowSurvival High (mo.)Survival Low (mo.)Cox HR95% CIp
SMA H-score37.2304727.248.91.461.06–2.020.02
Ki 67 H-score7.06106723.155.81.711.35–2.17<0.001
Caspase H-Score0.742502737.233.11.461.14–1.88<0.001
SMA Density1861.339383937.241.51.190.89–1.600.24
Ki 67 Density1107.625383927.241.51.271.04–1.570.02
Caspase Density174.07374036.437.21.251.02–1.530.03
SMA-Ki67 nK(25)1.156373836.437.21.050.76–1.450.76
SMA-Caspase nK(25)1.1313534.237.21.030.76–1.420.84
Caspase-Ki67 nK(25)1.164353333.138.40.920.65–1.310.65
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chakiryan, N.H.; Kimmel, G.J.; Kim, Y.; Johnson, J.O.; Clark, N.; Hajiran, A.; Chang, A.; Aydin, A.M.; Zemp, L.; Katende, E.; et al. Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma. Cancers 2021, 13, 3743. https://doi.org/10.3390/cancers13153743

AMA Style

Chakiryan NH, Kimmel GJ, Kim Y, Johnson JO, Clark N, Hajiran A, Chang A, Aydin AM, Zemp L, Katende E, et al. Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma. Cancers. 2021; 13(15):3743. https://doi.org/10.3390/cancers13153743

Chicago/Turabian Style

Chakiryan, Nicholas H., Gregory J. Kimmel, Youngchul Kim, Joseph O. Johnson, Noel Clark, Ali Hajiran, Andrew Chang, Ahmet M. Aydin, Logan Zemp, Esther Katende, and et al. 2021. "Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma" Cancers 13, no. 15: 3743. https://doi.org/10.3390/cancers13153743

APA Style

Chakiryan, N. H., Kimmel, G. J., Kim, Y., Johnson, J. O., Clark, N., Hajiran, A., Chang, A., Aydin, A. M., Zemp, L., Katende, E., Chahoud, J., Ferrall-Fairbanks, M. C., Spiess, P. E., Francis, N., Fournier, M., Dhillon, J., Park, J. Y., Wang, L., Mulé, J. J., ... Manley, B. J. (2021). Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma. Cancers, 13(15), 3743. https://doi.org/10.3390/cancers13153743

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