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.
What's new?
Mutations that affect chromatin remodeling can lead to cancer. In this paper, the authors investigated the impact of variations in chromatin remodeling genes and epigenetic modifier genes on neuroblastoma patients. They compared the frequency of these variations in NB cases with data from the Genome Aggregation Database (gnomAD). Neuroblastoma cases had a higher frequency of SMARCA4 and ATRX gene variations than the general population. Furthermore, NB patients with CRG/EMG mutations had poorer overall survival than NB cases without such mutations. These findings highlight the importance of chromatin remodeling in neuroblastoma as an avenue for new therapeutics.
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.
What's new?
The analysis of cell‐free DNA (cfDNA) in pediatric patients with kidney tumors is of interest for the molecular diagnosis and early identification of aggressive disease subtypes. In this study, significant amounts of circulating tumor DNA (ctDNA) were successfully isolated from pediatric patients with malignant kidney tumors. Accurate copy number profile and mutation calling to identify specific tumor cell genetic alterations was performed with whole‐exome sequencing using only low volumes of blood (≥100 μL) and low amounts of cfDNA (≥4 ng). The findings suggest that ctDNA analysis can facilitate the molecular diagnosis of specific subtypes of pediatric kidney tumors, potentially enabling earlier treatment.
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.
Inactivation of
(
), 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
deficiency. The PBRM1/PARP inhibitor synthetic lethality was recapitulated using several clinical PARP inhibitors in a series of
model systems and
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.
•RMH score is the most important variable associated with prognosis.•A low PD-L1 positive score (PD-L1 1+) is associated with a better prognosis.•CD3, CD8, FoxP3 and PD-L1 expression estimated with ...RNAseq related more to prognosis than with IHC.
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.
Abstract
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.
Supplementary information
Supplementary data are available at Bioinformatics online.
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.
The
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
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.
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.
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.
Display omitted
•DSRCT cell heterogeneity is driven, at least in part, by transcriptional plasticity•Metabolic and lineage programs are two major components of DSRCT heterogeneity•Both cell-intrinsic and -extrinsic mechanisms drive transcriptional variability•ScRNA-seq-derived signatures predict outcome, highlighting intertumor heterogeneity
DSRCT is a rare sarcoma subtype of dismal prognosis driven by the aberrant transcription factor EWSR1::WT1. Using single-cell RNA-sequencing and an in vitro modeling of EWSR1::WT1 activity, Henon et al. show that DSRCT cells are characterized by partly overlapping cell-dependent and microenvironment-dependent transcriptional programs, whose variable expression has prognostic significance.