- Open Access
Anaplastic oligodendrogliomas with 1p19q codeletion have a proneural gene expression profile
Molecular Cancervolume 7, Article number: 41 (2008)
In high grade gliomas, 1p19q codeletion and EGFR amplification are mutually exclusive and predictive of dramatically different outcomes. We performed a microarray gene expression study of four high grade gliomas with 1p19q codeletion and nine with EGFR amplification, identified by CGH-array.
The two groups of gliomas exhibited very different gene expression profiles and were consistently distinguished by unsupervised clustering analysis. One of the most striking differences was the expression of normal brain genes by oligodendrogliomas with 1p19q codeletion. These gliomas harbored a gene expression profile that partially resembled the gene expression of normal brain samples, whereas gliomas with EGFR amplification expressed many genes in common with glioblastoma cancer stem cells. The differences between the two types of gliomas and the expression of neuronal genes in gliomas with 1p19q codeletion were both validated in an independent series of 16 gliomas using real-time RT-PCR with a set of 22 genes differentiating the two groups of gliomas (AKR1C3, ATOH8, BMP2, C20orf42, CCNB1, CDK2, CHI3L1, CTTNBP2, DCX, EGFR, GALNT13, GBP1, IGFBP2, IQGAP1, L1CAM, NCAM1, NOG, OLIG2, PDPN, PLAT, POSTN, RNF135). Immunohistochemical study of the most differentially expressed neuronal gene, alpha-internexin, clearly differentiated the two groups of gliomas, with 1p19q codeletion gliomas showing specific staining in tumor cells.
These findings provide evidence for neuronal differentiation in oligodendrogliomas with 1p19q codeletion and support the hypothesis that the cell of origin for gliomas with 1p19q codeletion could be a bi-potential progenitor cell, able to give rise to both neurons and oligodendrocytes.
The 1p19q codeletion and EGFR amplification are mutually exclusive and related to dramatically different outcomes in high grade gliomas. The 1p19q codeletion is strongly associated with an oligodendroglial phenotype and favorable prognosis . It has recently been shown to be mediated by a specific t(1;19)(q10;p10) translocation . To date the efforts performed to identify the genes specifically involved in the breakpoint have failed, mostly because both 1p and 19q centromeric regions contain highly repeated sequences. As a consequence the molecular mechanisms underlying the particular phenotype and the favorable outcome of this subset of gliomas remain completely unknown. Reliable detection of 1p19q codeletion requires an appropriate technique, such as CGH-array. Indeed, the most widely used LOH studies may not distinguish this signature from partial distal 1p and 19q deletion or gain, which have radically different prognostic implications . On the other hand, EGFR amplification is tightly associated with chromosome 10 loss and gain of chromosome 7, representing another characteristic genomic signature . EGFR amplification is more frequent in glioblastomas, but it is also found in a subset of anaplastic oligodendrogliomas and, in this setting, is predictive of extremely poor prognosis . Recently, malignant gliomas have been separated into three expression profiles with distinct outcomes and histological correlations: 1) the proneural profile with a better prognosis, mostly corresponding to anaplastic gliomas (oligodendrogliomas and astrocytomas); 2) the proliferative and 3) mesenchymal profiles, corresponding mainly to glioblastomas . However, correlation with 1p19q codeletion is still missing. Based on a set of gliomas analyzed by CGH-array , we selected tumors displaying one of these two characteristics and mutually exclusive patterns -1p19q codeletion or EGFR amplification- and compared their gene expression profiles.
The microarray study was done on 13 gliomas selected from the Salpêtrière database, based on the following criteria: 1) CGH-array profile showing either whole 1p19q codeletion or EGFR amplification, 2) high quality RNA availability. The samples were provided as snap-frozen sections of areas immediately adjacent to the region used for the histopathological diagnosis according to the World Health Organization Classification (WHO 2000). This set included 4 grade III oligodendrogliomas with complete 1p19q codeletion and 9 gliomas with EGFR amplification (5 glioblastomas (GBM), 3 grade III oligodendrogliomas, 1 grade III oligoastrocytoma (OAIII)). Genomic characterization was performed using CGH array as previously described . Among the 9 tumors with EGFR amplification, 8 out of 9 had chromosome 10q loss and chromosome 7 gain; 4 had a gain of chromosome 1p, and 4 had a partial loss of chromosome 1p. Among the 4 tumors with complete 1p19q codeletion none had EGFR amplification, 10q loss or chromosome 7 gain, and 2 had complete chromosome 4 loss. In order to compare the gene expression profile of the gliomas with normal brain, we used the gene expression data of 5 samples of corpus callosum (GSM175855, GSM175856, GSM175857, GSM175858, GSM176050) and 5 samples of cortex (GSM176049, GSM176344, GSM176345, GSM176346, GSM176347), available in the Gene Expression Omnibus repository (GSE7307) . To compare the gene expression profile with glioblastomas cancer stem cells (CSC), we used the data of Beier et al. (GSE7181) . All raw and normalized data files for the microarray analysis have been deposited  at the European Bioinformatics Institute (Hinxton, UK), and are publicy available under accession number E-MEXP-1507.
RNA extraction and hybridization
Approximately 50 mg of tissue from each tumor were used for total RNA extraction using the RNeasy Lipid Tissue mini kit (Qiagen, CA), according to the manufacturer instructions. RNA quality was verified with the Bioanalyser System (Agilent Technologies, Paolo Alto, CA), using the RNA Nano Chips. One and half micrograms of RNA were processed for hybridization on the Genechip Human Genome U133 Plus 2.0 Expression array (Affymetrix, CA), which contains over 54.000 probe sets analyzing the expression level of over 47,000 transcripts and variants, including 38.500 well-characterized human genes. The processing was done according to the recommendations of the manufacturer.
Except as indicated, all transcriptome analysis was carried out using either R-system software (version 2.4.1) packages including those of Bioconductor or original R code (A. de R.). Normalization was performed using the RMA method . Clustering analysis was performed as previously reported . Class comparison using a univariate t-test was performed using BRB Array Tools developed by Dr. Richard Simon and Amy Peng Lam . A p-value < 0.001 was used to define differentially expressed genes. Gene set enrichment analysis was performed using GSEA v2.0 software  as described by Subramanian et al. . For enrichment analysis in specific gene ontology terms (GO terms), we used a hypergeometric test to measure the association between a gene (probe set) list and a GO term. To this end, we mapped both the gene list and GO term related proteins to non-redundant Entrez Gene identifiers. The mapping of a probe set list to Entrez Gene ids was done using the annotation file HG-U133_Plus_2.annot.csv . For each GO term, we obtained the list of non-redundant related protein identifiers – either directly associated with the GO term or with one of its descendants – and mapped it to a non-redundant list of Entrez Gene ids. GO terms and their relationships (parent/child) were downloaded from The Gene Ontology site . The list of proteins associated with GO terms (table gene_association.goa_human) was downloaded from . We designated a threshold significance level for the hypergeometric test of p < 0.01 and the criteria that a GO term was represented by at least 2 Entrez Gene identifiers. Enrichment analysis of genes located on specific chromosomes was performed using DAVID tools .
Quantitative RT-PCR was performed as previously described . TBP (Genbank accession no. NM_003194), which encodes the TATA box-binding protein (a component of the DNA-binding protein complex transcription factor II D), was selected as an endogenous control because the levels of its transcript did not change across various normal tissues and tumor samples. The expression of the following 22 genes was quantified: AKR1C3 (upper primer (UP) 5'-CGT ATT TCA ACC GGA GTA AAT TGC TA-3', lower primer (LP) 5'-GTT CGG GTC CAC CCA TCG T-3'); ATOH8 (UP 5'-CAC ACC ATC AGC GCA GCC TT-3', LP 5'-GAT GGC CAG TTT GGA CAG CTT CT-3'); BMP2 (UP 5'-CGC AGC TTC CAC CAT GAA GAA TC-3', LP 5'-GAA TCT CCG GGT TGT TTT CCC ACT-3'); C20orf42 (UP 5'-AAG GAA CTT GAA CAA GGA GAA CCA CT-3', LP 5'-GGC ACA ACT TCG CAG CCT CTA-3'); CCNB1 (UP 5'-TGG ATA ATG GTG AAT GGA CAC CAA-3', LP 5'-GCC AGG TGC TGC ATA ACT GGA-3'); CDK2 (UP 5'-GGA CGG AGC TTG TTA TCG CAA AT-3', LP 5'-CCT TGG CCG AAA TCC GCT T-3'); CHI3L1 (UP 5'-GAC CAC AGG CCA TCA CAG TCC-3', LP 5'-TGT ACC CCA CAG CAT AGT CAG TGT T-3'); CTTNBP2 (UP 5'-CCC TCT CCA TCC TTG AAG CAG T-3', LP 5'-GAA GCT TCT CCA TTT CCA GCT TCT-3'); DCX (UP 5'-AGC CAA GAG CCC TGG TCC TAT-3', LP 5'-TGG AGG TTC CGT TTG CTG AGT-3');EGFR (UP 5'-GGA GAA CTG CCA GAA ACT GAC C-3', LP 5'-GCC TGC AGC ACA CTG GTT G-3');GALNT13 (UP 5'-GTG GCC TAT TTT CTA TTG ACA GAA ACT-3', LP 5'-CCT CCA CAT TGC CAA ATC CTA A-3');GBP1 (UP 5'-CCT CGC TCT TAA ACT TCA GGA ACA-3', LP 5'-CCT TTC GTC GTC TCA TTT TCG T-3');IGFBP2 (UP 5'-GGC CCT CTG GAG CAC CTC TAC T-3', LP 5'-CCG TTC AGA GAC ATC TTG CAC TGT-3');IQGAP1 (UP 5'-AGA ACA GAC CAG ATA CAA GGC GA-3', LP 5'-CTT AGG CAA TCC AAT CTC ATC CA-3');L1CAM (UP 5'-CTG GTT GTC TTC CCC ACA GAT GA-3', LP 5'-TCG TCC AGC GGA ACT GCA CT-3');NCAM1 (UP 5'-CCA CAG CCA TCC CAG CCA A-3', LP 5'-GAC GAT GAG GAT GCC CAC GAT-3');NOG (UP 5'-AAG AAG CAG CGC CTA AGC AAG A-3', LP 5'-GTC GTT CCA CGC GTA CAG CA-3');OLIG2 (UP 5'-CGC CAG AGC CCG ATG ACC TT-3', LP 5'-GAC ACG GTG CCC CCA GTG AA-3');PDPN (UP 5'-GTG ACT CCA GGA ACC AGC GAA-3', LP 5'-TGA CAC TTG TTG CCA CCA GAG TT-3');PLAT (UP 5'-AGC AGG CCC TGT ACT TCT CAG ATT-3', LP 5'-ACG TGG CCC TGG TAT CTA TTT CA-3');POSTN (UP 5'-GTC CTA ATT CCT GAT TCT GCC AAA-3', LP 5'-GGG CCA CAA GAT CCG TGA A-3');RNF135 (UP 5'-TGC CTG ACC AGA GCC ACC C-3', LP 5'-GAT GGA TGG CCC ACT GAG CA-3'); TBP (UP 5'-TGC ACA GGA GCC AAG AGT GAA-3', LP 5'-CAC ATC ACA GCT CCC CAC CA-3'). RT-PCR validation was done in an independent set of 16 gliomas (8 grade III oligodendrogliomas with complete 1p19q codeletion and 8 glioblastomas with EGFR amplification) and on 3 samples of normal brain (temporal lobe obtained during surgery for epilepsy).
Immunohistochemical analyses were carried out on paraffin sections using antibodies directed against internexin neuronal intermediate filament protein alpha (INA) . Immunostaining was studied in five glioblastomas with EGFR amplification and five anaplastic oligodendrogliomas with complete 1p19q codeletion.
Transcriptomic differences partly reflect the underlying genomic alterations
To study the relationship between the differences in gene expression profile and the underlying genomic alterations, we looked at the genomic localization of the differentially expressed genes. Class comparison using a univariate t-test demonstrated that 4458 probe sets were differentially expressed between the two groups of gliomas with a p-value < 0.001 and a maximum false discovery rate of 1.2% (see Additional file 1). The set of overexpressed genes in gliomas with EGFR amplification was significantly enriched in genes located on chromosome 1 (305 genes, p < 10-4), chromosome 19 (151 genes, p < 10-4), chromosome 7 (109 genes, p < 10-4) and chromosome 4 (85 genes, p < 10-4). In contrast, the set of overexpressed genes in gliomas with 1p19q codeletion was significantly enriched in genes located on chromosome 10 (216 genes, p < 10-4). Thus, gene expression and copy number variation dynamics were tightly correlated in both tumor groups. This was even more obvious when we plotted the differentially expressed probe sets localized on chromosome 1, 19 and 10 according to their genomic localization and the log2 ratio of their geometric mean in both tumor groups (Figure 1). Indeed, almost all differentially expressed genes localized on 1p and 19q were underexpressed in gliomas with 1p19q codeletion, whereas differentially expressed genes localized to 1q and 19p were equally distributed (Figure 1).
Unsupervised analysis consistently distinguishes the two groups of gliomas
In order to study the differences in gene expression profile, we performed unsupervised hierarchical clustering analysis of the 13 glioma samples. This analysis was done using the 1366 probe sets whose expression varied the most across samples (this corresponds to the probe sets with a robust coefficient of variation (rCV) superior to the 97.5th rCV percentile). As shown in Figure 2, gliomas with 1p19q codeletion and gliomas with EGFR amplification segregated into two distinct groups. This clustering was extremely robust and was conserved across different gene lists and clustering methods. The genes were classified into three clusters: one cluster of genes overexpressed in gliomas with EGFR amplification (gene cluster A, 698 genes), one cluster of genes overexpressed in oligodendrogliomas with 1p19q codeletion (gene cluster C, 488 genes) and a smaller cluster of genes expressed by some samples of both groups (gene cluster B, 180 genes). Enrichment analysis was performed on these three gene clusters. First, chromosome enrichment analysis demonstrated that these gene clusters did not simply reflect the underlying genomic alterations. Indeed, gene cluster A was enriched in genes located on chromosome 4 (53 genes, p < 10-4) and on chromosome 7 (43 genes, p < 10-2), but not in genes located on chromosome 1 or chromosome 19. Neither gene cluster C nor gene cluster B was enriched in genes located on chromosome 10. In contrast, the three gene clusters were enriched in genes with specific ontologic classes. Gene cluster A was most significantly enriched in genes involved in: immune response (55 genes, p < 10-4), extracellular matrix (37 genes, p < 10-4), proliferation (45 genes, p < 10-4), blood vessel development (17 genes, p < 10-4) and embryonic development (14 genes, p < 10-4). This cluster contained several genes well-known to be overexpressed in glioblastomas (IGFBP2, CHI3L1) and markers of glioblastoma cancer stem cells as well (CD133, IQGAP1). In contrast, gene cluster C was significantly enriched in genes with different specific ontologic classes: nervous system development (42 genes, p < 10-4), synaptic transmission (26 genes, p < 10-4) and neurogenesis (13 genes, p < 10-4). Actually, most of the genes in gene cluster C were either related to neuronal function or neuronal development or known to be highly expressed in normal brain. Gene cluster B was also enriched in genes involved in synaptic transmission (17 genes, p < 10-4) and nervous system development (16 genes, p < 10-4). However, two ontologic classes specifically found in this cluster (neurofilament (3 genes, p < 10-4) and axon ensheathment (2 genes, p < 10-2)) suggested that this cluster contained genes expressed in more differentiated neural cells than gene cluster C.
Oligodendrogliomas with 1p19q codeletion express specific subsets of neuronal genes
In order to better characterize the expression of neuronal genes in gliomas with complete 1p19q codeletion, we performed a new hierarchical clustering analysis with samples of normal brain (GSE7307), including grey matter (cortex) and white matter (corpus callosum). As glioblastomas expressed genes of neural cancer stem cells, we also included samples of glioblastoma cancer stem cells (GSE7181)  in this analysis. Still using the most differentially expressed 1366 probe sets, the samples clustered into two major groups: one containing 1p19q codeleted gliomas, normal white matter and normal grey matter, and the other containing EGFR amplified gliomas and cancer stem cells (Figure 3). Again this clustering was found to be extremely robust and was conserved across different gene lists and clustering methods. Gene ontology enrichment analysis was performed on the main gene clusters (A to J). Despite the clustering of 1p19q codeleted gliomas with normal brain, they showed substantial differences. First, the genes characteristics of the corpus callosum (gene cluster A), enriched in myelination genes (5 genes, p < 10-4), were not overexpressed in gliomas with 1p19q codeletion. Second, gliomas with 1p19q only overexpressed one subset of the neuronal genes characteristic of the cortex samples (gene clusters B and C). For example, among the neurofilament genes (INA, NEFH, NEFM, NEFL) present in cluster C, alphainternexin (INA) was the only one overexpressed in both gliomas with 1p19q codeletion and cortex samples in comparison to gliomas with EGFR amplification. Third, 1p19q codeleted gliomas were characterized by one specific cluster (gene cluster D) enriched in genes involved in CNS development (6 genes, p < 10-4), which also contained genes known to be specifically expressed in neuronal cells in the physiological state, e.g. DCX and GALNT1 [21, 22]. Thus, gliomas with 1p19q codeletion had a specific gene expression pattern of neuronal genes different from the samples of normal brain. Gliomas with EGFR amplification and glioblastoma cancer stem cells were both characterized by a large gene cluster (G) most significantly enriched in genes involved in proliferation (34 genes, p < 10-4) and in CNS development (16 genes, p < 10-4). Gliomas with EGFR amplification segregated from the cancer stem cells by one main gene cluster (I) enriched in genes involved in immune response (11 genes, p < 10-4), extracellular matrix (7 genes, p < 10-4) and angiogenesis (4 genes, p < 10-4).
Most characteristic genes associated with oligodendrogliomas with 1p19q codeletion
To find the genes most specifically associated with oligodendrogliomas with 1p19q codeletion, we selected the probe sets that were consistently (> 2-fold) and significantly (t test p < 0.001) overexpressed in these gliomas when they were independently compared to each of the 4 other samples groups (i.e., gliomas with EGFR amplification, cortex samples, corpus callosum samples and glioblastoma cancer stem cells). Eighty-six probe sets corresponding to 39 well-annotated genes met these criteria (Table 1 and see additional file 1). Several genes on this list are known to be highly expressed in normal brain (CSMD3, C20orf42, CTTNBP2), and one is known to be specifically expressed by neuronal cells (GALNT13) . Two transcription factors that play a role in CNS development were also specifically overexpressed (ATOH8, NFIB). ATOH8 is a basic-helix-loop-helix transcription factor, whose homolog in mouse has been demonstrated to regulate neuronal versus glial fate . NFIB plays a role in brain development, and Nfib deficient mice exhibit callosal agenesis . Finally, an intriguing feature was the specific overexpression of both BMP2, which promotes astroglial differentiation, and its antagonist NOG, which has been shown to promote both neuronal and oligodendroglial differentiation [25–27].
Oligodendrogliomas with 1p19q have a proneural gene expression profile
As gliomas with 1p19q codeletion expressed neuronal genes, we performed gene set enrichment analysis (GSEA) to study the relationship between the gene expression profile of these gliomas and the "proneural" gene expression signatures that have been described in high grade gliomas with a good prognosis (Figure 4) [5, 14, 28]. This test determines the over-representation of a gene set (list of genes) at the extremes (top or bottom) of the ordered, non-redundant dataset (list of all of the genes being used to compare two groups of samples). This analysis demonstrated that gliomas with 1p19q codeletion in comparison to gliomas with EGFR amplification were significantly enriched in the "proneural" gene set associated with good prognosis reported by Phillips et al. and in the good prognosis neurogenesis-related gene set reported by Freije et al. (HC1A gene set) [5, 28]. They were also enriched in the HC1B gene set (neuronal genes) of Freije et al. In contrast, gliomas with EGFR amplification were enriched in gene sets associated with poor prognosis ("proliferation" and "mesenchymal" gene sets of Phillips et al., and HC2A (enriched in proliferation genes) and HC2B (enriched in extracellular matrix genes) gene sets of Freije et al.) [5, 28]. Next we used the 35 genes signature developed by Phillips et al. to distinguish the three groups of high grade gliomas (proneural, proliferative and mesenchymal) in order to perform unsupervised hierarchical clustering . As shown in Figure 5, the gliomas with 1p19q codeletion were classified as proneural. Thus, there was a clear association between the "proneural" gene expression profile and 1p19q codeletion.
Real-time RT-PCR validation
To validate these findings, we studied the expression of 22 selected genes differentially expressed (11 up and 11 down) between the two groups of gliomas in an independent data set of 16 gliomas (8 gliomas with EGFR amplification and 8 gliomas with 1p19q codeletion). This study was performed in comparison with 3 samples of normal brain obtained from epileptic surgery (Table 2, Figure 6 and see additional file 2). Using a univariate t-test, 21 out of the 22 genes studied were shown to be differentially expressed between the two groups of gliomas with a p-value < 0.05 (only NCAM1 was not validated). This confirmed that gliomas with 1p19q codeletion overexpressed neuronal/normal brain genes (AKR1C3, C20orf42, CTTNBP2, L1CAM, GALNT13) as well as genes implicated in gliogenesis and neurogenesis (OLIG2, BMP2, NOG, DCX, ATOH8). Except L1CAM and AKR1C3, all of these genes were also overexpressed in comparison to the normal brain samples, including two genes (DCX, GALNT13) known to be exclusively expressed in neuronal cells. BMP2, NOG, C20orf42, GALNT13 and OLIG2 belong to the list of proneural genes reported by Phillips .
In the gliomas with EGFR amplification, we confirmed the overexpression of genes implicated in proliferation (CCNB1, CDK2), extracellular matrix remodeling (PLAT, POSTN), immune response (GBP1), cancer stem cell signaling (IQGAP1) as well as several genes known to be highly expressed in glioblastomas (IGFBP2, CHI3L1, PDPN). CCNB1, CDK2 belong to the proliferative gene list, and CHI3L1 and PDPN, to the mesenchymal gene list of Phillips .
Finally, to validate the expression of neuronal genes in gliomas with 1p19q codeletion at the protein level, we studied the expression of the internexin neuronal intermediate filament protein alpha (INA) which was one of the neuronal genes most overexpressed in these gliomas in comparison to gliomas with EGFR amplification (FC = 15, p < 0.001). INA is a class-IV neuronal intermediate filament involved in the morphogenesis of neurons . Immunostaining for INA was positive in all five oligodendrogliomas with 1p19q codeletion examined. Immunopositivity was observed in some normal infiltrated neurons but was mostly seen in a specific cytoplasmic perinuclear staining pattern in tumor cells (Figure 7). Between 20 to 50% of tumor cells displayed this staining, which was different from the staining observed in the infiltrated normal neurons (Figure 7). Among the five glioblastomas with EGFR amplification, immunostaining was negative in four and positive in a scattered pattern in one, in a region displaying some features of oligodendroglial differentiation.
EGFR amplification and whole 1p19q codeletion are mutually exclusive and predictive of completely different outcomes [3, 4]. To date, no studies have compared the gene expression profile of these two types of gliomas. Indeed, among microarray studies of gliomas [5, 28, 30, 31], only a few have compared genetically well-defined tumors [32–35]. In addition, these studies were based on LOH or FISH [32–34], and not on CGH-array. Yet, there is a need when interpreting a difference in gene expression to analyze it in relation to the genomic profile. Our data reveals clearly distinct gene expression profiles in these 2 groups of gliomas: those with EGFR amplification express the proliferative and mesenchymal gene set defined by Phillips et al., while 1p19q codeleted gliomas express the proneural group . Moreover, gliomas with EGFR amplification clustered close to tumor stem cells. Indeed, the EGFR pathway is involved in the proliferation of normal neural stem cells and cancer stem cells . This result is consistent with the fact that several studies isolated stem like tumor cells from glioblastoma but not from oligodendroglioma. Such studies did not include genetic profiles of the tumors, but data from our group suggest indeed that the capacity of cell renewal (as reflected by the formation of spheroids derived from the tumor) in vitro is tightly correlated with the presence of EGFR amplification (unpublished results). EGFR activation upregulates genes involved in neural stem cell proliferation: one of these genes is ASPM (abnormal spindle-like microcephaly associated) that promotes neuroblast proliferation and symmetric division and is strongly upregulated in glioblastomas. Inhibition of ASPM inhibits glioblastoma cell growth and neural stem cell proliferation .
A proneural/normal brain gene expression profile is a factor related to good prognosis and correlates with younger age and grade III histology, with most anaplastic oligodendrogliomas being classified as proneural [5, 28]. As shown here, this gene expression profile can be determined by a simple, highly discriminating RT-PCR test, and this may be useful for clinical practice. Until now, a proneural gene expression profile has not been reported to be associated with 1p19q codeletion. In Freije's study the number of gliomas with 1p19q codeletion was too small (4 out 74 patients) to address this question . In Phillips' study the genomic/transcriptomic correlation was limited to patients with astrocytoma histology, and this may have limited the possibility of finding an association between 1p19q codeletion and the proneural gene expression profile . However, the authors noticed a negative correlation between the proneural gene expression profile and EGFR amplification, similar to the negative correlation between 1p19q codeletion and EGFR amplification [3, 5]. Our study demonstrates that there is a strong correlation between 1p19q codeletion and the expression of proneural genes, suggesting that gliomas with a 1p19q codeletion represent a subgroup of proneural gliomas. In addition, the expression of neuronal genes in 1p19q codeleted tumors is consistent with a previous study showing selective expression of neuronal genes in oligodendrogliomas with 1p loss . Whether there is a link between the good prognosis of proneural gliomas and the fact that gliomas with 1p19q codeletion display a proneural gene expression profile remains to be elucidated. We make the hypothesis that gliomas without 1p19q codeletion but with a gene expression profile similar to the 1p19q codeleted gliomas might also harbor a better prognosis.
The expression of "neuronal genes" in 1p19q codeleted gliomas can be interpreted in different ways. As advocated by some authors, this expression is probably due in part to the presence of infiltrated neurons in the tumor . Indeed, 1p19q codeletion has been suggested to be more frequent in tumors with indistinct, irregular borders, which therefore, are more likely to be contaminated with normal brain tissue . However, as shown here, this normal brain infiltration cannot completely explain the expression of neuronal genes by 1p19q codeleted gliomas. Indeed, these tumors only express a specific subset of neuronal genes (Figure 3). In addition, if the expression of neuronal genes was only due to infiltration of normal brain tissue, the expression pattern of the neuronal genes in these tumors would be similar to their expression in the normal brain samples, which was not the case. Furthermore, we have demonstrated that alpha-internexin (INA), a neuronal protein, was specifically expressed by 1p19q codeleted glioma tumor cells. Thus INA expression might be used as a simple surrogate marker of 1p19q codeletion. This hypothesis is currently being tested in a larger series of gliomas. Interestingly, recent ultrastructural analysis of oligodendrogliomas has shown neuronal structures such as synapses and neurosecretory granules . Thus, another hypothesis for the expression of neuronal genes in 1p19q codeleted glioma tumor cells is that the cell of origin of these tumors could be a progenitor cell giving rise to both neurons and oligodendrocytes [40, 41]. This progenitor has less capacity of self renewal than the more multipotent neural stem cells. This is consistent with the fact that 1p19q codeleted oligodendroglioma fails in our hands to form spheroids in vitro (unpublished data). In this setting it is interesting to note that concomitant overexpression of both BMPs and BMP antagonists, such as the concomitant overexpression of BMP2 and NOG observed in 1p19q codeleted gliomas in our study, has been demonstrated in white matter progenitor cells, which can give rise to both oligodendrocytes and neurons . Another non-exclusive explanation for the expression of "neuronal" genes in oligodendrogliomas could rely on the fact that some genes involved in neurogenesis and classified as "neuronal" may also play a role in oligodendroglial development, e.g. ASCL1/MASH1. This proneural gene specifies a population of telencephalic oligodendrocytes  and is also required for oligodendrocyte development in the spinal cord . On the other hand, Olig2 -implicated in oligodendroglial specification- is also involved in neurogenesis: during development, Olig2+ progenitors give rise to both motoneurons and oligodendrocytes in the ventral spinal cord, . Consistently with our results, these data, illustrating the tight connection that exists between neurons and oligodendrocytes fates, bring a new light on the pathogenesis of oligodendrogliomas with 1p19q codeletion. Finally it is important to remember that current WHO classification is only based on morphological similarity between normal cells and tumor cells, and the link between oligodendrocytes and oligodendrogliomas has never been demonstrated.
Idbaih A, Marie Y, Pierron G, Brennetot C, Hoang-Xuan K, Kujas M, Mokhtari K, Sanson M, Lejeune J, Aurias A, Delattre O, Delattre JY: Two types of chromosome 1p losses with opposite significance in gliomas. Ann Neurol. 2005, 58 (3): 483-487. 10.1002/ana.20607
Jenkins RB, Blair H, Ballman KV, Giannini C, Arusell RM, Law M, Flynn H, Passe S, Felten S, Brown PD, Shaw EG, Buckner JC: A t(1;19)(q10;p10) mediates the combined deletions of 1p and 19q and predicts a better prognosis of patients with oligodendroglioma. Cancer Res. 2006, 66 (20): 9852-9861. 10.1158/0008-5472.CAN-06-1796
Idbaih A, Marie Y, Lucchesi C, Pierron G, Manie E, Raynal V, Mosseri V, Hoang-Xuan K, Kujas M, Brito I, Mokhtari K, Sanson M, Barillot E, Aurias A, Delattre JY, Delattre O: BAC array CGH distinguishes mutually exclusive alterations that define clinicogenetic subtypes of gliomas. Int J Cancer. 2008, 122 (8): 1778-1786. 10.1002/ijc.23270
Dehais C, Laigle-Donadey F, Marie Y, Kujas M, Lejeune J, Benouaich-Amiel A, Pedretti M, Polivka M, Xuan KH, Thillet J, Delattre JY, Sanson M: Prognostic stratification of patients with anaplastic gliomas according to genetic profile. Cancer. 2006, 107 (8): 1891-1897. 10.1002/cncr.22211
Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, Misra A, Nigro JM, Colman H, Soroceanu L, Williams PM, Modrusan Z, Feuerstein BG, Aldape K: Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell. 2006, 9 (3): 157-173. 10.1016/j.ccr.2006.02.019
Gene Expression Omnibus repository (GSE7307) (http://www.ncbi.nlm.nih.gov/geo).
Beier D, Hau P, Proescholdt M, Lohmeier A, Wischhusen J, Oefner PJ, Aigner L, Brawanski A, Bogdahn U, Beier CP: CD133(+) and CD133(-) glioblastoma-derived cancer stem cells show differential growth characteristics and molecular profiles. Cancer Res. 2007, 67 (9): 4010-4015. 10.1158/0008-5472.CAN-06-4180
Bioconductor: version 2.0, http://www.bioconductor.org.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15- 10.1093/nar/gng015
Boyault S, Rickman DS, de Reynies A, Balabaud C, Rebouissou S, Jeannot E, Herault A, Saric J, Belghiti J, Franco D, Bioulac-Sage P, Laurent-Puig P, Zucman-Rossi J: Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets. Hepatology. 2007, 45 (1): 42-52. 10.1002/hep.21467
Simon R, Peng-Lam A: BRB Array Tools (version 3.5.0), http://linus.nci.nih.gov/BRB-ArrayTools.html.
GSEA: v2.0 software. http://www.broad.mit.edu/gsea.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102
Affymetrix: annotation file HG-U133_Plus_2.annot.csv. http://www.affymetrix.com. 2007.
TheGeneOntology: http://www.geneontology.org. 2007.
GOA: Gene Ontology Database. ftp://ftp.ebi.ac.uk/pub/databases/GO/goa.
DAVID: National Institute of Allergy and Infectious Diseases (NIAID), NIH. http://david.abcc.ncifcrf.gov/. 2008.
Bieche I, Parfait B, Le Doussal V, Olivi M, Rio MC, Lidereau R, Vidaud M: Identification of CGA as a novel estrogen receptor-responsive gene in breast cancer: an outstanding candidate marker to predict the response to endocrine therapy. Cancer Res. 2001, 61 (4): 1652-1658.
Novocastra: http://www.leica-microsystems.com/Histology_Systems, product code : NCL-A-INTER, clone 2E3.
Couillard-Despres S, Winner B, Schaubeck S, Aigner R, Vroemen M, Weidner N, Bogdahn U, Winkler J, Kuhn HG, Aigner L: Doublecortin expression levels in adult brain reflect neurogenesis. Eur J Neurosci. 2005, 21 (1): 1-14. 10.1111/j.1460-9568.2004.03813.x
Zhang Y, Iwasaki H, Wang H, Kudo T, Kalka TB, Hennet T, Kubota T, Cheng L, Inaba N, Gotoh M, Togayachi A, Guo J, Hisatomi H, Nakajima K, Nishihara S, Nakamura M, Marth JD, Narimatsu H: Cloning and characterization of a new human UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase, designated pp-GalNAc-T13, that is specifically expressed in neurons and synthesizes GalNAc alpha-serine/threonine antigen. J Biol Chem. 2003, 278 (1): 573-584. 10.1074/jbc.M203094200
Inoue C, Bae SK, Takatsuka K, Inoue T, Bessho Y, Kageyama R: Math6, a bHLH gene expressed in the developing nervous system, regulates neuronal versus glial differentiation. Genes Cells. 2001, 6 (11): 977-986. 10.1046/j.1365-2443.2001.00476.x
Steele-Perkins G, Plachez C, Butz KG, Yang G, Bachurski CJ, Kinsman SL, Litwack ED, Richards LJ, Gronostajski RM: The transcription factor gene Nfib is essential for both lung maturation and brain development. Mol Cell Biol. 2005, 25 (2): 685-698. 10.1128/MCB.25.2.685-698.2005
Kasai M, Satoh K, Akiyama T: Wnt signaling regulates the sequential onset of neurogenesis and gliogenesis via induction of BMPs. Genes Cells. 2005, 10 (8): 777-783. 10.1111/j.1365-2443.2005.00876.x
Lim DA, Tramontin AD, Trevejo JM, Herrera DG, Garcia-Verdugo JM, Alvarez-Buylla A: Noggin antagonizes BMP signaling to create a niche for adult neurogenesis. Neuron. 2000, 28 (3): 713-726. 10.1016/S0896-6273(00)00148-3
Mabie PC, Mehler MF, Marmur R, Papavasiliou A, Song Q, Kessler JA: Bone morphogenetic proteins induce astroglial differentiation of oligodendroglial-astroglial progenitor cells. J Neurosci. 1997, 17 (11): 4112-4120.
Freije WA, Castro-Vargas FE, Fang Z, Horvath S, Cloughesy T, Liau LM, Mischel PS, Nelson SF: Gene expression profiling of gliomas strongly predicts survival. Cancer Res. 2004, 64 (18): 6503-6510. 10.1158/0008-5472.CAN-04-0452
Chan SO, Chiu FC: Cloning and developmental expression of human 66 kd neurofilament protein. Brain Res Mol Brain Res. 1995, 29 (1): 177-184. 10.1016/0169-328X(94)00268-J
Liang Y, Diehn M, Watson N, Bollen AW, Aldape KD, Nicholas MK, Lamborn KR, Berger MS, Botstein D, Brown PO, Israel MA: Gene expression profiling reveals molecularly and clinically distinct subtypes of glioblastoma multiforme. Proc Natl Acad Sci U S A. 2005, 102 (16): 5814-5819. 10.1073/pnas.0402870102
Nutt CL, Mani DR, Betensky RA, Tamayo P, Cairncross JG, Ladd C, Pohl U, Hartmann C, McLaughlin ME, Batchelor TT, Black PM, von Deimling A, Pomeroy SL, Golub TR, Louis DN: Gene expression-based classification of malignant gliomas correlates better with survival than histological classification. Cancer Res. 2003, 63 (7): 1602-1607.
French PJ, Swagemakers SM, Nagel JH, Kouwenhoven MC, Brouwer E, van der Spek P, Luider TM, Kros JM, van den Bent MJ, Sillevis Smitt PA: Gene expression profiles associated with treatment response in oligodendrogliomas. Cancer Res. 2005, 65 (24): 11335-11344. 10.1158/0008-5472.CAN-05-1886
Mukasa A, Ueki K, Ge X, Ishikawa S, Ide T, Fujimaki T, Nishikawa R, Asai A, Kirino T, Aburatani H: Selective expression of a subset of neuronal genes in oligodendroglioma with chromosome 1p loss. Brain Pathol. 2004, 14 (1): 34-42.
Mukasa A, Ueki K, Matsumoto S, Tsutsumi S, Nishikawa R, Fujimaki T, Asai A, Kirino T, Aburatani H: Distinction in gene expression profiles of oligodendrogliomas with and without allelic loss of 1p. Oncogene. 2002, 21 (25): 3961-3968. 10.1038/sj.onc.1205495
Tews B, Felsberg J, Hartmann C, Kunitz A, Hahn M, Toedt G, Neben K, Hummerich L, von Deimling A, Reifenberger G, Lichter P: Identification of novel oligodendroglioma-associated candidate tumor suppressor genes in 1p36 and 19q13 using microarray-based expression profiling. Int J Cancer. 2006, 119 (4): 792-800. 10.1002/ijc.21901
Lee J, Kotliarova S, Kotliarov Y, Li A, Su Q, Donin NM, Pastorino S, Purow BW, Christopher N, Zhang W, Park JK, Fine HA: Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closelymirror the phenotype and genotype of primary tumors than do serum-cultured celllines. Cancer Cell. 2006, 9 (5): 391-403. 10.1016/j.ccr.2006.03.030
Horvath S, Zhang B, Carlson M, Lu KV, Zhu S, Felciano RM, Laurance MF, Zhao W, Qi S, Chen Z, Lee Y, Scheck AC, Liau LM, Wu H, Geschwind DH, Febbo PG, Kornblum HI, Cloughesy TF, Nelson SF, Mischel PS: Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target. Proc Natl Acad Sci U S A. 2006, 103 (46): 17402-17407. 10.1073/pnas.0608396103
Jenkinson MD, du Plessis DG, Smith TS, Joyce KA, Warnke PC, Walker C: Histological growth patterns and genotype in oligodendroglial tumours:correlation with MRI features. Brain. 2006, 129: 1884-1891. 10.1093/brain/awl108
Vyberg M, Ulhøi BP, Teglbjaerg PS: Neuronal features of oligodendrogliomas-an ultrastructural and immunohistochemical study. Histopathology. 2007, 50: 887-896. 10.1111/j.1365-2559.2007.02686.x
Menn B, Garcia-Verdugo JM, Yaschine C, Gonzalez-Perez O, Rowitch D, Alvarez-Buylla A: Origin of oligodendrocytes in the subventricular zone of the adult brain. J Neurosci. 2006, 26 (30): 7907-7918. 10.1523/JNEUROSCI.1299-06.2006
Nunes MC, Roy NS, Keyoung HM, Goodman RR, McKhann G, Jiang L, Kang J, Nedergaard M, Goldman SA: Identification and isolation of multipotential neural progenitor cells from the subcortical white matter of the adult human brain. Nat Med. 2003, 9 (4): 439-447. 10.1038/nm837
Sim FJ, Lang JK, Waldau B, Roy NS, Schwartz TE, Pilcher WH, Chandross KJ, Natesan S, Merrill JE, Goldman SA: Complementary patterns of gene expression by human oligodendrocyte progenitors and their environment predict determinants of progenitor maintenance and differentiation. Ann Neurol. 2006, 56 (1): 773-9. 10.1200/JCO.2005.02.2194
Parras CM, Hunt C, Sugimori M, Nakafuku M, Rowitch D, Guillemot F: The proneural gene Mash1 specifies an early population of telencephalic oligodendrocytes. J Neurosci. 2007, 27 (16): 4233-4242. 10.1523/JNEUROSCI.0126-07.2007
Sugimori M, Nagao M, Parras CM, Nakatani H, Lebel M, Guillemot F, Nakafuku M: Ascl1 is required for oligodendrocyte development in the spinal cord. Development. 2008, 135 (7): 1271-1281. 10.1242/dev.015370
Lu QR, Sun T, Zhu Z, Ma N, Garcia M, Stiles CD, Rowitch DH: Common developmental requirement for Olig function indicates a motor neuron/oligodendrocyte connection. Cell. 2002, 109 (1): 75-86. 10.1016/S0092-8674(02)00678-5
This study was supported by a grant of the Institut National du Cancer (INCA, PL032) and by the Ligue Nationale contre le Cancer.
The authors declare that they have no competing interests.
FD performed the major part of experiments and analysis. FD and MS drafted the manuscript. AI performed the CGH-array study and analysis. AR and SL provided the bioinformatic tools and participated to the analysis. IB and MV performed the real-time RT-PCR validation. YM, SP, JT and KM selected the samples, helped extracting the RNA and performed the immunohistochemistry validation. KHX, OD and JYD assisted with design of the study and with critical examination of the manuscript. MS conceived of and designed the study, participated in its experimental design and interpretation of results, and helped edit the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Complete list of genes differentially expressed (t-test, p < 0.001) between 1p19q codeleted oligodendrogliomas (n = 4) and 1) gliomas with EGFR amplification (n = 9), 2) cerebral cortex samples (n = 5), 3) corpus callosum samples (n = 5) 4) Beier's et al. gliomas cancer stem cells (n = 6) . (XLS 4 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.