GABA-A Channel Subunit Expression in Human Glioma Correlates with Tumor Histology and Clinical Outcome

GABA (γ-aminobutyric acid) is the main inhibitory neurotransmitter in the CNS and is present in high concentrations in presynaptic terminals of neuronal cells. More recently, GABA has been ascribed a more widespread role in the control of cell proliferation during development where low concentrations of extrasynaptic GABA induce a tonic activation of GABA receptors. The GABA-A receptor consists of a ligand-gated chloride channel, formed by five subunits that are selected from 19 different subunit isoforms. The functional and pharmacological properties of the GABA-A channels are dictated by their subunit composition. Here we used qRT-PCR to compare mRNA levels of all 19 GABA-A channel subunits in samples of human glioma (n = 29) and peri-tumoral tissue (n = 5). All subunits except the ρ1 and ρ3 subunit were consistently detected. Lowest mRNA levels were found in glioblastoma compared to gliomas of lower malignancy, except for the θ subunit. The expression and cellular distribution of the α1, γ1, ρ2 and θ subunit proteins was investigated by immunohistochemistry on tissue microarrays containing 87 gliomas grade II. We found a strong co-expression of ρ2 and θ subunits in both astrocytomas (r = 0.86, p<0.0001) and oligodendroglial tumors (r = 0.66, p<0.0001). Kaplan-Meier analysis and Cox proportional hazards modeling to estimate the impact of GABA-A channel subunit expression on survival identified the ρ2 subunit (p = 0.043) but not the θ subunit (p = 0.64) as an independent predictor of improved survival in astrocytomas, together with established prognostic factors. Our data give support for the presence of distinct GABA-A channel subtypes in gliomas and provide the first link between specific composition of the A-channel and patient survival.


Introduction
Gliomas are the most common form of primary brain tumor with an overall incidence of about 4-5 per 100.000 persons per year [1,2]. The majority of gliomas consist of glioblastomas, which are highly proliferative and invasive tumors characterized by remarkable biological heterogeneity and poor response to present treatments [3]. The prognosis for patients with diffuse low-grade gliomas is more favorable, but these tumors transform into malignant gliomas over time with eventually fatal outcome [4]. The etiology of gliomas is largely unknown, with exposure to highdose ionizing radiation as one of the few recognized risk factors [2].
GABA is the main inhibitory neurotransmitter in the central nervous system. The most common type of GABA receptors is the GABA-A channel, consisting of a ligand-gated pentameric chloride channel that is normally closed but can be opened by GABA [5]. Nineteen different GABA-A channel subunits have been cloned and are grouped into eight separate subfamilies (a1-6, b1-3, c1-3, d, e, h, p and r1-3). All nineteen subunits are expressed in the brain [6]. The GABA-A channel consists most often of three types of subunit isoforms: two as, two bs and a third type of subunit. All neurons contain GABA-A ion channels but the channel subtypes change during development, varying also between different brain regions and different types of neurons. The distinctive functional and pharmacological of the GABA-A channels are dictated by the composition of the subunits and can be modulated by intracellular proteins [7].
GABA is present in high concentrations in the presynaptic terminals of neuronal cells. In the postsynaptic terminal, brief exposure to high concentration of GABA results in opening of GABA-A channels and a subsequent increase in membrane conductance known as phasic inhibition. It has been shown that GABA-mediated signaling in the brain occurs also in a less spatially and temporally restricted manner, through low concentrations in the extracellular space that result in a persistent or tonic activation [5,8]. This tonic activation of GABA-A channels was first identified in voltage-clamp recordings in hippocampal and cerebellar neurons but is known to occur more generally in the mammalian brain [9]. Extrasynaptic pharmacologically active GABA-A channels in the mature brain have been found in nonneuronal cells such as astroglial cells [10,11], and also in T lymphocytes [12], suggesting an immunoregulatory role for extrasynaptic GABA [13].
There is growing evidence also for a widespread role of GABA in the growth regulation of many cell types, including neuronal stem cells and perhaps tumor stem cells [6]. In the postnatal subventricular zone along the lateral ventricle where adult neurogenesis occurs, GABAergic signaling is involved in the proliferative control of neuroblasts [6]. In primary brain tumors, GABA-A channel subunits have been detected in the neuronal component of gangliogliomas [14], but also in gliomas where the response to GABA correlated with the malignancy grade of the tumor [15,16]. GABA-evoked current response was restricted to low-grade gliomas and not recorded in glioblastoma, suggesting that the disappearance of GABA-A channels parallels the unlimited growth of malignant gliomas [16]. When co-cultured with neurons, an interaction between neurons and glioma cells was triggered resulting in functional expression of GABA-A channels by tumor cells within 24 hours [17].
These reports are all in favor of a role for extrasynaptic GABA in the growth control of gliomas, but so far this has not been substantiated by studies on clinical patient material. Also, a systematical analysis of the distribution of the different GABA-A channel subunit isoforms in human gliomas has not been performed. In the present work we compared mRNA levels of all 19 GABA-A channel subunits in gliomas of different malignancy grade and histological subtype. We then studied the distribution of the a1, c1, r2 and h subunit proteins in vivo in a clinical cohort of diffuse low-grade gliomas and correlated the presence of GABA-A channel subunit proteins with patient survival.

