Comprehensive genome wide analyses of single cells became increasingly important in cancer research, but remain to be a technically challenging task. Here, we provide a protocol for array comparative genomic hybridization (aCGH) of single cells. The protocol is based on an established adapter-linker PCR (WGAM) and allowed us to detect copy number alterations as small as 56 kb in single cells. In addition we report on factors influencing the success of single cell aCGH downstream of the amplification method, including the characteristics of the reference DNA, the labeling technique, the amount of input DNA, reamplification, the aCGH resolution, and data analysis. In comparison with two other commercially available non-linear single cell amplification methods, WGAM showed a very good performance in aCGH experiments. Finally, we demonstrate that cancer cells that were processed and identified by the CellSearch® System and that were subsequently isolated from the CellSearch® cartridge as single cells by fluorescence activated cell sorting (FACS) could be successfully analyzed using our WGAM-aCGH protocol. We believe that even in the era of next-generation sequencing, our single cell aCGH protocol will be a useful and (cost-) effective approach to study copy number alterations in single cells at resolution comparable to those reported currently for single cell digital karyotyping based on next generation sequencing data.

Resistance to standard cisplatin-based chemotherapies leads to worse survival outcomes for patients with esophageal squamous cell carcinoma (ESCC). Therefore, there is an urgent need to understand the aberrant mechanisms driving resistance in ESCC tumors. We hypothesized that ubiquilin-4 (UBQLN4), a protein that targets ubiquitinated proteins to the proteasome, regulates the expression of Meiotic Recombination 11 Homolog A (MRE11A), a critical component of the MRN complex and DNA damage repair pathways. Initially, immunohistochemistry analysis was conducted in specimens from patients with ESCC (n = 120). In endoscopic core ESCC biopsies taken from 61 patients who underwent neoadjuvant chemotherapy (NAC) (5-fluorouracil and cisplatin), low MRE11A and high UBQLN4 protein levels were associated with reduced pathological response to NAC (P < 0.001 and P < 0.001, respectively). Multivariable analysis of surgically resected ESCC tissues from 59 patients revealed low MRE11A and high UBLQN4 expression as independent factors that can predict shorter overall survival [P = 0.01, hazard ratio (HR) = 5.11, 95% confidence interval (CI), 1.45-18.03; P = 0.02, HR = 3.74, 95% CI, 1.19-11.76, respectively]. Suppression of MRE11A expression was associated with cisplatin resistance in ESCC cell lines. Additionally, MRE11A was found to be ubiquitinated after cisplatin treatment. We observed an amplification of UBQLN4 gene copy numbers and an increase in UBQLN4 protein levels in ESCC tissues. Binding of UBQLN4 to ubiquitinated-MRE11A increased MRE11A degradation, thereby regulating MRE11A protein levels following DNA damage and promoting cisplatin resistance. In summary, MRE11A and UBQLN4 protein levels can serve as predictors for NAC response and as prognostic markers in ESCC patients.

The majority of high-risk neuroblastomas harbor telomerase activity, and telomerase-interacting compounds, such as 6-thio-2'-deoxyguanosine (6-thio-dG), have been found to impair the growth of telomerase-positive neuroblastoma cell lines. It has remained unclear, however, how such drugs can be combined with other compounds used in current treatment concepts for neuroblastoma patients.

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

Replication competent vesicular stomatitis virus (VSV) is the basis of a vaccine against Ebola and VSV strains are developed as oncolytic viruses. Both functions depend on the ability of VSV to induce adequate amounts of interferon-α/β. It is therefore important to understand how VSV triggers interferon responses. VSV activates innate immunity via retinoic acid-inducible gene I (RIG-I), a sensor for viral RNA. Our results show that VSV needs to replicate for a robust interferon response. Analysis of RIG-I-associated RNA identified a copy-back defective-interfering (DI) genome and full-length viral genomes as main trigger of RIG-I. VSV stocks depleted of DI genomes lost most of their interferon-stimulating activity. The remaining full-length genome and leader-N-read-through sequences, however, still triggered RIG-I. Awareness for DI genomes as trigger of innate immune responses will help to standardize DI genome content and to purposefully deplete or use DI genomes as natural adjuvants in VSV-based therapeutics.

