Genomic ancestry in kidney cancer: Correlations with clinical and molecular features.

INTRODUCTION: Genetic ancestry (GA) refers to population hereditary patterns that contribute to phenotypic differences seen among race/ethnicity groups, and differences among GA groups may highlight unique biological determinants that add to our understanding of health care disparities. METHODS: A retrospective review of patients with renal cell carcinoma (RCC) was performed and correlated GA with clinicopathologic, somatic, and germline molecular data. All patients underwent next-generation sequencing of normal and tumor DNA using Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets, and contribution of African (AFR), East Asian (EAS), European (EUR), Native American, and South Asian (SAS) ancestry was inferred through supervised ADMIXTURE. Molecular data was compared across GA groups by Fisher exact test and Kruskal-Wallis test. RESULTS: In 953 patients with RCC, the GA distribution was: EUR (78%), AFR (4.9%), EAS (2.5%), SAS (2%), Native American (0.2%), and Admixed (12.2%). GA distribution varied by tumor histology and international metastatic RCC database consortium disease risk status (intermediate-poor: EUR 58%, AFR 88%, EAS 74%, and SAS 73%). Pathogenic/likely pathogenic germline variants in cancer-predisposition genes varied (16% EUR, 23% AFR, 8% EAS, and 0% SAS), and most occurred in CHEK2 in EUR (3.1%) and FH in AFR (15.4%). In patients with clear cell RCC, somatic alteration incidence varied with significant enrichment in BAP1 alterations (EUR 17%, AFR 50%, SAS 29%; p = .01). Comparing AFR and EUR groups within The Cancer Genome Atlas, significant differences were identified in angiogenesis and inflammatory pathways. CONCLUSION: Differences in clinical and molecular data by GA highlight population-specific variations in patients with RCC. Exploration of both genetic and nongenetic variables remains critical to optimize efforts to overcome health-related disparities.

Characterizing the Immune Microenvironment and Neoantigen Landscape of Hürthle Cell Carcinoma to Identify Potential Immunologic Vulnerabilities.

Hürthle cell carcinoma (HCC) is a rare type of thyroid cancer with high rates of distant metastasis and recurrence. Along with the scarcity of effective systemic therapies for HCC, these factors contribute to poor clinical outcomes. The immunologic features of HCC are poorly defined and response rates with immune checkpoint blockade have not been reported. A more comprehensive understanding of the immune landscape and factors that predict response to checkpoint inhibitors is needed. We performed RNA sequencing on 40 tumors to characterize the neoantigen landscape and immune microenvironment of HCC. We analyzed transcriptomic profiles, tumor-infiltrating immune cell populations, and measures of T-cell activation/dysfunction and correlated these to genetic features such as tumor mutation burden, neoantigen burden, mitochondrial mutations, and LOH from chromosomal uniparental disomy. Finally, immune profiles of patients with recurrence were compared with those of patients without recurrence. HCC tumors exhibited low levels of immune infiltration, with the more aggressive widely invasive phenotype associated with more immune depletion. There was a negative correlation between tumor mutation burden, neoantigen burden, programmed cell death ligand 1 (PD-L1) expression, and the immune infiltration score. HCC tumors that exhibited a global LOH from chromosomal uniparental disomy or haploidization had the lowest level of immune infiltration. HCC tumors that recurred displayed an immune-depleted microenvironment associated with global LOH and aerobic glycolysis. These findings offer new insights into the functional immune landscapes and immune microenvironment of HCC. Our data identify potential immunologic vulnerabilities for these understudied and often fatal cancers. Significance: The immune landscape of HCC is poorly defined and response rates to immunotherapy have not been reported. The authors found the immune microenvironment in HCC to be depleted. This immunosuppression is associated with a global LOH from haploidization and uniparental disomy, resulting in whole chromosome losses across the genome.

Nivolumab plus ipilimumab in advanced salivary gland cancer: a phase 2 trial.

