Showing papers by "Jo Vandesompele published in 2016"

Melanoma addiction to the long non-coding RNA SAMMSON.

Abstract: Focal amplifications of chromosome 3p13-3p14 occur in about 10% of melanomas and are associated with a poor prognosis. The melanoma-specific oncogene MITF resides at the epicentre of this amplicon. However, whether other loci present in this amplicon also contribute to melanomagenesis is unknown. Here we show that the recently annotated long non-coding RNA (lncRNA) gene SAMMSON is consistently co-gained with MITF. In addition, SAMMSON is a target of the lineage-specific transcription factor SOX10 and its expression is detectable in more than 90% of human melanomas. Whereas exogenous SAMMSON increases the clonogenic potential in trans, SAMMSON knockdown drastically decreases the viability of melanoma cells irrespective of their transcriptional cell state and BRAF, NRAS or TP53 mutational status. Moreover, SAMMSON targeting sensitizes melanoma to MAPK-targeting therapeutics both in vitro and in patient-derived xenograft models. Mechanistically, SAMMSON interacts with p32, a master regulator of mitochondrial homeostasis and metabolism, to increase its mitochondrial targeting and pro-oncogenic function. Our results indicate that silencing of the lineage addiction oncogene SAMMSON disrupts vital mitochondrial functions in a cancer-cell-specific manner; this silencing is therefore expected to deliver highly effective and tissue-restricted anti-melanoma therapeutic responses.

Asthma inflammatory phenotypes show differential microRNA expression in sputum

Abstract: Background Asthma is classified according to severity and inflammatory phenotype and is likely to be distinguished by specific microRNA (miRNA) expression profiles. Objective We sought to associate miRNA expression in sputum supernatants with the inflammatory cell profile and disease severity in asthmatic patients. Methods We investigated miRNA expression in sputum supernatants of 10 healthy subjects, 17 patients with mild-to-moderate asthma, and 9 patients with severe asthma using high-throughput, stem-loop, reverse transcriptase quantitative real-time PCR miRNA expression profiling (screening cohort, n = 36). Differentially expressed miRNAs were validated in an independent cohort (n = 60; 10 healthy subjects and 50 asthmatic patients). Cellular miRNA origin was examined by using in situ hybridization and reverse transcriptase quantitative real-time PCR. The functional role of miRNAs was assessed by using in silico analysis and in vitro transfecting miRNA mimics in human bronchial epithelial cells. Results In 2 independent cohorts expression of miR-629-3p, miR-223-3p, and miR-142-3p was significantly upregulated in sputum of patients with severe asthma compared with that in healthy control subjects and was highest in patients with neutrophilic asthma. Expression of the 3 miRNAs was associated with sputum neutrophilia, and miR-223-3p and miR-142-3p expression was associated also with airway obstruction (FEV 1 /forced vital capacity). Expression of miR-629-3p was localized in the bronchial epithelium, whereas miR-223-3p and miR-142-3p were expressed in neutrophils, monocytes, and macrophages. Transfecting human bronchial epithelial cells with miR-629-3p mimic induced epithelial IL-8 mRNA and protein expression. IL-1β and IL-8 protein levels were significantly increased in sputum of patients with severe asthma and were positively associated with sputum neutrophilia. Conclusions Expression of miR-223-3p, miR-142-3p, and miR-629-3p is increased in sputum of patients with severe asthma and is linked to neutrophilic airway inflammation, suggesting that these miRNAs contribute to this asthma inflammatory phenotype.

Benchmarking of RNA-sequencing analysis workflows using whole-transcriptome RT-qPCR expression data

