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.

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.

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.

Epigenetic regulation of neuroblastoma development.

In recent years, technological advances have enabled a detailed landscaping of the epigenome and the mechanisms of epigenetic regulation that drive normal cell function, development and cancer. Rather than merely a structural entity to support genome compaction, we now look at chromatin as a very dynamic and essential constellation that is actively participating in the tight orchestration of transcriptional regulation as well as DNA replication and repair. The unique feature of chromatin flexibility enabling fast switches towards more or less restricted epigenetic cellular states is, not surprisingly, intimately connected to cancer development and treatment resistance, and the central role of epigenetic alterations in cancer is illustrated by the finding that up to 50% of all mutations across cancer entities affect proteins controlling the chromatin status. We summarize recent insights into epigenetic rewiring underlying neuroblastoma (NB) tumor formation ranging from changes in DNA methylation patterns and mutations in epigenetic regulators to global effects on transcriptional regulatory circuits that involve key players in NB oncogenesis. Insights into the disruption of the homeostatic epigenetic balance contributing to developmental arrest of sympathetic progenitor cells and subsequent NB oncogenesis are rapidly growing and will be exploited towards the development of novel therapeutic strategies to increase current survival rates of patients with high-risk NB.

Correction to: Epigenetic regulation of neuroblastoma development.

The authorship names of this paper are incorrect.

Long non-coding RNAs in leukemia: biology and clinical impact.

PURPOSE OF REVIEW: Over the last years, long non-coding RNAs (lncRNAs) have emerged as putative regulators of malignant hematopoietic development. Here, we review recent literature on the involvement of lncRNAs in leukemia, including their role in driving or sustaining disease and their potential impact on diagnosis, classification, and prognosis. RECENT FINDINGS: Leukemogenesis is a complex process resulting from the accumulation of multiple genetic alterations. Over the last years, advances in high-throughput sequencing and transcriptome profiling have enabled the identification of lncRNAs involved in leukemia development. lncRNAs are able to distinguish different subtypes of human leukemia and several reports have identified specific patterns of lncRNA expression associated with clinical patient characteristics. Although functional studies on the actual role of these lncRNAs during transformation remain scarce, emerging evidence suggests that complex interactions between coding and non-coding transcript are truly involved in leukemia development. SUMMARY: Introduction of lncRNAs as an additional layer of complexity in human leukemia might provide new molecular genetic insights in the biology of this disease and could create unique opportunities for the identification of novel drug targets and diagnostic or prognostic biomarkers.

The H3K27me3 demethylase UTX is a gender-specific tumor suppressor in T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive form of leukemia that is mainly diagnosed in children and shows a skewed gender distribution toward males. In this study, we report somatic loss-of-function mutations in the X-linked histone H3K27me3 demethylase ubiquitously transcribed X (UTX) chromosome, in human T-ALL. Interestingly, UTX mutations were exclusively present in male T-ALL patients and allelic expression analysis revealed that UTX escapes X-inactivation in female T-ALL lymphoblasts and normal T cells. Notably, we demonstrate in vitro and in vivo that the H3K27me3 demethylase UTX functions as a bona fide tumor suppressor in T-ALL. Moreover, T-ALL driven by UTX inactivation exhibits collateral sensitivity to pharmacologic H3K27me3 inhibition. All together, our results show how a gender-specific and therapeutically relevant defect in balancing H3K27 methylation contributes to T-cell leukemogenesis.

The need for transparency and good practices in the qPCR literature.

Two surveys of over 1,700 publications whose authors use quantitative real-time PCR (qPCR) reveal a lack of transparent and comprehensive reporting of essential technical information. Reporting standards are significantly improved in publications that cite the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines, although such publications are still vastly outnumbered by those that do not.

A nanobody modulates the p53 transcriptional program without perturbing its functional architecture.

