Phase and context shape the function of composite oncogenic mutations.

Cancers develop as a result of driver mutations1,2 that lead to clonal outgrowth and the evolution of disease3,4. The discovery and functional characterization of individual driver mutations are central aims of cancer research, and have elucidated myriad phenotypes5 and therapeutic vulnerabilities6. However, the serial genetic evolution of mutant cancer genes7,8 and the allelic context in which they arise is poorly understood in both common and rare cancer genes and tumour types. Here we find that nearly one in four human tumours contains a composite mutation of a cancer-associated gene, defined as two or more nonsynonymous somatic mutations in the same gene and tumour. Composite mutations are enriched in specific genes, have an elevated rate of use of less-common hotspot mutations acquired in a chronology driven in part by oncogenic fitness, and arise in an allelic configuration that reflects context-specific selective pressures. cis-acting composite mutations are hypermorphic in some genes in which dosage effects predominate (such as TERT), whereas they lead to selection of function in other genes (such as TP53). Collectively, composite mutations are driver alterations that arise from context- and allele-specific selective pressures that are dependent in part on gene and mutation function, and which lead to complex-often neomorphic-functions of biological and therapeutic importance.

Lymphatic Invasion of Plakoglobin-Dependent Tumor Cell Clusters Drives Formation of Polyclonal Lung Metastases in Colon Cancer.

BACKGROUND & AIMS: Patients with colon cancer with liver metastases may be cured with surgery, but the presence of additional lung metastases often precludes curative treatment. Little is known about the processes driving lung metastasis. This study aimed to elucidate the mechanisms governing lung vs liver metastasis formation. METHODS: Patient-derived organoid (PDO) cultures were established from colon tumors with distinct patterns of metastasis. Mouse models recapitulating metastatic organotropism were created by implanting PDOs into the cecum wall. Optical barcoding was applied to trace the origin and clonal composition of liver and lung metastases. RNA sequencing and immunohistochemistry were used to identify candidate determinants of metastatic organotropism. Genetic, pharmacologic, in vitro, and in vivo modeling strategies identified essential steps in lung metastasis formation. Validation was performed by analyzing patient-derived tissues. RESULTS: Cecum transplantation of 3 distinct PDOs yielded models with distinct metastatic organotropism: liver only, lung only, and liver and lung. Liver metastases were seeded by single cells derived from select clones. Lung metastases were seeded by polyclonal clusters of tumor cells entering the lymphatic vasculature with very limited clonal selection. Lung-specific metastasis was associated with high expression of desmosome markers, including plakoglobin. Plakoglobin deletion abrogated tumor cell cluster formation, lymphatic invasion, and lung metastasis formation. Pharmacologic inhibition of lymphangiogenesis attenuated lung metastasis formation. Primary human colon, rectum, esophagus, and stomach tumors with lung metastases had a higher N-stage and more plakoglobin-expressing intra-lymphatic tumor cell clusters than those without lung metastases. CONCLUSIONS: Lung and liver metastasis formation are fundamentally distinct processes with different evolutionary bottlenecks, seeding entities, and anatomic routing. Polyclonal lung metastases originate from plakoglobin-dependent tumor cell clusters entering the lymphatic vasculature at the primary tumor site.

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.

AKT mutant allele-specific activation dictates pharmacologic sensitivities.

