The repertoire of mutational signatures in human cancer

Discovering the Mutational Profile of Early Colorectal Lesions: A Translational Impact

Simple Summary

Colorectal cancer (CRC) is one of the most common malignancies worldwide. Next-generation sequencing technologies have identified new candidate genes and deepened the knowledge of the molecular mechanisms underlying the progression of colonic adenomas towards CRC. The main genetic, epigenetic, and molecular alterations driving the onset and progression of CRC in both hereditary and sporadic settings have also been investigated. The evaluation of the CRC risk based on the molecular characterization of early pre-cancerous lesions may contribute to the development of targeted preventive strategies development, help define specific risk profiles, and identify patients who will benefit from targeted endoscopic surveillance.

Abstract

Colorectal cancer (CRC) develops through a multi-step process characterized by the acquisition of multiple somatic mutations in oncogenes and tumor-suppressor genes, epigenetic alterations and genomic instability. These events lead to the progression from precancerous lesions to advanced carcinomas. This process requires several years in a sporadic setting, while occurring at an early age and or faster in patients affected by hereditary CRC-predisposing syndromes. Since advanced CRC is largely untreatable or unresponsive to standard or targeted therapies, the endoscopic treatment of colonic lesions remains the most efficient CRC-preventive strategy. In this review, we discuss recent studies that have assessed the genetic alterations in early colorectal lesions in both hereditary and sporadic settings. Establishing the genetic profile of early colorectal lesions is a critical goal in the development of risk-based preventive strategies.

Genetic interactions among Brca1, Brca2, Palb2, and Trp53 in mammary tumor development

Inherited mutations in BRCA1, BRCA2, and PALB2 cause a high risk of breast cancer. Here, we conducted parallel conditional knockout (CKO) of Brca1, Palb2, and Brca2, individually and in combination, along with one copy of Trp53, in the mammary gland of nulliparous female mice. We observed a functional equivalence of the three genes in their basic tumor-suppressive activity, a linear epistasis of Palb2 and Brca2, but complementary roles of Brca1 and Palb2 in mammary tumor suppression, as combined ablation of either Palb2 or Brca2 with Brca1 led to delayed tumor formation. Whole-exome sequencing (WES) revealed both similarities and differences between Brca1 and Palb2 or Brca2 null tumors. Analyses of mouse mammary glands and cultured human cells showed that combined loss of BRCA1 and PALB2 led to high levels of reactive oxygen species (ROS) and increased apoptosis, implicating oxidative stress in the delayed tumor development in Brca1;Palb2 double CKO mice. The functional complementarity between BRCA1 and PALB2/BRCA2 and the role of ROS in tumorigenesis require further investigation.

The Quandary of DNA-Based Treatment Assessment in De Novo Metastatic Prostate Cancer in the Era of Precision Oncology

Guidelines for genetic testing have been established for multiple tumor types, frequently indicating the most confident molecularly targeted treatment options. However, considering the often-complex presentation of individual cancer patients, in addition to the combinatorial complexity and inherent uncertainties of molecular findings, deriving optimal treatment strategies frequently becomes very challenging. Here, we report a comprehensive analysis of a 68-year-old male with metastatic prostate cancer, encompassing pathology and MRI findings, transcriptomic results, and key genomics findings from whole-exome sequencing, both somatic aberrations and germline variants. We identify multiple somatic aberrations that are known to be enriched in prostate cancer, including a deletion of PTEN and a fusion transcript involving BRCA2. The gene expression patterns in the tumor biopsy were also strikingly similar to prostate tumor samples from TCGA. Furthermore, we detected multiple lines of evidence for homologous recombination repair deficiency (HRD), including a dominant contribution by mutational signature SBS3, which is specifically attributed to HRD. On the basis of the genomic and transcriptomic findings, and in light of the clinical case presentation, we discussed the personalized treatment options that exist for this patient and the various challenges that one faces in the process of translating high-throughput sequencing data towards treatment regimens.

Characteristics of Pan-Cancer Patients With Ultrahigh Tumor Mutation Burden

Background

Tumor mutation burden has been proven to be a good predictor for the efficacy of immunotherapy, especially in patients with hypermutation. However, most research focused on the analysis of hypermutation in individual tumors, and there is a lack of integrated research on the hypermutation across different cancers. This study aimed to characterize hypermutated patients to distinguish between these patients and non-hypermutated patients.

Methods

A total of 5,980 tumor samples involving 23 types of solid tumors from the in-house database were included in the study. Based on the cutoff value of tumor mutation burden (TMB), all samples were divided into hypermutated or non-hypermutated groups. Microsatellite instability status, PD-L1 expression and other mutation-related indicators were analyzed.

Results

Among the 5,980 tumor samples, 1,164 were selected as samples with hypermutation. Compared with the non-hypermutated group, a significant increase in the mutation rates of DNA mismatch repair genes and polymerase genes was detected in the hypermutated group, and there was an overlap between high TMB and high microsatellite instability or high PD-L1. In addition, we found that EGFR, KRAS and PIK3CA had a high frequency of both single nucleotide variation and copy number variation mutations. These identified mutant genes were enriched in the oncogenic signaling pathway and the DNA damage repair pathway. At the same time, the somatic cell characteristics and distribution of the two groups were significantly different.

Conclusions

This study identified genetic and phenotypic characteristics of hypermutated tumors and demonstrated that DNA damage repair is critically involved in hypermutation.

p53 Expression, Programmed Death Ligand 1, and Risk Factors in Urinary Tract Small Cell Carcinoma

Introduction: Small cell carcinoma of the urinary tract (SCCUT) is a rare finding with poor clinical course. This study sheds light on the molecular subtype and identifies risk factors in patients diagnosed with SCCUT.

Methods: Immunohistochemical expression of immunotherapy target programmed death ligand 1 (PD-L1) and luminal (GATA3), basal (p63), and p53 markers are assessed in patients diagnosed with SCCUT. Univariate analysis identified risk factors. Overall survival (OS) is computed using the Kaplan–Meier method.

Results: Tissue was available for 70.2% (33/47). All showed a high PD-L1 expression phenotype. p53 is seen in 93.9% (31/33), mostly as overexpression, GATA3 in 45.5% (15/33), and p63 in 57.6% (19/33). For the entire cohort (n = 47), 1-year survival was 59.6%, and the median OS was 17 months. Univariate analysis shows that chemotherapy [hazard ratio (HR) = 0.29, 95% confidence interval (CI) = 0.14–0.61, p = 0.001], radical surgery (HR = 0.37, 95% CI = 0.18–0.76, p = 0.007), and diagnosis of non-pure SCCUT (HR = 0.44, 95% CI = 0.22–0.86, p = 0.02) are favorable prognostic features. Metastasis had negative associations with survival (HR = 2.1, 95% CI = 1.1–4.2, p = 0.03).

Conclusions: In this series, pure and mixed SCCUT are characterized by p53 overexpression and a high PD-L1 phenotype. Histology of non-pure SCCUT is a positive prognosticator, and radical cystectomy or chemotherapy can improve OS. These findings demonstrate that SCCUT may be eligible for PD-L1 immunotherapy.

Genome instability and loss of protein homeostasis: converging paths to neurodegeneration?

Genome instability and loss of protein homeostasis are hallmark events of age-related diseases that include neurodegeneration. Several neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, Huntington's disease and amyotrophic lateral sclerosis are characterized by protein aggregation, while an impaired DNA damage response (DDR) as in many genetic DNA repair disorders leads to pronounced neuropathological features. It remains unclear to what degree these cellular events interconnect with each other in the development of neurological diseases. This review highlights how the loss of protein homeostasis and genome instability influence one other. We will discuss studies that illustrate this connection. DNA damage contributes to many neurodegenerative diseases, as shown by an increased level of DNA damage in patients, possibly due to the effects of protein aggregates on chromatin, the sequestration of DNA repair proteins and novel putative DNA repair functions. Conversely, genome stability is also important for protein homeostasis. For example, gene copy number variations and the loss of key DDR components can lead to marked proteotoxic stress. An improved understanding of how protein homeostasis and genome stability are mechanistically connected is needed and promises to lead to the development of novel therapeutic interventions.

