Loss of the DNA Repair Gene RNase H2 Identifies a Unique Subset of DDR-Deficient Leiomyosarcomas.

Targeting the DNA damage response (DDR) pathway is an emerging therapeutic approach for leiomyosarcoma (LMS), and loss of RNase H2, a DDR pathway member, is a potentially actionable alteration for DDR targeted treatments. Therefore, we designed a protein and genomic based RNase H2 screening assay to determine its prevalence and prognostic significance. Using a selective RNase H2 antibody on a pan-tumor tissue microarray (TMA), RNase H2 loss was more common in LMS (11.5%, 9/78) than across all tumors (3.8%, 32/843). In a separate LMS cohort, RNase H2 deficiency was confirmed in uterine LMS (U-LMS, 21%, 23/108) and soft-tissue LMS (ST-LMS) (30%, 39/102). In the TCGA database, RNASEH2B homozygous deletions (HomDels) were found in 6% (5/80) of LMS cases, with a higher proportion in U-LMS (15%; 4/27) compared to ST-LMS (2%; 1/53). Using the SNiPDx targeted-NGS sequencing assay to detect biallelic loss of function in select DDR related genes, we found RNASEH2B HomDels in 54% (19/35) of U-LMS cases with RNase H2 loss by IHC, and 7% (3/43) HomDels in RNase H2 intact cases. No RNASEH2B HomDels were detected in ST-LMS. In U-LMS patient cohort (n = 109), no significant overall survival difference was seen in patients with RNase H2 loss versus intact, or RNASEH2B HomDel (n=12) vs Non-HomDel (n=37). The overall diagnostic accuracy, sensitivity, and specificity of RNase H2 IHC for detecting RNASEH2B HomDels in U-LMS was 76%, 93% and 71% respectively, and it is being developed for future predictive biomarker driven clinical trials targeting DDR in U-LMS.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma.

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease. SIGNIFICANCE: The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Detection of Biallelic Loss of DNA Repair Genes in Formalin-Fixed, Paraffin-Embedded Tumor Samples Using a Novel Tumor-Only Sequencing Panel.

Patient selection for synthetic lethal-based cancer therapy may be improved by assessment of gene-specific loss of heterozygosity (LOH) and biallelic loss of function (LOF). This report describes SyNthetic lethal Interactions for Precision Diagnostics (SNiPDx), a targeted next-generation sequencing (NGS) panel for detection of LOH and biallelic LOF alterations in 26 target genes focused on DNA damage response pathways, in tumor-only formalin-fixed, paraffin-embedded (FFPE) samples. NGS was performed across all exons of these 26 genes and encompassed a total of 7632 genome-wide single-nucleotide polymorphisms on genomic DNA from 80 FFPE solid tumor samples. The Fraction and Allele-Specific Copy Number Estimates from Tumor Sequencing algorithm was optimized to assess tumor purity and copy number based on heterozygous single-nucleotide polymorphisms. SNiPDx demonstrated high sensitivity (95%) and specificity (91%) for LOH detection compared with whole genome sequencing. Positive agreement with local NGS-based testing in the detection of genetic alterations was 95%. SNiPDx detected 93% of biallelic ATM LOF mutations, 100% of ATM single-nucleotide variants and small insertions/deletions, and 100% of all ATM LOH status events identified by orthogonal NGS-based testing. SNiPDx is a novel, clinically feasible test for analysis of allelic status in FFPE tumor samples, which demonstrated high accuracy when compared with other NGS-based approaches in clinical use.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma.

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Functional patient-derived organoid screenings identify MCLA-158 as a therapeutic EGFR × LGR5 bispecific antibody with efficacy in epithelial tumors.

Patient-derived organoids (PDOs) recapitulate tumor architecture, contain cancer stem cells and have predictive value supporting personalized medicine. Here we describe a large-scale functional screen of dual-targeting bispecific antibodies (bAbs) on a heterogeneous colorectal cancer PDO biobank and paired healthy colonic mucosa samples. More than 500 therapeutic bAbs generated against Wingless-related integration site (WNT) and receptor tyrosine kinase (RTK) targets were functionally evaluated by high-content imaging to capture the complexity of PDO responses. Our drug discovery strategy resulted in the generation of MCLA-158, a bAb that specifically triggers epidermal growth factor receptor degradation in leucine-rich repeat-containing G-protein-coupled receptor 5-positive (LGR5+) cancer stem cells but shows minimal toxicity toward healthy LGR5+ colon stem cells. MCLA-158 exhibits therapeutic properties such as growth inhibition of KRAS-mutant colorectal cancers, blockade of metastasis initiation and suppression of tumor outgrowth in preclinical models for several epithelial cancer types.

