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.

Background Telomere maintenance is increasingly recognised as being fundamental to glioma oncogenesis with longer leucocyte 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 GWAS 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 analysed gene expression and DNA methylation data. Results Genetically increased LTL was significantly associated with increased glioma risk, IVW-RE ORSD 4.79 (95% CI: 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. Whilst complicated by patterns of linkage disequilibrium, genetic variation involving PARP1, PRRC2A, CARMIL1, POT1 and ATM-NPAT1 was implicated in the aetiology of glioma. Conclusions These observations extend the role of telomere-related genes in the development of glioma.

This paper introduces methodology for performing Bayesian inference sequentially on a sequence of posteriors on spaces of different dimensions. We show how this may be achieved through the use of sequential Monte Carlo (SMC) samplers (Del Moral et al., 2006, making use of the full flexibility of this framework in order that the method is computationally efficient. In particular, we introduce the innovation of using deterministic transformations to move particles effectively between target distributions with different dimensions. This approach, combined with adaptive methods, yields an extremely flexible and general algorithm for Bayesian model comparison that is suitable for use in applications where the acceptance rate in reversible jump Markov chain Monte Carlo (RJMCMC) is low. We demonstrate this approach on the well-studied problem of model comparison for mixture models, and for the novel application of inferring coalescent trees sequentially, as data arrives. arXiv:1612.06468v2 [stat.CO]

To characterise the somatic alterations in colorectal cancer (CRC), we conducted whole-genome sequencing analysis of 2,023 tumours. We provide the most detailed high-resolution map to date of somatic mutations in CRC, and demonstrate associations with clinicopathological features, in particular location in the large bowel. We refined the mutational processes and signatures acting in colorectal tumorigenesis. In analyses across the sample set or restricted to molecular subtypes, we identified 185 CRC driver genes, of which 117 were previously unreported. New drivers acted in various molecular pathways, including Wnt (CTNND1, AXIN1, TCF3), TGF-beta/BMP (TGFBR1) and MAP kinase (RASGRF1, RASA1, RAF1, and several MAP2K and MAP3K loci). Non-coding drivers included intronic neo-splice site alterations in APC and SMAD4. Whilst there was evidence of an excess of mutations in functionally active regions of the non-coding genome, no specific drivers were called with high confidence. Novel recurrent copy number changes included deletions of PIK3R1 and PWRN1, as well as amplification of CCND3 and NEDD9. Putative driver structural variants included BRD4 and SOX9 regulatory elements, and ACVR2A and ANKRD11 hotspot deletions. The frequencies of many driver mutations, including somatic Wnt and Ras pathway variants, showed a gradient along the colorectum. The Pks-pathogenic E. coli signature and TP53 mutations were primarily associated with rectal cancer. A set of unreported immune escape driver genes was found, primarily in hypermutated CRCs, most of which showed evidence of genetic evasion of the anti-cancer immune response. About 25% of cancers had a potentially actionable mutation for a known therapy. Thirty-three of the new driver genes were predicted to be essential, 17 possessed a druggable structure, and nine had a bioactive compound available. Our findings provide further insight into the genetics and biology of CRC, especially tumour subtypes defined by genomic instability or clinicopathological features.

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.

Identifying cancer driver genes is key for delivering the vision of precision oncology. The falling cost makes whole genome sequencing (WGS) attractive as a single all-encompassing test to identify uncommon cancer drivers in a patient, which may not be captured by panel testing but are potentially targeted by small molecules. We analysed 10,478 patient tumour genomes spanning 35 cancer types recruited to the UK 100,000 Genomes Project. We identified 330 driver genes, including 74 which are novel to any cancer. Across all cancer types 16% of the patients are eligible for a currently approved therapy. Computational chemogenomic analysis of cancer mutations identified 96 additional targets of compounds that are potentially active and represent candidates for future clinical trials, expanding opportunities for improved patient care.

To characterise the somatic alterations in colorectal cancer (CRC), we conducted whole-genome sequencing analysis of 2,023 tumours. We provide the most detailed high-resolution map to date of somatic mutations in CRC, and demonstrate associations with clinicopathological features, in particular location in the large bowel. We refined the mutational processes and signatures acting in colorectal tumorigenesis. In analyses across the sample set or restricted to molecular subtypes, we identified 185 CRC driver genes, of which 117 were previously unreported. New drivers acted in various molecular pathways, including Wnt (CTNND1, AXIN1, TCF3), TGF-β/BMP (TGFBR1) and MAP kinase (RASGRF1, RASA1, RAF1, and several MAP2K and MAP3K loci). Non-coding drivers included intronic neo-splice site alterations in APC and SMAD4. Whilst there was evidence of an excess of mutations in functionally active regions of the non-coding genome, no specific drivers were called with high confidence. Novel recurrent copy number changes included deletions of PIK3R1 and PWRN1, as well as amplification of CCND3 and NEDD9. Putative driver structural variants included BRD4 and SOX9 regulatory elements, and ACVR2A and ANKRD11 hotspot deletions. The frequencies of many driver mutations, including somatic Wnt and Ras pathway variants, showed a gradient along the colorectum. The Pks-pathogenic E. coli signature and TP53 mutations were primarily associated with rectal cancer. A set of unreported immune escape driver genes was found, primarily in hypermutated CRCs, most of which showed evidence of genetic evasion of the anti-cancer immune response. About 25% of cancers had a potentially actionable mutation for a known therapy. Thirty-three of the new driver genes were predicted to be essential, 17 possessed a druggable structure, and nine had a bioactive compound available. Our findings provide further insight into the genetics and biology of CRC, especially tumour subtypes defined by genomic instability or clinicopathological features.

For many cancers there are few well-established risk factors. Summary data from genome-wide association studies (GWAS) can be used in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to identify causal relationships. We performed a MR-PheWAS of breast, prostate, colorectal, lung, endometrial, oesophageal, renal, and ovarian cancers, comprising 378,142 cases and 485,715 controls. To derive a more comprehensive insight into disease aetiology we systematically mined the literature space for supporting evidence. We evaluated causal relationships for over 3,000 potential risk factors. In addition to identifying well-established risk factors (smoking, alcohol, obesity, lack of physical activity), we provide evidence for specific factors, including dietary intake, sex steroid hormones, plasma lipids and telomere length as determinants of cancer risk. We also implicate molecular factors including plasma levels of IL-18, LAG-3, IGF-1, CT-1, and PRDX1 as risk factors. Our analyses highlight the importance of risk factors that are common to many cancer types but also reveal aetiological differences. A number of the molecular factors we identify have the potential to be biomarkers. Our findings should aid public health prevention strategies to reduce cancer burden. We provide a R/Shiny app (https://mrcancer.shinyapps.io/mrcan/) to visualise findings.