Abstract: RNA-sequencing has become the gold standard for whole-transcriptome gene expression quantification. Multiple algorithms have been developed to derive gene counts from sequencing reads. While a number of benchmarking studies have been conducted, the question remains how individual methods perform at accurately quantifying gene expression levels from RNA-sequencing reads. We performed an independent benchmarking study using RNA-sequencing data from the well established MAQCA and MAQCB reference samples. RNA-sequencing reads were processed using five workflows (Tophat-HTSeq, Tophat-Cufflinks, STAR-HTSeq, Kallisto and Salmon) and resulting gene expression measurements were compared to expression data generated by wet-lab validated qPCR assays for all protein coding genes. All methods showed high gene expression correlations with qPCR data. When comparing gene expression fold changes between MAQCA and MAQCB samples, about 85% of the genes showed consistent results between RNA-sequencing and qPCR data. Of note, each method revealed a small but specific gene set with inconsistent expression measurements. A significant proportion of these method-specific inconsistent genes were reproducibly identified in independent datasets. These genes were typically smaller, had fewer exons, and were lower expressed compared to genes with consistent expression measurements. We propose that careful validation is warranted when evaluating RNA-seq based expression profiles for this specific gene set.

Renal microRNA- and RNA-profiles in progressive chronic kidney disease

Abstract: Background MicroRNAs (miRNAs) contribute to chronic kidney disease (CKD) progression via regulating mRNAs involved in renal homeostasis. However, their association with clinical outcome remains poorly understood. Materials and methods We performed miRNA and mRNA expression profiling on renal biopsy sections by qPCR (miRNA) and microarrays (mRNA) in a discovery (n = 43) and in a validation (n = 29) cohort. miRNAs differentiating stable and progressive cases were inversely correlated with putative target mRNAs, which were further characterized by pathway analysis using KEGG pathways. Results miR-30d, miR-140-3p, miR-532-3p, miR-194, miR-190, miR-204 and miR-206 were downregulated in progressive cases. These seven miRNAs correlated with upregulated 29 target mRNAs involved in inflammatory response, cell–cell interaction, apoptosis and intra-cellular signalling. In particular, miR-206 and miR-532-3p were associated with distinct biological processes via the expression of their target mRNAs: Reduced expression of miR-206 in progressive disease correlated with the upregulation of target mRNAs participating in inflammatory pathways (CCL19, CXCL1, IFNAR2, NCK2, PTK2B, PTPRC, RASGRP1 and TNFRSF25). Progressive cases also showed a lower expression of miR-532-3p and an increased expression of target transcripts involved in apoptosis pathways (MAP3K14, TNFRSF10B/TRAIL-R2, TRADD and TRAF2). In the validation cohort, we confirmed the decreased expression of miR-206 and miR-532-3p, and the inverse correlation of these miRNAs with the expression of nine of the 12 target genes. The levels of the identified miRNAs and the target mRNAs correlated with clinical parameters and histological damage indices. Conclusions These results suggest the involvement of specific miRNAs and mRNAs in biological pathways associated with the progression of CKD.

Straightforward and sensitive RT-qPCR based gene expression analysis of FFPE samples

Abstract: Fragmented RNA from formalin-fixed paraffin-embedded (FFPE) tissue is a known obstacle to gene expression analysis. In this study, the impact of RNA integrity, gene-specific reverse transcription and targeted cDNA preamplification was quantified in terms of reverse transcription polymerase chain reaction (RT-qPCR) sensitivity by measuring 48 protein coding genes on eight duplicate cultured cancer cell pellet FFPE samples and twenty cancer tissue FFPE samples. More intact RNA modestly increased gene detection sensitivity by 1.6 fold (earlier detection by 0.7 PCR cycles, 95% CI = 0.593–0.850). Application of gene-specific priming instead of whole transcriptome priming during reverse transcription further improved RT-qPCR sensitivity by a considerable 4.0 fold increase (earlier detection by 2.0 PCR cycles, 95% CI = 1.73–2.32). Targeted cDNA preamplification resulted in the strongest increase of RT-qPCR sensitivity and enabled earlier detection by an average of 172.4 fold (7.43 PCR cycles, 95% CI = 6.83–7.05). We conclude that gene-specific reverse transcription and targeted cDNA preamplification are adequate methods for accurate and sensitive RT-qPCR based gene expression analysis of FFPE material. The presented methods do not involve expensive or complex procedures and can be easily implemented in any routine RT-qPCR practice.

