Evaluation of the similarity of gene expression data estimated with SAGE and Affymetrix GeneChips

BACKGROUND: Serial Analysis of Gene Expression (SAGE) and microarrays have found awidespread application, but much ambiguity exists regarding the evaluation of these technologies. Cross-platform utilization of gene expression data from the SAGE and microarray technology could reduce the need for duplicate experiments and facilitate a more extensive exchange of data within the research community. This requires a measure for the correspondence of the different gene expression platforms. To date, a number of cross-platform evaluations (including a few studies using SAGE and Affymetrix GeneChips) have been conducted showing a variable, but overall low, concordance. This study evaluates these overall measures and introduces the between-ratio difference as a concordance measure pergene. RESULTS: In this study, gene expression measurements of Unigene clusters represented by both Affymetrix GeneChips HG-U133A and SAGE were compared using two independent RNA samples. After matching of the data sets the final comparison contains a small data set of 1094 unique Unigene clusters, which is unbiased with respect to expression level. Different overall correlation approaches, like Up/Down classification, contingency tables and correlation coefficients were used to compare both platforms. In addition, we introduce a novel approach to compare two platforms based on the calculation of differences between expression ratios observed in each platform for each individual transcript. This approach results in a concordance measure per gene (with statistical probability value), as opposed to the commonly used overall concordance measures between platforms. CONCLUSION: We can conclude that intra-platform correlations are generally good, but that overall agreement between the two platforms is modest. This might be due to the binomially distributed sampling variation in SAGE tag counts, SAGE annotation errors and the intensity variation between probe sets of a single gene in Affymetrix GeneChips. We cannot identify or advice which platform performs better since both have their (dis)-advantages. Therefore it is strongly recommended to perform follow-up studies of interesting genes using additional techniques. The newly introduced between-ratio difference is a filtering-independent measure for between-platform concordance. Moreover, the between-ratio difference per gene can be used to detect transcripts with similar regulation on both platforms

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors

Sequencing of neuroblastoma identifies chromothripsis and defects in neuritogenesis genes

Neuroblastoma is a childhood tumour of the peripheral sympathetic nervous system. The pathogenesis has for a long time been quite enigmatic, as only very few gene defects were identified in this often lethal tumour. Frequently detected gene alterations are limited to MYCN amplification (20%) and ALK activations (7%). Here we present a whole-genome sequence analysis of 87 neuroblastoma of all stages. Few recurrent amino-acid-changing mutations were found. In contrast, analysis of structural defects identified a local shredding of chromosomes, known as chromothripsis, in 18% of high-stage neuroblastoma. These tumours are associated with a poor outcome. Structural alterations recurrently affected ODZ3, PTPRD and CSMD1, which are involved in neuronal growth cone stabilization. In addition, ATRX, TIAM1 and a series of regulators of the Rac/Rho pathway were mutated, further implicating defects in neuritogenesis in neuroblastoma. Most tumours with defects in these genes were aggressive high-stage neuroblastomas, but did not carry MYCN amplifications. The genomic landscape of neuroblastoma therefore reveals two novel molecular defects, chromothripsis and neuritogenesis gene alterations, which frequently occur in high-risk tumours

Joint Binding of OTX2 and MYC in Promotor Regions Is Associated with High Gene Expression in Medulloblastoma

Both OTX2 and MYC are important oncogenes in medulloblastoma, the most common malignant brain tumor in childhood. Much is known about MYC binding to promoter regions, but OTX2 binding is hardly investigated. We used ChIP-on-chip data to analyze the binding patterns of both transcription factors in D425 medulloblastoma cells. When combining the data for all promoter regions in the genome, OTX2 binding showed a remarkable bi-modal distribution pattern with peaks around −250 bp upstream and +650 bp downstream of the transcription start sites (TSSs). Indeed, 40.2% of all OTX2-bound TSSs had more than one significant OTX2-binding peak. This OTX2-binding pattern was very different from the TSS-centered single peak binding pattern observed for MYC and other known transcription factors. However, in individual promoter regions, OTX2 and MYC have a strong tendency to bind in proximity of each other. OTX2-binding sequences are depleted near TSSs in the genome, providing an explanation for the observed bi-modal distribution of OTX2 binding. This contrasts to the enrichment of E-box sequences at TSSs. Both OTX2 and MYC binding independently correlated with higher gene expression. Interestingly, genes of promoter regions with multiple OTX2 binding as well as MYC binding showed the highest expression levels in D425 cells and in primary medulloblastomas. Genes within this class of promoter regions were enriched for medulloblastoma and stem cell specific genes. Our data suggest an important functional interaction between OTX2 and MYC in regulating gene expression in medulloblastoma

Iron-responsive element of Divalent metal transporter 1 (Dmt1) controls Notch-mediated cell fates

Notch receptor activation is regulated by the intramembrane protease ?-secretase, which cleaves and liberates the Notch intracellular domain (Nicd) that regulates gene transcription. While ?-secretase cleavage is necessary, we demonstrate it is insufficient for Notch activation and requires vesicular trafficking. Here, we report Divalent metal transporter 1 (Dmt1, Slc11A2) as a novel and essential regulator of Notch signalling. Dmt1-deficient cells are defective in Notch signalling and have perturbed endolysosomal trafficking and function. Dmt1 encodes for two isoforms, with and without an iron response element (ire). We show that isoform-specific silencing of Dmt1-ire and Dmt1+ire have opposite consequences on Notch-dependent cell fates in cell lines and intestinal organoids. Loss of Dmt1-ire suppresses Notch activation and promotes differentiation, whereas loss of Dmt1+ire causes Notch activation and maintains stem-progenitor cell fates. Dmt1 isoform expression correlates with Notch and Wnt signalling in Apc-deficient intestinal organoids and human colorectal cancers. Consistently Dmt1-ire silencing induces Notch-dependent differentiation in colorectal cancer cells. These data identify Dmt1 isoforms as binary switches controlling Notch cell fate decisions in normal and tumour cells

