Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing for a detailed description of the somatic mutational landscape of ccRCC. We identify candidate driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for therapeutic interventions. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The observations that higher T-cell infiltration is associated with better overall survival and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

Multi-State Gene Cluster Switches Determine the Adaptive Mitochondrial And Metabolic Landscape of Breast Cancer.

Adaptive metabolic switches are proposed to underlie conversions between cellular states during normal development as well as in cancer evolution. Metabolic adaptations represent important therapeutic targets in tumors, highlighting the need to characterize the full spectrum, characteristics, and regulation of the metabolic switches. To investigate the hypothesis that metabolic switches associated with specific metabolic states can be recognized by locating large alternating gene expression patterns, we developed a method to identify interspersed gene sets by massive correlated biclustering (MCbiclust) and to predict their metabolic wiring. Testing the method on breast cancer transcriptome datasets revealed a series of gene sets with switch-like behavior that could be used to predict mitochondrial content, metabolic activity, and central carbon flux in tumors. The predictions were experimentally validated by bioenergetic profiling and metabolic flux analysis of 13C-labelled substrates. The metabolic switch positions also distinguished between cellular states, correlating with tumor pathology, prognosis, and chemosensitivity. The method is applicable to any large and heterogeneous transcriptome dataset to discover metabolic and associated pathophysiological states.

Genetic landscape of interval and screen detected breast cancer.

Interval breast cancers (IBCs) are cancers diagnosed between screening episodes. Understanding the biological differences between IBCs and screen-detected breast-cancers (SDBCs) has the potential to improve mammographic screening and patient management. We analysed and compared the genomic landscape of 288 IBCs and 473 SDBCs by whole genome sequencing of paired tumour-normal patient samples collected as part of the UK 100,000 Genomes Project. Compared to SDBCs, IBCs were more likely to be lobular, higher grade, and triple negative. A more aggressive clinical phenotype was reflected in IBCs displaying features of genomic instability including a higher mutation rate and number of chromosomal structural abnormalities, defective homologous recombination and TP53 mutations. We did not however, find evidence to indicate that IBCs are associated with a significantly different immune response. While IBCs do not represent a unique molecular class of invasive breast cancer they exhibit a more aggressive phenotype, which is likely to be a consequence of the timing of tumour initiation. This information is relevant both with respect to treatment as well as informing the screening interval for mammography.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing the most detailed somatic mutational landscape to date. We identify new driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for drug repurposing. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The twin observations that higher T-cell infiltration is associated with better outcome and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

The evolution of lung cancer and impact of subclonal selection in TRACERx.

Lung cancer is the leading cause of cancer-associated mortality worldwide1. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource.

Using DNA sequencing data to quantify T cell fraction and therapy response.

The immune microenvironment influences tumour evolution and can be both prognostic and predict response to immunotherapy1,2. However, measurements of tumour infiltrating lymphocytes (TILs) are limited by a shortage of appropriate data. Whole-exome sequencing (WES) of DNA is frequently performed to calculate tumour mutational burden and identify actionable mutations. Here we develop T cell exome TREC tool (T cell ExTRECT), a method for estimation of T cell fraction from WES samples using a signal from T cell receptor excision circle (TREC) loss during V(D)J recombination of the T cell receptor-α gene (TCRA (also known as TRA)). TCRA T cell fraction correlates with orthogonal TIL estimates and is agnostic to sample type. Blood TCRA T cell fraction is higher in females than in males and correlates with both tumour immune infiltrate and presence of bacterial sequencing reads. Tumour TCRA T cell fraction is prognostic in lung adenocarcinoma. Using a meta-analysis of tumours treated with immunotherapy, we show that tumour TCRA T cell fraction predicts immunotherapy response, providing value beyond measuring tumour mutational burden. Applying T cell ExTRECT to a multi-sample pan-cancer cohort reveals a high diversity of the degree of immune infiltration within tumours. Subclonal loss of 12q24.31-32, encompassing SPPL3, is associated with reduced TCRA T cell fraction. T cell ExTRECT provides a cost-effective technique to characterize immune infiltrate alongside somatic changes.

Meta-analysis of tumor- and T cell-intrinsic mechanisms of sensitization to checkpoint inhibition.

Checkpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

The breast cancer oncogene IKKε coordinates mitochondrial function and serine metabolism.