Long noncoding RNA signatures define oncogenic subtypes in T-cell acute lymphoblastic leukemia.

Abstract: Long noncoding RNA signatures define oncogenic subtypes in T-cell acute lymphoblastic leukemia

miSTAR: miRNA target prediction through modeling quantitative and qualitative miRNA binding site information in a stacked model structure.

Abstract: In microRNA (miRNA) target prediction, typically two levels of information need to be modeled: the number of potential miRNA binding sites present in a target mRNA and the genomic context of each individual site. Single model structures insufficiently cope with this complex training data structure, consisting of feature vectors of unequal length as a consequence of the varying number of miRNA binding sites in different mRNAs. To circumvent this problem, we developed a two-layered, stacked model, in which the influence of binding site context is separately modeled. Using logistic regression and random forests, we applied the stacked model approach to a unique data set of 7990 probed miRNA-mRNA interactions, hereby including the largest number of miRNAs in model training to date. Compared to lower-complexity models, a particular stacked model, named miSTAR (miRNA stacked model target prediction; www.mi-star.org), displays a higher general performance and precision on top scoring predictions. More importantly, our model outperforms published and widely used miRNA target prediction algorithms. Finally, we highlight flaws in cross-validation schemes for evaluation of miRNA target prediction models and adopt a more fair and stringent approach.

Differential expression of lncRNAs during the HIV replication cycle: an underestimated layer in the HIV-host interplay.

Abstract: Studying the effects of HIV infection on the host transcriptome has typically focused on protein-coding genes However, recent advances in the field of RNA sequencing revealed that long non-coding RNAs (lncRNAs) add an extensive additional layer to the cell’s molecular network Here, we performed transcriptome profiling throughout a primary HIV infection in vitro to investigate lncRNA expression at the different HIV replication cycle processes (reverse transcription, integration and particle production) Subsequently, guilt-by-association, transcription factor and co-expression analysis were performed to infer biological roles for the lncRNAs identified in the HIV-host interplay Many lncRNAs were suggested to play a role in mechanisms relying on proteasomal and ubiquitination pathways, apoptosis, DNA damage responses and cell cycle regulation Through transcription factor binding analysis, we found that lncRNAs display a distinct transcriptional regulation profile as compared to protein coding mRNAs, suggesting that mRNAs and lncRNAs are independently modulated In addition, we identified five differentially expressed lncRNA-mRNA pairs with mRNA involvement in HIV pathogenesis with possible cis regulatory lncRNAs that control nearby mRNA expression and function Altogether, the present study demonstrates that lncRNAs add a new dimension to the HIV-host interplay and should be further investigated as they may represent targets for controlling HIV replication

Methyl-CpG-binding domain sequencing reveals a prognostic methylation signature in neuroblastoma.

Abstract: Accurate assessment of neuroblastoma outcome prediction remains challenging. Therefore, this study aims at establishing novel prognostic tumor DNA methylation biomarkers. In total, 396 low- and high-risk primary tumors were analyzed, of which 87 were profiled using methyl-CpG-binding domain (MBD) sequencing for differential methylation analysis between prognostic patient groups. Subsequently, methylation-specific PCR (MSP) assays were developed for 78 top-ranking differentially methylated regions and tested on two independent cohorts of 132 and 177 samples, respectively. Further, a new statistical framework was used to identify a robust set of MSP assays of which the methylation score (i.e. the percentage of methylated assays) allows accurate outcome prediction. Survival analyses were performed on the individual target level, as well as on the combined multimarker signature. As a result of the differential DNA methylation assessment by MBD sequencing, 58 of the 78 MSP assays were designed in regions previously unexplored in neuroblastoma, and 36 are located in non-promoter or non-coding regions. In total, 5 individual MSP assays (located in CCDC177, NXPH1, lnc-MRPL3-2, lnc-TREX1-1 and one on a region from chromosome 8 with no further annotation) predict event-free survival and 4 additional assays (located in SPRED3, TNFAIP2, NPM2 and CYYR1) also predict overall survival. Furthermore, a robust 58-marker methylation signature predicting overall and event-free survival was established. In conclusion, this study encompasses the largest DNA methylation biomarker study in neuroblastoma so far. We identified and independently validated several novel prognostic biomarkers, as well as a prognostic 58-marker methylation signature.

