HTSplotter: An end-to-end data processing, analysis and visualisation tool for chemical and genetic in vitro perturbation screening.

In biomedical research, high-throughput screening is often applied as it comes with automatization, higher-efficiency, and more and faster results. High-throughput screening experiments encompass drug, drug combination, genetic perturbagen or a combination of genetic and chemical perturbagen screens. These experiments are conducted in real-time assays over time or in an endpoint assay. The data analysis consists of data cleaning and structuring, as well as further data processing and visualisation, which, due to the amount of data, can easily become laborious, time-consuming and error-prone. Therefore, several tools have been developed to aid researchers in this process, but these typically focus on specific experimental set-ups and are unable to process data of several time points and genetic-chemical perturbagen screens. To meet these needs, we developed HTSplotter, a web tool and Python module that performs automatic data analysis and visualization of visualization of eitherendpoint or real-time assays from different high-throughput screening experiments: drug, drug combination, genetic perturbagen and genetic-chemical perturbagen screens. HTSplotter implements an algorithm based on conditional statements to identify experiment types and controls. After appropriate data normalization, including growth rate normalization, HTSplotter executes downstream analyses such as dose-response relationship and drug synergism assessment by the Bliss independence (BI), Zero Interaction Potency (ZIP) and Highest Single Agent (HSA) methods. All results are exported as a text file and plots are saved in a PDF file. The main advantage of HTSplotter over other available tools is the automatic analysis of genetic-chemical perturbagen screens and real-time assays where growth rate and perturbagen effect results are plotted over time. In conclusion, HTSplotter allows for the automatic end-to-end data processing, analysis and visualisation of various high-throughput in vitro cell culture screens, offering major improvements in terms of versatility, efficiency and time over existing tools.

Targeted AURKA degradation: Towards new therapeutic agents for neuroblastoma.

Aurora kinase A (AURKA) is a well-established target in neuroblastoma (NB) due to both its catalytic functions during mitosis and its kinase-independent functions, including stabilization of the key oncoprotein MYCN. We present a structure-activity relationship (SAR) study of MK-5108-derived PROTACs against AURKA by exploring different linker lengths and exit vectors on the thalidomide moiety. PROTAC SK2188 induces the most potent AURKA degradation (DC50,24h 3.9 nM, Dmax,24h 89%) and shows an excellent binding and degradation selectivity profile. Treatment of NGP neuroblastoma cells with SK2188 induced concomitant MYCN degradation, high replication stress/DNA damage levels and apoptosis. Moreover, SK2188 significantly outperforms the parent inhibitor MK-5108 in a cell proliferation screen and patient-derived organoids. Furthermore, altering the attachment point of the PEG linker to the 5-position of thalidomide allowed us to identify a potent AURKA degrader with a linker as short as 2 PEG units. With this, our SAR-study provides interesting lead structures for further optimization and validation of AURKA degradation as a potential therapeutic strategy in neuroblastoma.

Whole transcriptome profiling of liquid biopsies from tumour xenografted mouse models enables specific monitoring of tumour-derived extracellular RNA.

While cell-free DNA (cfDNA) is widely being investigated, free circulating RNA (extracellular RNA, exRNA) has the potential to improve cancer therapy response monitoring and detection due to its dynamic nature. However, it remains unclear in which blood subcompartment tumour-derived exRNAs primarily reside. We developed a host-xenograft deconvolution framework, exRNAxeno, with mapping strategies to either a combined human-mouse reference genome or both species genomes in parallel, applicable to exRNA sequencing data from liquid biopsies of human xenograft mouse models. The tool enables to distinguish (human) tumoural RNA from (murine) host RNA, to specifically analyse tumour-derived exRNA. We applied the combined pipeline to total exRNA sequencing data from 95 blood-derived liquid biopsy samples from 30 mice, xenografted with 11 different tumours. Tumoural exRNA concentrations are not determined by plasma platelet levels, while host exRNA concentrations increase with platelet content. Furthermore, a large variability in exRNA abundance and transcript content across individual mice is observed. The tumoural gene detectability in plasma is largely correlated with the RNA expression levels in the tumour tissue or cell line. These findings unravel new aspects of tumour-derived exRNA biology in xenograft models and open new avenues to further investigate the role of exRNA in cancer.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition.

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.

From DNA Copy Number Gains and Tumor Dependencies to Novel Therapeutic Targets for High-Risk Neuroblastoma.