The mechanism for generating double minutes chromosomes (dmin) and homogeneously staining regions (hsr) in cancer is still poorly understood. Through an integrated approach combining next-generation sequencing, single nucleotide polymorphism array, fluorescent in situ hybridization and polymerase chain reaction-based techniques, we inferred the fine structure of MYC-containing dmin/hsr amplicons harboring sequences from several different chromosomes in seven tumor cell lines, and characterized an unprecedented number of hsr insertion sites. Local chromosome shattering involving a single-step catastrophic event (chromothripsis) was recently proposed to explain clustered chromosomal rearrangements and genomic amplifications in cancer. Our bioinformatics analyses based on the listed criteria to define chromothripsis led us to exclude it as the driving force underlying amplicon genesis in our samples. Instead, the finding of coexisting heterogeneous amplicons, differing in their complexity and chromosome content, in cell lines derived from the same tumor indicated the occurrence of a multi-step evolutionary process in the genesis of dmin/hsr. Our integrated approach allowed us to gather a complete view of the complex chromosome rearrangements occurring within MYC amplicons, suggesting that more than one model may be invoked to explain the origin of dmin/hsr in cancer. Finally, we identified PVT1 as a target of fusion events, confirming its role as breakpoint hotspot in MYC amplification.

20% of children suffering from high hyperdiploid acute lymphoblastic leukemia develop recurrent disease. The molecular mechanisms are largely unknown. Here, we analyzed the genetic landscape of five patients at relapse, who developed recurrent disease without prior high-risk indication using whole-exome- and whole-genome-sequencing. Oncogenic mutations of RAS pathway genes (NRAS, KRAS, FLT3, n=4) and deactivating mutations of major epigenetic regulators (CREBBP, EP300, each n=2 and ARID4B, EZH2, MACROD2, MLL2, each n=1) were prominent in these cases and virtually absent in non-recurrent cases (n=6) or other pediatric acute lymphoblastic leukemia cases (n=18). In relapse nucleotide variations were detected in cell fate determining transcription factors (GLIS1, AKNA). Structural genomic alterations affected genes regulating B-cell development (IKZF1, PBX1, RUNX1). Eleven novel translocations involved the genes ART4, C12orf60, MACROD2, TBL1XR1, LRRN4, KIAA1467, and ELMO1/MIR1200. Typically, patients harbored only single structural variations, except for one patient who displayed massive rearrangements in the context of a germline tumor suppressor TP53 mutation and a Li-Fraumeni syndrome-like family history. Another patient harbored a germline mutation in the DNA repair factor ATM. In summary, the relapse patients of our cohort were characterized by somatic mutations affecting the RAS pathway, epigenetic and developmental programs and germline mutations in DNA repair pathways.

Long intergenic non-coding RNAs (lincRNAs) are emerging as integral components of signaling pathways in various cancer types. In neuroblastoma, only a handful of lincRNAs are known as upstream regulators or downstream effectors of oncogenes. Here, we exploit RNA sequencing data of primary neuroblastoma tumors, neuroblast precursor cells, neuroblastoma cell lines and various cellular perturbation model systems to define the neuroblastoma lincRNome and map lincRNAs up- and downstream of neuroblastoma driver genes MYCN, ALK and PHOX2B. Each of these driver genes controls the expression of a particular subset of lincRNAs, several of which are associated with poor survival and are differentially expressed in neuroblastoma tumors compared to neuroblasts. By integrating RNA sequencing data from both primary tumor tissue and cancer cell lines, we demonstrate that several of these lincRNAs are expressed in stromal cells. Deconvolution of primary tumor gene expression data revealed a strong association between stromal cell composition and driver gene status, resulting in differential expression of these lincRNAs. We also explored lincRNAs that putatively act upstream of neuroblastoma driver genes, either as presumed modulators of driver gene activity, or as modulators of effectors regulating driver gene expression. This analysis revealed strong associations between the neuroblastoma lincRNAs MIAT and MEG3 and MYCN and PHOX2B activity or expression. Together, our results provide a comprehensive catalogue of the neuroblastoma lincRNome, highlighting lincRNAs up- and downstream of key neuroblastoma driver genes. This catalogue forms a solid basis for further functional validation of candidate neuroblastoma lincRNAs.

Neuroblastoma is an aggressive childhood cancer arising from sympatho-adrenergic neuronal progenitors. The low survival rates for high-risk disease point to an urgent need for novel targeted therapeutic approaches. Detailed molecular characterization of the neuroblastoma genomic landscape indicates that ALK-activating mutations are present in 10% of primary tumours. Together with other mutations causing RAS/MAPK pathway activation, ALK mutations are also enriched in relapsed cases and ALK activation was shown to accelerate MYCN-driven tumour formation through hitherto unknown ALK-driven target genes. To gain further insight into how ALK contributes to neuroblastoma aggressiveness, we searched for known oncogenes in our previously reported ALK-driven gene signature. We identified ETV5, a bona fide oncogene in prostate cancer, as robustly upregulated in neuroblastoma cells harbouring ALK mutations, and show high ETV5 levels downstream of the RAS/MAPK axis. Increased ETV5 expression significantly impacted migration, invasion and colony formation in vitro, and ETV5 knockdown reduced proliferation in a murine xenograft model. We also established a gene signature associated with ETV5 knockdown that correlates with poor patient survival. Taken together, our data highlight ETV5 as an intrinsic component of oncogenic ALK-driven signalling through the MAPK axis and propose that ETV5 upregulation in neuroblastoma may contribute to tumour aggressiveness.

