Analysis of 10,478 cancer genomes identifies candidate driver genes and opportunities for precision oncology.

Tumor genomic profiling is increasingly seen as a prerequisite to guide the treatment of patients with cancer. To explore the value of whole-genome sequencing (WGS) in broadening the scope of cancers potentially amenable to a precision therapy, we analysed whole-genome sequencing data on 10,478 patients spanning 35 cancer types recruited to the UK 100,000 Genomes Project. We identified 330 candidate driver genes, including 74 that are new to any cancer. We estimate that approximately 55% of patients studied harbor at least one clinically relevant mutation, predicting either sensitivity or resistance to certain treatments or clinical trial eligibility. By performing computational chemogenomic analysis of cancer mutations we identify additional targets for compounds that represent attractive candidates for future clinical trials. This study represents one of the most comprehensive efforts thus far to identify cancer driver genes in the real world setting and assess their impact on informing precision oncology.

Genetic landscape of interval and screen detected breast cancer.

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

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

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

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features.

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

Frequency of pathogenic germline variants in cancer susceptibility genes in 1336 renal cell carcinoma cases.

Renal cell carcinoma (RCC) occurs in a number of cancer predisposition syndromes, but the genetic architecture of susceptibility to RCC is not well defined. We investigated the frequency of pathogenic and likely pathogenic (P/LP) germline variants in cancer susceptibility genes (CSGs) within a large series of unselected RCC participants. Whole-genome sequencing data on 1336 RCC participants and 5834 controls recruited to the UK 100 000 Genomes Project, a nationwide multicentre study, was analyzed to identify rare P/LP short variants (single nucleotide variants and insertions/deletions ranging from 1 to 50 base pairs) and structural variants in 121 CSGs. Among 1336 RCC participants [mean: 61.3 years (±12 SD), range: 13-88 years; 64% male], 85 participants [6.4%; 95% CI (5.1, 7.8)] had one or more P/LP germline variant in a wider range of CSGs than previously recognized. A further 64 intragenic variants in CSGs previously associated with RCC were classified as a variant of uncertain significance (VUS) (24 'hot VUSs') and were considered to be of potential clinical relevance as further evaluation might results in their reclassification. Most patients with P variants in well-established CSGs known to predispose to renal cell carcinoma (RCC-CSGs) were aged <50 years. Burden test analysis for filtered variants in CSGs demonstrated a significant excess of CHEK2 variants in European RCC participants compared with the healthy European controls (P = 0.0019). Approximately, 6% of the patients with RCC unselected for family history have a germline variant requiring additional follow-up analysis. To improve diagnostic yield, we suggest expanding the panel of RCC-CSGs tested to include CHEK2 and all SDHx subunits and raising the eligibility criteria for age-based testing.

Elongin C (ELOC/TCEB1)-associated von Hippel-Lindau disease.

Around 95% of patients with clinical features that meet the diagnostic criteria for von Hippel-Lindau disease (VHL) have a detectable inactivating germline variant in VHL. The VHL protein (pVHL) functions as part of the E3 ubiquitin ligase complex comprising pVHL, elongin C, elongin B, cullin 2 and ring box 1 (VCB-CR complex), which plays a key role in oxygen sensing and degradation of hypoxia-inducible factors. To date, only variants in VHL have been shown to cause VHL disease. We undertook trio analysis by whole-exome sequencing in a proband with VHL disease but without a detectable VHL mutation. Molecular studies were also performed on paired DNA extracted from the proband's kidney tumour and blood and bioinformatics analysis of sporadic renal cell carcinoma (RCC) dataset was undertaken. A de novo pathogenic variant in ELOC NM_005648.4(ELOC):c.236A>G (p.Tyr79Cys) gene was identified in the proband. ELOC encodes elongin C, a key component [C] of the VCB-CR complex. The p.Tyr79Cys substitution is a mutational hotspot in sporadic VHL-competent RCC and has previously been shown to mimic the effects of pVHL deficiency on hypoxic signalling. Analysis of an RCC from the proband showed similar findings to that in somatically ELOC-mutated RCC (expression of hypoxia-responsive proteins, no somatic VHL variants and chromosome 8 loss). These findings are consistent with pathogenic ELOC variants being a novel cause for VHL disease and suggest that genetic testing for ELOC variants should be performed in individuals with suspected VHL disease with no detectable VHL variant.