Salivary gland cancers (SGCs) are rare, aggressive cancers without effective treatments when metastasized. We conducted a phase 2 trial evaluating nivolumab (nivo, anti-PD-1) and ipilimumab (ipi, anti-CTLA-4) in 64 patients with metastatic SGC enrolled in two histology-based cohorts (32 patients each): adenoid cystic carcinoma (ACC; cohort 1) and other SGCs (cohort 2). The primary efficacy endpoint (≥4 objective responses) was met in cohort 2 (5/32, 16%) but not in cohort 1 (2/32, 6%). Treatment safety/tolerability and progression-free survival (PFS) were secondary endpoints. Treatment-related adverse events grade ≥3 occurred in 24 of 64 (38%) patients across both cohorts, and median PFS was 4.4 months (95% confidence interval (CI): 2.4, 8.3) and 2.2 months (95% CI: 1.8, 5.3) for cohorts 1 and 2, respectively. We present whole-exome, RNA and T cell receptor (TCR) sequencing data from pre-treatment and on-treatment tumors and immune cell flow cytometry and TCR sequencing from peripheral blood at serial timepoints. Responding tumors universally demonstrated clonal expansion of pre-existing T cells and mutational contraction. Responding ACCs harbored neoantigens, including fusion-derived neoepitopes, that induced T cell responses ex vivo. This study shows that nivo+ipi has limited efficacy in ACC, albeit with infrequent, exceptional responses, and that it could be promising for non-ACC SGCs, particularly salivary duct carcinomas. ClinicalTrials.gov identifier: NCT03172624 .

Immunometabolic coevolution defines unique microenvironmental niches in ccRCC.

Tumor cell phenotypes and anti-tumor immune responses are shaped by local metabolite availability, but intratumoral metabolite heterogeneity (IMH) and its phenotypic consequences remain poorly understood. To study IMH, we profiled tumor/normal regions from clear cell renal cell carcinoma (ccRCC) patients. A common pattern of IMH transcended all patients, characterized by correlated fluctuations in the abundance of metabolites and processes associated with ferroptosis. Analysis of intratumoral metabolite-RNA covariation revealed that the immune composition of the microenvironment, especially the abundance of myeloid cells, drove intratumoral metabolite variation. Motivated by the strength of RNA-metabolite covariation and the clinical significance of RNA biomarkers in ccRCC, we inferred metabolomic profiles from the RNA sequencing data of ccRCC patients enrolled in 7 clinical trials, and we ultimately identifyied metabolite biomarkers associated with response to anti-angiogenic agents. Local metabolic phenotypes, therefore, emerge in tandem with the immune microenvironment, influence ongoing tumor evolution, and are associated with therapeutic sensitivity.

Activation of the EGFR/PI3K/AKT pathway limits the efficacy of trametinib treatment in head and neck cancer.

Blocking the mitogen-activated protein kinase (MAPK) pathway with the MEK1/2 inhibitor trametinib has produced promising results in patients with head and neck squamous cell carcinoma (HNSCC). In the current study, we showed that trametinib treatment leads to overexpression and activation of the epidermal growth factor receptor (EGFR) in HNSCC cell lines and patient-derived xenografts. Knockdown of EGFR improved trametinib treatment efficacy both in vitro and in vivo. Mechanistically, we demonstrated that trametinib-induced EGFR overexpression hyperactivates the phosphatidylinositol 3-kinase (PI3K)/AKT pathway. In vitro, blocking the PI3K pathway with GDC-0941 (pictilisib), or BYL719 (alpelisib), prevented AKT pathway hyperactivation and enhanced the efficacy of trametinib in a synergistic manner. In vivo, a combination of trametinib and BYL719 showed superior antitumor efficacy vs. the single agents, leading to tumor growth arrest. We confirmed our findings in a syngeneic murine head and neck cancer cell line in vitro and in vivo. Taken together, our findings show that trametinib treatment induces hyperactivation of EGFR/PI3K/AKT; thus, blocking of the EGFR/PI3K pathway is required to improve trametinib efficacy in HNSCC.

