Amplification or overexpression of MYCN is associated with poor prognosis of human neuroblastoma. We have recently defined a MYCN-dependent transcriptional signature, including protein arginine methyltransferase 1 (PRMT1), which identifies a subgroup of patients with high-risk disease. Here we provide several lines of evidence demonstrating PRMT1 as a novel regulator of MYCN and implicating PRMT1 as a potential therapeutic target in neuroblastoma pathogenesis. First, we observed a strong correlation between MYCN and PRMT1 protein levels in primary neuroblastoma tumors. Second, MYCN physically associates with PRMT1 by direct protein-protein interaction. Third, depletion of PRMT1 through siRNA knockdown reduced neuroblastoma cell viability and MYCN expression. Fourth, we showed that PRMT1 regulates MYCN stability and identified MYCN as a novel substrate of PRMT1. Finally, we demonstrated that mutation of putatively methylated arginine R65 to alanine decreased MYCN stability by altering phosphorylation at residues serine 62 and threonine 58. These results provide mechanistic insights into the modulation of MYCN oncoprotein by PRMT1, and suggest that targeting PRMT1 may have a therapeutic impact on MYCN-driven oncogenesis.

Proliferation of the self-renewing epithelium of the gastric corpus occurs almost exclusively in the isthmus of the glands, from where cells migrate bidirectionally toward pit and base. The isthmus is therefore generally viewed as the stem cell zone. We find that the stem cell marker Troy is expressed at the gland base by a small subpopulation of fully differentiated chief cells. By lineage tracing with a Troy-eGFP-ires-CreERT2 allele, single marked chief cells are shown to generate entirely labeled gastric units over periods of months. This phenomenon accelerates upon tissue damage. Troy(+) chief cells can be cultured to generate long-lived gastric organoids. Troy marks a specific subset of chief cells that display plasticity in that they are capable of replenishing entire gastric units, essentially serving as quiescent "reserve" stem cells. These observations challenge the notion that stem cell hierarchies represent a "one-way street."

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

Neuroblastoma (NB) is the most common extracranial pediatric tumor. NB patients over 18 months of age at the time of diagnosis are often in the later stages of the disease, present with widespread dissemination, and often possess MYCN tumor gene amplification. MYCN is a transcription factor that regulates the expression of a number of genes including ornithine decarboxylase (ODC), a rate-limiting enzyme in the biosynthesis of polyamines. Inhibiting ODC in NB cells produces many deleterious effects including G(1) cell cycle arrest, inhibition of cell proliferation, and decreased tumor growth, making ODC a promising target for drug interference. DFMO treatment leads to the accumulation of the cyclin-dependent kinase inhibitor p27(Kip1) protein and causes p27(Kip1)/Rb-coupled G(1) cell cycle arrest in MYCN-amplified NB tumor cells through a process that involves p27(Kip1) phosphorylation at residues Ser10 and Thr198. While p27(Kip1) is well known for its role as a cyclin-dependent kinase inhibitor, recent studies have revealed a novel function of p27(Kip1) as a regulator of cell migration and invasion. In the present study we found that p27(Kip1) regulates the migration and invasion in NB and that these events are dependent on the state of phosphorylation of p27(Kip1). DFMO treatments induced MYCN protein downregulation and phosphorylation of Akt/PKB (Ser473) and GSK3-β (Ser9), and polyamine supplementation alleviated the DFMO-induced effects. Importantly, we provide strong evidence that p27(Kip1) mRNA correlates with clinical features and the survival probability of NB patients.

The neurotrophin-3 (NT-3) receptor tropomyosin receptor kinase C (TrkC/NTRK3) has been described as a dependence receptor and, as such, triggers apoptosis in the absence of its ligand NT-3. This proapoptotic activity has been proposed to confer a tumor suppressor activity to this classic tyrosine kinase receptor (RTK). By investigating interacting partners that might facilitate TrkC-induced cell death, we have identified the basic helix-loop-helix (bHLH) transcription factor Hey1 and importin-α3 (karyopherin alpha 4 [KPNA4]) as direct interactors of TrkC intracellular domain, and we show that Hey1 is required for TrkC-induced apoptosis. We propose here that the cleaved proapoptotic portion of TrkC intracellular domain (called TrkC killer-fragment [TrkC-KF]) is translocated to the nucleus by importins and interacts there with Hey1. We also demonstrate that Hey1 and TrkC-KF transcriptionally silence mouse double minute 2 homolog (MDM2), thus contributing to p53 stabilization. p53 transcriptionally regulates the expression of TrkC-KF cytoplasmic and mitochondrial interactors cofactor of breast cancer 1 (COBRA1) and B cell lymphoma 2-associated X (BAX), which will subsequently trigger the intrinsic pathway of apoptosis. Of interest, TrkC was proposed to constrain tumor progression in neuroblastoma (NB), and we demonstrate in an avian model that TrkC tumor suppressor activity requires Hey1 and p53.