Signatures of TOP1 transcription-associated mutagenesis in cancer and germline.

The mutational landscape is shaped by many processes. Genic regions are vulnerable to mutation but are preferentially protected by transcription-coupled repair1. In microorganisms, transcription has been demonstrated to be mutagenic2,3; however, the impact of transcription-associated mutagenesis remains to be established in higher eukaryotes4. Here we show that ID4-a cancer insertion-deletion (indel) mutation signature of unknown aetiology5 characterized by short (2 to 5 base pair) deletions -is due to a transcription-associated mutagenesis process. We demonstrate that defective ribonucleotide excision repair in mammals is associated with the ID4 signature, with mutations occurring at a TNT sequence motif, implicating topoisomerase 1 (TOP1) activity at sites of genome-embedded ribonucleotides as a mechanistic basis. Such TOP1-mediated deletions occur somatically in cancer, and the ID-TOP1 signature is also found in physiological settings, contributing to genic de novo indel mutations in the germline. Thus, although topoisomerases protect against genome instability by relieving topological stress6, their activity may also be an important source of mutations in the human genome.

Relationship between genetically determined telomere length and glioma risk.

BACKGROUND: Telomere maintenance is increasingly recognized as being fundamental to glioma oncogenesis with longer leukocyte telomere length (LTL) reported to increase risk of glioma. To gain further insight into the relationship between telomere genetics and risk of glioma, we conducted several complementary analyses, using genome-wide association studies data on LTL (78 592 individuals) and glioma (12 488 cases and 18 169 controls). METHODS: We performed both classical and summary Mendelian randomization (SMR), coupled with heterogeneity in dependent instruments tests, at genome-wide significant LTL loci to examine if an association was mediated by the same causal variant in glioma. To prioritize genes underscoring glioma-LTL associations, we analyzed gene expression and DNA methylation data. RESULTS: Genetically increased LTL was significantly associated with increased glioma risk, random-effects inverse variance weighted ORs per 1 SD unit increase in the putative risk factor (odds ratio [OR]SD) 4.79 (95% confidence interval: 2.11-10.85; P = 1.76 × 10-4). SMR confirmed the previously reported LTL associations at 3q26.2 (TERC; PSMR = 1.33 × 10-5), 5p15.33 (TERT; PSMR = 9.80 × 10-27), 10q24.33 (STN1 alias OBFC1; PSMR = 4.31 × 10-5), and 20q13.3 (STMN3/RTEL1; PSMR = 2.47 × 10-4) glioma risk loci. Our analysis implicates variation at 1q42.12 (PSMR = 1.55 × 10-2), 6p21.3 (PSMR = 9.76 × 10-3), 6p22.2 (PSMR = 5.45 × 10-3), 7q31.33 (PSMR = 6.52 × 10-3), and 11q22.3 (PSMR = 8.89 × 10-4) as risk factors for glioma risk. While complicated by patterns of linkage disequilibrium, genetic variation involving PARP1, PRRC2A, CARMIL1, POT1, and ATM-NPAT1 was implicated in the etiology of glioma. CONCLUSIONS: These observations extend the role of telomere-related genes in the development of glioma.

Cancer drivers and clonal dynamics in acute lymphoblastic leukaemia subtypes.

To obtain a comprehensive picture of composite genetic driver events and clonal dynamics in subtypes of paediatric acute lymphoblastic leukaemia (ALL) we analysed tumour-normal whole genome sequencing and expression data from 361 newly diagnosed patients. We report the identification of both structural drivers, as well as recurrent non-coding variation in promoters. Additionally we found the transcriptional profile of histone gene cluster 1 and CTCF altered tumours shared hallmarks of hyperdiploid ALL suggesting a 'hyperdiploid like' subtype. ALL subtypes are driven by distinct mutational processes with AID mutagenesis being confined to ETV6-RUNX1 tumours. Subclonality is a ubiquitous feature of ALL, consistent with Darwinian evolution driving selection and expansion of tumours. Driver mutations in B-cell developmental genes (IKZF1, PAX5, ZEB2) tend to be clonal and RAS/RTK mutations subclonal. In addition to identifying new avenues for therapeutic exploitation, this analysis highlights that targeted therapies should take into account composite mutational profile and clonality.