Stage 4S neuroblastoma tumors show a characteristic DNA methylation portrait

Abstract: Stage 4S neuroblastoma (NB) is a special type of NB found in infants with metastases at diagnosis and is associated with an excellent outcome due to its remarkable capacity to undergo spontaneous regression. As genomics have not been able to explain this intriguing clinical presentation, we here aimed at profiling the DNA methylome of stage 4S NB to better understand this phenomenon. To this purpose, differential methylation analyses between International Neuroblastoma Staging System (INSS) stage 4S, stage 4 and stage 1/2 were performed, using methyl-CpG-binding domain (MBD) sequencing data of 14 stage 4S, 14 stage 4, and 13 stage 1/2 primary NB tumors (all MYCN non-amplified in order not to confound results). Stage 4S-specific hyper- and hypomethylated promoters were determined and further characterized for genomic localization and function by cytogenetic band enrichment, gene set enrichment, transcription factor target enrichment and differential RNA expression analyses. We show that specific chromosomal locations are enriched for stage 4S differentially methylated promoters and that stage 4S tumors show characteristic hypermethylation of specific subtelomeric promoters. Furthermore, genes involved in important oncogenic pathways, in neural crest development and differentiation, and in epigenetic processes are differentially methylated and expressed in stage 4S tumors. Based on these findings, we describe new biological mechanisms possibly contributing to the stage 4S-specific tumor biology and spontaneous regression. In conclusion, this study is the first to describe the highly characteristic stage 4S DNA methylome. These findings will open new avenues to further unravel the NB pathology in general and stage 4S disease specifically.

Depletion of tRNA-halves enables effective small RNA sequencing of low-input murine serum samples

Abstract: The ongoing ascent of sequencing technologies has enabled researchers to gain unprecedented insights into the RNA content of biological samples. MiRNAs, a class of small non-coding RNAs, play a pivotal role in regulating gene expression. The discovery that miRNAs are stably present in circulation has spiked interest in their potential use as minimally-invasive biomarkers. However, sequencing of blood-derived samples (serum, plasma) is challenging due to the often low RNA concentration, poor RNA quality and the presence of highly abundant RNAs that dominate sequencing libraries. In murine serum for example, the high abundance of tRNA-derived small RNAs called 5′ tRNA halves hampers the detection of other small RNAs, like miRNAs. We therefore evaluated two complementary approaches for targeted depletion of 5′ tRNA halves in murine serum samples. Using a protocol based on biotinylated DNA probes and streptavidin coated magnetic beads we were able to selectively deplete 95% of the targeted 5′ tRNA half molecules. This allowed an unbiased enrichment of the miRNA fraction resulting in a 6-fold increase of mapped miRNA reads and 60% more unique miRNAs detected. Moreover, when comparing miRNA levels in tumor-carrying versus tumor-free mice, we observed a three-fold increase in differentially expressed miRNAs.

DNA methylation profiling of primary neuroblastoma tumors using methyl-CpG-binding domain sequencing

Abstract: Comprehensive genome-wide DNA methylation studies in neuroblastoma (NB), a childhood tumor that originates from precursor cells of the sympathetic nervous system, are scarce. Recently, we profiled the DNA methylome of 102 well-annotated primary NB tumors by methyl-CpG-binding domain (MBD) sequencing, in order to identify prognostic biomarker candidates. In this data descriptor, we give details on how this data set was generated and which bioinformatics analyses were applied during data processing. Through a series of technical validations, we illustrate that the data are of high quality and that the sequenced fragments represent methylated genomic regions. Furthermore, genes previously described to be methylated in NB are confirmed. As such, these MBD sequencing data are a valuable resource to further study the association of NB risk factors with the NB methylome, and offer the opportunity to integrate methylome data with other -omic data sets on the same tumor samples such as gene copy number and gene expression, also publically available.