Quantification of GABA-A channel subunit mRNAs levels
The histological classification of the tumors samples used for quantitative real-time PCR (qRT-PCR) is shown in Table 1. Mean expression levels of GABA-A channel subunit mRNAs were compared group-wise between gliomas grade II (n = 12), gliomas grade III (n = 10) and glioblastomas (n = 7) according to the World Health Organization (WHO) classification of brain tumors [1]. Seventeen GABA-A channel subunit mRNAs were consistently detected, whereas the r1 and r3 subunits were not expressed. The relative mRNA levels of the 17 detectable GABA-A subunits in gliomas of different malignancy grade are shown in Figure 1. Statistical analysis of quantitative differences between gliomas grade II, gliomas grade III and glioblastomas showed significantly higher mRNA levels of a1, a6, c1 and c2 GABA-A channel subunits in gliomas grade II compared to glioblastomas, and of a3, a6, b3, c1, c2 and p subunits in gliomas grade III compared to glioblastomas (p,0.05). In contrast, mRNA levels of the h subunit were 5-10 fold higher in glioblastomas than in gliomas with lower malignancy grade (p,0.05) ( Figure 2).
For 15 of the 17 detectable GABA-A subunits, no statistically significant differences in mRNA levels were found between astrocytomas and oligodendrogliomas grade II ( Figure S1) or between astrocytomas and oligodendrogliomas grade III ( Figure  S2). The two exceptions were the b1 subunit, for which expression was lower in astrocytomas grade II than in oligodendrogliomas grade II ( Figure S1), and the a3 subunit, showing lower expression in astrocytomas grade III than in oligodendrogliomas grade III ( Figure S2) (Mann-Whitney Rank Sum Test).
In peri-tumoral brain tissue, mRNAs for the 17 GABA-A channel subunits were consistently detected ( Figure S3). For most subunits belonging to the a, b and c subfamilies, as well as for the d subunit, mRNA levels were 3-5 fold higher in peri-tumoral brain tissue compared to tumor samples. Other subunits were expressed at approximately similar levels as those found in tumors ( Figure S3). Table 2 summarizes patient characteristics and tumor characteristics of the 91 gliomas that were used for immunohistochemical analysis. The results of the immunostainings are shown in Table 3. Tumor cells with immunoreactivity for a1, c1, r2 and h GABA-A channel subunits were identified in all histological tumor types. The subcellular distribution of the subunit proteins was found to be cytoplasmic and membranous. Of the total number of 91 samples, 17 showed immunoreactivity for the a1 subunit (19%), 24 for the c1 subunit (26%), 29 for the r2 subunit (32%), and 35 for the h subunit (38%) (Figure 3). The mean fraction of immunoreactive tumor cells denoted as + was between 4-8% for the a1, c1, r2 and h subunits ( Table 3). The mean fraction of immunoreactive tumor cells in astrocytomas denoted as ++ was 15% for the a1 subunit (n = 1), 20% for the c1 subunit (n = 2), 55% for the r2 subunit (n = 5) and 40% for the h subunit (n = 9). In oligodendrogliomas, the mean fraction of immunoreactive tumor cells denoted as ++ was 30% for the r2 subunit (n = 1) and 50% for the h subunit (n = 2), while in oligoastrocytomas these percentages were 50% for the c1 subunit (n = 3), 35% for the r2 subunit (n = 2), and 40% for the h subunit (n = 3).

Distribution of GABA-A channel subunit proteins in gliomas
Three astrocytomas grade II showed immunoreactivity for all four subunits, while none of the oligodendrogliomas and oligoastrocytomas expressed all four subunits. Interestingly, more than 80% of the gemistocytic astrocytomas expressed the r2 and h subunits, but lacked immunoreactivity for the a1 and c1 subunit proteins ( Figure 4). A similar pattern of predominantly r2 and h subunit expression was found in minigemistocytic oligoastrocytomas (Table 3).

