A pan-cancer analysis of the microbiome in metastatic cancer.

Microbial communities are resident to multiple niches of the human body and are important modulators of the host immune system and responses to anticancer therapies. Recent studies have shown that complex microbial communities are present within primary tumors. To investigate the presence and relevance of the microbiome in metastases, we integrated mapping and assembly-based metagenomics, genomics, transcriptomics, and clinical data of 4,160 metastatic tumor biopsies. We identified organ-specific tropisms of microbes, enrichments of anaerobic bacteria in hypoxic tumors, associations between microbial diversity and tumor-infiltrating neutrophils, and the association of Fusobacterium with resistance to immune checkpoint blockade (ICB) in lung cancer. Furthermore, longitudinal tumor sampling revealed temporal evolution of the microbial communities and identified bacteria depleted upon ICB. Together, we generated a pan-cancer resource of the metastatic tumor microbiome that may contribute to advancing treatment strategies.

EXO1 protects BRCA1-deficient cells against toxic DNA lesions.

Inactivating mutations in the BRCA1 and BRCA2 genes impair DNA double-strand break (DSB) repair by homologous recombination (HR), leading to chromosomal instability and cancer. Importantly, BRCA1/2 deficiency also causes therapeutically targetable vulnerabilities. Here, we identify the dependency on the end resection factor EXO1 as a key vulnerability of BRCA1-deficient cells. EXO1 deficiency generates poly(ADP-ribose)-decorated DNA lesions during S phase that associate with unresolved DSBs and genomic instability in BRCA1-deficient but not in wild-type or BRCA2-deficient cells. Our data indicate that BRCA1/EXO1 double-deficient cells accumulate DSBs due to impaired repair by single-strand annealing (SSA) on top of their HR defect. In contrast, BRCA2-deficient cells retain SSA activity in the absence of EXO1 and hence tolerate EXO1 loss. Consistent with a dependency on EXO1-mediated SSA, we find that BRCA1-mutated tumors show elevated EXO1 expression and increased SSA-associated genomic scars compared with BRCA1-proficient tumors. Overall, our findings uncover EXO1 as a promising therapeutic target for BRCA1-deficient tumors.

Pan-cancer whole-genome comparison of primary and metastatic solid tumours.

Metastatic cancer remains an almost inevitably lethal disease1-3. A better understanding of disease progression and response to therapies therefore remains of utmost importance. Here we characterize the genomic differences between early-stage untreated primary tumours and late-stage treated metastatic tumours using a harmonized pan-cancer analysis (or reanalysis) of two unpaired primary4 and metastatic5 cohorts of 7,108 whole-genome-sequenced tumours. Metastatic tumours in general have a lower intratumour heterogeneity and a conserved karyotype, displaying only a modest increase in mutations, although frequencies of structural variants are elevated overall. Furthermore, highly variable tumour-specific contributions of mutational footprints of endogenous (for example, SBS1 and APOBEC) and exogenous mutational processes (for example, platinum treatment) are present. The majority of cancer types had either moderate genomic differences (for example, lung adenocarcinoma) or highly consistent genomic portraits (for example, ovarian serous carcinoma) when comparing early-stage and late-stage disease. Breast, prostate, thyroid and kidney renal clear cell carcinomas and pancreatic neuroendocrine tumours are clear exceptions to the rule, displaying an extensive transformation of their genomic landscape in advanced stages. Exposure to treatment further scars the tumour genome and introduces an evolutionary bottleneck that selects for known therapy-resistant drivers in approximately half of treated patients. Our data showcase the potential of pan-cancer whole-genome analysis to identify distinctive features of late-stage tumours and provide a valuable resource to further investigate the biological basis of cancer and resistance to therapies.

The genome-wide mutational consequences of DNA hypomethylation.

