The pluripotency factor NANOG contributes to mesenchymal plasticity and is predictive for outcome in esophageal adenocarcinoma.

BACKGROUND: Despite the advent of neoadjuvant chemoradiotherapy (CRT), overall survival rates of esophageal adenocarcinoma (EAC) remain low. A readily induced mesenchymal transition of EAC cells contributes to resistance to CRT. METHODS: In this study, we aimed to chart the heterogeneity in cell state transition after CRT and to identify its underpinnings. A panel of 12 esophageal cultures were treated with CRT and ranked by their relative epithelial-mesenchymal plasticity. RNA-sequencing was performed on 100 pre-treatment biopsies. After RNA-sequencing, Ridge regression analysis was applied to correlate gene expression to ranked plasticity, and models were developed to predict mesenchymal transitions in patients. Plasticity score predictions of the three highest significant predictive models were projected on the pre-treatment biopsies and related to clinical outcome data. Motif enrichment analysis of the genes associated with all three models was performed. RESULTS: This study reveals NANOG as the key associated transcription factor predicting mesenchymal plasticity in EAC. Expression of NANOG in pre-treatment biopsies is highly associated with poor response to neoadjuvant chemoradiation, the occurrence of recurrences, and median overall survival difference in EAC patients (>48 months). Perturbation of NANOG reduces plasticity and resensitizes cell lines, organoid cultures, and patient-derived in vivo grafts. CONCLUSIONS: In conclusion, NANOG is a key transcription factor in mesenchymal plasticity in EAC and a promising predictive marker for outcome.

Esophageal cancer is the sixth most common cause of cancer-related death worldwide. Although chemotherapy combined with radiotherapy (chemoradiotherapy) followed by surgery has improved survival, tumor recurrence and metastatic disease (that has spread to other parts of the body) are often observed after several months. In this study, we assessed the effect of chemoradiotherapy on esophageal cells in the lab to predict the effect in patients with esophageal cancer. To investigate this, genes were assessed from 12 different cell lines and 100 patient tissues. We revealed that levels of one of the genes, NANOG, associates with poor response in patients. NANOG could be a promising marker to predict outcome in patients with esophageal cancer. This knowledge might help clinicians to treat patients with esophageal cancer appropriately, or may lead to new or optimized treatments.

Transcriptomic Differences Underlying the Activin-A Induced Large Osteoclast Formation in Both Healthy Control and Fibrodysplasia Ossificans Progressiva Osteoclasts.

Fibrodysplasia Ossificans Progressiva (FOP) is a very rare genetic disease characterized by progressive heterotopic ossification (HO) of soft tissues, leading to immobility and premature death. FOP is caused by a mutation in the Activin receptor Type 1 (ACVR1) gene, resulting in altered responsiveness to Activin-A. We recently revealed that Activin-A induces fewer, but larger and more active, osteoclasts regardless of the presence of the mutated ACVR1 receptor. The underlying mechanism of Activin-A-induced changes in osteoclastogenesis at the gene expression level remains unknown. Transcriptomic changes induced by Activin-A during osteoclast formation from healthy controls and patient-derived CD14-positive monocytes were studied using RNA sequencing. CD14-positive monocytes from six FOP patients and six age- and sex-matched healthy controls were differentiated into osteoclasts in the absence or presence of Activin-A. RNA samples were isolated after 14 days of culturing and analyzed by RNA sequencing. Non-supervised principal component analysis (PCA) showed that samples from the same culture conditions (e.g., without or with Activin-A) tended to cluster, indicating that the variability induced by Activin-A treatment was larger than the variability between the control and FOP samples. RNA sequencing analysis revealed 1480 differentially expressed genes induced by Activin-A in healthy control and FOP osteoclasts with p(adj) < 0.01 and a Log2 fold change of ≥±2. Pathway and gene ontology enrichment analysis revealed several significantly enriched pathways for genes upregulated by Activin-A that could be linked to the differentiation or function of osteoclasts, cell fusion or inflammation. Our data showed that Activin-A has a substantial effect on gene expression during osteoclast formation and that this effect occurred regardless of the presence of the mutated ACVR1 receptor causing FOP.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation.

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Implementation of paediatric precision oncology into clinical practice: The Individualized Therapies for Children with cancer program 'iTHER'.

