Circulating tumor DNA analysis enables molecular characterization of pediatric renal tumors at diagnosis.

Circulating tumor DNA (ctDNA) is a powerful tool for the molecular characterization of cancer. The most frequent pediatric kidney tumors (KT) are Wilms' tumors (WT), but other diagnoses may occur. According to the SIOP strategy, in most countries pediatric KT have a presumptive diagnosis of WT if they are clinically and radiologically compatible. The histologic confirmation is established after post-chemotherapy nephrectomy. Thus, there is a risk for a small fraction of patients to receive neoadjuvant chemotherapy that is not adapted to the disease. The aim of this work is to perform molecular diagnosis of pediatric KT by tumor genetic characterization based on the analysis of ctDNA. We analyzed ctDNA extracted from plasma samples of 18 pediatric patients with KT by whole-exome sequencing and compared the results to their matched tumor and germline DNA. Copy number alterations (CNAs) and single nucleotide variations (SNVs) were analyzed. We were able to detect tumor cell specific genetic alterations-CNAs, SNVs or both-in ctDNA in all patients except in one (for whom the plasma sample was obtained long after nephrectomy). These results open the door to new applications for the study of ctDNA with regards to the molecular diagnosis of KT, with a possibility of its usefulness for adapting the treatment early after diagnosis, but also for disease monitoring and follow up.

Study of chromatin remodeling genes implicates SMARCA4 as a putative player in oncogenesis in neuroblastoma.

In neuroblastoma (NB), genetic alterations in chromatin remodeling (CRGs) and epigenetic modifier genes (EMGs) have been described. We sought to determine their frequency and clinical impact. Whole exome (WES)/whole genome sequencing (WGS) data and targeted sequencing (TSCA®) of exonic regions of 33 CRGs/EMGs were analyzed in tumor samples from 283 NB patients, with constitutional material available for 55 patients. The frequency of CRG/EMG variations in NB cases was then compared to the Genome Aggregation Database (gnomAD). The sequencing revealed SNVs/small InDels or focal CNAs of CRGs/EMGs in 20% (56/283) of all cases, occurring at a somatic level in 4 (7.2%), at a germline level in 12 (22%) cases, whereas for the remaining cases, only tumor material could be analyzed. The most frequently altered genes were ATRX (5%), SMARCA4 (2.5%), MLL3 (2.5%) and ARID1B (2.5%). Double events (SNVs/small InDels/CNAs associated with LOH) were observed in SMARCA4 (n = 3), ATRX (n = 1) and PBRM1 (n = 1). Among the 60 variations, 24 (8.4%) targeted domains of functional importance for chromatin remodeling or highly conserved domains but of unknown function. Variations in SMARCA4 and ATRX occurred more frequently in the NB as compared to the gnomAD control cohort (OR = 4.49, 95%CI: 1.63-9.97, p = 0.038; OR 3.44, 95%CI: 1.46-6.91, p = 0.043, respectively). Cases with CRG/EMG variations showed a poorer overall survival compared to cases without variations. Genetic variations of CRGs/EMGs with likely functional impact were observed in 8.4% (24/283) of NB. Our case-control approach suggests a role of SMARCA4 as a player of NB oncogenesis.

QuantumClone: clonal assessment of functional mutations in cancer based on a genotype-aware method for clonal reconstruction.