DNA methylation is important for establishing and maintaining cell identity and for genomic stability. This is achieved by regulating the accessibility of regulatory and transcriptional elements and the compaction of subtelomeric, centromeric, and other inactive genomic regions. Carcinogenesis is accompanied by a global loss in DNA methylation, which facilitates the transformation of cells. Cancer hypomethylation may also cause genomic instability, for example through interference with the protective function of telomeres and centromeres. However, understanding the role(s) of hypomethylation in tumor evolution is incomplete because the precise mutational consequences of global hypomethylation have thus far not been systematically assessed. Here we made genome-wide inventories of all possible genetic variation that accumulates in single cells upon the long-term global hypomethylation by CRISPR interference-mediated conditional knockdown of DNMT1. Depletion of DNMT1 resulted in a genomewide reduction in DNA methylation. The degree of DNA methylation loss was similar to that observed in many cancer types. Hypomethylated cells showed reduced proliferation rates, increased transcription of genes, reactivation of the inactive X-chromosome and abnormal nuclear morphologies. Prolonged hypomethylation was accompanied by increased chromosomal instability. However, there was no increase in mutational burden, enrichment for certain mutational signatures or accumulation of structural variation to the genome. In conclusion, the primary consequence of hypomethylation is genomic instability, which in cancer leads to increased tumor heterogeneity and thereby fuels cancer evolution.

Common anti-cancer therapies induce somatic mutations in stem cells of healthy tissue.

Genome-wide mutation analyses have revealed that specific anti-cancer drugs are highly mutagenic to cancer cells, but the mutational impact of anti-cancer therapies on normal cells is not known. Here, we examine genome-wide somatic mutation patterns in 42 healthy adult stem cells (ASCs) of the colon or the liver from 14 cancer patients (mean of 3.2 ASC per donor) that received systemic chemotherapy and/or local radiotherapy. The platinum-based chemo-drug Oxaliplatin induces on average 535 ± 260 mutations in colon ASC, while 5-FU shows a complete mutagenic absence in most, but not all colon ASCs. In contrast with the colon, normal liver ASCs escape mutagenesis from systemic treatment with Oxaliplatin and 5-FU. Thus, while chemotherapies are highly effective at killing cancer cells, their systemic use also increases the mutational burden of long-lived normal stem cells responsible for tissue renewal thereby increasing the risk for developing second cancers.

Homologous Recombination Deficiency Scar: Mutations and Beyond-Implications for Precision Oncology.

Homologous recombination deficiency (HRD) is a prevalent in approximately 17% of tumors and is associated with enhanced sensitivity to anticancer therapies inducing double-strand DNA breaks. Accurate detection of HRD would therefore allow improved patient selection and outcome of conventional and targeted anticancer therapies. However, current clinical assessment of HRD mainly relies on determining germline BRCA1/2 mutational status and is insufficient for adequate patient stratification as mechanisms of HRD occurrence extend beyond functional BRCA1/2 loss. HRD, regardless of BRCA1/2 status, is associated with specific forms of genomic and mutational signatures termed HRD scar. Detection of this HRD scar might therefore be a more reliable biomarker for HRD. This review discusses and compares different methods of assessing HRD and HRD scar, their advances into the clinic, and their potential implications for precision oncology.

Machine learning-based tissue of origin classification for cancer of unknown primary diagnostics using genome-wide mutation features.

Cancers of unknown primary (CUP) origin account for ∼3% of all cancer diagnoses, whereby the tumor tissue of origin (TOO) cannot be determined. Using a uniformly processed dataset encompassing 6756 whole-genome sequenced primary and metastatic tumors, we develop Cancer of Unknown Primary Location Resolver (CUPLR), a random forest TOO classifier that employs 511 features based on simple and complex somatic driver and passenger mutations. CUPLR distinguishes 35 cancer (sub)types with ∼90% recall and ∼90% precision based on cross-validation and test set predictions. We find that structural variant derived features increase the performance and utility for classifying specific cancer types. With CUPLR, we could determine the TOO for 82/141 (58%) of CUP patients. Although CUPLR is based on machine learning, it provides a human interpretable graphical report with detailed feature explanations. The comprehensive output of CUPLR complements existing histopathological procedures and can enable improved diagnostics for CUP patients.

Functional RECAP (REpair CAPacity) assay identifies homologous recombination deficiency undetected by DNA-based BRCAness tests.

