Search

Summary: Neuroblastoma, a rare embryonic tumor arising from neural crest development, is responsible for 15% of pediatric cancer-related deaths. Recently, several single-cell transcriptome studies were performed on neuroblastoma patient samples to investigate the cell of origin and tumor heterogeneity. However, these individual studies involved a small number of tumors and cells, limiting the conclusions that could be drawn. To overcome this limitation, we integrated seven single-cell or single-nucleus datasets into a harmonized cell atlas covering 362,991 cells across 61 patients. We use this atlas to decipher the transcriptional landscape of neuroblastoma at single-cell resolution, revealing associations between transcriptomic profiles and clinical outcomes within the tumor compartment. In addition, we characterize the complex immune-cell landscape and uncover considerable heterogeneity among tumor-associated macrophages. Finally, we showcase the utility of our atlas as a resource by expanding it with additional data and using it as a reference for data-driven cell-type annotation.

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.

Summary: Studies defining normal and disrupted human neural crest cell development have been challenging given its early timing and intricacy of development. Consequently, insight into the early disruptive events causing a neural crest related disease such as pediatric cancer neuroblastoma is limited. To overcome this problem, we developed an in vitro differentiation model to recapitulate the normal in vivo developmental process of the sympathoadrenal lineage which gives rise to neuroblastoma. We used human in vitro pluripotent stem cells and single-cell RNA sequencing to recapitulate the molecular events during sympathoadrenal development. We provide a detailed map of dynamically regulated transcriptomes during sympathoblast formation and illustrate the power of this model to study early events of the development of human neuroblastoma, identifying a distinct subpopulation of cell marked by SOX2 expression in developing sympathoblast obtained from patient derived iPSC cells harboring a germline activating mutation in the anaplastic lymphoma kinase (ALK) gene.

The development of neuroblastoma (NB) is regulated by multiple core transcription factors. Here, SOX11 is identified as a potential epigenetic master regulator upstream of the core regulatory circuitry in adrenergic high-risk neuroblastoma.

Abstract 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.

Abstract 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.

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.

In high-risk neuroblastoma cases, gains in chromosome 17q are common. Here, the authors investigate the epigenomics and transcriptomics of neuroblastoma, identifying TBX2 as a core regulatory circuitry component enhancing the reactivation of DREAM targets by MYCN/FOXM1.

Chromosome 17q gains are a common alteration in high-risk neuroblastomas with unknown functional significance. We identified a 17q super-enhancer regulated T-box Transcription Factor 2 (TBX2) as constituent of a core regulatory circuitry driving proliferation through enhancing V-myc myelocytomatosis viral-related oncogene, neuroblastoma derived (avian) (MYCN)/Forkhead box protein M1(FOXM1) reactivation of dimerization partner, RB-like, E2F and multi-vulval class B (DREAM) targets, which can be affected synergistically by combined cyclin-dependent kinase 7 and Bromo-domain inhibition.

Abstract T-cell acute lymphoblastic leukemia (T-ALL) is a genetically heterogeneous disease that can be classified into different molecular genetic subtypes according to their mRNA gene expression profile. In this study, we applied RNA sequencing to investigate the full spectrum of miRNA expression in primary T-ALL patient samples, T-ALL leukemia cell lines and healthy donor thymocytes. Notably, this analysis revealed that genetic subtypes of human T-ALL also display unique miRNA expression signatures, which are largely conserved in human T-ALL cell lines with corresponding genetic background. Furthermore, small RNA-sequencing also unraveled the variety of isoforms that are expressed for each miRNA in T-ALL and showed that a significant number of miRNAs are actually represented by an alternative isomiR. Finally, comparison of CD34+ and CD4+CD8+ healthy donor thymocytes and T-ALL miRNA profiles allowed identifying several novel miRNAs with putative oncogenic or tumor suppressor functions in T-ALL. Altogether, this study provides a comprehensive overview of miRNA expression in normal and malignant T-cells and sets the stage for functional evaluation of novel miRNAs in T-ALL disease biology.