Statistical analysis
As a first step, we dichotomized the immunohistochemistry results of the four GABA-A channel subunits into two categories; no immunoreactivity ( = 0) versus positive immunoreactivity (+ or ++). These data are illustrated in Figure 5, where each number in the figure represents one tumor sample consisting of astrocytomas (n = 42), oligodendrogliomas (n = 27) or oligoastrocytomas (n = 18) grade II, with either no (no bar) or positive (colored bar) immunoreactivity for the four subunits. As demonstrated, there was a frequent co-expression of r2 and h subunits in astrocytomas grade II (r = 0.86, p,0.0001) and in oligodendrodendrogliomas/ oligoastrocytomas grade II (r = 0.66, p,0.0001). No statistically significant correlation between the a1 and c1 subunit or between any of the other four subunits was found.

Survival analysis
The mean postoperative survival for patients with gliomas grade II (n = 87) was 7,9 y. As expected, mean survival was longer amongst patients with oligodendrocytomas compared to oligoastrocytomas and astrocytomas (11.6 y, 6.9 y and 5.9 y respectively). We hypothesized that r2 expression and h subunit expression in grade II gliomas were prognostic factors for survival and tested these two parameter as dichotomized variables in the survival analysis for astrocytomas grade II (n = 42) and for the combined group of oligodendrogliomas and oligoastrocytomas (n = 45).
Both r2 subunit expression and h subunit expression were associated with statistically significant longer median survival in the histological subgroup of astrocytomas (Table 4). Of the established clinical and radiological prognostic factors for grade gliomas grade II (age, epilepsy, KPS, tumor size, tumor location, contrast enhancement) [18], the parameters seizures at onset, KPS$90 and tumor diameter #6 cm were identified as prognostic factors in the univariate analysis (Table 4). Young age (#40 years), non-enhancing tumor and non-central tumor localization were not associated with favorable prognosis. The presence of mutated IDH1 R132H protein showed a trend towards longer survival in the univariate analysis and was entered as a variable in the multivariate model [19]. Since all patients but one were treated by radiotherapy, this parameter was not tested. Likewise, the LOH 1p/19q status was positive for only two patients and could not be entered into the model. As shown in Table 4, expression of r2 subunits contributed to significantly longer survival in the multivariate model, together with seizures at onset, KPS$90, and mutated IDH1 protein. Expression of h subunits was not identified as an independent prognostic factor in the multivariate model. Small tumor size was a significantly favorable prognostic factor entered as a single variable (p = 0.019), but lost its significance in the multivariate model.
In the group of oligodendrogliomas and oligoastrocytomas (n = 45), the parameters oligoastrocytoma histology (versus oligodendroglioma) and tumors crossing midline structures (versus noncentral tumor location) were identified as unfavorable prognostic factors in both the univariate and multivariate model ( Table 4). The parameters expression of r2 subunit and expression of h subunit were not associated with significantly longer survival in these tumors.