Neuroblastoma is a pediatric tumor arising from the sympatho-adrenal lineage and a worldwide leading cause of childhood cancer-related deaths. About half of high-risk patients die from the disease while survivors suffer from multiple therapy-related side-effects. While neuroblastomas present with a low mutational burden, focal and large segmental DNA copy number aberrations are highly recurrent and associated with poor survival. It can be assumed that the affected chromosomal regions contain critical genes implicated in neuroblastoma biology and behavior. More specifically, evidence has emerged that several of these genes are implicated in tumor dependencies thus potentially providing novel therapeutic entry points. In this review, we briefly review the current status of recurrent DNA copy number aberrations in neuroblastoma and provide an overview of the genes affected by these genomic variants for which a direct role in neuroblastoma has been established. Several of these genes are implicated in networks that positively regulate MYCN expression or stability as well as cell cycle control and apoptosis. Finally, we summarize alternative approaches to identify and prioritize candidate copy-number driven dependency genes for neuroblastoma offering novel therapeutic opportunities.

MEIS2 Is an Adrenergic Core Regulatory Transcription Factor Involved in Early Initiation of TH-MYCN-Driven Neuroblastoma Formation.

Roughly half of all high-risk neuroblastoma patients present with MYCN amplification. The molecular consequences of MYCN overexpression in this aggressive pediatric tumor have been studied for decades, but thus far, our understanding of the early initiating steps of MYCN-driven tumor formation is still enigmatic. We performed a detailed transcriptome landscaping during murine TH-MYCN-driven neuroblastoma tumor formation at different time points. The neuroblastoma dependency factor MEIS2, together with ASCL1, was identified as a candidate tumor-initiating factor and shown to be a novel core regulatory circuit member in adrenergic neuroblastomas. Of further interest, we found a KEOPS complex member (gm6890), implicated in homologous double-strand break repair and telomere maintenance, to be strongly upregulated during tumor formation, as well as the checkpoint adaptor Claspin (CLSPN) and three chromosome 17q loci CBX2, GJC1 and LIMD2. Finally, cross-species master regulator analysis identified FOXM1, together with additional hubs controlling transcriptome profiles of MYCN-driven neuroblastoma. In conclusion, time-resolved transcriptome analysis of early hyperplastic lesions and full-blown MYCN-driven neuroblastomas yielded novel components implicated in both tumor initiation and maintenance, providing putative novel drug targets for MYCN-driven neuroblastoma.

MYCN-induced nucleolar stress drives an early senescence-like transcriptional program in hTERT-immortalized RPE cells.

MYCN is an oncogenic driver in neural crest-derived neuroblastoma and medulloblastoma. To better understand the early effects of MYCN activation in a neural-crest lineage context, we profiled the transcriptome of immortalized human retina pigment epithelial cells with inducible MYCN activation. Gene signatures associated with elevated MYC/MYCN activity were induced after 24 h of MYCN activation, which attenuated but sustained at later time points. Unexpectedly, MYCN activation was accompanied by reduced cell growth. Gene set enrichment analysis revealed a senescence-like signature with strong induction of p53 and p21 but in the absence of canonical hallmarks of senescence such as ß-galactosidase positivity, suggesting incomplete cell fate commitment. When scrutinizing the putative drivers of this growth attenuation, differential gene expression analysis identified several regulators of nucleolar stress. This process was also reflected by phenotypic correlates such as cytoplasmic granule accrual and nucleolar coalescence. Hence, we propose that the induction of MYCN congests the translational machinery, causing nucleolar stress and driving cells into a transient pre-senescent state. Our findings shed new light on the early events induced by MYCN activation and may help unravelling which factors are required for cells to tolerate unscheduled MYCN overexpression during early malignant transformation.

Recurrent chromosomal imbalances provide selective advantage to human embryonic stem cells under enhanced replicative stress conditions.

Human embryonic stem cells (hESCs) and embryonal tumors share a number of common features, including a compromised G1/S checkpoint. Consequently, these rapidly dividing hESCs and cancer cells undergo elevated levels of replicative stress, inducing genomic instability that drives chromosomal imbalances. In this context, it is of interest that long-term in vitro cultured hESCs exhibit a remarkable high incidence of segmental DNA copy number gains, some of which are also highly recurrent in certain malignancies such as 17q gain (17q+). The selective advantage of DNA copy number changes in these cells has been attributed to several underlying processes including enhanced proliferation. We hypothesized that these recurrent chromosomal imbalances become rapidly embedded in the cultured hESCs through a replicative stress driven Darwinian selection process. To this end, we compared the effect of hydroxyurea-induced replicative stress vs normal growth conditions in an equally mixed cell population of isogenic euploid and 17q + hESCs. We could show that 17q + hESCs rapidly overtook normal hESCs. Our data suggest that recurrent chromosomal segmental gains provide a proliferative advantage to hESCs under increased replicative stress, a process that may also explain the highly recurrent nature of certain imbalances in cancer.

Large-scale circular RNA deregulation in T-ALL: unlocking unique ectopic expression of molecular subtypes.