Germline BRCA1/2 mutation status is predictive for response to Poly-[ADP-Ribose]-Polymerase (PARP) inhibitors in breast cancer (BC) patients. However, non-germline BRCA1/2 mutated and homologous recombination repair deficient (HRD) tumors are likely also PARP-inhibitor sensitive. Clinical validity and utility of various HRD biomarkers are under investigation. The REpair CAPacity (RECAP) test is a functional method to select HRD tumors based on their inability to form RAD51 foci. We investigated whether this functional test defines a similar group of HRD tumors as DNA-based tests. An HRD enriched cohort (n = 71; 52 primary and 19 metastatic BCs) selected based on the RECAP test (26 RECAP-HRD; 37%), was subjected to DNA-based HRD tests (i.e., Classifier of HOmologous Recombination Deficiency (CHORD) and BRCA1/2-like classifier). Whole genome sequencing (WGS) was carried out for 38 primary and 19 metastatic BCs. The RECAP test identified all bi-allelic BRCA deficient samples (n = 15) in this cohort. RECAP status partially correlated with DNA-based HRD test outcomes (70% concordance for both RECAP-CHORD and RECAP-BRCA1/2-like classifier). RECAP selected additional samples unable to form RAD51 foci, suggesting that this functional assay identified deficiencies in other DNA repair genes, which could also result in PARP-inhibitor sensitivity. Direct comparison of these HRD tests in clinical trials will be required to evaluate the optimal predictive test for clinical decision making.

Distinct Genomic Profiles Are Associated with Treatment Response and Survival in Ovarian Cancer.

The majority of patients with ovarian cancer ultimately develop recurrent chemotherapy-resistant disease. Treatment stratification is mainly based on histological subtype and stage, prior response to platinum-based chemotherapy, and time to recurrent disease. Here, we integrated clinical treatment, treatment response, and survival data with whole-genome sequencing profiles of 132 solid tumor biopsies of metastatic epithelial ovarian cancer to explore genome-informed stratification opportunities. Samples from primary and recurrent disease harbored comparable numbers of single nucleotide variants and structural variants. Mutational signatures represented platinum exposure, homologous recombination deficiency, and aging. Unsupervised hierarchical clustering based on genomic input data identified specific ovarian cancer subgroups, characterized by homologous recombination deficiency, genome stability, and duplications. The clusters exhibited distinct response rates and survival probabilities which could thus potentially be used for genome-informed therapy stratification for more personalized ovarian cancer treatment.

MutationalPatterns: the one stop shop for the analysis of mutational processes.

BACKGROUND: The collective of somatic mutations in a genome represents a record of mutational processes that have been operative in a cell. These processes can be investigated by extracting relevant mutational patterns from sequencing data. RESULTS: Here, we present the next version of MutationalPatterns, an R/Bioconductor package, which allows in-depth mutational analysis of catalogues of single and double base substitutions as well as small insertions and deletions. Major features of the package include the possibility to perform regional mutation spectra analyses and the possibility to detect strand asymmetry phenomena, such as lesion segregation. On top of this, the package also contains functions to determine how likely it is that a signature can cause damaging mutations (i.e., mutations that affect protein function). This updated package supports stricter signature refitting on known signatures in order to prevent overfitting. Using simulated mutation matrices containing varied signature contributions, we showed that reliable refitting can be achieved even when only 50 mutations are present per signature. Additionally, we incorporated bootstrapped signature refitting to assess the robustness of the signature analyses. Finally, we applied the package on genome mutation data of cell lines in which we deleted specific DNA repair processes and on large cancer datasets, to show how the package can be used to generate novel biological insights. CONCLUSIONS: This novel version of MutationalPatterns allows for more comprehensive analyses and visualization of mutational patterns in order to study the underlying processes. Ultimately, in-depth mutational analyses may contribute to improved biological insights in mechanisms of mutation accumulation as well as aid cancer diagnostics. MutationalPatterns is freely available at http://bioconductor.org/packages/MutationalPatterns .

The Mutagenic Impact of Environmental Exposures in Human Cells and Cancer: Imprints Through Time.