Discussion
The aims of the present study were 1) to explore the distribution of GABA-A channel subunits in human gliomas of various histological subtypes, and 2) to search for correlations between GABA-A channel subunits and patient survival. We showed that 17 of the 19 different GABA-A channel subunit mRNAs were  present, only r1 and r3 subunits could not be detected. While no major differences in subunit mRNA levels were found between astrocytomas and oligodendrogliomas, there was a general downregulation of GABA-A channel subunits in glioblastoma. Thus, highest mRNA levels were recorded in gliomas grade II and grade III, except for the h subunit that showed 5-10 fold higher levels in glioblastomas.
Our findings are in agreement with previous reports demonstrating functional bindings sites for the GABA-A channel receptor in low-grade and anaplastic gliomas but not in glioblastomas [16,17]. The lack of responsiveness to GABA in glioblastomas has been explained by either the loss of GABA-A channels or by defective channels [16]. Thus, tumor cells equipped with functional GABA-A channels may be able to maintain low or intermediate proliferative activity upon physiological levels of GABA in the extracellular fluid. GABA triggered a depolarization in the majority of glioma cells and a concomitant increase of intracellular Ca 2+ [16]. Since the GABA-A channel itself is not permeable for Ca 2+ ions, the Ca 2+ influx has been associated with an activation of voltage-gated Ca 2+ channels in the tumor cell membrane. The central role of intracellular Ca 2+ as a downstream effector of GABA receptor signaling in cell proliferation has been confirmed [20].
The other main finding of this study is the strong co-expression of r2 and h subunit proteins in gliomas grade II and the association between r2 subunit expression and longer patient survival in diffuse astrocytomas. In spite of its strong co-expression with the r2 subunit, the h subunit did not emerge as an independent predictor of survival in the multivariate model. These findings are in agreement with the qRT-PCR results, showing an opposite trend for the h subunit with up-regulation in glioblastomas. The prognostic impact of r2 subunits in diffuse astrocytomas is intriguing but is not necessarily unique for the GABA A-channel subunit families. It is possible that we failed to identify associations for other subunits in oligodendrogliomas and oligoastrocytomas, which were combined into one group in the present study due to the limited number of tumor samples. Neither can we exclude that there has been a selection bias influencing our results, since many tumor biopsies from the original clinical cohort were too small to be included in the TMA. Thus, our study provides a first link between GABA-A channel composition and survival in gliomas and supports a presumed functional role of GABA in gliomas, but the findings need to be confirmed in larger studies.
The GABA-A channel subunits exist as a family of subtypes with distinct temporal and spatial patterns of expression and distinct properties that presumably underlie a precise role for each subtype [21,9]. The differential sensitivity to GABA is dependent on the specific subunit composition of the A-channel, and may further be influenced by intracellular proteins interacting with the channel complex and by post-translational modifications [22]. Whilst the a1 subunit has frequently been the focus of studies, due to its abundance in the brain and its reliability of measurement, considerably less is known on the other subunits expressed at a lower level such as h and the r2 subunits. In the brain, functional GABA-A channel receptors containing the h subunit are formed together with a, b and c subunits [7]. Synaptic channels are activated by high GABA concentrations, whereas extrasynaptic channels can be activated by extremely low concentrations of GABA. Tumors such as gliomas are probably not exposed to high GABA concentrations and one would expect mostly tonic currents to be active in these tumors. It seems that tonic currents more or less can contain any subunit in their channel complex, whereas synaptic channels are more restricted to certain types of subunits. As such, the c2 subunit is a major component of synaptic channels, while the a4 a5 a6 and d subunits are part of extrasynaptic GABA A-channels and thus involved in tonic currents [7,23]. The a1, a2 and a3 subunits can be located both at and outside of synapses and the e and h subunits are probably mostly extrasynaptic [7,23].
In the present study, r2 subunit expression was found in a relatively large number of glioma samples while no detectable mRNA levels were found by qRT-PCR for the r1 and r3 subunits. The r1 subunit is predominantly expressed in the retina and visual pathways, while r2 expression has also been found in hippocampus and amygdala [24,25]. Rho subunits form functional homomeric GABA-A-r receptors (previously known as GABA C channels) [7,26], but can potentially also participate in heteromeric GABA-A channel formation [7,27]. The GABA-A-r receptor has a number of distinctive and unique functional properties, such as a long mean opening time of the channel and slow desensitization, suggesting a more widespread function than previously thought [28]. In addition, the pharmacological properties of the GABA-A-r receptor, including the lack of response to benzodiazepines and barbiturates, set apart this class of receptor.  We found a strong positive correlation between the expression of h and r2 subunits. Expression of the h and r2 subunits in brain tissue is relatively uncommon, and there is no evidence to suggest that they form a functional channel together [29,30]. However, normal transcriptional regulation in tumor cells is disrupted and may not follow the normal ''rules'' for GABA-A channel subunit expression. It is therefore possible that they could have been expressed without participating in forming functional channels. Interestingly, the strongest correlation between the expression of r2 and h subunits was found in gemistocytic astrocytomas. The presence of gemistocytes in a tumor has been correlated with worse prognosis [31,32], although the proportion of gemistocytes in the tumor does not seem to have an effect [32]. The reasons for development of gemistocytes in tumors are still poorly understood, but it has been suggested that they may develop due to losing the competition for substrates with adjacent cells [31,32].
It is not known whether tumor cells are the source of GABA. Only neurons contain glutamic acid decarboxylase (GAD), the enzyme that catalyzes the decarboxylation of glutamate to GABA, and are able to produce GABA. However, astrocytes can take up GABA and there is a minor pathway not involving GAD that allows synthesis of GABA [33]. Both neurons and astrocytes can release GABA, but only neurons can release high concentrations of GABA within a short time-span from synaptic vesicles at the presynaptic terminal [33].