Motivation: In cancer, clonal evolution is assessed based on information coming from single nucleotide variants and copy number alterations. Nonetheless, existing methods often fail to accurately combine information from both sources to truthfully reconstruct clonal populations in a given tumor sample or in a set of tumor samples coming from the same patient. Moreover, previously published methods detect clones from a single set of variants. As a result, compromises have to be done between stringent variant filtering [reducing dispersion in variant allele frequency estimates (VAFs)] and using all biologically relevant variants. Results: We present a framework for defining cancer clones using most reliable variants of high depth of coverage and assigning functional mutations to the detected clones. The key element of our framework is QuantumClone, a method for variant clustering into clones based on VAFs, genotypes of corresponding regions and information about tumor purity. We validated QuantumClone and our framework on simulated data. We then applied our framework to whole genome sequencing data for 19 neuroblastoma trios each including constitutional, diagnosis and relapse samples. We confirmed an enrichment of damaging variants within such pathways as MAPK (mitogen-activated protein kinases), neuritogenesis, epithelial-mesenchymal transition, cell survival and DNA repair. Most pathways had more damaging variants in the expanding clones compared to shrinking ones, which can be explained by the increased total number of variants between these two populations. Functional mutational rate varied for ancestral clones and clones shrinking or expanding upon treatment, suggesting changes in clone selection mechanisms at different time points of tumor evolution. Availability and implementation: Source code and binaries of the QuantumClone R package are freely available for download at https://CRAN.R-project.org/package=QuantumClone. Contact: gudrun.schleiermacher@curie.fr or valentina.boeva@inserm.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Molecular profiling of biliary tract cancers reveals distinct genomic landscapes between circulating and tissue tumor DNA.

Biliary tract cancers (BTCs) are heterogeneous malignancies with dismal prognosis due to tumor aggressiveness and poor response to limited current therapeutic options. Tumor exome profiling has allowed to successfully establish targeted therapeutic strategies in the clinical management of cholangiocarcinoma (CCA). Still, whether liquid biopsy profiling could inform on BTC biology and patient management is unknown. In order to test this and generate novel insight into BTC biology, we analyzed the molecular landscape of 128 CCA patients, using a 394-gene NGS panel (Foundation Medicine). Among them, 32 patients had matched circulating tumor (ct) DNA and tumor DNA samples, where both samples were profiled. In both tumor and liquid biopsies, we identified an increased frequency of alterations in genes involved in genome integrity or chromatin remodeling, including ARID1A (15%), PBRM1 (9%), and BAP1 (14%), which were validated using an in-house-developed immunohistochemistry panel. ctDNA and tumor DNA showed variable concordance, with a significant correlation in the total number of detected variants, but some heterogeneity in the detection of actionable mutations. FGFR2 mutations were more frequently identified in liquid biopsies, whereas KRAS alterations were mostly found in tumors. All IDH1 mutations detected in tumor DNA were also identified in liquid biopsies. These findings provide novel insights in the concordance between the tumor and liquid biopsies genomic landscape in a large cohort of patients with BTC and highlight the complementarity of both analyses when guiding therapeutic prescription.

Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in DSRCTs.

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive sarcoma driven by the EWSR1::WT1 chimeric transcription factor. Despite this unique oncogenic driver, DSRCT displays a polyphenotypic differentiation of unknown causality. Using single-cell multi-omics on 12 samples from five patients, we find that DSRCT tumor cells cluster into consistent subpopulations with partially overlapping lineage- and metabolism-related transcriptional programs. In vitro modeling shows that high EWSR1::WT1 DNA-binding activity associates with most lineage-related states, in contrast to glycolytic and profibrotic states. Single-cell chromatin accessibility analysis suggests that EWSR1::WT1 binding site variability may drive distinct lineage-related transcriptional programs, supporting some level of cell-intrinsic plasticity. Spatial transcriptomics reveals that glycolytic and profibrotic states specifically localize within hypoxic niches at the periphery of tumor cell islets, suggesting an additional role of tumor cell-extrinsic microenvironmental cues. We finally identify a single-cell transcriptomics-derived epithelial signature associated with improved patient survival, highlighting the clinical relevance of our findings.

Replication Stress is an Actionable Genetic Vulnerability in Desmoplastic Small Round Cell Tumors.