The p53 transcription factor plays an important role in genome integrity. To perform this task, p53 regulates the transcription of genes promoting various cellular outcomes including cell cycle arrest, apoptosis or senescence. The precise regulation of this activity remains elusive as numerous mechanisms, e.g. posttranslational modifications of p53 and (non-)covalent p53 binding partners, influence the p53 transcriptional program. We developed a novel, non-invasive tool to manipulate endogenous p53. Nanobodies (Nb), raised against the DNA-binding domain of p53, allow us to distinctively target both wild type and mutant p53 with great specificity. Nb3 preferentially binds 'structural' mutant p53, i.e. R175H and R282W, while a second but distinct nanobody, Nb139, binds both mutant and wild type p53. The co-crystal structure of the p53 DNA-binding domain in complex with Nb139 (1.9 Å resolution) reveals that Nb139 binds opposite the DNA-binding surface. Furthermore, we demonstrate that Nb139 does not disturb the functional architecture of the p53 DNA-binding domain using conformation-specific p53 antibody immunoprecipitations, glutaraldehyde crosslinking assays and chromatin immunoprecipitation. Functionally, the binding of Nb139 to p53 allows us to perturb the transactivation of p53 target genes. We propose that reduced recruitment of transcriptional co-activators or modulation of selected post-transcriptional modifications account for these observations.

Comprehensive analysis of transcriptome variation uncovers known and novel driver events in T-cell acute lymphoblastic leukemia.

RNA-seq is a promising technology to re-sequence protein coding genes for the identification of single nucleotide variants (SNV), while simultaneously obtaining information on structural variations and gene expression perturbations. We asked whether RNA-seq is suitable for the detection of driver mutations in T-cell acute lymphoblastic leukemia (T-ALL). These leukemias are caused by a combination of gene fusions, over-expression of transcription factors and cooperative point mutations in oncogenes and tumor suppressor genes. We analyzed 31 T-ALL patient samples and 18 T-ALL cell lines by high-coverage paired-end RNA-seq. First, we optimized the detection of SNVs in RNA-seq data by comparing the results with exome re-sequencing data. We identified known driver genes with recurrent protein altering variations, as well as several new candidates including H3F3A, PTK2B, and STAT5B. Next, we determined accurate gene expression levels from the RNA-seq data through normalizations and batch effect removal, and used these to classify patients into T-ALL subtypes. Finally, we detected gene fusions, of which several can explain the over-expression of key driver genes such as TLX1, PLAG1, LMO1, or NKX2-1; and others result in novel fusion transcripts encoding activated kinases (SSBP2-FER and TPM3-JAK2) or involving MLLT10. In conclusion, we present novel analysis pipelines for variant calling, variant filtering, and expression normalization on RNA-seq data, and successfully applied these for the detection of translocations, point mutations, INDELs, exon-skipping events, and expression perturbations in T-ALL.

The Notch driven long non-coding RNA repertoire in T-cell acute lymphoblastic leukemia.

Genetic studies in T-cell acute lymphoblastic leukemia have uncovered a remarkable complexity of oncogenic and loss-of-function mutations. Amongst this plethora of genetic changes, NOTCH1 activating mutations stand out as the most frequently occurring genetic defect, identified in more than 50% of T-cell acute lymphoblastic leukemias, supporting a role as an essential driver for this gene in T-cell acute lymphoblastic leukemia oncogenesis. In this study, we aimed to establish a comprehensive compendium of the long non-coding RNA transcriptome under control of Notch signaling. For this purpose, we measured the transcriptional response of all protein coding genes and long non-coding RNAs upon pharmacological Notch inhibition in the human T-cell acute lymphoblastic leukemia cell line CUTLL1 using RNA-sequencing. Similar Notch dependent profiles were established for normal human CD34(+) thymic T-cell progenitors exposed to Notch signaling activity in vivo. In addition, we generated long non-coding RNA expression profiles (array data) from ex vivo isolated Notch active CD34(+) and Notch inactive CD4(+)CD8(+) thymocytes and from a primary cohort of 15 T-cell acute lymphoblastic leukemia patients with known NOTCH1 mutation status. Integration of these expression datasets with publicly available Notch1 ChIP-sequencing data resulted in the identification of long non-coding RNAs directly regulated by Notch activity in normal and malignant T cells. Given the central role of Notch in T-cell acute lymphoblastic leukemia oncogenesis, these data pave the way for the development of novel therapeutic strategies that target hyperactive Notch signaling in human T-cell acute lymphoblastic leukemia.

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.

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.

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.