Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas
Genome Biology volume 15, Article number: 471 (2014)
The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behavior and can thereby complement the prognostically favorable 1p/19q co-deletion.
Genome-wide, 50 base pair single-end sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analyzed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors.
CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.
Diffuse low-grade gliomas (LGGs) are regarded as slow growing malignant brain tumors and patients can live up to 30 years with this disease. In a subset of patients the tumor exerts a more aggressive behavior and survival can be as short as two years . Personalized timing of postoperative treatment is crucial to forestall progression in the latter group whilst preventing long-term side-effects for patients with more favorable prospects . The disease course of patients with LGGs is correlated with gene mutations, such as in p53 and IDH1, hypermethylation of MGMT as well as chromosomal copy number aberrations (CNAs). Regarding the latter, assessment of combined loss of 1p and 19q currently is implemented in routine clinical care in specific glioma subgroups given its favorable prognostic and predictive value ,. Other CNAs, such as losses of chromosomes 10 and 11p, have been reported to be prognostically unfavorable, but have not been introduced into clinical practice yet, possibly due the limited number of samples included in the studies and/or lack of validation -. Unfavorable events might go undetected as a consequence of intratumoral heterogeneity in gliomas , which is particularly salient if they are only present in the more malignant subclones of LGGs that seed outgrowth of a recurrent tumor , thereby promoting a large extent of resection. As current knowledge on the temporal and spatial evolution of CNAs in LGGs is limited, we evaluated CNAs in a clinically and histologically representative cohort of formalin-fixed archival samples using shallow whole genome sequencing (shallow WGS). We demonstrate that loss of part or whole chromosome 10q is prognostically unfavorable and often present in a subclonal manner.
Clinical and histological data
We studied 173 formalin-fixed paraffin-embedded (FFPE) samples from 98 LGG patients, including spatially distinct regions and paired recurrent tumors, by shallow WGS. Patients had either deceased or had passed the median survival time of six years for LGGs. Other inclusion criteria and patient characteristics of this discovery cohort are summarized in Figure 1 and Table 1. Age at diagnosis, overall survival and postoperative treatment (type and timing) varied extensively between patients, but not between the five participating hospitals, which contributed nearly equal numbers of cases. Comparison of overall survival between patients treated immediately after surgery and those for whom postoperative treatment was withheld did not reveal statistically significant differences. Characteristics of the LGG patients of the French validation (n =126)  and confirmation cohorts (n =184), the latter from publicly available data from The Cancer Genome Atlas (TCGA) , are also listed in Table 1. Due to the retrospective character of this study, the cohorts are not matched; there is considerable variation in duration of follow-up and the percentage of patients deceased and for which information on overall survival is available, as well as in the distribution of histological subgroups.
Copy number detection by shallow WGS in LGGs
To obtain genome-wide copy numbers from the FFPE samples of our discovery cohort, we evaluated the use of shallow WGS. First, for sample LGG284 a paired-end 100 (PE100) sequence run was performed. Copy number profiles were produced by counting the unique sequence tags per 15 kb bin of the paired-end 100 bp reads from both ends (PE100 in Figure S1A in Additional file 1), the single 100 bp read from one end (SR100 in Figure S1B in Additional file 1) and the trimmed first single 50 bp read from the same end (SR50 in Figure S1C in Additional file 1). The noise (measured as variance) of the different profiles is very similar and CNAs observed are indistinguishable from each other, which implies that the uniqueness of the 50 bp sequence tags suffices to infer copy number levels, and longer reads are not necessary. Array comparative genomic hybridization (array CGH) was performed on the same DNA sample, which confirmed the CNAs detected (Figure S1D in Additional file 1). For an additional eight samples both 50 bp single-read (SR50) shallow WGS and array CGH were applied as technical validation. Shallow WGS and array analysis invariably yielded the same CNA profiles (Figure S2 in Additional file 1). Based on this information, all subsequent analyses were performed using 50 bp single-read (SR50) shallow WGS since it is more cost-effective and allows the use of samples with short DNA fragments, which are frequently obtained with FFPE materials. The most frequent CNAs, detected in more than 10% of cases, are whole or partial loss of chromosomal arms 9p, 10q, 12p, 13 and 14, as well as gain of chromosomal arms 7q, 8q, 10p and 11q. The most frequent CNAs in this cohort are co-deletion of 1p and 19q often accompanied by loss of whole chromosome 4, all commensurate with previous reports  (Figure 2).
The prognostic value of CNAs in discovery, validation and confirmation cohorts
Association of survival with CNAs detected in the discovery cohort was tested. In addition to the known prognostically favorable 1p/19q co-deletion, five further chromosomal losses at chromosomes 9p, distal 10q, 11p, 13q, and 22q presented with statistical significance (Table 2). No associations were observed with gains. Significant regions were verified in the French validation cohort of 126 diffuse LGG patients (Table 1). Loss of distal 10q was an unfavorable CNA in both cohorts, whereas losses of chromosomal regions at 9p, 11p, 13q and 22q were not substantiated in the validation cohort (Table 2). In the discovery cohort, median overall survival for patients with or without loss of whole or distal 10q was respectively 6.6 years versus 16.7 years (18/98, P-value =0.009). The size of chromosome 10 deletion varies from whole chromosome loss (5/18) to 22.5 Mbp distal loss (10q25.2- 10qter). An association between loss of this region with overall survival was finally tested in the TCGA dataset of LGG. Despite the limited number of patients deceased in this cohort (Table 1), a significant association with overall survival was observed (P-value =0.0018) (Figure 3), which confirms that distal 10q is a prognostically unfavorable chromosomal aberration.
In the discovery cohort, absence of IDH1 or IDH2 mutation (11/98) was overrepresented in patients with distal 10q loss (7/11; five with whole chromosome 10 loss and two with distal 10q loss). After splitting the cohort by IDH status, a trend for distal 10q loss was observed; in the IDH mutant subgroup (n =87) the log rank test for loss of 10q (n =11) yielded a P-value of 0.077, and in the IDH wild-type subgroup (n =11), a similar P-value of 0.068 was yielded through the test for distal loss of 10q (n =7).
Co-deletion of 1p/19q was predominantly detected in LGGs with oligodendroglial histological features while loss of distal 10q was more frequently identified in astrocytic LGGs. However, there was no one-to-one relationship between histological features and these CNAs (Figure 2). Co-deletion of 1p/19q combined with distal 10q loss was observed in three LGGs of the discovery cohort and four LGGs of the validation cohort (all with oligodendroglial features) and none in the confirmation cohort. This limited number of patients does not allow for proper statistical survival analysis, but median survival of the patients in these cohorts combined (13.4 years) suggests that loss of 1p/19q and distal 10q counteracts overall survival. Simultaneous testing of both CNAs classified LGG patients with a favorable (1p/19q co-deletion), unfavorable (distal 10q loss), or intermediate (both) prognosis in all three cohorts (Figure 3A,B,C). In the discovery cohort, hazard ratios of 1p/19q co-deletion without distal 10q loss and of distal 10q loss without 1p/19q co-deletion were 0.30 (95% confidence interval 0.15 to 0.58), and 2.91 (95% confidence interval 1.53 to 5.55), respectively (Table 3).
Intratumoral heterogeneity of CNAs in LGGs
In addition to the above-mentioned CNAs detected in single samples, we studied intratumoral heterogeneity by shallow WGS of multiple, spatially distinct regions obtained during the same surgery. Among other CNAs, distribution of 10q loss was assessed, illustrated for LGG240 in Figure 4 (more examples are provided in Additional file 2). In the original chromosomal profile of LGG240, marginal deflection of 10q was observed (a smaller distance from the 0-line than for the losses in 1p, 4, or 19q (Figure 4A) . Assuming clonality of the 1p/19q co-deletion , this difference in extent of copy number loss suggests that 10q loss would only be present in about 30 to 35% of the tumor cells (Figure S3 in Additional file 1).
To further delineate intratumoral heterogeneity of CNAs in this sample, the originally outlined area was divided into three sub-regions and an additional tumor region within the same paraffin section was included (Figure 4B). The 1p/19q co-deletion as well as chromosome 4 loss were present in all sub-regions and assumed to be clonally present. Losses of chromosomes 9, 10, 13, 15, 18 and gain of chromosome 11 were present in one or few sub-regions and assumed to be heterogeneously present (Figure 4C). To technically validate the intratumoral copy number heterogeneity observed in LGG240, array CGH was performed for all but one (insufficient amount of DNA) of the spatially distinct regions, which confirmed either clonality of 1p/19q and chromosome 4 losses or heterogeneity of all six chromosomally aberrant regions (Figure S4 in Additional file 1). Three additional samples with a clonal type of deflection and four with a marginal deflection of (distal) 10q were technically validated by array CGH (Figure S2 in Additional file 1). Intratumoral heterogeneity was detected in 15 out of 17 LGGs analyzed for this purpose; 68% of the CNAs (84/124) were not homogeneously present in spatially distinct regions obtained during the same surgery, such as loss of chromosomal arm 5q, chromosome 13 and gain of 11p (Figure 5A). Co-deletion of 1p and 19q was the only CNA that was consistently present in all spatially distinct regions of LGGs with this combination of CNAs; others, such as gain of chromosomal arm 7q, were most often, but not always, clonal. Loss of 10q was heterogeneously present in seven out of eight patients (Figure 5B). Histological variability did not correspond to the extent of heterogeneity.
Temporal evolution of CNAs in LGGs
Forty-seven out of 98 patients were subjected to a second surgery because of tumor progression. Of 20 patients, 24 recurrent tumors could be retrieved from medical archives. Almost 50% of CNAs (99/207) in the initial and paired recurrent tumors were shared and 15% (31/207) were uniquely detected in the initial tumor. A substantial proportion (37%, 77/207) of CNAs was uniquely identified in the recurrent tumor, such as loss of genomic regions at chromosomes 4, 10 and 15 (Figures 5C and 6). 1p/19q co-deletion was consistently identified in initial as well as recurrent tumors and there were no cases with new 1p/19q co-deletion. In four patients, de novo loss of 10q (including distal 10q losses) surfaced in the recurrence. In hindsight, marginal deflection of 10q was observed in the initial tumor of one of these four patients, and was not detected by the calling algorithm . In two out of four patients with new loss of 10q a higher malignancy grade (WHO grade III or IV) had been assigned to the recurrent tumor (Figure 5D). In one of the three patients for which both spatially distinct regions of the initial tumor and recurrences were analyzed, subclonal 10q loss was present in one of the regions of the initial tumor, but undetectable in the recurrence (Additional file 3). In the other two patients, 10q loss was detected in both the initial and paired recurrent tumor.
Intratumoral heterogeneity at the genomic level has been observed in numerous types of cancer, although its implications for treatment often remain undetermined. In the present study, archival LGG material and matched clinical data provide insight into spatial and temporal evolution of prognostically relevant CNAs. All LGG samples collected for this study were included for copy number profiling by shallow WGS combined with a depth of coverage approach, yielding high quality data without technical dropouts. This approach proved to be particularly beneficial for our study, since no matched normal DNA is required, which is a major advantage when analyzing long-term archived FFPE tumor samples. While shallow WGS cannot detect copy-neutral loss of heterozygosity or rearrangements, it is cost-effective, with a quality comparable to array CGH and applicable to DNA isolated from the FFPE samples. This allowed us to include samples that had been archived for over 30 years and collate a representative cohort, including LGG patients with long survival.
The relatively low incidence of LGGs and relatively long overall survival of patients necessitated this retrospective, multi-center approach. The observed variability in postoperative treatment can be attributed to the lack of a standard of treatment for these patients. Despite these variations, distal 10q loss (including whole chromosome losses) was significantly associated with an unfavorable prognosis in the discovery, validation and confirmation cohorts. Previously, some studies with smaller cohorts of specific histological subgroups of LGG have reported a correlation between 10q and survival ,. Furthermore, a high prevalence of whole chromosome 10 loss and strong negative correlation with survival have been reported for grade III and IV gliomas -. Partial loss of 10q is much more frequently detected in histological grade II diffuse gliomas compared with grade III and IV gliomas. In each of the previously published studies the whole of chromosome 10 or the entire 10q arm was taken into account. Here we demonstrate that, different from higher grade gliomas, the distal end of 10q is frequently lost and associates with overall survival in three cohorts. Spatial as well as temporal analyses suggest that subclones with distal loss of 10q are involved in tumor progression, since the loss surfaces in paired recurrent tumors with a higher malignancy grade.
Identification of the genes and their proteins affected by CNAs may elucidate the biological underpinnings of their clinical relevance but can be challenging, especially when a genomic region is large. Only after many years were mutations in CIC and FUBP1 associated with co-deletion of chromosomal arms 1p and 19q ,. In total 148 genes are located on 10q25.2-qter, including MGMT, DMBT1 and ERCC6,,, while the usual suspect, PTEN, is located more proximal to the centromere  and is preferentially lost in higher grade gliomas .
Based on our findings, we suggest that patients with an LGG should be simultaneously tested for both 1p/19q co-deletion and distal loss of 10q, since these two phenomena seem to have counteractive effects on survival. Introduction of heterogeneous CNAs, such as distal 10q loss, in daily clinical practice requires a robust diagnostic test. The well-known clonal features of 1p/19q co-deletion have been helpful to interpret these intermediate copy number levels in LGGs . Analysis of multiple spatially distinct regions could reveal subclones. We currently favor genome-wide analysis, which visualizes chromosomes 1p, distal 10q and 19q simultaneously, and at the same time may provide insight into intratumoral heterogeneity within one region. Both extent of resection as well as the subclonal character of important markers for progression command alternative diagnostic procedures to assess their presence in a postsurgical situation, which may in the future be offered through peripheral blood screening . Meanwhile, detailed registration of the positions of samples from different regions within a tumor obtained during the same surgery may provide more accurate insight into biologically relevant topics such as the physical distance and direction of growth of tumor subclones as well as the overall extent of heterogeneity of a tumor .
Copy number analysis by shallow WGS is a robust approach for archival clinical LGG specimens. For a large proportion of LGG patients, analysis of CNAs with prognostic value may improve personalized timing of therapy. Thereby, loss of distal 10q without 1p/19q co-deletion is indicative for urgent postoperative treatment, while in LGGs without loss of 10q and with 1p/19q co-deletion a wait-and-scan policy should be considered. The subclonal character of whole or distal 10q loss in a subset of samples emphasizes the need for maximal extent of resection, illustrates that single sample diagnostics may be insufficient for LGG and favors future studies on genome-wide analysis of multiple spatially distinct samples to map tumor progression.
Materials and methods
Clinical data and sample collection for discovery, validation and confirmation cohorts
Approval for collection of clinical data and FFPE tumor samples for the 98 patients of the discovery cohort was obtained from the institutional review boards of all five Dutch hospitals, namely the Medical Ethical Committee (in Dutch: Medisch-Ethische Toetsingsingscommissie or METc) of the Academical Medical Center (AMC), the METc of the Isala klinieken in Zwolle, the METc of the VU University Medical Center (VUmc) in Amsterdam, the METc of the St Elisabeth Hospital in Tilburg, and the METc of the Arnhem - Nijmegen Region for samples from Radboud University Medical Center in Nijmegen (CMO). Experimental methods in this manuscript are in compliance with the Helsinki Declaration. Inclusion criteria and characteristics of the discovery cohort are summarized in Figure 1 and Table 1. Clinical features of the validation cohort, from a French hospital, can be found in Table 1; materials and methods are presented in more detail by Alentorn et al. . For the confirmation cohort, copy number data of 531 lower grade glioma patients from the TCGA database were downloaded on 12 June 2014 via the Cancer Browser at UCSC . Clinical data were available for 373 of these patients, including grade and overall survival; 184 of these samples were categorized as `diffuse glioma histological grade 2' and selected as a confirmation cohort . Presumably as a consequence of the fact that only fresh frozen samples were included in the TCGA cohort, a limited number of patients are contained in the dataset that had deceased during follow-up.
Histological revision of samples in the discovery cohort was performed by two experienced neuropathologists (EA and PW). For all samples in the discovery cohort, including paired recurrent tumors, areas containing >60% tumor cells were outlined on hematoxylin and eosin stained slides, and tumor cell percentage estimated and registered for each sample (Additional files 2 and 3 and Table S1 in Additional file 4) and 10 adjacent sections were used for DNA isolation . For the assessment of intratumoral heterogeneity, spatially distinct regions were selected based on histological variability and/or plain physical distance (Additional file 4). These samples were obtained either from one FFPE block, or individual blocks from the same surgery (Table S1 in Additional file 4). DNA (500 ng) was fragmented by sonication (Covaris™ S2, Woburn, MA, USA), and sequenced using a 50 bp single-read (50 bp SR) modus (Illumina TruSeq DNA-kit and HiSeq 2000, San Diego, CA, USA). The 100 bp paired-end (100 bp PE) sequencing modus and array CGH were used for comparison and technical validation. Array CGH was performed as described previously . IDH1 and IDH2 mutation analysis was performed as described previously .
Copy number data from shallow WGS were analyzed using a novel Bioconductor script called QDNAseq . QDNAseq infers copy numbers through depth of coverage by binning reads uniquely aligned to the human reference genome build GRCh37/hg19 with Burrow's Wheeler Alignment (BWA) . PCR duplicates and reads with mapping qualities below 37 (highest value returned by BWA) were filtered. Copy numbers were inferred from the number of sequence reads per 15 kb bin. A simultaneous Loess correction for sequence mappability and GC content is applied within QDNAseq, which reduces noise of the copy number profiles, particularly for those with more degraded DNA. Problematic genome regions were furthermore filtered by applying our procedures to sequence data from the 1000 Genomes Project  to obtain a blacklist that eliminates problematic regions and the most common copy number variants of germ-line origin. Sequence data as well as all array CGH data have been uploaded to the European Genome-phenome Archive (EGA; accession number EGAS00001000643).
Calling of CNAs into discreet categories (normal, gain or loss) for the discovery and validation set was performed with the Bioconductor/R-package CGHcall . A weighted hierarchical clustering of the CNAs was performed using call probabilities to assess similarity of chromosomal profiles . Association with survival was tested using a log rank test with significance estimated over 10,000 permutations. After discovery of regions of interest for survival, consecutive regions in the same chromosome with P-values <0.05 were fused together to final regions and the log rank test was repeated. Chromosomal regions that were still significant in the discovery set after multiple testing correction according to Benjamini-Hochberg were verified in the independent French validation cohort. Therefore, genomic coordinates were converted to the NCBI35/hg17 genome build using the UCSC liftOver tool  and P-values were calculated with the log rank test and adjusted with the more stringent Holm-Bonferroni method. The statistical significance of a CNA was calculated compared to the rest of the cohort not bearing this CNA - for example, samples with loss of distal 10q versus samples without loss of distal 10q.
Regions significant in both the discovery and validation cohorts were tested in the TCGA confirmation cohort for which CNA data were generated with Affymetrix SNP 6.0 arrays (Santa Clara, CA, USA). TCGA level 3 copy number data were publicly available at the time of download and mapped to NCBI36/hg18. These level 3 data involve beginning and end positions of chromosomal segments with deflection values, resulting from TCGA preprocessing (for level definitions and preprocessing see ). Segment values were converted to CNA discreet categories by setting thresholds whereby a log2 ratio of >0.20 is gain, < -0.23 is loss and all other values are normal copy number; these values correspond to 30% of the tumor cells with that CNA. Those segments overlapping for at least 90% of the 1p, 19q (excluding centromeres) or 10q25.2-qter region (corresponding to NCBI35/hg17: chr10: 112939991-135323881 and NCBI36/hg18: chr10:112939991-135284990) were taken into consideration. At this setting 14 patients had a distal 10q loss, of which 4 deceased during follow-up and no patients had both distal 10q and 1p/19q loss. P-value calculations were performed as described above without corrections since only one region was tested for confirmation. The threshold settings were selected based on the fact that a deletion in 30% of the tumor cells, as observed for the chromosome 10 loss in LGG sample 240 (Figure S3 in Additional file 1), should not be missed. All other threshold values for calling losses of these regions were stepwise tested as well as the percentage of overlap with the 10q25.2-qter region and are presented in Table S2 in Additional file 4. Significance for the 10q loss remained for many different settings, but the number of patients with this loss substantially decreased with widening margins for calling CNAs to lower than should be expected based on the discovery and validation cohorts.
To assess the presence of CNAs between spatially distinct regions and/or recurrences from the same patient, common regions were detected with CGHregions , and regions smaller than 5 Mbp were excluded. Clinical parameters were analyzed with log rank test. All reported P-values were two-sided, and <0.05 was considered statistically significant.
Both array CGH and sequence data have been uploaded to the European Genome-phenome Archive (EGA; accession number EGAS00001000643).
HFT: study conceptualization and design, sample and clinical data collection, DNA isolation, statistical analysis and manuscript writing. IS: study conceptualization and design, sequence data generation and interpretation, bioinformatics analysis, manuscript writing. HFE: DNA isolation, sequence library preparations, and manuscript review. AA: validation cohort, data collection, analysis, interpretation and manuscript review. DS and MC: sequence data generation and bioinformatics, manuscript review. RF: sample collection, histological analysis and manuscript review. AMG: clinical data collection and manuscript review. GB: clinical data collection and manuscript review. WAB: clinical data collection and manuscript review. GAM: histological analysis and manuscript review. MH: sample collection, histological analysis and manuscript review. AI: validation cohort, clinical data collection, study design, data interpretation and manuscript review. KHX: validation cohort clinical data collection and manuscript review. KM: histological analysis and manuscript review. RGV: confirmation cohort analysis and manuscript review. PV: study conceptualization, histological analysis and manuscript review. MAW: statistical analysis, bioinformatics and manuscript review. JJH: study conceptualization and design, manuscript writing. EA: histological analysis and revision, manuscript review. JCR: study conceptualization and design, clinical data collection and analysis, manuscript writing. PW: study conceptualization and design, sample collection, histological analysis and revision, manuscript writing. BY: study conceptualization and design, sequence data generation, sample collection, bioinformatics, manuscript writing. All authors read and approved the final manuscript.
The present study was performed by a wide variety of professionals. HFT is a PhD student in neuro-oncology and resident in neurology; IS, DS and MC are PhD students specialized in bioinformatics; HE is a research technician specialized in DNA isolation from FFPE, microarrays and next generation sequencing; AA, AMG, AI, KHX, JJH and JCR are neuro-oncologists providing clinical care for brain tumor patients and are involved in several research projects in this field; GB and WAB are neurosurgeons that provide daily clinical care for brain tumor patients; RF, MH, EA, GAM, PV and PW are pathologists with a special interest in neuro-oncology; MW is a mathematician and statistical expert in genome-wide data analysis; RGV and BY are molecular biologists specialized in tumor genome analysis.
comparative genomic hybridization
copy number aberration
diffuse low-grade glioma
polymerase chain reaction
The Cancer Genome Atlas
whole genome sequencing
Claus EB, Black PM: Survival rates and patterns of care for patients diagnosed with supratentorial low-grade gliomas: data from the SEER program, 1973-2001. Cancer. 2006, 106: 1358-1363. 10.1002/cncr.21733.
Douw L, Klein M, Fagel SS, van den Heuvel J, Taphoorn MJ, Aaronson NK, Postma TJ, Vandertop WP, Mooij JJ, Boerman RH, Beute GN, Sluimer JD, Slotman BJ, Reijneveld JC, Heimans JJ: Cognitive and radiological effects of radiotherapy in patients with low-grade glioma: long-term follow-up. Lancet Neurol. 2009, 8: 810-818. 10.1016/S1474-4422(09)70204-2.
Erdem-Eraslan L, Gravendeel LA, de Rooi J, Eilers PH, Idbaih A, Spliet WG, den Dunnen WF, Teepen JL, Wesseling P, Sillevis Smit PA, Kros JM, Gorlia T, van den Bent MJ, French PJ: Intrinsic molecular subtypes of glioma are prognostic and predict benefit from adjuvant procarbazine, lomustine, and vincristine chemotherapy in combination with other prognostic factors in anaplastic oligodendroglial brain tumors: a report from EORTC study 26951. J Clin Oncol. 2013, 31: 328-336. 10.1200/JCO.2012.44.1444.
van Thuijl HF, Ylstra B, Wurdinger T, van Nieuwenhuizen D, Heimans JJ, Wesseling P, Reijneveld JC: Genetics and pharmacogenomics of diffuse gliomas. Pharmacol Ther. 2013, 137: 78-88. 10.1016/j.pharmthera.2012.09.002.
Alentorn A, van Thuijl HF, Marie Y, Alshehhi H, Carpentier C, Boisselier B, Laigle-Donadey F, Mokhtari K, Scheinin I, Wesseling P, Ylstra B, Capelle L, Hoang-Xuan K, Delattre JY, Reijneveld JC, Idbaih A: Clinical value of chromosome arms 19q and 11p losses in low-grade gliomas. Neuro Oncol. 2014, 16: 400-408. 10.1093/neuonc/not227.
Bissola L, Eoli M, Pollo B, Merciai BM, Silvani A, Salsano E, Maccagnano C, Bruzzone MG, Fruhman Conti AM, Solero CL, Giombini S, Broggi G, Boiardi A, Finocchiaro G: Association of chromosome 10 losses and negative prognosis in oligoastrocytomas. Ann Neurol. 2002, 52: 842-845. 10.1002/ana.10405.
Houillier C, Mokhtari K, Carpentier C, Criniere E, Marie Y, Rousseau A, Kaloshi G, Dehais C, Laffaire J, Laigle-Donadey F, Hoang-Xuan K, Sanson M, Delattre JY: Chromosome 9p and 10q losses predict unfavorable outcome in low-grade gliomas. Neuro Oncol. 2010, 12: 2-6. 10.1093/neuonc/nop002.
Snuderl M, Fazlollahi L, Le LP, Nitta M, Zhelyazkova BH, Davidson CJ, Akhavanfard S, Cahill DP, Aldape KD, Betensky RA, Louis DN, Iafrate AJ: Mosaic amplification of multiple receptor tyrosine kinase genes in glioblastoma. Cancer Cell. 2011, 20: 810-817. 10.1016/j.ccr.2011.11.005.
Johnson BE, Mazor T, Hong C, Barnes M, Aihara K, McLean CY, Fouse SD, Yamamoto S, Ueda H, Tatsuno K, Asthana S, Jalbert LE, Nelson SJ, Bollen AW, Gustafson WC, Charron E, Weiss WA, Smirnov IV, Song JS, Olshen AB, Cha S, Zhao Y, Moore RA, Mungall AJ, Jones SJ, Hirst M, Marra MA, Saito N, Aburatani H, Mukasa A, et al: Mutational analysis reveals the origin and therapy-driven evolution of recurrent glioma. Science. 2014, 343: 189-193. 10.1126/science.1239947.
Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, Zheng S, Chakravarty D, Sanborn JZ, Berman SH, Beroukhim R, Bernard B, Wu CJ, Genovese G, Shmulevich I, Barnholtz-Sloan J, Zou L, Vegesna R, Shukla SA, Ciriello G, Yung WK, Zhang W, Sougnez C, Mikkelsen T, Aldape K, Bigner DD, van Meir EG, Prados M, Sloan A, Black KL, et al: The somatic genomic landscape of glioblastoma. Cell. 2013, 155: 462-477. 10.1016/j.cell.2013.09.034.
Idbaih A, Criniere E, Ligon KL, Delattre O, Delattre JY: Array-based genomics in glioma research. Brain Pathol. 2010, 20: 28-38. 10.1111/j.1750-3639.2009.00274.x.
van de Wiel MA, Kim KI, Vosse SJ, van Wieringen WN, Wilting SM, Ylstra B: CGHcall: calling aberrations for array CGH tumor profiles. Bioinformatics. 2007, 23: 892-894. 10.1093/bioinformatics/btm030.
Jeuken JW, Sijben A, Bleeker FE, Boots-Sprenger SH, Rijntjes J, Gijtenbeek JM, Mueller W, Wesseling P: The nature and timing of specific copy number changes in the course of molecular progression in diffuse gliomas: further elucidation of their genetic `life story'. Brain Pathol. 2011, 21: 308-320. 10.1111/j.1750-3639.2010.00447.x.
Horbinski C, Nikiforova MN, Hobbs J, Bortoluzzi S, Cieply K, Dacic S, Hamilton RL: The importance of 10q status in an outcomes-based comparison between 1p/19q fluorescence in situ hybridization and polymerase chain reaction-based microsatellite loss of heterozygosity analysis of oligodendrogliomas. J Neuropathol Exp Neurol. 2012, 71: 73-82. 10.1097/NEN.0b013e318240fa65.
Ramirez C, Bowman C, Maurage CA, Dubois F, Blond S, Porchet N, Escande F: Loss of 1p, 19q, and 10q heterozygosity prospectively predicts prognosis of oligodendroglial tumors–towards individualized tumor treatment?. Neuro Oncol. 2010, 12: 490-499.
Sasaki H, Zlatescu MC, Betensky RA, Ino Y, Cairncross JG, Louis DN: PTEN is a target of chromosome 10q loss in anaplastic oligodendrogliomas and PTEN alterations are associated with poor prognosis. Am J Pathol. 2001, 159: 359-367. 10.1016/S0002-9440(10)61702-6.
Thiessen B, Maguire JA, McNeil K, Huntsman D, Martin MA, Horsman D: Loss of heterozygosity for loci on chromosome arms 1p and 10q in oligodendroglial tumors: relationship to outcome and chemosensitivity. J Neurooncol. 2003, 64: 271-278. 10.1023/A:1025689004046.
Bettegowda C, Agrawal N, Jiao Y, Sausen M, Wood LD, Hruban RH, Rodriguez FJ, Cahill DP, McLendon R, Riggins G, Velculescu VE, Oba-Shinjo SM, Marie SK, Vogelstein B, Bigner D, Yan H, Papadopoulos N, Kinzler KW: Mutations in CIC and FUBP1 contribute to human oligodendroglioma. Science. 2011, 333: 1453-1455. 10.1126/science.1210557.
Yip S, Butterfield YS, Morozova O, Chittaranjan S, Blough MD, An J, Birol I, Chesnelong C, Chiu R, Chuah E, Corbett R, Docking R, Firme M, Hirst M, Jackman S, Karsan A, Li H, Louis DN, Maslova A, Moore R, Moradian A, Mungall KL, Perizzolo M, Qian J, Roldan G, Smith EE, Tamura-Wells J, Thiessen N, Varhol R, Weiss S, et al: Concurrent CIC mutations, IDH mutations, and 1p/19q loss distinguish oligodendrogliomas from other cancers. J Pathol. 2012, 226: 7-16. 10.1002/path.2995.
Boulay JL, Ionescu MC, Sivasankaran B, Labuhn M, Dolder-Schlienger B, Taylor E, Morin P, Hemmings BA, Lino MM, Jones G, Maier D, Merlo A: The 10q25.3-26.1 G protein-coupled receptor gene GPR26 is epigenetically silenced in human gliomas. Int J Oncol. 2009, 35: 1123-1131. 10.3892/ijo_00000428.
Ramalho-Carvalho J, Pires M, Lisboa S, Graca I, Rocha P, Barros-Silva JD, Sawa-Bordalo J, Mauricio J, Resende M, Teixeira MR, Honavar M, Henrique R, Jeronimo C: Altered expression of MGMT in high-grade gliomas results from the combined effect of epigenetic and genetic aberrations. PLoS One. 2013, 8: e58206-10.1371/journal.pone.0058206.
Sottoriva A, Spiteri I, Piccirillo SG, Touloumis A, Collins VP, Marioni JC, Curtis C, Watts C, Tavare S: Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Proc Natl Acad Sci U S A. 2013, 110: 4009-4014. 10.1073/pnas.1219747110.
Zhu J, Sanborn JZ, Benz S, Szeto C, Hsu F, Kuhn RM, Karolchnik D, Archi J, Lenburg ME, Esserman LJ, Kent WJ, Haussler D, Wangt T: The UCSC Cancer Genomics Browser. Nat Methods. 2009, 6: 239-240. 10.1038/nmeth0409-239.
van Essen HF, Ylstra B: High-resolution copy number profiling by array CGH using DNA isolated from formalin-fixed, paraffin-embedded tissues. Methods Mol Biol. 2012, 838: 329-341. 10.1007/978-1-61779-507-7_16.
Krijgsman O, Israeli D, van Essen HF, Eijk PP, Berens ML, Mellink CH, Nieuwint AW, Weiss MM, Steenbergen RD, Meijer GA, Ylstra B: Detection limits of DNA copy number alterations in heterogeneous cell populations. Cell Oncol (Dordr). 2013, 36: 27-36. 10.1007/s13402-012-0108-2.
Scheinin I, Sie D, Bengtsson H, van de Wiel MA, Olshen AB, van Thuijl HF, van Essen HF, Eijk PP, Rustenburg F, Meijer GA, Reijneveld JC, Wesseling P, Pinkel D, Albertson DG, Ylstra B: DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.Genome Res 2014, gr.175141.114.,
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.
Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, Kang HM, Marth GT, McVean GA: An integrated map of genetic variation from 1,092 human genomes. Nature. 2012, 491: 56-65. 10.1038/nature11632.
van Wieringen WN, van de Wiel MA, Ylstra B: Weighted clustering of called array CGH data. Biostatistics. 2008, 9: 484-500. 10.1093/biostatistics/kxm048.
Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, Hillman-Jackson J, Kuhn RM, Pedersen JS, Pohl A, Raney BJ, Rosenbloom KR, Siepel A, Smith KE, Sugnet CW, Sultan-Qurraie A, Thomas DJ, Trumbower H, Weber RJ, Weirauch M, Zweig AS, Haussler D, Kent WJ: The UCSC Genome Browser Database: Update 2006. Nucleids Acids Res. 2006, 1: D590-8. 10.1093/nar/gkj144.
The Cancer Genome Atlas Data Portal , [https://tcga-data.nci.nih.gov/tcga]
van de Wiel MA, Wieringen WN: CGHregions: dimension reduction for array CGH data with minimal information loss. Cancer Inform. 2007, 3: 55-63.
This work was supported by grants from the Dutch Cancer Society (KWF kankerbestrijding, grant nr 2009-4470), foundation STOPHersentumoren.nl and Edli foundation. The validation cohort is part of the national program Cartes d’Identité des Tumeurs® (CIT). The research leading to these results received funding from the program `Investissements d’avenir’ ANR-10-IAIHU-06.
AI: Novartis, Hoffmann-La Roche, IntselChimos and BetaInnov, no conflict of interest related to the present work. KHX: LOC network (primary CNS lymphoma) supported by the INCa (National Institute of Cancer). AA: Obra Social `la Caixa' and ARTC. The other authors report no competing interests.
Electronic supplementary material
Additional file 1: Figure S1.: Evaluation of shallow WGS for genome-wide copy number analysis for read lengths of sample LGG 284. Copy number profiles were produced by counting the number of uniquely mapped sequence tags per 15 kb bin for different settings. (A) Paired-end 100 (PE100), 100 bp reads from both ends. (B) Single-end 100 (SR100), 100 bp reads. (C) Single-end 50 (SR50), 50 bp reads from one end. (D) Technical validation by array CGH of the same sample. Y-axis, normalized log 2 sequence read counts per bin; x-axis, 15 kb bins ordered by genomic position from chromosomes 1 to22. Figure S2. Technical validation of shallow WGS by array CGH CNA analysis of eight LGGs tested on both platforms. The level of deflection of CNAs was comparable. Also, heterogeneously present CNAs with only marginal deflection could be reproduced; LGGs 168, 184 and 193 show clonal distal 10q loss, while LGGs 187, 189, 210, 211 and 224 show subclonal loss of 10q. Figure S3. Minimal tumor cell percentage required for detection of CNAs. Calculation of the percentage of tumor cells with a chromosome 10 deletion in LGG240. To simulate a limited tumor percentage, a virtual tumor-normal mixture experiment was performed with LGG240 (diffuse low-grade glioma) and NA18960 ((1000 Genomes Project Consortium). Figure S4. Technical validation of shallow WGS by array CGH. Technical validation of LGG240 (original) and three spatially distinct regions. (PDF 10 MB)
Additional file 2: Copy number profiles generated with shallow WGS of original LGGs of the discovery cohort and spatially distinct regions of this tumor obtained during the same surgery. A representative picture of histology (hematoxylin and eosin staining, original magnification × 200) is depicted in the top left corner of each profile. (PDF 14 MB)
Additional file 3: Copy number profiles generated with shallow WGS of initial LGGs and paired recurrent tumors. A representative picture of histology (hematoxylin and eosin staining, original magnification × 200) is depicted in the top left corner of each profile. (PDF 10 MB)
Additional file 4: Table S1.: Overview of all samples of the discovery cohort, including paired recurrent tumors and spatially distinct regions. Columns show ID number, relative spatial distance (one FFPE block or individual FFPE block), histology, prognostically relevant CNAs, distal 10q loss, 1p/19q co-deletion, read count, sequence depth and tumor cell percentage as estimated by neuropathologists. Table S2A-C. Matrix with P-values for overall survival of the TCGA LGG (grade 2) confirmation cohort at various thresholds. Vertical rows: log2 ratio thresholds for calling CNA discrete categories (loss, normal or gain) from the log2 ratio level 3 data downloaded from the Cancer Browser at UCSC . Horizontal rows: thresholds for overlap with the 10q25.2-qter region. Grey, non-significant settings; shades of beige, significant settings (darker indicates higher signioficance); bold, threshold settings used for Figure 3. (A) P-values rounded to four decimal places for samples with or without a distal chromosome 10q loss. (B) Number of grade 2 LGG patients with clinical data (n =184) from the lower grade glioma TCGA dataset with a distal 10q deletion. (C) Number of grade 2 LGG patients (n =184) from the lower grade glioma TCGA dataset with a distal 10q deletion that had deceased during follow-up. Overall survival is significant for various settings of distal 10q loss. The number of patients with distal 10q loss substantially decreased with widening margins for calling CNAs, to lower than should be expected based on the French and Dutch cohorts. (XLSX 34 KB)
Authors’ original submitted files for images
About this article
Cite this article
van Thuijl, H.F., Scheinin, I., Sie, D. et al. Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas. Genome Biol 15, 471 (2014). https://doi.org/10.1186/s13059-014-0471-6