AKT- a key molecular regulator of PI-3K signaling pathway, is somatically mutated in diverse solid cancer types, and aberrant AKT activation promotes altered cancer cell growth, survival, and metabolism1-8. The most common of AKT mutations (AKT1 E17K) sensitizes affected solid tumors to AKT inhibitor therapy7,8. However, the pathway dependence and inhibitor sensitivity of the long tail of potentially activating mutations in AKT is poorly understood, limiting our ability to act clinically in prospectively characterized cancer patients. Here we show, through population-scale driver mutation discovery combined with functional, biological, and therapeutic studies that some but not all missense mutations activate downstream AKT effector pathways in a growth factor-independent manner and sensitize tumor cells to diverse AKT inhibitors. A distinct class of small in-frame duplications paralogous across AKT isoforms induce structural changes different than those of activating missense mutations, leading to a greater degree of membrane affinity, AKT activation, and cell proliferation as well as pathway dependence and hyper-sensitivity to ATP-competitive, but not allosteric AKT inhibitors. Assessing these mutations clinically, we conducted a phase II clinical trial testing the AKT inhibitor capivasertib (AZD5363) in patients with solid tumors harboring AKT alterations (NCT03310541). Twelve patients were enrolled, out of which six harbored AKT1-3 non-E17K mutations. The median progression free survival (PFS) of capivasertib therapy was 84 days (95% CI 50-not reached) with an objective response rate of 25% (n = 3 of 12) and clinical benefit rate of 42% (n = 5 of 12). Collectively, our data indicate that the degree and mechanism of activation of oncogenic AKT mutants vary, thereby dictating allele-specific pharmacological sensitivities to AKT inhibition.

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.

Respiratory complex and tissue lineage drive recurrent mutations in tumour mtDNA.

Mitochondrial DNA (mtDNA) encodes protein subunits and translational machinery required for oxidative phosphorylation (OXPHOS). Using repurposed whole-exome sequencing data, in the present study we demonstrate that pathogenic mtDNA mutations arise in tumours at a rate comparable to those in the most common cancer driver genes. We identify OXPHOS complexes as critical determinants shaping somatic mtDNA mutation patterns across tumour lineages. Loss-of-function mutations accumulate at an elevated rate specifically in complex I and often arise at specific homopolymeric hotspots. In contrast, complex V is depleted of all non-synonymous mutations, suggesting that impairment of ATP synthesis and mitochondrial membrane potential dissipation are under negative selection. Common truncating mutations and rarer missense alleles are both associated with a pan-lineage transcriptional programme, even in cancer types where mtDNA mutations are comparatively rare. Pathogenic mutations of mtDNA are associated with substantial increases in overall survival of colorectal cancer patients, demonstrating a clear functional relationship between genotype and phenotype. The mitochondrial genome is therefore frequently and functionally disrupted across many cancers, with major implications for patient stratification, prognosis and therapeutic development.

Putative Drivers of Aggressiveness in TCEB1-mutant Renal Cell Carcinoma: An Emerging Entity with Variable Clinical Course.

BACKGROUND: TCEB1-mutant renal cell carcinoma (RCC) is a rare variant of RCC with clear-cell features. Owing to its unique morphological and molecular features it has recently been proposed as a separate entity. Initial series suggested an indolent, early-stage phenotype. Here we expand our clinical cohort and describe a more detailed genomic analysis looking for potential drivers of aggressiveness. DESIGN, SETTING, AND PARTICIPANTS: We identified five new cases in our institutional sequencing cohort, four of whom were found to have high-stage disease (American Joint Committee on Cancer stage III/IV). Twelve previously reported cases were pooled for comparison purposes (Sato, The Cancer Genome Atlas, TRACERx Renal). OUTCOME MEASURES AND STATISTICAL ANALYSIS: We used our previously validated pipeline to analyze somatic mutations and copy number alterations (CNAs) in seven tumor samples with available data and estimated the number of cancer cells bearing each somatic mutation. The oncogenic potential of mutations was assessed using OncoKB and two other algorithms. Mann-Whitney U tests were used to evaluate differences in genomic markers between stage groups. RESULTS AND LIMITATIONS: All tumors showed biallelic inactivation of the TCEB1 gene according to a combination of somatic mutation and CNA analyses. Mutations were always found in residues involved in hydrophobic interactions with VHL. We found that high-stage tumors had additional oncogenic mutations (median 1, interquartile range [IQR] 1-1 vs 2, IQR 2-2; median difference 1, 95% confidence interval [CI] 1-1; p= 0.002) and showed whole-genome doubling events. They also seemed to have a higher burden of somatic CNAs (median fraction CNA genome 0.10, IQR 0.10-0.15 vs 0.63, IQR 0.58-0.68), however, this finding did not reach statistical significance (median difference 0.49, 95% CI 0.33-0.63; p=0.052). CONCLUSIONS: TCEB1-mutant RCC can show variable behavior ranging from very indolent to aggressive. Specific molecular events leading to high genomic instability seem to be associated with aggressiveness. This study expands the clinical spectrum of TCEB1-mutant RCC. PATIENT SUMMARY: We present four cases of aggressive TCEB1-mutant renal cell carcinoma, a rare type of kidney cancer. In-depth analysis of the genomes of these tumors revealed certain abnormalities that might explain this aggressive behavior.

