Abstract
Keywords
INTRODUCTION
Alzheimer’s disease (AD) is a chronic neurodegenerative disease and is characterized by cognitive impairment including short-term memory loss and psychological symptoms such as depression, apathy, delusion, and other symptoms that are due to aggregation of amyloid-β (Aβ) plaques and neurofibrillary tangles. Because of reduced levels of acetylcholine in the brain of AD patients, acetylcholinesterase inhibitors (AChEIs) have been used as a first-line treatment for AD [1, 2]. However, the progression of AD cannot be stopped even if treatment with AChEIs continues. Therefore, a new, fundamental therapeutic strategy based on the original causes of AD is needed.
Several studies identified certain genotypes of apolipoprotein E (APOE) as a genetic risk factor for AD [3]. Large-scale genome-wide association studies with European AD subjects also identified various candidate genes for AD, including CR1, CLU, PICALM, BIN1, CD2AP, CD33, EPHA1, MS4A6A/MS4E4, and ABCA7 (ATP-binding cassette subfamily A member 7 gene) [4–8]. ABCA7 encodes a member of the superfamily of ATP-binding cassette transporters expressed on both microglia and neurons in human brain. ABCA7 is implicated in transport of phospholipids and enhances phagocytosis by macrophages [9–13]. Deletion of ABCA7 in amyloidogenic mice increases insoluble Aβ levels and thioflavine-S-positive plaques in the brain regardless of changes in amyloid-β protein precursor processing [12]. ABCA7 deficiency facilitates the processing of amyloid-β protein precursor and Aβ production by increasing the levels of β-secretase 1 in primary neurons and mouse brain [14]. Moreover, phagocytic clearance of Aβ oligomers in the hippocampus is reduced in ABCA7 deletion mice [15]. These results indicate that dysfunction of ABCA7 may cause accumulation of Aβ through impaired phagocytic function and increase the susceptibility to AD.
A single nucleotide polymorphism (SNP) in ABCA7, rs3764650, is associated with development of AD, and meta-analyses of all data in European studies conclude that ABCA7 rs3764650 is a susceptibility loci for AD (OR = 1.22) [4]. Although association studies investigating a link between rs3764650 and AD in Chinese (n = 633 and n = 1,224), Japanese (n = 1,735), Korean (n = 844), African American (n = 5,896), and Canadian (n = 1,104) individuals reported a weak or negligible association [16], a meta-analysis of these studies identified a significant association [17]. Moreover, ABCA7 expression is increased in AD brain, although the rs3764650 major allele T that decreases AD risk is associated with increased ABCA7 expression [18]. Additionally, methylation rates of three CpG sites in AD brain are associated with greater odds for AD diagnosis [19]. However, these sites are far (more than thousands of nucleotides) from the ABCA7 transcription start site (TSS), and the correlation between methylation rates of these CpG sites and ABCA7 expression is weak(r range, –0.2 to 0.2).
Although these results from mouse and human brain suggest that ABCA7 is related to AD pathogenesis, whether the gene polymorphism, expression, and methylation of ABCA7 are related to the development and psychological symptoms of AD has not been well investigated. In transgenic AD model mice, neutrophil deletion or inhibition of neutrophil trafficking via lymphocyte function-associated antigen-1 block reduces AD-like neuropathology and alleviates cognitive impairment [20]. This result demonstrates that neutrophils contribute to AD pathogenesis and cognitive impairment. Therefore, investigation of peripheral blood from AD patients may increase our understanding of AD pathogenesis and identify convenient clinical biomarkers [21–24].
The aim of this study was to examine the blood gene expression and methylation of ATP-binding cassette sub-family A member 7 gene (ABCA7) as a biological marker of AD.
METHODS
Subjects
Demographic data for each group of participants are shown in Table 1. We enrolled 50 AD subjects [11 males and 39 females, mean age±standard deviation (SD) = 77.7±6.05 years] who visited Ehime University Hospital, Zaidan Niihama Hospital, or Kuroda Hospital, Ehime, Japan from August 2013 to January 2016. AD subjects were diagnosed with and classified as having probable AD dementia according to criteria by the National Institute on Aging and the Alzheimer’s Association [25]. Of the 50 participants, 25 AD subjects took donepezil, and the other 25 did not take AChEIs at the time of blood sampling. AD subjects were evaluated with the Mini-Mental State Examination (MMSE) to assess their cognitive function [26], Clinical Dementia Rating (CDR) as measured by family caregivers [27], and Neuropsychiatric Inventory (NPI) to assess their psychological symptoms [28]. Control participants were 50 elderly individuals (11 males and 39 females, mean age±SD = 76.3±6.02 years) without cognitive impairment, psychiatric signs, or a past history of mental disorders as determined by at least two certified psychiatrists based on clinical interviews. All participants were unrelated, of Japanese origin, and provided written informed consent forms that were approved by the institutional ethics committees of Ehime University Hospital, Zaidan Niihama Hospital, and Kuroda Hospital.
Blood sample collection, extraction of genomic DNA (gDNA), and synthesis of complementary DNA (cDNA)
gDNA was obtained from whole peripheral blood samples collected in potassium EDTA tubes and extracted using the QIAamp DNA Blood Mini Kit (Qiagen, Tokyo, Japan) according to the standard protocol. Total RNA was isolated from whole peripheral blood samples using PaxGene Blood RNA Systems tubes (BD, Tokyo, Japan) according to the standard protocol. The RNA concentration and purity were measured with a NanoDrop–1000 (Thermo Fisher Scientific, Yokohama, Japan), and the 260/280 ratio was between 1.8 and 2.0. RNA (1.0 μg) in a 40-μl total reaction volume per sample was used to synthesize cDNA by using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA).
Genotyping analysis
SNP analysis of rs3764650 was performed with real-time quantitative PCR. A specific TaqMan probe for the SNP analysis was selected for NM_019112.3; assay ID: C__27478162_10 (TaqMan Assays, Applied Biosystems). The minor allele frequency of rs3764650 was based on the current International HapMap project database (http://hapmap.ncbi.nlm.nih.gov/index.html.en). The additional SNPs, rs429358 (NG_007084.2; assayID: C__3084793_20) and rs7412 (NG_007084.2; assay ID: C__904973_10, TaqMan Assays, Applied Biosystems), were used to determine the APOE isoform. PCR was performed with 1×TaqMan PCR Master Mix, 1×TaqMan SNP genotyping assay, 10 ng gDNA, and ultrapure water in a final reaction volume of 6 μl in each well of an optical plate. Allelic discrimination was determined using the StepOnePlus Real-Time PCR System and analyzed using StepOnePlussoftware.
Bisulfite conversion and pyrosequencing
The methylation status was analyzed in the region of ABCA7 that is considered the promotor region (Fig. 1). This region included five CpG sites that bind several transcription factors (TFs) investigated by JASPAR 2016 (http://jaspar.genereg.net). JASPAR is an open-access database storing curated, non-redundant TF binding profiles representing TF binding preferences as position frequency matrices for multiple species in six taxonomic groups [29, 30]. Therefore, we anticipated that methylation of these CpG sites will change the ABCA7 mRNA expression level. gDNA extracted from peripheral blood was treated with bisulfate using the EpiTect Plus DNA Bisulfite Kit (Qiagen, Hilden, Germany) and then amplified with PCR using the forward primer 5′-GGGTTTTCCTCAAAATCAGGGTAGCCACTA-3′ and the reverse primer 5′-[Biotin]-GACTAAACTAGAGGGAGCCTGGCCAGC-3′. PCR was performed in a 40-μl final volume per sample and included 2.5 μl gDNA, 0.2 μM each primer, AmpliTaq gold (Applied Biosystems), 10×PCR buffer with 15 mM MgCl2, and 2 mM dNTPs. Cycling conditions were as follows: denaturation at 94°C for 10 min, followed by 45 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 1 min, with a final extension of 10 min at 72°C. The PCR products were sequenced with PyroMark Q24 Advanced (Qiagen, Hilden, Germany) using the sequencing primer 5′-CCTCCTGCCATCCTCTGCAGAAGC-3′. Methylation rates at each CpG were quantified in duplicate using PyroMark Q24 Advanced Software (Qiagen, Tokyo, Japan).
Expression analysis
For mRNA expression analysis, quantitative reverse transcription-PCR was performed using the StepOnePlus Real-Time PCR System (Applied Biosystems). Specific TaqMan probes were Hs01105117_m1 for ABCA7 and Hs99999905_m1 for GAPDH. We used GAPDH as a housekeeping gene because previous studies including ours consistently identified GAPDH as a suitable housekeeping gene for peripheral blood gene expression analysis using the Paxgene blood RNA system [31–34]. The final volume of the reactions was 10 μl and included the TaqMan Universal Master Mix (Applied Biosystems). The expression levels were determined in triplicate. The ΔΔCt method was used to determine relative expression levels using StepOne software (Applied Biosystems) [35]. To correct for observational errors, PCR was performed in all plates with gDNA from the same control subjects in this study.
Statistical analysis
Statistical analysis was performed using SPSS Statistics version 22.0 (IBM Corp., Tokyo, Japan) and EZR [36]. For analysis of the effects of SNPs, all participants were divided into two groups according to the presence or absence of the rs3764650G allele in ABCA7 as in a previous study [18] and the APOE risk variant ɛ4. Additionally, all participants were divided into three groups according to the presence of the G allele of rs3764650 in a “dose-dependent” manner; T/T (dose = 0), G/T (dose = 1), and G/G (dose = 2). Linear regression analysis was performed between the ABCA7 mRNA expression level and the age, sex, rs3764650 allele, methylation rate of each CpG site, and APOE allele for AD patients and controls. The age and duration of illness in AD patients were compared between treated and not treated with AChEIs by the Student’s t-test or Mann-Whitney U test after the Shapiro-Wilk test. The ABCA7 mRNA expression level, mean methylation rate of all CpG sites were compared between AD and control subjects or AD patients treated or not treated with AChEIs using the Mann-Whitney U test with the Shapiro-Wilk test after Bonferroni correction. Sex differences and the distributions of ABCA7 alleles were compared with the χ2 test. The difference in mRNA expression levels between the two groups according to the rs3764650 allele was compared with the Mann-Whitney U test in AD and control subjects. Correlations among age, duration of illness, MMSE total score, NPI total score, rs3764650 allele, ABCA7 mRNA expression level, and methylation rate of each CpG site were analyzed with Spearman’s correlation test after Bonferroni correction. Relationships among CDR scores or MMSE recall, the ABCA7 mRNA expression level, and the methylation rate of each CpG site were analyzed with the Kruskal-Wallis test after the Steel-Dwass test.
RESULTS
Participant characteristics and genotyping
The demographic data of the participants are shown in Table 1. AD and control subjects did not differ in gender or age. Clinical characteristics and the results of psychological tests in AD subjects are also shown in Table 1. The distribution of the rs3764650 alleles is shown in Table 2. The ratio of each allele in AD was not significantly different from those in control subjects (number with the G allele, p = 0.548; presence or absence of the G allele, p = 0.603).
Methylation status
The sequence and the position of the five CpG sites are depicted in Fig. 1. Methylation rates at each CpG site are shown in Fig. 2. Methylation rates at all CpG sites were correlated with each other in AD subjects as well as control subjects (Fig. 1). Methylation rates at each CpG site were not associated with sex or age at blood sampling. The methylation rates at CpG3 were significantly correlated with the number of G alleles at rs3764650. We found no correlations between methylation rates at each CpG site and education, age of onset, or duration of illness, and no significant differences between methylation rates at each CpG site in AD subjects and those in control subjects (Table 3 and Fig. 2).
ABCA7 mRNA expression level
Linear regression analysis between the ABCA7 mRNA expression level and diagnosis, age, sex, rs3764650 allele, methylation rate of each CpG site, and APOE allele revealed that only diagnosis was significantly associated with the ABCA7 mRNA expression level (p = 0.003). The ABCA7 mRNA expression level was not associated with sex (AD subjects, p = 0.582; control subjects, p = 0.433) or age at blood sampling (AD, p = 0.349; control, p = 0.983). The ABCA7 mRNA expression level among the groups divided according to the rs3764650 allele was not different in AD or control subjects (AD, p = 0.727; control, p = 0.894). The ABCA7 mRNA expression level in AD subjects was significantly higher than that in controls (Fig. 3, p = 0.017). Moreover, we found a significant negative correlation between the duration of illness and the ABCA7 mRNA expression level (r = –0.433, p = 0.002, Table 4 and Fig. 4). Even excluding 5 AD patients used in the period beyond 10 years of illness, the ABCA7 mRNA expression level was negatively correlated with duration of illness (Spearman’s correlation coefficients, r = –0.391, p = 0.007). The ABCA7 mRNA expression level in AD subjects treated with AChEIs was significantly lower than in those not treated with AChEIs (Supplementary Figure 1, p = 0.013, *p≤0.025 after Bonferroni correction was considered significant), and the duration of illness in AD subjects treated with AChEIs was longer than in those not treated with AChEIs (p = 0.015).
Relationships among the ABCA7 mRNA expression level, methylation rate, rs3764650 allele, and scores on the MMSE, CDR, and NPI
Correlations between the ABCA7 mRNA expression level and each parameter including psychological tests in AD subjects are shown in Table 4. The ABCA7 mRNA expression level was significantly associated with the ADAS total score and CDR score (ADAS total score, p = 0.005; CDR score, Kruskal-Wallis test, p = 0.032, Steel-Dwass, score 0.5 : 2 points, p = 0.01, Table 4 and Fig. 5, Supplementary Figure 2). Although the ABCA7 mRNA expression level was not correlated with the MMSE or NPI score, the recall score in the MMSE subscale was significantly associated with the ABCA7 mRNA expression level (Kruskal-Wallis test, p = 0.0068, Steel-Dwass, score 0 : 2 points, p = 0.049, Table 4 and Supplementary Figure 3). Other than CpG3, methylation rates in AD subjects were not correlated with the ABCA7 mRNA expression level, rs3764650, MMSE total score, CDR score, or NPI total score (Table 3). We found no allele effects of rs3764650 on cognition assessed with MMSE, ADAS, or CDR (data not shown).
DISCUSSION
Our study revealed three major findings. First, the frequencies of the rs3764650 alleles in AD subjects were not different from those in control subjects in our study, although this SNP was significantly associated with AD in a genome-wide association study [4]. Another study using AD brains reported that the rs3764650 T allele is associated with increased ABCA7 expression and decreased AD risk [18]. Although a meta-analysis including Japanese individuals identified a significant association between AD and rs3764650 [17], a study with only Japanese patients reported a negligible association [16]. However, we performed power calculations with Fisher’s exact test using EZR, and the sample size was not sufficient to detect differences in the distribution of alleles with rs3764650 (n = 2184). Therefore, we cannot conclude that rs3764650 is not associated with AD. We found no allele effects on ABCA7 mRNA expression in peripheral blood or cognition assessed with MMSE, ADAS, or CDR. Our results are consistent with previous reports that showed no association between the MMSE score and rs3764650 [16, 37]. Differences in ethnicity, sample size, and experimental setting (peripheral blood vs. brain tissue) may influence the results.
Second, the ABCA7 mRNA expression level in peripheral blood from AD patients was significantly higher than that in control subjects. ABCA7 is expressed throughout the mouse brain with particularly high expression in CA1 hippocampal neurons [38], as well as in the human brain [39, 40]. Additional analysis in human brain cells showed that the ABCA7 mRNA expression level is the highest in microglia [10]. ABCA7 likely plays roles in transporting phospholipids and enhancing phagocytosis by macrophages [9–13]. Peripheral lymphocyte infiltration was observed in capillaries of AD brain tissue [41]. Recent study showed that in a mouse model of AD, Aβ-specific lymphocyte occurred in the brain after a single Aβ immunization provided that small amounts of IFN-γ were expressed in the brain, and induced microglia activation and efficient clearance of Aβ [42]. These findings indicated that peripheral lymphocyte could be influenced on the change of brain in AD. In another study, the ABCA7 mRNA expression level in peripheral blood was negatively correlated with hippocampal atrophy as assessed with magnetic resonance imaging in 50 normal individuals and 98 with mild cognitive impairment [37]. Therefore, the ABCA7 mRNA expression level may be associated with cognitive impairment in AD. The ABCA7 mRNA level in AD subjects treated with AChEIs was significantly lower than in those not treated with AChEIs in this study, and the duration of illness in AD subjects treated with AChEIs was significantly longer than in those not treated with AChEIs. Because the ABCA7 mRNA expression level in AD subjects showed a significant negative correlation with the duration of illness, we conclude that the ABCA7 mRNA expression level was more influenced on longer duration of illness rather than on whether medicated with AChEIs or not. Another study reported that the ABCA7 expression level in the parietal lobes of European Americans with AD was positively associated with the CDR [43]. The difference between this finding and ours may be influenced by the timing of brain or blood sampling. In the future, we should analyze ABCA7 expression levels in human brain and peripheral blood during the same period. Additionally, many neurodegenerative diseases including AD share chronic immune activation as a common feature, although these diseases are probably triggered by many different initiating events in the early stages [44]. Higher ABCA7 expression levels in peripheral blood may reflect neuroprotective effects in early stages of AD and a decrease in ABCA7 mRNA may reflect disease progression and cognitive decline.
Third, a previous study reported that three CpG sites were significantly associated with pathological AD diagnosis. After adjusting for multiple testing, cg02308560 and cg04587220 (30 base pairs apart) were observed in a polycomb-repressed region in HMHA1 proximal to the 3’-untranslated region of ABCA7. However, this study used human brain, and these CpG sites were far (about 30,000 base pairs) from the TSS of ABCA7. We set the target CpG sites within 500 base pairs of the TSS of ABCA7. Although this region has five CpG sites and is predicted to bind several major TFs [29, 30], we found no association between methylation rates of ABCA7 and AD or the ABCA7 expression level. Although the methylation rate at only CpG3 was significantly correlated with the number of G alleles at rs3764650, we found no correlations between the methylation rate at each CpG site and any other clinical parameters in AD or control subjects. Because this is the first study showing a relationship between AD and the methylation rates of ABCA7 in peripheral blood, further studies with other tissues, such as cerebral cortices, are warranted to confirm the finding.
Our study has several limitations. The first limitation is the small sample size. Second, this study was cross-sectional. We recruited AD patients in various stages and found an inverse correlation between ABCA7 mRNA expression and duration of illness. We should examine ABCA7 mRNA prospectively from the time of mild cognitive impairment through late stages of AD to confirm the stage-dependent expression changes. Third, we selected target CpG sites that may theoretically bind major TFs, and a verification study in vitro such as a luciferase assay is needed to confirm that these target CpG sites are actually bound by TFs and regulate gene expression. Finally, although there was no evidence of cognitive decline from a previous level or mental retardation in healthy controls, no education data for healthy controls is the limitation.
In summary, our results suggest that a higher ABCA7 mRNA expression level in peripheral blood may be a marker for neuroprotective effects against AD pathology, and a decrease in ABCA7 mRNA may reflect disease progression and cognitive decline regardless of rs3764650 and the methylation rates.
