37Integrative multi-omic characterization of the immune landscape in renal cell carcinoma
Bibliographic record
Abstract
Abstract Background Renal Cell Carcinoma (RCC) has been characterized as being amongst the most immune infiltrated solid tumors with a highly heterogenous immune landscape. Within spatially organized cellular networks (CNs) of the tumor immune microenvironment (TIME), key immune cell-cell interactions (CCIs) impact immune cell function and organization ultimately impacting the patient’s overall response. Multiple studies have observed an association of tertiary lymphoid structures, a commonly observed spatial CN, with positive clinical outcomes in RCC, however additional CNs and the CCIs that control these networks need to be identified to better harness and potentially reprogram immune responses to improve patient outcomes. The recent explosion of interest in the heterogeneity of the immune landscape in RCC has led to numerous publications using the latest technologies in spatial transcriptomics, proteomics, and metabolomics. However, many of these studies use this data in isolation and therefore, may be hindered by the technological biases inherent in each method. Here, we have developed a novel approach to integrate spatial and single cell multi-omic data harnessing the strengths of each technology to better interrogate CNs that exist in the RCC TIME. Methods Fresh surgical samples were procured at the University Health Network (Toronto, Canada) through the REnal cancer MicroEnvironment DiscoverY (REMEDY) study. Bulk RNA sequencing (RNA-seq), single cell RNA sequencing (scRNA-seq), single cell suspension mass cytometry (SMC), whole transcriptome digital spatial profiling (DSP), and imaging mass cytometry (IMC) was performed on spatially concordant tumor regions across 54 patients. scRNA-seq enabled the identification of immune, stromal, and malignant high-resolution cell states, which informed a tailored antibody panel design for SMC and IMC and served as a reference framework for harmonized cell-type annotation across modalities. This enabled the integration of our transcriptomic and proteomic data to delineate RCC-specific CNs enriched for defined CCIs across unique biological pathways. Results Using this integrative approach, we identified seven high-resolution patient immunophenotypes. To evaluate their prognostic and predictive relevance, we derived representative gene signatures using a linear mixed model (Flash-MM) to interrogate publicly available bulk RNA-seq datasets, including TCGA, JAVELIN, and IMmotion151. This analysis revealed immunophenotype-specific associations with survival following surgery or systemic therapy in univariate models. To explore potential biological mechanisms associated with these divergent clinical outcomes, we incorporated spatial information into traditional CCI analyses and performed pathway analyses to define functional relationships. This revealed that patient subtypes with high lymphoid infiltration exhibit greater spatial heterogeneity, potentially reflecting the coexistence of multiple activated immune pathways. In contrast, patient subtypes with low immune infiltration were enriched in more developmental signaling pathways. In addition to our biological observations, we were also able to assess the technological differences between patient matched samples and compare the ability of each technology to capture inter-patient and intra-patient heterogeneity. Conclusions Collectively, this work identifies clinically distinct subgroups defined by CNs and details the interpatient cellular heterogeneity that exists in RCC, providing the foundation for future personalized therapeutic interventions against this disease.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".