Cancer therapy shapes the fitness landscape of clonal hematopoiesis.

Acquired mutations are pervasive across normal tissues. However, understanding of the processes that drive transformation of certain clones to cancer is limited. Here we study this phenomenon in the context of clonal hematopoiesis (CH) and the development of therapy-related myeloid neoplasms (tMNs). We find that mutations are selected differentially based on exposures. Mutations in ASXL1 are enriched in current or former smokers, whereas cancer therapy with radiation, platinum and topoisomerase II inhibitors preferentially selects for mutations in DNA damage response genes (TP53, PPM1D, CHEK2). Sequential sampling provides definitive evidence that DNA damage response clones outcompete other clones when exposed to certain therapies. Among cases in which CH was previously detected, the CH mutation was present at tMN diagnosis. We identify the molecular characteristics of CH that increase risk of tMN. The increasing implementation of clinical sequencing at diagnosis provides an opportunity to identify patients at risk of tMN for prevention strategies.

Comprehensive molecular comparison of BRCA1 hypermethylated and BRCA1 mutated triple negative breast cancers.

Homologous recombination deficiency (HRD) is a defining characteristic in BRCA-deficient breast tumors caused by genetic or epigenetic alterations in key pathway genes. We investigated the frequency of BRCA1 promoter hypermethylation in 237 triple-negative breast cancers (TNBCs) from a population-based study using reported whole genome and RNA sequencing data, complemented with analyses of genetic, epigenetic, transcriptomic and immune infiltration phenotypes. We demonstrate that BRCA1 promoter hypermethylation is twice as frequent as BRCA1 pathogenic variants in early-stage TNBC and that hypermethylated and mutated cases have similarly improved prognosis after adjuvant chemotherapy. BRCA1 hypermethylation confers an HRD, immune cell type, genome-wide DNA methylation, and transcriptional phenotype similar to TNBC tumors with BRCA1-inactivating variants, and it can be observed in matched peripheral blood of patients with tumor hypermethylation. Hypermethylation may be an early event in tumor development that progress along a common pathway with BRCA1-mutated disease, representing a promising DNA-based biomarker for early-stage TNBC.

11p15.5 epimutations in children with Wilms tumor and hepatoblastoma detected in peripheral blood.

BACKGROUND: Constitutional or somatic mosaic epimutations are increasingly recognized as a mechanism of gene dysregulation resulting in cancer susceptibility. Beckwith-Wiedemann syndrome is the cancer predisposition syndrome most commonly associated with epimutation and is extremely variable in its phenotypic presentation, which can include isolated tumors. Because to the authors' knowledge large-scale germline DNA sequencing studies have not included methylation analysis, the percentage of pediatric cancer predisposition that is due to epimutations is unknown. METHODS: Germline methylation testing at the 11p15.5 locus was performed in blood for 24 consecutive patients presenting with hepatoblastoma (3 patients) or Wilms tumor (21 patients). RESULTS: Six individuals with Wilms tumor and 1 patient with hepatoblastoma were found to have low-level gain of methylation at imprinting control 1, and a child with hepatoblastoma was found to have loss of methylation at imprinting control 2. The loss of methylation at imprinting control 2 was found to be maternally inherited, despite not being associated with any detectable genomic alteration. CONCLUSIONS: Overall, 33% of patients (8 of 24 patients) with Wilms tumor or hepatoblastoma were found to have an epigenetic susceptibility that was detectable in the blood. It is interesting to note that low-level gain of methylation at imprinting control 1 predominantly was detected in females with bilateral Wilms tumors. Further studies in larger cohorts are needed to determine the efficacy of testing all patients with Wilms tumor or hepatoblastoma for 11p15.5 epimutations in the blood as part of DNA analysis because this hallmark of predisposition will not be detected by sequencing-based approaches and detecting a cancer predisposition may modify treatment.

A practical framework and online tool for mutational signature analyses show inter-tissue variation and driver dependencies.