RT-qPCR gene expression analysis in zebrafish: Preanalytical precautions and use of expressed repetitive elements for normalization.

Abstract: Gene expression analysis is increasingly important in many fields of biological research. Understanding patterns of expressed genes is assumed to provide insight into complex regulatory networks and can lead to the identification of genes relevant to specific biological processes, including disease. Among different techniques, reverse transcription quantitative polymerase chain reaction (RT-qPCR) is currently regarded as the gold standard for targeted quantification of RNA gene expression, especially because of its high sensitivity, specificity, accuracy, and precision, and also because of its practical simplicity and processing speed. However, different critical factors can influence the outcome of RT-qPCR studies, including isolation of RNA, reverse transcription to cDNA, and data analysis. These factors need to be addressed in order to obtain biologically meaningful results. In this chapter, we describe how RT-qPCR can be used in a reliable way to successfully study differential gene expression in zebrafish. Hereby, we especially focus on how expressed repetitive elements can be employed as reference targets in zebrafish RT-qPCR studies and how they can further improve the quality of the data.

Long non-coding RNA expression profiling in the NCI60 cancer cell line panel using high-throughput RT-qPCR.

Abstract: Long non-coding RNAs (lncRNAs) form a new class of RNA molecules implicated in various aspects of protein coding gene expression regulation. To study lncRNAs in cancer, we generated expression profiles for 1707 human lncRNAs in the NCI60 cancer cell line panel using a high-throughput nanowell RT-qPCR platform. We describe how qPCR assays were designed and validated and provide processed and normalized expression data for further analysis. Data quality is demonstrated by matching the lncRNA expression profiles with phenotypic and genomic characteristics of the cancer cell lines. This data set can be integrated with publicly available omics and pharmacological data sets to uncover novel associations between lncRNA expression and mRNA expression, miRNA expression, DNA copy number, protein coding gene mutation status or drug response.

Quantitative and Functional Requirements for Bioluminescent Cancer Models.

Abstract: Background: Bioluminescent cancer models are widely used but detailed quantification of the luciferase signal and functional comparison with a non-transfected control cell line are generally lacking. In the present study, we provide quantitative and functional tests for luciferase-transfected cells. Materials and Methods: We quantified the luciferase expression in BLM and HCT8/E11 transfected cancer cells, and examined the effect of long-term luciferin exposure. The present study also investigated functional differences between parental and transfected cancer cells. Results: Our results showed that quantification of different single-cell-derived populations are superior with droplet digital polymerase chain reaction. Quantification of luciferase protein level and luciferase bioluminescent activity is only useful when there is a significant difference in copy number. Continuous exposure of cell cultures to luciferin leads to inhibitory effects on mitochondrial activity, cell growth and bioluminescence. These inhibitory effects correlate with luciferase copy number. Cell culture and mouse xenograft assays showed no significant functional differences between luciferase-transfected and parental cells. Conclusion: Luciferase-transfected cells should be validated by quantitative and functional assays before starting large-scale experiments.

Abstract A091: Specific myelomonocytic cells heavily infiltrate orthotopic lung tumors and display a hypoxia-driven miRNA expression signature that directs tumor-supporting functions and negatively impacts on clinical outcome

