Francesco Favero

Clonal status of actionable driver events and the timing of mutational processes in cancer evolution

Pan-cancer analysis of tumor subclonal structures uncovers widespread intratumor heterogeneity in actionable driver events and reveals mutational processes fueling their selection. Uncovering a tumor’s family tree In the era of targeted anticancer drugs, correctly identifying the mutations in a tumor becomes an essential part of optimizing cancer treatment. This is not necessarily straightforward because tumors can contain both “driver” mutations, which control tumor growth and therefore should be blocked with specific drugs, and “passenger” mutations, which, as their name suggests, may not contribute to the progression of a tumor and are unlikely to be useful therapeutic targets. McGranahan et al. identified patterns of driver event evolution in a wide variety of tumor types, revealing specific patterns of mutations that will be important in the design of future therapeutic regimens for cancer. Deciphering whether actionable driver mutations are found in all or a subset of tumor cells will likely be required to improve drug development and precision medicine strategies. We analyzed nine cancer types to determine the subclonal frequencies of driver events, to time mutational processes during cancer evolution, and to identify drivers of subclonal expansions. Although mutations in known driver genes typically occurred early in cancer evolution, we also identified later subclonal “actionable” mutations, including BRAF (V600E), IDH1 (R132H), PIK3CA (E545K), EGFR (L858R), and KRAS (G12D), which may compromise the efficacy of targeted therapy approaches. More than 20% of IDH1 mutations in glioblastomas, and 15% of mutations in genes in the PI3K (phosphatidylinositol 3-kinase)–AKT–mTOR (mammalian target of rapamycin) signaling axis across all tumor types were subclonal. Mutations in the RAS–MEK (mitogen-activated protein kinase kinase) signaling axis were less likely to be subclonal than mutations in genes associated with PI3K-AKT-mTOR signaling. Analysis of late mutations revealed a link between APOBEC-mediated mutagenesis and the acquisition of subclonal driver mutations and uncovered putative cancer genes involved in subclonal expansions, including CTNNA2 and ATXN1. Our results provide a pan-cancer census of driver events within the context of intratumor heterogeneity and reveal patterns of tumor evolution across cancers. The frequent presence of subclonal driver mutations suggests the need to stratify targeted therapy response according to the proportion of tumor cells in which the driver is identified.

MicroRNA Signature of Malignant Mesothelioma with Potential Diagnostic and Prognostic Implications

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of target genes, and may behave as oncogenes or tumor suppressors. Human malignant mesothelioma is an asbestos-related cancer, with poor prognosis and low median survival. Here we report, for the first time, a cross-evaluation of miRNA expression in mesothelioma (MPP-89, REN) and human mesothelial cells (HMC-telomerase reverse transcriptase). Microarray profiling, confirmed by real-time quantitative RT-PCR, revealed a differential expression of miRNAs between mesothelioma and mesothelial cells. In addition, a computational analysis combining miRNA and gene expression profiles allowed the accurate prediction of genes potentially targeted by dysregulated miRNAs. Several predicted genes belong to terms of Gene Ontology (GO) that are associated with the development and progression of mesothelioma. This suggests that miRNAs may be key players in mesothelioma oncogenesis. We further investigated miRNA expression on a panel of 24 mesothelioma specimens, representative of the three histotypes (epithelioid, biphasic, and sarcomatoid), by quantitative RT-PCR. The expression of miR-17-5p, miR-21, miR-29a, miR-30c, miR-30e-5p, miR-106a, and miR-143 was significantly associated with the histopathological subtypes. Notably, the reduced expression of two miRNAs (miR-17-5p and miR-30c) correlated with better survival of patients with sarcomatoid subtype. Our preliminary analysis points at miRNAs as potential diagnostic and prognostic markers of mesothelioma, and suggests novel tools for the therapy of this malignancy.

Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data

We describe our algorithm and software for determining copy number profiles from tumor genome sequencing data, and find that it compares favorably to existing algorithms for the same purpose.

Targeting Human Cancer by a Glycosaminoglycan Binding Malaria Protein.

