Distinct Genomic Aberrations between Low-Grade and High-Grade Gliomas of Chinese Patients

Abstract

Background

Glioma is a type of tumor that develops in the central nerve system, mainly the brain. Alterations of genomic sequence and sequence segments (such as copy number variations or CNV and copy neutral loss of heterozygosities or cnLOH) are thought to be a major determinant of the tumor grade.

Methods

We mapped genomic variations between low-grade and high-grade gliomas (LGG and HGG) in Chinese population based on Illumina’s Beadchip and validated the results using real-time qPCR.

Results

At the cytoband level, we discovered: (1) unique losses in LGG on 5q, 8p and 11q, and in HGG on 6q, 11p, 13q and 19q; (2) unique gains in the LGG on 1p and in HGG at 5p, 7p, 7q and 20q; and (3) cnLOH in HGG only on 3q, 8q, 10p, 14q, 15q, 17p, 17q, 18q and 21q. Subsequently, we confirmed well-characterized oncogenes among tumor-related loci (such as EGFR and KIT) and detected novel genes that gained chromosome sequences (such as AASS, HYAL4, NDUFA5 and SPAM1) in both LGG and HGG. In addition, we found gains, losses, and cnLOH in several genes, including VN1R2, VN1R4, and ZNF677, in multiple samples. Mapping grade-associated pathways and their related gene ontology (GO) terms, we classified LGG-associated functions as “arachidonic acid metabolism”, “DNA binding” and “regulation of DNA-dependent transcription” and the HGG-associated as “neuroactive ligand-receptor interaction”, “neuronal cell body” and “defense response to bacterium”.

Conclusion

LGG and HGG appear to have different molecular signatures in genomic variations and our results provide invaluable information for the diagnosis and treatment of gliomas in patients with variable duration or diverse tumor differentiation.

About the Authors:

Yunbo Li

Contributed equally to this work with: Yunbo Li, Dapeng Wang, Lei Wang

Affiliations First Hospital of Jilin University, Changchun, PR China, Jilin University, Changchun, PR China

Dapeng Wang

Contributed equally to this work with: Yunbo Li, Dapeng Wang, Lei Wang

* E-mail: [email protected] (DW); [email protected] (FZ); [email protected] (SF)

Affiliation: CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China

Lei Wang

Contributed equally to this work with: Yunbo Li, Dapeng Wang, Lei Wang

Current address: Life Technologies Corporation, Beijing, PR China

Affiliation: CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China

Jinhai Yu

Affiliation: First Hospital of Jilin University, Changchun, PR China

Danhua Du

Affiliation: First Hospital of Jilin University, Changchun, PR China

Ye Chen

Affiliation: Norman Bethune College of Medicine, Jilin University, Changchun, PR China

Peng Gao

Affiliation: First Hospital of Jilin University, Changchun, PR China

Duen-Mei Wang

Affiliation: CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China

Jun Yu

Affiliation: CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China

Feng Zhang

* E-mail: [email protected] (DW); [email protected] (FZ); [email protected] (SF)

Affiliations CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China, Northeast Normal University, Changchun, PR China

Shuanglin Fu

* E-mail: [email protected] (DW); [email protected] (FZ); [email protected] (SF)

Affiliation: First Hospital of Jilin University, Changchun, PR China

Introduction

Glioma is a major tumor type that derives from glial cells of the central nerve system, including of spine and the brain. Gliomas are classified into four grades, from I to IV, with increasing exacerbation according to histology [1]. Alternatively, they are categorized into ependymomas, astrocytomas, oligodendrogliomas, and mixed gliomas, according to the cell types that are anatomically defined from different brain regions based on the brain topology [2]. Glioma is considered as low prevalence but an increasing detection probability due to enhanced early detection techniques and procedures, which include X-ray, magnetic resonance imaging (MRI), and computed tomography (CT) [3]. Patients with gliomas often have a morbid state in various degrees, including pain, epilepsy, mental disorder, visual disturbance, hearing impairment, insomnia, and nausea [4]. Traditional therapies for the tumors, such as surgery, radiotherapy, chemotherapy, and oral medication, are often individually or selectively cooperated to remit or to cure the symptoms according to patient’s condition and tumor grade [5], [6].