Conclusions
The present study shows a down-regulation of subunit mRNA levels in glioblastomas, except for the h subunit, and the presence of distinct GABA-A channel subunit proteins in gliomas grade II. The correlation between r2 subunits and favorable survival in diffuse astrocytomas suggests that specific subunit compositions affect clinical outcome in glioma. Taken the distinct pharmacological properties of the different GABA-A channel subtypes, this may open up for new therapeutic strategies for glioma. Our study also highlights the complexity of GABA signaling in gliomas and stresses the need for functional studies in this field.

Ethics Statement
Informed consent was obtained for the use of human brain tissue and for access to medical records for research purposes, and all material was obtained in a manner compliant with the Declaration of Helsinki. The studies involving human tissue samples were approved by the Ethics Committee of Uppsala University (Application Dnr Ups 02-330) and the Ethics Committee of Karolinska Institutet (Application Dnr Ki 02-254). Written informed consent was obtained prior to sample collection. Tissue Samples for quantitative real-time quantitative PCR Human tissue samples for quantitative real-time quantitative PCR (qRT-PCR) were obtained from the files of the department of neuropathology of the Academic Medical Center, University of Amsterdam. Histological classification of tumors according to the WHO is shown in Table 1. Twenty-nine gliomas (12 gliomas grade II, 10 gliomas grade III, 7 primary glioblastomas), all with temporal lobe location, were included. Five samples of peritumoral brain derived from patients with primary glioblastomas in the temporal lobe were included. These samples consisted of tumor-free cortex on the basis of absence of immunoreactivity for Ki67 (1:500; Rabbit polyclonal, DAKO), IDH1 R132H (1:100; Clone H09, Dianova) and p53 (1:2000; Clone DO-7+BP53-12, Neomarkers). Samples were quickly removed at surgery and immediately divided into two parts; one part was fixed in 4% paraformaldehyde for 24 hours, paraffin embedded and used for histopathological diagnosis, the remaining part was snap frozen in liquid nitrogen and maintained at 280uC until used for RNA isolation.

Quantitative real-time PCR
Total RNA was isolated with TRIzol reagent according to the manufacture's instruction and further transcribed to cDNA using oligodT primers (Table S1). qRT-PCR was performed in a 10 ml reaction mixture containing 4 ml cDNA (0.5 ng/ml), 16 PCR    (Table S1). Amplification was performed in 96-well optical plates using the ABI PRISM 7900HT Sequence Detection System (Applied Biosystems) with an initial denaturation of 5 min at 95uC, followed by 45 cycles of 95uC for 15 s, 60uC for 30 s and 72uC for 30 s. A melting curve was determined at the end of cycling to ensure the amplification of a single PCR product. Each reaction was run in duplicate. Cycle threshold values (Ct) were determined with the SDS 2.3 software supplied with the instrument. The expression of each target gene relative to a reference gene TATA-binding protein (TBP) was calculated using the 2 2DCt method as previously described [34].

Tissue samples for immunohistochemistry
The distribution of GABA-A channel subunit proteins was explored by immunohistochemistry in a tissue microarray (TMA) of glioma samples. The TMA was prepared from paraffinembedded tissue blocks as previously described [35]. Samples were included from patients aged $16 years old operated for low-grade gliomas between 1984-2001 at the Uppsala University Hospital, as part of a previously described clinical cohort of diffuse lowgrade gliomas [36]. Paraffin-embedded tissue blocks were identified and re-used for preparation of the TMA. All tumors included in the TMA were re-evaluated for the purpose of the study by one of the neuropathologists (IA). Of the original cohort of 116 grade II gliomas, 3 tumors were re-classified as non-gliomas, 18 blocks contained too little tumor for TMA preparation (mainly biopsies), and in 4 blocks areas with infiltrative tumor edge were left without representative areas of tumor bulk. Thus, 91 glioma samples remained and were included in the TMA presented in this study. Clinical and histopathological tumor characteristics are presented in Table 3. Clinical data collection A retrospective chart review was performed for all tumor cases as previously described [36]. The following clinical data were  Tissue microarray preparation TMAs were prepared by the Swedish Proteome Resource Centre (HPR) facilities at the department of pathology, the Rudbeck laboratory, Uppsala University Hospital [37]. In brief, donor blocks were sectioned, hematoxylin-eosin stained and immunostained with immunohistochemical markers to select appropriate areas for the TMA construct. Doublets of representative areas were included for each tumor. Using an automated tissue arrayer (Beecher Instruments, Silver Spring, MD), two 1.0mm diameter punches were taken from each donor block and transferred to the recipient TMA block. In total, two TMA blocks were prepared.