Lymph node metastases (LNM) are an important prognostic factor for patients with intrahepatic cholangiocarcinoma, but underlying genetic alterations are poorly understood. Whole genome array comparative genomic hybridization (aCGH) was performed in 37 tumors and 14 matched LNM. Genomic analyses of tumors confirmed known and identified new (gains in 19q) copy number alterations (CNA). Tumors with LNM (N1) had more alterations and exclusive gains (3p, 4q, 5p, 13q) and losses (17p and 20p). LNM shared most alterations with their matched tumors (86%), but 79% acquired new isolated gains [12q14 (36%); 1p13, 2p23, 7p22, 7q11, 11q12, 13q13 and 14q12 (>20%)]. Unsupervised clustering revealed a poor prognosis subclass with increased alterations significantly associated to tumor differentiation and survival. TP53 and KRAS mutations occurred in 19% of tumors and 6% of metastases. Pathway analyses revealed association to cancer-associated pathways. Advanced tumor stage, microvascular/perineural invasion, and microscopic positive resection margin (R1) were significantly correlated to metastases, while N1-status, R1-resection, and poor tumor differentiation were significantly correlated to survival. ACGH identified clear differences between N0 (no LNM) and N1 tumors, while N1 tumors and matched LNM displayed high clonality with exclusive gains in the metastases. A novel subclass with increased CNAs and poor tumor differentiation was significantly correlated to survival.

Trait-associated loci often map to genomic regions encoding long noncoding RNAs (lncRNAs), but the role of these lncRNAs in disease etiology is largely unexplored. We show that a pair of sense/antisense lncRNA (6p22lncRNAs) encoded by CASC15 and NBAT1 located at the neuroblastoma (NB) risk-associated 6p22.3 locus are tumor suppressors and show reduced expression in high-risk NBs. Loss of functional synergy between 6p22lncRNAs results in an undifferentiated state that is maintained by a gene-regulatory network, including SOX9 located on 17q, a region frequently gained in NB. 6p22lncRNAs regulate SOX9 expression by controlling CHD7 stability via modulating the cellular localization of USP36, encoded by another 17q gene. This regulatory nexus between 6p22.3 and 17q regions may lead to potential NB treatment strategies.

Ferroptosis has emerged as an attractive strategy in cancer therapy. Understanding the operational networks regulating ferroptosis may unravel vulnerabilities that could be harnessed for therapeutic benefit. Using CRISPR-activation screens in ferroptosis hypersensitive cells, we identify the selenoprotein P (SELENOP) receptor, LRP8, as a key determinant protecting MYCN-amplified neuroblastoma cells from ferroptosis. Genetic deletion of LRP8 leads to ferroptosis as a result of an insufficient supply of selenocysteine, which is required for the translation of the antiferroptotic selenoprotein GPX4. This dependency is caused by low expression of alternative selenium uptake pathways such as system Xc- . The identification of LRP8 as a specific vulnerability of MYCN-amplified neuroblastoma cells was confirmed in constitutive and inducible LRP8 knockout orthotopic xenografts. These findings disclose a yet-unaccounted mechanism of selective ferroptosis induction that might be explored as a therapeutic strategy for high-risk neuroblastoma and potentially other MYCN-amplified entities.

The consistent focus of 'Advances in Neuroblastoma Research' congresses on the topic neuroblastoma sets it as a model for a circumscribed scientific community.

Telomere maintenance mechanisms (TMM) are a hallmark of high-risk neuroblastoma, and are conferred by activation of telomerase or alternative lengthening of telomeres (ALT). However, detection of TMM is not yet part of the clinical routine, and consensus on TMM detection, especially on ALT assessment, remains to be achieved.

TERT gene rearrangement with transcriptional superenhancers leads to TERT overexpression and neuroblastoma. No targeted therapy is available for clinical trials in patients with TERT-rearranged neuroblastoma.