Although the degree of copy number variations (CNVs) among healthy human populations do vary, abnormal increase of CNVs and cnLOHs are thought to be at least one of the important causative factors for gliomas and other cancers [7], [8], [9]. Large-scale genomic aberration and gene expression studies on gliomas have lead to identifications of genes, which are involved in numerous cellular functions and metabolic pathways (PCDH9, CXCL12, MYC, PDGFRA, PARK2, DMBT1, TOP2A, PTEN, ARF, TP53, P16, CDKN2B, RB1, EGFR, and NF1 [9], [10], [11], [12], [13], [14], [15], [16], [17]), as well as those related to neural development, cell signaling (RAS/RAF, RTK, MAPK, PI3K, and ROCK), and tumor suppression (p53 and RB) [1], [6], [11], [16], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28]. These genes serve as biomarkers useful for tumor classification, prognosis, and targeted therapies [17], [23], [29], [30], [31], [32], [33], [34].

Previous studies have been adding complexity for the understanding of gliomas at molecular levels. First, both nuclear and mitochondrial genomes show numerous structural variations [35]. Second, there are heterogeneous patterns of structural variations and/or expressions depending on the ways of classifying and grading gliomas, even between primary and secondary gliomas and between pediatric and adult gliomas [14], [20], [36], [37]. Third, consistency of the experimental results at DNA, RNA, and protein levels is poor due to complications of transcriptional and post-translational regulations, as well as enzymatic modifications [6], [12], [38], [39]. Fourth, relationship between sequence variations and prognosis are influenced by certain patient-specific parameters. For instance, older patients tend to have shorter survival time than younger patients [1], [40]. Fifth, co-variations can occur in different chromosomes and one gene may have multiple variations in different combinations [41], [42]. Last, relatively homogenous DNA variations have been found when topological locations (e.g., central vs. peripheral) are compared, where cross-contamination between normal and tumor tissues may contribute to the incorrect detection of structural variations [6], [43], [44]. Nevertheless, the existing literature related to genomic alterations and their molecular mechanisms in gliomas lacks information from Chinese populations.

Together with increasing throughput and efficiency of sequencing and microarray technologies for SNP and CNV discovery and typing, in-depth study on special tumor types become feasible, especially when two novel parameters–physical (Log R ratio) and genetic (B allele frequency) alterations/variability–are used to increase sensitivity and accuracy [45]. We report the first study in which we compared chromosomal aberrations of low-grade and high-grade gliomas (LGG and HGG) in Chinese populations to unveil differences of tumor-related sequence variations, including chromosomes and cytobands, and individual genes in the context of pathways and function categories.

Materials and Methods

Ethics Statement

This study was approved by the ethics committee at First Hospital of Jilin University, and all patients who agreed to this study signed their informed consents. In particular, some of the written informed consents were also obtained from the next of kin, caretakers, or guardians on the behalf of the minors/children participants involved in this study.

Sample Preparation and Genotyping

18 gliomas, including nine of each for both low-grade (grade II) and high-grade (6 grade III and 3 grade IV) tumors (as LGG and HGG), as well as 18 paired blood samples were collected at First Hospital of Jilin University (Table 1). In general, patients with LGG (average age of 41) and HGG (average age of 47) have no prior history of chemotherapy or radiation treatment, and the tumors and matched peripheral blood samples were obtained at the time of tumor resection. The blood samples were immediately snap-frozen in liquid nitrogen and stored at −80°C before DNA extraction. Portions of the solid tumors were treated in the same way as the blood samples and others were paraffin-embedded for histopathology. The genomic DNA was isolated by using QIAamp DNA Mini Kit (QIAGEN Inc). DNA quality and purity were determined based on agarose gel electrophoresis and spectrophotometry. Four additional gliomas for RNA RT-qPCR validation were also collected at the same hospital and treated according to the same protocol as used for the other samples (Table 1). Genotyping was carried out by using Illumina Sentrix Human CNV370-Quad v3 BeadChip according to the manufacturer’s instruction.

[Figure omitted. See PDF.]

Table 1. Overview of samples used in this study.

https://doi.org/10.1371/journal.pone.0057168.t001

Data Analysis and Functional Annotation

Illumina’s KaryoStudio (V1.0.3), together with the cnvPartition algorithm, was adopted to identify CNV regions. Illumina’s GenomeStudio software (V2009.1) with LOHscore plug-in was used to discover copy neutral loss of heterozygosity (cnLOH) regions. The samples were divided into LGG and HGG, and the analysis was limited to autosomal regions only due to the experiment design. The data used in this study were submitted to GEO with an accession of GSE34888.

GeneCoDis3 was used to relate genes to the corresponding Gene Ontology (GO) terms with cut-off = 0.05 [46], and WebGestalt V2 was used to map genes onto significant KEGG pathways with cut-off = 0.1 [47]. OMIM (Online Mendelian Inheritance in Man) was used to find relationship between genes and diseases [48].

Real-time qPCR-based Validation Experiment

