Research article | Open | Published:
Rheumatoid arthritis-associated polymorphisms are not protective against Alzheimer's disease
Molecular Neurodegenerationvolume 6, Article number: 33 (2011)
Rheumatoid arthritis (RA) and Alzheimer's disease (AD) are inversely associated. To test the hypothesis that genetic elements associated with increased RA risk are associated with decreased AD risk, we evaluated RA genetic risk factors recently identified in genome-wide association studies (GWAS) for their association with AD in a two-stage, case-control analysis.
In our Stage 1 analysis of ~800 AD and ~1,200 non-AD individuals, three of seventeen RA-associated SNPs were nominally associated with AD (p < 0.05) with one SNP, rs2837960, retaining significance after correction for multiple testing (p = 0.03). The rs2837960_G (minor) allele, which is associated with increased RA risk, was associated with increased AD risk. Analysis of these three SNPs in a Stage 2 population, consisting of ~1,100 AD and ~2,600 non-AD individuals, did not confirm their association with AD. Analysis of Stage 1 and 2 combined suggested that rs2837960 shows a trend for association with AD. When the Stage 2 population was age-matched for the Stage 1 population, rs2837960 exhibited a non-significant trend with AD. Combined analysis of Stage 1 and the age-matched Stage 2 subset showed a significant association of rs2837960 with AD (p = 0.002, OR 1.24) that retained significance following correction for age, sex and APOE (p = 0.02, OR = 1.20). Rs2837960 is near BACE2, which encodes an aspartic protease capable of processing the AD-associated amyloid precursor protein. Testing for an association between rs2837960 and the expression of BACE2 isoforms in human brain, we observed a trend between rs2837960 and the total expression of BACE2 and the expression of a BACE2 transcript lacking exon 7 (p = 0.07 and 0.10, respectively).
RA-associated SNPs are generally not associated with AD. Moreover, rs2837960_G is associated with increased risk of both RA and, in individuals less than 80 years of age, with AD. Overall, these results contest the hypothesis that genetic variants associated with RA confer protection against AD. Further investigation of rs2837960 is necessary to elucidate the mechanism by which rs2837960 contributes to both AD and RA risk, likely via modulation of BACE2 expression.
There is a long-standing, inverse relationship between the prevalence of Alzheimer's disease (AD) and of rheumatoid arthritis (RA). Jenkinson and colleagues first described the decreased prevalence of RA in patients suffering from senile dementia of the Alzheimer's type as compared to cognitively intact individuals . Further retrospective studies of clinical and autopsy data revealed that patients with RA exhibit a reduced prevalence of AD . A study by Myllykangas-Luosujarvi and colleagues evaluating AD pathology in patients with and without RA revealed that AD-associated neuropathology occurred four times less often in patients with RA as compared to the general population .
The basis of this inverse relationship is unclear but may include both genetic and environmental factors. RA and AD each have a strong genetic component, i.e., 50% of RA risk and 60% of AD risk is attributable to genetic factors, supporting the original hypothesis of Jenkinson and colleagues that genetics might explain the relationship between AD and RA [4, 5]. Alternatively, anti-inflammatory medications used therapeutically for the treatment of RA could decrease AD risk by reducing AD-associated inflammation or via other mechanisms, .e.g., modulation of APP processing [6, 7]. Supporting this possibility, an initial double-blind, placebo-controlled study by Rogers et al. provided evidence that indomethacin slowed cognitive decline in patients with AD relative to placebo . These findings were further supported by Breitner and colleagues who found that multiple anti-inflammatory medications slow disease progression and delay disease onset . However, there has been little success replicating these findings in larger, randomized clinical trials [10–13]. Hence, whether anti-inflammatory agents delay the onset of AD remains unclear.
The recent advent of RA genome wide association studies (GWAS) has identified single nucleotide polymorphisms (SNP)s associated with RA that provide a foundation for evaluating the initial hypothesis of Jenkinson et al. that genetic variants that increase the risk of RA also decrease the risk of AD. To this end, we tested whether seventeen RA-associated SNPs with genome-wide significance were associated with AD in a two-stage analysis using separate AD case-control populations. We found that none of the seventeen alleles associated with increased RA risk were also associated with reduced AD risk. Rather, we found three RA-associated SNPs that were nominally associated with AD (p < 0.05). One of these SNPs, rs2837960, was found to be significantly associated with AD in a combined analysis of our Stage 1 and Stage 2 populations when the Stage 2 population was restricted to individuals of similar age as Stage 1. The gene closest to rs2837960 is BACE2, the product of which has been implicated in amyloid protein precursor (APP) processing [14, 15]. When we evaluated the expression of BACE2 isoforms as a function of rs2837960, we found a trend for BACE2 expression with rs2837960. In summary, genetic variants that increase RA risk do not decrease AD risk. The inverse relationship between RA and AD may thus be better explained by environmental factors such as the use of anti-inflammatory medications. Further functional investigation of rs2837960 is needed to elucidate the mechanism by which this SNP may modulate AD and RA.
RA-associated SNPs are generally not associated with AD
To evaluate whether RA-associated SNPs are also associated with AD, we began by identifying SNPs that are robustly associated with RA risk and then evaluated these SNPs for their association with AD in an exploratory Stage 1 case-control population of ~800 AD and ~1200 non-AD individuals. Contrary to the hypothesis that alleles associated with increased RA risk are also associated with reduced AD risk, only three of the seventeen RA-associated SNPs in our Stage 1 study were nominally significant for association with AD (p < 0.05, Table 1). Moreover, for each these SNPs, the allele associated with increased RA risk was also associated with increased AD risk, further refuting the hypothesis that genetics underlies the inverse epidemiologic relationship between RA and AD prevalence. The AD-associated SNPs are located in or near the genes BACE2 (rs2837960; p = 0.002, OR = 1.29), TRAF1/C5 (rs3761847; p = 0.006, OR = 1.19) and SALL3 (rs2002842; p = 0.04, OR = 1.15). When a Bonferroni correction for multiple testing was applied to minimize false-positive associations between RA-associated SNPs and AD, only rs2837960 exhibited a significant association with AD (p = 0.03). Furthermore, rs2837960 remained significantly associated with AD after correcting for age, sex and APOE genotype per logistic regression (p = 0.012 OR = 1.26, Table 2). Rs3761847 remained significant following correction of Stage 1 analysis for covariates (p = 0.007 OR = 1.21).
Rs2837960 is significantly associated with AD in individuals
The three RA- and potentially AD-associated SNPs from Stage 1 were evaluated further for their association with AD in a Stage 2 case-control series that consisted of 2677 non-AD and 1102 AD subjects. Contrary to our initial results, these SNPs exhibited no association with AD in our Stage 2 population (Table 3). Correction for age, sex and APOE genotype had marginal, non-significant effects on the association between these three SNPs and AD. Our Stage 1 and Stage 2 populations were combined with the intention of clarifying the overall association between these three SNPs and AD. Analysis of this combined population, consisting of 3949 non-AD and 1965 AD subjects, suggested that rs2837960 is significantly associated with AD (p = 0.04, OR = 1.11, Table 4). When these results were corrected for age, sex and APOE genotype this association between rs2837960 and AD showed only a trend (p = 0.14).
Due to the large discrepancy in average age between our Stage 1 and Stage 2 populations, and the possibility that the impact of genetic risk factors may decline with age, we next evaluated a subset of Stage 2 individuals that had an age of AD diagnosis between 60 and 80 years of age along with non-AD individuals with the same age range . Thus, our age parameters and average population ages for Stage 2 mimicked those of Stage 1. This effort resulted in a sample population that consisted of 186 AD (average age 73) and 912 non-AD individuals (average age of 74). This population showed a trend in the OR of rs2837960 that was consistent with that of the Stage 1 analysis (Table 5), i.e., the minor rs2837960_G allele appeared to impart an increased, although non-significant, risk of AD. The SNPs rs3761847 and rs2002842, that were associated with AD in Stage 1, failed to show an association with AD in Stage 2. We note that there was significant overlap in the 95% confidence intervals between Stage 1 and Stage 2 for all three of these SNPs. Analysis of the combined data from Stage 1 and the Stage 2 subset revealed that only rs2837960 maintained a significant association with AD (p = 0.002, OR = 1.27, Table 6). This association between rs2837960 and AD remained significant following correction for the covariates age, sex and APOE (p = 0.02, OR = 1.20). Thus, when the Stage 1 and Stage 2 populations were matched for age, rs2837960 was significantly associated with AD risk.
In silico analyses suggest that rs2837960 is associated with BACE2 expression and splicing
To gain insights into the possible actions of rs2837960, we first identified genes within the vicinity of rs2837960 by using HapMap . This analysis found that rs2837960 resides within a haplotypic block that includes the promoter region and 5'-UTR of BACE2 (Figure 1A). Due to the presence of distinct proximal and distal BACE2 promoters, and the potential contribution of upstream regulatory elements to BACE2 expression, we considered the possibility that rs3846662 may affect BACE2 expression [18–20]. We then queried whether there was an association between rs2837960 and BACE2 expression in human brain by using the SNPexpress database, which includes data from 93 human brain samples . This analysis revealed a trend towards increased BACE2 expression in rs2837960 minor allele carriers (p = 0.08, Figure 1B). Since these data suggest that rs2837960 may modulate BACE2 expression and since BACE2 function in turn may modulate AD and cognition, we hypothesized that rs2837960 modulates AD risk by altering BACE2 expression in the human brain [22–24].
Relationship between rs2837960 and BACE2 isoforms in the human brain
To pursue this hypothesis, we sought to directly quantify BACE2 isoforms as a function of rs2837960 in a series of human brain samples. We began by confirming the presence of previously identified BACE2 isoforms that lack exons seven or eight (BACE2d7 and BACE2d8, respectively) as well as the full-length, nine-exon BACE2 transcript (BACE2_FL). We also identified a novel isoform lacking both exons 7 and 8 (BACE2d7/8) that was confirmed by direct sequencing (Figure 2). Considering the biological relevance of these isoforms, we note that the loss of exon 7 is predicted to result in an in-frame deletion of 50 amino acids, resulting in a 50.3 kDa peptide. Moreover, the BACE2 protein found in human brain appears to correspond to the BACE2d7 isoform based upon its size and pattern of epitope reactivity . In contrast, the loss of exon 8 or exons 7-8 results in a frameshift and prematurely truncated BACE2 proteins of 43.0 kDa and 37.5 kDa, respectively. Hence, we chose to use real-time PCR to quantify BACE2tot and BACE2d7 since they represent the expression of total BACE2 and functional BACE2, respectively. A total of 53 brain cDNA samples were analyzed for genotypic associations between rs2837960 and expression of BACE2tot and BACEd7. This effort revealed a trend between rs2837960 and the expression of BACE2tot (p = 0.07) and BACE2d7 (p = 0.10, Figure 3). It is noteworthy that the trend toward greater expression of BACE2tot and BACE2d7 in rs2837960 minor allele carriers in our results is similar to the results from SNPExpress.
The primary finding of this investigation is that the majority of seventeen SNPs that exhibit a genome-wide significant association with RA are not associated with AD. Furthermore, the minor allele of rs2837960, which was found to be significantly associated with AD risk after combined analysis of Stage 1 and age-matched Stage 2 data, was associated with an increased risk of both RA and AD. Hence, these results contest the hypothesis that genetics underlie the inverse relationship between RA and AD, i.e. that alleles associated with an increased risk of RA are protective against AD. A secondary finding is that we have pursued the role of rs2837960 in its possible regulation of the nearby BACE2 gene. We report the presence of multiple BACE2 isoforms in human brain and that rs2837960 shows a trend for association with BACE2tot and BACE2d7, which represent total BACE2 and functional BACE2, respectively . In summary, the genetic underpinnings of RA have negligible overlap with AD with the exception of rs2837960, which is associated with both RA and AD, possibly through its effects on BACE2 expression.
RA and AD each have a strong genetic component that accounts for approximately 50% and 60% of their risk, respectively [4, 5]. The remainder of RA and AD risk is likely derived from environmental influences. The vast majority of RA-associated SNPs implicate gene products involved in immune system processes. Chronic inflammation of the brain is a common feature of AD pathology, raising the possibility that RA-associated SNPs that influence immune system function could influence AD risk [25–27]. It is well established that some of the most strongly AD-associated genes, including CLU, CR1, TNF and CCR2, exhibit ontological association with immune system processes [28–36]. Hence, the impetus for pursuing genetic overlap between RA and AD is greater than that provided by their epidemiologic relationship alone. However, our results indicate that RA-associated SNPs, which pertain largely to gene products involved in immune system processes, are not associated with AD.
There are several possible interpretations of our primary findings. The lack of overlap between RA-associated SNPs and AD could be due in part to the tissue-specific expression of DNA and RNA binding proteins required to interact with these SNPs to manifest effects on gene expression . However, if any of the seventeen RA-associated SNPs included in this study are capable of modulating peripheral immune system activity, either alone or in combination with each other, then it is probable that their peripheral effects on the immune system would indirectly affect immune system activity within the CNS. Evidence supporting the ability of peripheral inflammation to modulate CNS inflammation has been reported previously . Therefore, if RA-associated SNPs are only functional in the periphery then their effects on immune system function and inflammation should manifest in the CNS, even if the same SNPs do not modulate endogenous immune system function within the brain.
What is yet unclear is whether RA-associated alleles actually propagate inflammation and, if so, why they would not be expected to increase, rather than decrease, AD risk. In fact, the results of our study suggest that alleles that increase RA risk may likewise increase AD risk, i.e., rs2738960_G increases risk of both RA and AD, while rs3761847_G and rs2002842_A show a similar trend. If these observations are replicated in future studies, alleles that are pro-inflammatory may emerge as risk factors for both RA and AD. More explicitly, considering the role of genetics and environment in RA and AD, these results suggest that RA genetics alone may enhance rather than reduce AD risk. Hence, the inverse epidemiologic relationship between RA and AD is likely explained by an environmental RA-associated influence. In this regard, McGeer et al. postulated that the reduced prevalence of AD in RA patients is related to the use of anti-inflammatory drugs for the treatment of RA . Multiple studies of anti-inflammatory agents have since been performed to test for their ability to modify AD risk and cognitive decline in AD patients, yielding mixed results [8, 39–42]. To some extent, variability in study outcome may be explained by the additional ability of a subset of anti-inflammatory medications to reduce production of the neurotoxic Aβ1-42 peptide . Further investigation is required to clarify the functional genetics of RA- and AD-associated SNPs and the role of anti-inflammatory medications in AD.
In pursuit of the functional genetics of rs2837960, which is associated with an increased risk of RA and AD, we investigated its association with the expression of BACE2 isoforms in human brain . Thus, our secondary finding is that the minor allele of rs2837960 showed a strong trend for association with increased expression of BACE2tot and BACE2d7, the latter of which may represent the majority of functional BACE2 in human brain . BACE2 encodes a transmembrane aspartic protease and is ~75% homologous with BACE1 with regard to amino acid sequence . Although the function of BACE2 is disputed, it appears to possess both β-secretase and α-secretase-like activities . Data obtained from the study of BACE1/BACE2 double-knockout mice suggest that BACE2 expressed in glia contributes significantly to Aβ production . This glial-specific expression is likely due to the more distal of the two distinct BACE2 promoters, neither of which share similarity with the BACE1 promoter [18–20].
Several factors are consistent with the possibility that rs2837960, or SNPs in tight linkage with rs2837960 (LD of r2 >0.8), are functional in modulating BACE2 expression. This evidence includes the observation that (i) rs2837960 resides within a haplotypic block that spans the region containing both the proximal and distal BACE2 promoters as well as the 5'UTR and first exon of BACE2, (ii) the region surrounding rs2837960 and its proxy SNPs (r2 = 1.0, ~4 kb window) is well conserved in primates per rVISTA analysis (data not shown), and (iii) the alleles of rs2837960 and its proxy SNPs are predicted to differentially affect transcription factor binding per PROMO 3.0 analysis of the TRANSFAC database (data not shown) [45–47].
Other studies that examined the association between BACE2 polymorphisms and AD risk have yielded mixed results [48–54]. These studies differ with our study in that (i) they have focused on SNPs much more proximal to BACE2 that are not in strong linkage disequilibrium with rs2837960 and (ii) they generally utilized smaller populations than those utilized in our present study. Future analyses of the association between BACE2 SNPs and AD should therefore take into account SNPs that are more distal to BACE2, such as rs2837960, as well as utilize larger population sizes that are sufficiently powered to detect associations with AD. Thus, in future studies rs2837960 may emerge as a risk factor for both RA and AD that functionally modulates BACE2 expression. Elucidation of the precise mechanism by which rs2837960, or a SNP that is proxy to it, modulates BACE2 expression may contribute to a better understanding of the role of BACE2 in both AD and RA pathology.
In summary, we have provided evidence that RA genetics do not underlie the inverse relationship between RA and AD but rather may promote AD. Furthermore, we have found that rs2837960 is associated with both RA and AD and that it, or one of its proxy SNPs, may modulate the expression of BACE2. As we learn more about the pathologic processes behind both RA and AD, including the contribution of BACE2 to each disease, a greater understanding of the factors underlying the inverse relationship between these two diseases may be obtained.
Materials and methods
The Human Genome (HuGE) Navigator (http://www.hugenavigator.net) was queried using the search term "rheumatoid arthritis" to identify RA-associated SNPs of genome-wide significance . Six available studies utilizing individuals of European decent were chosen to mimic the AD MAYO GWAS demographics (Table 7). Sample sizes ranged from ~1,600 (810 RA, 794 non-RA) to ~25,500 (7,322 RA, 18,262 non-RA). Thus, we identified twenty-eight candidate SNPs for study from the literature. SNPs which appeared more than once or that were in tight linkage disequilibrium with each other, i.e. r2 >0.8 (according to the CEU HapMap population), were considered to be redundant and only those with the lowest RA-associated p-value were retained for further analysis . This effort reduced the number of candidate SNPs to twenty-two. If a candidate RA-associated SNP was not available within the Mayo Clinic AD GWAS, an appropriate proxy SNP (LD of r2 >0.8) was selected by using the HapMap-based SNAP proxy search (http://www.broadinstitute.org/mpg/snap/) . Ultimately, seventeen of the candidate SNPs or their proxies were present in our AD GWAS dataset.
Case and Control Samples
The Mayo Clinic case-control samples used for the Stage 1 analysis have been described in detail in a prior GWAS publication . The Mayo Clinic case-control series used for the Stage 2 study have also been previously described . Briefly, clinical diagnoses of probable AD were made according to NINCDS-ADRDA criteria for samples from Jacksonville, FL (JS) and Rochester, MN (RS); age-matched controls had a score of 0 on the Clinical Dementia Rating scale. Additional samples were obtained from the Mayo Clinic brain bank (AUT); autopsy-confirmed diagnosis of AD (NINCDS-ADRDA, Braak score >4.0) was utilized for AD samples while non-AD samples exhibited limited AD pathology (Braak <2.5, not including other unrelated pathology).
AD Association Testing
Association testing of RA-associated SNPs for AD risk was carried out in two stages by using PLINK software (http://pngu.mgh.harvard.edu/purcell/plink/) . All genotyped samples were subject to strict quality control including elimination of samples with call rates <90%, MAF <0.01, HW p < 0.001, discrepancy between reported and genotyped sex, cryptic relatedness and discordant genotype clustering. Stage 1 consisted of 1264 non-AD and 843 AD subjects with average ages of 74.3 ± 4.5 (age at last assessment, mean ± SD) and 72.4 ± 4.6 years (age at diagnosis), respectively. The non-AD and AD groups in this series consisted of 51.7% and 57.5% female individuals, respectively. Stage 1 samples were genotyped by using HumanHap300-Duo Genotyping BeadChips processed with an Illumina BeadLab station (Illumina, San Diego, CA) at the Mayo Clinic Genotyping Shared Resource center (Rochester, MN).
We proceeded to test for an association between the seventeen RA-associated SNPs and AD in this Stage 1 case-control population. Stage 1 association testing was performed by using PLINK to generate allelic models that included odds ratios (OR), 95% confidence intervals (CI)s and uncorrected p-values. Logistic regression was also performed using the covariates age, sex and APOE genotype. With regards to multiple testing we expected to obtain approximately one false positive result given α = 0.05 (seventeen unique SNPs; 17 tests × 0.05 = 0.85). Bonferroni correction for multiple testing was also applied to data generated using allelic models.
Stage 2 samples were genotyped by using SEQUENOM MassARRAY iPLEX Platform (Sequenom, San Diego, CA). Overall, Stage 2 consisted of 2677 non-AD and 1102 AD subjects with average ages of 81.0 ± 6.2 and 83.5 ± 6.6 years of age, respectively. The non-AD and AD groups were composed of 55.0% and 64.0% female individuals, respectively. Stage 2 AD-SNP association testing was performed using only the three SNPs identified in Stage 1 as being associated with both RA and AD. PLINK software was used to generate odds ratios, 95% CIs and p-values per allelic modeling. Logistic regression including the covariates age, sex and APOE genotype was also performed. To evaluate the overall significance of Stage 1 and 2 data, they were combined and examined collectively.
Due to the considerable difference in mean age between Stage 1 and Stage 2 individuals, and due to our interest in focusing on genetic, rather than environmental factors, we also chose to examine only Stage 2 individuals between 60 and 80 years of age. Hence, when Stage 2 was limited to individuals between 60 and 80 years of age, our analysis included 912 non-AD and 186 AD subjects with average ages of 73.9 ± 3.8 and 72.8 ± 5.1 years. The non-AD and AD groups consisted of 49.9% and 57.0% female individuals, respectively. Similar to our analysis of our initial Stage 2 population, logistic regression of this modified Stage 2 population was also performed to test for an association between the three AD-associated SNPs from Stage 1. Furthermore, we evaluated the overall significance of RA-associated SNP associations with AD in combined Stage 1 and Stage 2 individuals between 60 and 80 years of age.
Human anterior cingulate brain specimens were generously provided by the Sanders-Brown AD Center Neuropathology Core and have been described elsewhere . The samples were from deceased individuals with an average age at death for females of 82 ± 7 years (mean ± SD, n = 29) and for males of 81 ± 8 (n = 24). The average postmortem interval (PMI) for females and males was 3.2 ± 0.8 h and 3.0 ± 0.8 h, respectively.
Evaluation of BACE2 isoforms in vivo
To gain insights into the functionality of rs2837960 we tested for an association between rs2837960 and BACE2 expression in human brain. We first screened anterior cingulate samples for the presence of BACE2 and its known alternatively spliced isoforms that lack exons 7 and 8, respectively. Total RNA and genomic DNA were prepared from human tissue samples; the RNA was reverse transcribed as we have reported elsewhere [61, 62]. Conventional PCR using Platinum Taq (Invitrogen, Carlsbad, CA) was used to amplify the region of BACE2 spanning exons 6-9 (Table 8). Thermal cycling conditions consisted of denaturation at 95°C for 5 min followed by 32 cycles of 95°C for 30 s, 60°C for 30 s, 72°C for 1 min and a final extension at 72°C for 2 min. PCR products were separated using 8% TBE-PAGE gel electrophoresis and visualized using SYBR-gold fluorescent stain (Invitrogen) and a fluorescence imager (FLA-2000, Fuji). To confirm the identities of the BACE2 splice variants, bands were excised, purified and directly sequenced (Davis Sequencing, Davis, CA).
To quantify total BACE2 expression (BACE2tot) and expression of the BACE2 isoform lacking exon 7 (BACE2d7) we designed separate primer sets. BACE2tot expression was measured by amplification of a product spanning a non-alternatively spliced region of BACE2 (exons 5-6). Isoform-specific primers designed to amplify BACE2d7 consisted of a forward primer specific to BACE2 exon 6 and a reverse primer specific to the junction of exons 6-8 (Table 8). DNA samples were genotyped using a TaqMan SNP Genotyping Assay (ID # C_2688271_10; ABI, Carlsbad, CA).
Quantitative real-time PCR reactions contained ~20 ng of sample cDNA together with 10 μl of PerfeCTA SYBR green SuperMix (Quanta Biosciences, Gaithersburg, MD), 10 μl of ddH2O and 20 pmol of forward and reverse primers. Cycling conditions included a 3 minute denaturation step at 95°C followed by 40 cycles of denaturation for 15 seconds at 95°C and annealing/extension for 45 seconds at 60°C using an MJ Opticon 4 thermal cycler (Biorad, Hercules, CA). A melting curve was generated following cycling to assess the purity of amplification product. Fidelity of amplification was also assessed via visual inspection of PCR products on 8% TBE-PAGE gel stained with SYBR gold. Standard curves were generated from purified PCR products that were quantified by A260/280 spectrophotometric analysis. Standard curves were then used to calculate the copy number for each BACE2 isoform measured.
Hypoxanthine-guanine phosphoribosyltransferase (HPRT) and ribosomal protein L32 (RPL32) were used as housekeeping genes per the analysis of geNorm software as described previously [63–65]. Expression levels of each of these genes were measured by using quantitative real-time PCR and gene specific primers under conditions identical to cycling conditions for BACE2. Standard curves were used to generate exact copy numbers, which in turn were used to calculate the sample-specific geometric mean of HPRT and RPL32 expression. The geometric mean was in turn used to normalize subsequent BACE2 expression data. Analysis of the association between BACE2 isoforms and rs2837960 genotype was performed using non-parametric Jonckheere-Terpstra testing (PASW Statistics, v.18, IBM, Somers, NY).
Jenkinson ML, Bliss MR, Brain AT, Scott DL: Rheumatoid arthritis and senile dementia of the Alzheimer's type. Br J Rheumatol. 1989, 28: 86-88.
McGeer PL, McGeer E, Rogers J, Sibley J: Anti-inflammatory drugs and Alzheimer disease. Lancet. 1990, 335: 1037-
Myllykangas-Luosujarvi R, Isomaki H: Alzheimer's disease and rheumatoid arthritis. Br J Rheumatol. 1994, 33: 501-502. 10.1093/rheumatology/33.5.501.
Scott DL, Wolfe F, Huizinga TW: Rheumatoid arthritis. Lancet. 2010, 376: 1094-1108. 10.1016/S0140-6736(10)60826-4.
Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA, Berg S, Fiske A, Pedersen NL: Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry. 2006, 63: 168-174. 10.1001/archpsyc.63.2.168.
Weggen S, Eriksen JL, Das P, Sagi SA, Wang R, Pietrzik CU, Findlay KA, Smith TE, Murphy MP, Bulter T, et al: A subset of NSAIDs lower amyloidogenic Abeta42 independently of cyclooxygenase activity. Nature. 2001, 414: 212-216. 10.1038/35102591.
Eriksen JL, Sagi SA, Smith TE, Weggen S, Das P, McLendon DC, Ozols VV, Jessing KW, Zavitz KH, Koo EH, Golde TE: NSAIDs and enantiomers of flurbiprofen target gamma-secretase and lower Abeta 42 in vivo. J Clin Invest. 2003, 112: 440-449.
Rogers J, Kirby LC, Hempelman SR, Berry DL, McGeer PL, Kaszniak AW, Zalinski J, Cofield M, Mansukhani L, Willson P, et al: Clinical trial of indomethacin in Alzheimer's disease. Neurology. 1993, 43: 1609-1611.
Breitner JC, Gau BA, Welsh KA, Plassman BL, McDonald WM, Helms MJ, Anthony JC: Inverse association of anti-inflammatory treatments and Alzheimer's disease: initial results of a co-twin control study. Neurology. 1994, 44: 227-232.
Szekely CA, Green RC, Breitner JC, Ostbye T, Beiser AS, Corrada MM, Dodge HH, Ganguli M, Kawas CH, Kuller LH, et al: No advantage of A beta 42-lowering NSAIDs for prevention of Alzheimer dementia in six pooled cohort studies. Neurology. 2008, 70: 2291-2298. 10.1212/01.wnl.0000313933.17796.f6.
Szekely CA, Breitner JC, Fitzpatrick AL, Rea TD, Psaty BM, Kuller LH, Zandi PP: NSAID use and dementia risk in the Cardiovascular Health Study: role of APOE and NSAID type. Neurology. 2008, 70: 17-24. 10.1212/01.wnl.0000284596.95156.48.
Breitner JC, Haneuse SJ, Walker R, Dublin S, Crane PK, Gray SL, Larson EB: Risk of dementia and AD with prior exposure to NSAIDs in an elderly community-based cohort. Neurology. 2009, 72: 1899-1905. 10.1212/WNL.0b013e3181a18691.
Sonnen JA, Larson EB, Walker RL, Haneuse S, Crane PK, Gray SL, Breitner JC, Montine TJ: Nonsteroidal anti-inflammatory drugs are associated with increased neuritic plaques. Neurology. 2010, 75: 1203-1210. 10.1212/WNL.0b013e3181f52db1.
Ahmed RR, Holler CJ, Webb RL, Li F, Beckett TL, Murphy MP: BACE1 and BACE2 enzymatic activities in Alzheimer's disease. J Neurochem. 2010, 112: 1045-1053. 10.1111/j.1471-4159.2009.06528.x.
Bettegazzi B, Mihailovich M, Di Cesare A, Consonni A, Macco R, Pelizzoni I, Codazzi F, Grohovaz F, Zacchetti D: beta-Secretase activity in rat astrocytes: translational block of BACE1 and modulation of BACE2 expression. The European journal of neuroscience. 2011, 33: 236-243. 10.1111/j.1460-9568.2010.07482.x.
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.278.16.1349.
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.
Lahiri DK, Maloney B, Ge YW: Functional domains of the BACE1 and BACE2 promoters and mechanisms of transcriptional suppression of the BACE2 promoter in normal neuronal cells. J Mol Neurosci. 2006, 29: 65-80. 10.1385/JMN:29:1:65.
Maloney B, Ge YW, Greig NH, Lahiri DK: Characterization of the human beta-secretase 2 (BACE2) 5'-flanking region: identification of a 268-bp region as the basal BACE2 promoter. J Mol Neurosci. 2006, 29: 81-99. 10.1385/JMN:29:1:81.
Sun X, Wang Y, Qing H, Christensen MA, Liu Y, Zhou W, Tong Y, Xiao C, Huang Y, Zhang S, et al: Distinct transcriptional regulation and function of the human BACE2 and BACE1 genes. Faseb J. 2005, 19: 739-749. 10.1096/fj.04-3426com.
Sanders MA, Verhaak RG, Geertsma-Kleinekoort WM, Abbas S, Horsman S, van der Spek PJ, Lowenberg B, Valk PJ: SNPExpress: integrated visualization of genome-wide genotypes, copy numbers and gene expression levels. BMC Genomics. 2008, 9: 41-10.1186/1471-2164-9-41.
Farzan M, Schnitzler CE, Vasilieva N, Leung D, Choe H: BACE2, a beta-secretase homolog, cleaves at the beta site and within the amyloid-beta region of the amyloid-beta precursor protein. Proc Natl Acad Sci USA. 2000, 97: 9712-9717.
Stockley JH, Ravid R, O'Neill C: Altered beta-secretase enzyme kinetics and levels of both BACE1 and BACE2 in the Alzheimer's disease brain. FEBS Lett. 2006, 580: 6550-6560. 10.1016/j.febslet.2006.10.076.
Sloan CD, Shen L, West JD, Wishart HA, Flashman LA, Rabin LA, Santulli RB, Guerin SJ, Rhodes CH, Tsongalis GJ, et al: Genetic pathway-based hierarchical clustering analysis of older adults with cognitive complaints and amnestic mild cognitive impairment using clinical and neuroimaging phenotypes. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 1060-1069. 10.1002/ajmg.b.31078.
Amor S, Puentes F, Baker D, van der Valk P: Inflammation in neurodegenerative diseases. Immunology. 2010, 129: 154-169. 10.1111/j.1365-2567.2009.03225.x.
Perry VH, Nicoll JA, Holmes C: Microglia in neurodegenerative disease. Nat Rev Neurol. 2010, 6: 193-201. 10.1038/nrneurol.2010.17.
Metcalfe MJ, Figueiredo-Pereira ME: Relationship between tau pathology and neuroinflammation in Alzheimer's disease. Mt Sinai J Med. 2010, 77: 50-58. 10.1002/msj.20163.
Bertram L, McQueen MB, Mullin K, Blacker D, Tanzi RE: Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet. 2007, 39: 17-23. 10.1038/ng1934.
Xie Z, Harris-White ME, Wals PA, Frautschy SA, Finch CE, Morgan TE: Apolipoprotein J (clusterin) activates rodent microglia in vivo and in vitro. J Neurochem. 2005, 93: 1038-1046. 10.1111/j.1471-4159.2005.03065.x.
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.
Zanjani H, Finch CE, Kemper C, Atkinson J, McKeel D, Morris JC, Price JL: Complement activation in very early Alzheimer disease. Alzheimer Dis Assoc Disord. 2005, 19: 55-66. 10.1097/01.wad.0000165506.60370.94.
McGeer EG, McGeer PL: Neuroinflammation in Alzheimer's disease and mild cognitive impairment: a field in its infancy. J Alzheimers Dis. 2010, 19: 355-361.
McCoy MK, Tansey MG: TNF signaling inhibition in the CNS: implications for normal brain function and neurodegenerative disease. J Neuroinflammation. 2008, 5: 45-10.1186/1742-2094-5-45.
Conductier G, Blondeau N, Guyon A, Nahon JL, Rovere C: The role of monocyte chemoattractant protein MCP1/CCL2 in neuroinflammatory diseases. J Neuroimmunol. 2010
El Khoury J, Toft M, Hickman SE, Means TK, Terada K, Geula C, Luster AD: Ccr2 deficiency impairs microglial accumulation and accelerates progression of Alzheimer-like disease. Nat Med. 2007, 13: 432-438. 10.1038/nm1555.
Britschgi M, Wyss-Coray T: Immune cells may fend off Alzheimer disease. Nat Med. 2007, 13: 408-409. 10.1038/nm0407-408.
Heinzen EL, Ge D, Cronin KD, Maia JM, Shianna KV, Gabriel WN, Welsh-Bohmer KA, Hulette CM, Denny TN, Goldstein DB: Tissue-specific genetic control of splicing: implications for the study of complex traits. PLoS Biol. 2008, 6: e1-
Akiyama H, Barger S, Barnum S, Bradt B, Bauer J, Cole GM, Cooper NR, Eikelenboom P, Emmerling M, Fiebich BL, et al: Inflammation and Alzheimer's disease. Neurobiol Aging. 2000, 21: 383-421. 10.1016/S0197-4580(00)00124-X.
Aisen PS: Anti-inflammatory therapy for Alzheimer's disease. Neurobiol Aging. 2000, 21: 447-448. 10.1016/S0197-4580(00)00140-8. discussion 451-443
Zandi PP, Anthony JC, Hayden KM, Mehta K, Mayer L, Breitner JC: Reduced incidence of AD with NSAID but not H2 receptor antagonists: the Cache County Study. Neurology. 2002, 59: 880-886.
Aisen PS, Schafer KA, Grundman M, Pfeiffer E, Sano M, Davis KL, Farlow MR, Jin S, Thomas RG, Thal LJ: Effects of rofecoxib or naproxen vs placebo on Alzheimer disease progression: a randomized controlled trial. JAMA. 2003, 289: 2819-2826. 10.1001/jama.289.21.2819.
de Jong D, Jansen R, Hoefnagels W, Jellesma-Eggenkamp M, Verbeek M, Borm G, Kremer B: No effect of one-year treatment with indomethacin on Alzheimer's disease progression: a randomized controlled trial. PLoS One. 2008, 3: e1475-10.1371/journal.pone.0001475.
Solans A, Estivill X, de La Luna S: A new aspartyl protease on 21q22.3, BACE2, is highly similar to Alzheimer's amyloid precursor protein beta-secretase. Cytogenet Cell Genet. 2000, 89: 177-184. 10.1159/000015608.
Dominguez D, Tournoy J, Hartmann D, Huth T, Cryns K, Deforce S, Serneels L, Camacho IE, Marjaux E, Craessaerts K, et al: Phenotypic and biochemical analyses of BACE1- and BACE2-deficient mice. J Biol Chem. 2005, 280: 30797-30806. 10.1074/jbc.M505249200.
Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I: VISTA: computational tools for comparative genomics. Nucleic Acids Res. 2004, 32: W273-279. 10.1093/nar/gkh458.
Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM: PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002, 18: 333-334. 10.1093/bioinformatics/18.2.333.
Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, Messeguer X: Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003, 31: 3651-3653. 10.1093/nar/gkg605.
Giedraitis V, Kilander L, Degerman-Gunnarsson M, Sundelof J, Axelsson T, Syvanen AC, Lannfelt L, Glaser A: Genetic analysis of Alzheimer's disease in the Uppsala Longitudinal Study of Adult Men. Dementia and geriatric cognitive disorders. 2009, 27: 59-68. 10.1159/000191203.
Gold G, Blouin JL, Herrmann FR, Michon A, Mulligan R, Duriaux Sail G, Bouras C, Giannakopoulos P, Antonarakis SE: Specific BACE1 genotypes provide additional risk for late-onset Alzheimer disease in APOE epsilon 4 carriers. American journal of medical genetics Part B, Neuropsychiatric genetics: the official publication of the International Society of Psychiatric Genetics. 2003, 119B: 44-47.
Li H, Wetten S, Li L, St Jean PL, Upmanyu R, Surh L, Hosford D, Barnes MR, Briley JD, Borrie M, et al: Candidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer disease. Archives of neurology. 2008, 65: 45-53. 10.1001/archneurol.2007.3.
Myllykangas L, Wavrant-De Vrieze F, Polvikoski T, Notkola IL, Sulkava R, Niinisto L, Edland SD, Arepalli S, Adighibe O, Compton D, et al: Chromosome 21 BACE2 haplotype associates with Alzheimer's disease: a two-stage study. J Neurol Sci. 2005, 236: 17-24. 10.1016/j.jns.2005.04.008.
Nowotny P, Kwon JM, Chakraverty S, Nowotny V, Morris JC, Goate AM: Association studies using novel polymorphisms in BACE1 and BACE2. Neuroreport. 2001, 12: 1799-1802. 10.1097/00001756-200107030-00008.
Reiman EM, Webster JA, Myers AJ, Hardy J, Dunckley T, Zismann VL, Joshipura KD, Pearson JV, Hu-Lince D, Huentelman MJ, et al: GAB2 alleles modify Alzheimer's risk in APOE epsilon4 carriers. Neuron. 2007, 54: 713-720. 10.1016/j.neuron.2007.05.022.
Yu Y, Jia J: Lack of association between the polymorphisms of beta-site APP-cleaving enzyme 2 (BACE2) 5'-flanking region and sporadic Alzheimer's disease. Brain Res. 2009, 1257: 10-15.
Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ: A navigator for human genome epidemiology. Nature genetics. 2008, 40: 124-125. 10.1038/ng0208-124.
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.
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.
Carrasquillo MM, Belbin O, Hunter TA, Ma L, Bisceglio GD, Zou F, Crook JE, Pankratz VS, Dickson DW, Graff-Radford NR, et al: Replication of CLU, CR1, and PICALM associations with alzheimer disease. Arch Neurol. 2010, 67: 961-964. 10.1001/archneurol.2010.147.
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.
Zou F, Gopalraj RK, Lok J, Zhu H, Ling IF, Simpson JF, Tucker HM, Kelly JF, Younkin SG, Dickson DW, et al: Sex-dependent association of a common low-density lipoprotein receptor polymorphism with RNA splicing efficiency in the brain and Alzheimer's disease. Hum Mol Genet. 2008, 17: 929-935.
Zhu H, Tucker HM, Grear KE, Simpson JF, Manning AK, Cupples LA, Estus S: A common polymorphism decreases low-density lipoprotein receptor exon 12 splicing efficiency and associates with increased cholesterol. Hum Mol Genet. 2007, 16: 1765-1772. 10.1093/hmg/ddm124.
Ling IF, Gopalraj RK, Simpson JF, Estus S: Expression and regulation of a low-density lipoprotein receptor exon 12 splice variant. J Neurochem. 2010
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-
Ling IF, Estus S: Role of SFRS13A in low-density lipoprotein receptor splicing. Hum Mutat. 2010, 31: 702-709. 10.1002/humu.21244.
Zhang X, Ding L, Sandford AJ: Selection of reference genes for gene expression studies in human neutrophils by real-time PCR. BMC Mol Biol. 2005, 6: 4-10.1186/1471-2199-6-4.
Acknowledgements and Funding
This work was supported by grants from the National Institutes of Health/National Institute on Aging P01-AG030128 (SE), RO1-AG18023 (SGY), T32-AG000242 (CRS) and Mayo Alzheimer's Disease Research Center P50-AG16574 (SGY). Human tissue was generously provided by the University of Kentucky AD Center, which is supported by P30-AG028383.
The authors declare that they have no competing interests.
CRS and SE contributed to experimental design and wrote the manuscript. CRS performed RT-PCR of BACE2 isoforms, rs2837960 genotyping of human brain tissue and association testing between BACE2 isoform expression and rs2837960. FZ performed SEQUENOM genotyping of samples for Stage 2 analyses. SGY, FZ and CRS performed statistical analysis of RA-associated SNPs and AD GWAS data. All authors have read and approved the final manuscript.