Multidisciplinary management of locally advanced pancreatic adenocarcinoma: Biology is King

The annual incidence of pancreatic cancer is nearly 50,000 patients. The 5-year overall survival is only 9%, and there remains a great need for better therapy. A subset of these patients presents with locally advanced disease. Multidisciplinary therapy has evolved to include some combination of systemic chemotherapy, locoregional radiation, and surgery in select patients with excellent biology. This review will address the thoughtful evidence-based and individualized approach to these patients.

Transcription-coupled nucleotide excision repair: New insights revealed by genomic approaches

Elongation of RNA polymerase II (Pol II) is affected by many factors including DNA damage. Bulky damage, such as lesions caused by ultraviolet (UV) radiation, arrests Pol II and inhibits gene transcription, and may lead to genome instability and cell death. Cells activate transcription-coupled nucleotide excision repair (TC-NER) to remove Pol II-impeding damage and allow transcription resumption. TC-NER initiation in humans is mediated by Cockayne syndrome group B (CSB) protein, which binds to the stalled Pol II and promotes assembly of the repair machinery. Given the complex nature of the TC-NER pathway and its unique function at the interface between transcription and repair, new approaches are required to gain in-depth understanding of the mechanism. Advances in genomic approaches provide an important opportunity to investigate how TC-NER is initiated upon damage-induced Pol II stalling and what factors are involved in this process. In this Review, we discuss new mechanisms of TC-NER revealed by genome-wide DNA damage mapping and new TC-NER factors identified by high-throughput screening. As TC-NER conducts strand-specific repair of mutagenic damage, we also discuss how this repair pathway causes mutational strand asymmetry in the cancer genome.

TERT promoter hotspot mutations and gene amplification in metaplastic breast cancer

Metaplastic breast cancers (MBCs) are characterized by complex genomes, which seem to vary according to their histologic subtype. TERT promoter hotspot mutations and gene amplification are rare in common forms of breast cancer, but present in a subset of phyllodes tumors. Here, we sought to determine the frequency of genetic alterations affecting TERT in a cohort of 60 MBCs with distinct predominant metaplastic components (squamous, 23%; spindle, 27%; osseous, 8%; chondroid, 42%), and to compare the repertoire of genetic alterations of MBCs according to the presence of TERT promoter hotspot mutations or gene amplification. Forty-four MBCs were subjected to: whole-exome sequencing (WES; n = 27) or targeted sequencing of 341-468 cancer-related genes (n = 17); 16 MBCs were subjected to Sanger sequencing of the TERT promoter, TP53 and selected exons of PIK3CA, HRAS, and BRAF. TERT promoter hotspot mutations (n = 9) and TERT gene amplification (n = 1) were found in 10 of the 60 MBCs analyzed, respectively. These TERT alterations were less frequently found in MBCs with predominant chondroid differentiation than in other MBC subtypes (p = 0.01, Fisher's exact test) and were mutually exclusive with TP53 mutations (p < 0.001, CoMEt). In addition, a comparative analysis of the MBCs subjected to WES or targeted cancer gene sequencing (n = 44) revealed that MBCs harboring TERT promoter hotspot mutations or gene amplification (n = 6) more frequently harbored PIK3CA than TERT wild-type MBCs (n = 38; p = 0.001; Fisher's exact test). In conclusion, TERT somatic genetic alterations are found in a subset of TP53 wild-type MBCs with squamous/spindle differentiation, highlighting the genetic diversity of these cancers.

Aust D, B. Baretton G. Routine Molecular Pathology Diagnostics in Precision Oncology

BACKGROUND

Technical advances in the field of molecular genetics permit precise genomic characterization of malignant tumors. This has not only improved our understanding of tumor biology but also paved the way for molecularly stratified treatment strategies in routine clinical practice.

METHODS

A selective search of PubMed to identify literature on molecular pathology methods, their indications, the challenges associated with molecular findings, and future developments.

RESULTS

Tumors can be characterized with the aid of immunohistochemistry, in-situ hybridization, and sequencing of DNA or RNA. The benefits of molecularly stratified tumor treatment have been demonstrated by randomized clinical trials on numerous tumor entities, e.g., non-small-cell lung cancer, colorectal cancer, and breast cancer. Therefore, initiation of specific treatment for these entities should be preceded by molecular pathology biomarker analyses, generally carried out on tumor tissue. Randomized controlled trials and non-controlled studies show that enhanced progression-free survival ensues if the pharmacological treatment is oriented on the findings of molecular pathology diagnostics. In next-generation sequencing, numerous relevant gene sequences or even whole genes can be sequenced in parallel, dispensing with complex staged diagnostics and reducing the use of biomaterials. These new methods also complement the currently relevant predictive biomarkers by permitting the investigation of genetic alterations presently of interest in the context of clinical studies. Prior to widespread routine clinical application, however, sequencing of large gene panels or whole genomes or exomes need to be even more stringently validated.

CONCLUSION

Quality-assured molecular pathology assays are universally available for the determination of currently relevant predictive biomarkers. However, the integration of extensive genomic analyses into routine molecular pathology diagnostics represents a future challenge in precision oncology.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Predicting clinical outcomes using cancer progression associated signatures

Somatic mutation signatures are an informative facet of cancer aetiology, however they are rarely useful for predicting patient outcome. The aim of this study is to evaluate the utility of a panel of 142 mutation-signature–associated metrics (P142) for predicting cancer progression in patients from a ‘TCGA PanCancer Atlas’ cohort. The P142 metrics are comprised of AID/APOBEC and ADAR deaminase associated SNVs analyzed for codon context, strand bias, and transitions/transversions. TCGA tumor-normal mutation data was obtained for 10,437 patients, representing 31 of the most prevalent forms of cancer. Stratified random sampling was used to split patients into training, tuning and validation cohorts for each cancer type. Cancer specific machine learning (XGBoost) models were built using the output from the P142 panel to predict patient Progression Free Survival (PFS) status as either “High PFS” or “Low PFS”. Predictive performance of each model was evaluated using the validation cohort. Models accurately predicted PFS status for several cancer types, including adrenocortical carcinoma, glioma, mesothelioma, and sarcoma. In conclusion, the P142 panel of metrics successfully predicted cancer progression status in patients with some, but not all cancer types analyzed. These results pave the way for future studies on cancer progression associated signatures.

Moss enables high sensitivity single-nucleotide variant calling from multiple bulk DNA tumor samples

Intra-tumor heterogeneity renders the identification of somatic single-nucleotide variants (SNVs) a challenging problem. In particular, low-frequency SNVs are hard to distinguish from sequencing artifacts. While the increasing availability of multi-sample tumor DNA sequencing data holds the potential for more accurate variant calling, there is a lack of high-sensitivity multi-sample SNV callers that utilize these data. Here we report Moss, a method to identify low-frequency SNVs that recur in multiple sequencing samples from the same tumor. Moss provides any existing single-sample SNV caller the ability to support multiple samples with little additional time overhead. We demonstrate that Moss improves recall while maintaining high precision in a simulated dataset. On multi-sample hepatocellular carcinoma, acute myeloid leukemia and colorectal cancer datasets, Moss identifies new low-frequency variants that meet manual review criteria and are consistent with the tumor's mutational signature profile. In addition, Moss detects the presence of variants in more samples of the same tumor than reported by the single-sample caller. Moss' improved sensitivity in SNV calling will enable more detailed downstream analyses in cancer genomics.

Evolutionary Trajectories and Genomic Divergence in Localized Breast Cancers after Ipsilateral Breast Tumor Recurrence

Simple Summary