Four RNA samples, reverse-transcribed (Table 1), were subjected to real-time qPCR in a final reaction volume of 10 µl, which contains: 1 µl of 10×buffer, 0.25 µl of dNTP, 2 µl of gene-specific forward and reverse primer mix (Table 2), 1 µl cDNA (100 ng/µl), 0.1 µl of rTaq, 0.5 µl of fluorochrome (TIANGEN BIOTECH) and 5.15 µl of ddH2O.

[Figure omitted. See PDF.]

Table 2. Primers and results of real-time qPCRs.

https://doi.org/10.1371/journal.pone.0057168.t002

Results

Case and Control

Using 18 solid tumors from patients (9 as LGG and 9 as HGG) as the case sample and the patients’ DNA isolated from their peripheral blood as the control, we compared chromosome aberrations between the case and control based on the Illumina’s BeadChip platform that defines both copy number variation (CNVs) and copy neutral loss of heterozygosity (cnLOH). This technology allows us to find sequence differences between tumor and blood samples at cytogenetic and molecular levels, which include both gains and losses. We only discovered sequence amplifications in the tumor not in the corresponding periphery blood. In other three aberration categories of CNVs, there are always more variations in the tumor than in the control, and the significance of variations between the paired samples can be evaluated statistically (P-value<0.05, paired t-test) when the total length influenced by CNVs in each sample was considered.

CNVs among Autosomal Arms

After the survey of different CNVs in 44 autosomal arms in LGG and HGG, we removed those shared by tumors and their controls to identify tumor-specific variations. We made a few interesting observations. First, 13q harbors a major high-frequency genomic aberration class. For instance, there are 3, 3, 3, 2, and 2 samples found to have homozygous deletion, hemizygous deletion, cnLOHs, duplication and amplification, respectively. Second, we found that 22, 29, 38, 38, and 6 chromosomal arms have variations in the five categories, respectively. Obviously, cnLOH is one of the most frequent categories, and our results agreed with a previous report [43]. In addition, the skewed distribution of sequence amplification showed that there were more copy number gains found on only a few chromosomal arms. Third, we paid more attention to CNVs between LGG and HGG and found several homozygous deletions of 13q in HGG but not in LGG. There were also cnLOHs in three HGG, such as those on 15q and 17q, but none was found in LGG. In addition, four HGG had duplications on 3q but none was found in LGG.

Comparison of LGG and HGG at the Cytoband Level

To pinpoint the position of variations on chromosomes, we investigated CNV and cnLOH in LGG and HGG. To reduce false positives, we only referred to those that occurred in at least two tumor samples (four samples for “duplications” category) after removing what appeared in the controls, i.e., corresponding blood samples (Table 3). Taking homozygous deletion as an example, we found several events on 8p11.23 in LGG, and 13q12.11 in HGG. In both hemizygous deletion and cnLOH, we found more cytobands in HGG than in LGG. Surprisingly, we found more duplication events in LGG than in HGG. Taken together, we conjectured that losses and gains occur more often in HGG and LGG, respectively. Furthermore, we found that gains usually occurred on 7q and 1p in HGG and LGG, respectively. Decreasing copy number is usually discovered on 6q, 13q and 19q in HGG. We also noticed that there were overlaps between the two tumor grades in the cnLOH category, especially the fact that cnLOHs were spread out more broadly in cross-chromosome cytobanding than the rest of copy number variation categories.

[Figure omitted. See PDF.]

Table 3. Comparison of sample numbers of genomic aberration in gliomas at the cytoband level.

https://doi.org/10.1371/journal.pone.0057168.t003

Genes Found to Associate with Gliomas