iTHER is a Dutch prospective national precision oncology program aiming to define tumour molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory paediatric tumours. Between April 2017 and April 2021, 302 samples from 253 patients were included. Comprehensive molecular profiling including low-coverage whole genome sequencing (lcWGS), whole exome sequencing (WES), RNA sequencing (RNA-seq), Affymetrix, and/or 850k methylation profiling was successfully performed for 226 samples with at least 20% tumour content. Germline pathogenic variants were identified in 16% of patients (35/219), of which 22 variants were judged causative for a cancer predisposition syndrome. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with clinical priority very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%) priority. iTHER led to revision or refinement of diagnosis in 8 patients (3.5%). Temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analyses. Of 137 patients with follow-up beyond twelve months, 21 molecularly matched treatments were applied in 19 patients (13.9%), with clinical benefit in few. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates the feasibility of comprehensive molecular profiling across all ages, tumour types and stages in paediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants. Therefore, WES and RNA-seq is nowadays standard clinical care at the Princess Máxima Center for all children with cancer, including patients at primary diagnosis. Improved access to innovative treatments within biology-driven combination trials is required to ultimately improve survival.

Fibroblast activation protein identifies Consensus Molecular Subtype 4 in colorectal cancer and allows its detection by 68Ga-FAPI-PET imaging.

BACKGROUND: In colorectal cancer (CRC), the consensus molecular subtype 4 (CMS4) is associated with therapy resistance and poor prognosis. Clinical diagnosis of CMS4 is hampered by locoregional and temporal variables influencing CMS classification. Diagnostic tools that comprehensively detect CMS4 are therefore urgently needed. METHODS: To identify targets for molecular CMS4 imaging, RNA sequencing data of 3232 primary CRC patients were explored. Heterogeneity of marker expression in relation to CMS4 status was assessed by analysing 3-5 tumour regions and 91.103 single-tumour cells (7 and 29 tumours, respectively). Candidate marker expression was validated in CMS4 peritoneal metastases (PM; n = 59). Molecular imaging was performed using the 68Ga-DOTA-FAPI-46 PET tracer. RESULTS: Fibroblast activation protein (FAP) mRNA identified CMS4 with very high sensitivity and specificity (AUROC > 0.91), and was associated with significantly shorter relapse-free survival (P = 0.0038). Heterogeneous expression of FAP among and within tumour lesions correlated with CMS4 heterogeneity (AUROC = 1.00). FAP expression was homogeneously high in PM, a near-homogeneous CMS4 entity. FAPI-PET identified focal and diffuse PM that were missed using conventional imaging. Extra-peritoneal metastases displayed extensive heterogeneity of tracer uptake. CONCLUSION: FAP expression identifies CMS4 CRC. FAPI-PET may have value in the comprehensive detection of CMS4 tumours in CRC. This is especially relevant in patients with PM, for whom effective imaging tools are currently lacking.

Peritoneal metastases from colorectal cancer belong to Consensus Molecular Subtype 4 and are sensitised to oxaliplatin by inhibiting reducing capacity.

BACKGROUND: Peritoneal metastases (PM) in colorectal cancer (CRC) are associated with therapy resistance and poor survival. Oxaliplatin monotherapy is widely applied in the intraperitoneal treatment of PM, but fails to yield clinical benefit. We aimed to identify the mechanism(s) underlying PM resistance to oxaliplatin and to develop strategies overcoming such resistance. EXPERIMENTAL DESIGN: We generated a biobank consisting of 35 primary tumour regions and 59 paired PM from 12 patients. All samples were analysed by RNA sequencing. We also generated a series of PM-derived organoid (PMDO) cultures and used these to design and test strategies to overcome resistance to oxaliplatin. RESULTS: PM displayed various hallmarks of aggressive CRC biology. The vast majority of PM and paired primary tumours belonged to the Consensus Molecular Subtype 4 (CMS4). PMDO cultures were resistant to oxaliplatin and expressed high levels of glutamate-cysteine ligase (GCLC) causing detoxification of oxaliplatin through glutathione synthesis. Genetic or pharmacological targeting of GCLC sensitised PMDOs to a 1-h exposure to oxaliplatin, through increased platinum-DNA adduct formation. CONCLUSIONS: These results link oxaliplatin resistance of colorectal PM to their CMS4 status and high reducing capacity. Inhibiting the reducing capacity of PM may be an effective strategy to overcome PM resistance to oxaliplatin.