Desmoplastic small round cell tumor (DSRCT) is an aggressive sarcoma subtype that is driven by the EWS-WT1 chimeric transcription factor. The prognosis for DSRCT is poor, and major advances in treating DSCRT have not occurred for over two decades. To identify effective therapeutic approaches to target DSRCT, we conducted a high-throughput drug sensitivity screen in a DSRCT cell line assessing chemosensitivity profiles for 79 small-molecule inhibitors. DSRCT cells were sensitive to PARP and ATR inhibitors (PARPi, ATRi), as monotherapies and in combination. These effects were recapitulated using multiple clinical PARPi and ATRi in three biologically distinct, clinically-relevant models of DSRCT, including cell lines, a patient-derived xenograft (PDX)-derived organoid model, and a cell line-derived xenograft mouse model. Mechanistically, exposure to a combination of PARPi and ATRi caused increased DNA damage, G2/M checkpoint activation, micronuclei accumulation, replication stress, and R-loop formation. EWS-WT1 silencing abrogated these phenotypes and was epistatic with exogenous expression of the R-loop resolution enzyme RNase H1 in reversing the sensitivity to PARPi and ATRi monotherapies. The combination of PARPi and ATRi also induced EWS-WT1-dependent cell-autonomous activation of the cGAS/STING innate immune pathway and cell surface expression of PD-L1. Taken together, these findings point towards a role for EWS-WT1 in generating R-loop-dependent replication stress that leads to a targetable vulnerability, providing a rationale for the clinical assessment of PARPi and ATRi in DSRCT.

Prognostic value of tumor immune biomarkers in biopsies from patients with refractory solid cancers.

PD-L1 and tumor-infiltrating lymphocytes play a key role in the immune escape of cancer, although their prognostic value remains unknown in patients with refractory solid cancer compared to other known prognostic estimation methods. In this ancillary study, we assessed the prognostic value of previously-defined prognostic scores (such as the Royal Marsden Hospital (RMH) score) and of PD-L1, CD3, CD8 and FOXP3 expressions based on immunohistochemistry (IHC) and RNA sequencing (RNAseq) of tumor samples from patients included in the personalized-medicine MOSCATO-02 trial. We collected biopsies with successful IHC analysis from 266 patients treated between April 2016 and September 2017, among whom 170 (63.9%) also had a matched RNAseq. We used a Random Forest model to identify the best prognostic factor, and a Lasso-penalized Cox model to validate the findings. We found that the RMH score was the strongest prognostic factor, with high scores associated with a higher risk of death (Hazard Ratio (HR)=1.29; CI95%[1.19-1.21]). The PD-L1 expression score obtained from IHC analyses was the second-best performing predictor, with the 1+ score (low expression) linked to a lower risk of death (HR=0.564; CI95%[0.539-0.580]). Other tested variables, including primary tumor type and subsequent treatments received following biopsy, were not found significantly linked to prognosis. We found modest correlation between IHC and RNAseq expressions of immune genes, but RNAseq related better to prognosis. Overall, our study supports the use of the RMH score and the assessment of PD-L1 expression in IHC to estimate prognosis in patients with advanced cancer.

Epigenetic gene alterations in metastatic solid tumours: results from the prospective precision medicine MOSCATO and MATCH-R trials.

INTRODUCTION: Although the role of epigenetic alterations in oncogenesis has been well studied, their prevalence in metastatic solid tumours is still poorly described. We therefore aimed at: (i) describing the presence of epigenetic gene alterations (EGA) - defined by an alteration in a gene encoding an epigenetic regulator; and (ii) evaluating their relationship with clinical characteristics and outcome in patients (pts) included in prospective molecular profiling trials. MATERIALS AND METHODS: On-purpose tumour biopsies from pts with metastatic solid tumours enrolled in the Gustave Roussy-sponsored MOSCATO (NCT01566019) and MATCHR (NCT02517892) trials were molecularly profiled using whole exome sequencing (WES). Alterations in 176 epigenetic genes were assessed and classified as pathogenic variants (PV) or non-pathogenic variants by a molecular tumour board. Clinical characteristics and outcome were collected. RESULTS: Between Dec 2011 and Oct 2016, WES was successfully performed in 292 pts presenting various solid tumours. We found 496 epigenetic gene alterations in 134 patients (49%), including 237 pathogenic variants in 86 patients; 63 tumour samples (47%) presented ≥3 EGAs. The median number of previous treatment lines was 3 (1-10). The most frequently altered genes were KMT2D and KMT2C (16% each), ARID1A and SETD2 (10% each) and KMT2A (8%).; 31% of EGA co-occurred with a driver gene alteration (p < 0.001). Outcome was not correlated with the presence of EGA. CONCLUSIONS: Epigenetic alterations occur frequently in metastatic solid tumours. With the current development of epigenetic modifiers, they increasingly represent actionable targets. Such genes should now be systematically analysed in molecular profiling studies.