Neuroblastomas are childhood tumors with frequent fatal relapses after induction treatment, which is related to tumor evolution with additional genomic events. Our whole-genome sequencing data analysis revealed a high frequency of somatic cytosine > adenine (C > A) substitutions in primary neuroblastoma tumors, which was associated with poor survival. We showed that increased levels of C > A substitutions correlate with copy number loss (CNL) of OGG1 or MUTYH Both genes encode DNA glycosylases that recognize 8-oxo-guanine (8-oxoG) lesions as a first step of 8-oxoG repair. Tumor organoid models with CNL of OGG1 or MUTYH show increased 8-oxoG levels compared to wild-type cells. We used CRISPR-Cas9 genome editing to create knockout clones of MUTYH and OGG1 in neuroblastoma cells. Whole-genome sequencing of single-cell OGG1 and MUTYH knockout clones identified an increased accumulation of C > A substitutions. Mutational signature analysis of these OGG1 and MUTYH knockout clones revealed enrichment for C > A signatures 18 and 36, respectively. Clustering analysis showed that the knockout clones group together with tumors containing OGG1 or MUTYH CNL. In conclusion, we demonstrate that defects in 8-oxoG repair cause accumulation of C > A substitutions in neuroblastoma, which contributes to mutagenesis and tumor evolution.

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.

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.

In the past decades, the incidence of esophageal adenocarcinoma has increased dramatically in Western populations. Better understanding of disease etiology along with the identification of novel prognostic and predictive biomarkers are urgently needed to improve the dismal survival probabilities. Here, we performed comprehensive RNA (coding and non-coding) profiling in various samples from 17 patients diagnosed with esophageal adenocarcinoma, high-grade dysplastic or non-dysplastic Barrett's esophagus. Per patient, a blood plasma sample, and a healthy and disease esophageal tissue sample were included. In total, this comprehensive dataset consists of 102 sequenced libraries from 51 samples. Based on this data, 119 expression profiles are available for three biotypes, including miRNA (51), mRNA (51) and circRNA (17). This unique resource allows for discovery of novel biomarkers and disease mechanisms, comparison of tissue and liquid biopsy profiles, integration of coding and non-coding RNA patterns, and can serve as a validation dataset in other RNA landscaping studies. Moreover, structural RNA differences can be identified in this dataset, including protein coding mutations, fusion genes, and circular RNAs.

SP140 is a bromodomain-containing protein expressed predominantly in immune cells. Genetic polymorphisms and epigenetic modifications in the SP140 locus have been linked to Crohn's disease (CD), suggesting a role in inflammation.

Mutations affecting the RAS-MAPK pathway occur frequently in relapsed neuroblastoma tumors and are associated with response to MEK inhibition in vitro. However, these inhibitors alone do not lead to tumor regression in vivo, indicating the need for combination therapy.

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

In principle, germline cells possess the capability to transmit a nearly unaltered set of genetic material to infinite future generations, whereas somatic cells are limited by strict growth constraints necessary to assure an organism's physical structure and eventual mortality. As the potential to replicate indefinitely is a key feature of cancer, we hypothesized that the activation of a "germline program" in somatic cells can contribute to oncogenesis. Our group recently described over one thousand germline specific genes that can be ectopically expressed in cancer, yet how germline specific processes contribute to the malignant properties of cancer is poorly understood. We here show that the expression of germ cell/cancer (GC) genes correlates with malignancy in lung adenocarcinoma (LUAD). We found that LUAD cells expressing more GC genes can repair DNA double strand breaks more rapidly, show higher rates of proliferation and are more resistant to ionizing radiation, compared to LUAD cells that express fewer GC genes. In particular, we identified the HORMA domain protein regulator TRIP13 to be predominantly responsible for this malignant phenotype, and that TRIP13 inhibition or expression levels affect the response to ionizing radiation and subsequent DNA repair. Our results demonstrate that GC genes are viable targets in oncology, as they induce increased radiation resistance and increased propagation in cancer cells. Because their expression is normally restricted to germline cells, we anticipate that GC gene directed therapeutic options will effectively target cancer, with limited side effects besides (temporary) infertility.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

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.