RaS–MAPK pathway-driven tumor progression is associated with loss of CIC and other genomic aberrations in neuroblastoma

Mutations affecting the RAS–MAPK pathway frequently occur in relapsed neuroblastoma tumors, which suggests that activation of this pathway is associated with a more aggressive phenotype. To explore this hypothesis, we generated several model systems to define a neuroblastoma RAS–MAPK pathway signature. Activation of this pathway in primary tumors indeed correlated with poor survival and was associated with known activating mutations in ALK and other RAS–MAPK pathway genes. Integrative analysis showed that mutations in PHOX2B, CIC, and DMD were also associated with an activated RAS–MAPK pathway. Mutation of PHOX2B and deletion of CIC in neuroblastoma cell lines induced activation of the RAS–MAPK pathway. This activation was independent of phosphorylated ERK in CIC knockout systems. Furthermore, deletion of CIC caused a significant increase in tumor growth in vivo. These results show that the RAS–MAPK pathway is involved in tumor progression and establish CIC as a powerful tumor suppressor that functions downstream of this pathway in neuroblastoma. Significance: This work identifies CIC as a powerful tumor suppressor affecting the RAS-MAPK pathway in neuroblastoma and reinforces the importance of mutation-driven activation of this pathway in cancer

TERT rearrangements are frequent in neuroblastoma and identify aggressive tumors

Whole-genome sequencing detected structural rearrangements of TERT in 17 of 75 high-stage neuroblastomas, with five cases resulting from chromothripsis. Rearrangements were associated with increased TERT expression and targeted regions immediately up-and downstream of TERT, positioning a super-enhancer close to the breakpoints in seven cases. TERT rearrangements (23%), ATRX deletions (11%) and MYCN amplifications (37%) identify three almost non-overlapping groups of high-stage neuroblastoma, each associated with very poor prognosi

Full transcriptome analysis of rhabdomyosarcoma, normal, and fetal skeletal muscle: statistical comparison of multiple SAGE libraries

Rhabdomyosarcoma (RMS) is the most frequent soft tissue sarcoma in children. Improved treatment strategies have increased overall survival, but the response of approximately one-third of the patients is still poor. To increase the knowledge of RMS pathogenesis, we performed the first full transcriptome analysis of RMS using serial analysis of gene expression (SAGE). With a G-test for the simultaneous comparison of subsets of SAGE libraries of normal skeletal muscle, embryonal (ERMS) and alveolar (ARMS) RMS, we identified 251 differentially expressed genes. A literature-mining procedure demonstrated that 158 of these genes have not previously been associated with RMS or normal muscle. Gene Ontology (GO) analysis assigned 198 of the 251 genes to muscle-specific classes, including those involved in normal myogenic development, as well as tumor-related classes. Prominent GO classes were those associated with proliferation and actin reorganization, which are processes that play roles during early muscle development, muscle function, and tumor progression. Using custom microarrays, we confirmed the (up- or down-) regulation of 80% of 98 differentially expressed genes. Another SAGE library of 19- to 22-week-old fetal skeletal muscle was compared with the RMS and normal muscle transcriptomes. Cluster analysis showed that the RMS and fetal muscle SAGE libraries formed one cluster distinct from normal muscle samples. Moreover, the expression profile of 86% of the differentially expressed genes between normal muscle and RMS was highly similar in fetal muscle and RMS. In conclusion, the G-test is a robust tool for analyzing groups of SAGE libraries and correctly identifies genes marking the difference between fully differentiated skeletal muscle and RMS. This study not only substantiates the close association between embryonic myogenesis and RMS development but also provides a rich source of candidate genes to further elucidate the etiology of RMS or to identify diagnostic and/or prognostic marker

Specific and Sensitive Detection of Neuroblastoma mRNA Markers by Multiplex RT-qPCR

mRNA RT-qPCR is shown to be a very sensitive technique to detect minimal residual disease (MRD) in patients with neuroblastoma. Multiple mRNA markers are known to detect heterogeneous neuroblastoma cells in bone marrow (BM) or blood from patients. However, the limited volumes of BM and blood available can hamper the detection of multiple markers. To make optimal use of these samples, we developed a multiplex RT-qPCR for the detection of MRD in neuroblastoma. GUSB and PHOX2B were tested as single markers. The adrenergic markers TH, GAP43, CHRNA3 and DBH and mesenchymal markers POSTN, PRRX1 and FMO3 were tested in multiplex. Using control blood and BM, we established new thresholds for positivity. Comparison of multiplex and singleplex RT-qPCR results from 21 blood and 24 BM samples from neuroblastoma patients demonstrated a comparable sensitivity. With this multiplex RT-qPCR, we are able to test seven different neuroblastoma mRNA markers, which overcomes tumor heterogeneity and improves sensitivity of MRD detection, even in those samples of low RNA quantity. With resources and time being saved, reduction in sample volume and consumables can assist in the introduction of MRD by RT-qPCR into clinical practice