Patterns of immune equilibrium and escape in indolent and progressing tumors.

By comparing indolent/slowly progressing with aggressive/rapidly progressing tumor types, Pandey et al. identify human evidence of immune equilibrium in indolent tumors and immune escape in progressing tumors, suggesting a link between these mechanisms and the epidemiologic phenomenon of overdiagnosis.

A multimodal atlas of tumour metabolism reveals the architecture of gene-metabolite covariation.

Tumour metabolism is controlled by coordinated changes in metabolite abundance and gene expression, but simultaneous quantification of metabolites and transcripts in primary tissue is rare. To overcome this limitation and to study gene-metabolite covariation in cancer, we assemble the Cancer Atlas of Metabolic Profiles of metabolomic and transcriptomic data from 988 tumour and control specimens spanning 11 cancer types in published and newly generated datasets. Meta-analysis of the Cancer Atlas of Metabolic Profiles reveals two classes of gene-metabolite covariation that transcend cancer types. The first corresponds to gene-metabolite pairs engaged in direct enzyme-substrate interactions, identifying putative genes controlling metabolite pool sizes. A second class of gene-metabolite covariation represents a small number of hub metabolites, including quinolinate and nicotinamide adenine dinucleotide, which correlate to many genes specifically expressed in immune cell populations. These results provide evidence that gene-metabolite covariation in cellularly heterogeneous tissue arises, in part, from both mechanistic interactions between genes and metabolites, and from remodelling of the bulk metabolome in specific immune microenvironments.

Geospatial characterization of immune cell distributions and dynamics across the microenvironment in clear cell renal cell carcinoma.

INTRODUCTION: In clear cell renal cell carcinoma (ccRCC), tumor-associated macrophage (TAM) induction of CD8+T cells into a terminally exhausted state has been implicated as a major mechanism of immunotherapy resistance, but a deeper biological understanding is necessary. METHODS: Primary ccRCC tumor samples were obtained from 97 patients between 2004 and 2018. Multiplex immunofluorescence using lymphoid and myeloid markers was performed in seven regions of interest per patient across three predefined zones, and geospatial analysis was performed using Ripley's K analysis, a methodology adapted from ecology. RESULTS: Clustering of CD163+M2 like TAMs into the stromal compartment at the tumor-stroma interface was associated with worse clinical stage (tumor/CD163+nK(75): stage I/II: 4.4 (IQR -0.5 to 5.1); stage III: 1.4 (IQR -0.3 to 3.5); stage IV: 0.6 (IQR -2.1 to 2.1); p=0.04 between stage I/II and stage IV), and worse overall survival (OS) and cancer-specific survival (CSS) (tumor/CD163+nK(75): median OS-hi=149 months, lo=86 months, false-discovery rate (FDR)-adj. Cox p<0.001; median CSS-hi=174 months, lo=85 months; FDR-adj. Cox p<0.001). An RNA-seq differential gene expression score was developed using this geospatial metric, and was externally validated in multiple independent cohorts of patients with ccRCC including: TCGA KIRC, and the IMmotion151, IMmotion150, and JAVELIN Renal 101 clinical trials. In addition, this CD163+ geospatial pattern was found to be associated with a higher TIM-3+ proportion of CD8+T cells, indicative of terminal exhaustion (tumor-core: 0.07 (IQR 0.04-0.14) vs 0.40 (IQR 0.15-0.66), p=0.05). CONCLUSIONS: Geospatial clustering of CD163+M2 like TAMs into the stromal compartment at the tumor-stromal interface was associated with poor clinical outcomes and CD8+T cell terminal exhaustion.

Novel Genetic Subtypes of Urothelial Carcinoma With Differential Outcomes on Immune Checkpoint Blockade.