Ipsilateral breast tumor relapse (IBTR) occurs in 5–10% of localized breast cancers (BCs) within 10 years of incidence, despite proper treatment of the primary lesion. However, the clinical consequences of evolutionary trajectories of BC cells and their impact on IBTR remain poorly understood. Here, we conducted a longitudinal genomic analysis of 10 matched localized BC patients with IBTR. Overall, we identified the differences in homologous recombination deficiency, chromosomal instability, and somatic mutation drivers between primary and relapsed lesions. Our analyses highlighted three clonal architectures that shape by distinct mutagenic processes and subclonal diversification during relapse progression. Finally, this study provided a framework, which integrated actionable biomarkers with clonal architectures, towards improvement of future treatment decisions.

Abstract

The evolutionary trajectories that drive clinical and therapeutic consequences in localized breast cancers (BCs) with ipsilateral breast tumor relapse (IBTR) remain largely unknown. Analyses of longitudinal paired whole-exome sequencing data from 10 localized BC patients with IBTR reveal that, compared to primary breast tumors, homologous recombination (HR) deficiency, inactivation of the HR pathway, chromosomal instability, and somatic driver mutations are more frequent. Furthermore, three major models of evolution in IBTR are summarized, through which relative contributions of mutational signatures shift, and the subclonal diversity expansions are shown. Optimal treatment regimens are suggested by the clinically relevant molecular features, such as HR deficiency (20%) or specific alterations (30%) with sensitivity to available FDA-approved drugs. Finally, a rationale for the development of the therapeutic management framework is provided. This study sheds light on the complicated evolution patterns in IBTR and has significant clinical implications for future improvement of treatment decisions.

Genomic profiling of multiple tissues in two patients with multiple endocrine neoplasia type 1

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant tumor syndrome. This hereditary cancer is caused by germline variants in MEN1. Two patients with MEN1 were identified via whole exome sequencing and gene expression profile analysis, conducted for 5,063 patients with various types of cancers. We obtained multiple tumors from each patient; tumors derived from these two MEN1 patients had a loss of the normal MEN1 allele and frequently chromosomal copy number changes. Thus, we investigated whether structural variants were present in the MEN1 patient genomes. Whole-genome sequencing revealed no catastrophic rearrangements, and the tumor samples had very low somatic variants. The two patients had germline variants in MEN1 and some chromosomal copy number changes including on chromosome 11. The only pathogenic variant detected was the MEN1 germline variant, and chromosomal rearrangements led to tumorigenesis in somatic cells. Furthermore, the MEN1 tumor samples displayed a specific signature characterized by T:A>C:G transition. Studies of multiple tumors obtained from single patients are rare in hereditary cancer syndromes, and our results provide insights that the second hit of the tumor suppressor gene MEN1 may be caused by a gross genome rearrangement, not a small insertion and deletion, nor a change in epigenetic regulation.

Ppp6c haploinsufficiency accelerates UV-induced BRAF(V600E)-initiated melanomagenesis

According to TCGA database, mutations in PPP6C (encoding phosphatase PP6) are found in c. 10% of tumors from melanoma patients, in which they coexist with BRAF and NRAS mutations. To assess PP6 function in melanoma carcinogenesis, we generated mice in which we could specifically induce BRAF(V600E) expression and delete Ppp6c in melanocytes. In these mice, melanoma susceptibility following UVB irradiation exhibited the following pattern: Ppp6c semi‐deficient (heterozygous) > Ppp6c wild‐type > Ppp6c‐deficient (homozygous) tumor types. Next‐generation sequencing of Ppp6c heterozygous and wild‐type melanoma tumors revealed that all harbored Trp53 mutations. However, Ppp6c heterozygous tumors showed a higher Signature 1 (mitotic/mitotic clock) mutation index compared with Ppp6c wild‐type tumors, suggesting increased cell division. Analysis of cell lines derived from either Ppp6c heterozygous or wild‐type melanoma tissues showed that both formed tumors in nude mice, but Ppp6c heterozygous tumors grew faster compared with those from the wild‐type line. Ppp6c knockdown via siRNA in the Ppp6c heterozygous line promoted the accumulation of genomic damage and enhanced apoptosis relative to siRNA controls. We conclude that in the presence of BRAF(V600E) expression and UV‐induced Trp53 mutation, Ppp6c haploinsufficiency promotes tumorigenesis.

Integrated approaches for precision oncology in colorectal cancer: The more you know, the better

Colorectal cancer (CRC) is one of the most common human malignancies accounting for approximately 10 % of worldwide cancer incidence and mortality. While early-stage CRC is mainly a preventable and curable disease, metastatic colorectal cancer (mCRC) remains an unmet clinical need. Moreover, about 25 % of CRC cases are diagnosed only at the metastatic stage. Despite the extensive molecular and functional knowledge on this disease, systemic therapy for mCRC still relies on traditional 5-fluorouracil (5-FU)-based chemotherapy regimens. On the other hand, targeted therapies and immunotherapy have shown effectiveness only in a limited subset of patients. For these reasons, there is a growing need to define the molecular and biological landscape of individual patients to implement novel, rationally driven, tailored therapies. In this review, we explore current and emerging approaches for CRC management such as genomic, transcriptomic and metabolomic analysis, the use of liquid biopsies and the implementation of patients' preclinical avatars. In particular, we discuss the contribution of each of these tools in elucidating patient specific features, with the aim of improving our ability in advancing the diagnosis and treatment of colorectal tumors.

MutViz 2.0: visual analysis of somatic mutations and the impact of mutational signatures on selected genomic regions

