Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes.

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Obesity shapes selection for driver mutations in cancer.

Obesity is a leading risk factor for cancer, but whether obesity is linked to specific genomic subtypes of cancer is unknown. Here, we examined the relationship between obesity and tumor genotype in two large clinicogenomic corpora. Obesity was associated with specific driver mutations in lung adenocarcinoma, endometrial carcinoma, and cancers of unknown primary, independent of clinical covariates and genetic ancestry. Obesity is therefore a putative driver of etiologic heterogeneity across cancers.

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.

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.

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.

Clinical sequencing of soft tissue and bone sarcomas delineates diverse genomic landscapes and potential therapeutic targets.

The genetic, biologic, and clinical heterogeneity of sarcomas poses a challenge for the identification of therapeutic targets, clinical research, and advancing patient care. Because there are > 100 sarcoma subtypes, in-depth genetic studies have focused on one or a few subtypes. Herein, we report a comparative genetic analysis of 2,138 sarcomas representing 45 pathological entities. This cohort is prospectively analyzed using targeted sequencing to characterize subtype-specific somatic alterations in targetable pathways, rates of whole genome doubling, mutational signatures, and subtype-agnostic genomic clusters. The most common alterations are in cell cycle control and TP53, receptor tyrosine kinases/PI3K/RAS, and epigenetic regulators. Subtype-specific associations include TERT amplification in intimal sarcoma and SWI/SNF alterations in uterine adenosarcoma. Tumor mutational burden, while low compared to other cancers, varies between and within subtypes. This resource will improve sarcoma models, motivate studies of subtype-specific alterations, and inform investigations of genetic factors and their correlations with treatment response.

Genomic and Metabolic Hallmarks of SDH- and FH-deficient Renal Cell Carcinomas.

BACKGROUND: Succinate dehydrogenase-deficient and fumarate hydratase-deficient renal cell carcinomas (SDHRCC and FHRCC) are rare kidney cancers driven by loss of TCA cycle enzymes. OBJECTIVE: To define and compare the genomic and metabolomic hallmarks of SDHRCC and FHRCC. DESIGN, SETTING, AND PARTICIPANTS: We analyzed SDHRCC and FHRCC tumors with either immunohistochemical evidence of loss of protein expression or genomically confirmed biallelic inactivation of SDHA/B/C/D/AF2 or FH. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Somatic alterations were identified using clinical pipelines, with allele-specific copy number alterations (CNAs) identified using FACETS. Mass spectrometry-based metabolomic profiling was performed on available SDHRCC and FHRCC tumors. RESULTS AND LIMITATIONS: Tumors were analyzed for 42 patients (25 FHRCC, 17 SDHRCC). In the germline analysis, 16/17 SDHRCCs harbored a germline alteration in SDHB, whereas only 17/22 FHRCCs had pathogenic germline FH variants. SDHRCCs had a lower mutation burden (p = 0.02) and CNA burden (p = 0.0002) than FHRCCs. All SDHRCCs presented with deletion of chromosome 1p (overlapping SDHB), whereas FHRCCs demonstrated high but not ubiquitous loss of 1q (FH locus). Both SDHRCCs and FHRCCs exhibited significant idiopathic accumulation of the metabolite guanine. FHRCC tumors had elevated levels of urea cycle metabolites (argininosuccinate, citrulline, and fumarate), whereas SDHRCC tumors had elevation of numerous acylcarnitines. These characteristic metabolic changes allowed identification of a previously unrecognized SDH-deficient RCC. CONCLUSIONS: Despite sharing similar genetic etiology, SDHRCC and FHRCC represent distinct molecular entities with unique genetic and metabolic abnormalities. PATIENT SUMMARY: Kidney cancers driven by loss of the gene encoding either the succinate dehydrogenase or fumarate hydratase enzyme are rare. We sought to define and compare the genetic and metabolic features of these cancer entities.

A Single Session of Aerobic Exercise Mediates Plasticity-Related Phosphorylation in both the Rat Motor Cortex and Hippocampus.

A single session of aerobic exercise may offer one means to "prime" motor regions to be more receptive to the acquisition of a motor skill; however, the mechanisms whereby this priming may occur are not clear. One possible explanation may be related to the post-translational modification of plasticity-related receptors and their associated intracellular signaling molecules, given that these proteins are integral to the development of synaptic plasticity. In particular, phosphorylation governs the biophysical properties (e.g., Ca2+ conductance) and the migratory patterns (i.e., trafficking) of plasticity-related receptors by altering the relative density of specific receptor subunits at synapses. We hypothesized that a single session of exercise would alter the subunit phosphorylation of plasticity-related receptors (AMPA receptors, NMDA receptors) and signaling molecules (PKA, CaMKII) in a manner that would serve to prime motor cortex. Young, male Sprague-Dawley rats (n = 24) were assigned to either exercise (Moderate, Exhaustion), or non-exercising (Sedentary) groups. Immediately following a single session of treadmill exercise, whole tissue homogenates were prepared from both the motor cortex and hippocampus. We observed a robust (1.2-2.0× greater than sedentary) increase in tyrosine phosphorylation of AMPA (GluA1,2) and NMDA (GluN2A,B) receptor subunits, and a clear indication that exercise preferentially affects pPKA over pCaMKII. The changes were found, specifically, following moderate, but not maximal, acute aerobic exercise in both motor cortex and hippocampus. Given the requirement for these proteins during the early phases of plasticity induction, the possibility exists that exercise-induced priming may occur by altering the phosphorylation of plasticity-related proteins.