IκB kinase ε (IKKε) is a key molecule at the crossroads of inflammation and cancer. Known to regulate cytokine secretion via NFκB and IRF3, the kinase is also a breast cancer oncogene, overexpressed in a variety of tumours. However, to what extent IKKε remodels cellular metabolism is currently unknown. Here, we used metabolic tracer analysis to show that IKKε orchestrates a complex metabolic reprogramming that affects mitochondrial metabolism and consequently serine biosynthesis independently of its canonical signalling role. We found that IKKε upregulates the serine biosynthesis pathway (SBP) indirectly, by limiting glucose-derived pyruvate utilisation in the TCA cycle, inhibiting oxidative phosphorylation. Inhibition of mitochondrial function induces activating transcription factor 4 (ATF4), which in turn drives upregulation of the expression of SBP genes. Importantly, pharmacological reversal of the IKKε-induced metabolic phenotype reduces proliferation of breast cancer cells. Finally, we show that in a highly proliferative set of ER negative, basal breast tumours, IKKε and PSAT1 are both overexpressed, corroborating the link between IKKε and the SBP in the clinical context.

Macrophages induce malignant traits in mammary epithelium via IKKε/TBK1 kinases and the serine biosynthesis pathway.

During obesity, macrophages infiltrate the breast tissue leading to low-grade chronic inflammation, a factor considered responsible for the higher risk of breast cancer associated with obesity. Here, we formally demonstrate that breast epithelial cells acquire malignant properties when exposed to medium conditioned by macrophages derived from human healthy donors. These effects were mediated by the breast cancer oncogene IKKε and its downstream target-the serine biosynthesis pathway as demonstrated by genetic or pharmacological tools. Furthermore, amlexanox, an FDA-approved drug targeting IKKε and its homologue TBK1, delayed in vivo tumour formation in a combined genetic mouse model of breast cancer and high-fat diet-induced obesity/inflammation. Finally, in human breast cancer tissues, we validated the link between inflammation-IKKε and alteration of cellular metabolism. Altogether, we identified a pathway connecting obesity-driven inflammation to breast cancer and a potential therapeutic strategy to reduce the risk of breast cancer associated with obesity.

MCbiclust: a novel algorithm to discover large-scale functionally related gene sets from massive transcriptomics data collections.

The potential to understand fundamental biological processes from gene expression data has grown in parallel with the recent explosion of the size of data collections. However, to exploit this potential, novel analytical methods are required, capable of discovering large co-regulated gene networks. We found current methods limited in the size of correlated gene sets they could discover within biologically heterogeneous data collections, hampering the identification of multi-gene controlled fundamental cellular processes such as energy metabolism, organelle biogenesis and stress responses. Here we describe a novel biclustering algorithm called Massively Correlated Biclustering (MCbiclust) that selects samples and genes from large datasets with maximal correlated gene expression, allowing regulation of complex networks to be examined. The method has been evaluated using synthetic data and applied to large bacterial and cancer cell datasets. We show that the large biclusters discovered, so far elusive to identification by existing techniques, are biologically relevant and thus MCbiclust has great potential in the analysis of transcriptomics data to identify large-scale unknown effects hidden within the data. The identified massive biclusters can be used to develop improved transcriptomics based diagnosis tools for diseases caused by altered gene expression, or used for further network analysis to understand genotype-phenotype correlations.

The mitochondrial calcium uniporter regulates breast cancer progression via HIF-1α.

Triple-negative breast cancer (TNBC) represents the most aggressive breast tumor subtype. However, the molecular determinants responsible for the metastatic TNBC phenotype are only partially understood. We here show that expression of the mitochondrial calcium uniporter (MCU), the selective channel responsible for mitochondrial Ca(2+) uptake, correlates with tumor size and lymph node infiltration, suggesting that mitochondrial Ca(2+) uptake might be instrumental for tumor growth and metastatic formation. Accordingly, MCU downregulation hampered cell motility and invasiveness and reduced tumor growth, lymph node infiltration, and lung metastasis in TNBC xenografts. In MCU-silenced cells, production of mitochondrial reactive oxygen species (mROS) is blunted and expression of the hypoxia-inducible factor-1α (HIF-1α) is reduced, suggesting a signaling role for mROS and HIF-1α, downstream of mitochondrial Ca(2+) Finally, in breast cancer mRNA samples, a positive correlation of MCU expression with HIF-1α signaling route is present. Our results indicate that MCU plays a central role in TNBC growth and metastasis formation and suggest that mitochondrial Ca(2+) uptake is a potential novel therapeutic target for clinical intervention.