Fluorescent in situ hybridization
All samples included in the TMA were studied by fluorescent in situ hybridization analysis (FISH) to identify losses of chromosomal arms 1p and 19q (LOH 1p/19q), as part of our previous study [36]. The commercially purchased probes used for hybridization were Zytolight SPEC 1p36/1q25 and 19q13/19p13 dual color probes (Nordic BioSite, Sweden). Slides were assessed under a fluorescence microscope (Olympus BX 50 Deutschland GmbH), and a minimum of 100 non-overlapping nuclei was calculated for each hybridization. A tumor was considered deleted if .50% of the nuclei harbored two signals of the reference probe but only one signal of the target probe.

Evaluation of immunohistochemistry
The neuropathologists in this study (EA, IA) evaluated all imunostainings and identified immunoreactive cell types as well as the cellular distribution of the proteins. For each sample, the entire piece of micro tissue was examined through light microscopy at magnification 20-406. The percentage of entrapped neurons, identified by immunoreactivity for neuronal nuclear protein, was estimated for all samples in the TMA. Representative areas containing the highest density of immunoreactive cells for the GABA-A channel subunits were used for counting $200 tumor cells per section. The fraction of immunoreactive tumor cells labeled with GABA-A subunit antibodies was then calculated and graded as 0-2 (0 = no positive cells; 1 (+) = few (,10%) positive cells; 2 (++) = several ($10%) positive cells).

Characterization of the TMA
Clinical characteristics and histological and molecular characterization of the tumors included in the TMA are shown in Table 3. Of the 91 tumors included in the present study, 4 tumors were re-classified as high-grade gliomas (3 oligoastrocytomas grade III, 1 astrocytoma grade III). The remaining tumors consisted of astrocytoma grad II (n = 42), oligodendroglioma grade II (n = 27) or oligoastrocytoma grade II (n = 18). Histological analysis revealed a fraction of 0-5% of all cells to consist of entrapped neurons (mean percentage in all sections 2.7%). Molecular characterization confirmed the presence of mutated IDH1 R132H protein in 76/91 tumor samples (84%) and LOH 1p/ 19q in 26/91 tumor samples (29%) ( Table 3).

Statistical analysis
Statistical analysis comparing mRNA levels of individual GABA-A channel subunits between different types of gliomas was carried out using SigmaPlot v11 (Systat Software Inc., USA), and assessed by Kruskal-Wallis ANOVA on ranks, with significance level set to p,0.05. The spearman correlation test was used to calculate the coefficients between a1, c1, r2 and h subunit protein expression and was performed in JMP statistical software, version 5.0.1a (SAS Institute Inc., Cary, North Carolina, USA).
Survival curves were plotted according to the Kaplan-Meier method (product-limit method) and the log-rank probability test (Mantel-Cox) estimated the prognostic value of each specific clinical, histological, molecular or radiological parameter in the univariate analysis. The Cox proportional hazard model was used to calculate the impact of each prognostic factor in the multivariate analysis. Stepwise exclusion of variables was used to achieve a model with as few variables as possible. The level for confounders to be removed from the model when adjusted for variables already in the model (P-to-remove) was set to .0.1. The natural logarithm (ln) of the cumulative hazard plots was used to confirm the assumption of the proportional hazard functions. Statistical calculations were performed in JMP, version 5.0.1a (SAS Institute Inc., Cary, North Carolina, USA). Figure S1 Detailed presentation of the RT-PCR results for each of the 17 subunits, showing quantitative mRNA levels between astrocytomas (n = 6) and oligodendrogliomas grade II (n = 6). The normalized mRNA expression of each target gene relative to a reference gene TATA-binding protein (TBP) was calculated using the 2 2DCt method. (TIF) Figure S2 Detailed presentation of the RT-PCR results for each of the 17 subunits, showing quantitative mRNA levels between astrocytomas (n = 5) and oligodendrogliomas grade III (n = 5). The normalized mRNA expression of each target gene relative to a reference gene TATA-binding protein (TBP) was calculated using the 2 2DCt method. (TIF) Figure S3 Relative mRNA levels for 17 GABAR A subunits in respectively peri-tumoral tissue (n = 5), gliomas grade II (n = 12), gliomas grade III (n = 10), and glioblastomas (GBM) (n = 7). (TIF)