TRK xDFG Mutations Trigger a Sensitivity Switch from Type I to II Kinase Inhibitors.

On-target resistance to next-generation TRK inhibitors in TRK fusion-positive cancers is largely uncharacterized. In patients with these tumors, we found that TRK xDFG mutations confer resistance to type I next-generation TRK inhibitors designed to maintain potency against several kinase domain mutations. Computational modeling and biochemical assays showed that TRKAG667 and TRKCG696 xDFG substitutions reduce drug binding by generating steric hindrance. Concurrently, these mutations stabilize the inactive (DFG-out) conformations of the kinases, thus sensitizing these kinases to type II TRK inhibitors. Consistently, type II inhibitors impede the growth and TRK-mediated signaling of xDFG-mutant isogenic and patient-derived models. Collectively, these data demonstrate that adaptive conformational resistance can be abrogated by shifting kinase engagement modes. Given the prior identification of paralogous xDFG resistance mutations in other oncogene-addicted cancers, these findings provide insights into rational type II drug design by leveraging inhibitor class affinity switching to address recalcitrant resistant alterations. SIGNIFICANCE: In TRK fusion-positive cancers, TRK xDFG substitutions represent a shared liability for type I TRK inhibitors. In contrast, they represent a potential biomarker of type II TRK inhibitor activity. As all currently available type II agents are multikinase inhibitors, rational drug design should focus on selective type II inhibitor creation.This article is highlighted in the In This Issue feature, p. 1.

Leveraging Systematic Functional Analysis to Benchmark an In Silico Framework Distinguishes Driver from Passenger MEK Mutants in Cancer.

Despite significant advances in cancer precision medicine, a significant hurdle to its broader adoption remains the multitude of variants of unknown significance identified by clinical tumor sequencing and the lack of biologically validated methods to distinguish between functional and benign variants. Here we used functional data on MAP2K1 and MAP2K2 mutations generated in real-time within a co-clinical trial framework to benchmark the predictive value of a three-part in silico methodology. Our computational approach to variant classification incorporated hotspot analysis, three-dimensional molecular dynamics simulation, and sequence paralogy. In silico prediction accurately distinguished functional from benign MAP2K1 and MAP2K2 mutants, yet drug sensitivity varied widely among activating mutant alleles. These results suggest that multifaceted in silico modeling can inform patient accrual to MEK/ERK inhibitor clinical trials, but computational methods need to be paired with laboratory- and clinic-based efforts designed to unravel variabilities in drug response. SIGNIFICANCE: Leveraging prospective functional characterization of MEK1/2 mutants, it was found that hotspot analysis, molecular dynamics simulation, and sequence paralogy are complementary tools that can robustly prioritize variants for biologic, therapeutic, and clinical validation.See related commentary by Whitehead and Sebolt-Leopold, p. 4042.

Efficacy and Determinants of Response to HER Kinase Inhibition in HER2-Mutant Metastatic Breast Cancer.

