Kuo, Fengshen

Abstract Background 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. In vitro mechanistic studies have demonstrated that the anti-tumor activity of lymphoid and myeloid cell populations is mediated by metabolite availability and signaling in the TME, raising the possibility that the immune response and metabolism of ccRCC tumors coevolve and jointly influence the likelihood that a patient responds to therapy.. However, both the broad patterns of coordination between metabolite abundance and TME cellular composition, as well as the precise cell populations producing metabolic phenotypes of interest, remain unknown. Methods To study IMH, we multiregionally profiled the metabolome, transcriptome, and genome of 187 tumor/normal regions from 31 clear cell renal cell carcinoma (ccRCC) patients. Using these measurements and additional multimodal metabolomic/transcriptomic profiling of ccRCC and other diseases, we developed computational models that can be used to understand RNA-metabolite covariation and ultimately impute metabolite levels from RNA sequencing data. Results Analysis of intratumoral metabolite-RNA covariation revealed that the immune composition of the microenvironment, and 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 deployed and benchmarked a machine learning method (MIRTH) to impute metabolite levels directly from RNA sequencing data of primary and metastatic ccRCC tumors. We inferred metabolomic profiles from RNA sequencing data of ccRCC patients enrolled in 6 clinical trials, ultimately identifying specific metabolite biomarkers associated with response to anti-angiogenic agents. Conclusions Local metabolic phenotypes therefore emerge in tandem with the immune microenvironment and associate with therapeutic sensitivity. CDMRP DOD Funding: yes

Abstract 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.

Abstract Background Patients with human papillomavirus–related oropharyngeal cancers have excellent outcomes but experience clinically significant toxicities when treated with standard chemoradiotherapy (70 Gy). We hypothesized that functional imaging could identify patients who could be safely deescalated to 30 Gy of radiotherapy. Methods In 19 patients, pre- and intratreatment dynamic fluorine-18-labeled fluoromisonidazole positron emission tomography (PET) was used to assess tumor hypoxia. Patients without hypoxia at baseline or intratreatment received 30 Gy; patients with persistent hypoxia received 70 Gy. Neck dissection was performed at 4 months in deescalated patients to assess pathologic response. Magnetic resonance imaging (weekly), circulating plasma cell-free DNA, RNA-sequencing, and whole-genome sequencing (WGS) were performed to identify potential molecular determinants of response. Samples from an independent prospective study were obtained to reproduce molecular findings. All statistical tests were 2-sided. Results Fifteen of 19 patients had no hypoxia on baseline PET or resolution on intratreatment PET and were deescalated to 30 Gy. Of these 15 patients, 11 had a pathologic complete response. Two-year locoregional control and overall survival were 94.4% (95% confidence interval = 84.4% to 100%) and 94.7% (95% confidence interval = 85.2% to 100%), respectively. No acute grade 3 radiation–related toxicities were observed. Microenvironmental features on serial imaging correlated better with pathologic response than tumor burden metrics or circulating plasma cell-free DNA. A WGS-based DNA repair defect was associated with response (P = .02) and was reproduced in an independent cohort (P = .03). Conclusions Deescalation of radiotherapy to 30 Gy on the basis of intratreatment hypoxia imaging was feasible, safe, and associated with minimal toxicity. A DNA repair defect identified by WGS was predictive of response. Intratherapy personalization of chemoradiotherapy may facilitate marked deescalation of radiotherapy.