In childhood cancer, the frequency of cancer-associated germline variants and their inheritance patterns are not thoroughly investigated. Moreover, the identification of children carrying a genetic predisposition by clinical means remains challenging. In this single-center study, we performed trio whole-exome sequencing and comprehensive clinical evaluation of a prospectively enrolled cohort of 160 children with cancer and their parents. We identified in 11/160 patients a pathogenic germline variant predisposing to cancer and a further eleven patients carried a prioritized VUS with a strong association to the cancerogenesis of the patient. Through clinical screening, 51 patients (31.3%) were identified as suspicious for an underlying cancer predisposition syndrome (CPS), but only in ten of those patients a pathogenic variant could be identified. In contrast, one patient with a classical CPS and ten patients with prioritized VUS were classified as unremarkable in the clinical work-up. Taken together, a monogenetic causative variant was detected in 13.8% of our patients using WES. Nevertheless, the still unclarified clinical suspicious cases emphasize the need to consider other genetic mechanisms including new target genes, structural variants, or polygenic interactions not previously associated with cancer predisposition.

Genomic instability can be a hallmark of both human genetic disease and cancer. We identify a deleterious UBQLN4 mutation in families with an autosomal recessive syndrome reminiscent of genome instability disorders. UBQLN4 deficiency leads to increased sensitivity to genotoxic stress and delayed DNA double-strand break (DSB) repair. The proteasomal shuttle factor UBQLN4 is phosphorylated by ATM and interacts with ubiquitylated MRE11 to mediate early steps of homologous recombination-mediated DSB repair (HRR). Loss of UBQLN4 leads to chromatin retention of MRE11, promoting non-physiological HRR activity in vitro and in vivo. Conversely, UBQLN4 overexpression represses HRR and favors non-homologous end joining. Moreover, we find UBQLN4 overexpressed in aggressive tumors. In line with an HRR defect in these tumors, UBQLN4 overexpression is associated with PARP1 inhibitor sensitivity. UBQLN4 therefore curtails HRR activity through removal of MRE11 from damaged chromatin and thus offers a therapeutic window for PARP1 inhibitor treatment in UBQLN4-overexpressing tumors.

Current risk stratification systems for neuroblastoma patients consider clinical, histopathological, and genetic variables, and additional prognostic markers have been proposed in recent years. We here sought to select highly informative covariates in a multistep strategy based on consecutive Cox regression models, resulting in a risk score that integrates hazard ratios of prognostic variables.

The prerequisite to prevent childhood B-cell acute lymphoblastic leukemia (B-ALL) is to decipher its etiology. The current model suggests that infection triggers B-ALL development through induction of activation-induced cytidine deaminase (AID; also known as AICDA) in precursor B-cells. This evidence has been largely acquired through the use of ex vivo functional studies. However, whether this mechanism governs native non-transplant B-ALL development is unknown. Here we show that, surprisingly, AID genetic deletion does not affect B-ALL development in Pax5-haploinsufficient mice prone to B-ALL upon natural infection exposure. We next test the effect of premature AID expression from earliest pro-B-cell stages in B-cell transformation. The generation of AID off-target mutagenic activity in precursor B-cells does not promote B-ALL. Likewise, known drivers of human B-ALL are not preferentially targeted by AID. Overall these results suggest that infections promote B-ALL through AID-independent mechanisms, providing evidence for a new model of childhood B-ALL development.

The discovery of frequent 8p11-p12 amplifications in squamous cell lung cancer (SQLC) has fueled hopes that FGFR1, located inside this amplicon, might be a therapeutic target. In a clinical trial, only 11% of patients with 8p11 amplification (detected by FISH) responded to FGFR kinase inhibitor treatment. To understand the mechanism of FGFR1 dependency, we performed deep genomic characterization of 52 SQLCs with 8p11-p12 amplification, including 10 tumors obtained from patients who had been treated with FGFR inhibitors. We discovered somatically altered variants of FGFR1 with deletion of exons 1-8 that resulted from intragenic tail-to-tail rearrangements. These ectodomain-deficient FGFR1 variants (ΔEC-FGFR1) were expressed in the affected tumors and were tumorigenic in both in vitro and in vivo models of lung cancer. Mechanistically, breakage-fusion-bridges were the source of 8p11-p12 amplification, resulting from frequent head-to-head and tail-to-tail rearrangements. Generally, tail-to-tail rearrangements within or in close proximity upstream of FGFR1 were associated with FGFR1 dependency. Thus, the genomic events shaping the architecture of the 8p11-p12 amplicon provide a mechanistic explanation for the emergence of FGFR1-driven SQLC. Specifically, we believe that FGFR1 ectodomain-deficient and FGFR1-centered amplifications caused by tail-to-tail rearrangements are a novel somatic genomic event that might be predictive of therapeutically relevant FGFR1 dependency.