PBRM1 Deficiency Confers Synthetic Lethality to DNA Repair Inhibitors in Cancer.

Inactivation of Polybromo 1 (PBRM1), a specific subunit of the PBAF chromatin remodeling complex, occurs frequently in cancer, including 40% of clear cell renal cell carcinomas (ccRCC). To identify novel therapeutic approaches to targeting PBRM1-defective cancers, we used a series of orthogonal functional genomic screens that identified PARP and ATR inhibitors as being synthetic lethal with PBRM1 deficiency. The PBRM1/PARP inhibitor synthetic lethality was recapitulated using several clinical PARP inhibitors in a series of in vitro model systems and in vivo in a xenograft model of ccRCC. In the absence of exogenous DNA damage, PBRM1-defective cells exhibited elevated levels of replication stress, micronuclei, and R-loops. PARP inhibitor exposure exacerbated these phenotypes. Quantitative mass spectrometry revealed that multiple R-loop processing factors were downregulated in PBRM1-defective tumor cells. Exogenous expression of the R-loop resolution enzyme RNase H1 reversed the sensitivity of PBRM1-deficient cells to PARP inhibitors, suggesting that excessive levels of R-loops could be a cause of this synthetic lethality. PARP and ATR inhibitors also induced cyclic GMP-AMP synthase/stimulator of interferon genes (cGAS/STING) innate immune signaling in PBRM1-defective tumor cells. Overall, these findings provide the preclinical basis for using PARP inhibitors in PBRM1-defective cancers. SIGNIFICANCE: This study demonstrates that PARP and ATR inhibitors are synthetic lethal with the loss of PBRM1, a PBAF-specific subunit, thus providing the rationale for assessing these inhibitors in patients with PBRM1-defective cancer. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/11/2888/F1.large.jpg.

ALK mutation dynamics and clonal evolution in a neuroblastoma model exhibiting two ALK mutations.

The ALK gene is a major oncogene of neuroblastoma cases exhibiting ALK activating mutations. Here, we characterized two neuroblastoma cell lines established from a stage 4 patient at diagnosis either from the primary tumor (PT) or from the bone marrow (BM). Both cell lines exhibited similar genomic profiles. All cells in the BM-derived cell line exhibited an ALK F1174L mutation, whereas this mutation was present in only 5% of the cells in the earliest passages of the PT-derived cell line. The BM-derived cell line presented with a higher proliferation rate in vitro and injections in Nude mice resulted in tumor formation only for the BM-derived cell line. Next, we observed that the F1174L mutation frequency in the PT-derived cell line increased with successive passages. Further Whole Exome Sequencing revealed a second ALK mutation, L1196M, in this cell line. Digital droplet PCR documented that the allele fractions of both mutations changed upon passages, and that the F1174L mutation reached 50% in late passages, indicating clonal evolution. In vitro treatment of the PT-derived cell line exhibiting the F1174L and L1196M mutations with the alectinib inhibitor resulted in an enrichment of the L1196M mutation. Using xenografts, we documented a better efficacy of alectinib compared to crizotinib on tumor growth and an enrichment of the L1196M mutation at the end of both treatments. Finally, single-cell RNA-seq analysis was consistent with both mutations resulting in ALK activation. Altogether, this study provides novel insights into ALK mutation dynamics in a neuroblastoma model harbouring two ALK mutations.

Whole-Exome Sequencing of Cell-Free DNA Reveals Temporo-spatial Heterogeneity and Identifies Treatment-Resistant Clones in Neuroblastoma.