HER2 mutations define a subset of metastatic breast cancers with a unique mechanism of oncogenic addiction to HER2 signaling. We explored activity of the irreversible pan-HER kinase inhibitor neratinib, alone or with fulvestrant, in 81 patients with HER2-mutant metastatic breast cancer. Overall response rate was similar with or without estrogen receptor (ER) blockade. By comparison, progression-free survival and duration of response appeared longer in ER+ patients receiving combination therapy, although the study was not designed for direct comparison. Preexistent concurrent activating HER2 or HER3 alterations were associated with poor treatment outcome. Similarly, acquisition of multiple HER2-activating events, as well as gatekeeper alterations, were observed at disease progression in a high proportion of patients deriving clinical benefit from neratinib. Collectively, these data define HER2 mutations as a therapeutic target in breast cancer and suggest that coexistence of additional HER signaling alterations may promote both de novo and acquired resistance to neratinib. SIGNIFICANCE: HER2 mutations define a targetable breast cancer subset, although sensitivity to irreversible HER kinase inhibition appears to be modified by the presence of concurrent activating genomic events in the pathway. These findings have implications for potential future combinatorial approaches and broader therapeutic development for this genomically defined subset of breast cancer.This article is highlighted in the In This Issue feature, p. 161.

V211D Mutation in MEK1 Causes Resistance to MEK Inhibitors in Colon Cancer.

We report the emergence of the novel MEK1 V211D gatekeeper mutation in a patient with BRAF K601E colon cancer treated with the allosteric MEK inhibitor binimetinib and the anti-EGFR antibody panitumumab. The MEK1 V211D mutation concurrently occurs in the same cell with BRAF K601E and leads to RAF-independent activity but remains regulated by RAF. The V211D mutation causes resistance to binimetinib by both increasing the catalytic activity of MEK1 and reducing its affinity for the drug. Moreover, the mutant exhibits reduced sensitivity to all the allosteric MEK inhibitors tested. Thus, this mutation serves as a general resistance mutation for current MEK inhibitors; however, it is sensitive to a newly reported ATP-competitive MEK inhibitor, which therefore could be used to overcome drug resistance. SIGNIFICANCE: We report a resistance mechanism to allosteric MEK inhibitors in the clinic. A MEK1 V211D mutation developed in a patient with BRAF K601E colon cancer on MEK and EGFR inhibitors. This mutant increases the catalytic activity of MEK1 and reduces its affinity for binimetinib, but remains sensitive to ATP-competitive MEK inhibitors.This article is highlighted in the In This Issue feature, p. 1143.

Altered Nuclear Export Signal Recognition as a Driver of Oncogenesis.

Altered expression of XPO1, the main nuclear export receptor in eukaryotic cells, has been observed in cancer, and XPO1 has been a focus of anticancer drug development. However, mechanistic evidence for cancer-specific alterations in XPO1 function is lacking. Here, genomic analysis of 42,793 cancers identified recurrent and previously unrecognized mutational hotspots in XPO1. XPO1 mutations exhibited striking lineage specificity, with enrichment in a variety of B-cell malignancies, and introduction of single amino acid substitutions in XPO1 initiated clonal, B-cell malignancy in vivo. Proteomic characterization identified that mutant XPO1 altered the nucleocytoplasmic distribution of hundreds of proteins in a sequence-specific manner that promoted oncogenesis. XPO1 mutations preferentially sensitized cells to inhibitors of nuclear export, providing a biomarker of response to this family of drugs. These data reveal a new class of oncogenic alteration based on change-of-function mutations in nuclear export signal recognition and identify therapeutic targets based on altered nucleocytoplasmic trafficking. SIGNIFICANCE: Here, we identify that heterozygous mutations in the main nuclear exporter in eukaryotic cells, XPO1, are positively selected in cancer and promote the initiation of clonal B-cell malignancies. XPO1 mutations alter nuclear export signal recognition in a sequence-specific manner and sensitize cells to compounds in clinical development inhibiting XPO1 function.This article is highlighted in the In This Issue feature, p. 1325.