Tumor-educated blood platelets (TEPs) are implicated as central players in the systemic and local responses to tumor growth, thereby altering their RNA profile. We determined the diagnostic potential of TEPs by mRNA sequencing of 283 platelet samples. We distinguished 228 patients with localized and metastasized tumors from 55 healthy individuals with 96% accuracy. Across six different tumor types, the location of the primary tumor was correctly identified with 71% accuracy. Also, MET or HER2-positive, and mutant KRAS, EGFR, or PIK3CA tumors were accurately distinguished using surrogate TEP mRNA profiles. Our results indicate that blood platelets provide a valuable platform for pan-cancer, multiclass cancer, and companion diagnostics, possibly enabling clinical advances in blood-based "liquid biopsies".

Medulloblastoma is the most common malignant brain tumor in children. Despite recent improvements in cure rates, prediction of disease outcome remains a major challenge and survivors suffer from serious therapy-related side-effects. Recent data showed that patients with WNT-activated tumors have a favorable prognosis, suggesting that these patients could be treated less intensively, thereby reducing the side-effects. This illustrates the potential benefits of a robust classification of medulloblastoma patients and a detailed knowledge of associated biological mechanisms.

The Hox clusters play a crucial role in body patterning during animal development. They encode both Hox transcription factor and micro-RNA genes that are activated in a precise temporal and spatial sequence that follows their chromosomal order. These remarkable collinear properties confer functional unit status for Hox clusters. We developed the TranscriptView platform to establish high resolution transcriptional profiling and report here that transcription in the Hox clusters is far more complex than previously described in both human and mouse. Unannotated transcripts can represent up to 60% of the total transcriptional output of a cluster. In particular, we identified 14 non-coding Transcriptional Units antisense to Hox genes, 10 of which (70%) have a detectable mouse homolog. Most of these Transcriptional Units in both human and mouse present conserved sizeable sequences (>40 bp) overlapping Hox transcripts, suggesting that these Hox antisense transcripts are functional. Hox clusters also display at least seven polycistronic clusters, i.e., different genes being co-transcribed on long isoforms (up to 30 kb). This work provides a reevaluated framework for understanding Hox gene function and dys-function. Such extensive transcriptions may provide a structural explanation for Hox clustering.

Somatic mutations in cancer can result in neoantigens against which patients can be vaccinated. The quest for tumor specific neoantigens has yielded no targets that are common to all tumors, yet foreign to healthy cells. Single base pair substitutions (SNVs) at best can alter 1 amino acid which can result in a neoantigen; with the exception of rare site-specific oncogenic driver mutations (such as RAS) such mutations are private. Here, we describe a source of common neoantigens induced by frame shift mutations, based on analysis of 10,186 TCGA tumor samples. We find that these frame shift mutations can produce long neoantigens. These are completely new to the body, and indeed recent evidence suggests that frame shifts can be highly immunogenic. We report that many different frame shift mutations converge to the same small set of 3' neo open reading frame peptides (NOPs), all encoded by the Neo-ORFeome. We find that a fixed set of only 1,244 neo-peptides in as much as 30% of all TCGA cancer patients. For some tumor classes this is higher; e.g. for colon and cervical cancer, peptides derived from only ten genes (saturated at 90 peptides) can be applied to 39% of all patients. 50% of all TCGA patients can be achieved at saturation (using all those peptides in the library found more than once). A pre-fabricated library of vaccines (peptide, RNA or DNA) based on this set can provide off the shelf, quality certified, 'personalized' vaccines within hours, saving months of vaccine preparation. This is crucial for critically ill cancer patients with short average survival expectancy after diagnosis.

The memory CD8 T-cell pool must select for clones that bind immunodominant epitopes with high affinity to efficiently counter reinfection. At the same time, it must retain a level of clonal diversity to allow recognition of pathogens with mutated epitopes. How the level of diversity within the memory pool is controlled is unclear, especially in the context of a selective drive for antigen affinity. We find that preservation of clones that bind the activating antigen with low affinity depends on expression of the transcription factor Eomes in the first days after antigen encounter. Eomes is induced at low activating signal strength and directly drives transcription of the prosurvival protein Bcl-2. At higher signal intensity, T-bet is induced which suppresses Bcl-2 and causes a relative survival advantage for cells of low affinity. Clones activated with high-affinity antigen form memory largely independent of Eomes and have a proliferative advantage over clones that bind the same antigen with low affinity. This causes high-affinity clones to prevail in the memory pool, despite their relative survival deficit. Genetic or therapeutic targeting of the Eomes/Bcl-2 axis reduces the clonal diversity of the memory pool, which diminishes its ability to respond to pathogens carrying mutations in immunodominant epitopes. Thus, we demonstrate on a molecular level how sufficient diversity of the memory pool is established in an environment of affinity-based selection.