During life, the DNA of our cells is continuously exposed to external damaging processes. Despite the activity of various repair mechanisms, DNA damage eventually results in the accumulation of mutations in the genomes of our cells. Oncogenic mutations are at the root of carcinogenesis, and carcinogenic agents are often highly mutagenic. Over the past decade, whole genome sequencing data of healthy and tumor tissues have revealed how cells in our body gradually accumulate mutations because of exposure to various mutagenic processes. Dissection of mutation profiles based on the type and context specificities of the altered bases has revealed a variety of signatures that reflect past exposure to environmental mutagens, ranging from chemotherapeutic drugs to genotoxic gut bacteria. In this review, we discuss the latest knowledge on somatic mutation accumulation in human cells, and how environmental mutagenic factors further shape the mutation landscapes of tissues. In addition, not all carcinogenic agents induce mutations, which may point to alternative tumor-promoting mechanisms, such as altered clonal selection dynamics. In short, we provide an overview of how environmental factors induce mutations in the DNA of our healthy cells and how this contributes to carcinogenesis. A better understanding of how environmental mutagens shape the genomes of our cells can help to identify potential preventable causes of cancer.

GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing.

GRIDSS2 is the first structural variant caller to explicitly report single breakends-breakpoints in which only one side can be unambiguously determined. By treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints, GRIDSS2 can explain 47% of somatic centromere copy number changes using single breakends to non-centromere sequence. On a cohort of 3782 deeply sequenced metastatic cancers, GRIDSS2 achieves an unprecedented 3.1% false negative rate and 3.3% false discovery rate and identifies a novel 32-100 bp duplication signature. GRIDSS2 simplifies complex rearrangement interpretation through phasing of structural variants with 16% of somatic calls phasable using paired-end sequencing.

Pan-cancer landscape of homologous recombination deficiency.

Homologous recombination deficiency (HRD) results in impaired double strand break repair and is a frequent driver of tumorigenesis. Here, we develop a genome-wide mutational scar-based pan-cancer Classifier of HOmologous Recombination Deficiency (CHORD) that can discriminate BRCA1- and BRCA2-subtypes. Analysis of a metastatic (n = 3,504) and primary (n = 1,854) pan-cancer cohort reveals that HRD is most frequent in ovarian and breast cancer, followed by pancreatic and prostate cancer. We identify biallelic inactivation of BRCA1, BRCA2, RAD51C or PALB2 as the most common genetic cause of HRD, with RAD51C and PALB2 inactivation resulting in BRCA2-type HRD. We find that while the specific genetic cause of HRD is cancer type specific, biallelic inactivation is predominantly associated with loss-of-heterozygosity (LOH), with increased contribution of deep deletions in prostate cancer. Our results demonstrate the value of pan-cancer genomics-based HRD testing and its potential diagnostic value for patient stratification towards treatment with e.g. poly ADP-ribose polymerase inhibitors (PARPi).

Single-cell derived tumor organoids display diversity in HLA class I peptide presentation.

Tumor heterogeneity is a major cause of therapeutic resistance. Immunotherapy may exploit alternative vulnerabilities of drug-resistant cells, where tumor-specific human leukocyte antigen (HLA) peptide ligands are promising leads to invoke targeted anti-tumor responses. Here, we investigate the variability in HLA class I peptide presentation between different clonal cells of the same colorectal cancer patient, using an organoid system. While clone-specific differences in HLA peptide presentation were observed, broad inter-clone variability was even more prevalent (15-25%). By coupling organoid proteomics and HLA peptide ligandomics, we also found that tumor-specific ligands from DNA damage control and tumor suppressor source proteins were prominently presented by tumor cells, coinciding likely with the silencing of such cytoprotective functions. Collectively, these data illustrate the heterogeneous HLA peptide presentation landscape even within one individual, and hint that a multi-peptide vaccination approach against highly conserved tumor suppressors may be a viable option in patients with low tumor-mutational burden.

Publisher Correction: The mutational impact of culturing human pluripotent and adult stem cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The mutational impact of culturing human pluripotent and adult stem cells.