PURPOSE: Immune checkpoint blockade (ICB) therapy has significantly improved clinical outcomes in bladder cancer. Identification of correlates of benefit is critical to select appropriate therapy for individual patients. METHODS: To reveal genetic variables associated with benefit from ICB, we performed whole-exome sequencing on tumor specimens from 88 patients with advanced bladder cancer treated with ICB. RESULTS: We identified several genetic factors that correlated with progression-free and overall survival after ICB therapy including ARID1A mutation, tumor mutational burden, intratumoral heterogeneity, the ratio of nonsynonymous to synonymous mutations in the immunopeptidome (immune dN/dS), and tumor cell purity. In addition, we noted that neutrophil-to-lymphocyte ratio and smoking history were negatively associated with overall survival. These genetic characteristics define four molecular subtypes demonstrating differential sensitivity to ICB. We validated the association of these four subtypes with clinical benefit from ICB in an independent cohort (IMvigor210). Finally, we showed that these molecular subtypes also correlate with outcome, although with distinct relationships, among patients not treated with ICB using The Cancer Genome Atlas (TCGA) bladder cancer cohort. Using parallel RNA sequencing data, the subtypes were also shown to correlate with immune infiltration and inflammation, respectively, in the IMvigor210 and TCGA cohorts. CONCLUSION: Together, our study defines molecular subgroups of bladder cancer that influence benefit from ICB.

Genomic and Transcriptomic Characteristics of Metastatic Thyroid Cancers with Exceptional Responses to Radioactive Iodine Therapy.

PURPOSE: The determinants of response or resistance to radioiodine (RAI) are unknown. We aimed to identify genomic and transcriptomic factors associated with structural responses to RAI treatment of metastatic thyroid cancer, which occur infrequently, and to test whether high MAPK pathway output was associated with RAI refractoriness. EXPERIMENTAL DESIGN: Exceptional response to RAI was defined as reduction of tumor volume based on RECIST v1.1. We performed a retrospective case-control study of genomic and transcriptomic characteristics of exceptional responders (ER; n = 8) versus nonresponders (NR; n = 16) matched by histologic type and stage at presentation on a 1:2 ratio. RESULTS: ER are enriched for mutations that activate MAPK through RAF dimerization (RAS, class 2 BRAF, RTK fusions), whereas NR are associated with BRAFV600E, which signals as a monomer and is unresponsive to negative feedback. ER have a lower MAPK transcriptional output and a higher thyroid differentiation score (TDS) than NR (P < 0.05). NR are enriched for 1q-gain (P < 0.05) and mutations of genes regulating mRNA splicing and the PI3K pathway. BRAFV600E tumors with 1q-gain have a lower TDS than BRAFV600E/1q-quiet tumors and transcriptomic signatures associated with metastatic propensity. CONCLUSIONS: ER tumors have a lower MAPK output and higher TDS than NR, whereas NR have a high frequency of BRAFV600E and 1q-gain. Molecular profiling of thyroid cancers and further functional validation of the key findings discriminating ER from NR may help predict response to RAI therapy.

Spatiotemporal evolution of the clear cell renal cell carcinoma microenvironment links intra-tumoral heterogeneity to immune escape.

