Growing evidence suggests a key role for RNA binding proteins (RBPs) in genome stability programs. Additionally, recent developments in RNA sequencing technologies, as well as mass-spectrometry techniques, have greatly expanded our knowledge on protein-RNA interactions. We here use full transcriptome sequencing and label-free LC/MS/MS to identify global changes in protein-RNA interactions in response to etoposide-induced genotoxic stress. We show that RBPs have distinct binding patterns in response to genotoxic stress and that inactivation of the RBP regulator module, p38/MK2, can affect the entire spectrum of protein-RNA interactions that take place in response to stress. In addition to validating the role of known RBPs like Srsf1, Srsf2, Elavl1 in the genotoxic stress response, we add a new collection of RBPs to the DNA damage response. We identify Khsrp as a highly regulated RBP in response to genotoxic stress and further validate its role as a driver of the G(1/)S transition through the suppression of Cdkn1a(P21) transcripts. Finally, we identify KHSRP as an indicator of overall survival, as well as disease free survival in glioblastoma multiforme.

Gene expression profiling is being widely applied in cancer research to identify biomarkers for clinical endpoint prediction. Since RNA-seq provides a powerful tool for transcriptome-based applications beyond the limitations of microarrays, we sought to systematically evaluate the performance of RNA-seq-based and microarray-based classifiers in this MAQC-III/SEQC study for clinical endpoint prediction using neuroblastoma as a model.

Human protein-coding genes are often accompanied by divergently transcribed non-coding RNAs whose functions, especially in cell fate decisions, are poorly understood. Using an hESC-based cardiac differentiation model, we define a class of divergent lncRNAs, termed yin yang lncRNAs (yylncRNAs), that mirror the cell-type-specific expression pattern of their protein-coding counterparts. yylncRNAs are preferentially encoded from the genomic loci of key developmental cell fate regulators. Most yylncRNAs are spliced polyadenylated transcripts showing comparable expression patterns in vivo in mouse and in human embryos. Signifying their developmental function, the key mesoderm specifier BRACHYURY (T) is accompanied by yylncT, which localizes to the active T locus during mesoderm commitment. yylncT binds the de novo DNA methyltransferase DNMT3B, and its transcript is required for activation of the T locus, with yylncT depletion specifically abolishing mesodermal commitment. Collectively, we report a lncRNA-mediated regulatory layer safeguarding embryonic cell fate transitions.

Kinase inhibitors represent the backbone of targeted cancer therapy, yet only a limited number of oncogenic drivers are directly druggable. By interrogating the activity of 1,505 kinase inhibitors, we found that BRD4-NUT-rearranged NUT midline carcinoma (NMC) cells are specifically killed by CDK9 inhibition (CDK9i) and depend on CDK9 and Cyclin-T1 expression. We show that CDK9i leads to robust induction of apoptosis and of markers of DNA damage response in NMC cells. While both CDK9i and bromodomain inhibition over time result in reduced Myc protein expression, only bromodomain inhibition induces cell differentiation and a p21-induced cell-cycle arrest in these cells. Finally, RNA-seq and ChIP-based analyses reveal a BRD4-NUT-specific CDK9i-induced perturbation of transcriptional elongation. Thus, our data provide a mechanistic basis for the genotype-dependent vulnerability of NMC cells to CDK9i that may be of relevance for the development of targeted therapies for NMC patients.

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.

Despite the successful introduction of potent anti-cancer therapeutics, most of these drugs lead to only modest tumor-shrinkage or transient responses, followed by re-growth of tumors. Combining different compounds has resulted in enhanced tumor control and prolonged survival. However, methods querying the efficacy of such combinations have been hampered by limited scalability, analytical resolution, statistical feasibility, or a combination thereof. We have developed a theoretical framework modeling cellular viability as a stochastic lifetime process to determine synergistic compound combinations from high-throughput cellular screens. We apply our method to data derived from chemical perturbations of 65 cancer cell lines with two inhibitors. Our analysis revealed synergy for the combination of both compounds in subsets of cell lines. By contrast, in cell lines in which inhibition of one of both targets was sufficient to induce cell death, no synergy was detected, compatible with the topology of the oncogenically activated signaling network. In summary, we provide a tool for the measurement of synergy strength for combination perturbation experiments that might help define pathway topologies and direct clinical trials.

The IGF receptor type 1 (IGF-1R) pathway is frequently deregulated in human tumors and has become a target of interest for anti-cancer therapy.

Deciphering the evolution of cancer cells under therapeutic pressure is a crucial step to understand the mechanisms that lead to treatment resistance. To this end, we analyzed whole-exome sequencing data of eight chronic lymphocytic leukemia (CLL) patients that developed resistance upon BCL2-inhibition by venetoclax. Here, we report recurrent mutations in BTG1 (2 patients) and homozygous deletions affecting CDKN2A/B (3 patients) that developed during treatment, as well as a mutation in BRAF and a high-level focal amplification of CD274 (PD-L1) that might pinpoint molecular aberrations offering structures for further therapeutic interventions.

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.

High resolution molecular studies have demonstrated that the clonal acquisition of gene mutations is an important mechanism that may promote rapid disease progression and drug resistance in chronic lymphocytic leukemia (CLL). Therefore, the early and sensitive detection of such mutations is an important prerequisite for future predictive CLL diagnostics in the clinical setting.

Genomic translocation events frequently underlie cancer development through generation of gene fusions with oncogenic properties. Identification of such fusion transcripts by transcriptome sequencing might help to discover new potential therapeutic targets. We developed TRUP (Tumor-specimen suited RNA-seq Unified Pipeline) (https://github.com/ruping/TRUP), a computational approach that combines split-read and read-pair analysis with de novo assembly for the identification of chimeric transcripts in cancer specimens. We apply TRUP to RNA-seq data of different tumor types, and find it to be more sensitive than alternative tools in detecting chimeric transcripts, such as secondary rearrangements in EML4-ALK-positive lung tumors, or recurrent inactivating rearrangements affecting RASSF8.

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.

Loss of the tumor suppressors RB1 and TP53 and MYC amplification are frequent oncogenic events in small cell lung cancer (SCLC). We show that Myc expression cooperates with Rb1 and Trp53 loss in the mouse lung to promote aggressive, highly metastatic tumors, that are initially sensitive to chemotherapy followed by relapse, similar to human SCLC. Importantly, MYC drives a neuroendocrine-low "variant" subset of SCLC with high NEUROD1 expression corresponding to transcriptional profiles of human SCLC. Targeted drug screening reveals that SCLC with high MYC expression is vulnerable to Aurora kinase inhibition, which, combined with chemotherapy, strongly suppresses tumor progression and increases survival. These data identify molecular features for patient stratification and uncover a potential targeted treatment approach for MYC-driven SCLC.

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.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.