Genetic changes acquired during in vitro culture pose a risk for the successful application of stem cells in regenerative medicine. To assess the genetic risks induced by culturing, we determined all mutations in individual human stem cells by whole genome sequencing. Individual pluripotent, intestinal, and liver stem cells accumulate 3.5 ± 0.5, 7.2 ± 1.1 and 8.3 ± 3.6 base substitutions per population doubling, respectively. The annual in vitro mutation accumulation rate of adult stem cells is nearly 40-fold higher than the in vivo mutation accumulation rate. Mutational signature analysis reveals that in vitro induced mutations are caused by oxidative stress. Reducing oxygen tension in culture lowers the mutational load. We use the mutation rates, spectra, and genomic distribution to model the accumulation of oncogenic mutations during typical in vitro expansion, manipulation or screening experiments using human stem cells. Our study provides empirically defined parameters to assess the mutational risk of stem cell based therapies.

Mutational signature in colorectal cancer caused by genotoxic pks+ E. coli.

Various species of the intestinal microbiota have been associated with the development of colorectal cancer1,2, but it has not been demonstrated that bacteria have a direct role in the occurrence of oncogenic mutations. Escherichia coli can carry the pathogenicity island pks, which encodes a set of enzymes that synthesize colibactin3. This compound is believed to alkylate DNA on adenine residues4,5 and induces double-strand breaks in cultured cells3. Here we expose human intestinal organoids to genotoxic pks+ E. coli by repeated luminal injection over five months. Whole-genome sequencing of clonal organoids before and after this exposure revealed a distinct mutational signature that was absent from organoids injected with isogenic pks-mutant bacteria. The same mutational signature was detected in a subset of 5,876 human cancer genomes from two independent cohorts, predominantly in colorectal cancer. Our study describes a distinct mutational signature in colorectal cancer and implies that the underlying mutational process results directly from past exposure to bacteria carrying the colibactin-producing pks pathogenicity island.

Pancreatic cancer organoids recapitulate disease and allow personalized drug screening.

We report the derivation of 30 patient-derived organoid lines (PDOs) from tumors arising in the pancreas and distal bile duct. PDOs recapitulate tumor histology and contain genetic alterations typical of pancreatic cancer. In vitro testing of a panel of 76 therapeutic agents revealed sensitivities currently not exploited in the clinic, and underscores the importance of personalized approaches for effective cancer treatment. The PRMT5 inhibitor EZP015556, shown to target MTAP (a gene commonly lost in pancreatic cancer)-negative tumors, was validated as such, but also appeared to constitute an effective therapy for a subset of MTAP-positive tumors. Taken together, the work presented here provides a platform to identify novel therapeutics to target pancreatic tumor cells using PDOs.

The genomic landscape of metastatic breast cancer highlights changes in mutation and signature frequencies.

The whole-genome sequencing of prospectively collected tissue biopsies from 442 patients with metastatic breast cancer reveals that, compared to primary breast cancer, tumor mutational burden doubles, the relative contributions of mutational signatures shift and the mutation frequency of six known driver genes increases in metastatic breast cancer. Significant associations with pretreatment are also observed. The contribution of mutational signature 17 is significantly enriched in patients pretreated with fluorouracil, taxanes, platinum and/or eribulin, whereas the de novo mutational signature I identified in this study is significantly associated with pretreatment containing platinum-based chemotherapy. Clinically relevant subgroups of tumors are identified, exhibiting either homologous recombination deficiency (13%), high tumor mutational burden (11%) or specific alterations (24%) linked to sensitivity to FDA-approved drugs. This study provides insights into the biology of metastatic breast cancer and identifies clinically useful genomic features for the future improvement of patient management.

5-Fluorouracil treatment induces characteristic T>G mutations in human cancer.

5-Fluorouracil (5-FU) is a chemotherapeutic drug commonly used for the treatment of solid cancers. It is proposed that 5-FU interferes with nucleotide synthesis and incorporates into DNA, which may have a mutational impact on both surviving tumor and healthy cells. Here, we treat intestinal organoids with 5-FU and find a highly characteristic mutational pattern that is dominated by T>G substitutions in a CTT context. Tumor whole genome sequencing data confirms that this signature is also identified in vivo in colorectal and breast cancer patients who have received 5-FU treatment. Taken together, our results demonstrate that 5-FU is mutagenic and may drive tumor evolution and increase the risk of secondary malignancies. Furthermore, the identified signature shows a strong resemblance to COSMIC signature 17, the hallmark signature of treatment-naive esophageal and gastric tumors, which indicates that distinct endogenous and exogenous triggers can converge onto highly similar mutational signatures.