Lack of an association between gallstone disease and bilirubin levels with risk of colorectal cancer: a Mendelian randomisation analysis.

BACKGROUND: Epidemiological studies of the relationship between gallstone disease and circulating levels of bilirubin with risk of developing colorectal cancer (CRC) have been inconsistent. To address possible confounding and reverse causation, we examine the relationship between these potential risk factors and CRC using Mendelian randomisation (MR). METHODS: We used two-sample MR to examine the relationship between genetic liability to gallstone disease and circulating levels of bilirubin with CRC in 26,397 patients and 41,481 controls. We calculated the odds ratio per genetically predicted SD unit increase in log bilirubin levels (ORSD) for CRC and tested for a non-zero causal effect of gallstones on CRC. Sensitivity analysis was applied to identify violations of estimator assumptions. RESULTS: No association between either gallstone disease (P value = 0.60) or circulating levels of bilirubin (ORSD = 1.00, 95% confidence interval (CI) = 0.96-1.03, P value = 0.90) with CRC was shown. CONCLUSIONS: Despite the large scale of this study, we found no evidence for a causal relationship between either circulating levels of bilirubin or gallstone disease with risk of developing CRC. While the magnitude of effect suggested by some observational studies can confidently be excluded, we cannot exclude the possibility of smaller effect sizes and non-linear relationships.

Searching for causal relationships of glioma: a phenome-wide Mendelian randomisation study.

BACKGROUND: The aetiology of glioma is poorly understood. Summary data from genome-wide association studies (GWAS) can be used in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to search for glioma risk factors. METHODS: We performed an MR-PheWAS analysing 316 phenotypes, proxied by 8387 genetic variants, and summary genetic data from a GWAS of 12,488 glioma cases and 18,169 controls. Causal effects were estimated under a random-effects inverse-variance-weighted (IVW-RE) model, with robust adjusted profile score (MR-RAPS), weighted median and mode-based estimates computed to assess the robustness of findings. Odds ratios per one standard deviation increase in each phenotype were calculated for all glioma, glioblastoma (GBM) and non-GBM tumours. RESULTS: No significant associations (P < 1.58 × 10-4) were observed between phenotypes and glioma under the IVW-RE model. Suggestive associations (1.58 × 10-4 < P < 0.05) were observed between leukocyte telomere length (LTL) with all glioma (ORSD = 3.91, P = 9.24 × 10-3) and GBM (ORSD = 4.86, P = 3.23 × 10-2), but the association was primarily driven by the TERT variant rs2736100. Serum low-density lipoprotein cholesterol and plasma HbA1C showed suggestive associations with glioma (ORSD = 1.11, P = 1.39 × 10-2 and ORSD = 1.28, P = 1.73 × 10-2, respectively), both associations being reliant on single genetic variants. CONCLUSIONS: Our study provides further insight into the aetiological basis of glioma for which published data have been mixed.

An enhanced genetic model of relapsed IGH-translocated multiple myeloma evolutionary dynamics.

Most patients with multiple myeloma (MM) die from progressive disease after relapse. To advance our understanding of MM evolution mechanisms, we performed whole-genome sequencing of 80 IGH-translocated tumour-normal newly diagnosed pairs and 24 matched relapsed tumours from the Myeloma XI trial. We identify multiple events as potentially important for survival and therapy-resistance at relapse including driver point mutations (e.g., TET2), translocations (MAP3K14), lengthened telomeres, and increased genomic instability (e.g., 17p deletions). Despite heterogeneous mutational processes contributing to relapsed mutations across MM subtypes, increased AID/APOBEC activity is particularly associated with shorter progression time to relapse, and contributes to higher mutational burden at relapse. In addition, we identify three enhanced major clonal evolution patterns of MM relapse, independent of treatment strategies and molecular karyotypes, questioning the viability of "evolutionary herding" approach in treating drug-resistant MM. Our data show that MM relapse is associated with acquisition of new mutations and clonal selection, and suggest APOBEC enzymes among potential targets for therapy-resistant MM.