Circular RNAs (circRNAs) are stable RNA molecules that can drive cancer through interactions with microRNAs and proteins and by the expression of circRNA encoded peptides. The aim of the study was to define the circRNA landscape and potential impact in T-cell acute lymphoblastic leukemia (T-ALL). Analysis by CirComPara of RNA-sequencing data from 25 T-ALL patients, immature, HOXA overexpressing, TLX1, TLX3, TAL1, or LMO2 rearranged, and from thymocyte populations of human healthy donors disclosed 68 554 circRNAs. Study of the top 3447 highly expressed circRNAs identified 944 circRNAs with significant differential expression between malignant T cells and normal counterparts, with most circRNAs displaying increased expression in T-ALL. Next, we defined subtype-specific circRNA signatures in molecular genetic subgroups of human T-ALL. In particular, circZNF609, circPSEN1, circKPNA5, and circCEP70 were upregulated in immature, circTASP1, circZBTB44, and circBACH1 in TLX3, circHACD1, and circSTAM in HOXA, circCAMSAP1 in TLX1, and circCASC15 in TAL-LMO. Backsplice sequences of 14 circRNAs ectopically expressed in T-ALL were confirmed, and overexpression of circRNAs in T-ALL with specific oncogenic lesions was substantiated by quantification in a panel of 13 human cell lines. An oncogenic role of circZNF609 in T-ALL was indicated by decreased cell viability upon silencing in vitro. Furthermore, functional predictions identified circRNA-microRNA gene axes informing modes of circRNA impact in molecular subtypes of human T-ALL.

PHF6 Expression Levels Impact Human Hematopoietic Stem Cell Differentiation.

Transcriptional control of hematopoiesis involves complex regulatory networks and functional perturbations in one of these components often results in malignancies. Loss-of-function mutations in PHF6, encoding a presumed epigenetic regulator, have been primarily described in T cell acute lymphoblastic leukemia (T-ALL) and the first insights into its function in normal hematopoiesis only recently emerged from mouse modeling experiments. Here, we investigated the role of PHF6 in human blood cell development by performing knockdown studies in cord blood and thymus-derived hematopoietic precursors to evaluate the impact on lineage differentiation in well-established in vitro models. Our findings reveal that PHF6 levels differentially impact the differentiation of human hematopoietic progenitor cells into various blood cell lineages, with prominent effects on lymphoid and erythroid differentiation. We show that loss of PHF6 results in accelerated human T cell development through reduced expression of NOTCH1 and its downstream target genes. This functional interaction in developing thymocytes was confirmed in vivo using a phf6-deficient zebrafish model that also displayed accelerated developmental kinetics upon reduced phf6 or notch1 activation. In summary, our work reveals that appropriate control of PHF6 expression is important for normal human hematopoiesis and provides clues towards the role of PHF6 in T-ALL development.

Distinct Notch1 and BCL11B requirements mediate human γδ/αß T cell development.

γδ and αß T cells have unique roles in immunity and both originate in the thymus from T-lineage committed precursors through distinct but unclear mechanisms. Here, we show that Notch1 activation is more stringently required for human γδ development compared to αß-lineage differentiation and performed paired mRNA and miRNA profiling across 11 discrete developmental stages of human T cell development in an effort to identify the potential Notch1 downstream mechanism. Our data suggest that the miR-17-92 cluster is a Notch1 target in immature thymocytes and that miR-17 can restrict BCL11B expression in these Notch-dependent T cell precursors. We show that enforced miR-17 expression promotes human γδ T cell development and, consistently, that BCL11B is absolutely required for αß but less for γδ T cell development. This study suggests that human γδ T cell development is mediated by a stage-specific Notch-driven negative feedback loop through which miR-17 temporally restricts BCL11B expression and provides functional insights into the developmental role of the disease-associated genes BCL11B and the miR-17-92 cluster in a human context.

The ETS transcription factor ETV5 is a target of activated ALK in neuroblastoma contributing to increased tumour aggressiveness.

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.

Publisher Correction: In silico discovery of a FOXM1 driven embryonal signaling pathway in therapy resistant neuroblastoma tumors.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Purification of high-quality RNA from a small number of fluorescence activated cell sorted zebrafish cells for RNA sequencing purposes.