BACKGROUND: Intratumoral heterogeneity (ITH) is a hallmark of clear cell renal cell carcinoma (ccRCC) that reflects the trajectory of evolution and influences clinical prognosis. Here, we seek to elucidate how ITH and tumor evolution during immune checkpoint inhibitor (ICI) treatment can lead to therapy resistance. METHODS: Here, we completed a single-arm pilot study to examine the safety and feasibility of neoadjuvant nivolumab in patients with localized RCC. Primary endpoints were safety and feasibility of neoadjuvant nivolumab. Then, we spatiotemporally profiled the genomic and immunophenotypic characteristics of 29 ccRCC patients, including pre- and post-therapy samples from 17 ICI-treated patients. Deep multi-regional whole-exome and transcriptome sequencing were performed on 29 patients at different time points before and after ICI therapy. T cell repertoire was also monitored from tissue and peripheral blood collected from a subset of patients to study T cell clonal expansion during ICI therapy. RESULTS: Angiogenesis, lymphocytic infiltration, and myeloid infiltration varied significantly across regions of the same patient, potentially confounding their utility as biomarkers of ICI response. Elevated ITH associated with a constellation of both genomic features (HLA LOH, CDKN2A/B loss) and microenvironmental features, including elevated myeloid expression, reduced peripheral T cell receptor (TCR) diversity, and putative neoantigen depletion. Hypothesizing that ITH may itself play a role in shaping ICI response, we derived a transcriptomic signature associated with neoantigen depletion that strongly associated with response to ICI and targeted therapy treatment in several independent clinical trial cohorts. CONCLUSIONS: These results argue that genetic and immune heterogeneity jointly co-evolve and influence response to ICI in ccRCC. Our findings have implications for future biomarker development for ICI response across ccRCC and other solid tumors and highlight important features of tumor evolution under ICI treatment. TRIAL REGISTRATION: The study was registered on ClinicalTrial.gov (NCT02595918) on November 4, 2015.

FOXA1 repression drives lineage plasticity and immune heterogeneity in bladder cancers with squamous differentiation.

Cancers arising from the bladder urothelium often exhibit lineage plasticity with regions of urothelial carcinoma adjacent to or admixed with regions of divergent histomorphology, most commonly squamous differentiation. To define the biologic basis for and clinical significance of this morphologic heterogeneity, here we perform integrated genomic analyses of mixed histology bladder cancers with separable regions of urothelial and squamous differentiation. We find that squamous differentiation is a marker of intratumoral genomic and immunologic heterogeneity in patients with bladder cancer and a biomarker of intrinsic immunotherapy resistance. Phylogenetic analysis confirms that in all cases the urothelial and squamous regions are derived from a common shared precursor. Despite the presence of marked genomic heterogeneity between co-existent urothelial and squamous differentiated regions, no recurrent genomic alteration exclusive to the urothelial or squamous morphologies is identified. Rather, lineage plasticity in bladder cancers with squamous differentiation is associated with loss of expression of FOXA1, GATA3, and PPARG, transcription factors critical for maintenance of urothelial cell identity. Of clinical significance, lineage plasticity and PD-L1 expression is coordinately dysregulated via FOXA1, with patients exhibiting morphologic heterogeneity pre-treatment significantly less likely to respond to immune checkpoint inhibitors.

Associations between Pretreatment Body Composition Features and Clinical Outcomes among Patients with Metastatic Clear Cell Renal Cell Carcinoma Treated with Immune Checkpoint Blockade.

PURPOSE: High body mass index (BMI) may lead to improved immune-checkpoint blockade (ICB) outcomes in metastatic clear cell renal cell carcinoma (mccRCC). However, BMI is a crude body size measure. We investigated BMI and radiographically assessed body composition (BC) parameters association with mccRCC ICB outcomes. EXPERIMENTAL DESIGN: Retrospective study of ICB-treated patients with mccRCC. BMI and BC variables [skeletal muscle index (SMI) and multiple adiposity indexes] were determined using pretreatment CT scans. We examined the associations between BMI and BC variables with ICB outcomes. Therapeutic responses per RECIST v1.1 were determined. We compared whole-transcriptomic patterns with BC variables in a separate cohort of 62 primary tumor samples. RESULTS: 205 patients with mccRCC were included in the cohort (74% were male, 71% were overweight/obese, and 53% were classified as low SMI). High-BMI patients experienced longer overall survival (OS) than normal-weight patients [unadjusted HR, 0.66; 95% confidence interval (CI), 0.45-0.97; P = 0.035]. The only BC variable associated with OS was SMI [unadjusted HR comparing low vs. high SMI 1.65 (95% CI: 1.13-2.43); P = 0.009]. However, this OS association became nonsignificant after adjusting for International Metastatic Renal Cell Carcinoma Database Consortium score and line of therapy. No OS association was seen for adiposity and no BC variable was associated with progression-free survival or radiological responses. Tumors from patients with low SMI displayed increased angiogenic, inflammatory, and myeloid signals. CONCLUSIONS: Our findings highlight the relevance of skeletal muscle in the BMI paradox. Future studies should investigate if addressing low skeletal muscle in metastatic patients treated with ICB can improve survival.