Purpose: Neuroblastoma displays important clinical and genetic heterogeneity, with emergence of new mutations at tumor progression.Experimental Design: To study clonal evolution during treatment and follow-up, an innovative method based on circulating cell-free DNA (cfDNA) analysis by whole-exome sequencing (WES) paired with target sequencing was realized in sequential liquid biopsy samples of 19 neuroblastoma patients.Results: WES of the primary tumor and cfDNA at diagnosis showed overlap of single-nucleotide variants (SNV) and copy number alterations, with 41% and 93% of all detected alterations common to the primary neuroblastoma and cfDNA. CfDNA WES at a second time point indicated a mean of 22 new SNVs for patients with progressive disease. Relapse-specific alterations included genes of the MAPK pathway and targeted the protein kinase A signaling pathway. Deep coverage target sequencing of intermediate time points during treatment and follow-up identified distinct subclones. For 17 seemingly relapse-specific SNVs detected by cfDNA WES at relapse but not tumor or cfDNA WES at diagnosis, deep coverage target sequencing detected these alterations in minor subclones, with relapse-emerging SNVs targeting genes of neuritogenesis and cell cycle. Furthermore a persisting, resistant clone with concomitant disappearance of other clones was identified by a mutation in the ubiquitin protein ligase HERC2Conclusions: Modelization of mutated allele fractions in cfDNA indicated distinct patterns of clonal evolution, with either a minor, treatment-resistant clone expanding to a major clone at relapse, or minor clones collaborating toward tumor progression. Identification of treatment-resistant clones will enable development of more efficient treatment strategies. Clin Cancer Res; 24(4); 939-49. ©2017 AACR.

Genomic Copy Number Profiling Using Circulating Free Tumor DNA Highlights Heterogeneity in Neuroblastoma.

PURPOSE: The tumor genomic copy number profile is of prognostic significance in neuroblastoma patients. We have studied the genomic copy number profile of cell-free DNA (cfDNA) and compared this with primary tumor arrayCGH (aCGH) at diagnosis. EXPERIMENTAL DESIGN: In 70 patients, cfDNA genomic copy number profiling was performed using the OncoScan platform. The profiles were classified according to the overall pattern, including numerical chromosome alterations (NCA), segmental chromosome alterations (SCA), and MYCN amplification (MNA). RESULTS: Interpretable and dynamic cfDNA profiles were obtained in 66 of 70 and 52 of 70 cases, respectively. An overall identical genomic profile between tumor aCGH and cfDNA was observed in 47 cases (3 NCAs, 22 SCAs, 22 MNAs). In one case, cfDNA showed an additional SCA not detected by tumor aCGH. In 4 of 8 cases with a silent tumor aCGH profile, cfDNA analysis revealed a dynamic profile (3 SCAs, 1 NCA). In 14 cases, cfDNA analysis did not reveal any copy number changes. A total of 378 breakpoints common to the primary tumor and cfDNA of any given patient were identified, 27 breakpoints were seen by tumor aCGH, and 54 breakpoints were seen in cfDNA only, including two cases with interstitial IGFR1 gains and two alterations targeting TERT CONCLUSIONS: These results demonstrate the feasibility of cfDNA copy number profiling in neuroblastoma patients, with a concordance of the overall genomic profile in aCGH and cfDNA dynamic cases of 97% and a sensitivity of 77%, respectively. Furthermore, neuroblastoma heterogeneity is highlighted, suggesting that cfDNA might reflect genetic alterations of more aggressive cell clones. Clin Cancer Res; 22(22); 5564-73. ©2016 AACRSee related commentary by Janku and Kurzrock, p. 5400.

Relapsed neuroblastomas show frequent RAS-MAPK pathway mutations.

The majority of patients with neuroblastoma have tumors that initially respond to chemotherapy, but a large proportion will experience therapy-resistant relapses. The molecular basis of this aggressive phenotype is unknown. Whole-genome sequencing of 23 paired diagnostic and relapse neuroblastomas showed clonal evolution from the diagnostic tumor, with a median of 29 somatic mutations unique to the relapse sample. Eighteen of the 23 relapse tumors (78%) showed mutations predicted to activate the RAS-MAPK pathway. Seven of these events were detected only in the relapse tumor, whereas the others showed clonal enrichment. In neuroblastoma cell lines, we also detected a high frequency of activating mutations in the RAS-MAPK pathway (11/18; 61%), and these lesions predicted sensitivity to MEK inhibition in vitro and in vivo. Our findings provide a rationale for genetic characterization of relapse neuroblastomas and show that RAS-MAPK pathway mutations may function as a biomarker for new therapeutic approaches to refractory disease.