BACKGROUND: Transgenic zebrafish lines with the expression of a fluorescent reporter under the control of a cell-type specific promoter, enable transcriptome analysis of FACS sorted cell populations. RNA quality and yield are key determinant factors for accurate expression profiling. Limited cell number and FACS induced cellular stress make RNA isolation of sorted zebrafish cells a delicate process. We aimed to optimize a workflow to extract sufficient amounts of high-quality RNA from a limited number of FACS sorted cells from Tg(fli1a:GFP) zebrafish embryos, which can be used for accurate gene expression analysis. RESULTS: We evaluated two suitable RNA isolation kits (the RNAqueous micro and the RNeasy plus micro kit) and determined that sorting cells directly into lysis buffer is a critical step for success. For low cell numbers, this ensures direct cell lysis, protects RNA from degradation and results in a higher RNA quality and yield. We showed that this works well up to 0.5× dilution of the lysis buffer with sorted cells. In our sort settings, this corresponded to 30,000 and 75,000 cells for the RNAqueous micro kit and RNeasy plus micro kit respectively. Sorting more cells dilutes the lysis buffer too much and requires the use of a collection buffer. We also demonstrated that an additional genomic DNA removal step after RNA isolation is required to completely clear the RNA from any contaminating genomic DNA. For cDNA synthesis and library preparation, we combined SmartSeq v4 full length cDNA library amplification, Nextera XT tagmentation and sample barcoding. Using this workflow, we were able to generate highly reproducible RNA sequencing results. CONCLUSIONS: The presented optimized workflow enables to generate high quality RNA and allows accurate transcriptome profiling of small populations of sorted zebrafish cells.

ALK positively regulates MYCN activity through repression of HBP1 expression.

ALK mutations occur in 10% of primary neuroblastomas and represent a major target for precision treatment. In combination with MYCN amplification, ALK mutations infer an ultra-high-risk phenotype resulting in very poor patient prognosis. To open up opportunities for future precision drugging, a deeper understanding of the molecular consequences of constitutive ALK signaling and its relationship to MYCN activity in this aggressive pediatric tumor entity will be essential. We show that mutant ALK downregulates the 'HMG-box transcription factor 1' (HBP1) through the PI3K-AKT-FOXO3a signaling axis. HBP1 inhibits both the transcriptional activating and repressing activity of MYCN, the latter being mediated through PRC2 activity. HBP1 itself is under negative control of MYCN through miR-17~92. Combined targeting of HBP1 by PI3K antagonists and MYCN signaling by BET- or HDAC-inhibitors blocks MYCN activity and significantly reduces tumor growth, suggesting a novel targeted therapy option for high-risk neuroblastoma.

In silico discovery of a FOXM1 driven embryonal signaling pathway in therapy resistant neuroblastoma tumors.

Chemotherapy resistance is responsible for high mortality rates in neuroblastoma. MYCN, an oncogenic driver in neuroblastoma, controls pluripotency genes including LIN28B. We hypothesized that enhanced embryonic stem cell (ESC) gene regulatory programs could mark tumors with high pluripotency capacity and subsequently increased risk for therapy failure. An ESC miRNA signature was established based on publicly available data. In addition, an ESC mRNA signature was generated including the 500 protein coding genes with the highest positive expression correlation with the ESC miRNA signature score in 200 neuroblastomas. High ESC m(i)RNA expression signature scores were significantly correlated with poor neuroblastoma patient outcome specifically in the subgroup of MYCN amplified tumors and stage 4 nonamplified tumors. Further data-mining identified FOXM1, as the major predicted driver of this ESC signature, controlling a large set of genes implicated in cell cycle control and DNA damage response. Of further interest, re-analysis of published data showed that MYCN transcriptionally activates FOXM1 in neuroblastoma cells. In conclusion, a novel ESC m(i)RNA signature stratifies neuroblastomas with poor prognosis, enabling the identification of therapy-resistant tumors. The finding that this signature is strongly FOXM1 driven, warrants for drug design targeted at FOXM1 or key components controlling this pathway.

TBX2 is a neuroblastoma core regulatory circuitry component enhancing MYCN/FOXM1 reactivation of DREAM targets.

Chromosome 17q gains are almost invariably present in high-risk neuroblastoma cases. Here, we perform an integrative epigenomics search for dosage-sensitive transcription factors on 17q marked by H3K27ac defined super-enhancers and identify TBX2 as top candidate gene. We show that TBX2 is a constituent of the recently established core regulatory circuitry in neuroblastoma with features of a cell identity transcription factor, driving proliferation through activation of p21-DREAM repressed FOXM1 target genes. Combined MYCN/TBX2 knockdown enforces cell growth arrest suggesting that TBX2 enhances MYCN sustained activation of FOXM1 targets. Targeting transcriptional addiction by combined CDK7 and BET bromodomain inhibition shows synergistic effects on cell viability with strong repressive effects on CRC gene expression and p53 pathway response as well as several genes implicated in transcriptional regulation. In conclusion, we provide insight into the role of the TBX2 CRC gene in transcriptional dependency of neuroblastoma cells warranting clinical trials using BET and CDK7 inhibitors.