Urinary 3-Methoxytyramine Is a Biomarker for MYC Activity in Patients With Neuroblastoma.

PURPOSE: Elevated urinary 3-methoxytyramine (3MT) level at diagnosis was recently put forward as independent risk factor for poor prognosis in neuroblastoma. Here, we investigated the biologic basis underlying the putative association between elevated 3MT levels and poor prognosis. METHODS: Urinary 3MT levels and prognosis were investigated in both retrospective Italian (N = 90) and prospective Dutch (N = 95) cohorts. From the Dutch Cancer Oncology Group cohort (N = 122), patients with available urinary 3MT and gene expression data (n = 90) were used to generate a 3MT gene signature. The 3MT gene signature score was then used to predict survival outcome in the Children's Oncology Group (N = 247) and German Pediatric Oncology Group (N = 498) cohorts and compared with other known gene signatures. Immunohistochemistry of MYCN and dopamine ß-hydroxylase proteins was performed on primary tumors. RESULTS: Elevated urinary 3MT levels were associated with poor prognosis in a retrospective cohort and a prospective cohort. Moreover, elevated urinary 3MT levels were associated with eight differentially expressed genes, providing a 3MT gene signature that successfully predicted poor clinical outcome. Even among low-risk patients, high 3MT signature score was associated with poor 5-year overall survival (72% v 99% among low-risk patients with a low 3MT signature score), and the 3MT signature score was correlated with MYC activity in the tumor (R = 82%, P < .0001). Finally, a strong MYCN and weak dopamine ß-hydroxylase staining of tumors derived from patients with elevated urinary 3MT levels was observed, linking MYC activity in the tumor to both catecholamine biosynthesis and elevated urinary 3MT levels. CONCLUSION: Elevated urinary 3MT is a promising biomarker for poor prognosis and reflects increased MYC activity in the tumor. Therefore, urinary 3MT levels should be measured at diagnosis and may assist in assessing risk.

Mesenchymal-Type Neuroblastoma Cells Escape ALK Inhibitors.

Cancer therapy frequently fails due to the emergence of resistance. Many tumors include phenotypically immature tumor cells, which have been implicated in therapy resistance. Neuroblastoma cells can adopt a lineage-committed adrenergic (ADRN) or an immature mesenchymal (MES) state. They differ in epigenetic landscape and transcription factors, and MES cells are more resistant to chemotherapy. Here we analyzed the response of MES cells to targeted drugs. Activating anaplastic lymphoma kinase (ALK) mutations are frequently found in neuroblastoma and ALK inhibitors (ALKi) are in clinical trials. ALKi treatment of ADRN neuroblastoma cells with a tumor-driving ALK mutation induced cell death. Conversely, MES cells did not express either mutant or wild-type ALK and were resistant to ALKi, and MES cells formed tumors that progressed under ALKi therapy. In assessing the role of MES cells in relapse development, TRAIL was identified to specifically induce apoptosis in MES cells and to suppress MES tumor growth. Addition of TRAIL to ALKi treatment of neuroblastoma xenografts delayed relapses in a subset of the animals, suggesting a role for MES cells in relapse formation. While ADRN cells resembled normal embryonal neuroblasts, MES cells resembled immature precursor cells, which also lacked ALK expression. Resistance to targeted drugs can therefore be an intrinsic property of immature cancer cells based on their resemblance to developmental precursors. SIGNIFICANCE: In neuroblastoma, mesenchymal tumor cells lack expression of the tumor-driving ALK oncogene and are resistant to ALKi, but dual treatment with ALKi and mesenchymal cell-targeting TRAIL delays tumor relapse.

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome.

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

Author Correction: c-Kit-positive ILC2s exhibit an ILC3-like signature that may contribute to IL-17-mediated pathologies.

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

c-Kit-positive ILC2s exhibit an ILC3-like signature that may contribute to IL-17-mediated pathologies.