Plasmodium falciparum engineer infected erythrocytes to present the malarial protein, VAR2CSA, which binds a distinct type chondroitin sulfate (CS) exclusively expressed in the placenta. Here, we show that the same CS modification is present on a high proportion of malignant cells and that it can be specifically targeted by recombinant VAR2CSA (rVAR2). In tumors, placental-like CS chains are linked to a limited repertoire of cancer-associated proteoglycans including CD44 and CSPG4. The rVAR2 protein localizes to tumors in vivo and rVAR2 fused to diphtheria toxin or conjugated to hemiasterlin compounds strongly inhibits in vivo tumor cell growth and metastasis. Our data demonstrate how an evolutionarily refined parasite-derived protein can be exploited to target a common, but complex, malignancy-associated glycosaminoglycan modification.

Reduced Expression of the ROCK Inhibitor Rnd3 Is Associated with Increased Invasiveness and Metastatic Potential in Mesenchymal Tumor Cells

Background Mesenchymal and amoeboid movements are two important mechanisms adopted by cancer cells to invade the surrounding environment. Mesenchymal movement depends on extracellular matrix protease activity, amoeboid movement on the RhoA-dependent kinase ROCK. Cancer cells can switch from one mechanism to the other in response to different stimuli, limiting the efficacy of antimetastatic therapies. Methodology and Principal Findings We investigated the acquisition and molecular regulation of the invasion capacity of neoplastically transformed human fibroblasts, which were able to induce sarcomas and metastases when injected into immunocompromised mice. We found that neoplastic transformation was associated with a change in cell morphology (from fibroblastic to polygonal), a reorganization of the actin cytoskeleton, a decrease in the expression of several matrix metalloproteases and increases in cell motility and invasiveness. In a three-dimensional environment, sarcomagenic cells showed a spherical morphology with cortical actin rings, suggesting a switch from mesenchymal to amoeboid movement. Accordingly, cell invasion decreased after treatment with the ROCK inhibitor Y27632, but not with the matrix protease inhibitor Ro 28-2653. The increased invasiveness of tumorigenic cells was associated with reduced expression of Rnd3 (also known as RhoE), a cellular inhibitor of ROCK. Indeed, ectopic Rnd3 expression reduced their invasive ability in vitro and their metastatic potential in vivo. Conclusions These results indicate that, during neoplastic transformation, cells of mesenchymal origin can switch from a mesenchymal mode of movement to an amoeboid one. In addition, they point to Rnd3 as a possible regulator of mesenchymal tumor cell invasion and to ROCK as a potential therapeutic target for sarcomas.

Development of synchronous VHL syndrome tumors reveals contingencies and constraints to tumor evolution

BackgroundGenomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution.ResultsWe perform whole exome sequencing on four clear cell renal cell carcinomas removed from both kidneys of a patient with a germline VHL mutation. We report that tumors arising in this context are clonally independent and harbour distinct secondary events exemplified by loss of chromosome 3p, despite an identical genetic background and tissue microenvironment. We propose that divergent mutational and copy number anomalies are contingent upon the nature of 3p loss of heterozygosity occurring early in tumorigenesis. However, despite distinct 3p events, genomic, proteomic and immunohistochemical analyses reveal evidence for convergence upon the PI3K-AKT-mTOR signaling pathway. Four germline tumors in this young patient, and in a second, older patient with VHL syndrome demonstrate minimal intra-tumor heterogeneity and mutational burden, and evaluable tumors appear to follow a linear evolutionary route, compared to tumors from patients with sporadic clear cell renal cell carcinoma.ConclusionsIn tumors developing from a germline VHL mutation, the evolutionary principles of contingency and convergence in tumor development are complementary. In this small set of patients with early stage VHL-associated tumors, there is reduced mutation burden and limited evidence of intra-tumor heterogeneity.

Pan-cancer analysis of genomic scar signatures associated with homologous recombination deficiency suggests novel indications for existing cancer drugs.

