Research article | Open | Published:
Glutathione S-transferase omega genes in Alzheimer and Parkinson disease risk, age-at-diagnosis and brain gene expression: an association study with mechanistic implications
Molecular Neurodegenerationvolume 7, Article number: 13 (2012)
Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects).
We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10-11-1.9 × 10-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003).
These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.
Glutathione S-Transferase (GST) family of genes have been implicated in multiple neuropsychiatric [1–4] and neurodegenerative diseases [5–11]; where altered levels or function of these enzymes is thought to impact levels of oxidative stress and/or inflammation in a way that contributes to disease susceptibility. A linkage locus on chromosome 10q that has been implicated in both Alzheimer's (AD)[11–13] and Parkinson's disease (PD) harbors two GST genes of the omega class: GSTO1 and GSTO2, which are approximately 75 kb apart.
GSTOs have enzymatic activities as thioltransferases and dehydroascorbate reductases that promote antioxidant activity and can also function in metabolism of drugs and toxins. Additionally, GSTO1 was shown to promote activation of the pro-inflammatory cytokine, interleukin-1β (IL-1β) by post-translational processing. Given their location and function, they have been studied as candidate genes in AD and PD[5, 6, 9, 11, 14, 16–18]. Li et al. compared hippocampal gene expression levels in 6 AD vs. 2 control brains and identified significantly lower GSTO1 levels in the AD hippocampi. This group studied AD and PD families that showed linkage to chromosome 10, using the age-at-onset phenotype  and identified association of multiple SNPs at the GSTO locus with delayed age-at-onset of both diseases, with the strongest effects observed for GSTO1 rs4925 and GSTO2 rs2297235 SNPs that are in tight linkage disequilibrium (LD). No significant influence was detected for either AD or PD risk in this study.
Since this initial report, several follow-up studies have been published with mixed approaches and results. Kölsch et al. reported association of rs4925 with earlier age-at-onset of AD, thus in opposite direction to the original report, and no effect on AD risk. Lee et al. found modest association of rs4925 with AD risk in Carribean-Hispanic families that show linkage to chromosome 10q as did Capurso et al. in an Italian case-control series, though neither study detected an age-at-onset effect. A case-control study by Wahner et al. was the only report for an effect of GSTO locus on PD risk, with both rs4925 and rs2297235 conferring protection, especially in those with smoking history. Additionally, several studies reported lack of association with age-at-onset or risk of AD[17, 18] or PD.
Additional investigation of the GSTO locus is needed to further elucidate the role of these genetic variants in AD and PD, especially given the potential to establish the glutathione metabolism as a molecular pathway that is common to multiple, chronic neurologic diseases. An important shortcoming of most prior reports on the GSTO locus is the modest sample sizes, which could underlie the inconsistent results likely due to lack of power, sample or locus heterogeneity or a combination of these factors. Both AD[19, 20] and PD are complex diseases with substantial genetic component. Some of the genetic risk for these diseases has been identified via linkage and association studies and shown to influence age-at-onset[19–23]. More recently, genome-wide association studies (GWAS) of AD[24–28] and PD[29, 30], with sample sizes exceeding 10,000 subjects provide considerably greater power for detection of susceptibility loci. Despite their advantages, GWAS do not explain all of the underlying genetic component of these and other complex diseases, thus necessitating alternative approaches, including analysis of quantitative phenotypes.
In this study, we assessed the GSTO locus for its role in AD and PD, using an in-depth approach aimed at surmounting these challenges. Given the original report of association with delayed age-at-onset of AD and PD [5, 7], and with risk of AD in some follow-up studies [9, 11], we postulate that GSTO locus variants confer risk of LOAD in older age. We have a collection of > 8,000 late-onset AD (LOAD) case-controls, which includes a large series of older subjects ≥ 80 years of age-at-diagnosis/death (clinical/autopsied LOADs) or evaluation (controls). We analyzed two previously reported, coding SNPs in GSTO1 and GSTO2 for association both with disease risk and age-at-diagnosis in the LOAD series, as well as a large PD series. Reduced expression levels of GSTO1  and other glutathione metabolism genes [10, 32] have been reported in AD. We therefore analyzed the GSTO SNPs for association with brain GSTO1 and GSTO2 levels in > 750 brain samples from autopsied subjects with AD and other brain pathologies to determine whether they influenced disease risk by affecting brain gene expression. In an expression GWAS (eGWAS) testing association of 24,526 transcript levels measured in these brain samples with 213,528 cisSNPs within ± 100 kb of the tested transcript, we identified 686 genes that have significant cisSNP/transcript associations (in-press, PLoS Genetics). We analyzed these genes to discover molecular pathways that are enriched for genes with significant brain cisSNPs, and identified glutathione metabolism as one of the top pathways. Our results suggest that GSTO locus variants influence brain GSTO2 levels and confer AD risk at older age. These findings have mechanistic implications for the GSTO locus and glutathione metabolism genes, which should be explored further in AD and other chronic, neurologic diseases to identify functional variants that influence disease risk by altering brain gene expression levels.
Association of GSTOlocus SNPs with LOAD and PD risk
GSTO locus SNPs rs156697 and rs4925 were tested for association with disease risk in an older series (> 80 years) of 1,368 LOADs vs. 1,623 controls; in a younger series of 2,193 LOADs vs. 3,060 controls (60-80 years) and in 678 PDs vs. 712 controls (Table 1 Table 2 Table 3 and Table 4), using logistic regression analysis. GSTO2 rs156697 was significantly associated with LOAD in the older series (p = 0.038), with the minor allele conferring increased risk (OR = 1.14, 95% confidence interval = 95%CI = 1.01-1.30) (Table 5). There was no other significant disease risk association in the combined younger LOAD series (Table 6) or LOAD series of all ages (Table 7), although there were trends for increased LOAD risk with the minor allele of rs156697 in the LOAD series of all ages (p = 0.18, OR = 1.06) and of rs4925 in the older LOAD series (p = 0.15, OR = 1.10). Analysis of the six individual LOAD series of the older age group revealed consistently increased risk estimates for rs156697 (Table 5), which is also evident from the meta-analysis (Figure 1a), where there is no evidence of between-series heterogeneity (Breslow-Day p = 0.97). In contrast, the younger LOAD series is significantly heterogeneous (Breslow-Day p = 0.004), with three series showing increased (JS, NW, PS) and the other three (RS, AUT, NCRAD) with protective risk estimates (Table 6 Figure 1b).
PD series were composed of those with (PD FAM) or without (PD-SPO) family history of PD (Table 4). Assessment of these individual series revealed significant association of rs4925 with lower PD risk in the familial PD series (p = 0.034, OR = 0.78) and a trend for decreased risk with rs156697 (p = 0.116, OR = 0.83) in this series (Table 8). There was no association with PD risk in the sporadic PD or the combined series.
Association of GSTOlocus SNPs with AD and PD age-at-diagnosis
We employed age-at-diagnosis (clinical) or death (autopsy) as the surrogate quantitative variable for age-at-onset in our LOAD subjects (3,561 LOADs), who had an age range of 61-105 (mean age = 78, Table 1 Table 2 and Table 3). Analysis of rs156697 and rs4925 did not identify any significant association with age-at-diagnosis in the combined or individual LOAD series, although there was a trend for rs156697 for delayed age-at-diagnosis in the combined series (p = 0.098) (Additional File 1: Table S1). There was no significant age-at-diagnosis association when the older and younger LOADs were assessed separately (data not shown). The combined PD series had an age range of 25-94 (mean age = 64, Table 4). There was no significant association of either GSTO SNP with PD age-at-diagnosis (Additional File 1: Table S1).
Association of GSTOlocus SNPs with brain gene expression
In an eGWAS using 773 brain samples, we measured levels of 24,526 transcripts from the cerebellum and temporal cortex of autopsied subjects with and without AD pathology (in-press, PLoS Genetics). Control subjects without AD pathology often had other brain pathologies. We determined the association of rs156697 and rs4925 with brain GSTO1 and GSTO2 levels in these series, using linear regression analysis. Both SNPs had significant association with brain GSTO2 (Table 9, 10), but not with GSTO1 levels (data not shown). This association was significant for GSTO2 transcript levels measured from both the cerebellum and the temporal cortex; and in both the AD and control autopsy subjects, although the effect size estimates appeared to be bigger for the ADs and the cerebellum. The minor alleles of both SNPs were associated with lower brain GSTO2 levels in all analyses, with an additive pattern of association (Figure 2).
Twenty other cisSNPs at the GSTO locus were tested for association with brain GSTO1 and GSTO2 levels in our eGWAS. Although rs156697 had the strongest association, many of the additional cisSNPs also showed significant association with brain GSTO2 (Additional File 1: Table S2, Additional File 2: Figure S1), but not GSTO1 (data not shown) levels. The strongest cisSNPs were in an LD block encompassing GSTO2 (Additional File 2: Figure S1.)
Discovery of glutathione metabolism pathway in a brain gene expression GWAS
In our brain eGWAS we identified 686 genes with cerebellar transcript levels that are significantly influenced by cisSNPs, which were submitted to pathway analysis to discover molecular pathways that are significantly regulated in the brain. Glutathione metabolism was identified as one of the most significant pathways (p = 0.0035), where six genes from this pathway out of the thirty that existed within our eGWAS had significant cisSNPs that influence both the cerebellar and temporal cortex levels of these genes (Additional Files 3: Figure S2a and 4: Figure S2b). Five of the significant genes are enzymes that directly catalyze the binding of reduced glutathione to substrates (GSTO2, GSTT1, GSTT2, GSTM3, GSTM5) and GCTG is involved in amino acid metabolism, including glutamate (Additional File 1).
GSTO1 and GSTO2, which are evolutionarily conserved genes, previously implicated in AD[5, 6, 9, 11] and PD[5, 16], have diverse attributed functions including antioxidant activity via generation of ascorbate (Vitamin C) [14, 34, 35]; biotransformation of inorganic arsenic[14, 34]; modulation of ryanodine receptors and thus calcium release and apoptosis; and post-translational processing of the pro-inflammatory cytokine, IL-1β. Given their functions which are relevant for the pathophysiology of neurodegenerative diseases and their location in linkage regions for AD[11–13] and PD, GSTO locus variants have previously been studied for their association with risk and age-at-onset of AD and PD with mixed results[5, 6, 9, 11, 14, 16–18].
In this study, we assessed two coding polymorphisms, rs4925 (Ala140Asp) in GSTO1 and rs156697 (Asn142Asp) in GSTO2 in a large LOAD series of > 8,000 subjects, ~3,000 of whom were older (> 80 years) and in a PD series of > 1,300 subjects including both familial and sporadic cases. We found significant LOAD risk association for the minor allele of rs156697 in older subjects and a suggestive trend for delayed age-at-diagnosis. These results are consistent with the original and some of the follow-up reports on this locus[9, 11], and suggest that the reported delay in age-at-onset is likely to be due to an increased risk conferred in older subjects. Given the age-dependent decline in key glutathione metabolism components and their role in mitigating oxidative stress, the postulate that risky GSTO variants lead to increased risk in older LOADs due to accumulation of oxidative damage with increasing age, is a plausible scenario. It should be emphasized that our study utilized age-at-diagnosis as a surrogate for age-at-onset and unrelated case-controls, rather than family-based series. These differences could underlie the marginal age association in our study, in comparison to the original study.
Given the tight LD (r2 = 0.73, D' = 0.94 in HapMap3) between the two coding SNPs tested for AD and PD risk association in this study, we did not correct for multiple testing. If corrected, the AD association in the older ADs would no longer be significant (p = 0.076). Furthermore, GSTO locus variants were not reported to have significant or suggestive association with AD risk in the recent, large GWAS[26–28]. Although, these findings could collectively suggest that the AD risk association in our study is a false positive, there are alternative explanations: First, the effect conferred by GSTO2 rs156697 is age-specific based on our results, and others. Additionally, unlike the older series in our study, the younger LOAD series had significantly heterogeneous results for the rs156697 SNP. Thus, the large LOAD GWAS need to be re-analyzed focusing on the different age groups and also for age-at-onset or diagnosis association. Second, the effect of the GSTO2 variant is likely modest for LOAD risk, despite strong effects on brain gene expression. Third, although GSTO2 rs156697 has the strongest effect on brain expression of this gene in that locus, it may still not be the functional variant, thus leading to weak or heterogeneous effects on LOAD risk. Our results in LOAD risk and brain gene expression provide support for functional variant discovery efforts in the GSTO2 region and screening of such variants for their effects in transcriptional assays.
There was no significant association of GSTO SNPs with disease risk in the combined PD series. This may not be surprising given the difference in sample size and therefore power between the LOAD and PD series. Whereas our older LOAD series (1,338 LOAD vs. 1,604 controls) have ~61% power to detect the effect of the GSTO2 rs156697 SNP (OR = 1.14), the combined PD series (661 PDs vs. 702 controls) tested for this SNP, has ~32% power to detect this effect at α = 0.05. There was, however, association with decreased risk in the familial PD cases for the GSTO1 rs4925 minor allele. Although consistent with one other study in PD, this finding requires further replication. It is intriguing to note that this variant also conferred a protective effect in the LEAPS-PD GWAS, which assessed PD sib-pairs in its first stage[38, 39]. The opposite direction of association in the familial PDs (and some of the younger LOAD series) vs. the older LOAD series could have several explanations including the tested SNPs not being functional themselves but marking different functional variants of opposing effects; heterogeneity due to different gene-gene or gene-environment interactions in different groups; and false positivity in some of the tested series.
Although both GSTO SNPs are in coding regions, they do not lead to any change in the enzymatic activities of GSTO1[34, 35] or GSTO2. While their effects on LOAD and PD could be due to other, untested alterations in protein function, another potential mechanism of action is influencing levels of gene expression. Indeed, both SNPs had highly significant effects on brain gene expression levels of GSTO2, but not GSTO1. Amongst the 22 cisSNPs tested for association with GSTO levels in our brain eGWAS, rs156697 had the strongest effect, where the risky minor allele was associated with lower brain GSTO2 levels. These results strongly suggest that the risk conferred by the GSTO locus is most likely due to variants which influence GSTO2 levels in the brain. These findings are biologically compatible with the very high antioxidant function of GSTO2, where its dehydroascorbate reductase activity was found to be 70-100% greater than that of GSTO1.
Brain expression levels of other key enzymes of glutathione metabolism are also significantly influenced by genetic variants, as was identified from pathway analysis of our significant brain eGWAS results. Given our findings with GSTO2 and other studies implicating glutathione metabolism genes in neurodegenerative diseases[8, 10, 32], it will be important to analyze these additional glutathione metabolism genes with high brain regulation, for variants that influence risk of AD and other neurodegenerative diseases.
In summary, our results support GSTO2 as a risk gene for older LOAD subjects, where risky genotypes reduce brain levels of this gene, which likely leads to accumulation of oxidative damage worsening with increasing age. These findings have implications for disease mechanism, as well as the search for genetic risk variants in AD and other neurodegenerative diseases. First, it will be important to analyze the existing large LOAD and PD risk GWAS by different age-strata and also using age-at-onset as the outcome, where available. Second, GSTO2 should be sequenced for variants that may influence gene expression and thereby disease risk. Third, association with expression levels provides a unique opportunity to identify the actual disease gene at the linkage or association locus. Fourth, individual or combined assessment of glutathione pathway genes that are regulated in the brain, may uncover additional neurodegenerative risk variants. Further establishment of GSTO2 and other glutathione metabolism genes in AD and PD awaits discovery and mechanistic studies of functional genetic variants.
Subjects and samples
LOAD and PD DNA samples
Unrelated subjects from six independent LOAD case-control series, consisting of Caucasians with an age-at-diagnosis (LOAD), evaluation (elderly controls) or death (autopsy series) ≥ 60 years, were utilized in this study (3,561 LOAD vs. 4,683 controls; Table 1 Table 2 and Table 3). Subjects with younger (60-80) and older (> 80) ages were assessed both separately and jointly, as per our prior reports[40–42] and given the age-specific effects observed for many LOAD risk variants, including APOE. Four case-control series were comprised of Caucasian subjects collected in the United States, with three series collected at Mayo Clinic in Jacksonville, Florida (JS: 864 LOADs, 972 controls), Rochester, Minnesota, (RS: 597 LOADs, 2,406 controls) and an autopsy-confirmed series from the Brain Bank at Mayo Clinic Florida (AUT: 581 LOADs, 360 controls). The fourth Caucasian-American series was from the National Cell Repository for Alzheimer's Disease (NCRAD: 695 LOADs, 208 controls). These series were previously described in detail. Two additional Caucasian series were from Poland (PS: 479 LOADs, 186 controls) and Norway (NW: 345 LOADs, 551 controls) were also included in this study. All clinical LOAD subjects had a diagnosis of probable or possible AD and all autopsied LOAD subjects of definite AD made according to NINCDS-ADRDA criteria. All controls from the clinical Caucasian-American series had a clinical dementia rating score of 0. All autopsied LOAD brains had Braak scores of ≥ 4.0. Brains employed as controls had Braak scores of ≤ 2.5 but often had pathologies unrelated to AD.
We utilized age 80 as the arbitrary cutoff to define the older vs. younger LOAD series. This decision is partly based on the knowledge that the strongest genetic factor conferring LOAD risk, APOE ε4 has age-specific effects with highest effect sizes in younger ages (60-75) (reviewed in ). Age-specific LOAD risk association has also been demonstrated for other genetic factors both by others[48, 49] and by our group. These results suggest that older vs. younger LOAD subjects may be heterogeneous and harbor different genetic risk factors. Consequently, we have divided our LOAD case-control series into older vs. younger age groups using the 80 year cutoff and analyze these series both separately and jointly in all of our studies assessing LOAD genetic risk, including the Mayo LOAD GWAS which was focused on the 60-80 year group.
Caucasian-American, unrelated PD patients and controls (PD: 678 PDs, 712 controls) were recruited and diagnosed as described[52, 53], by a neurologist according to published criteria. Control subjects lacked any history suggestive of parkinsonism. PDs with family history of parkinsonism (familial, PD FAM) and those without (sporadic, PD-SPO) were analyzed both separately and jointly. All DNA samples were isolated from peripheral blood, with the exception of samples in the autopsy series where DNA was isolated from donated brain tissue, as described in previous publications[40, 44, 45, 52]. This study was approved by the appropriate institutional review board and appropriate informed consent was obtained from all participants.
Brain RNA for gene expression studies was obtained from the Mayo Clinic Autopsy (AUT) series, described above. These subjects were part of a larger expression GWAS (eGWAS) (in-press Neurology and PloSGenetics). AUT subjects with an age-at-death of 60-80 years were included in the Mayo LOAD GWAS. RNA was extracted from the frozen cerebellum and temporal cortex samples of these autopsied subjects, where available, using the AB(Applied Biosystems) RNA was extracted from the frozen cerebellum and temporal cortex samples of these autopsied subjects, where available, using the Ambion RNAqueous kit according to the manufacturer’s instructions. The quantity and quality of the RNA samples were determined by the Agilent 2100 Bioanalyzer using the Agilent RNA 6000 Nano Chip. In total, 399 temporal cortex (202 LOADs, 197 Controls) and 374 cerebellar samples (197 LOADs, 177 Controls) were assessed.
This study initially focused on four known coding variants from dbSNP within GSTO1 (rs4925, rs11509438) and GSTO2 (rs156697, rs34400162). We determined that rs11509438 was below the required minor allele frequency (MAF) cutoff of 5% (MAF = 3.3%) and rs34400162 was monomorphic in our series. We therefore focused on rs4925 and r156697 in all downstream analyses. Taqman (Applied Biosystems) was used to genotype rs4925 and rs155697 in all case-control series.
The genotypes for the eGWAS were obtained as part of the Mayo LOAD GWAS using Illumina (San Diego, CA) HumanHap300-Duo Genotyping BeadChips, and were analyzed with an Illumina BeadLab Station (Illumina, San Diego, CA), followed by quality control (QC), as previously described.
Transcript levels were measured using the Whole Genome DASL assay (Illumina, San Diego, CA). The RNA samples were randomized across the chips and plates using a stratified approach to ensure balance with respect to diagnosis, age, gender, RINs (RNA integrity numbers) and APOE genotype. Replicate samples were utilized for QC. Raw probe level mRNA expression data were exported from GenomeStudio software (Illumina Inc.) for preprocessing with background correction, variance stabilizing transformation, quantile normalization and probe filtering using the lumi package of BioConductor[55, 56].
Disease-risk association analysis
SNPs rs4925 and rs156697 were assessed for association with LOAD by multivariate logistic regression analysis using an allelic dosage model, adjusted for the following covariates: APOE ε4 dosage (0, 1, 2), age at diagnosis/evaluation/death, and gender. Analysis was executed for each of the six individual case control series and for the series combined, with a series covariate included in the model. All analyses were conducted separately on subjects of ages 60-80 years and those > 80 years, separately, and jointly.
Analysis for PD risk was done in a similar fashion. The PDs were classified as familial vs. sporadic based on the presence of family history of parkinsonism, in this cohort. We have therefore assessed the familial and sporadic PDs against the common PD control group both separately and jointly.
Meta analysis was also performed for rs4925 and rs156997 association with LOAD risk using the DerSimonian-Laird random effects model. The younger (60-80) and older (> 80) age groups were analyzed separately, in addition to combined ages. Breslow-Day test for non-compatibility was used to test for series heterogeneity. Test statistics are reported for each series as well as the pooled test statistics from the random effects model.
Age-at-diagnosis association analysis
We employed age-at-diagnosis or death as the surrogate quantitative variable for age-at-onset in our LOAD subjects (3,561 LOADs). An additive model was assumed for the GSTO locus SNPs, with the minor allele dosage (0, 1, 2) as the independent variable, and APOE ε4 dosage (0, 1, 2), and gender as covariates. Analysis was executed for the LOADs both individually for the six series and with all series combined, where a series covariate was included in the model.
PD age-at-diagnosis analysis was done similarly, for the familial and sporadic PDs both individually and jointly. Both the disease risk and age-at-diagnosis associations for the two GSTO SNPs were conducted in StatsDirect (v.2.5.8).
Gene expression level association analysis
All GSTO1 and GSTO2 trancript measurements and GSTO locus cisSNP/transcript associations were conducted as part of our brain eGWAS. Pathway analysis described below was also conducted using the cisSNP/transcript association results from this eGWAS. GSTO locus SNPs rs4925 and rs156697 were tested for association with brain levels of GSTO1 and GSTO2 transcripts measured in the cerebellum and temporal cortex of autopsied subjects as part of our eGWAS. Preprocessed probe transcript levels (GSTO1 = ILMN_2227573 and GSTO2 = ILMN_1740234) were used as the quantitative traits and the analyses were conducted as described (in-press Neurology, PLoS Genetics). An additive model was assumed, with the minor allele dosage (0, 1, 2) as the independent variable, and APOE ε4 dosage (0, 1, 2), age-at-death, gender, PCR plate, RIN, adjusted RIN2 (defined as (RIN-RINmean)2) as covariates. The cerebellar and temporal cortex results were analyzed separately. The autopsied LOADs and controls without AD pathology were analyzed both separately and jointly for rs4925 and rs156697, with the joint analysis including diagnosis as an additional covariate. Linear regression analysis to test for SNP/transcript associations were done with PLINK. Box plots depicting transcript levels by SNP genotype were generated in R, for the residuals of the cerebellar and temporal cortex associations from the multivariate linear regression analysis described above.
Our brain eGWAS assessed association of brain transcript levels with their nearby cisSNPs described as those residing within the gene or its ± 100 kb flanking region. We determined that there were 20 additional cisSNPs tested for GSTO1 and GSTO2 in our eGWAS. The brain transcript level associations with these 20 cisSNPs were also assessed using the same analytical approach.
Linkage disequilibrium analysis
Linkage disequilibrium for the 20 GSTO locus cisSNPs, rs4925 and rs156697 was evaluated using the HapMap Caucasian (CEU) data and assessed in HaploView version 4.1 with the solid spine algorithm.
In our brain eGWAS, we measured expression levels of 24,526 transcripts in 773 brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (control, cerebellar n = 177, temporal cortex n = 197) (in-press, PLoS Genetics). Association studies were carried out for the transcripts that could be detected in brain tissue, which is ~70% of all tested transcripts and for the 213,528 cisSNPs which reside within ± 100 kb of the genomic region of the transcript. We identified 686 genes with cerebellar transcript levels that are significantly influenced by cisSNPs in both the AD and control samples. Importantly, 625 of these genes could be tested in the temporal cortex, of which 471 were also significant for the transcript associations in this other brain region.
To discover the molecular pathways which harbor the top genes with significant cisSNP associations, we performed pathway analysis using MetaCore. The total number of tested genes and the number of significant genes that belong to a MetaCore pathway were used to determine the molecular pathways that are significantly enriched for genes influenced by cisSNPs in the brain. Out of 11,897 tested genes with eGWAS data, 3,316 belong to at least one MetaCore pathway. Out of 686 genes with significant cerebellar transcript/cisSNP associations, 188 belong to at least one MetaCore pathway. Pathways with less than five tested genes were excluded from analysis. Fisher 2 × 2 test with mid-p value was calculated in R to determine significance of enrichment.
Gawryluk JW, Wang JF, Andreazza AC, Shao L, Yatham LN, Young LT: Prefrontal cortex glutathione S-transferase levels in patients with bipolar disorder, major depression and schizophrenia. Int J Neuropsychopharmacol. 2011, 14: 1069-1074. 10.1017/S1461145711000617.
Andreazza AC, Kapczinski F: Kauer-Sant'Anna M, Walz JC, Bond DJ, Goncalves CA, Young LT, Yatham LN: 3-Nitrotyrosine and glutathione antioxidant system in patients in the early and late stages of bipolar disorder. J Psychiatry Neurosci. 2009, 34: 263-271.
Gravina P, Spoletini I, Masini S, Valentini A, Vanni D, Paladini E, Bossu P, Caltagirone C, Federici G, Spalletta G, Bernardini S: Genetic polymorphisms of glutathione S-transferases GSTM1, GSTT1, GSTP1 and GSTA1 as risk factors for schizophrenia. Psychiatry Res. 2011, 187: 454-456. 10.1016/j.psychres.2010.10.008.
Rodriguez-Santiago B, Brunet A, Sobrino B, Serra-Juhe C, Flores R, Armengol L, Vilella E, Gabau E, Guitart M, Guillamat R, et al: Association of common copy number variants at the glutathione S-transferase genes and rare novel genomic changes with schizophrenia. Mol Psychiatry. 2010, 15: 1023-1033. 10.1038/mp.2009.53.
Li YJ, Oliveira SA, Xu P, Martin ER, Stenger JE, Scherzer CR, Hauser MA, Scott WK, Small GW, Nance MA, et al: Glutathione S-transferase omega-1 modifies age-at-onset of Alzheimer disease and Parkinson disease. Hum Mol Genet. 2003, 12: 3259-3267. 10.1093/hmg/ddg357.
Kolsch H, Linnebank M, Lutjohann D, Jessen F, Wullner U, Harbrecht U, Thelen KM, Kreis M, Hentschel F, Schulz A, et al: Polymorphisms in glutathione S-transferase omega-1 and AD, vascular dementia, and stroke. Neurology. 2004, 63: 2255-2260. 10.1212/01.WNL.0000147294.29309.47.
Li YJ, Scott WK, Zhang L, Lin PI, Oliveira SA, Skelly T, Doraiswamy MP, Welsh-Bohmer KA, Martin ER, Haines JL, et al: Revealing the role of glutathione S-transferase omega in age-at-onset of Alzheimer and Parkinson diseases. Neurobiol Aging. 2006, 27: 1087-1093. 10.1016/j.neurobiolaging.2005.05.013.
van de Giessen E, Fogh I, Gopinath S, Smith B, Hu X, Powell J, Andersen P, Nicholson G: Al Chalabi A, Shaw CE: Association study on glutathione S-transferase omega 1 and 2 and familial ALS. Amyotroph Lateral Scler. 2008, 9: 81-84. 10.1080/17482960701702553.
Capurso C, Panza F, Seripa D, Frisardi V, Imbimbo BP, Verdile G, Vendemiale G, Pilotto A, Solfrizzi V: Polymorphisms in glutathione S-transferase omega-1 gene and increased risk of sporadic Alzheimer disease. Rejuvenation Res. 2010, 13: 645-652. 10.1089/rej.2010.1052.
Maes OC, Schipper HM, Chong G, Chertkow HM, Wang E: A GSTM3 polymorphism associated with an etiopathogenetic mechanism in Alzheimer disease. Neurobiol Aging. 2010, 31: 34-45. 10.1016/j.neurobiolaging.2008.03.007.
Lee JH, Mayeux R, Mayo D, Mo J, Santana V, Williamson J, Flaquer A, Ciappa A, Rondon H, Estevez P, et al: Fine mapping of 10q and 18q for familial Alzheimer's disease in Caribbean Hispanics. Mol Psychiatry. 2004, 9: 1042-1051. 10.1038/sj.mp.4001538.
Blacker D, Bertram L, Saunders AJ, Moscarillo TJ, Albert MS, Wiener H, Perry RT, Collins JS, Harrell LE, Go RC, et al: Results of a high-resolution genome screen of 437 Alzheimer's Disease families. Hum Mol Genet. 2003, 12: 23-32. 10.1093/hmg/ddg007.
Li YJ, Scott WK, Hedges DJ, Zhang F, Gaskell PC, Nance MA, Watts RL, Hubble JP, Koller WC, Pahwa R, et al: Age at onset in two common neurodegenerative diseases is genetically controlled. Am J Hum Genet. 2002, 70: 985-993. 10.1086/339815.
Whitbread AK, Masoumi A, Tetlow N, Schmuck E, Coggan M, Board PG: Characterization of the omega class of glutathione transferases. Methods Enzymol. 2005, 401: 78-99.
Laliberte RE, Perregaux DG, Hoth LR, Rosner PJ, Jordan CK, Peese KM, Eggler JF, Dombroski MA, Geoghegan KF, Gabel CA: Glutathione s-transferase omega 1-1 is a target of cytokine release inhibitory drugs and may be responsible for their effect on interleukin-1beta posttranslational processing. J Biol Chem. 2003, 278: 16567-16578. 10.1074/jbc.M211596200.
Wahner AD, Glatt CE, Bronstein JM, Ritz B: Glutathione S-transferase mu, omega, pi, and theta class variants and smoking in Parkinson's disease. Neurosci Lett. 2007, 413: 274-278. 10.1016/j.neulet.2006.11.053.
Nishimura M, Sakamoto T, Kaji R, Kawakami H: Influence of polymorphisms in the genes for cytokines and glutathione S-transferase omega on sporadic Alzheimer's disease. Neurosci Lett. 2004, 368: 140-143. 10.1016/j.neulet.2004.06.076.
Ozturk A, Desai PP, Minster RL, Dekosky ST, Kamboh MI: Three SNPs in the GSTO1, GSTO2 and PRSS11 genes on chromosome 10 are not associated with age-at-onset of Alzheimer's disease. Neurobiol Aging. 2005, 26: 1161-1165. 10.1016/j.neurobiolaging.2004.11.001.
Ertekin-Taner N: Genetics of Alzheimer's disease: a centennial review. Neurol Clin. 2007, 25: 611-667. 10.1016/j.ncl.2007.03.009.
Ertekin-Taner N: Genetics of Alzheimer disease in the pre- and post-GWAS era. Alzheimers Res Ther. 2010, 2: 3-10.1186/alzrt26.
Wider C, Wszolek ZK: Clinical genetics of Parkinson's disease and related disorders. Parkinsonism Relat Disord. 2007, 13 (Suppl 3): S229-S232.
Golub Y, Berg D, Calne DB, Pfeiffer RF, Uitti RJ, Stoessl AJ, Wszolek ZK, Farrer MJ, Mueller JC, Gasser T, Fuchs J: Genetic factors influencing age at onset in LRRK2-linked Parkinson disease. Parkinsonism Relat Disord. 2009, 15: 539-541. 10.1016/j.parkreldis.2008.10.008.
Haugarvoll K, Wszolek ZK: Clinical features of LRRK2 parkinsonism. Parkinsonism Relat Disord. 2009, 15 (Suppl 3): S205-S208.
Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, Hamshere ML, Pahwa JS, Moskvina V, Dowzell K, Williams A, et al: Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat Genet. 2009, 41: 1088-1093. 10.1038/ng.440.
Lambert JC, Heath S, Even G, Campion D, Sleegers K, Hiltunen M, Combarros O, Zelenika D, Bullido MJ, Tavernier B, et al: Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet. 2009, 41: 1094-1099. 10.1038/ng.439.
Seshadri S, Fitzpatrick AL, Ikram MA, DeStefano AL, Gudnason V, Boada M, Bis JC, Smith AV, Carassquillo MM, Lambert JC, et al: Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA. 2010, 303: 1832-1840. 10.1001/jama.2010.574.
Naj AC, Jun G, Beecham GW, Wang LS, Vardarajan BN, Buros J, Gallins PJ, Buxbaum JD, Jarvik GP, Crane PK, et al: Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer's disease. Nat Genet. 2011, 43: 436-441. 10.1038/ng.801.
Hollingworth P, Harold D, Sims R, Gerrish A, Lambert JC, Carrasquillo MM, Abraham R, Hamshere ML, Pahwa JS, Moskvina V, et al: Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease. Nat Genet. 2011, 43: 429-435. 10.1038/ng.803.
Satake W, Nakabayashi Y, Mizuta I, Hirota Y, Ito C, Kubo M, Kawaguchi T, Tsunoda T, Watanabe M, Takeda A, et al: Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson's disease. Nat Genet. 2009, 41: 1303-1307. 10.1038/ng.485.
Simon-Sanchez J, Schulte C, Bras JM, Sharma M, Gibbs JR, Berg D, Paisan-Ruiz C, Lichtner P, Scholz SW, Hernandez DG, et al: Genome-wide association study reveals genetic risk underlying Parkinson's disease. Nat Genet. 2009, 41: 1308-1312. 10.1038/ng.487.
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 2009, 461: 747-753. 10.1038/nature08494.
Liu H, Wang H, Shenvi S, Hagen TM, Liu RM: Glutathione metabolism during aging and in Alzheimer disease. Ann N Y Acad Sci. 2004, 1019: 346-349. 10.1196/annals.1297.059.
Ekins S, Nikolsky Y, Bugrim A, Kirillov E, Nikolskaya T: Pathway mapping tools for analysis of high content data. Methods Mol Biol. 2007, 356: 319-350.
Schmuck EM, Board PG, Whitbread AK, Tetlow N, Cavanaugh JA, Blackburn AC, Masoumi A: Characterization of the monomethylarsonate reductase and dehydroascorbate reductase activities of Omega class glutathione transferase variants: implications for arsenic metabolism and the age-at-onset of Alzheimer's and Parkinson's diseases. Pharmacogenet Genomics. 2005, 15: 493-501. 10.1097/01.fpc.0000165725.81559.e3.
Whitbread AK, Tetlow N, Eyre HJ, Sutherland GR, Board PG: Characterization of the human Omega class glutathione transferase genes and associated polymorphisms. Pharmacogenetics. 2003, 13: 131-144. 10.1097/00008571-200303000-00003.
Dulhunty A, Gage P, Curtis S, Chelvanayagam G, Board P: The glutathione transferase structural family includes a nuclear chloride channel and a ryanodine receptor calcium release channel modulator. J Biol Chem. 2001, 276: 3319-3323. 10.1074/jbc.M007874200.
Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, de Bakker PI: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008, 24: 2938-2939. 10.1093/bioinformatics/btn564.
Maraganore DM, de Andrade M, Lesnick TG, Strain KJ, Farrer MJ, Rocca WA, Pant PV, Frazer KA, Cox DR, Ballinger DG: High-resolution whole-genome association study of Parkinson disease. Am J Hum Genet. 2005, 77: 685-693. 10.1086/496902.
Lill CM, Roehr JT, McQueen MB, Kavvoura F, Bagade S, Schjeide BMM, Schjeide L, Meissner E, Zauft U, Allen NC, et al: The PDGene Database. Alzheimer Research Forum. Available at: http://www.pdgene.org/. Accessed [03-04-2012]
Zou F, Carrasquillo MM, Pankratz VS, Belbin O, Morgan K, Allen M, Wilcox SL, Ma L, Walker LP, Kouri N, et al: Gene expression levels as endophenotypes in genome-wide association studies of Alzheimer disease. Neurology. 2010, 74: 480-486. 10.1212/WNL.0b013e3181d07654.
Burgess JD, Pedraza O, Graff-Radford NR, Hirpa M, Zou F, Miles R, Nguyen T, Li M, Lucas JA, Ivnik RJ, et al: Association of common KIBRA variants with episodic memory and AD risk. Neurobiol Aging. 2011, 32: 557 e551-559.
Allen M, Cox C, Belbin O, Ma L, Bisceglio GD, Wilcox SL, Howell CC, Hunter TA, Culley O, Walker LP, et al: Association and heterogeneity at the GAPDH locus in Alzheimer's disease. Neurobiol Aging. 2012, 33: 203 e225-233.
Khachaturian AS, Corcoran CD, Mayer LS, Zandi PP, Breitner JC: Apolipoprotein E epsilon4 count affects age at onset of Alzheimer disease, but not lifetime susceptibility: The Cache County Study. Arch Gen Psychiatry. 2004, 61: 518-524. 10.1001/archpsyc.61.5.518.
Klimkowicz-Mrowiec A, Marona M, Wolkow P, Maruszak A, Styczynska M, Barcikowska M, Zekanowski C, Szczudlik A, Slowik A: Interleukin-1 gene -511 CT polymorphism and the risk of Alzheimer's disease in a Polish population. Dement Geriatr Cogn Disord. 2009, 28: 461-464. 10.1159/000259460.
Sando SB, Melquist S, Cannon A, Hutton ML, Sletvold O, Saltvedt I, White LR, Lydersen S, Aasly JO: APOE epsilon 4 lowers age at onset and is a high risk factor for Alzheimer's disease; a case control study from central Norway. BMC Neurol. 2008, 8: 9-10.1186/1471-2377-8-9.
McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan EM: Clinical diagnosis of Alzheimer's disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer's Disease. Neurology. 1984, 34: 939-944. 10.1212/WNL.34.7.939.
Farrer LA, Cupples LA, Haines JL, Hyman B, Kukull WA, Mayeux R, Myers RH, Pericak-Vance MA, Risch N, van Duijn CM: Effects of age, sex, and ethnicity on the association between apolipoprotein E genotype and Alzheimer disease. A meta-analysis. APOE and Alzheimer Disease Meta Analysis Consortium. JAMA. 1997, 278: 1349-1356. 10.1001/jama.1997.03550160069041.
Li Y, Nowotny P, Holmans P, Smemo S, Kauwe JS, Hinrichs AL, Tacey K, Doil L, van Luchene R, Garcia V, et al: Association of late-onset Alzheimer's disease with genetic variation in multiple members of the GAPD gene family. Proc Natl Acad Sci USA. 2004, 101: 15688-15693. 10.1073/pnas.0403535101.
Lin PI, Martin ER, Bronson PG, Browning-Large C, Small GW, Schmechel DE, Welsh-Bohmer KA, Haines JL, Gilbert JR, Pericak-Vance MA: Exploring the association of glyceraldehyde-3-phosphate dehydrogenase gene and Alzheimer disease. Neurology. 2006, 67: 64-68. 10.1212/01.wnl.0000223438.90113.4e.
Burgess JD, Pedraza O, Graff-Radford NR, Hirpa M, Zou F, Miles R, Nguyen T, Li M, Lucas JA, Ivnik RJ, et al: Association of common KIBRA variants with episodic memory and AD risk. Neurobiol Aging. 2011, 557: e1-9.
Carrasquillo MM, Zou F, Pankratz VS, Wilcox SL, Ma L, Walker LP, Younkin SG, Younkin CS, Younkin LH, Bisceglio GD, et al: Genetic variation in PCDH11X is associated with susceptibility to late-onset Alzheimer's disease. Nat Genet. 2009, 41: 192-198. 10.1038/ng.305.
Wider C, Vilarino-Guell C, Heckman MG, Jasinska-Myga B, Ortolaza-Soto AI, Diehl NN, Crook JE, Cobb SA, Bacon JA, Aasly JO, et al: SNCA, MAPT, and GSK3B in Parkinson disease: a gene-gene interaction study. Eur J Neurol. 2011, 18: 876-881. 10.1111/j.1468-1331.2010.03297.x.
Uitti RJ, Baba Y, Wszolek ZK, Putzke DJ: Defining the Parkinson's disease phenotype: initial symptoms and baseline characteristics in a clinical cohort. Parkinsonism Relat Disord. 2005, 11: 139-145. 10.1016/j.parkreldis.2004.10.007.
Gelb DJ, Oliver E, Gilman S: Diagnostic criteria for Parkinson disease. Arch Neurol. 1999, 56: 33-39. 10.1001/archneur.56.1.33.
Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.
Lin SM, Du P, Huber W, Kibbe WA: Model-based variance-stabilizing transformation for Illumina microarray data. Nucleic Acids Res. 2008, 36: e11-
DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7: 177-188. 10.1016/0197-2456(86)90046-2.
Allen M, Zou F, Chai HS, Younkin CS, Crook J, Pankratz VS, Carrasquillo MM, Rowley CN, Nair AA, Middha S, et al: Novel late-onset Alzheimer's disease loci variants associate with brain gene expression. Neurology. 2012,
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.
Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449: 851-861. 10.1038/nature06258.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
Support for this research was provided by the National Institutes of Health grants: National Institute on Aging (R01 032990 to NET and R01 AG018023 to NRG-R and SGY); Mayo Alzheimer's Disease Research Center: (P50 AG016574 to RCP, DWD, NRG-R, SGY, and NET); Mayo Alzheimer's Disease Patient Registry: (U01 AG006576 to RCP); National Institute on Aging (AG025711, AG017216, AG003949 to DWD). This project was also generously supported by the Robert and Clarice Smith and Abigail Van Buren Alzheimer's Disease Research Program (to RCP, DWD, NRG-R, and SGY), and by the Palumbo Professorship in Alzheimer's Disease Research (to SGY). N.E.T. is the recipient of National Institutes of Health (KL2 RR024151), and Siragusa Foundation grants. RJU and ZKW are partially supported by the NIH/NINDS P50 NS072187-01 S2 and NS057567, and Mayo Clinic Florida (MCF) Research Committee CR program. ZKW is also partially supported by NIH/NINDS 1RC2NS070276 and Dystonia Medical Research Foundation.
N. Graff-Radford, M.D. has served as a consultant to Codman and received grant support from Elan Pharmaceutical Research, Pfizer Pharmaceuticals, Medivation, and Forrest. R.C. Petersen, M.D., Ph.D. has been a consultant to GE Healthcare and Elan Pharmaceuticals, has served on a data safety monitoring board in a clinical trial sponsored by Elan Pharmaceuticals, and a safety monitoring board for Wyeth Pharmaceuticals.
MA carried out SNP genotyping assays, executed the disease risk and age at onset statistical analysis and helped to draft the manuscript. FZ participated in the organization and design of the gene expression study. HSC, CSY and AAN significantly contributed to the analysis of the gene expression data. RM participated in the collection of the gene expression data and carried out SNP genotyping assays. JEC and VSP significantly contributed to the design of the gene expression study and statistical analysis. MMC led the design, collection and analysis of SNP genotyping data used for the gene expression study. CNR and RP carried out SNP genotyping assays. TN, LM, KGM, GB and AIO provided technical expertise in the preparation, organization and maintenance of the RNA and DNA samples. SM, SM and CG participated in the analysis of the gene expression data. DS and FR participated in the collection of gene expression data. CPK and JJ supervised the collection of gene expression data. SBS, JOA, MB, RJU, ZKW, OAR, RCP, NRG and DWD collected and provided the case-control samples used in these studies. SGY significantly contributed to the conception and organization of the study. NET conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.