Abstract In the last decade, the role for noncoding RNAs in disease was clearly established, starting with microRNAs and later expanded towards long noncoding RNAs. This was also the case for T cell acute lymphoblastic leukemia, which is a malignant blood disorder arising from oncogenic events during normal T cell development in the thymus. By studying the transcriptomic profile of protein-coding genes, several oncogenic events leading to T cell acute lymphoblastic leukemia (T-ALL) could be identified. In recent years, it became apparent that several of these oncogenes function via microRNAs and long noncoding RNAs. In this review, we give a detailed overview of the studies that describe the noncoding RNAome in T-ALL oncogenesis and normal T cell development.

Strong Notch signalling promotes initial T cell lineage specification of lymphoid progenitors but is also permissive for thymic natural killer (NK) cell development. Here the authors show that GATA3 directs human T-lineage commitment by modulating Notch activity and repressing the NK programme.

Despite the established contribution of deregulated microRNA (miRNA) function to carcinogenesis, relatively few miRNA-cancer gene interactions have been validated, making it difficult to appreciate the true complexity of miRNA-cancer gene regulatory networks.In this effort, we identify miRNA interactomes of 17 well-established cancer genes, involved in various cancer types, through a miRNome-wide 3' UTR reporter screening. Using a novel and performant strategy for high-throughput screening data analysis, we identify 390 interactions, quadrupling the size of the known miRNA interactome for the cancer genes under investigation. Clear enrichments of established and predicted interactions underscore the validity of the interactome data set. Interactomes appear to be primarily driven by canonical binding site interactions. Nonetheless, non-canonical binding sites, such as offset 6mer and seed-mismatched or G:U wobble sites, also have regulatory activity, albeit clearly less pronounced. Furthermore, we observe enhanced regulation in the presence of 3' supplementary pairing for both canonical and non-canonical binding sites.Altogether, the cancer gene-miRNA interactome data set represents a unique resource that will aid in the unraveling of regulatory miRNA networks and the dynamic regulation of key protein-coding cancer genes. In addition, it uncovers aspects of the functional miRNA binding site's architecture and the relative contributions of different binding site types.

Effective inhibitors of invasion and metastasis represent a serious unmet clinical need. We have recently identified 4-fluoro-3',4',5'-trimethoxychalcone or C16 as a potent anti-invasive molecule. In this paper, we report on the development of an optimized vehicle for oral administration of C16. We also explore its pharmacokinetic and toxicity profile in rodents as a prelude to a broad-scope evaluation as a pharmacological tool in animal models of disease. C16 showed suboptimal pharmacokinetics with limited oral bioavailability and whole blood stability. Rapid metabolism with elimination via glutathione conjugation was observed. An oral dosing routine using medicated gels was developed to overcome bioavailability issues and yielded sustained whole blood levels above the half maximal effective concentration (EC50) in a 7-day study. The compound proved well-tolerated in acute and chronic experiments at 300 mg/kg PO dosing. The medicated gel formulation is highly suitable for evaluation of C16 in animal models of disease.