Here we identify a group 2 innate lymphoid cell (ILC2) subpopulation that can convert into interleukin-17 (IL-17)-producing NKp44- ILC3-like cells. c-Kit and CCR6 define this ILC2 subpopulation that exhibits ILC3 features, including RORγt, enabling the conversion into IL-17-producing cells in response to IL-1ß and IL-23. We also report a role for transforming growth factor-ß in promoting the conversion of c-Kit- ILC2s into RORγt-expressing cells by inducing the upregulation of IL23R, CCR6 and KIT messenger RNA in these cells. This switch was dependent on RORγt and the downregulation of GATA-3. IL-4 was able to reverse this event, supporting a role for this cytokine in maintaining ILC2 identity. Notably, this plasticity has physiological relevance because a subset of RORγt+ ILC2s express the skin-homing receptor CCR10, and the frequencies of IL-17-producing ILC3s are increased at the expense of ILC2s within the lesional skin of patients with psoriasis.

A NOTCH feed-forward loop drives reprogramming from adrenergic to mesenchymal state in neuroblastoma.

Transition between differentiation states in development occurs swift but the mechanisms leading to epigenetic and transcriptional reprogramming are poorly understood. The pediatric cancer neuroblastoma includes adrenergic (ADRN) and mesenchymal (MES) tumor cell types, which differ in phenotype, super-enhancers (SEs) and core regulatory circuitries. These cell types can spontaneously interconvert, but the mechanism remains largely unknown. Here, we unravel how a NOTCH3 intracellular domain reprogrammed the ADRN transcriptional landscape towards a MES state. A transcriptional feed-forward circuitry of NOTCH-family transcription factors amplifies the NOTCH signaling levels, explaining the swift transition between two semi-stable cellular states. This transition induces genome-wide remodeling of the H3K27ac landscape and a switch from ADRN SEs to MES SEs. Once established, the NOTCH feed-forward loop maintains the induced MES state. In vivo reprogramming of ADRN cells shows that MES and ADRN cells are equally oncogenic. Our results elucidate a swift transdifferentiation between two semi-stable epigenetic cellular states.

RAS-MAPK Pathway-Driven Tumor Progression Is Associated with Loss of CIC and Other Genomic Aberrations in Neuroblastoma.

Mutations affecting the RAS-MAPK pathway frequently occur in relapsed neuroblastoma tumors, which suggests that activation of this pathway is associated with a more aggressive phenotype. To explore this hypothesis, we generated several model systems to define a neuroblastoma RAS-MAPK pathway signature. Activation of this pathway in primary tumors indeed correlated with poor survival and was associated with known activating mutations in ALK and other RAS-MAPK pathway genes. Integrative analysis showed that mutations in PHOX2B, CIC, and DMD were also associated with an activated RAS-MAPK pathway. Mutation of PHOX2B and deletion of CIC in neuroblastoma cell lines induced activation of the RAS-MAPK pathway. This activation was independent of phosphorylated ERK in CIC knockout systems. Furthermore, deletion of CIC caused a significant increase in tumor growth in vivo These results show that the RAS-MAPK pathway is involved in tumor progression and establish CIC as a powerful tumor suppressor that functions downstream of this pathway in neuroblastoma.Significance: This work identifies CIC as a powerful tumor suppressor affecting the RAS-MAPK pathway in neuroblastoma and reinforces the importance of mutation-driven activation of this pathway in cancer. Cancer Res; 78(21); 6297-307. ©2018 AACR.

Neuroblastoma is composed of two super-enhancer-associated differentiation states.

Neuroblastoma and other pediatric tumors show a paucity of gene mutations, which has sparked an interest in their epigenetic regulation. Several tumor types include phenotypically divergent cells, resembling cells from different lineage development stages. It has been proposed that super-enhancer-associated transcription factor (TF) networks underlie lineage identity, but the role of these enhancers in intratumoral heterogeneity is unknown. Here we show that most neuroblastomas include two types of tumor cells with divergent gene expression profiles. Undifferentiated mesenchymal cells and committed adrenergic cells can interconvert and resemble cells from different lineage differentiation stages. ChIP-seq analysis of isogenic pairs of mesenchymal and adrenergic cells identified a distinct super-enhancer landscape and super-enhancer-associated TF network for each cell type. Expression of the mesenchymal TF PRRX1 could reprogram the super-enhancer and mRNA landscapes of adrenergic cells toward a mesenchymal state. Mesenchymal cells were more chemoresistant in vitro and were enriched in post-therapy and relapse tumors. Two super-enhancer-associated TF networks, which probably mediate lineage control in normal development, thus dominate epigenetic control of neuroblastoma and shape intratumoral heterogeneity.