Genomic landscape of platinum resistant and sensitive testicular cancers.

While most testicular germ cell tumours (TGCTs) exhibit exquisite sensitivity to platinum chemotherapy, ~10% are platinum resistant. To gain insight into the underlying mechanisms, we undertake whole exome sequencing and copy number analysis in 40 tumours from 26 cases with platinum-resistant TGCT, and combine this with published genomic data on an additional 624 TGCTs. We integrate analyses for driver mutations, mutational burden, global, arm-level and focal copy number (CN) events, and SNV and CN signatures. Albeit preliminary and observational in nature, these analyses provide support for a possible mechanistic link between early driver mutations in RAS and KIT and the widespread copy number events by which TGCT is characterised.

Search for multiple myeloma risk factors using Mendelian randomization.

The etiology of multiple myeloma (MM) is poorly understood. Summary data from genome-wide association studies (GWASs) of multiple phenotypes can be exploited in a Mendelian randomization (MR) phenome-wide association study (PheWAS) to search for factors influencing MM risk. We performed an MR-PheWAS analyzing 249 phenotypes, proxied by 10 225 genetic variants, and summary genetic data from a GWAS of 7717 MM cases and 29 304 controls. Odds ratios (ORs) per 1 standard deviation increase in each phenotype were estimated under an inverse variance weighted random effects model. A Bonferroni-corrected threshold of P = 2 × 10-4 was considered significant, whereas P < .05 was considered suggestive of an association. Although no significant associations with MM risk were observed among the 249 phenotypes, 28 phenotypes showed evidence suggestive of association, including increased levels of serum vitamin B6 and blood carnitine (P = 1.1 × 10-3) with greater MM risk and ω-3 fatty acids (P = 5.4 × 10-4) with reduced MM risk. A suggestive association between increased telomere length and reduced MM risk was also noted; however, this association was primarily driven by the previously identified risk variant rs10936599 at 3q26 (TERC). Although not statistically significant, increased body mass index was associated with increased risk (OR, 1.10; 95% confidence interval, 0.99-1.22), supporting findings from a previous meta-analysis of prospective observational studies. Our study did not provide evidence supporting any modifiable factors examined as having a major influence on MM risk; however, it provides insight into factors for which the evidence has previously been mixed.

Lack of association between modifiable exposures and glioma risk: a Mendelian randomization analysis.

BACKGROUND: The etiological basis of glioma is poorly understood. We have used genetic markers in a Mendelian randomization (MR) framework to examine if lifestyle, cardiometabolic, and inflammatory factors influence the risk of glioma. This methodology reduces bias from confounding and is not affected by reverse causation. METHODS: We identified genetic instruments for 37 potentially modifiable risk factors and evaluated their association with glioma risk using data from a genome-wide association study of 12 488 glioma patients and 18 169 controls. We used the estimated odds ratio of glioma associated with each of the genetically defined traits to infer evidence for a causal relationship with the following exposures:Lifestyle and dietary factors-height, plasma insulin-like growth factor 1, blood carnitine, blood methionine, blood selenium, blood zinc, circulating adiponectin, circulating carotenoids, iron status, serum calcium, vitamins (A1, B12, B6, E, and 25-hydroxyvitamin D), fatty acid levels (monounsaturated, omega-3, and omega-6) and circulating fetuin-A;Cardiometabolic factors-birth weight, high density lipoprotein cholesterol, low density lipoprotein cholesterol, total cholesterol, total triglycerides, basal metabolic rate, body fat percentage, body mass index, fasting glucose, fasting proinsulin, glycated hemoglobin levels, diastolic and systolic blood pressure, waist circumference, waist-to-hip ratio; andInflammatory factors- C-reactive protein, plasma interleukin-6 receptor subunit alpha and serum immunoglobulin E. RESULTS: After correction for the testing of multiple potential risk factors and excluding associations driven by one single nucleotide polymorphism, no significant association with glioma risk was observed (ie, PCorrected > 0.05). CONCLUSIONS: This study did not provide evidence supporting any of the 37 factors examined as having a significant influence on glioma risk.