DNA copy number analysis has been instrumental for the identification of genetic alterations in B-cell precursor acute lymphoblastic leukemia. Notably, some of these genetic defects have been associated with poor treatment outcome and might be relevant for future risk stratification. In this study, we characterized recurrent deletions of CD200 and BTLA genes, mediated by recombination-activating genes, and used breakpoint-specific polymerase chain reaction assay to screen a cohort of 1154 cases of B-cell precursor acute lymphoblastic leukemia uniformly treated according to the EORTC-CLG 58951 protocol. CD200/BTLA deletions were identified in 56 of the patients (4.8%) and were associated with an inferior 8-year event free survival in this treatment protocol [70.2% ± 1.2% for patients with deletions versus 83.5% ± 6.4% for non-deleted cases (hazard ratio 2.02; 95% confidence interval 1.23–3.32; P=0.005)]. Genetically, CD200/BTLA deletions were strongly associated with ETV6-RUNX1-positive leukemias (P

It is still unclear why absence of the fragile X protein (FMRP) leads to mental retardation and specific behavioral problems. In neurons, the protein transports specific mRNAs towards the actively translating ribosomes near the synapses.To unravel the mechanism leading to the disorder, we performed global gene expression analysis by means of the differential display method using the fragile X mouse model. To verify differential expression, we used microarray technology and real-time PCR. Three differentially expressed cDNAs showed consistent underexpression in the fragile X knockout mouse, including a GABAA receptor subunit δ, a Rho guanine exchange factor 12 and an EST BU563433. In addition, we identified 5 genes that showed differential expression dependent on the sample of RNA analysis. We consider their differential expression as provisional. It is possible that these differentially expressed genes play an important role in the cognitive and behavioral problems observed in the fragile X syndrome.

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.

T-cell acute lymphoblastic leukemia arises from the leukemic transformation of developing thymocytes and results from cooperative genetic lesions. Inactivation of the PHF6 gene is frequently observed in T-cell acute lymphoblastic leukemia, suggesting an important tumor suppressive role for PHF6 in the pathobiology of this leukemia. Although the precise function of PHF6 is still unknown, this gene is most likely involved in chromatin regulation, a strongly emerging theme in T-cell acute lymphoblastic leukemia. In this context, our previous description of a cooperative microRNA regulatory network controlling several well-known T-cell acute lymphoblastic leukemia tumor suppressor genes, including PHF6, is of great importance. Given the high frequency of PHF6 lesions in T-cell acute lymphoblastic leukemia and the integration of PHF6 in this microRNA regulatory network, we aimed to identify novel oncogenic microRNAs in T-cell acute lymphoblastic leukemia which suppress PHF6. To this end, we performed an unbiased PHF6 3′UTR-microRNA library screen and combined the results with microRNA profiling data of samples from patients with T-cell acute lymphoblastic leukemia and normal thymocyte subsets. We selected miR-128-3p as a candidate PHF6-targeting, oncogenic microRNA and demonstrated regulation of PHF6 expression upon modulation of this microRNA in T-cell acute lymphoblastic leukemia cell lines. In vivo evidence of an oncogenic role of this microRNA in T-cell acute lymphoblastic leukemia was obtained through accelerated leukemia onset in a NOTCH1-induced T-cell acute lymphoblastic leukemia mouse model upon miR-128-3p over-expression. We conclude that miR-128-3p is a strong novel candidate oncogenic microRNA in T-cell acute lymphoblastic leukemia which targets the PHF6 tumor suppressor gene.

Structural genomic variations play an important role in human disease and phenotypic diversity. With the rise of high-throughput sequencing tools, mate-pair/paired-end/single-read sequencing has become an important technique for the detection and exploration of structural variation. Several analysis tools exist to handle different parts and aspects of such sequencing based structural variation analyses pipelines. A comprehensive analysis platform to handle all steps, from processing the sequencing data, to the discovery and visualization of structural variants, is missing. The ViVar platform is built to handle the discovery of structural variants, from Depth Of Coverage analysis, aberrant read pair clustering to split read analysis. ViVar provides you with powerful visualization options, enables easy reporting of results and better usability and data management. The platform facilitates the processing, analysis and visualization, of structural variation based on massive parallel sequencing data, enabling the rapid identification of disease loci or genes. ViVar allows you to scale your analysis with your work load over multiple (cloud) servers, has user access control to keep your data safe and is easy expandable as analysis techniques advance. URL: https://www.cmgg.be/vivar/

The selection and validation of stably expressed reference genes is a critical issue for proper RT-qPCR data normalization. In zebrafish expression studies, many commonly used reference genes are not generally applicable given their variability in expression levels under a variety of experimental conditions. Inappropriate use of these reference genes may lead to false interpretation of expression data and unreliable conclusions. In this study, we evaluated a novel normalization method in zebrafish using expressed repetitive elements (ERE) as reference targets, instead of specific protein coding mRNA targets. We assessed and compared the expression stability of a number of EREs to that of commonly used zebrafish reference genes in a diverse set of experimental conditions including a developmental time series, a set of different organs from adult fish and different treatments of zebrafish embryos including morpholino injections and administration of chemicals. Using geNorm and rank aggregation analysis we demonstrated that EREs have a higher overall expression stability compared to the commonly used reference genes. Moreover, we propose a limited set of ERE reference targets (hatn10, dna15ta1 and loopern4), that show stable expression throughout the wide range of experiments in this study, as strong candidates for inclusion as reference targets for qPCR normalization in future zebrafish expression studies. Our applied strategy to find and evaluate candidate expressed repeat elements for RT-qPCR data normalization has high potential to be used also for other species.

BACKGROUND:Neuroblastoma is a pediatric cancer that exhibits a wide clinical spectrum ranging from spontaneous regression in low-risk patients to fatal disease in high-risk patients. The identification of single nucleotide polymorphisms (SNPs) may help explain the heterogeneity of neuroblastoma and assist in identifying patients at higher risk for poor survival. SNPs in the TP53 pathway are of special importance, as several studies have reported associations between TP53 pathway SNPs and cancer. Of note, less than 2% of neuroblastoma tumors have a TP53 mutation at diagnosis. PATIENTS AND METHODS:We selected 21 of the most frequently studied SNPs in the TP53 pathway and evaluated their association with outcome in 500 neuroblastoma patients using TaqMan allelic discrimination assays. RESULTS AND CONCLUSION:We investigated the impact of 21 SNPs on overall survival, event-free survival, age at diagnosis, MYCN status, and stage of the disease in 500 neuroblastoma patients. A missense SNP in exon 10 of the CASP8 gene SNP D302H was associated with worse overall and event-free survival in patients with MYCN-amplified neuroblastoma tumors.

BACKGROUND: Measuring messenger RNA (mRNA) levels using the reverse transcription quantitative polymerase chain reaction (RT-qPCR) is common practice in many laboratories. A specific set of mRNAs as internal control reference genes is considered as the preferred strategy to normalize RT-qPCR data. Proper selection of reference genes is a critical issue, especially in cancer cells that are subjected to different in vitro manipulations. These manipulations may result in dramatic alterations in gene expression levels, even of assumed reference genes. In this study, we evaluated the expression levels of 11 commonly used reference genes as internal controls for normalization of 19 experiments that include neuroblastoma, T-ALL, melanoma, breast cancer, non small cell lung cancer (NSCL), acute myeloid leukemia (AML), prostate cancer, colorectal cancer, and cervical cancer cell lines subjected to various perturbations. RESULTS: The geNorm algorithm in the software package qbase+ was used to rank the candidate reference genes according to their expression stability. We observed that the stability of most of the candidate reference genes varies greatly in perturbation experiments. Expressed Alu repeats show relatively stable expression regardless of experimental condition. These Alu repeats are ranked among the best reference assays in all perturbation experiments and display acceptable average expression stability values (M

Chronic lymphocytic leukemia (CLL) is a disease with variable clinical outcome. Several prognostic factors such as the immunoglobulin heavy chain variable genes (IGHV) mutation status are linked to the B-cell receptor (BCR) complex, supporting a role for triggering the BCR in vivo in the pathogenesis. The miRNA profile upon stimulation and correlation with IGHV mutation status is however unknown. To evaluate the transcriptional response of peripheral blood CLL cells upon BCR stimulation in vitro, miRNA and mRNA expression was measured using hybridization arrays and qPCR. We found both IGHV mutated and unmutated CLL cells to respond with increased expression of MYC and other genes associated with BCR activation, and a phenotype of cell cycle progression. Genome-wide expression studies showed hsa-miR-132-3p/hsa-miR-212 miRNA cluster induction associated with a set of downregulated genes, enriched for genes modulated by BCR activation and amplified by Myc. We conclude that BCR triggering of CLL cells induces a transcriptional response of genes associated with BCR activation, enhanced cell cycle entry and progression and suggest that part of the transcriptional profiles linked to IGHV mutation status observed in isolated peripheral blood are not cell intrinsic but rather secondary to in vivo BCR stimulation.

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.

Chemotherapy induces apoptosis and tumor regression primarily through activation of p53-mediated transcription. Neuroblastoma is a p53 wild type malignancy at diagnosis and repression of p53 signaling plays an important role in its pathogenesis. Recently developed small molecule inhibitors of the MDM2-p53 interaction are able to overcome this repression and potently activate p53 dependent apoptosis in malignancies with intact p53 downstream signaling. We used the small molecule MDM2 inhibitor, Nutlin-3a, to determine the p53 drug response signature in neuroblastoma cells. In addition to p53 mediated apoptotic signatures, GSEA and pathway analysis identified a set of p53-repressed genes that were reciprocally over-expressed in neuroblastoma patients with the worst overall outcome in multiple clinical cohorts. Multifactorial regression analysis identified a subset of four genes (CHAF1A, RRM2, MCM3, and MCM6) whose expression together strongly predicted overall and event-free survival (p

Neuroblastoma is an embryonic tumor arising from immature sympathetic nervous system cells. Recurrent genomic alterations include MYCN and ALK amplification as well as recurrent patterns of gains and losses of whole or large partial chromosome segments. A recent whole genome sequencing effort yielded no frequently recurring mutations in genes other than those affecting ALK. However, the study further stresses the importance of DNA copy number alterations in this disease, in particular for genes implicated in neuritogenesis. Here we provide additional evidence for the importance of focal DNA copy number gains and losses, which are predominantly observed in MYCN amplified tumors. A focal 5 kb gain encompassing the MYCN regulated miR-17~92 cluster as sole gene was detected in a neuroblastoma cell line and further analyses of the array CGH data set demonstrated enrichment for other MYCN target genes in focal gains and amplifications. Next we applied an integrated genomics analysis to prioritize MYCN down regulated genes mediated by MYCN driven miRNAs within regions of focal heterozygous or homozygous deletion. We identified RGS5, a negative regulator of G-protein signaling implicated in vascular normalization, invasion and metastasis, targeted by a focal homozygous deletion, as a new MYCN target gene, down regulated through MYCN activated miRNAs. In addition, we expand the miR-17~92 regulatory network controlling TGFß signaling in neuroblastoma with the ring finger protein 11 encoding gene RNF11, which was previously shown to be targeted by the miR-17~92 member miR-19b. Taken together, our data indicate that focal DNA copy number imbalances in neuroblastoma (1) target genes that are implicated in MYCN signaling, possibly selected to reinforce MYCN oncogene addiction and (2) serve as a resource for identifying new molecular targets for treatment.

Several microRNAs (miRNAs) that are either specifically enriched or highly expressed in neurons and glia have been described, but the identification of miRNAs modulating neural stem cell (NSC) biology remains elusive. In this study, we exploited high throughput miRNA expression profiling to identify candidate miRNAs enriched in NSC/early progenitors derived from the murine subventricular zone (SVZ). Then, we used lentiviral miRNA sensor vectors (LV.miRT) to monitor the activity of shortlisted miRNAs with cellular and temporal resolution during NSC differentiation, taking advantage of in vitro and in vivo models that recapitulate physiological neurogenesis and gliogenesis and using known neuronal- and glial-specific miRNAs as reference. The LV.miRT platform allowed us monitoring endogenous miRNA activity in low represented cell populations within a bulk culture or within the complexity of CNS tissue, with high sensitivity and specificity. In this way we validated and extended previous results on the neuronal-specific miR-124 and the astroglial-specific miR-23a. Importantly, we describe for the first time a cell type- and differentiation stage-specific modulation of miR-93 and miR-125b in SVZ-derived NSC cultures and in the SVZ neurogenic niche in vivo, suggesting key roles of these miRNAs in regulating NSC function.

One of the first and most important steps in the metastatic cascade is the loss of cell-cell and cell-matrix interactions. N-cadherin, a crucial mediator of homotypic and heterotypic cell-cell interactions, might play a central role in the metastasis of neuroblastoma (NB), a solid tumor of neuroectodermal origin. Using Reverse Transcription Quantitative PCR (RT-qPCR), Western blot, immunocytochemistry and Tissue MicroArrays (TMA) we demonstrate the expression of N-cadherin in neuroblastoma tumors and cell lines. All neuroblastic tumors (n = 356) and cell lines (n = 10) expressed various levels of the adhesion protein. The N-cadherin mRNA expression was significantly lower in tumor samples from patients suffering metastatic disease. Treatment of NB cell lines with the N-cadherin blocking peptide ADH-1 (Exherin, Adherex Technologies Inc.), strongly inhibited tumor cell proliferation in vitro by inducing apoptosis. Our results suggest that N-cadherin signaling may play a role in neuroblastoma disease, marking involvement of metastasis and determining neuroblastoma cell viability.

BACKGROUND: Although the throughput of next generation sequencing is increasing and at the same time the cost is substantially reduced, for the majority of laboratories whole genome sequencing of large cohorts of cancer samples is still not feasible. In addition, the low number of genomes that are being sequenced is often problematic for the downstream interpretation of the significance of the variants. Targeted resequencing can partially circumvent this problem; by focusing on a limited number of candidate cancer genes to sequence, more samples can be included in the screening, hence resulting in substantial improvement of the statistical power. In this study, a successful strategy for prioritizing candidate genes for targeted resequencing of cancer genomes is presented. RESULTS: Four prioritization strategies were evaluated on six different cancer types: genes were ranked using these strategies, and the positive predictive value (PPV) or mutation rate within the top-ranked genes was compared to the baseline mutation rate in each tumor type. Successful strategies generate gene lists in which the top is enriched for known mutated genes, as evidenced by an increase in PPV. A clear example of such an improvement is seen in colon cancer, where the PPV is increased by 2.3 fold compared to the baseline level when 100 top fitSNP genes are sequenced. CONCLUSIONS: A gene prioritization strategy based on the fitSNP scores appears to be most successful in identifying mutated cancer genes across different tumor entities, with variance of gene expression levels as a good second best.

Fluorescence in situ hybridization was performed to characterize 81 cases of myeloid and lymphoid malignancies with cytogenetic 1p36 alterations not affecting the PRDM16 locus. In total, three subgroups were identified: balanced translocations (N = 27) and telomeric rearrangements (N = 15), both mainly observed in myeloid disorders; and unbalanced non-telomeric rearrangements (N = 39), mainly observed in lymphoid proliferations and frequently associated with a highly complex karyotype. The 1p36 rearrangement was isolated in 12 cases, mainly myeloid disorders. The breakpoints on 1p36 were more widely distributed than previously reported, but with identifiable rare breakpoint cluster regions, such as the TP73 locus. We also found novel partner loci on 1p36 for the known multi-partner genes HMGA2 and RUNX1. We precised the common terminal 1p36 deletion, which has been suggested to have an adverse prognosis, in B-cell lymphomas [follicular lymphomas and diffuse large B-cell lymphomas with t(14;18)(q32;q21) as well as follicular lymphomas without t(14;18)]. Intrachromosomal telomeric repetitive sequences were detected in at least half the cases of telomeric rearrangements. It is unclear how the latter rearrangements occurred and whether they represent oncogenic events or result from chromosomal instability during oncogenesis.

Chromosomal rearrangements involving the EVI1 proto-oncogene are a recurrent finding in myeloid leukemias and are indicative of a poor prognosis. Rearrangements of the EVI1 locus are often associated with monosomy 7 or cytogenetic detectable deletions of part of 7q. As EVI1 overexpression alone is not sufficient to induce leukemia, loss of a 7q tumour suppressor gene might be a required cooperating event. To test this hypothesis, we performed high-resolution array comparative genomic hybridization analysis of twelve EVI1 overexpressing patients and three EVI1 deregulated cell lines to search for 7q submicroscopic deletions. This analysis lead to the delineation of two critical regions, one of 0.39 Mb on 7q35 containing the CNTNAP2 gene and one of 1.33 Mb on chromosome bands 7q35-q36 comprising nine genes in EVI1 deregulated cell lines. These findings open the way to further studies aimed at identifying the culprit EVI1 implicated tumour suppressor genes on 7q.

MiRNAs regulate gene expression at a post-transcriptional level and their dysregulation can play major roles in the pathogenesis of many different forms of cancer, including neuroblastoma, an often fatal paediatric cancer originating from precursor cells of the sympathetic nervous system. We have analyzed a set of neuroblastoma (n = 145) that is broadly representative of the genetic subtypes of this disease for miRNA expression (430 loci by stem-loop RT qPCR) and for DNA copy number alterations (array CGH) to assess miRNA involvement in disease pathogenesis. The tumors were stratified and then randomly split into a training set (n = 96) and a validation set (n = 49) for data analysis. Thirty-seven miRNAs were significantly over- or under-expressed in MYCN amplified tumors relative to MYCN single copy tumors, indicating a potential role for the MYCN transcription factor in either the direct or indirect dysregulation of these loci. In addition, we also determined that there was a highly significant correlation between miRNA expression levels and DNA copy number, indicating a role for large-scale genomic imbalances in the dysregulation of miRNA expression. In order to directly assess whether miRNA expression was predictive of clinical outcome, we used the Random Forest classifier to identify miRNAs that were most significantly associated with poor overall patient survival and developed a 15 miRNA signature that was predictive of overall survival with 72.7% sensitivity and 86.5% specificity in the validation set of tumors. We conclude that there is widespread dysregulation of miRNA expression in neuroblastoma tumors caused by both over-expression of the MYCN transcription factor and by large-scale chromosomal imbalances. MiRNA expression patterns are also predicative of clinical outcome, highlighting the potential for miRNA mediated diagnostics and therapeutics.

To date, the contribution of disrupted potentially cis-regulatory conserved non-coding sequences (CNCs) to human disease is most likely underestimated, as no systematic screens for putative deleterious variations in CNCs have been conducted. As a model for monogenic disease we studied the involvement of genetic changes of CNCs in the cis-regulatory domain of FOXL2 in blepharophimosis syndrome (BPES). Fifty-seven molecularly unsolved BPES patients underwent high-resolution copy number screening and targeted sequencing of CNCs. Apart from three larger distant deletions, a de novo deletion as small as 7.4 kb was found at 283 kb 5' to FOXL2. The deletion appeared to be triggered by an H-DNA-induced double-stranded break (DSB). In addition, it disrupts a novel long non-coding RNA (ncRNA) PISRT1 and 8 CNCs. The regulatory potential of the deleted CNCs was substantiated by in vitro luciferase assays. Interestingly, Chromosome Conformation Capture (3C) of a 625 kb region surrounding FOXL2 in expressing cellular systems revealed physical interactions of three upstream fragments and the FOXL2 core promoter. Importantly, one of these contains the 7.4 kb deleted fragment. Overall, this study revealed the smallest distant deletion causing monogenic disease and impacts upon the concept of mutation screening in human disease and developmental disorders in particular.

Chromosomal translocations involving the EVI1 locus are a recurrent finding in myeloid leukemia and are associated with poor prognosis. In this study, we performed a detailed molecular characterization of the recurrent translocation t(3;17)(q26;q22) in 13 hematologic malignancies. The EVI1 gene locus was rearranged in all 13 patients and was associated with EVI1 overexpression. In 9 out of 13 patients, the 17q breakpoints clustered in a 250 kb region on band 17q22 encompassing the MSI2 (musashi homologue 2) gene. Expression analyses failed to demonstrate ectopic MSI2 expression or the presence of an MSI2/EVI1 fusion gene. In conclusion, we show for the first time that the t(3;17) is indeed a recurrent chromosomal aberration in myeloid malignancies. In keeping with findings in other recurrent 3q26 rearrangements, overexpression of the EVI1 gene appears to be the major contributor to leukemogenesis in patients with a t(3;17).

The human 1p36 region is deleted in many different types of tumors, and so it probably harbors one or more tumor suppressor genes. In a Belgian neuroblastoma patient, a constitutional balanced translocation t(1;17)(p36.2;q11.2) may have led to the development of the tumor by disrupting or activating a gene. Here, we report the cloning of both translocation breakpoints and the identification of a novel gene that is disrupted by this translocation. This gene, named NBPF1 for Neuroblastoma BreakPoint Family member 1, belongs to a recently described gene family encoding highly similar proteins, the functions of which are unknown. The translocation truncates NBPF1 and gives rise to two chimeric transcripts of NBPF1 sequences fused to sequences derived from chromosome 17. On chromosome 17, the translocation disrupts one of the isoforms of ACCN1, a potential glioma tumor suppressor gene. Expression of the NBPF family in neuroblastoma cell lines is highly variable, but it is decreased in cell lines that have a deletion of chromosome 1p. More importantly, expression profiling of the NBPF1 gene showed that its expression is significantly lower in cell lines with heterozygous NBPF1 loss than in cell lines with a normal 1p chromosome. Meta-analysis of the expression of NBPF and ACCN1 in neuroblastoma tumors indicates a role for the NBPF genes and for ACCN1 in tumor aggressiveness. Additionally, DLD1 cells with inducible NBPF1 expression showed a marked decrease of clonal growth in a soft agar assay. The disruption of both NBPF1 and ACCN1 genes in this neuroblastoma patient indicates that these genes might suppress development of neuroblastoma and possibly other tumor types.

Background and Objectives Classical Hodgkin’s lymphoma (cHL) is a common malignant lymphoma characterized by the presence of large, usually multinucleated malignant Hodgkin and Reed Sternberg (HRS) cells which are thought to be derived from germinal center B-cells. In cHL, the HRS cells constitute less than 1% of the tumor volume; consequently the profile of genetic aberrations in cHL is still poorly understood.Design and Methods In this study, we subjected four commonly used cHL cell lines to array comparative genomic hybridization (aCGH) in order to delineate known chromosomal aberrations in more detail and to search for small hitherto undetected genomic imbalances.Results The aCGH profiles of the four cell lines tested confirmed the complex patterns of rearrangements previously demonstrated with multicolor fluorescence in situ hybridization and chromosomal CGH (cCGH). Importantly, aCGH allowed a much more accurate delineation of imbalances as compared to previous studies performed at a chromosomal level of resolution. Furthermore, we detected 35 hitherto undetected aberrations including a homozygous deletion of chromosomal region 15q26.2 in the cell line HDLM2 encompasing RGMA and CHD2 and an amplification of the STAT6 gene in cell line L1236 leading to STAT6 overexpression. Finally, in cell line KM-H2 we found a 2.35 Mb deletion at 16q12.1 putatively defining a small critical region for the recurrent 16q deletion in cHL. This region contains the CYLD gene, a known suppressor gene of the NF-κB pathway.Interpretation and Conclusions aCGH was performed on four cHL cell lines leading to the improved delineation of known chromosomal imbalances and the detection of 35 hitherto undetected aberrations. More specifically, our results highlight STAT6 as a potential transcriptional target and identified RGMA, CHD2 and CYLD as candidate tumor suppressors in cHL.