Abstract: Tumor-derived micro-RNA expression patterns are known to convey prognostic information in lung cancer. However, the cellular source of these signatures is unclear. In parallel, immunological parameters such as the presence and activation status of dendritic cells (DCs) emerge as important prognostic factors in this disease as well, presumably by affecting the quality of immune surveillance.We hypothesized that tumor-associated DCs are involved in the way tumor-derived miRNA expression patterns impact on disease outcome.Using a preclinical model of lung cancer featuring orthotopic tumor growth in syngeneic, immunocompetent hosts, we found that lung tumors were preferentially infiltrated by CD11c+/MHCII+/CD11b+/CD103- DCs (TIDCs). Compared to their peritumoral counterparts, displayed a dramatic overexpression of PD-L1 along with markers associated with tumor-associated macrophages (TAM). Transcriptome analysis of these CD11b+ TIDCs vs peritumoral CD11b+ DCs using GSEA and ISMARA indicated loss of immunogenicity, upregulation of features related to immunosuppressive TAMs, and exposure to hypoxia. Global miRNA profiling of CD11b+TIDCs revealed a specific signature dominated by a strong overexpression of mir-31, a prototypical lung cancer oncomir. Exposure to hypoxia appeared as a powerful driver of endogenous miR-31 expression in DCs. Confocal microscopy of intrapulmonary tumors confirmed that CD11c+ dendritic-shaped leukocytes colocalize within hypoxic regions. We found that miR-31 overexpression and hypoxia alike cause DCs to secrete the tumor-supporting factors VEGF, S100A8 and S100A9. Conditioned medium of miR-31-overexpression DCs induced pro-invasive cell morphology changes in lung carcinoma cells cultured on a collagen matrix, a phenomenon that could be reproduced by exposure to recombinant S100A8 protein. By re-processing TCGA expression data for non-coding RNA sequences, we discovered that the lung TIDC miRNA signature imparts a strong negative impact on overall survival in non-squamous lung cancer. Moreover, the TIDC-derived miRNA signature induced a profound difference in outcome among lung cancer patients with regional lymph node invasion.In summary, lung tumors predominantly recruit inflammatory-type DCs that are reprogrammed at the level of immuno-phenotype and micro-RNA repertoire. The miRNA signature of lung tumor-infiltrating DCs is dominated by the miR-31, a hypoxia-induced oncomir stimulating the release of pro-metastatic proteins by DCs. Finally, we show for the first time that a miRNA signature derived from a stromal immune cell has a dramatic impact on long term outcome in lung cancer. Citation Format: Elisabeth Brabants, Lotte Pyfferoen, Celine Everaert, Simon Tavernier, Kelly Heyns, Nancy De Cabooter, Glenn Wagemans, Kim Deswarte, Hamida Hammad, Olivier De Wever, Jo Vandesompele, Bart Lambrecht, Pieter Mestdagh, Karim Vermaelen. Specific myelomonocytic cells heavily infiltrate orthotopic lung tumors and display a hypoxia-driven miRNA expression signature that directs tumor-supporting functions and negatively impacts on clinical outcome [abstract]. In: Proceedings of the Second CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; 2016 Sept 25-28; New York, NY. Philadelphia (PA): AACR; Cancer Immunol Res 2016;4(11 Suppl):Abstract nr A091.

Abstract 2664: Eradication of neuroblastoma by suppressing the expression of a single long noncoding RNA

Abstract: N-Myc gene amplification occurs in one quarter of human neuroblastoma tissues, and is a marker for poor patient prognosis. We performed RNA sequencing experiments, and identified 5 transcripts, including RP1NB1, which were most considerably differentially expressed between N-Myc gene amplified and nonamplified human neuroblastoma cell lines. Affymetrix microarray studies revealed that DEPD was one of the few genes considerably downregulated in neuroblastoma cells after RP1NB1 depletion. Chromatin immunoprecipitation assays showed that knocking down RP1NB1 expression reduced histone H3 lysine 4 trimethylation, a marker for active gene transcription, at the DEPD gene promoter. Luciferase assays demonstrated that knocking down RP1NB1 decreased DEPD gene promoter activity. Depletion of RP1NB1 or DEPD with two independent siRNAs or shRNAs significantly reduced ERK protein phosphorylation, N-Myc protein phosphorylation at Serine 62, N-Myc protein stabilization, neuroblastoma cell proliferation and survival. Clonogenic assays showed that knocking down RP1NB1 with doxycycline completely abolished colony formation capacity of neuroblastoma cells stably transfected with doxycycline-inducible RP1NB1 shRNAs. Importantly, treatment with doxycycline in mice xenografted with neuroblastoma cells stably transfected with doxycycline-inducible RP1NB1 shRNA led to tumor eradication. In human neuroblastoma tissues from 600 neuroblastoma patients, high levels of RP1NB1 gene expression correlated with DEPD gene expression and poor patient prognosis. In conclusion, this study identifies the novel long noncoding RNA RP1NB1 as an important regulator of N-Myc protein stability and neuroblastoma tumorigenesis. Citation Format: Andrew E. Tee, Pei Y. Liu, Giorgio Milazzo, Kate M. Hannan, Jesper Maag, Nenad Bartonicek, Renhua Song, Chen C. Jiang, Xu D. Zhang, Murray D. Norris, Michelle Haber, Glenn M. Marshall, Jinyan Li, Jo Vandesompele, John S. Mattick, Pieter Mestdagh, Giovanni Perini, Ross D. Hannan, Marcel E. Dinger, Tao Liu. Eradication of neuroblastoma by suppressing the expression of a single noncoding RNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2453.