Modifiable pathways for colorectal cancer: a mendelian randomisation analysis.

BACKGROUND: Epidemiological studies have linked lifestyle, cardiometabolic, reproductive, developmental, and inflammatory factors to the risk of colorectal cancer. However, which specific factors affect risk and the strength of these effects are unknown. We aimed to examine the relationship between potentially modifiable risk factors and colorectal cancer. METHODS: We used a random-effects model to examine the relationship between 39 potentially modifiable risk factors and colorectal cancer in 26 397 patients with colorectal cancer and 41 481 controls (ie, people without colorectal cancer). These population data came from a genome-wide association study of people of European ancestry, which was amended to exclude UK BioBank data. In the model, we used genetic variants as instruments via two-sample mendelian randomisation to limit bias from confounding and reverse causation. We calculated odds ratios per genetically predicted SD unit increase in each putative risk factor (ORSD) for colorectal cancer risk. We did mendelian randomisation Egger regressions to identify evidence of potential violations of mendelian randomisation assumptions. A Bonferroni-corrected threshold of p=1·3 × 10-3 was considered significant, and p values less than 0·05 were considered to be suggestive of an association. FINDINGS: No putative risk factors were significantly associated with colorectal cancer risk after correction for multiple testing. However, suggestive associations with increased risk were noted for genetically predicted body fat percentage (ORSD 1·14 [95% CI 1·03-1·25]; p=0·0086), body-mass index (1·09 [1·01-1·17]; p=0·023), waist circumference (1·13 [1·02-1·26]; p=0·018), basal metabolic rate (1·10 [1·03-1·18]; p=0·0079), and concentrations of LDL cholesterol (1·14 [1·04-1·25]; p=0·0056), total cholesterol (1·09 [1·01-1·18]; p=0·025), circulating serum iron (1·17 [1·00-1·36]; p=0·049), and serum vitamin B12 (1·21 [1·04-1·42]; p=0·016), although potential pleiotropy among genetic variants used as instruments for vitamin B12 constrains the finding. A suggestive association was also noted between adult height and increased risk of colorectal cancer (ORSD 1·04 [95% CI 1·00-1·08]; p=0·032). Low blood selenium concentration had a suggestive association with decreased risk of colorectal cancer (ORSD 0·85 [95% CI 0·75-0·96]; p=0·0078) based on a single variant, as did plasma concentrations of interleukin-6 receptor subunit α (also based on a single variant; 0·98 [0·96-1·00]; p=0·035). Risk of colorectal cancer was not associated with any sex hormone or reproductive factor, serum calcium, or circulating 25-hydroxyvitamin D concentrations. INTERPRETATION: This analysis identified several modifiable targets for primary prevention of colorectal cancer, including lifestyle, obesity, and cardiometabolic factors, that should inform public health policy. FUNDING: Cancer Research UK, UK Medical Research Council Human Genetics Unit Centre, DJ Fielding Medical Research Trust, EU COST Action, and the US National Cancer Institute.

Mutational processes contributing to the development of multiple myeloma.

To gain insight into multiple myeloma (MM) tumorigenesis, we analyzed the mutational signatures in 874 whole-exome and 850 whole-genome data from the CoMMpass Study. We identified that coding and non-coding regions are differentially dominated by distinct single-nucleotide variant (SNV) mutational signatures, as well as five de novo structural rearrangement signatures. Mutational signatures reflective of different principle mutational processes-aging, defective DNA repair, and apolipoprotein B editing complex (APOBEC)/activation-induced deaminase activity-characterize MM. These mutational signatures show evidence of subgroup specificity-APOBEC-attributed signatures associated with MAF translocation t(14;16) and t(14;20) MM; potentially DNA repair deficiency with t(11;14) and t(4;14); and aging with hyperdiploidy. Mutational signatures beyond that associated with APOBEC are independent of established prognostic markers and appear to have relevance to predicting high-risk MM.