Tumor-associated macrophages expressing the transcription factor IRF8 promote T cell exhaustion in cancer.

Tumors are populated by antigen-presenting cells (APCs) including macrophage subsets with distinct origins and functions. Here, we examined how cancer impacts mononuclear phagocytic APCs in a murine model of breast cancer. Tumors induced the expansion of monocyte-derived tumor-associated macrophages (TAMs) and the activation of type 1 dendritic cells (DC1s), both of which expressed and required the transcription factor interferon regulatory factor-8 (IRF8). Although DC1s mediated cytotoxic T lymphocyte (CTL) priming in tumor-draining lymph nodes, TAMs promoted CTL exhaustion in the tumor, and IRF8 was required for TAMs' ability to present cancer cell antigens. TAM-specific IRF8 deletion prevented exhaustion of cancer-cell-reactive CTLs and suppressed tumor growth. Tumors from patients with immune-infiltrated renal cell carcinoma had abundant TAMs that expressed IRF8 and were enriched for an IRF8 gene expression signature. Furthermore, the TAM-IRF8 signature co-segregated with CTL exhaustion signatures across multiple cancer types. Thus, CTL exhaustion is promoted by TAMs via IRF8.

Immune Checkpoint Blockade Outcome in Small-Cell Lung Cancer and Its Relationship With Retinoblastoma Mutation Status and Function.

PURPOSE: Immune checkpoint blockade (ICB) in conjunction with chemotherapy is approved for the treatment of extensive-stage small-cell lung cancer (SCLC). Although specific genomic abnormalities such as KEAP1 and STK11 gene mutations are associated with resistance to ICB in non-SCLC, no genomic abnormality has been found in association with resistance to ICB in SCLC. MATERIALS AND METHODS: We first analyzed a retrospective cohort of 42 patients with SCLC treated with single-agent ICB or ICB combination (data set A). We then validated our results in a large prospective clinical trial of 460 patients (CheckMate 032, data set B). DNA and RNA sequencing were performed. RESULTS: In data set A, patients treated with ICB with RB1 wild-type (WT) had a median overall survival (OS) of 23.1 months (95% CI, 9 to 37.5), whereas the RB1 mutant OS was 5 months (95% CI, 2.5 to 26; P = .04). Differentially expressed gene analysis between RB1 mutant and RB1 WT samples indicated the enrichment of downregulated immune-related genes and an immune exclusion phenotype among RB1 mutant but not in the RB1 WT tumor samples. We then assessed results from 460 patients enrolled in CheckMate 032, a trial of nivolumab (NIVO) or NIVO + ipilimumab only in SCLC. In this large cohort, RB1 WT patients had significantly improved outcome with NIVO therapy compared with mutant patients (hazard ratio, 1.41; 95% CI, 1.02 to 2.01; P = .041). High RB1 loss-of-function (LOF) signature scores significantly associated with neuroendocrine subtypes (ASCL1 and NeuroD1). However, neuroendocrine subtypes did not associate with OS. Remarkably, patients with lower RB1 LOF scores had longer OS following treatment with NIVO. CONCLUSION: SCLC patients with RB1 WT status or lower RB1 LOF signature scores by transcriptomics have better outcomes with ICB monotherapy.

Functional landscapes of POLE and POLD1 mutations in checkpoint blockade-dependent antitumor immunity.