Mutational signatures are patterns of mutations that arise during tumorigenesis. We present an enhanced, practical framework for mutational signature analyses. Applying these methods on 3,107 whole genome sequenced (WGS) primary cancers of 21 organs reveals known signatures and nine previously undescribed rearrangement signatures. We highlight inter-organ variability of signatures and present a way of visualizing that diversity, reinforcing our findings in an independent analysis of 3,096 WGS metastatic cancers. Signatures with a high level of genomic instability are dependent on TP53 dysregulation. We illustrate how uncertainty in mutational signature identification and assignment to samples affects tumor classification, reinforcing that using multiple orthogonal mutational signature data is not only beneficial, it is essential for accurate tumor stratification. Finally, we present a reference web-based tool for cancer and experimentally-generated mutational signatures, called Signal (https://signal.mutationalsignatures.com), that also supports performing mutational signature analyses.

Revealing the impact of structural variants in multiple myeloma.

The landscape of structural variants (SVs) in multiple myeloma remains poorly understood. Here, we performed comprehensive analysis of SVs in a large cohort of 752 multiple myeloma patients by low coverage long-insert whole genome sequencing. We identified 68 SV hotspots involving 17 new candidate driver genes, including the therapeutic targets BCMA (TNFRSF17), SLAMF and MCL1. Catastrophic complex rearrangements termed chromothripsis were present in 24% of patients and independently associated with poor clinical outcomes. Templated insertions were the second most frequent complex event (19%), mostly involved in super-enhancer hijacking and activation of oncogenes such as CCND1 and MYC. Importantly, in 31% of patients two or more seemingly independent putative driver events were caused by a single structural event, demonstrating that the complex genomic landscape of multiple myeloma can be acquired through few key events during tumor evolutionary history. Overall, this study reveals the critical role of SVs in multiple myeloma pathogenesis.

Whole-genome sequencing of triple-negative breast cancers in a population-based clinical study.

Whole-genome sequencing (WGS) brings comprehensive insights to cancer genome interpretation. To explore the clinical value of WGS, we sequenced 254 triple-negative breast cancers (TNBCs) for which associated treatment and outcome data were collected between 2010 and 2015 via the population-based Sweden Cancerome Analysis Network-Breast (SCAN-B) project (ClinicalTrials.gov ID:NCT02306096). Applying the HRDetect mutational-signature-based algorithm to classify tumors, 59% were predicted to have homologous-recombination-repair deficiency (HRDetect-high): 67% explained by germline/somatic mutations of BRCA1/BRCA2, BRCA1 promoter hypermethylation, RAD51C hypermethylation or biallelic loss of PALB2. A novel mechanism of BRCA1 abrogation was discovered via germline SINE-VNTR-Alu retrotransposition. HRDetect provided independent prognostic information, with HRDetect-high patients having better outcome on adjuvant chemotherapy for invasive disease-free survival (hazard ratio (HR) = 0.42; 95% confidence interval (CI) = 0.2-0.87) and distant relapse-free interval (HR = 0.31, CI = 0.13-0.76) compared to HRDetect-low, regardless of whether a genetic/epigenetic cause was identified. HRDetect-intermediate, some possessing potentially targetable biological abnormalities, had the poorest outcomes. HRDetect-low cancers also had inadequate outcomes: ~4.7% were mismatch-repair-deficient (another targetable defect, not typically sought) and they were enriched for (but not restricted to) PIK3CA/AKT1 pathway abnormalities. New treatment options need to be considered for now-discernible HRDetect-intermediate and HRDetect-low categories. This population-based study advocates for WGS of TNBC to better inform trial stratification and improve clinical decision-making.

Genomic landscape and chronological reconstruction of driver events in multiple myeloma.

The multiple myeloma (MM) genome is heterogeneous and evolves through preclinical and post-diagnosis phases. Here we report a catalog and hierarchy of driver lesions using sequences from 67 MM genomes serially collected from 30 patients together with public exome datasets. Bayesian clustering defines at least 7 genomic subgroups with distinct sets of co-operating events. Focusing on whole genome sequencing data, complex structural events emerge as major drivers, including chromothripsis and a novel replication-based mechanism of templated insertions, which typically occur early. Hyperdiploidy also occurs early, with individual trisomies often acquired in different chronological windows during evolution, and with a preferred order of acquisition. Conversely, positively selected point mutations, whole genome duplication and chromoplexy events occur in later disease phases. Thus, initiating driver events, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of evolution that are biologically relevant but heterogeneous across patients.

Partially methylated domains are hypervariable in breast cancer and fuel widespread CpG island hypermethylation.

Global loss of DNA methylation and CpG island (CGI) hypermethylation are key epigenomic aberrations in cancer. Global loss manifests itself in partially methylated domains (PMDs) which extend up to megabases. However, the distribution of PMDs within and between tumor types, and their effects on key functional genomic elements including CGIs are poorly defined. We comprehensively show that loss of methylation in PMDs occurs in a large fraction of the genome and represents the prime source of DNA methylation variation. PMDs are hypervariable in methylation level, size and distribution, and display elevated mutation rates. They impose intermediate DNA methylation levels incognizant of functional genomic elements including CGIs, underpinning a CGI methylator phenotype (CIMP). Repression effects on tumor suppressor genes are negligible as they are generally excluded from PMDs. The genomic distribution of PMDs reports tissue-of-origin and may represent tissue-specific silent regions which tolerate instability at the epigenetic, transcriptomic and genetic level.

Author Correction: Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

Corrigendum: A somatic-mutational process recurrently duplicates germline susceptibility loci and tissue-specific super-enhancers in breast cancers.

This corrects the article DOI: 10.1038/ng.3771.

Short inverted repeats contribute to localized mutability in human somatic cells.

Selected repetitive sequences termed short inverted repeats (SIRs) have the propensity to form secondary DNA structures called hairpins. SIRs comprise palindromic arm sequences separated by short spacer sequences that form the hairpin stem and loop respectively. Here, we show that SIRs confer an increase in localized mutability in breast cancer, which is domain-dependent with the greatest mutability observed within spacer sequences (∼1.35-fold above background). Mutability is influenced by factors that increase the likelihood of formation of hairpins such as loop lengths (of 4-5 bp) and stem lengths (of 7-15 bp). Increased mutability is an intrinsic property of SIRs as evidenced by how almost all mutational processes demonstrate a higher rate of mutagenesis of spacer sequences. We further identified 88 spacer sequences showing enrichment from 1.8- to 90-fold of local mutability distributed across 283 sites in the genome that intriguingly, can be used to inform the biological status of a tumor.

Whole-Genome Sequencing Reveals Breast Cancers with Mismatch Repair Deficiency.

Mismatch repair (MMR)-deficient cancers have been discovered to be highly responsive to immune therapies such as PD-1 checkpoint blockade, making their definition in patients, where they may be relatively rare, paramount for treatment decisions. In this study, we utilized patterns of mutagenesis known as mutational signatures, which are imprints of the mutagenic processes associated with MMR deficiency, to identify MMR-deficient breast tumors from a whole-genome sequencing dataset comprising a cohort of 640 patients. We identified 11 of 640 tumors as MMR deficient, but only 2 of 11 exhibited germline mutations in MMR genes or Lynch Syndrome. Two additional tumors had a substantially reduced proportion of mutations attributed to MMR deficiency, where the predominant mutational signatures were related to APOBEC enzymatic activity. Overall, 6 of 11 of the MMR-deficient cases in this cohort were confirmed genetically or epigenetically as having abrogation of MMR genes. However, IHC analysis of MMR-related proteins revealed all but one of 10 samples available for testing as MMR deficient. Thus, the mutational signatures more faithfully reported MMR deficiency than sequencing of MMR genes, because they represent a direct pathophysiologic readout of repair pathway abnormalities. As whole-genome sequencing continues to become more affordable, it could be used to expose individually abnormal tumors in tissue types where MMR deficiency has been rarely detected, but also rarely sought. Cancer Res; 77(18); 4755-62. ©2017 AACR.

HRDetect is a predictor of BRCA1 and BRCA2 deficiency based on mutational signatures.

Approximately 1-5% of breast cancers are attributed to inherited mutations in BRCA1 or BRCA2 and are selectively sensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. In other cancer types, germline and/or somatic mutations in BRCA1 and/or BRCA2 (BRCA1/BRCA2) also confer selective sensitivity to PARP inhibitors. Thus, assays to detect BRCA1/BRCA2-deficient tumors have been sought. Recently, somatic substitution, insertion/deletion and rearrangement patterns, or 'mutational signatures', were associated with BRCA1/BRCA2 dysfunction. Herein we used a lasso logistic regression model to identify six distinguishing mutational signatures predictive of BRCA1/BRCA2 deficiency. A weighted model called HRDetect was developed to accurately detect BRCA1/BRCA2-deficient samples. HRDetect identifies BRCA1/BRCA2-deficient tumors with 98.7% sensitivity (area under the curve (AUC) = 0.98). Application of this model in a cohort of 560 individuals with breast cancer, of whom 22 were known to carry a germline BRCA1 or BRCA2 mutation, allowed us to identify an additional 22 tumors with somatic loss of BRCA1 or BRCA2 and 47 tumors with functional BRCA1/BRCA2 deficiency where no mutation was detected. We validated HRDetect on independent cohorts of breast, ovarian and pancreatic cancers and demonstrated its efficacy in alternative sequencing strategies. Integrating all of the classes of mutational signatures thus reveals a larger proportion of individuals with breast cancer harboring BRCA1/BRCA2 deficiency (up to 22%) than hitherto appreciated (∼1-5%) who could have selective therapeutic sensitivity to PARP inhibition.