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

Kidney cancer: Links between hereditary syndromes and sporadic tumorigenesis

Multiple hereditary syndromes predispose to kidney cancer, including Von Hippel-Lindau syndrome, BAP1-Tumor Predisposition Syndrome, Hereditary Papillary Renal Cell Carcinoma, Tuberous Sclerosis Complex, Birt-Hogg-Dubé syndrome, Hereditary Paraganglioma-Pheochromocytoma Syndrome, Fumarate Hydratase Tumor Predisposition Syndrome, and Cowden syndrome. In some cases, mutations in the genes that cause hereditary kidney cancer are tightly linked to similar histologic features in sporadic RCC. For example, clear cell RCC occurs in the hereditary syndrome VHL, and sporadic ccRCC usually has inactivation of the VHL gene. In contrast, mutations in FLCN, the causative gene for Birt-Hogg-Dube syndrome, are rarely found in sporadic RCC. Here, we focus on the genes and pathways that link hereditary and sporadic RCC.

von-Hippel Lindau and Hypoxia-Inducible Factor at the Center of Renal Cell Carcinoma Biology

The most common form of kidney cancer is clear cell renal cell carcinoma (ccRCC). Biallelic VHL tumor suppressor gene inactivation is the usual initiating event in both hereditary (VHL Disease) and sporadic ccRCCs. The VHL protein, pVHL, earmarks the alpha subunits of the HIF transcription factor for destruction in an oxygen-dependent manner. Deregulation of HIF2 drives ccRCC pathogenesis. Drugs inhibiting the HIF2-responsive growth factor VEGF are now mainstays of ccRCC treatment. A first-in-class allosteric HIF2 inhibitor was recently approved for treating VHL Disease-associated neoplasms and appears active against sporadic ccRCC in early clinical trials.

Inter- and intra-tumor heterogeneity of genetic and immune profiles in inherited renal cell carcinoma

Patients with von Hippel-Lindau disease (vHL) are at risk of developing spatially and temporally multiple clear cell renal cell carcinomas (ccRCCs), which offers a valuable opportunity to analyze inter- and intra-tumor heterogeneity of genetic and immune profiles within the same patient. Here, we perform whole-exome and RNA sequencing, digital gene expression, and immunohistochemical analyses for 81 samples from 51 ccRCCs of 10 patients with vHL. Inherited ccRCCs are clonally independent and have less genomic alterations than sporadic ccRCCs. Hierarchical clustering of transcriptome profiles shows two clusters with distinct immune signatures: immune hot and cold clusters. Interestingly, not only samples from the same tumors but also different tumors from the same patients tend to show a similar immune signature, whereas samples from different patients frequently exhibit different signatures. Our findings reveal the genetic and immune landscape of inherited ccRCCs, demonstrating the relevance of host factors in shaping anti-tumor immunity.

Unraveling the Drivers of Tumorigenesis in the Context of Evolution: Theoretical Models and Bioinformatics Tools

Cancer originates from somatic cells that have accumulated mutations. These mutations alter the phenotype of the cells, allowing them to escape homeostatic regulation that maintains normal cell numbers. The emergence of malignancies is an evolutionary process in which the random accumulation of somatic mutations and sequential selection of dominant clones cause cancer cells to proliferate. The development of technologies such as high-throughput sequencing has provided a powerful means to measure subclonal evolutionary dynamics across space and time. Here, we review the patterns that may be observed in cancer evolution and the methods available for quantifying the evolutionary dynamics of cancer. An improved understanding of the evolutionary trajectories of cancer will enable us to explore the molecular mechanism of tumorigenesis and to design tailored treatment strategies.

Cell Adaptive Fitness and Cancer Evolutionary Dynamics

Genome instability of cancer cells translates into increased entropy and lower information processing capacity, leading to metabolic reprograming toward higher energy states, presumed to be aligned with a cancer growth imperative. Dubbed as the cell adaptive fitness, the proposition postulates that the coupling between cell signaling and metabolism constrains cancer evolutionary dynamics along trajectories privileged by the maintenance of metabolic sufficiency for survival. In particular, the conjecture postulates that clonal expansion becomes restricted when genetic alterations induce a sufficiently high level of disorder, that is, high entropy, in the regulatory signaling network, abrogating as a result the ability of cancer cells to successfully replicate, leading to a stage of clonal stagnation. The proposition is analyzed in the context of an in-silico model of tumor evolutionary dynamics to illustrate how cell-inherent adaptive fitness may predictably constrain clonal evolution of tumors, which would have significant implications for the design of adaptive cancer therapies.

Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues.

The Clinical and Molecular Features in the VHL Renal Cancers; Close or Distant Relatives with Sporadic Clear Cell Renal Cell Carcinoma?

Von Hippel-Lindau (VHL) disease is an autosomal dominant inherited cancer syndrome caused by germline mutations in the VHL tumor suppressor gene, characterized by the susceptibility to a wide array of benign and malign neoplasms, including clear-cell renal cell carcinoma. Moreover, VHL somatic inactivation is a crucial molecular event also in sporadic ccRCCs tumorigenesis. While systemic biomarkers in the VHL syndrome do not currently play a role in clinical practice, a new promising class of predictive biomarkers, microRNAs, has been increasingly studied. Lots of pan-genomic studies have deeply investigated the possible biological role of microRNAs in the development and progression of sporadic ccRCC; however, few studies have investigated the miRNA profile in VHL patients. Our review summarize all the new insights related to clinical and molecular features in VHL renal cancers, with a particular focus on the overlap with sporadic ccRCC.

Involvement of PBRM1 in VHL disease-associated clear cell renal cell carcinoma and its putative relationship with the HIF pathway

Von Hippel-Lindau (VHL) disease is the main cause of inherited clear-cell renal cell carcinoma (ccRCC) and is caused by germline mutations in the VHL tumor suppressor gene. Bi-allelic VHL alterations lead to inactivation of pVHL, which plays a major role by downstream activation of the hypoxia inducible factor (HIF) pathway. Somatic VHL mutations occur in 80% of sporadic ccRCC cases and the second most frequently mutated gene is polybromo 1 (PBRM1). As there is currently no data regarding PBRM1 involvement in VHL disease-associated ccRCC, the aim of the present study was to assess the PBRM1 mutational status, and PBRM1 and HIF expression in VHL disease-associated ccRCC series compared with a sporadic series. PBRM1 gene was screened by Sanger sequencing for 23 VHL-disease-associated ccRCC and 22 sporadic ccRCC cases. Immunohistochemical studies were performed to detect the expression of PBRM1, HIF1 and HIF2 for all cases. In VHL-associated tumors, 13.0% (n=3/23) had PBRM1 somatic mutations and 17.4% (n=4/23) had a loss of PBRM1 nuclear expression. In sporadic cases, 27.3% (n=6/22) showed PBRM1 somatic mutations and 45.5% (n=10/22) had a loss of PBRM1 nuclear expression. Loss of PBRM1 was associated with an advanced tumor stage. HIF1-positive tumors were observed more frequently in the VHL-associated ccRCC than in the sporadic series. Furthermore, in the VHL cohort, PBRM1 expression appeared to be associated more with HIF1 than with HIF2. Given that hereditary tumors tend to be less aggressive, these results would suggest that co-expression of PBRM1 and HIF1 may have a less oncogenic role in VHL-associated ccRCC.

Clonal Evolution of Multiple Myeloma-Clinical and Diagnostic Implications