Patterns of somatic single nucleotide variants observed in human cancers vary widely between different tumor types. They depend not only on the activity of diverse mutational processes, such as exposure to ultraviolet light and the deamination of methylated cytosines, but largely also on the sequence content of different genomic regions on which these processes act. With MutViz (http://gmql.eu/mutviz/), we have presented a user-friendly web tool for the identification of mutation enrichments that offers preloaded mutations from public datasets for a variety of cancer types, well organized within an effective database architecture. Somatic mutation patterns can be visually and statistically analyzed within arbitrary sets of small, user-provided genomic regions, such as promoters or collections of transcription factor binding sites. Here, we present MutViz 2.0, a largely extended and consolidated version of the tool: we took into account the immediate (trinucleotide) sequence context of mutations, improved the representation of clinical annotation of tumor samples and devised a method for signature refitting on limited genomic regions to infer the contribution of individual mutational processes to the mutation patterns observed in these regions. We described both the features of MutViz 2.0, concentrating on the novelties, and the substantial re-engineering of the cloud-based architecture.

Genomic landscape of extraordinary responses in metastatic breast cancer

Extreme responders to anticancer therapy are rare among advanced breast cancer patients. Researchers, however, have yet to investigate treatment responses therein on the whole exome level. We performed whole exome analysis to characterize the genomic landscape of extreme responders among metastatic breast cancer patients. Clinical samples were obtained from breast cancer patients who showed exceptional responses to anti-HER2 therapy or hormonal therapy and from those who did not. Matched breast tumor tissue (somatic DNA) and blood samples (germline DNA) were collected from a total of 30 responders and 15 non-responders. Whole exome sequencing using Illumina HiSeq2500 was performed for all 45 patients (90 samples). Somatic single nucleotide variants (SNVs), indels, and copy number variants (CNVs) were identified for the genomes of each patient. Group-specific somatic variants and mutational burden were statistically analyzed. Sequencing of cancer exomes for all patients revealed 1839 somatic SNVs (1661 missense, 120 nonsense, 43 splice-site, 15 start/stop-lost) and 368 insertions/deletions (273 frameshift, 95 in-frame), with a median of 0.7 mutations per megabase (range, 0.08 to 4.2 mutations per megabase). Responders harbored a significantly lower nonsynonymous mutational burden (median, 26 vs. 59, P = 0.02) and fewer CNVs (median 13.6 vs. 97.7, P = 0.05) than non-responders. Multivariate analyses of factors influencing progression-free survival showed that a high mutational burden and visceral metastases were significantly related with disease progression. Extreme responders to treatment for metastatic breast cancer are characterized by fewer nonsynonymous mutations and CNVs.

Folate Intake Alters Mutation Frequency and Profiles in a Tissue- and Dose-Specific Manner in MutaMouse Male Mice

BACKGROUND

While cancer is common, its incidence varies widely by tissue. These differences are attributable to variable risk factors, such as environmental exposure, genetic inheritance, and lifetime number of stem cell divisions in a tissue. Folate deficiency is generally associated with increased risk for colorectal cancer (CRC) and acute lymphocytic leukemia (ALL). Conversely, high folic acid (FA) intake has also been associated with higher CRC risk.

OBJECTIVE

Our objective was to compare the effect of folate intake on mutant frequency (MF) and types of mutations in the colon and bone marrow of mice.

METHODS

Five-week-old MutaMouse male mice were fed a deficient (0 mg FA/kg), control (2 mg FA/kg), or supplemented (8 mg FA/kg) diet for 20 wk. Tissue MF was assessed using the lacZ mutant assay and comparisons made by 2-factor ANOVA. LacZ mutant plaques were sequenced using next-generation sequencing, and diet-specific mutation profiles within each tissue were compared by Fisher's exact test.

RESULTS

In the colon, the MF was 1.5-fold and 1.3-fold higher in mice fed the supplemented diet compared with mice fed the control (P = 0.001) and deficient (P = 0.008) diets, respectively. This contrasted with the bone marrow MF in the same mice where the MF was 1.7-fold and 1.6-fold higher in mice fed the deficient diet compared with mice fed the control (P = 0.02) and supplemented (P = 0.03) diets, respectively. Mutation profiles and signatures (mutation context) were tissue-specific.

CONCLUSIONS

Our data indicate that dietary folate intake affects mutagenesis in a tissue- and dose-specific manner in mice. Mutation profiles were generally tissue- but not dose-specific, suggesting that altered cellular folate status appears to interact with endogenous mutagenic mechanisms in each tissue to create a permissive context in which specific mutation types accumulate. These data illuminate potential mechanisms underpinning differences in observed associations between folate intake/status and cancer.

Interpretation of allele-specific chromatin accessibility using cell state-aware deep learning

Genomic sequence variation within enhancers and promoters can have a significant impact on the cellular state and phenotype. However, sifting through the millions of candidate variants in a personal genome or a cancer genome, to identify those that impact cis-regulatory function, remains a major challenge. Interpretation of noncoding genome variation benefits from explainable artificial intelligence to predict and interpret the impact of a mutation on gene regulation. Here we generate phased whole genomes with matched chromatin accessibility, histone modifications, and gene expression for 10 melanoma cell lines. We find that training a specialized deep learning model, called DeepMEL2, on melanoma chromatin accessibility data can capture the various regulatory programs of the melanocytic and mesenchymal-like melanoma cell states. This model outperforms motif-based variant scoring, as well as more generic deep learning models. We detect hundreds to thousands of allele-specific chromatin accessibility variants (ASCAVs) in each melanoma genome, of which 15%-20% can be explained by gains or losses of transcription factor binding sites. A considerable fraction of ASCAVs are caused by changes in AP-1 binding, as confirmed by matched ChIP-seq data to identify allele-specific binding of JUN and FOSL1. Finally, by augmenting the DeepMEL2 model with ChIP-seq data for GABPA, the TERT promoter mutation, as well as additional ETS motif gains, can be identified with high confidence. In conclusion, we present a new integrative genomics approach and a deep learning model to identify and interpret functional enhancer mutations with allelic imbalance of chromatin accessibility and gene expression.

Integrated mutational landscape analysis of uterine leiomyosarcomas

Uterine leiomyosarcomas (uLMS) are aggressive tumors arising from the smooth muscle layer of the uterus. We analyzed 83 uLMS sample genetics, including 56 from Yale and 27 from The Cancer Genome Atlas (TCGA). Among them, a total of 55 Yale samples including two patient-derived xenografts (PDXs) and 27 TCGA samples have whole-exome sequencing (WES) data; 10 Yale and 27 TCGA samples have RNA-sequencing (RNA-Seq) data; and 11 Yale and 10 TCGA samples have whole-genome sequencing (WGS) data. We found recurrent somatic mutations in TP53, MED12, and PTEN genes. Top somatic mutated genes included TP53, ATRX, PTEN, and MEN1 genes. Somatic copy number variation (CNV) analysis identified 8 copy-number gains, including 5p15.33 (TERT), 8q24.21 (C-MYC), and 17p11.2 (MYOCD, MAP2K4) amplifications and 29 copy-number losses. Fusions involving tumor suppressors or oncogenes were deetected, with most fusions disrupting RB1, TP53, and ATRX/DAXX, and one fusion (ACTG2-ALK) being potentially targetable. WGS results demonstrated that 76% (16 of 21) of the samples harbored chromoplexy and/or chromothripsis. Clinically actionable mutational signatures of homologous-recombination DNA-repair deficiency (HRD) and microsatellite instability (MSI) were identified in 25% (12 of 48) and 2% (1 of 48) of fresh frozen uLMS, respectively. Finally, we found olaparib (PARPi; P = 0.002), GS-626510 (C-MYC/BETi; P < 0.000001 and P = 0.0005), and copanlisib (PIK3CAi; P = 0.0001) monotherapy to significantly inhibit uLMS-PDXs harboring derangements in C-MYC and PTEN/PIK3CA/AKT genes (LEY11) and/or HRD signatures (LEY16) compared to vehicle-treated mice. These findings define the genetic landscape of uLMS and suggest that a subset of uLMS may benefit from existing PARP-, PIK3CA-, and C-MYC/BET-targeted drugs.

Distinct patterns of somatic genomic alterations and mutational signatures in central and peripheral-type small-cell lung cancer

Background

Generally, small-cell lung cancer (SCLC) can be classified into central and peripheral-type. Genomic landscape and genetic heterogeneity between central and peripheral-type SCLCs is scarce.

Methods

We conducted whole-exome sequencing (WES) of 41 tumor/control pairs of SCLC samples. In all cases, 8,186 genes with non-synonymous mutations were identified, as well as significantly mutated genes (SMGs), tumor mutation burden (TMB), weighted genome instability index (wGII) and somatic copy number alteration (SCNA), the driver recurrent SCNA, and mutational signatures between central and peripheral-type SCLCs.

Results

TMB of peripheral-type SCLCs were higher than central-type. Smoking patients had significantly higher wGII and CNA burden than the non-smokers in SCLCs, these alterations were more obviously in central-types. Furthermore, the driver SCNA regions of central and peripheral-type SCLCs were different. Central and peripheral-type SCLCs had no common recurrent amplification (AMP) cytobands or genes, but had collective recurrent deletion (DEL) including four cytobands and five genes. Interestingly, the AMP 12q24.31 was negative correlated with overall survival (OS) in central but not peripheral-type SCLCs. Moreover, a de novo mutational signature A was found significantly higher in peripheral than central-type SCLCs. In parallel, signature D was predictive of poor outcome, which was mainly in peripheral SCLC patients. COSMIC signatures analysis revealed that the association with signature D and OS was similar to the positive relationship between signature 13 and outcome.

Conclusions

Central and peripheral-type SCLCs were different in immunotherapy response, genome instability, the driver SCNAs and mutational signatures, which leaded to differences of prognosis.

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

Identifying key questions in the ecology and evolution of cancer

The application of evolutionary and ecological principles to cancer prevention and treatment, as well as recognizing cancer as a selection force in nature, has gained impetus over the last 50 years. Following the initial theoretical approaches that combined knowledge from interdisciplinary fields, it became clear that using the eco-evolutionary framework is of key importance to understand cancer. We are now at a pivotal point where accumulating evidence starts to steer the future directions of the discipline and allows us to underpin the key challenges that remain to be addressed. Here, we aim to assess current advancements in the field and to suggest future directions for research. First, we summarize cancer research areas that, so far, have assimilated ecological and evolutionary principles into their approaches and illustrate their key importance. Then, we assembled 33 experts and identified 84 key questions, organized around nine major themes, to pave the foundations for research to come. We highlight the urgent need for broadening the portfolio of research directions to stimulate novel approaches at the interface of oncology and ecological and evolutionary sciences. We conclude that progressive and efficient cross-disciplinary collaborations that draw on the expertise of the fields of ecology, evolution and cancer are essential in order to efficiently address current and future questions about cancer.

Inherent mosaicism and extensive mutation of human placentas

Placentas can exhibit chromosomal aberrations that are absent from the fetus1. The basis of this genetic segregation, which is known as confined placental mosaicism, remains unknown. Here we investigated the phylogeny of human placental cells as reconstructed from somatic mutations, using whole-genome sequencing of 86 bulk placental samples (with a median weight of 28 mg) and of 106 microdissections of placental tissue. We found that every bulk placental sample represents a clonal expansion that is genetically distinct, and exhibits a genomic landscape akin to that of childhood cancer in terms of mutation burden and mutational imprints. To our knowledge, unlike any other healthy human tissue studied so far, the placental genomes often contained changes in copy number. We reconstructed phylogenetic relationships between tissues from the same pregnancy, which revealed that developmental bottlenecks genetically isolate placental tissues by separating trophectodermal lineages from lineages derived from the inner cell mass. Notably, there were some cases with full segregation-within a few cell divisions of the zygote-of placental lineages and lineages derived from the inner cell mass. Such early embryonic bottlenecks may enable the normalization of zygotic aneuploidy. We observed direct evidence for this in a case of mosaic trisomic rescue. Our findings reveal extensive mutagenesis in placental tissues and suggest that mosaicism is a typical feature of placental development.

Clonal expansion in non-cancer tissues

Cancer is a clonal disorder derived from a single ancestor cell and its progenies that are positively selected by acquisition of 'driver mutations'. However, the evolution of positively selected clones does not necessarily imply the presence of cancer. On the contrary, it has become clear that expansion of these clones in phenotypically normal or non-cancer tissues is commonly seen in association with ageing and/or in response to environmental insults and chronic inflammation. Recent studies have reported expansion of clones harbouring mutations in cancer driver genes in the blood, skin, oesophagus, bronchus, liver, endometrium and bladder, where the expansion could be so extensive that tissues undergo remodelling of an almost entire tissue. The presence of common cancer driver mutations in normal tissues suggests a strong link to cancer development, providing an opportunity to understand early carcinogenic processes. Nevertheless, some driver mutations are unique to normal tissues or have a mutation frequency that is much higher in normal tissue than in cancer, indicating that the respective clones may not necessarily be destined for evolution to cancer but even negatively selected for carcinogenesis depending on the mutated gene. Moreover, tissues that are remodelled by genetically altered clones might define functionalities of aged tissues or modified inflammatory processes. In this Review, we provide an overview of major findings on clonal expansion in phenotypically normal or non-cancer tissues and discuss their biological significance not only in cancer development but also in ageing and inflammatory diseases.

mmsig: a fitting approach to accurately identify somatic mutational signatures in hematological malignancies

Mutational signatures have emerged as powerful biomarkers in cancer patients, with prognostic and therapeutic implications. Wider clinical utility requires access to reproducible algorithms, which allow characterization of mutational signatures in a given tumor sample. Here, we show how mutational signature fitting can be applied to hematological cancer genomes to identify biologically and clinically important mutational processes, highlighting the importance of careful interpretation in light of biological knowledge. Our newly released R package mmsig comes with a dynamic error-suppression procedure that improves specificity in important clinical and biological settings. In particular, mmsig allows accurate detection of mutational signatures with low abundance, such as those introduced by APOBEC cytidine deaminases. This is particularly important in the most recent mutational signature reference (COSMIC v3.1) where each signature is more clearly defined. Our mutational signature fitting algorithm mmsig is a robust tool that can be implemented immediately in the clinic.

Rustad et al. present a software package for the R statistical environment for the accurately quantify of somatic mutational signatures in hematological malignancies

Unraveling the genomic landscape of colorectal cancer through mutational signatures

Colorectal cancer, along with most other cancer types, is driven by somatic mutations. Characteristic patterns of somatic mutations, known as mutational signatures, arise as a result of the activities of different mutational processes. Mutational signatures have diverse origins, including exogenous and endogenous sources. In the case of colorectal cancer, the analysis of mutational signatures has elucidated specific signatures for classically associated DNA repair deficiencies, namely mismatch repair (leading to microsatellite instability), base excision repair (due to MUTYH or NTHL1 mutations), and polymerase proofreading (due to POLE and POLD1 exonuclease domain mutations). Additional signatures also play a role in colorectal cancer, including those related to normal aging and those associated with gut microbiota, as well as a number of signatures with unknown etiologies. This chapter provides an overview of the current knowledge of mutational signatures, with a focus on colorectal cancer and on the recently reported signatures in physiologically normal and inflammatory bowel disease-affected somatic colon tissues.

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

Recent Advances in Our Knowledge of mCRC Tumor Biology and Genetics: A Focus on Targeted Therapy Development

Metastatic colorectal cancer (mCRC) remains a highly lethal malignancy although considerable progress has resulted from characterizing molecular alterations such as RAS mutation status and extent of microsatellite instability (MSI) to guide optimal use of available therapies. The availability of gene expression profiling, next generation sequencing technologies, proteomics analysis and other technologies provides high resolution information on individual tumors, including metastatic lesions to better define intra-tumor and inter-tumor heterogeneity. Recent literature applying this information to further customize personalized therapies is reviewed. Current biomarker-based stratification used to select optimal therapy that is personalized to the mutation profile of individual tumors is described. Recent literature using whole exome sequencing of metastatic lesions and primary CRC tumors and other advanced technologies to more fully elucidate the tumor biology specific to mCRC sub-types and to develop more precise therapies that improve outcomes is also reviewed.

Engaging the Innate and Adaptive Antitumor Immune Response in Lymphoma

Immunotherapy has emerged as a powerful therapeutic strategy for many malignancies, including lymphoma. As in solid tumors, early clinical trials have revealed that immunotherapy is not equally efficacious across all lymphoma subtypes. For example, immune checkpoint inhibition has a higher overall response rate and leads to more durable outcomes in Hodgkin lymphomas compared to non-Hodgkin lymphomas. These observations, combined with a growing understanding of tumor biology, have implicated the tumor microenvironment as a major determinant of treatment response and prognosis. Interactions between lymphoma cells and their microenvironment facilitate several mechanisms that impair the antitumor immune response, including loss of major histocompatibility complexes, expression of immunosuppressive ligands, secretion of immunosuppressive cytokines, and the recruitment, expansion, and skewing of suppressive cell populations. Accordingly, treatments to overcome these barriers are being rapidly developed and translated into clinical trials. This review will discuss the mechanisms of immune evasion, current avenues for optimizing the antitumor immune response, clinical successes and failures of lymphoma immunotherapy, and outstanding hurdles that remain to be addressed.

The Inherited and Familial Component of Early-Onset Colorectal Cancer

Early-onset colorectal cancer (EOCRC), defined as that diagnosed before the age of 50, accounts for 10–12% of all new colorectal cancer (CRC) diagnoses. Epidemiological data indicate that EOCRC incidence is increasing, despite the observed heterogeneity among countries. Although the cause for such increase remains obscure, ≈13% (range: 9–26%) of EOCRC patients carry pathogenic germline variants in known cancer predisposition genes, including 2.5% of patients with germline pathogenic variants in hereditary cancer genes traditionally not associated with CRC predisposition. Approximately 28% of EOCRC patients have family history of the disease. This article recapitulates current evidence on the inherited syndromes that predispose to EOCRC and its familial component. The evidence gathered support that all patients diagnosed with an EOCRC should be referred to a specialized genetic counseling service and offered somatic and germline pancancer multigene panel testing. The identification of a germline pathogenic variant in a known hereditary cancer gene has relevant implications for the clinical management of the patient and his/her relatives, and it may guide surgical and therapeutic decisions. The relative high prevalence of hereditary cancer syndromes and familial component among EOCRC patients supports further research that helps understand the genetic background, either monogenic or polygenic, behind this increasingly common disease.

The origins and genetic interactions of KRAS mutations are allele- and tissue-specific

Mutational activation of KRAS promotes the initiation and progression of cancers, especially in the colorectum, pancreas, lung, and blood plasma, with varying prevalence of specific activating missense mutations. Although epidemiological studies connect specific alleles to clinical outcomes, the mechanisms underlying the distinct clinical characteristics of mutant KRAS alleles are unclear. Here, we analyze 13,492 samples from these four tumor types to examine allele- and tissue-specific genetic properties associated with oncogenic KRAS mutations. The prevalence of known mutagenic mechanisms partially explains the observed spectrum of KRAS activating mutations. However, there are substantial differences between the observed and predicted frequencies for many alleles, suggesting that biological selection underlies the tissue-specific frequencies of mutant alleles. Consistent with experimental studies that have identified distinct signaling properties associated with each mutant form of KRAS, our genetic analysis reveals that each KRAS allele is associated with a distinct tissue-specific comutation network. Moreover, we identify tissue-specific genetic dependencies associated with specific mutant KRAS alleles. Overall, this analysis demonstrates that the genetic interactions of oncogenic KRAS mutations are allele- and tissue-specific, underscoring the complexity that drives their clinical consequences.

Clonal architecture in mesothelioma is prognostic and shapes the tumour microenvironment

Malignant Pleural Mesothelioma (MPM) is typically diagnosed 20-50 years after exposure to asbestos and evolves along an unknown evolutionary trajectory. To elucidate this path, we conducted multi-regional exome sequencing of 90 tumour samples from 22 MPMs acquired at surgery. Here we show that exomic intratumour heterogeneity varies widely across the cohort. Phylogenetic tree topology ranges from linear to highly branched, reflecting a steep gradient of genomic instability. Using transfer learning, we detect repeated evolution, resolving 5 clusters that are prognostic, with temporally ordered clonal drivers. BAP1/-3p21 and FBXW7/-chr4 events are always early clonal. In contrast, NF2/-22q events, leading to Hippo pathway inactivation are predominantly late clonal, positively selected, and when subclonal, exhibit parallel evolution indicating an evolutionary constraint. Very late somatic alteration of NF2/22q occurred in one patient 12 years after surgery. Clonal architecture and evolutionary clusters dictate MPM inflammation and immune evasion. These results reveal potentially drugable evolutionary bottlenecking in MPM, and an impact of clonal architecture on shaping the immune landscape, with potential to dictate the clinical response to immune checkpoint inhibition.

Development of a system combining comprehensive genotyping and organoid cultures for identifying and testing genotype-oriented personalised medicine for pancreatobiliary cancers

BACKGROUND

Pancreatobiliary cancer is a highly aggressive tumour with a dismal prognosis. Personalised medicine represents a promising and effective therapeutic approach for this intractable disease. In this study, we aimed to establish a system for identifying and testing genotype-oriented targeted drugs for pancreatobiliary cancers by combining exome sequencing and organoid culture of primary tumours.

METHODS

Tumour cells isolated from resected tumours were subjected to organoid cultures based on published protocols with modifications. Exome sequencing was performed on the primary tumours. Histopathological and molecular features of the primary tumours were validated in the corresponding organoids. Genotype-oriented candidate targeted drugs were identified from exome sequencing, and their efficacies were tested in the organoids.

RESULTS

Organoid cultures succeeded in 30 of 54 (55.6%) cases. Six primary cancers of the biliary tract and gall bladder were subjected to exome sequencing, which revealed a variety of somatic mutations of genes involved in signalling pathways, epigenetic modifiers, genome maintenance and metabolic enzymes. Most of the organoids of these 6 cases showed identical histopathological features and genomic aberrations as those of the primary tumours. Some of the aberrations were candidates for targeted therapies. Integrin-linked kinase (ILK) was one such candidate target, and an ILK inhibitor was confirmed to suppress proliferation of patient-derived organoids.

CONCLUSIONS

By combining exome sequencing and organoid culture, our model enabled to identify genotype-oriented targets for personalised medicine and to test efficacies of candidate targeted drugs in the organoids. The current proof-of-concept approach could increase therapeutic opportunities for patients with pancreatobiliary cancers.

Integrated Genomic and Transcriptomic Analysis reveals key genes for predicting dual-phenotype Hepatocellular Carcinoma Prognosis

Dual-phenotype hepatocellular carcinoma (DPHCC) expresses both hepatocyte and cholangiocyte markers, and is characterized by high recurrence and low survival rates. The underlying molecular mechanisms of DPHCC pathogenesis are unclear.

We performed whole exome sequencing and RNA sequencing of three subtypes of HCC (10 DPHCC, 10 CK19-positive HCC, and 14 CK19-negative HCC), followed by integrated bioinformatics analysis, including somatic mutation analysis, mutation signal analysis, differential gene expression analysis, and pathway enrichment analysis. Cox proportional hazard regression analyses were applied for exploring survival related characteristics.

We found that mutated genes in DPHCC patients were associated with carcinogenesis and immunity, and the up-regulated genes were mainly enriched in transcription-related and cancer-related pathways, and the down-regulated genes were mainly enriched in immune-related pathways. CXCL9 was selected as the hub gene, which is associated with immune cells and survival prognosis. Our results showed that low CXCL9 expression was significantly associated with poor prognosis, and its expression was significantly reduced in DPHCC samples.

In conclusion, we explored the molecular mechanisms governing DPHCC development and progression and identified CXCL9, which influences the immune microenvironment and prognosis of DPHCC and might be new clinically significant biomarkers for predicting prognosis.

Translational Research in the Era of Precision Medicine: Where We Are and Where We Will Go

The advent of Precision Medicine has globally revolutionized the approach of translational research suggesting a patient-centric vision with therapeutic choices driven by the identification of specific predictive biomarkers of response to avoid ineffective therapies and reduce adverse effects. The spread of "multi-omics" analysis and the use of sensors, together with the ability to acquire clinical, behavioral, and environmental information on a large scale, will allow the digitization of the state of health or disease of each person, and the creation of a global health management system capable of generating real-time knowledge and new opportunities for prevention and therapy in the individual person (high-definition medicine). Real world data-based translational applications represent a promising alternative to the traditional evidence-based medicine (EBM) approaches that are based on the use of randomized clinical trials to test the selected hypothesis. Multi-modality data integration is necessary for example in precision oncology where an Avatar interface allows several simulations in order to define the best therapeutic scheme for each cancer patient.

Mutational signatures associated with exposure to carcinogenic microplastic compounds bisphenol A and styrene oxide

Microplastic pollutants in oceans and food chains are concerning to public health. Common plasticizing compounds Bisphenol-A (BPA) and Styrene-7,8-Oxide (SO) are now labeled as carcinogens. We show that BPA and SO cause deoxyribonucleic acid damage and mutagenesis in human cells, and analyze the genome-wide point mutation and genomic rearrangement patterns associated with BPA and SO exposure. A subset of the single- and doublet base substitutions shows mutagenesis near or at guanine, consistent with these compounds' preferences to form guanosine adducts. Presence of other mutational signatures suggest additional mutagenesis probably due to complex effects of BPA and SO on diverse cellular processes. Analyzing data for 19 cancer cohorts, we find that tumors of digestive and urinary organs show relatively high similarity in mutational profiles, and the burden of such mutations increases with age. Even within the same cancer type, proportions of corresponding mutational patterns vary among the cohorts from different countries, as does the amount of microplastic waste in ocean waters. BPA and SO are relatively mild mutagens, and other environmental agents can also potentially generate similar, complex mutational patterns in cancer genomes. Nonetheless, our findings call for systematic evaluation of public health consequences of microplastic exposure worldwide.

High tumor mutation burden fails to predict immune checkpoint blockade response across all cancer types

BACKGROUND

High tumor mutation burden (TMB-H) has been proposed as a predictive biomarker for response to immune checkpoint blockade (ICB), largely due to the potential for tumor mutations to generate immunogenic neoantigens. Despite recent pan-cancer approval of ICB treatment for any TMB-H tumor, as assessed by the targeted FoundationOne CDx assay in nine tumor types, the utility of this biomarker has not been fully demonstrated across all cancers.

PATIENTS AND METHODS

Data from over 10 000 patient tumors included in The Cancer Genome Atlas were used to compare approaches to determine TMB and identify the correlation between predicted neoantigen load and CD8 T cells. Association of TMB with ICB treatment outcomes was analyzed by both objective response rates (ORRs, N = 1551) and overall survival (OS, N = 1936).

RESULTS

In cancer types where CD8 T-cell levels positively correlated with neoantigen load, such as melanoma, lung, and bladder cancers, TMB-H tumors exhibited a 39.8% ORR to ICB [95% confidence interval (CI) 34.9-44.8], which was significantly higher than that observed in low TMB (TMB-L) tumors [odds ratio (OR) = 4.1, 95% CI 2.9-5.8, P < 2 × 10^-16]. In cancer types that showed no relationship between CD8 T-cell levels and neoantigen load, such as breast cancer, prostate cancer, and glioma, TMB-H tumors failed to achieve a 20% ORR (ORR = 15.3%, 95% CI 9.2-23.4, P = 0.95), and exhibited a significantly lower ORR relative to TMB-L tumors (OR = 0.46, 95% CI 0.24-0.88, P = 0.02). Bulk ORRs were not significantly different between the two categories of tumors (P = 0.10) for patient cohorts assessed. Equivalent results were obtained by analyzing OS and by treating TMB as a continuous variable.

CONCLUSIONS

Our analysis failed to support application of TMB-H as a biomarker for treatment with ICB in all solid cancer types. Further tumor type-specific studies are warranted.

Genomic Instability Profiles at the Single Cell Level in Mouse Colorectal Cancers of Defined Genotypes

The genomes of many human CRCs have been sequenced, revealing a large number of genetic alterations. However, the molecular mechanisms underlying the accumulation of these alterations are still being debated. In this study, we examined colorectal tumours that developed in mice with Apc^lox/lox, LSL-Kras^G12D, and Tp53^lox/lox targetable alleles. Organoids were derived from single cells and the spectrum of mutations was determined by exome sequencing. The number of single nucleotide substitutions (SNSs) correlated with the age of the tumour, but was unaffected by the number of targeted cancer-driver genes. Thus, tumours that expressed mutant Apc, Kras, and Tp53 alleles had as many SNSs as tumours that expressed only mutant Apc. In contrast, the presence of large-scale (>10 Mb) copy number alterations (CNAs) correlated strongly with Tp53 inactivation. Comparison of the SNSs and CNAs present in organoids derived from the same tumour revealed intratumoural heterogeneity consistent with genomic lesions accumulating at significantly higher rates in tumour cells compared to normal cells. The rate of acquisition of SNSs increased from the early stages of cancer development, whereas large-scale CNAs accumulated later, after Tp53 inactivation. Thus, a significant fraction of the genomic instability present in cancer cells cannot be explained by aging processes occurring in normal cells before oncogenic transformation.

RRM2B Is Frequently Amplified Across Multiple Tumor Types: Implications for DNA Repair, Cellular Survival, and Cancer Therapy

RRM2B plays a crucial role in DNA replication, repair and oxidative stress. While germline RRM2B mutations have been implicated in mitochondrial disorders, its relevance to cancer has not been established. Here, using TCGA studies, we investigated RRM2B alterations in cancer. We found that RRM2B is highly amplified in multiple tumor types, particularly in MYC-amplified tumors, and is associated with increased RRM2B mRNA expression. We also observed that the chromosomal region 8q22.3–8q24, is amplified in multiple tumors, and includes RRM2B, MYC along with several other cancer-associated genes. An analysis of genes within this 8q-amplicon showed that cancers that have both RRM2B-amplified along with MYC have a distinct pattern of amplification compared to cancers that are unaltered or those that have amplifications in RRM2B or MYC only. Investigation of curated biological interactions revealed that gene products of the amplified 8q22.3–8q24 region have important roles in DNA repair, DNA damage response, oxygen sensing, and apoptosis pathways and interact functionally. Notably, RRM2B-amplified cancers are characterized by mutation signatures of defective DNA repair and oxidative stress, and at least RRM2B-amplified breast cancers are associated with poor clinical outcome. These data suggest alterations in RR2MB and possibly the interacting 8q-proteins could have a profound effect on regulatory pathways such as DNA repair and cellular survival, highlighting therapeutic opportunities in these cancers.

Comprehensive Genomic Characterization of Fifteen Early-Onset Lynch-Like Syndrome Colorectal Cancers

Simple Summary

The most prevalent type of hereditary colorectal cancer is called Lynch syndrome and it is characterized by a tumor phenotype called microsatellite instability (MSI). This disease is a consequence of germline (inheritable) variants in any of the four mismatch repair (MMR) DNA genes, being their identification essential to ensure their appropriate diagnosis and implementation of preventive measurements. Nevertheless, only 50% of patients with MSI and suspected Lynch syndrome actually carry a germline pathogenic variant in an MMR gene that explains the clinical entity. The remaining 50% are termed Lynch-like syndrome, and their causes remain unknown. In this work, we tried to elucidate the molecular mechanisms that underlie this rare entity in a group of early-onset Lynch-like syndrome colorectal cancer, through whole-exome sequencing of germline and tumor samples. We observed that one-third of these patients have somatic alterations in genes associated with the MMR system and that these could be the mechanism causing their unexplained MSI. Furthermore, we found that patients who showed biallelic somatic alterations also carried germline variants in new candidate genes associated with DNA repair functions and that this could be, partly, the cause of the early onset in this cohort.

Abstract

Lynch-like syndrome (LLS) is an increasingly common clinical challenge with an underlying molecular basis mostly unknown. To shed light onto it, we focused on a very young LLS early-onset colorectal cancer (CRC) cohort (diagnosis ≤ 40 y.o.), performing germline and tumor whole-exome sequencing (WES) of 15 patients, and additionally analyzing their corresponding tumor mutational burden (TMB) and mutational signatures. We identified four cases (27%) with double somatic putative variants in mismatch repair (MMR) core genes, as well as three additional cases (20%) with double MSH3 somatic alterations in tumors with unexplained MSH2/MSH6 loss of expression, and two cases (13%) with POLD1 potential biallelic alterations. Average TMB was significantly higher for LLS cases with double somatic alterations. Lastly, nine predicted deleterious variants in genes involved in the DNA repair functions and/or previously associated with CRC were found in nine probands, four of which also showed MMR biallelic somatic inactivation. In conclusion, we contribute new insights into LLS CRC, postulating MSH3 and POLD1 double somatic alterations as an underlying cause of a microsatellite instability (MSI) phenotype, proposing intrinsic biological differences between LLS with and without somatic alterations, and suggesting new predisposing candidate genes in this scenario.

APOBEC3A catalyzes mutation and drives carcinogenesis in vivo

The APOBEC3 family of antiviral DNA cytosine deaminases is implicated as the second largest source of mutation in cancer. This mutational process may be a causal driver or inconsequential passenger to the overall tumor phenotype. We show that human APOBEC3A expression in murine colon and liver tissues increases tumorigenesis. All other APOBEC3 family members, including APOBEC3B, fail to promote liver tumor formation. Tumor DNA sequences from APOBEC3A-expressing animals display hallmark APOBEC signature mutations in TCA/T motifs. Bioinformatic comparisons of the observed APOBEC3A mutation signature in murine tumors, previously reported APOBEC3A and APOBEC3B mutation signatures in yeast, and reanalyzed APOBEC mutation signatures in human tumor datasets support cause-and-effect relationships for APOBEC3A-catalyzed deamination and mutagenesis in driving multiple human cancers.

Meta-Analysis and Systematic Review of the Genomics of Mucosal Melanoma

Mucosal melanoma is a rare subtype of melanoma. To date, there has been no comprehensive systematic collation and statistical analysis of the aberrations and aggregated frequency of driver events across multiple studies. Published studies using whole genome, whole exome, targeted gene panel, or individual gene sequencing were identified. Datasets from these studies were collated to summarize mutations, structural variants, and regions of copy-number alteration. Studies using next-generation sequencing were divided into the "main" cohort (n = 173; fresh-frozen samples), "validation" cohort (n = 48; formalin-fixed, paraffin-embedded samples) and a second "validation" cohort comprised 104 tumors sequenced using a targeted panel. Studies assessing mutations in BRAF, KIT, and NRAS were summarized to assess hotspot mutations. Statistical analysis of the main cohort variant data revealed KIT, NF1, BRAF, NRAS, SF3B1, and SPRED1 as significantly mutated genes. ATRX and SF3B1 mutations occurred more commonly in lower anatomy melanomas and CTNNB1 in the upper anatomy. NF1, PTEN, CDKN2A, SPRED1, ATM, CHEK2, and ARID1B were commonly affected by chromosomal copy loss, while TERT, KIT, BRAF, YAP1, CDK4, CCND1, GAB2, MDM2, SKP2, and MITF were commonly amplified. Further notable genomic alterations occurring at lower frequencies indicated commonality of signaling networks in tumorigenesis, including MAPK, PI3K, Notch, Wnt/β-catenin, cell cycle, DNA repair, and telomere maintenance pathways. This analysis identified genomic aberrations that provide some insight to the way in which specific pathways may be disrupted. IMPLICATIONS: Our analysis has shown that mucosal melanomas have a diverse range of genomic alterations in several biological pathways. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/19/6/991/F1.large.jpg.

Comparison of Structural and Short Variants Detected by Linked-Read and Whole-Exome Sequencing in Multiple Myeloma

Linked-read sequencing was developed to aid the detection of large structural variants (SVs) from short-read sequencing efforts. We performed a systematic evaluation to determine if linked-read exome sequencing provides more comprehensive and clinically relevant information than whole-exome sequencing (WES) when applied to the same set of multiple myeloma patient samples. We report that linked-read sequencing detected a higher number of SVs (n = 18,455) than WES (n = 4065). However, linked-read predictions were dominated by inversions (92.4%), leading to poor detection of other types of SVs. In contrast, WES detected 56.3% deletions, 32.6% insertions, 6.7% translocations, 3.3% duplications and 1.2% inversions. Surprisingly, the quantitative performance assessment suggested a higher performance for WES (AUC = 0.791) compared to linked-read sequencing (AUC = 0.766) for detecting clinically validated cytogenetic alterations. We also found that linked-read sequencing detected more short variants (n = 704) compared to WES (n = 109). WES detected somatic mutations in all MM-related genes while linked-read sequencing failed to detect certain mutations. The comparison of somatic mutations detected using linked-read, WES and RNA-seq revealed that WES and RNA-seq detected more mutations than linked-read sequencing. These data indicate that WES outperforms and is more efficient than linked-read sequencing for detecting clinically relevant SVs and MM-specific short variants.

Immune-Related Mutational Landscape and Gene Signatures: Prognostic Value and Therapeutic Impact for Head and Neck Cancer

Simple Summary

Immunotherapy has emerged as a standard-of-care for most human malignancies, including head and neck cancer, but only a limited number of patients exhibit a durable clinical benefit. An urgent medical need is the establishment of accurate response predictors, which is handicapped by the growing body of molecular, cellular and clinical variables that modify the complex nature of an effective anti-tumor immune response. This review summarizes more recent efforts to elucidate immune-related mutational landscapes and gene expression signatures by integrative analysis of multi-omics data, and highlights their potential therapeutic impact for head and neck cancer. A better knowledge of the underlying principles and relevant interactions could pave the way for rational therapeutic combinations to improve the efficacy of immunotherapy, in particular for those cancer patients at a higher risk for treatment failure.

Abstract

Immunotherapy by immune checkpoint inhibition has become a main pillar in the armamentarium to treat head and neck cancer and is based on the premise that the host immune system can be reactivated to successfully eliminate cancer cells. However, the response rate remains low and only a small subset of head and neck cancer patients achieves a durable clinical benefit. The availability of multi-omics data and emerging computational technologies facilitate not only a deeper understanding of the cellular composition in the tumor immune microenvironment but also enables the study of molecular principles in the complex regulation of immune surveillance versus tolerance. These knowledges will pave the way to apply immunotherapy more precisely and effectively. This review aims to provide a holistic view on how the immune landscape dictates the tumor fate and vice versa, and how integrative analysis of multi-omics data contribute to our current knowledge on the accuracy of predictive biomarkers and on a broad range of factors influencing the response to immunotherapy in head and neck cancer.

Transitioning the Molecular Tumor Board from Proof of Concept to Clinical Routine: A German Single-Center Analysis

Molecular precision oncology faces two major challenges: first, to identify relevant and actionable molecular variants in a rapidly changing field and second, to provide access to a broad patient population. Here, we report a four-year experience of the Molecular Tumor Board (MTB) of the Comprehensive Cancer Center Freiburg (Germany) including workflows and process optimizations. This retrospective single-center study includes data on 488 patients enrolled in the MTB from February 2015 through December 2018. Recommendations include individual molecular diagnostics, molecular stratified therapies, assessment of treatment adherence and patient outcomes including overall survival. The majority of MTB patients presented with stage IV oncologic malignancies (90.6%) and underwent an average of 2.1 previous lines of therapy. Individual diagnostic recommendations were given to 487 patients (99.8%). A treatment recommendation was given in 264 of all cases (54.1%) which included a molecularly matched treatment in 212 patients (43.4%). The 264 treatment recommendations were implemented in 76 patients (28.8%). Stable disease was observed in 19 patients (25.0%), 17 had partial response (22.4%) and five showed a complete remission (6.6%). An objective response was achieved in 28.9% of cases with implemented recommendations and for 4.5% of the total population (22 of 488 patients). By optimizing the MTB workflow, case-discussions per session increased significantly while treatment adherence and outcome remained stable over time. Our data demonstrate the feasibility and effectiveness of molecular-guided personalized therapy for cancer patients in a clinical routine setting showing a low but robust and durable disease control rate over time.

Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling

Although melanoma is initiated by acquisition of point mutations and limited focal copy number alterations in melanocytes-of-origin, the nature of genetic changes that characterise lethal metastatic disease is poorly understood. Here, we analyze the evolution of human melanoma progressing from early to late disease in 13 patients by sampling their tumours at multiple sites and times. Whole exome and genome sequencing data from 88 tumour samples reveals only limited gain of point mutations generally, with net mutational loss in some metastases. In contrast, melanoma evolution is dominated by whole genome doubling and large-scale aneuploidy, in which widespread loss of heterozygosity sculpts the burden of point mutations, neoantigens and structural variants even in treatment-naïve and primary cutaneous melanomas in some patients. These results imply that dysregulation of genomic integrity is a key driver of selective clonal advantage during melanoma progression.

Reconstructing the Lineage Histories and Differentiation Trajectories of Individual Cancer Cells in Myeloproliferative Neoplasms

Some cancers originate from a single mutation event in a single cell. Blood cancers known as myeloproliferative neoplasms (MPNs) are thought to originate when a driver mutation is acquired by a hematopoietic stem cell (HSC). However, when the mutation first occurs in individuals and how it affects the behavior of HSCs in their native context is not known. Here we quantified the effect of the JAK2-V617F mutation on the self-renewal and differentiation dynamics of HSCs in treatment-naive individuals with MPNs and reconstructed lineage histories of individual HSCs using somatic mutation patterns. We found that JAK2-V617F mutations occurred in a single HSC several decades before MPN diagnosis-at age 9 ± 2 years in a 34-year-old individual and at age 19 ± 3 years in a 63-year-old individual-and found that mutant HSCs have a selective advantage in both individuals. These results highlight the potential of harnessing somatic mutations to reconstruct cancer lineages.