Defects in pathways governing genomic fidelity have been linked to improved response to immune checkpoint blockade therapy (ICB). Pathogenic POLE/POLD1 mutations can cause hypermutation, yet how diverse mutations in POLE/POLD1 influence antitumor immunity following ICB is unclear. Here, we comprehensively determined the effect of POLE/POLD1 mutations in ICB and elucidated the mechanistic impact of these mutations on tumor immunity. Murine syngeneic tumors harboring Pole/Pold1 functional mutations displayed enhanced antitumor immunity and were sensitive to ICB. Patients with POLE/POLD1 mutated tumors harboring telltale mutational signatures respond better to ICB than patients harboring wild-type or signature-negative tumors. A mutant POLE/D1 function-associated signature-based model outperformed several traditional approaches for identifying POLE/POLD1 mutated patients that benefit from ICB. Strikingly, the spectrum of mutational signatures correlates with the biochemical features of neoantigens. Alterations that cause POLE/POLD1 function-associated signatures generate T cell receptor (TCR)-contact residues with increased hydrophobicity, potentially facilitating T cell recognition. Altogether, the functional landscapes of POLE/POLD1 mutations shape immunotherapy efficacy.

A Targetable Myeloid Inflammatory State Governs Disease Recurrence in Clear-Cell Renal Cell Carcinoma.

It is poorly understood how the tumor immune microenvironment influences disease recurrence in localized clear-cell renal cell carcinoma (ccRCC). Here we performed whole-transcriptomic profiling of 236 tumors from patients assigned to the placebo-only arm of a randomized, adjuvant clinical trial for high-risk localized ccRCC. Unbiased pathway analysis identified myeloid-derived IL6 as a key mediator. Furthermore, a novel myeloid gene signature strongly correlated with disease recurrence and overall survival on uni- and multivariate analyses and is linked to TP53 inactivation across multiple data sets. Strikingly, effector T-cell gene signatures, infiltration patterns, and exhaustion markers were not associated with disease recurrence. Targeting immunosuppressive myeloid inflammation with an adenosine A2A receptor antagonist in a novel, immunocompetent, Tp53-inactivated mouse model significantly reduced metastatic development. Our findings suggest that myeloid inflammation promotes disease recurrence in ccRCC and is targetable as well as provide a potential biomarker-based framework for the design of future immuno-oncology trials in ccRCC. SIGNIFICANCE: Improved understanding of factors that influence metastatic development in localized ccRCC is greatly needed to aid accurate prediction of disease recurrence, clinical decision-making, and future adjuvant clinical trial design. Our analysis implicates intratumoral myeloid inflammation as a key driver of metastasis in patients and a novel immunocompetent mouse model. This article is highlighted in the In This Issue feature, p. 2221.

Mitonuclear genotype remodels the metabolic and microenvironmental landscape of Hürthle cell carcinoma.

Hürthle cell carcinomas (HCCs) display two exceptional genotypes: near-homoplasmic mutation of mitochondrial DNA (mtDNA) and genome-wide loss of heterozygosity (gLOH). To understand the phenotypic consequences of these genetic alterations, we analyzed genomic, metabolomic, and immunophenotypic data of HCC and other thyroid cancers. Both mtDNA mutations and profound depletion of citrate pools are common in HCC and other thyroid malignancies, suggesting that thyroid cancers are broadly equipped to survive tricarboxylic acid cycle impairment, whereas metabolites in the reduced form of NADH-dependent lysine degradation pathway were elevated exclusively in HCC. The presence of gLOH was not associated with metabolic phenotypes but rather with reduced immune infiltration, indicating that gLOH confers a selective advantage partially through immunosuppression. Unsupervised multimodal clustering revealed four clusters of HCC with distinct clinical, metabolomic, and microenvironmental phenotypes but overlapping genotypes. These findings chart the metabolic and microenvironmental landscape of HCC and shed light on the interaction between genotype, metabolism, and the microenvironment in cancer.