Plasma cell dyscrasias are a heterogeneous group of diseases characterized by the expansion of bone marrow plasma cells. Malignant transformation of plasma cells depends on the continuity of events resulting in a sequence of well-defined disease stages, from monoclonal gammopathy of undetermined significance (MGUS) through smoldering myeloma (SMM) to symptomatic multiple myeloma (MM). Evolution of a pre-malignant cell into a malignant cell, as well as further tumor progression, dissemination, and relapse, require development of multiple driver lesions conferring selective advantage of the dominant clone and allowing subsequent evolution under selective pressure of microenvironment and treatment. This process of natural selection facilitates tumor plasticity leading to the formation of genetically complex and heterogenous tumors that are notoriously difficult to treat. Better understanding of the mechanisms underlying tumor evolution in MM and identification of lesions driving the evolution from the premalignant clone is therefore a key to development of effective treatment and long-term disease control. Here, we review recent advances in clonal evolution patterns and genomic landscape dynamics of MM, focusing on their clinical implications.

The Therapeutic Landscape of Renal Cell Carcinoma: From the Dark Age to the Golden Age

Oncologic treatments for renal cell carcinoma (RCC) have undergone a major revolution in the past 2 decades, moving away from the pre-2004 Dark Age during which interleukin 2 and interferon-α were the only therapeutic options and induced treatment responses in only 5% to 10% of patients with metastatic disease. The development of anti-angiogenic tyrosine kinase inhibitors against vascular endothelial growth factor receptor 2 and inhibitors of mammalian target of rapamycin complex 1 in 2005 introduced the Modern Age with better overall and progression-free survival and a greater number of patients (30%-40%) responding to and (∼80%) benefiting from these targeted therapeutic agents. The coming of age of the immuno-oncology era with the use of immune checkpoint inhibitors (ICIs) have ushered us into the Golden Age of metastatic RCC care, in which combined administrations of two ICIs (anti-programmed cell death protein 1/programmed death-ligand 1 and anti-cytotoxic T-lymphocyte-associated protein 4 or one tyrosine kinase inhibitor plus one ICI (anti-programmed cell death protein 1/programmed death-ligand 1) have recast the treatment landscape of clear cell RCC, the most common RCC subtype, with an approximately 60% response rate and an approximately 90% disease control rate that further improves metastatic RCC survival. Exciting clinical trials are in the pipeline investigating complementary/synergistic molecular mechanisms, based on studies investigating the biology, pathology, and genomics of renal carcinoma and the respective treatment outcome. This will enable us to enter the Diamond Age of precision medicine in which a specific treatment can be tailored to the specific biological and pathologic circumstance of an individual kidney tumor to offer more effective yet less toxic therapy.

Molecular and Metabolic Subtypes in Sporadic and Inherited Clear Cell Renal Cell Carcinoma

The promise of personalized medicine is a therapeutic advance where tumor signatures obtained from different omics platforms, such as genomics, transcriptomics, proteomics, and metabolomics, in addition to environmental factors including metals and metalloids, are used to guide the treatments. Clear cell renal carcinoma (ccRCC), the most common type of kidney cancer, can be sporadic (frequently) or genetic (rare), both characterized by loss of the von Hippel-Lindau (VHL) gene that controls hypoxia inducible factors. Recently, several genomic subtypes were identified with different prognoses. Transcriptomics, proteomics, metabolomics and metallomic data converge on altered metabolism as the principal feature of the disease. However, in view of multiple biochemical alterations and high level of tumor heterogeneity, identification of clearly defined subtypes is necessary for further improvement of treatments. In the future, single-cell combined multi-omics approaches will be the next generation of analyses gaining deeper insights into ccRCC progression and allowing for design of specific signatures, with better prognostic/predictive clinical applications.

HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures

Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles.

Genetic Alterations in Renal Cancers: Identification of The Mechanisms Underlying Cancer Initiation and Progression and of Therapeutic Targets

Renal cell cancer (RCC) involves three most recurrent sporadic types: clear-cell RCC (70-75%, CCRCC), papillary RCCC (10-15%, PRCC), and chromophobe RCC (5%, CHRCC). Hereditary cases account for about 5% of all cases of RCC and are caused by germline pathogenic variants. Herein, we review how a better understanding of the molecular biology of RCCs has driven the inception of new diagnostic and therapeutic approaches. Genomic research has identified relevant genetic alterations associated with each RCC subtype. Molecular studies have clearly shown that CCRCC is universally initiated by Von Hippel Lindau (VHL) gene dysregulation, followed by different types of additional genetic events involving epigenetic regulatory genes, dictating disease progression, aggressiveness, and differential response to treatments. The understanding of the molecular mechanisms that underlie the development and progression of RCC has considerably expanded treatment options; genomic data might guide treatment options by enabling patients to be matched with therapeutics that specifically target the genetic alterations present in their tumors. These new targeted treatments have led to a moderate improvement of the survival of metastatic RCC patients. Ongoing studies based on the combination of immunotherapeutic agents (immune check inhibitors) with VEGF inhibitors are expected to further improve the survival of these patients.

Exploiting the circuit breaker cancer evolution model in human clear cell renal cell carcinoma

The incessant interactions between susceptible humans and their respective macro/microenvironments registered throughout their lifetime result in the ultimate manifestation of individual cancers. With the average lifespan exceeding 50 years of age in humans since the beginning of 20th century, aging - the "time" factor - has played an ever-increasing role alongside host and environmental factors in cancer incidences. Cancer is a genetic/epigenetic disease due to gain-of-function mutations in cancer-causing genes (oncogene; OG) and/or loss-of-function mutations in tumor-suppressing genes (tumor suppressor genes; TSG). In addition to their integral relationship with cancer, a timely deployment of specific OG and/or TSG is in fact needed for higher organisms like human to cope with respective physiological and pathological conditions. Over the past decade, extensive human kidney cancer genomics have been performed and novel mouse models recapitulating human kidney cancer pathobiology have been generated. With new genomic, genetic, mechanistic, clinical and therapeutic insights accumulated from studying clear cell renal cell carcinoma (ccRCC)-the most common type of kidney cancer, we conceived a cancer evolution model built upon the OG-TSG signaling pair analogous to the electrical circuit breaker (CB) that permits necessary signaling output and at the same time prevent detrimental signaling overdrive. Hence, this viewpoint aims at providing a step-by-step mechanistic explanation/illustration concerning how inherent OG-TSG CBs intricately operate in concert for the organism's wellbeing; and how somatic mutations, the essential component for genetic adaptability, inadvertently triggers a sequential outage of specific sets of CBs that normally function to maintain and protect and individual tissue homeostasis.

Significance of PI3K signalling pathway in clear cell renal cell carcinoma in relation to VHL and HIF status

Renal cell carcinoma (RCC) includes diverse tumour types characterised by various genetic abnormalities. The genetic changes, like mutations, deletions and epigenetic alterations, play a crucial role in the modification of signalling networks, tumour pathogenesis and prognosis. The most prevalent RCC type, clear cell RCC (ccRCC), is asymptomatic in the early stages and has a poorer prognosis compared with the papillary and the chromophobe types RCCs. Generally, ccRCC is refractory to chemotherapy and radiation therapy. Loss of von Hippel-Lindau (VHL) gene and upregulation of hypoxia-inducible factors (HIF), the signature of most sporadic ccRCC, promote multiple growth factors. Hence, VHL/HIF and a variety of pathways, including phosphatase and TEnsin homolog on chromosome 10/phosphatidylinositol-3-kinase (PI3K)/AKT, are closely connected and contribute to the ontogeny of ccRCC. In the recent decade, multiple targeting agents have been developed based on blocking major signalling pathways directly or indirectly involved in ccRCC tumour progression, metastasis, angiogenesis and survival. However, most of these drugs have limitations; either metastatic ccRCC develops resistance to these agents, or despite blocking receptors, tumour cells use alternate signalling pathways. This review compiles the state of knowledge about the PI3K/AKT signalling pathway confined to ccRCC and its cross-talks with VHL/HIF pathway.

[Hereditary kidney cancers: The pathologist's view in 2020]

Hereditary predispositions to adult kidney tumors involve around 5% of tumors and include a dozen of autosomal dominant syndromes. The most frequent tumors encountered in these setting are clear cell renal cell carcinomas, papillary renal cell carcinomas, chromophobe renal cell carcinomas and angiomyolipomas. Their detection is essential in order to adapt individual care and perform genetic screening of at-risk relatives, especially in the national french network PREDIR, labeled by the National Cancer Institute and dedicated to hereditary predispositions to kidney tumors. Targeted genetic analysis, which was guided in particular by the renal tumor subtype, has recently evolved into genetic analysis using panels of genes. Pathologist contribution's remains however central in the diagnosis of hereditary forms since we currently have immunohistochemical biomarkers that allow us to diagnose two specifically hereditary entities: hereditary leiomyomatosis and renal cell carcinoma associated-renal cell carcinoma, associated with a loss of fumarate hydratase and succinate dehydrogenase-deficient renal cell carcinoma associated with a loss of succinate deshydrogenase B expression. These diagnoses must however be confirmed by the identification of pathogenic germline variation in the corresponding genes. Improvement of kidney tumors characterization has also lead to identify new subtypes, expanding the algorithm of renal tumors associated with hereditary setting. Here we aim to review all subtypes of adult renal tumors encountered in predisposition syndromes.

A Pilot Study of Baseline Spatial Genomic Heterogeneity in Primary Gastric Cancers Using Multi-Region Endoscopic Sampling

Intertumoral heterogeneity among actionable biomarkers including ERBB2, FGFR2 and EGFR has been observed to occur under therapeutic pressure in advanced gastric cancer. However, baseline intratumoral heterogeneity at diagnosis is understudied and may impact clinical outcomes. We sought to explore intratumoral heterogeneity in primary advanced gastric cancers via DNA sequencing from multi-region endoscopic sampling at diagnosis. Patients with newly diagnosed advanced gastric adenocarcinoma underwent endoscopic mapping and pre-determined 8-sector biopsy of the primary tumor with concurrent plasma cfDNA sampling. Biopsy samples were subjected to targeted next generation sequencing and plasma cfDNA was analyzed via a 28-gene cfDNA assay. Expectedly, we observed that the majority of genetic alterations were shared among multi-sector biopsies within the same gastric primary tumor. However, all samples contained private subclonal alterations between biopsy sectors, including actionable alterations in GNAS and STK11. Cell free DNA analyses also exhibited both shared and non-shared alterations between mutations detected in cfDNA and tumor tissue biopsies confirming baseline intertumoral heterogeneity. This is the first dataset to confirm baseline intratumoral heterogeneity and confirms that multi-sector endoscopic biopsy is feasible and capable of capturing intratumoral heterogeneity among relevant genomic alterations in gastric cancer. Both multi-sector endoscopic biopsies and cfDNA analyses are complementary in capturing the diverse mutational landscape at disease presentation.

The circular RNA hsa-circ-0072309 plays anti-tumour roles by sponging miR-100 through the deactivation of PI3K/AKT and mTOR pathways in the renal carcinoma cell lines

Aims: To explore the roles and regulatory mechanisms of the circular RNA (circRNA)-hsa-circ-0072309 in CAKI-1 and ACHN cells. Methods: CAKI-1 and ACHN cells were transfected with hsa-circ-0072309 overproduction vector (circRNA) and microRNA-100 (miR-100) mimic or the corresponding controls. Cell viability was detected with the CCK-8. The protein expression levels of p53, c-Myc, cleaved-caspase-3/9, matrix metalloproteinase (MMP)-2/9, vimentin, AKT, PI3K and mTOR were individually determined through western blot. qRT-PCR was used to examine the expressions of hsa-circ-0072309 and miR-100. The apoptotic rate and the migration or invasion rates were separately determined by the annexin v-FITC/PI with a flow cytometer and modified two-chamber migration assay or millicell hanging cell culture. Results: The hsa-circ-0072309 was poorly expressed in tumor tissue. Abundant hsa-circ-0072309 induced the inhibitions of cell proliferation, migration and invasion, as well as the PI3K/AKT and the mTOR cascades but enhancement of apoptosis. circRNA stimulated the down-regulation of miR-100, which was low-expressed in tumour tissue and whose overproduction abolished the impacts of circRNA on these elements mentioned above. Conclusion: The hsa-circ-0072309 played anti-tumour roles by targeting miR-100 by blocking the PI3K/AKT and mTOR cascades in the CAKI-1 and ACHN cell lines.

The incidence of consecutive manifestations in Von Hippel-Lindau disease

Von Hippel-Lindau (VHL) disease is an autosomal dominant rare tumor syndrome characterized by high penetrance. VHL mutation carriers develop numerous manifestations in multiple organs during life. The natural course of development of new and growth of existing VHL-related manifestations is still unclear. In this study we aimed to gain insight into the development of subsequent manifestations in VHL disease. We retrospectively scored each new VHL-related manifestation as detected by standard follow-up (retina, central nervous system, kidneys and pancreas, excluding adrenal and endolymfatic sac manifestations) in 75 VHL mutation carriers. The Kaplan–Meier method was used to plot the cumulative proportions of all consecutive manifestations in each organ against age. The cumulative average number of manifestations in all organs during life was calculated by summating these cumulative proportions. Poisson model parameters were used to calculate average time to the detection of consecutive VHL manifestations in each organ. Consecutive VHL-related kidney and retina manifestations during life occur linearly according to Poisson distribution model. The total number of VHL manifestations rises linearly, with an average of seven VHL-related lesions at age 60 years. The incidence of consecutive VHL-related manifestations is constant during life in VHL mutation carriers. Our data is consistent with the notion that somatic inactivation of the remaining allele (Knudson’s “two-hit” hypothesis) is the determining factor in developing new VHL-related manifestations.

Perspective: Potential Impact and Therapeutic Implications of Oncogenic PI3K Activation on Chromosomal Instability

Genetic activation of the class I PI3K pathway is very common in cancer. This mostly results from oncogenic mutations in PIK3CA, the gene encoding the ubiquitously expressed PI3Kα catalytic subunit, or from inactivation of the PTEN tumour suppressor, a lipid phosphatase that opposes class I PI3K signalling. The clinical impact of PI3K inhibitors in solid tumours, aimed at dampening cancer-cell-intrinsic PI3K activity, has thus far been limited. Challenges include poor drug tolerance, incomplete pathway inhibition and pre-existing or inhibitor-induced resistance. The principle of pharmacologically targeting cancer-cell-intrinsic PI3K activity also assumes that all cancer-promoting effects of PI3K activation are reversible, which might not be the case. Emerging evidence suggests that genetic PI3K pathway activation can induce and/or allow cells to tolerate chromosomal instability, which-even if occurring in a low fraction of the cell population-might help to facilitate and/or drive tumour evolution. While it is clear that such genomic events cannot be reverted pharmacologically, a role for PI3K in the regulation of chromosomal instability could be exploited by using PI3K pathway inhibitors to prevent those genomic events from happening and/or reduce the pace at which they are occurring, thereby dampening cancer development or progression. Such an impact might be most effective in tumours with clonal PI3K activation and achievable at lower drug doses than the maximum-tolerated doses of PI3K inhibitors currently used in the clinic.

SEECancer: a resource for somatic events in evolution of cancer genome

Cancer cells progressively evolve from a premalignant to a malignant state, which is driven by accumulating somatic alterations that confer normal cells a fitness advantage. Improvements in high-throughput sequencing techniques have led to an increase in construction of tumor phylogenetics and identification of somatic driver events that specifically occurred in different tumor progression stages. Here, we developed the SEECancer database (http://biocc.hrbmu.edu.cn/SEECancer), which aims to present the comprehensive cancer evolutionary stage-specific somatic events (including early-specific, late-specific, relapse-specific, metastasis-specific, drug-resistant and drug-induced genomic events) and their temporal orders. By manually curating over 10 000 published articles, 1231 evolutionary stage-specific genomic events and 5772 temporal orders involving 82 human cancers and 23 tissue origins were collected and deposited in the SEECancer database. Each entry contains the somatic event, evolutionary stage, cancer type, detection approach and relevant evidence. SEECancer provides a user-friendly interface for browsing, searching and downloading evolutionary stage-specific somatic events and temporal relationships in various cancers. With increasing attention on cancer genome evolution, the necessary information in SEECancer will facilitate understanding of cancer etiology and development of evolutionary therapeutics, and help clinicians to discover biomarkers for monitoring tumor progression.

Estimating the predictability of cancer evolution

MOTIVATION

How predictable is the evolution of cancer? This fundamental question is of immense relevance for the diagnosis, prognosis and treatment of cancer. Evolutionary biologists have approached the question of predictability based on the underlying fitness landscape. However, empirical fitness landscapes of tumor cells are impossible to determine in vivo. Thus, in order to quantify the predictability of cancer evolution, alternative approaches are required that circumvent the need for fitness landscapes.

RESULTS

We developed a computational method based on conjunctive Bayesian networks (CBNs) to quantify the predictability of cancer evolution directly from mutational data, without the need for measuring or estimating fitness. Using simulated data derived from >200 different fitness landscapes, we show that our CBN-based notion of evolutionary predictability strongly correlates with the classical notion of predictability based on fitness landscapes under the strong selection weak mutation assumption. The statistical framework enables robust and scalable quantification of evolutionary predictability. We applied our approach to driver mutation data from the TCGA and the MSK-IMPACT clinical cohorts to systematically compare the predictability of 15 different cancer types. We found that cancer evolution is remarkably predictable as only a small fraction of evolutionary trajectories are feasible during cancer progression.

AVAILABILITY AND IMPLEMENTATION

https://github.com/cbg-ethz/predictability\_of\_cancer\_evolution.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Revealing clonality and subclonality of driver genes for clinical survival benefits in breast cancer

PURPOSE

Genomic studies have revealed that genomic aberrations play important roles in the progression of this disease. The aim of this study was to evaluate the associations between clinical survival outcomes of the clonality and subclonality status of driver genes in breast cancer.

METHODS

We performed an integrated analysis to infer the clonal status of 55 driver genes in breast cancer data from TCGA. We used the chi-squared test to assess the relations between clonality of driver gene mutations and clinicopathological factors. The Kaplan-Meier method was performed for the visualization and the differences between survival curves were calculated by log-rank test. Univariate and multivariate Cox proportional hazards regression models were used to adjust for clinicopathological factors.

RESULTS

We identified a high proportion of clonal mutations in these driver genes. Among them, there were 17 genes showing significant associations between their clonality and multiple clinicopathologic factors. Performing survival analysis on BRCA patients with clonal or subclonal driver gene mutations, we found that clonal ERBB2, FOXA1, and KMT2C mutations and subclonal GATA3 and RB1 mutations predicted shorter overall survival compared with those with wild type. Furthermore, clonal ERBB2 and FOXA1 mutations and subclonal GATA3 and RB1 mutations independently predicted for shorter overall survival after adjusting for clinicopathological factors. By longitudinal analysis, the clonality of ERBB2, FOXA1, GATA3, and RB1 significantly predicted patients' outcome within some specific BRCA tumor stages and histological subtypes.

CONCLUSIONS

In summary, these clonal or subclonal mutations of driver genes have implications for diagnosis, prognosis, and treatment with BRCA patients.

Intratumor genetic heterogeneity in squamous cell carcinoma of the oral cavity

BACKGROUND

We sought to evaluate intratumor heterogeneity in squamous cell carcinoma of the oral cavity (OCC) and specifically determine the effect of physical separation and histologic differentiation within the same tumor.

METHODS

We performed whole exome sequencing on five biopsy sites-two from well-differentiated, two from poorly differentiated regions, and one from normal parenchyma-from five primary OCC specimens.

RESULTS

We found high levels of intratumor heterogeneity and, in four primary tumors, identified only 0 to 2 identical mutations in all subsites. We found that the heterogeneity inversely correlated with physical separation and that pairs of well-differentiated samples were more similar to each other than analogous poorly differentiated specimens. Only TP53 mutations, but not other purported "driver mutations" in head and neck squamous cell carcinoma, were found in multiple biopsy sites.

CONCLUSION

These data highlight the challenges to characterization of the mutational landscape of OCC with single site biopsy and have implications for personalized medicine.

Reprogramming normal human epithelial tissues to a common, lethal neuroendocrine cancer lineage

The use of potent therapies inhibiting critical oncogenic pathways active in epithelial cancers has led to multiple resistance mechanisms, including the development of highly aggressive, small cell neuroendocrine carcinoma (SCNC). SCNC patients have a dismal prognosis due in part to a limited understanding of the molecular mechanisms driving this malignancy and the lack of effective treatments. Here, we demonstrate that a common set of defined oncogenic drivers reproducibly reprograms normal human prostate and lung epithelial cells to small cell prostate cancer (SCPC) and small cell lung cancer (SCLC), respectively. We identify shared active transcription factor binding regions in the reprogrammed prostate and lung SCNCs by integrative analyses of epigenetic and transcriptional landscapes. These results suggest that neuroendocrine cancers arising from distinct epithelial tissues may share common vulnerabilities that could be exploited for the development of drugs targeting SCNCs.

Chromosome 3p Loss-Orchestrated VHL, HIF, and Epigenetic Deregulation in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common renal cell carcinoma subtype, and metastatic ccRCC is associated with 5-year survival rates of 10% to 20%. Genetically, ccRCC originates from sequential losses of multiple tumor suppressor genes. Remarkably, chromosome 3p loss occurs in more than 90% of sporadic ccRCCs. This results in concurrent one-copy loss of four tumor suppressor genes that are also mutated individually at high frequency in ccRCC (ie, VHL, 80%; PBRM1, 29% to 46%; BAP1, 6% to 19%; and SETD2, 8% to 30%). Pathogenically, 3p loss probably represents the first genetic event that occurs in sporadic ccRCC and the second genetic event in VHL-mutated hereditary ccRCC. VHL constitutes the substrate recognition module of the VCB-Cul2 E3 ligase that degrades HIF1/2α, whereas PBRM1, BAP1, and SETD2 are epigenetic modulators that regulate gene transcription. Because 3p loss and VHL inactivation are nearly universal truncal events in ccRCC, the resulting HIF1/2 signaling overdrive and accompanied tumor hypervascularization probably underlie the therapeutic benefits observed with vascular endothelial growth factor receptor inhibitors, including sorafenib, sunitinib, pazopanib, axitinib, bevacizumab, cabozantinib, and lenvatinib. Furthermore, recent marked advances in ccRCC genomics, transcriptomics, proteomics, metabolomics, molecular mechanisms, mouse models, prognostic and predictive biomarkers, and clinical trials have rendered invaluable translational insights concerning precision kidney cancer therapeutics. With an armamentarium encompassing 13 drugs that exploit seven unique therapeutic mechanisms (ie, cytokines, vascular endothelial growth factor receptor, mTORC1, cMET/AXL, fibroblast growth factor receptor, programmed cell death-1 and programmed death-ligand 1, and cytotoxic T-cell lymphocyte associated-4) to treat metastatic renal cell carcinoma, one of the imminent clinical questions concerning care of patients with metastatic ccRCC is how a personalized treatment strategy, through rationally combining and sequencing different therapeutic modalities, can be formulated to offer the best clinical outcome for individual patients. Here, we attempt to integrate recent discoveries of immediate translational impacts and discuss future translational challenges and opportunities.

Structural Alterations in Human Fibroblast Growth Factor Receptors in Carcinogenesis

Fibroblast growth factor (FGF) plays an important role in human embryogenesis, angiogenesis, cell proliferation, and differentiation. Carcinogenesis is accompanied by aberrant constitutive activation of FGF receptors (FGFRs) resulting from missense mutation in the FGFR1-4 genes, generation of chimeric oncogenes, FGFR1-4 gene amplification, alternative splicing shift toward formation of mesenchymal FGFR isoforms, and FGFR overexpression. Altogether, these alterations contribute to auto- and paracrine stimulation of cancer cells and neoangiogenesis. Certain missense mutations are found at a high rate in urinary bladder cancer and can be used for non-invasive cancer recurrence diagnostics by analyzing urine cell pellet DNA. Chimeric FGFR1/3 and amplified FGFR1/2 genes can predict cell response to the targeted therapy in various oncological diseases. In recent years, high-throughput sequencing has been used to analyze exomes of virtually all human tumors, which allowed to construct phylogenetic trees of clonal cancer evolution with special emphasis on driver mutations in FGFR1-4 genes. At present, FGFR blockers, such as multi-kinase inhibitors, specific FGFR inhibitors, and FGF ligand traps are being tested in clinical trials. In this review, we discuss current data on the functioning of the FGFR family proteins in both normal and cancer cells, mutations in the FGFR1-4 genes, and mechanisms underlying their oncogenic potential, which might be interesting to a broad range of scientists searching for specific tumor markers and targeted anti-cancer drugs.

Integrative analysis of dysregulated microRNAs and mRNAs in multiple recurrent synchronized renal tumors from patients with von Hippel-Lindau disease

Von Hippel-Lindau (VHL) disease is a rare auto-somal dominant syndrome that is the main cause of inherited clear-cell renal cell carcinoma (ccRCC), which generally occurs in the form of multiple recurrent synchronized tumors. Affected patients are carriers of a germline mutation in the VHL tumor suppressor gene. Somatic mutations of this gene are also found in sporadic ccRCC and numerous pan-genomic studies have reported a dysregulation of microRNA (miRNA) expression in these sporadic tumors. In order to investigate the molecular mechanisms underlying the pathogenesis of VHL-associated ccRCC, particularly in the context of multiple tumors, the present study characterized the mRNA and miRNA transcriptome through an integrative analysis compared with sporadic renal tumors. In the present study, two series of ccRCC samples were used. The first set consisted of several samples from different tumors occurring in the same patient, for two independent patients affected with VHL disease. The second set consisted of 12 VHL-associated tumors and 22 sporadic ccRCC tumors compared with a pool of normal renal tissue. For each sample series, an expression analysis of miRNAs and mRNAs was conducted using microarrays. The results indicated that multiple tumors within the kidney of a patient with VHL disease featured a similar pattern of miRNA and gene expression. In addition, the expression levels of miRNA were able to distinguish VHL-associated tumors from sporadic ccRCC, and it was identified that 103 miRNAs and 2,474 genes were differentially expressed in the ccRCC series compared with in normal renal tissue. The majority of dysregulated genes were implicated in 'immunity' and 'metabolism' pathways. Taken together, these results allow a better understanding of the occurrence of ccRCC in patients with VHL disease, by providing insights into dysregulated miRNA and mRNA. In the set of patients with VHL disease, there were few differences in miRNA and mRNA expression, thus indicating a similar molecular evolution of these synchronous tumors and suggesting that the same molecular mechanisms underlie the pathogenesis of these hereditary tumors.

Genomic classifications of renal cell carcinoma: a critical step towards the future application of personalized kidney cancer care with pan-omics precision

Over the past 20 years, classifications of kidney cancer have undergone major revisions based on morphological refinements and molecular characterizations. The 2016 WHO classification of renal tumors recognizes more than ten different renal cell carcinoma (RCC) subtypes. Furthermore, the marked inter- and intra-tumor heterogeneity of RCC is now well appreciated. Nevertheless, contemporary multi-omics studies of RCC, encompassing genomics, transcriptomics, proteomics, and metabolomics, not only highlight apparent diversity but also showcase and underline commonality. Here, we wish to provide an integrated perspective concerning the future 'functional' classification of renal cancer by bridging gaps among morphology, biology, multi-omics, and therapeutics. This review focuses on recent progress and elaborates the potential value of contemporary pan-omics approaches with a special emphasis on cancer genomics unveiled through next-generation sequencing technology, and how an integrated multi-omics approach might impact precision-based personalized kidney cancer care in the near future. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The von Hippel–Lindau Tumor Suppressor Protein

The von Hippel–Lindau (VHL) gene is a two-hit tumor suppressor gene and is linked to the development of the most common form of kidney cancer, clear cell renal carcinoma; blood vessel tumors of the retina, cerebellum, and spinal cord called hemangioblastomas; and tumors of the sympathoadrenal nervous system called paragangliomas. The VHL gene product, pVHL, is the substrate recognition subunit of a cullin-dependent ubiquitin ligase that targets the α subunits of hypoxia-inducible factor (HIF) for destruction when oxygen is plentiful. Mounting evidence implicates HIF2 in the pathogenesis of pVHL-defective tumors and has provided a conceptual foundation for the development of drugs to treat them that inhibit HIF2-responsive gene products such as VEGF and, more recently, HIF2 itself. pVHL has additional, noncanonical functions that are cancer relevant, including roles related to the primary cilium, chromosome stability, extracellular matrix formation, and survival signaling.

Dynamic epigenetic changes to VHL occur with sunitinib in metastatic clear cell renal cancer

Background

Genetic intratumoral heterogeneity (ITH) hinders biomarker development in metastatic clear cell renal cancer (mccRCC). Epigenetic relative to genetic ITH or the presence of consistent epigenetic changes following targeted therapy in mccRCC have not been evaluated. The aim of this study was to determine methylome/genetic ITH and to evaluate specific epigenetic and genetic changes associated with sunitinib therapy.

Patients and methods

Multi-region DNA sampling performed on sequential frozen pairs of primary tumor tissue from 14 metastatic ccRCC patients, in the Upfront Sunitinib (SU011248) Therapy Followed by Surgery in Patients with Metastatic Renal Cancer: a Pilot Phase II Study (SuMR; ClinicalTrials.gov identifier: NCT01024205), at presentation (biopsy) and after 3-cycles of 50mg sunitinib (nephrectomy). Untreated biopsy and nephrectomy samples before and after renal artery ligation were controls. Ion Proton sequencing of 48 key ccRCC genes, and MethylCap-seq DNA methylation analysis was performed, data was analysed using the statistical computing environment R.

Results

Unsupervised hierarchical clustering revealed complete methylome clustering of biopsy and three nephrectomy samples for each patient (14/14 patients). For mutational status, untreated biopsy and all treated nephrectomy samples clustered together in 8/13 (61.5%) patients. The only methylation target significantly altered following sunitinib therapy was VHL promoter region 7896829 which was hypermethylated with treatment (FDR=0.077, P<0.001) and consistent for all patients (pre-treatment 50% patients had VHL mutations, 14% patients VHL hypermethylation). Renal artery ligation did not affect this result. No significant differences in driver or private mutation count was found with sunitinib treatment.

Conclusions

Demonstration of relative methylome homogeneity and consistent VHL hypermethylation, after sunitinib, may overcome the hurdle of ITH present at other molecular levels for biomarker research.

Renal Cell Carcinoma in von Hippel-Lindau Disease-From Tumor Genetics to Novel Therapeutic Strategies

von Hippel-Lindau (VHL) disease is an autosomal dominant syndrome caused by mutations in the VHL tumor-suppressor gene, leading to the dysregulation of many hypoxia-induced genes. Affected individuals are at increased risk of developing recurrent and bilateral kidney cysts and dysplastic lesions which may progress to clear cell renal cell carcinoma (ccRCC). Following the eponymous VHL gene inactivation, ccRCCs evolve through additional genetic alterations, resulting in both intratumor and intertumor heterogeneity. Genomic studies have identified frequent mutations in genes involved in epigenetic regulation and phosphoinositide 3-kinase-AKT-mechanistic target of rapamycin (mTOR) pathway activation. Currently, local therapeutic options include nephron-sparing surgery and alternative ablative procedures. For advanced metastatic disease, systemic treatment, including inhibition of vascular endothelial growth factor pathways and mTOR pathways, as well as immunotherapy are available. Multimodal therapy, targeting multiple signaling pathways and/or enhancing the immune response, is currently being investigated. A deeper understanding of the fundamental biology of ccRCC development and progression, as well as the development of novel and targeted therapies will be accelerated by new preclinical models, which will greatly inform the search for clinical biomarkers for diagnosis, prognosis, and response to treatment.

The SWI/SNF Protein PBRM1 Restrains VHL-Loss-Driven Clear Cell Renal Cell Carcinoma

PBRM1 is the second most commonly mutated gene after VHL in clear cell renal cell carcinoma (ccRCC). However, the biological consequences of PBRM1 mutations for kidney tumorigenesis are unknown. Here, we find that kidney-specific deletion of Vhl and Pbrm1, but not either gene alone, results in bilateral, multifocal, transplantable clear cell kidney cancers. PBRM1 loss amplified the transcriptional outputs of HIF1 and STAT3 incurred by Vhl deficiency. Analysis of mouse and human ccRCC revealed convergence on mTOR activation, representing the third driver event after genetic inactivation of VHL and PBRM1. Our study reports a physiological preclinical ccRCC mouse model that recapitulates somatic mutations in human ccRCC and provides mechanistic and therapeutic insights into PBRM1 mutated subtypes of human ccRCC.

Simultaneous evolutionary expansion and constraint of genomic heterogeneity in multifocal lung cancer

Recent genomic analyses have revealed substantial tumor heterogeneity across various cancers. However, it remains unclear whether and how genomic heterogeneity is constrained during tumor evolution. Here, we sequence a unique cohort of multiple synchronous lung cancers (MSLCs) to determine the relative diversity and uniformity of genetic drivers upon identical germline and environmental background. We find that each multicentric primary tumor harbors distinct oncogenic alterations, including novel mutations that are experimentally demonstrated to be functional and therapeutically targetable. However, functional studies show a strikingly constrained tumorigenic pathway underlying heterogeneous genetic variants. These results suggest that although the mutation-specific routes that cells take during oncogenesis are stochastic, genetic trajectories may be constrained by selection for functional convergence on key signaling pathways. Our findings highlight the robust evolutionary pressures that simultaneously shape the expansion and constraint of genomic diversity, a principle that holds important implications for understanding tumor evolution and optimizing therapeutic strategies.Across cancer types tumor heterogeneity has been observed, but how this relates to tumor evolution is unclear. Here, the authors sequence multiple synchronous lung cancers, highlighting the evolutionary pressures that simultaneously shape the expansion and constraint of genomic heterogeneity.

Epigenome Aberrations: Emerging Driving Factors of the Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC), the most common form of Kidney cancer, is characterized by frequent mutations of the von Hippel-Lindau (VHL) tumor suppressor gene in ~85% of sporadic cases. Loss of pVHL function affects multiple cellular processes, among which the activation of hypoxia inducible factor (HIF) pathway is the best-known function. Constitutive activation of HIF signaling in turn activates hundreds of genes involved in numerous oncogenic pathways, which contribute to the development or progression of ccRCC. Although VHL mutations are considered as drivers of ccRCC, they are not sufficient to cause the disease. Recent genome-wide sequencing studies of ccRCC have revealed that mutations of genes coding for epigenome modifiers and chromatin remodelers, including PBRM1, SETD2 and BAP1, are the most common somatic genetic abnormalities after VHL mutations in these tumors. Moreover, recent research has shed light on the extent of abnormal epigenome alterations in ccRCC tumors, including aberrant DNA methylation patterns, abnormal histone modifications and deregulated expression of non-coding RNAs. In this review, we discuss the epigenetic modifiers that are commonly mutated in ccRCC, and our growing knowledge of the cellular processes that are impacted by them. Furthermore, we explore new avenues for developing therapeutic approaches based on our knowledge of epigenome aberrations of ccRCC.

Molecular Classification of Renal Cell Carcinoma and Its Implication in Future Clinical Practice

Renal cell carcinoma (RCC) encompasses a wide spectrum of morphologically and molecularly distinct (>10) cancer subtypes originated from the kidney epithelium. Metastatic RCC (mRCC) is lethal and refractory to conventional chemotherapeutic agents. The incorporation of targeted therapies and immune checkpoint inhibitors into the current practice of mRCC has markedly improved the median overall survival of clear cell RCC (ccRCC) patients, the most common subtype, but not rare kidney cancer (RKC or non-ccRCC, nccRCC). Varied treatment response in mRCC patients is observed, which presents clinical challenges/opportunities at the modern mRCC therapeutic landscape consisting of 12 approved drugs representing 6 different effective mechanisms. Key contributing factors include inter- and intra-RCC heterogeneity. With the advances in pan-omics technologies, we now have a better understanding of the molecular pathobiology of individual RCC subtype. Here, we attempt to classify ccRCC based on contemporary molecular features with emphasis on their respective potential significance in clinical practice.

HIF activation causes synthetic lethality between the VHL tumor suppressor and the EZH1 histone methyltransferase

Inactivation of the von Hippel-Lindau tumor suppressor protein (pVHL) is the signature lesion in the most common form of kidney cancer, clear cell renal cell carcinoma (ccRCC). pVHL loss causes the transcriptional activation of hypoxia-inducible factor (HIF) target genes, including many genes that encode histone lysine demethylases. Moreover, chromatin regulators are frequently mutated in this disease. We found that ccRCC displays increased H3K27 acetylation and a shift toward mono- or unmethylated H3K27 caused by an HIF-dependent increase in H3K27 demethylase activity. Using a focused short hairpin RNA library, as well as CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9) and a pharmacological inhibitor, we discovered that pVHL-defective ccRCC cells are hyperdependent on the H3K27 methyltransferase EZH1 for survival. Therefore, targeting EZH1 could be therapeutically useful in ccRCC.

Cousins not twins: intratumoural and intertumoural heterogeneity in syndromic neuroendocrine tumours

Hereditary endocrine neoplasias, including phaeochromocytoma/paraganglioma and medullary thyroid cancer, are caused by autosomal dominant mutations in several familial cancer genes. A common feature of these diseases is the presentation of multiple primary tumours, or multifocal disease representing independent tumour clones that have arisen from the same initiating genetic lesion, but have undergone independent clonal evolution. Such tumours provide an opportunity to discover common cooperative changes required for tumourigenesis, while controlling for the genetic background of the individual. We performed genomic analysis of synchronous and metachronous tumours from five patients bearing germline mutations in the genes SDHB, RET, and MAX. Using whole exome sequencing and high-density single-nucleotide polymorphism arrays, we analysed two to four primary tumours from each patient. We also applied multi-region sampling, to assess intratumoural heterogeneity and clonal evolution, in two cases involving paraganglioma and medullary thyroid cancer, respectively. Heterogeneous patterns of genomic change existed between synchronous or metachronous tumours, with evidence of branching evolution. We observed striking examples of evolutionary convergence involving the same rare somatic copy-number events in synchronous primary phaeochromocytoma/paraganglioma. Convergent events also occurred during clonal evolution of metastatic medullary thyroid cancer. These observations suggest that genetic or epigenetic changes acquired early within precursor cells, or pre-existing within the genetic background of the individual, create contingencies that determine the evolutionary trajectory of the tumour. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Overcome tumor heterogeneity-imposed therapeutic barriers through convergent genomic biomarker discovery: A braided cancer river model of kidney cancer

Tumor heterogeneity, encompassing genetic, epigenetic, and microenvironmental variables, is extremely complex and presents challenges to cancer diagnosis and therapy. Genomic efforts on genetic intratumor heterogeneity (G-ITH) confirm branched evolution, support the trunk-branch cancer model, and present a seemingly insurmountable obstacle to conquering cancers. G-ITH is conspicuous in clear cell renal cell carcinoma (ccRCC), where its presence complicates identification and validation of biomarkers and thwarts efforts in advancing precision cancer therapeutics. However, long-term clinical benefits on targeted therapy are not uncommon in metastatic ccRCC patients, implicating that there are underlying constraints during ccRCC evolution, which in turn force a nonrandom sequence of parallel gene/pathway/function/phenotype convergence within individual tumors. Accordingly, we proposed a "braided cancer river model" depicting ccRCC evolution, which deduces cancer development based on multiregion tumor genomics of exceptional mTOR inhibitor (mTORi) responders. Furthermore, we employ an outlier case to explore the river model and highlight the importance of "Five NGS Matters: Number, Frequency, Position, Site and Time" in assessing cancer genomics for precision medicine. This mutable cancer river model may capture clinically significant phenotype-convergent events, predict vulnerability/resistance mechanisms, and guide effective therapeutic strategies. Our model originates from studying exceptional responders in ccRCC, which warrants further refinement and future validation concerning its applicability to other cancer types. The goal of this review is employing kidney cancer as an example to illustrate critical issues concerning tumor heterogeneity.

Tissue-specific tumorigenesis: context matters

How can we treat cancer more effectively? Traditionally, tumours from the same anatomical site are treated as one tumour entity. This concept has been challenged by recent breakthroughs in cancer genomics and translational research that have enabled molecular tumour profiling. The identification and validation of cancer drivers that are shared between different tumour types, spurred the new paradigm to target driver pathways across anatomical sites by off-label drug use, or within so-called basket or umbrella trials which are designed to test whether molecular alterations in one tumour entity can be extrapolated to all others. However, recent clinical and preclinical studies suggest that there are tissue- and cell type-specific differences in tumorigenesis and the organization of oncogenic signalling pathways. In this Opinion article, we focus on the molecular, cellular, systemic and environmental determinants of organ-specific tumorigenesis and the mechanisms of context-specific oncogenic signalling outputs. Investigation, recognition and in-depth biological understanding of these differences will be vital for the design of next-generation clinical trials and the implementation of molecularly guided cancer therapies in the future.

A population genetics perspective on the determinants of intra-tumor heterogeneity

Cancer results from the acquisition of somatic alterations in a microevolutionary process that typically occurs over many years, much of which is occult. Understanding the evolutionary dynamics that are operative at different stages of progression in individual tumors might inform the earlier detection, diagnosis, and treatment of cancer. Although these processes cannot be directly observed, the resultant spatiotemporal patterns of genetic variation amongst tumor cells encode their evolutionary histories. Such intra-tumor heterogeneity is pervasive not only at the genomic level, but also at the transcriptomic, phenotypic, and cellular levels. Given the implications for precision medicine, the accurate quantification of heterogeneity within and between tumors has become a major focus of current research. In this review, we provide a population genetics perspective on the determinants of intra-tumor heterogeneity and approaches to quantify genetic diversity. We summarize evidence for different modes of evolution based on recent cancer genome sequencing studies and discuss emerging evolutionary strategies to therapeutically exploit tumor heterogeneity. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer?, edited by Dr. Robert A. Gatenby.

Clonal Heterogeneity and Tumor Evolution: Past, Present, and the Future

Intratumor heterogeneity, which fosters tumor evolution, is a key challenge in cancer medicine. Here, we review data and technologies that have revealed intra-tumor heterogeneity across cancer types and the dynamics, constraints, and contingencies inherent to tumor evolution. We emphasize the importance of macro-evolutionary leaps, often involving large-scale chromosomal alterations, in driving tumor evolution and metastasis and consider the role of the tumor microenvironment in engendering heterogeneity and drug resistance. We suggest that bold approaches to drug development, harnessing the adaptive properties of the immune-microenvironment while limiting those of the tumor, combined with advances in clinical trial-design, will improve patient outcome.

Constraints in cancer evolution

Next-generation deep genome sequencing has only recently allowed us to quantitatively dissect the extent of heterogeneity within a tumour, resolving patterns of cancer evolution. Intratumour heterogeneity and natural selection contribute to resistance to anticancer therapies in the advanced setting. Recent evidence has also revealed that cancer evolution might be constrained. In this review, we discuss the origins of intratumour heterogeneity and subsequently focus on constraints imposed upon cancer evolution. The presence of (1) parallel evolution, (2) convergent evolution and (3) the biological impact of acquiring mutations in specific orders suggest that cancer evolution may be exploitable. These constraints on cancer evolution may help us identify cancer evolutionary rule books, which could eventually inform both diagnostic and therapeutic approaches to improve survival outcomes.

Targeting HIF2 in Clear Cell Renal Cell Carcinoma

Inactivation of the von Hippel-Lindau tumor-suppressor protein (pVHL) is the signature "truncal" event in clear cell renal cell carcinoma, which is the most common form of kidney cancer. pVHL is part of a ubiquitin ligase the targets the α subunit of the hypoxia-inducible factor (HIF) transcription factor for destruction when oxygen is available. Preclinical studies strongly suggest that deregulation of HIF, and particularly HIF2, drives pVHL-defective renal carcinogenesis. Although HIF2α was classically considered undruggable, structural and chemical work by Rick Bruick and Kevin Gardner at University of Texas Southwestern laid the foundation for the development of small molecule direct HIF2α antagonists (PT2385 and the related tool compound PT2399) by Peloton Therapeutics that block the dimerization of HIF2α with its partner protein ARNT1. These compounds inhibit clear cell renal cell carcinoma growth in preclinical models, and PT2385 has now entered the clinic. Nonetheless, the availability of such compounds, together with clustered regularly interspaced short palindromic repeat (CRISPR)-based gene editing approaches, has revealed a previously unappreciated heterogeneity among clear cell renal carcinomas and patient-derived xenografts with respect to HIF2 dependence, suggesting that predictive biomarkers will be needed to optimize the use of such agents in the clinic.

Challenges in molecular testing in non-small-cell lung cancer patients with advanced disease

Lung cancer diagnostics have progressed greatly in the previous decade. Development of molecular testing to identify an increasing number of potentially clinically actionable genetic variants, using smaller samples obtained via minimally invasive techniques, is a huge challenge. Tumour heterogeneity and cancer evolution in response to therapy means that repeat biopsies or circulating biomarkers are likely to be increasingly useful to adapt treatment as resistance develops. We highlight some of the current challenges faced in clinical practice for molecular testing of EGFR, ALK, and new biomarkers such as PDL1. Implementation of next generation sequencing platforms for molecular diagnostics in non-small-cell lung cancer is increasingly common, allowing testing of multiple genetic variants from a single sample. The use of next generation sequencing to recruit for molecularly stratified clinical trials is discussed in the context of the UK Stratified Medicine Programme and The UK National Lung Matrix Trial.

Whole Exome Sequencing Identifies TSC1/TSC2 Biallelic Loss as the Primary and Sufficient Driver Event for Renal Angiomyolipoma Development

Renal angiomyolipoma is a kidney tumor in the perivascular epithelioid (PEComa) family that is common in patients with Tuberous Sclerosis Complex (TSC) and Lymphangioleiomyomatosis (LAM) but occurs rarely sporadically. Though histologically benign, renal angiomyolipoma can cause life-threatening hemorrhage and kidney failure. Both angiomyolipoma and LAM have mutations in TSC2 or TSC1. However, the frequency and contribution of other somatic events in tumor development is unknown. We performed whole exome sequencing in 32 resected tumor samples (n = 30 angiomyolipoma, n = 2 LAM) from 15 subjects, including three with TSC. Two germline and 22 somatic inactivating mutations in TSC2 were identified, and one germline TSC1 mutation. Twenty of 32 (62%) samples showed copy neutral LOH (CN-LOH) in TSC2 or TSC1 with at least 8 different LOH regions, and 30 of 32 (94%) had biallelic loss of either TSC2 or TSC1. Whole exome sequencing identified a median of 4 somatic non-synonymous coding region mutations (other than in TSC2/TSC1), a mutation rate lower than nearly all other cancer types. Three genes with mutations were known cancer associated genes (BAP1, ARHGAP35 and SPEN), but they were mutated in a single sample each, and were missense variants with uncertain functional effects. Analysis of sixteen angiomyolipomas from a TSC subject showed both second hit point mutations and CN-LOH in TSC2, many of which were distinct, indicating that they were of independent clonal origin. However, three tumors had two shared mutations in addition to private somatic mutations, suggesting a branching evolutionary pattern of tumor development following initiating loss of TSC2. Our results indicate that TSC2 and less commonly TSC1 alterations are the primary essential driver event in angiomyolipoma/LAM, whereas other somatic mutations are rare and likely do not contribute to tumor development.

We performed comprehensive genome analysis of a kidney tumor called angiomyolipoma. These tumors are known to develop in most individuals who have Tuberous Sclerosis Complex (TSC) and those who have sporadic lymphangioleiomyomatosis (LAM), and are seen rarely in the general population. In these angiomyolipomas, we found consistent involvement of the TSC2 and TSC1 genes that are known to cause TSC, but very few (<5 on average) mutations elsewhere in the protein-coding regions. This is in stark contrast to other adult solid tumours that typically harbor hundreds to thousands of such mutations. Our results indicate that genetic alterations in TSC2/TSC1 are the primary and essential driver genetic events for development and progression of renal angiomyolipoma. Analysis of multiple angiomyolipomas from a single patient showed distinct genetic aberrations in the majority of samples, indicating that most of the tumors had developed independently. Branched clonal evolution was evident from the observation of three tumors that shared 2 mutations in addition to mutations private to each. Our results indicate that therapeutic approaches for treatment of patients with angiomyolipoma should focus on the consequences of TSC2/TSC1 loss, including but not limited to mTOR activation.

Independent Tumor Origin in Two Cases of Synchronous Bilateral Clear Cell Renal Cell Carcinoma

Bilateral renal cell carcinomas (RCCs) pose a challenge for clinical treatment and management. Most bilateral RCCs are sporadic, and do not show a hereditary pattern indicative of VHL syndrome or other inherited cancers. The origin and evolution of these sporadic bilateral RCCs remains elusive. We obtained normal and tumor samples from two male patients suffering from early stage synchronous bilateral clear cell RCC (ccRCC), and analyzed genomic DNA using whole exome sequencing and bisulfite pyrosequencing. We detected distinct 3p loss of heterozygosity (LOH) in both tumors in each patient. Two tumors within the same patient harbored distinct driver mutations and different CpG hypermethylation sites in the VHL promoter. Moreover, tumors exhibit independent evolutionary trajectories. Therefore, distinct 3p LOH, combined with contingent driver gene mutations and independent VHL hypermethylation, led to independent tumor origin and parallel evolution of bilateral ccRCC in these two patients. Our results indicate that tumors in these two cases were not due to common germline oncogenic mutations. They were results of multiple de novo mutations in each kidney, rather than primary ccRCC with contralateral renal metastasis. Therefore, histopathologic and genetic profiling from single tumor specimen may underestimate the mutational burden and somatic heterogeneity of bilateral ccRCCs.

Patients with genetically heterogeneous synchronous colorectal cancer carry rare damaging germline mutations in immune-related genes

Synchronous colorectal cancers (syCRCs) are physically separated tumours that develop simultaneously. To understand how the genetic and environmental background influences the development of multiple tumours, here we conduct a comparative analysis of 20 syCRCs from 10 patients. We show that syCRCs have independent genetic origins, acquire dissimilar somatic alterations, and have different clone composition. This inter- and intratumour heterogeneity must be considered in the selection of therapy and in the monitoring of resistance. SyCRC patients show a higher occurrence of inherited damaging mutations in immune-related genes compared to patients with solitary colorectal cancer and to healthy individuals from the 1,000 Genomes Project. Moreover, they have a different composition of immune cell populations in tumour and normal mucosa, and transcriptional differences in immune-related biological processes. This suggests an environmental field effect that promotes multiple tumours likely in the background of inflammation.

Some individuals present with multiple synchronous colorectal tumours, but the genetic understanding of this is unclear. Here, the authors use a sequencing strategy to show that the synchronous tumours are genetically independent and the patients harbour rare germline damaging mutations in genes associated with the immune system.