BackgroundOvarian and triple-negative breast cancers with BRCA1 or BRCA2 loss are highly sensitive to treatment with PARP inhibitors and platinum-based cytotoxic agents and show an accumulation of genomic scars in the form of gross DNA copy number aberrations. Cancers without BRCA1 or BRCA2 loss but with accumulation of similar genomic scars also show increased sensitivity to platinum-based chemotherapy. Therefore, reliable biomarkers to identify DNA repair-deficient cancers prior to treatment may be useful for directing patients to platinum chemotherapy and possibly PARP inhibitors. Recently, three SNP array-based signatures of chromosomal instability were published that each quantitate a distinct type of genomic scar considered likely to be caused by improper DNA repair. They measure telomeric allelic imbalance (named NtAI), large scale transition (named LST), and loss of heterozygosity (named HRD-LOH), and it is suggested that these signatures may act as biomarkers for the state of DNA repair deficiency in a given cancer.ResultsWe explored the pan-cancer distribution of scores of the three signatures utilizing a panel of 5371 tumors representing 15 cancer types from The Cancer Genome Atlas, and found a good correlation between scores of the three signatures (Spearman’s ρ 0.73–0.87). In addition we found that cancer types ordinarily receiving platinum as standard of care have higher median scores of all three signatures. Interestingly, we also found that smaller subpopulations of high-scoring tumors exist in most cancer types, including those for which platinum chemotherapy is not standard therapy.ConclusionsWithin several cancer types that are not ordinarily treated with platinum chemotherapy, we identified tumors with high levels of the three genomic biomarkers. These tumors represent identifiable subtypes of patients which may be strong candidates for clinical trials with PARP inhibitors or platinum-based chemotherapeutic regimens.

BCL9L Dysfunction Impairs Caspase-2 Expression Permitting Aneuploidy Tolerance in Colorectal Cancer

Summary Chromosomal instability (CIN) contributes to cancer evolution, intratumor heterogeneity, and drug resistance. CIN is driven by chromosome segregation errors and a tolerance phenotype that permits the propagation of aneuploid genomes. Through genomic analysis of colorectal cancers and cell lines, we find frequent loss of heterozygosity and mutations in BCL9L in aneuploid tumors. BCL9L deficiency promoted tolerance of chromosome missegregation events, propagation of aneuploidy, and genetic heterogeneity in xenograft models likely through modulation of Wnt signaling. We find that BCL9L dysfunction contributes to aneuploidy tolerance in both TP53-WT and mutant cells by reducing basal caspase-2 levels and preventing cleavage of MDM2 and BID. Efforts to exploit aneuploidy tolerance mechanisms and the BCL9L/caspase-2/BID axis may limit cancer diversity and evolution.

Glioblastoma adaptation traced through decline of an IDH1 clonal driver and macro-evolution of a double-minute chromosome

In a glioblastoma tumour with multi-region sequencing before and after recurrence, we find an IDH1 mutation that is clonal in the primary but lost at recurrence. We also describe the evolution of a double-minute chromosome encoding regulators of the PI3K signalling axis that dominates at recurrence, emphasizing the challenges of an evolving and dynamic oncogenic landscape for precision medicine.

Mitochondrial mutations drive prostate cancer aggression

Nuclear mutations are well known to drive tumor incidence, aggression and response to therapy. By contrast, the frequency and roles of mutations in the maternally inherited mitochondrial genome are poorly understood. Here we sequence the mitochondrial genomes of 384 localized prostate cancer patients, and identify a median of one mitochondrial single-nucleotide variant (mtSNV) per patient. Some of these mtSNVs occur in recurrent mutational hotspots and associate with aggressive disease. Younger patients have fewer mtSNVs than those who diagnosed at an older age. We demonstrate strong links between mitochondrial and nuclear mutational profiles, with co-occurrence between specific mutations. For example, certain control region mtSNVs co-occur with gain of the MYC oncogene, and these mutations are jointly associated with patient survival. These data demonstrate frequent mitochondrial mutation in prostate cancer, and suggest interplay between nuclear and mitochondrial mutational profiles in prostate cancer.In prostate cancer, the role of mutations in the maternally-inherited mitochondrial genome are not well known. Here, the authors demonstrate frequent, age-dependent mitochondrial mutation in prostate cancer. Strong links between mitochondrial and nuclear mutational profiles are associated with clinical aggressivity.