Atypical Teratoid/Rhabdoid Tumors Are Comprised of Three Epigenetic Subgroups with Distinct Enhancer Landscapes.

Atypical teratoid/rhabdoid tumor (ATRT) is one of the most common brain tumors in infants. Although the prognosis of ATRT patients is poor, some patients respond favorably to current treatments, suggesting molecular inter-tumor heterogeneity. To investigate this further, we genetically and epigenetically analyzed 192 ATRTs. Three distinct molecular subgroups of ATRTs, associated with differences in demographics, tumor location, and type of SMARCB1 alterations, were identified. Whole-genome DNA and RNA sequencing found no recurrent mutations in addition to SMARCB1 that would explain the differences between subgroups. Whole-genome bisulfite sequencing and H3K27Ac chromatin-immunoprecipitation sequencing of primary tumors, however, revealed clear differences, leading to the identification of subgroup-specific regulatory networks and potential therapeutic targets.

TERT rearrangements are frequent in neuroblastoma and identify aggressive tumors.

Whole-genome sequencing detected structural rearrangements of TERT in 17 of 75 high-stage neuroblastomas, with five cases resulting from chromothripsis. Rearrangements were associated with increased TERT expression and targeted regions immediately up- and downstream of TERT, positioning a super-enhancer close to the breakpoints in seven cases. TERT rearrangements (23%), ATRX deletions (11%) and MYCN amplifications (37%) identify three almost non-overlapping groups of high-stage neuroblastoma, each associated with very poor prognosis.

Molecular Classification of Ependymal Tumors across All CNS Compartments, Histopathological Grades, and Age Groups.

Ependymal tumors across age groups are currently classified and graded solely by histopathology. It is, however, commonly accepted that this classification scheme has limited clinical utility based on its lack of reproducibility in predicting patients' outcome. We aimed at establishing a uniform molecular classification using DNA methylation profiling. Nine molecular subgroups were identified in a large cohort of 500 tumors, 3 in each anatomical compartment of the CNS, spine, posterior fossa, supratentorial. Two supratentorial subgroups are characterized by prototypic fusion genes involving RELA and YAP1, respectively. Regarding clinical associations, the molecular classification proposed herein outperforms the current histopathological classification and thus might serve as a basis for the next World Health Organization classification of CNS tumors.

Mutations in the SIX1/2 pathway and the DROSHA/DGCR8 miRNA microprocessor complex underlie high-risk blastemal type Wilms tumors.

Blastemal histology in chemotherapy-treated pediatric Wilms tumors (nephroblastoma) is associated with adverse prognosis. To uncover the underlying tumor biology and find therapeutic leads for this subgroup, we analyzed 58 blastemal type Wilms tumors by exome and transcriptome sequencing and validated our findings in a large replication cohort. Recurrent mutations included a hotspot mutation (Q177R) in the homeo-domain of SIX1 and SIX2 in tumors with high proliferative potential (18.1% of blastemal cases); mutations in the DROSHA/DGCR8 microprocessor genes (18.2% of blastemal cases); mutations in DICER1 and DIS3L2; and alterations in IGF2, MYCN, and TP53, the latter being strongly associated with dismal outcome. DROSHA and DGCR8 mutations strongly altered miRNA expression patterns in tumors, which was functionally validated in cell lines expressing mutant DROSHA.

Genome sequencing of SHH medulloblastoma predicts genotype-related response to smoothened inhibition.

Smoothened (SMO) inhibitors recently entered clinical trials for sonic-hedgehog-driven medulloblastoma (SHH-MB). Clinical response is highly variable. To understand the mechanism(s) of primary resistance and identify pathways cooperating with aberrant SHH signaling, we sequenced and profiled a large cohort of SHH-MBs (n = 133). SHH pathway mutations involved PTCH1 (across all age groups), SUFU (infants, including germline), and SMO (adults). Children >3 years old harbored an excess of downstream MYCN and GLI2 amplifications and frequent TP53 mutations, often in the germline, all of which were rare in infants and adults. Functional assays in different SHH-MB xenograft models demonstrated that SHH-MBs harboring a PTCH1 mutation were responsive to SMO inhibition, whereas tumors harboring an SUFU mutation or MYCN amplification were primarily resistant.