To identify specific genes associated with gliomas, we pooled all genomic aberrations occurred in at least six tumor samples. After filtered our data based on known variations found in the controls and Database of Genomic Variants (hg18.v8) (http://projects.tcag.ca/variation/), we had 24 genes and the related information was summarized (chromosome location, aberration category, tumor grading) (Table 4). These genes are all clustered on 1p, 7q, and 19q. Among them, 17 genes are only gains, and three of other genes, VN1R2, VN1R4 and ZNF677, have all three types of genomic aberrations–gain, loss and cnLOH. Referencing to the annotation of the OMIM Morbid Map (http://www.ncbi.nlm.nih.gov/omim), we found that AASS, TAS2R16 and TSPAN12 are previously identified to be disease-related and associated with “hyperlysinemia”, “alcohol dependence” and “exudative vitreoretinopathy”, respectively.

[Figure omitted. See PDF.]

Table 4. Selected genes involved in genomic aberration of gliomas.

https://doi.org/10.1371/journal.pone.0057168.t004

Pathway Enrichment and Functional Annotation of the Associated Genes

To decipher functional relevance of genes and gene classification related to gliomas, we established datasets based on the following two criteria: (i) elimination of genes shared by both the case and the control and (ii) removal of genes shared by less than four samples. Altogether, we collected 442 LGG and 111 HGG associated genes. After mapping the two groups of genes onto the KEGG pathways, we obtained 12 pathways for LGG, which belong to “lipid metabolism” (in the number of pathways: 6), “neurodegenerative diseases” (1), “endocrine system” (1), “nervous system” (1), “circulatory system” (1), “signal transduction” (1), and “immune system” (1) (Table 5). These genes were categorized in 10 pathways correlated to HGG, among which two were classified as “lipid metabolism” and others were distributed in diverse classifications such as “signaling molecules and interaction”, “signal transduction”, “endocrine system”, “carbohydrate metabolism”, “amino acid metabolism”, “glycan biosynthesis and metabolism”, and “infectious diseases”. The most four significant pathways identified in LGG are “arachidonic acid metabolism”, “linoleic acid metabolism”, “alpha-Linolenic acid metabolism”, and “ether lipid metabolism”. In HGG, the most four enriched terms are “metabolic pathways”, “neuroactive ligand-receptor interaction”, “calcium signaling pathway”, and “melanogenesis”. We found 12 and 4 enriched GO terms associated with LGG and HGG, respectively (Table 6). In LGG, the top two GO terms were “DNA binding”, and “regulation of transcription, DNA-dependent”. In HGG, the top two GO terms were “neuronal cell body”, and “defense response to bacterium”. To our current knowledge, tumors often alter cellular processes, such as proliferation, growth, programmed death, differentiation, division, mutation-induced DNA damage, and repair [1]. All these newly discovered distinct function categories may introduce features of primary (low-grade) or malignant (high-grade) tumors.

[Figure omitted. See PDF.]

Table 5. Pathway analysis of genes involved in genomic aberration in LGG (A) and HGG (B).

https://doi.org/10.1371/journal.pone.0057168.t005

[Figure omitted. See PDF.]

Table 6. GO enrichment analysis of genes involved in genomic aberration in LGG (A) and HGG (B).

https://doi.org/10.1371/journal.pone.0057168.t006

Correlating Copy Number Variation with Gene Expression

We used real-time qPCR to validate the expression of seven identified genes: AASS, CYP2J2, CYP4A11, PLA2G2A, PLA2G5, PTEN, and RB1. First, all genes showed increased expression when compared in the two tumor grades (HGG vs. LGG) (Table 2). CYP2J2, CYP4A11, PLA2G5 and PTEN exhibited even over 10 fold changes. We did not compare the gene expressions between the tumors and the controls, their corresponding blood cells, due to the tissue-specific nature of the gene expression. Of the seven genes, we noticed RB1 gains in HGG and losses in both grades. Other different aberrations include (1) CYP2J2, CYP4A11, PLA2G2A, and PLA2G5 gains in LGG, (2) AASS gain in both, and (3) PTEN gains in HGG only. The correlation between expression and copy number is complex, and has been demonstrated to be positive under most conditions with several exceptions of negative correlations [12]. There are a few interpretations for the complicated relations between gene (or segment) copy number and gene expression. First, the higher gene dosage leads to an increased transcript production and a positive correlation. Second, negative correlation has been proposed that some possible compensation or time-series feedback effects may act on the progression or cell response of tumors [12]. Last, a handful of CNVs can destroy regulatory regions to some extents and inhibit the expression of relevant genes. In any case, we have found significant differential gene expressions in selected genomic structures and genes between LGG and HGG.

Discussion

In accordance with Knudson’s two-hit hypothesis on tumor formation [49], tumorigenesis occurs concurrently with activation of oncogenes and inactivation of tumor suppressor genes (TSGs). In particular, tumor formation is a process where the potential for malignancy increases with mutation accumulation [50]. Our results support the concept that genomes in HGG tend to be deleted more intensively than those in LGG at both cytoband and molecular levels. In other words, there are gains of genes or partial sequences in the primary stage of the tumors and subsequently losses of tumor suppressor genes or TSGs appear to escape from normal cellular controls. For instance, cnLOHs in HGG happened in twelve chromosomes (2, 3, 6, 8, 9, 10, 14, 15, 17, 18, 21, and 22). In contrast, cnLOHs in LGG occurred only in four chromosomes (2, 6, 9, and 22) (Table 3). This particular clear observation in cnLOHs has not been found in other variation types, and the result indicates that disappearance of heterozygosity and allelic losses together with gene functions are major contributors to HGG.

We have effectively validated many locations, genes, pathways, and function categories in keeping with the previous studies at different levels. At the chromosome level, among the five variation classes, cnLOH is the major type, which is distributed over all chromosomes. Furthermore, we have validated many variations at the cytoband level (Table 3). For instance, we found duplications on 7q31.31, 7q31.32, 7q31.33, 7q33, 7q34, 7q35 and 7q36.1 in HGG, which is in agreement with gains of 7q reported in both grades [37]. We discovered hemizygous deletions on numerous cytobands (e.g., 6q12, 6q13, 6q14.1, 6q16.3, 6q21, 6q22.1, 6q23.1, 6q24.1, 13q12.2, 13q13.1, 13q14.12, 13q21.1, 19q12, and 19q13.11) in line with losses of 6q, 13q, and 19q in HGG [37], [44]. There have been some contradicting results in this study; for example, we only identified duplications on 1p in LGG and hemizygous deletions on 19q in HGG but LOHs on 1p-19q were reported in both grades by another group [41]. We also have some novel findings specific to Chinese populations and observed the complicated and altered roles of traditional tumor-related genes in our study. For example, KIT gains in 2 LGG and 2 HGG, MGMT gains in 1 LGG and 3 HGG, MXI1 gains in 1 LGG and 3 HGG, DMBT1 gains in 1 LGG and 3 HGG, IDH1 only gains in one LGG. Thus, TSGs and oncogenes are relative and condition-sensitive concept, i.e., TSGs in one sample may be oncogenes in another. We finally inferred that TSGs not only can be lost but also can be gained, whereas oncogenes appear generally being gained in these tumors.

We did not find any obvious major pathways shared by LGG and HGG but identified 6 and 2 subclasses from the major class “lipid metabolism” shared by the two grades, respectively. In HGG, we noticed some interesting variations, including AASS (4 gains) in “lysine degradation” and CHRM2 (4 gains), CYSLTR2 (2 gains and 2 losses), GRM8 (4 gains), HTR2A (2 gains and 2 losses), and MLNR (2 gains and 2 losses) in “neuroactive ligand-receptor interaction”. We further found variations of GRM8 (2 gains) and CYSLTR2 (1 loss), HTR2A (1 loss), MLNR (1 loss) in LGG. Gains in the four LGG are CYP2J2, CYP4A11, CYP4A22, PLA2G2A, PLA2G2C, PLA2G2D, PLA2G2E, PLA2G2F, and PLA2G5, and they are involved in “arachidonic acid metabolism”. Looking into well-known oncogenes and TSGs, we identified gains of EGFR in “de novo glioma pathway” (2 in LGG and 3 in HGG), as well as MDM2 (1 in LGG and 3 in HGG) and PTEN (1 in HGG), and losses of CDKN2A (2 in HGG). In “secondary glioma pathway”, gains of PDGFA (1 in LGG and 2 in HGG), PDGFRA (2 in LGG and 2 in HGG) and CDK4 (1 in LGG and 2 in HGG) and RB1 (2 in HGG), and losses of RB1 (1 in LGG and 2 in HGG) are also identified.

Our study was designed to uncover molecular differences between LGG and HGG, which are largely based on pathology, morphology, and degree of malignancy, and for which we know that altered genomic regions of HGG are more severe than those of LGG and that patient survival time is shorter in HGG (data not shown). Based on our study, we hope to associate genetics to clinical outcomes in three crucial aspects of fighting the disease, including diagnosis, treatment, and prognosis. First, specific and dominant genomic variations in cytobands or genes may be used to classify the two tumor grades at molecular level. Second, different genes with DNA aberrations may be further pursued as potential drug targets, especially in the personalized or population-specific cancer treatment for Chinese populations. Third, the disease-associated genes from LGG and HGG may be used for inferring survival time and guidance for disease treat [51], [52], [53]. In addition, we proposed a theory that losses of genes or alleles may play greater role in the formation and progression of HGG although both gains and losses contribute. This phenomenon is only observed at the cytoband level and relevant to concentrated alternations of homologous or heterogeneous genome composition [54]. Therefore, this is not a decisive but a substantial measure for glioma morbidity. Further large-scale investigations are needed for further validations.

Ultimately, we should be able to build an intricate network that is composed of pathways involved in gliomas both identified in this study and other known in the previous literature as we exemplified in Figure 1. The first characteristic of the network is the difficulty to fully discriminate all three complex modules of the pathways, such as LGG, HGG and others. However, we are able to identify major hubs among the modules, such as MP and MAPK, in the categories of HGG and others, respectively; some may even be obviously LGG associated, such as AAM. In addition, the tissue-specificity of gliomas is blurred by AP and other five vertices including PD, ALS, AD, HUD and PAD, which are common for neurodegenerative diseases and other neuro-pathological conditions. As a final note, the patterns of the three categories may suggest that there are differences in disease-causing mechanisms between Chinese and other populations. It is undoubted that studies on larger samples, coupled with in-depth analyses and validations, are what we should do in the near future.

[Figure omitted. See PDF.]

Figure 1. A network of pathways related to LGG, HGG, and others.

The abbreviations for KEGG pathways used in this analysis are listed as follows. AAM: Arachidonic acid metabolism; LAM: Linoleic acid metabolism; ALAM: alpha-Linolenic acid metabolism; ELM: Ether lipid metabolism; GPM: Glycerophospholipid metabolism; PD: Prion diseases; GNRH: GnRH signaling pathway; LTD: Long-term depression; VSMC: Vascular smooth muscle contraction; VEGF: VEGF signaling pathway; FCER: Fc epsilon RI signaling pathway; FAM: Fatty acid metabolism; MP: Metabolic pathways; NLI: Neuroactive ligand-receptor interaction; CSP: Calcium signaling pathway; ME: Melanogenesis; FMM: Fructose and mannose metabolism; LD: Lysine degradation; AEM: Androgen and estrogen metabolism; GM: Glycerolipid metabolism; GD: Glycosaminoglycan degradation; VCI: Vibrio cholerae infection; AD: Alzheimer's disease; PAD: Parkinson's disease; ALS: Amyotrophic lateral sclerosis; HUD: Huntington's disease; PI3K: PI3K-Akt signaling pathway; HIF-1: HIF-1 signaling pathway; MAPK: MAPK signaling pathway; CC: Cell cycle; P53: p53 signaling pathway; AP: Apoptosis; GL: Glioma; ERBB: ErbB signaling pathway; WNT: Wnt signaling pathway; and OP: Oxidative phosphorylation.

https://doi.org/10.1371/journal.pone.0057168.g001

Author Contributions

Conceived and designed the experiments: SLF FZ JY DMW. Performed the experiments: YBL LW FZ. Analyzed the data: DPW LW. Contributed reagents/materials/analysis tools: YBL JHY DHD YC PG FZ. Wrote the paper: DPW JY.

Citation: Li Y, Wang D, Wang L, Yu J, Du D, Chen Y, et al. (2013) Distinct Genomic Aberrations between Low-Grade and High-Grade Gliomas of Chinese Patients. PLoS ONE 8(2): e57168. https://doi.org/10.1371/journal.pone.0057168

References

1. Gladson CL, Prayson RA, Liu WM (2010) The pathobiology of glioma tumors. Annu Rev Pathol 5: 33–50.

2. Larjavaara S, Mantyla R, Salminen T, Haapasalo H, Raitanen J, et al. (2007) Incidence of gliomas by anatomic location. Neuro Oncol 9: 319–325.

3. Radhakrishnan K, Mokri B, Parisi JE, O'Fallon WM, Sunku J, et al. (1995) The trends in incidence of primary brain tumors in the population of Rochester, Minnesota. Ann Neurol 37: 67–73.

4. Cheng JX, Liu BL, Zhang X, Lin W, Zhang YQ, et al. (2010) Health-related quality of life in glioma patients in China. BMC Cancer 10: 305.

5. Taylor LP (2010) Diagnosis, treatment, and prognosis of glioma: five new things. Neurology 75: S28–32.

6. Parsons DW, Jones S, Zhang X, Lin JC, Leary RJ, et al. (2008) An integrated genomic analysis of human glioblastoma multiforme. Science 321: 1807–1812.

7. Blauw HM, Veldink JH, van Es MA, van Vught PW, Saris CG, et al. (2008) Copy-number variation in sporadic amyotrophic lateral sclerosis: a genome-wide screen. Lancet Neurol 7: 319–326.

8. Stark M, Hayward N (2007) Genome-wide loss of heterozygosity and copy number analysis in melanoma using high-density single-nucleotide polymorphism arrays. Cancer Res 67: 2632–2642.

9. Yin D, Ogawa S, Kawamata N, Tunici P, Finocchiaro G, et al. (2009) High-resolution genomic copy number profiling of glioblastoma multiforme by single nucleotide polymorphism DNA microarray. Mol Cancer Res 7: 665–677.

10. van den Boom J, Wolter M, Kuick R, Misek DE, Youkilis AS, et al. (2003) Characterization of gene expression profiles associated with glioma progression using oligonucleotide-based microarray analysis and real-time reverse transcription-polymerase chain reaction. Am J Pathol 163: 1033–1043.

11. The Cancer Genome Atlas Research Network (2008) Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 455: 1061–1068.

12. Kotliarov Y, Kotliarova S, Charong N, Li A, Walling J, et al. (2009) Correlation analysis between single-nucleotide polymorphism and expression arrays in gliomas identifies potentially relevant target genes. Cancer Res 69: 1596–1603.

13. Hu J, Jiang C, Ng HK, Pang JC, Tong CY, et al. (2003) Genome-wide allelotype study of primary glioblastoma multiforme. Chin Med J (Engl) 116: 577–583.

14. Paugh BS, Qu C, Jones C, Liu Z, Adamowicz-Brice M, et al. (2010) Integrated molecular genetic profiling of pediatric high-grade gliomas reveals key differences with the adult disease. J Clin Oncol 28: 3061–3068.

15. de Tayrac M, Etcheverry A, Aubry M, Saikali S, Hamlat A, et al. (2009) Integrative genome-wide analysis reveals a robust genomic glioblastoma signature associated with copy number driving changes in gene expression. Genes Chromosomes Cancer 48: 55–68.

16. Ferrer-Luna R, Mata M, Nunez L, Calvar J, Dasi F, et al. (2009) Loss of heterozygosity at 1p-19q induces a global change in oligodendroglial tumor gene expression. J Neurooncol 95: 343–354.

17. Pfister S, Remke M, Benner A, Mendrzyk F, Toedt G, et al. (2009) Outcome prediction in pediatric medulloblastoma based on DNA copy-number aberrations of chromosomes 6q and 17q and the MYC and MYCN loci. J Clin Oncol 27: 1627–1636.

18. Freire P, Vilela M, Deus H, Kim YW, Koul D, et al. (2008) Exploratory analysis of the copy number alterations in glioblastoma multiforme. PLoS One 3: e4076.

19. Ohgaki H, Dessen P, Jourde B, Horstmann S, Nishikawa T, et al. (2004) Genetic pathways to glioblastoma: a population-based study. Cancer Res 64: 6892–6899.

20. Ohgaki H, Kleihues P (2007) Genetic pathways to primary and secondary glioblastoma. Am J Pathol 170: 1445–1453.

21. Liu F, Park PJ, Lai W, Maher E, Chakravarti A, et al. (2006) A genome-wide screen reveals functional gene clusters in the cancer genome and identifies EphA2 as a mitogen in glioblastoma. Cancer Res 66: 10815–10823.

22. Ruano Y, Mollejo M, Ribalta T, Fiano C, Camacho FI, et al. (2006) Identification of novel candidate target genes in amplicons of Glioblastoma multiforme tumors detected by expression and CGH microarray profiling. Mol Cancer 5: 39.

23. Dreyfuss JM, Johnson MD, Park PJ (2009) Meta-analysis of glioblastoma multiforme versus anaplastic astrocytoma identifies robust gene markers. Mol Cancer 8: 71.

24. Jiang R, Mircean C, Shmulevich I, Cogdell D, Jia Y, et al. (2006) Pathway alterations during glioma progression revealed by reverse phase protein lysate arrays. Proteomics 6: 2964–2971.

25. Chakravarti A, Zhai G, Suzuki Y, Sarkesh S, Black PM, et al. (2004) The prognostic significance of phosphatidylinositol 3-kinase pathway activation in human gliomas. J Clin Oncol 22: 1926–1933.

26. Jeuken J, van den Broecke C, Gijsen S, Boots-Sprenger S, Wesseling P (2007) RAS/RAF pathway activation in gliomas: the result of copy number gains rather than activating mutations. Acta Neuropathol 114: 121–133.

27. Zohrabian VM, Forzani B, Chau Z, Murali R, Jhanwar-Uniyal M (2009) Rho/ROCK and MAPK signaling pathways are involved in glioblastoma cell migration and proliferation. Anticancer Res 29: 119–123.

28. Rao SK, Edwards J, Joshi AD, Siu IM, Riggins GJ (2010) A survey of glioblastoma genomic amplifications and deletions. J Neurooncol 96: 169–179.

29. Misra A, Pellarin M, Nigro J, Smirnov I, Moore D, et al. (2005) Array comparative genomic hybridization identifies genetic subgroups in grade 4 human astrocytoma. Clin Cancer Res 11: 2907–2918.

30. Czernicki T, Zegarska J, Paczek L, Cukrowska B, Grajkowska W, et al. (2007) Gene expression profile as a prognostic factor in high-grade gliomas. Int J Oncol 30: 55–64.

31. Brennan C, Momota H, Hambardzumyan D, Ozawa T, Tandon A, et al. (2009) Glioblastoma subclasses can be defined by activity among signal transduction pathways and associated genomic alterations. PLoS One 4: e7752.

32. Hirose Y, Aldape KD, Chang S, Lamborn K, Berger MS, et al. (2003) Grade II astrocytomas are subgrouped by chromosome aberrations. Cancer Genet Cytogenet 142: 1–7.

33. Marko NF, Prayson RA, Barnett GH, Weil RJ (2010) Integrated molecular analysis suggests a three-class model for low-grade gliomas: a proof-of-concept study. Genomics 95: 16–24.

34. Qaddoumi I, Sultan I, Broniscer A (2009) Pediatric low-grade gliomas and the need for new options for therapy: Why and how? Cancer Biol Ther 8: 4–10.

35. Liang BC, Hays L (1996) Mitochondrial DNA copy number changes in human gliomas. Cancer Lett 105: 167–173.

36. Arslantas A, Artan S, Oner U, Muslumanoglu MH, Ozdemir M, et al. (2007) Genomic alterations in low-grade, anaplastic astrocytomas and glioblastomas. Pathol Oncol Res 13: 39–46.

37. Idbaih A, Carvalho Silva R, Criniere E, Marie Y, Carpentier C, et al. (2008) Genomic changes in progression of low-grade gliomas. J Neurooncol 90: 133–140.

38. Koschny R, Koschny T, Froster UG, Krupp W, Zuber MA (2002) Comparative genomic hybridization in glioma: a meta-analysis of 509 cases. Cancer Genet Cytogenet 135: 147–159.

39. Kotliarov Y, Steed ME, Christopher N, Walling J, Su Q, et al. (2006) High-resolution global genomic survey of 178 gliomas reveals novel regions of copy number alteration and allelic imbalances. Cancer Res 66: 9428–9436.

40. Batchelor TT, Betensky RA, Esposito JM, Pham LD, Dorfman MV, et al. (2004) Age-dependent prognostic effects of genetic alterations in glioblastoma. Clin Cancer Res 10: 228–233.

41. Campbell BA, Horsman DE, Maguire J, Young S, Curman D, et al. (2008) Chromosomal alterations in oligodendroglial tumours over multiple surgeries: is tumour progression associated with change in 1p/19q status? J Neurooncol 89: 37–45.

42. Wong KK, Tsang YT, Chang YM, Su J, Di Francesco AM, et al. (2006) Genome-wide allelic imbalance analysis of pediatric gliomas by single nucleotide polymorphic allele array. Cancer Res 66: 11172–11178.

43. Kuga D, Mizoguchi M, Guan Y, Hata N, Yoshimoto K, et al. (2008) Prevalence of copy-number neutral LOH in glioblastomas revealed by genomewide analysis of laser-microdissected tissues. Neuro Oncol 10: 995–1003.

44. Vranova V, Necesalova E, Kuglik P, Cejpek P, Pesakova M, et al. (2007) Screening of genomic imbalances in glioblastoma multiforme using high-resolution comparative genomic hybridization. Oncol Rep 17: 457–464.

45. Peiffer DA, Le JM, Steemers FJ, Chang W, Jenniges T, et al. (2006) High-resolution genomic profiling of chromosomal aberrations using Infinium whole-genome genotyping. Genome Res 16: 1136–1148.

46. Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A (2012) GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res 40: W478–483.

47. Zhang B, Kirov S, Snoddy J (2005) WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 33: W741–748.

48. Amberger J, Bocchini CA, Scott AF, Hamosh A (2009) McKusick's Online Mendelian Inheritance in Man (OMIM). Nucleic Acids Res 37: D793–796.

49. Knudson AG Jr (1971) Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci U S A 68: 820–823.

50. Bredel M, Bredel C, Juric D, Harsh GR, Vogel H, et al. (2005) High-resolution genome-wide mapping of genetic alterations in human glial brain tumors. Cancer Res 65: 4088–4096.

51. Stupp R, Tonn JC, Brada M, Pentheroudakis G (2010) High-grade malignant glioma: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol 21 Suppl 5v190–193.

52. Cavaliere R, Lopes MB, Schiff D (2005) Low-grade gliomas: an update on pathology and therapy. Lancet Neurol 4: 760–770.

53. Pouratian N, Schiff D (2010) Management of low-grade glioma. Curr Neurol Neurosci Rep 10: 224–231.

54. Costantini M, Clay O, Auletta F, Bernardi G (2006) An isochore map of human chromosomes. Genome Res 16: 536–541.