Methyl-CpG-binding domain sequencing reveals a prognostic methylation signature in neuroblastoma

Abstract: Neuroblastoma, a childhood tumor that originates from precursor cells of the sympathetic nervous system, is a heterogeneous disease with prognosis ranging from long-term survival for localized tumors to fatal outcome for metastatic disease. One of the main challenges in the management of neuroblastoma remains accurate outcome prediction, enabling the choice of risk-related therapy. Therefore, we aimed at establishing novel prognostic tumor DNA methylation biomarkers. In total, 396 low- and high-risk primary tumors were analyzed, of which 87 were profiled using methyl-CpG-binding domain (MBD) sequencing for differential methylation analysis between prognostic patient groups. Subsequently, methylation-specific PCR (MSP) assays were developed for 78 top-ranking differentially methylated regions and tested on two independent cohorts of 132 and 177 samples, respectively. Further, a new statistical framework was developed to identify a robust set of MSP assays of which the methylation score (i.e. the percentage of methylated assays) allows accurate outcome prediction. Survival analyses were performed on the individual target level, as well as on the combined multimarker signature. As a result of the unbiased differential DNA methylation assessment by MBD sequencing, 58 of the 78 MSP assays were designed in regions previously unexplored in neuroblastoma, and 36 are located in non-promoter or non-coding regions. In total, 5 individual MSP assays (located in CCDC177, NXPH1, lnc-MRPL3-2, lnc-TREX1-1 and one on a region from chromosome 8 with no further annotation) predict event-free survival and 4 additional assays (located in SPRED3, TNFAIP2, NPM2 and CYYR1) also predict overall survival. Furthermore, a robust 58-marker methylation signature predicting overall and event-free survival was established. In conclusion, this study encompasses the largest unbiased DNA methylation biomarker study in neuroblastoma so far. We identified and independently validated several novel prognostic biomarkers, as well as a prognostic 58-marker methylation signature.

Zipper plot: visualizing transcriptional activity of genomic regions

Abstract: Summary: Reconstructing transcript models from RNA-sequencing (RNA-seq) data and establishing these as independent transcriptional units can be a challenging task. The Zipper plot is an application that enables users to interrogate putative transcription start sites (TSSs) in relation to various features that are indicative for transcriptional activity. These features are obtained from publicly available datasets including CAGE-sequencing (CAGE-seq), ChIP-sequencing (ChIP-seq) for histone marks and DNase-sequencing (DNase-seq). The Zipper plot application requires three input fields (chromosome, genomic coordinate (hg19) of the TSS and strand) and generates a report that includes a detailed summary table, a Zipper plot and several statistics derived from this plot. Availability and Implementation: The Zipper plot is implemented using the statistical programming language R and is freely available at http://zipperplot.cmgg.be Contact: Pieter.Mestdagh@UGent.be; Katleen.DePreter@UGent.be; Francisco.AvilaCobos@UGent.be Supplementary information: Supplementary Methods available online.