Neuron-specific methylome analysis reveals epigenetic regulation and tau-related dysfunction of BRCA1 in Alzheimer's disease.

Alzheimer's disease (AD) is a chronic neurodegenerative disease characterized by pathology of accumulated amyloid β (Aβ) and phosphorylated tau proteins in the brain. Postmortem degradation and cellular complexity within the brain have limited approaches to molecularly define the causal relationship between pathological features and neuronal dysfunction in AD. To overcome these limitations, we analyzed the neuron-specific DNA methylome of postmortem brain samples from AD patients, which allowed differentially hypomethylated region of the BRCA1 promoter to be identified. Expression of BRCA1 was significantly up-regulated in AD brains, consistent with its hypomethylation. BRCA1 protein levels were also elevated in response to DNA damage induced by Aβ. BRCA1 became mislocalized to the cytoplasm and highly insoluble in a tau-dependent manner, resulting in DNA fragmentation in both in vitro cellular and in vivo mouse models. BRCA1 dysfunction under Aβ burden is consistent with concomitant deterioration of genomic integrity and synaptic plasticity. The Brca1 promoter region of AD model mice brain was similarly hypomethylated, indicating an epigenetic mechanism underlying BRCA1 regulation in AD. Our results suggest deterioration of DNA integrity as a central contributing factor in AD pathogenesis. Moreover, these data demonstrate the technical feasibility of using neuron-specific DNA methylome analysis to facilitate discovery of etiological candidates in sporadic neurodegenerative diseases.

Alzheimer’s disease (AD) is the most prevalent form of sporadic neurodegeneration. Although various studies have indicated that amyloid β (Aβ) and phosphorylated tau (pTau) are deeply involved in its pathogenesis, our current knowledge about these molecules is insufficient for the introduction of drugs that may ameliorate or cure the disease (1). There are several different strategies commonly used to investigate the pathogenesis of sporadic neurodegenerative diseases. The first approach has been to analyze genomic data of patients (i.e., SNPs) and compare with existing disease phenotypic signatures. This is a rather straightforward method since causal relation is often clear and also suitable when an underlying factor produces a particularly large effect. The discovery of the APOE ε4 allele as a risk factor is an important example of successful implementation of this approach in AD research (2). Despite this success, the methodology has not led to the discovery of additional factors largely involved in disease pathogenesis, since the effects of other genetic risk loci are often relatively small (3–5). A second approach has been to analyze postmortem brain samples from AD patients biochemically or histopathologically. While a more conventional method compared with genomic analysis, difficulties remain in extracting useful information from such samples. Unlike classic pathological analysis, high-throughput screening by transcriptome analysis have been proved to be powerful providing us with insights into the etiology of AD (6–9). However, many of these studies have analyzed bulk messenger RNA (mRNA) extracted from unprocessed brain samples, which are a mixture of mRNA originating from different cell types, including neurons, astrocytes, oligodendrocytes, microglia, and blood cells. Laser microdissection has been used to exclusively extract neuronal mRNA to overcome this drawback; however, the relatively small number of cells that can be analyzed (a maximum of ∼1 × 103 for each sample) and manually selected dissection areas can lead to problems of sampling bias. In addition, the vulnerability of mRNA molecules to postmortem degradation may also bias results (10, 11).

DNA methylation is an epigenetic gene regulation mark that occurs at the cytosine residue of CG dinucleotide sequences. In the genome, regions rich in CG doublet sequences, known as CpG islands, tend to be located at gene promoter regions. Methylation occurring at promoter regions inhibits downstream gene expression through interaction with methyl-CpG–binding proteins. In contrast, hypomethylated promoter CpG islands correspond to positive downstream gene expression (12). Compared with histone modification, another epigenetic mark that occurs in a more rapid manner, DNA methylation is thought to reflect constitutive regulation of gene expression (13). Although methylome does not provide direct information about gene expression, it serves as a rough representation of the entire transcriptome. Analysis of DNA methylation offers three significant benefits: (i) stability of DNA methylation in postmortem brains (14, 15), (ii) high-throughput isolation of up to 1 × 106–7 neuronal from nonneuronal nuclei by fluorescence-activated cell sorting (FACS) to eliminate any sampling bias (16), and (iii) implementation of high-throughput analysis using a DNA microarray system (17, 18).

Here, we present an analysis of the neuron-specific DNA methylome from human AD brains. Using this strategy, we identified three candidate molecules potentially related to AD but whose link to the disease has not been well studied. One of the molecules that emerged from the screen was BRCA1. We investigate the role of BRCA1 in human AD brains and model systems and demonstrate the importance of its interaction with the tau protein in the pathogenesis of AD.

Results

Genome-Wide Neuron-Specific DNA Methylome Analysis Reveals Dysregulated Genes in the AD Brain.

We examined inferior temporal gyrus samples from 30 age-matched normal control (NC) and 30 AD patients (Table 1). Neuronal and nonneuronal nuclei were separated by FACS using an anti-NeuN antibody (). Starting material of 100–400 mg of brain tissue resulted in yields of ∼1 × 106–7 nuclei and 600–2,000 ng of genomic DNA. Extracted DNA was bisulfite converted and subjected to genome-wide DNA methylation analysis using an Infinium 450k methylation array. The data were analyzed using a filtering procedure described in . The Infinium 450k contains 485,512 probes in total. Data with detection P values ≥ 0.01 were excluded for quality control purposes, along with probes on the sex chromosomes, to obtain 414,222 probes. Principal-component analysis of the top 1,000 most variable probes showed no evidence of separation between NC and AD, suggesting that there was little difference between the two methylomes (). Subsequently, probes that met the following criteria of statistical difference between the two groups were further analyzed: (i) mean β-value difference > 0.05 (ii), value of P < 0.05 by two one-sided t tests. Using this procedure, the number of probes of interest was reduced to 278 (). Since transcriptional regulation through DNA methylation requires regional modification rather than a single isolated change, we attempted to further refine the results by testing if the methylation occurred at a regional level. We defined differentially methylated regions (DMRs) as those including three or more consecutive differentially methylated probes. By analyzing the 278 probes in this manner, we found eight DMRs consisting of 36 statistically significant probes ().

Table 1.

Clinical demographics of postmortem NC and AD brains used in neuron-specific methylome analysis

Demographic variable	NC	AD	Statistics
No. of samples	30	30
Male ratio, %	50	50
Age at death	76.7 ± 7.4	79.4 ± 7.4	P = 0.16
APOE ε4 ratio	26.7	66.7	P = 0.004

P values are calculated by t test or χ2 test.

Three of the detected DMRs were located at the promoter CpG islands of the BRCA1, ZNF714, and AURKC genes. Another DMR was found at the promoter CpG island of LOC441666; however, this locus is documented as a pseudogene. The remaining four DMRs at DUSP5P1, PCDHB7, SSR1, and ERICH1 were located at intergenic regions or the gene body, where the biological significance was considered less significant and so were not pursued further. Our analysis focused on the DMRs at BRCA1, ZNF714, and AURKC. The DMRs at BRCA1 and AURKC were hypomethylated in AD samples, while that of ZNF714 was hypermethylated (Fig. 1 and ). Methylation levels of the DMRs at each of these three genes were revealed by correlation analysis and found to be tightly regulated, implying that methylation levels correlate with downstream gene regulation ().

Fig. 1.

Details of DMR associated with BRCA1. The uppermost panel shows genomic structures from the University of California, Santa Cruz, genome browser including chromosomal ideogram, GC content (vertical bars), RefSeq gene map (blue lines and bars), and CpG islands (green bars). The second panel depicts β values of each sample at each CpG site. Dots are β values from each sample, and solid lines are mean β values of each group. Light green areas correspond to CpG islands. CG identifiers listed are the sites where methylation levels were validated using an additional pyrosequencing assay. The last two panels show β-value differences (Upper, βNC − βAD) and significance levels on a −log(P) scale at each CpG site. Orange dots represent differentially methylated probes (DMPs) discovered by Infinium screening.

Next, we analyzed the CpG methylation levels of the DMRs at BRCA1, AURKC, and ZNF714 by pyrosequencing to confirm the results of the aforementioned screening. Pyrosequencing and methylation array analysis were found to be consistent (). Methylation levels of the differentially methylated probes (DMPs) associated with these three genes did not correlate with age at death of the subjects (). The methylation levels of the ZNF714 and AURKC DMPs did demonstrate significant differences according to sex, whereas the BRCA1 DMPs showed no such difference (). Interestingly, methylation levels of the BRCA1 and AURKC DMRs were significantly correlated with the number of APOE ε4 alleles, but only in AD neuronal samples (). This suggests that the amount of Aβ deposition and subsequent neurodegeneration could be related to their altered methylation.

To test whether these results were specific to AD, we analyzed neuronal DNA from 16 inferior temporal gyrus samples of dementia with Lewy bodies (DLB) (). For BRCA1 and AURKC, the differences in methylation levels between the NC and DLB groups were not statistically significant, while for ZNF714 the changes observed for AD were also seen in DLB (). Next, we monitored differences in methylation levels of nonneuronal samples from the inferior temporal gyrus or neuronal samples from the cerebellum of the same subjects. Surprisingly, the differences remained at several CpGs of nonneuronal DNA (), and even in DNA from cerebellar neuronal cells (predominantly granular cells as Purkinje cells are NeuN negative) (). Thus, the aberrant DNA methylation pattern appeared to be a widespread change within AD brains. Based on this, we decided to perform transcriptional analysis from postmortem brains using bulk mRNA. For this experiment, mRNA extracted from the temporal lobe in previously untested samples was used to validate the results of the methylome analysis (). In AD brains, expression of BRCA1 and AURKC were clearly up-regulated; however, ZNF714 showed no statistically significant difference between the two groups (). This finding was mostly in agreement with the expression changes postulated from the methylation changes. Together, we hypothesized that the difference in promoter methylation levels of BRCA1 and AURKC would have a significant influence on expression of these genes, particularly in AD brains.

BRCA1 Is Mislocalized and Insoluble in the AD Brain.

We next examined the protein expression levels of BRCA1 and AURKC in postmortem brains, starting with NC and advanced AD subjects at Braak stages 5–6 (). Immunohistochemical analysis showed remarkable staining of BRCA1 in the cytoplasm of neuronal cells in only the AD brains. Staining was prominent in the hippocampal CA1 region and entorhinal cortex, less so in the parietal lobe and anterior cingulate gyrus. Localization of BRCA1 in the occipital lobe or the cerebellum was nuclear, and no significant cytoplasmic signal was detected (Fig. 2 ). Furthermore, BRCA1 was localized to the cytoplasm of neuronal cells that well colocalized with pTau (Fig. 2). The result was far less remarkable for AURKC, which showed a slightly increased staining in the hippocampal CA1 region and entorhinal cortex (). As changes in AURKC staining were less marked than those of BRCA1, further analysis was focused on BRCA1.

Fig. 2.

BRCA1 is mislocalized at the cytoplasm of AD brains, and it occurs at the insoluble fraction. (A) Immunohistochemical images of various regions from advanced-stage AD or NC by anti-BRCA1 antibody. Representative immunohistochemical images are shown from a total of n = 6 (hippocampus, entorhinal cortex, anterior cingulate gyrus, and parietal lobe) or n = 4 (occipital lobe and cerebellum) each. (B) Statistical analysis of A. Error bars represent means ± SEM. Significance was determined using one-way ANOVA followed by post hoc Holm–Sidak method. Mean ± SEM: NC: 3.72 ± 1.37/mm2 in the hippocampus, 1.35 ± 0.52/mm2 in the entorhinal cortex, 0/mm2 in the anterior cingulate gyrus, the parietal lobe, the occipital lobe, and the cerebellum. AD: 67.32 ± 21.27/mm2 in the hippocampus, 39.69 ± 8.05/mm2 in the entorhinal cortex, 15.33 ± 1.70/mm2 in the anterior cingulate gyrus, 8.46 ± 2.31/mm2 in the parietal lobe, 0.68 ± 0.48/mm2 in the occipital lobe, and 0/mm2 in the cerebellum. (C) Immunofluorescence images of advanced stage AD brains colabeled with DAPI, anti-BRCA1, and anti-pTau antibodies. Representative images from n = 4. (D) Immunohistochemical images of inferior temporal gyrus samples from NC or advanced-stage AD brains by anti–γ-H2ax antibody. Representative images are shown from a total of n = 3 samples of AD and NC. (E) Serial solubilization of postmortem brain samples using various detergents. The samples were treated with indicated detergents in a serial manner, each time saving the centrifuged supernatants as soluble fractions for each detergent. Three AD and NC samples were tested. (F) Immuno-EM of purified PHF by anti-BRCA1 antibody. *P < 0.05.

The distribution pattern of BRCA1 was similar to that of advanced AD pathology, consistent with a model where the AD pathology itself could be driving expression (Fig. 2 ). To test this hypothesis, the immunostaining pattern of BRCA1 was compared with those of pathological hallmarks of AD (e.g., deposition of Aβ or pTau) in subjects harboring earlier stage AD pathology. At Braak stage 3, Aβ- and pTau-positive pathology was already present in both the hippocampal CA1 region and the entorhinal cortex, while cytoplasmic BRCA1 was detected mainly in the hippocampal CA1 region (), especially in the neuronal cells having prominent pTau deposition (). This indicates that accumulation of cytoplasmic BRCA1 may be a consequence of pTau deposition, which well explained the observation that the region without tau pathology, that is, the occipital lobe and the cerebellum, is mostly free from cytoplasmic accumulation of BRCA1 despite decreased methylation.

BRCA1 is a nuclear protein with a central role in DNA repair. Among various DNA repair related proteins, only γ-H2ax, a well-established marker for initiation of DNA double-strand break (DSB) repair, was prominent in AD brains (Fig. 2 and ). This result suggested that Aβ induced DNA DSBs, of which BRCA1 specializes in their repair and possible failure of the downstream process. Considering that aberrant cytoplasmic localization of BRCA1 could be indicative of its dysfunction, we tested its solubility in AD brains. Extraction via a series of detergents revealed a substantial shift in BRCA1 solubility in AD brains, suggesting that the majority of BRCA1 in the AD brain is dysfunctional (Fig. 2). For further assessment, we performed immuno-electron microscopy (EM) of the paired helical filaments (PHF) fraction to find that PHFs were positive with anti-BRCA1 antibodies (Fig. 2). Immunoprecipitation of soluble tau from human AD brain lysates showed no association of BRCA1, indicating no direct binding of BRCA1 to soluble tau (). Taken together, we concluded that insolubilized tau as PHF sequesters BRCA1 by coaggregation.

Aβ Burden Induces DNA DSBs and Up-Regulation of BRCA1 Expression in Vitro and in Vivo.

We next focused on the mechanism of BRCA1 up-regulation. Since the BRCA1 DMR was also hypomethylated in nonneuronal cells that are essentially free of tau pathology, we reasoned that its up-regulation might be more closely related to Aβ than tau. To test whether Aβ itself was sufficient to up-regulate BRCA1, we analyzed N2a swe.10 cell line that expresses human amyloid precursor protein (APP) (KM670/671NL). This cell line produces significant levels of intracellular APP (Fig. 3 and ) and secretes extracellular Aβ40 and Aβ42 into the medium (). N2a swe.10 cells expressed more soluble BRCA1 compared with control cells (mean ± SEM: 1.00 ± 0.06 in N2a cells vs. 1.44 ± 0.15 in N2a swe.10 cells) (Fig. 3 ). Full-length BRCA1 was found in the soluble fractions of both control and N2a swe.10 cells (Fig. 3). To confirm that BRCA1 expression was caused by Aβ and not by APP overexpression, we treated N2a swe.10 cells with a potent γ-secretase inhibitor, compound E. The treatment significantly decreased BRCA1 expression as well as Aβ, without having any effect on APP expression in N2a swe.10 cells, while BRCA1 expression was not affected in N2a cells (Fig. 3 and ). To eliminate any off-target effect of compound E, we treated N2a cells with recombinant Aβ40 or Aβ42 and found that both of them increased BRCA1 but Aβ40 required much higher concentration compared with Aβ42, suggesting that their toxicity is driving BRCA1 expression (Fig. 3). Immunostaining revealed multiple positive nuclear γ-H2ax signals only in N2a swe.10 cells, consistent with the activation of DNA repair machinery (Fig. 3 ). To exclude the possibility that overexpressed APP may affect formation of nuclear γ-H2ax signals and BRCA1 up-regulation, N2a cells were cocultured with N2a swe.10 or N2a cells (Fig. 3). Aβ secreted from donor N2a swe.10 cells was able to induce up-regulation of BRCA1 and multiple positive nuclear γ-H2ax signals were detected in recipient cells (Fig. 3 ). We analyzed the methylation level of Brca1 CpG island, but there was no difference between the two cell lines (). To evaluate the amount of DNA damage in N2a swe.10 cells, single-cell gel electrophoresis (comet assay) was performed. Control experiment using etoposide resulted in significant DNA fragmentation (Fig. 4 ). We found no significant DNA fragmentation in N2a swe.10 cells compared with controls (Fig. 4 ), indicating that sufficient DNA repair continued to occur despite activated DSBs in these cells. Knockdown of BRCA1 by lentiviral shRNA impaired DNA repair to increase DNA damage in N2a swe.10 cells (Fig. 4 ). These data indicated that the presence of Aβ triggers DNA DSBs, with the recruitment of BRCA1 sufficient to neutralize Aβ toxicity.

Fig. 3.

Aβ confers BRCA1 overexpression. (A) Western blot of N2a and swe.10 cells using anti-mouse BRCA1 and anti-APP antibodies. Actin was detected as a loading control. (B) Quantitative measurement of the relative amount of APP and BRCA1 from experiments in A. n = 4. Relative quantity against actin expression level was normalized to N2a = 1.0. Mean ± SEM: APP: 1.00 ± 0.05 in N2a cells vs. 2.19 ± 0.10 in N2a swe.10 cells; BRCA1: 1.00 ± 0.06 in N2a cells vs. 1.44 ± 0.15 in N2a swe.10 cells. (C) Serial fractionation of N2a and swe.10 cells using various detergents. (D) Concentrations of Aβ40 and Aβ42 in the culture media supernatant of N2a swe.10 cells treated with dimethyl sulfoxide (DMSO) or 25 nM compound E were measured by ELISA. n = 6. Mean ± SEM: Aβ40: 23.00 ± 1.28 pM in DMSO vs. 5.08 ± 1.00 pM in compound E; Aβ42: 4.16 ± 0.29 pM in DMSO vs. 1.45 ± 0.04 pM in compound E. (E) Western blot of N2a swe.10 cells treated with DMSO or compound E by anti-mouse BRCA1 and anti-APP antibodies. (F) Quantitative measurement of the relative amount of BRCA1 and APP from E. n = 6. Relative expression level was normalized to DMSO = 1.0. Mean ± SEM: BRCA1: 1.00 ± 0.08 in DMSO vs. 0.41 ± 0.02 in compound E; APP: 1.00 ± 0.03 in DMSO vs. 0.94 ± 0.02 in compound E. (G) Effect of recombinant Aβ40 and Aβ42 on BRCA1 expression. Western blot of N2a cells by anti-mouse BRCA1 and anti-actin antibodies. Quantitative measurement of the relative amount of BRCA1 is shown below (n = 6 independent wells). Each relative expression level of BRCA1 was normalized to DMSO = 1. One-way ANOVA [Aβ40: F(2,15) = 11.63, P = 0.0009; Aβ42: F(3,20) = 7.981, P = 0.0011] with post hoc Turkey method. Mean ± SEM: Aβ40: 1.00 ± 0.09 in DMSO, 1.10 ± 0.11 in 100 nM, 1.80 ± 0.17 in 1 μM; Aβ42: 1.00 ± 0.10 in DMSO, 0.97 ± 0.12 in 100 pM, 1.50 ± 0.08 in 1 nM, 1.62 ± 0.16 in 10 nM. (H) Immunofluorescence images of N2a and swe.10 cells stained with DAPI, and anti–γ-H2ax antibody. Insets show single nuclei at high magnification. (I) Quantitative analysis of the number of cells with nuclear γ-H2ax foci in H. n = 9 visual fields (3 visual fields from 3 experiments). Mean ± SEM: 0.34 ± 0.20% in N2a cells and 3.02 ± 0.25% in N2a swe.10 cells. (J) Diagram of coculture system. Recipient cells were used for biochemical and immunohistochemical analysis. (K) Western blot of recipient N2a cells by anti-mouse BRCA1 and anti-actin antibodies. Quantitative measurement of the relative amount of BRCA1 is shown below (n = 4 independent wells). Each relative expression level was normalized to N2a-donor culture. Mean ± SEM: 1.00 ± 0.11 in N2a cells and 1.44 ± 0.14 in N2a swe.10 cells. (L) Quantitative analysis of the number of cells with nuclear γ-H2ax foci in the recipient N2a cells. n = 9 visual fields (3 visual fields from 3 experiments). Mean ± SEM: 0.31 ± 0.20% in N2a cells and 2.81 ± 0.35% in N2a swe.10 cells. *P < 0.05 and ****P < 0.0001. ns, not significant.

Fig. 4.

Functional relevance of BRCA1 dysfunction in an in vitro and in vivo neuronal model with Aβ burden. (A) Comet assay images of N2a treated with DMSO or 50 μM etoposide for 6 h. (B) Quantification of tail DNA%. Closed circles represent N2a cells treated with DMSO, and open circles, etoposide. n = 39 cells (DMSO) and 31 cells (etoposide). Mean ± SEM: 5.07 ± 0.21% in DMSO and 6.60 ± 0.73% in etoposide. (C) Comet assay images of N2a and N2a swe.10 cells. Representative images are shown. Insets are a high magnification of single nuclei. (D) Quantification of tail DNA%. Closed circles represent N2a cells, and open circles, swe.10. n = 150 cells. Gray bars represent means. Mean ± SEM: 5.72 ± 0.34% in N2a cells and 6.68 ± 0.42% in N2a swe.10 cells. Biological replicates are shown in . (E) Representative comet assay images are shown from a total of n = 80 N2a and swe.10 cells after lentiviral transduction of BRCA1 shRNA. (F) Statistical significance of BRCA1 knockdown on DNA fragmentation was determined using the Tukey method. Bars represent means. Mean ± SEM: N2a: 3.45 ± 0.24% in shRNA control, 3.59 ± 0.27% in shRNA #1, 4.40 ± 0.34% in shRNA #2; N2a swe.10: 3.67 ± 0.35% in shRNA control, 5.59 ± 1.37% in shRNA #1, 18.06 ± 3.10% in shRNA #2. Biological replicates are shown in . (G) Representative images of the differentiated N2a and swe.10 cells treated with shRNAs are shown (n = 8 visual fields). (H) Quantification of cells with neurite-like process. Statistical significance was calculated by the Tukey method. Mean ± SEM: N2a: 64.86 ± 3.65% in shRNA control, 59.64 ± 4.33% in shRNA #1, 61.62 ± 4.36% in shRNA #2; N2a swe.10: 68.98 ± 6.62% in shRNA control, 45.46 ± 5.81% in shRNA #1, 24.61 ± 4.98% in shRNA #2. (I) Representative images of primary neuronal cultures of cortical tissues from 3×Tg mice. Total length of neurites (J), maximal length of neurites (K), and neurite length distribution (L) were analyzed. Control (n = 12 neurons), shRNA #1 (n = 11 neurons), and shRNA #2 (n = 10 neurons). Statistical significance was determined by the Tukey method. Boxes extend from the 25th to 75th percentiles, and the lines in the boxes represent the median. The whiskers show the minimum and maximum values. (M) Knockdown of BRCA1 in vivo. Lentivirus expressing shRNA against BRCA1 was stereotactically injected into the dentate gyrus (DG) of APP/PS1 mice at 3 mo of age, and mice were killed 3 wk after the surgery. (N) Immunofluorescence images of the neuronal cell in DG stained with DAPI, and anti–γ-H2ax antibody. (O) Spine density of the neuronal cells in DG. n = 8 (control) and 7 (shRNA #1) dendrites. Representative images are shown. (P) Statistical analysis of O. Mean ± SEM: 0.76 ± 0.02/μm in shRNA control and 0.61 ± 0.03/μm in shRNA #1. (Q) Nucleus from inferior temporal gyrus of four NC and four AD brains were subjected to comet assay. Representative images are shown. (R) Statistical analysis of Q. Boxes extend from the 25th to 75th percentiles, and the lines in the boxes represent the median. The whiskers show the minimum and maximum values. n = 70–100 nuclei from each sample. Statistical significance was determined using two-way ANOVA [NC vs. AD: F(1,707) = 100.6, P < 0.0001; sample: F(3,707) = 52.30, P < 0.0001]. Mean ± SEM: NC: 3.54 ± 0.24%, 3.36 ± 0.23%, 3.13 ± 0.21%, and 2.58 ± 0.15%; AD: 7.45 ± 0.34%, 5.46 ± 0.19%, 3.64 ± 0.22%, and 2.80 ± 0.24%. **P < 0.01 and ****P < 0.0001.

To assess the functional relevance of BRCA1 function against Aβ in neuronal cells, we analyzed their morphological changes upon BRCA1 knockdown. BRCA1 knockdown suppressed formation of neurite-like processes in N2a swe.10 cells upon differentiation but not in N2a cells (Fig. 4 ). Furthermore, in primary cultures of neuronal cells from 3×Tg-AD mice harboring APP KM670/671NL and MAPT P301L transgenes, along with knocked-in PSEN1 M146V, which is a well-established model for both Aβ and pTau deposition (19), BRCA1 knockdown did not have significant effect on the global morphology (Fig. 4 ), but reduced the number of short neurites (mean ± SEM: for dendrites with 0–1 AU in length, 50.50 ± 7.16 in shRNA control, 39.18 ± 5.38 in shRNA #1, 22.60 ± 4.07 in shRNA #2) (Fig. 4). We further assessed in vivo function of BRCA1 in APP/PS1 mice that expresses APP KM670/671NL and PSEN1dE9 (20, 21). Knocking down BRCA1 using stereotactic lentivirus injection into the dentate gyrus (DG) induced γ-H2ax foci and reduced spine density (Fig. 4 ).

We next used APP/PS1 mice to assess the effect of Aβ toxicity on BRCA1 expression and DNA damage in vivo. These mice exhibit Aβ pathology, but tau pathology is absent (). Although γ-H2ax could be detected from the earliest stage (), suggesting DNA DSBs were triggered in the hippocampal neurons, comet assays of hippocampal cells in the mice showed no significant fragmentation of genomic DNA (). Immunohistochemistry showed no cytoplasmic BRCA1 in the hippocampal regions in these mice, and biochemical analysis demonstrated increased soluble BRCA1 [mean ± SEM: 1.00 ± 0.10 in wild type (WT) vs. 1.59 ± 0.09 in APP/PS1 in PBS fraction] (). Data from these mice indicated that DNA DSBs were induced by Aβ but sufficiently repaired by BRCA1. To see whether DNA damage was actually happening in human AD brain, temporal lobe samples were subjected to comet assay, which showed significant DNA damage in AD brains (Fig. 4 ).

If this is the case, what makes BRCA1 insoluble and dysfunctional in human AD brains? We hypothesized that aggregated tau could have a role in this phenomenon, since there were almost no glial cells positive for cytoplasmic BRCA1, despite a widespread methylation change in AD brains (Fig. 2 and ). To evaluate its role in BRCA1 inactivation, we first tested whether BRCA1 was insolubilized in in vitro seed-dependent tau aggregation model (22). Upon aggregation of P301L tau by adding in vitro-generated seed, endogenous BRCA1 was shifted to insoluble fraction (Fig. 5 ) and cytoplasmic tau aggregates colocalized with BRCA1, precluding its nuclear localization (Fig. 5). To investigate whether this was also true in in vivo model, we histopathologically and biochemically analyzed 3×Tg-AD mice. As previously reported, Aβ and pTau accumulate in an age-dependent manner in hippocampal CA1 neurons (19, 23) (Fig. 5 ). In these mice, positive γ-H2ax signal preceded the appearance of cytoplasmic BRCA1 and the number of positive cells increased in an age-dependent manner (Fig. 5). Cytoplasmic BRCA1 was also detected from the age of 6 mo, gradually increasing with age up to 12 mo (Fig. 5). Together with hippocampal tau deposition in the 3×Tg-AD mice brain, we postulated that the cytoplasmic BRCA1 observed was the result of neuronal tau deposition. Consistent with this hypothesis, cortical regions of 3×Tg-AD mice show only Aβ pathology without tau deposition (23), similar to the APP/PS1 mice. These same cortical regions also exhibited increased levels of soluble BRCA1 (Fig. 5). In contrast, the hippocampal region of 3×Tg-AD mice had increased levels of insoluble BRCA1 along with increased insoluble pTau at the age of 12 mo (Fig. 5).

Fig. 5.

Cytoplasmic tau is important for mislocalization and insolubilization of BRCA1. (A) SH-SY5Y cells transfected with wild type (WT) or P301L mutant 4-repeat tau were treated with in vitro-aggregated tau seeds [either full length (FL) or C-terminal fragment (251)]. After fractionation, supernatant (sup) or pellet (ppt) was subjected to Western blot. (B) Quantification of BRCA1 in the pellet fraction. Mean ± SEM: 1.00 ± 0.09 in WT, 1.08 ± 0.08 in P301L, 0.87 ± 0.11 in WT + FL, 0.95 ± 0.09 in WT + 251, and 1.74 ± 0.15 in P301L + FL. (C) Immunocytochemical images of seed-dependent tau aggregation colabeled with DAPI, pTau, and BRCA1 antibodies. (D) Representative images of Aβ staining of 3×Tg-AD mice hippocampal region. (E–G) Representative immunohistochemical images of 3×Tg-AD mice hippocampus by anti-pTau (E), anti-γ-H2ax (F), and BRCA1 (G) staining. Samples are from 3×Tg-AD mice at 3, 6, 9, and 12 mo of age. n = 10–15 visual fields from six to seven mice for 3×Tg-AD mice. Quantitation of positive cells for each protein are shown at the Right side of the images. Statistical significance was determined by the Tukey method. Mean ± SEM: pTau: 33.26 ± 2.57/mm in 3 mo, 37.30 ± 1.99/mm in 6 mo, 44.56 ± 3.19/mm in 9 mo, and 55.12 ± 2.69/mm in 12 mo; γ-H2ax: 1.28 ± 0.22/mm in 3 mo, 3.23 ± 0.90/mm in 6 mo, 4.14 ± 0.66/mm in 9 mo, and 7.97 ± 0.90/mm in 12 mo; BRCA1: 0.56 ± 0.21/mm in 3 mo, 1.07 ± 0.24/mm in 6 mo, 3.41 ± 0.33/mm in 9 mo, and 6.94 ± 0.77/mm in 12 mo. (H and I) Expression of BRCA1 in detergent fractionated the cortical (H) and hippocampal (I) regions of 3×Tg-AD mice at 12 mo of age. n = 2 mice at 12 mo of age. Representative blot was shown. (J) Microscopic images of comet assays. Representative images are shown from a total of n = 136 cells from WT and 3×Tg-AD animals. Insets are magnified nuclei. (K) Quantification of tail DNA% in E. Solid circles are WT, and open circles, 3×Tg-AD mice. Gray bars represent means. Statistical significance was determined using the Holm–Sidak method. Mean ± SEM: WT: 36.71 ± 1.03% in 3 mo, 32.11 ± 0.92% in 6 mo, 36.55 ± 0.92% in 9 mo, and 36.32 ± 0.98% in 12 mo; 3×Tg-AD: 34.10 ± 0.88% in 3 mo, 33.67 ± 0.92% in 6 mo, 41.24 ± 0.90% in 9 mo, and 40.87 ± 1.65% in 12 mo. *P < 0.05.

Hippocampal cells from 3×Tg mice were monitored for DNA damage by comet assay. DNA fragmentation increased from 9 mo in an age-dependent manner, suggesting that the balance between DNA damage and repair was significantly compromised along with tau accumulation in 3×Tg-AD mice (mean ± SEM: for 9 mo, 36.55 ± 0.92% in WT, vs. 41.24 ± 0.90% in 3×Tg; for 12 mo, 36.32 ± 0.98% in WT, vs. 40.87 ± 1.65% in 3×Tg) (Fig. 5 ). These data indicated that the presence of aggregated tau induced mislocalization and insolubility of BRCA1, leading to accelerated DNA fragmentation.

BRCA1 expression was up-regulated by promoter demethylation in human AD brains. To determine whether this phenomenon could be reproduced in the brain of an AD model mouse, we measured the methylation level of the Brca1 gene in 3×Tg-AD mice. Pyrosequencing clearly showed that the Brca1 promoter in hippocampal DNA from 3×Tg-AD mice was significantly hypomethylated compared with that of WT mice (), suggesting that Aβ-induced DNA damage leads to up-regulation of BRCA1 expression through demethylation of the promoter region.

Discussion

Genome-Wide Neuron-Specific DNA Methylation Analysis Reveals Previously Undiscovered Genes.

The methodology adopted in this study initially focused on genome-wide neuron-specific methylome analysis of the postmortem brain and the successful identification of biologically relevant molecules. Previous studies of AD epigenetics have used bulk DNA from postmortem brains and, in several instances, have demonstrated aberrant CpG methylation at genes related to familial AD and APP processing (24) or newly discovered genes (25, 26), which were not confirmed in our neuron-specific analysis (). This discrepancy could be attributed to the level of resolution reached by the Infinium microarray [the exact CpGs detected in the previous study (24) are not on the array] or to the advantage of our neuron-specific analysis over a bulk DNA analysis made up of various cell types. Our FACS purification strategy demonstrated that the DNA amount of neuronal/nonneuronal cell ratios varied among samples from 1:0.19 to 1:2.8 [mean, 1:1.20 ± 0.65 (SD)], which could easily skew the methylome analysis, given that patterns of methylation differ among cell types (16, 27).

The sample numbers used in our study are relatively small, which does limit our capacity to perform precise statistical analyses such as multiple comparisons. These analyses were nonetheless performed; however, none of our probes was able to meet statistical criteria. Recent studies of site-directed methylation editing have revealed that methylation changes over broader ranges (rather than a single mark) are responsible for the expression level of mRNA (28). Thus, strict correction of multiple comparisons with statistical methods that have been useful in detecting single, differentially methylated spots may not be as supportive for detecting broader differential regions reflecting the transcriptome. The goal of our study was to identify a biologically relevant molecule for AD through methylome analysis. Therefore, rather than correcting multiple comparisons by statistically stringent method, we attached more importance to the discovery of methylation signature related to transcriptional regulation. Additionally, we chose to directly confirm results at the mRNA or protein level rather than simply replicating the methylome analysis in a different sample sets, since we considered it as a screening method for intraneuronal environment alteration in AD brains. We believe that neuron-specific analysis can reduce the intragroup variation of methylation level and enhance the specificity and sensitivity in detecting DMRs.

Methylation levels were validated independently in the DMR of interest. Moreover, of three candidate genes extracted from the methylome analysis, two corresponding proteins were found to be differentially expressed in AD brains, strongly supporting the validity of our method. Methylome analysis is not a method to directly analyze protein expression in the postmortem brain. However, given the limitations of previously reported methodologies, such as postmortem degradation and sampling bias, we believe that this methodology produced a more precise representation of molecular neuronal pathology compared with previous studies. In doing so, it has the potential to be used for further in-depth analysis of the pathomechanisms underlying AD and other neurodegenerative diseases.

Methylome Analysis Reveals a Unique Pathological Pathway.

Using methylome analysis, we found that BRCA1 and AURKC were up-regulated in AD brains. We also demonstrated that the methylation levels at the promoter regions of these genes were, to some extent, correlated with the number of APOE ε4 alleles, whose gene product is known to impair Aβ clearance in the brain compared with other alleles, such as ε2 or ε3 (29). Polymorphic variants of TREM2 gene are also the known risk factor influencing for developing AD on the same level as one copy of the APOE ε4 allele (30); however, there was no polymorphic variants detected in our cases. Together with our analyses using postmortem human and mouse brains, this supports our hypothesis that increased Aβ burden induces up-regulation of BRCA1 and AURKC to repair DSBs. However, BRCA1 protein was mislocalized to the cytoplasm and predominantly found in the detergent-insoluble fraction of both human AD and 3×Tg-AD mice brains that suggested it might be dysfunctional (31). Indeed, significant DNA fragmentation could be seen in both human AD and mouse model brains. Similar finding has also been reported in a previous study (32); however, our data have suggested an alternative role for BRCA1. In addition to the precision and relevance of our results, we provide three pieces of information that also support our interpretation. First, the methylome, transcriptome, proteome, and neuropathological analyses reported here consistently demonstrate BRCA1 up-regulation during the pathological process of AD. Second, we and others have shown that the molecular weight of BRCA1 is always >250 kDa (33, 34) (), but the band indicated on immunoblots from a conflicting report (32) appears to be much smaller. In addition, knockdown of human and murine BRCA1 by multiple shRNAs confirmed that the band was BRCA1 (). Third, our data corroborated the previous result demonstrating that large number of neurons with neurofibrillary tangles were positive for BRCA1 (35).

Differential extraction demonstrated the presence of Sarkosyl-insoluble BRCA1 in the hippocampus of 3×Tg AD mice at 12 mo without SDS-insoluble BRCA1. SDS-insoluble tau represents a later and more extensive state of aggregation than Sarkosyl-insoluble tau (36, 37). In human AD, the disease time course is quite long, usually 10–15 y from the beginning of tau aggregation to autopsy, whereas mice have a life expectancy of only 2 y. Considering the results of in vitro seed-dependent tau aggregation model and immuno-EM, BRCA1 aggregates in a tau-dependent manner. The amount of insoluble tau is very small in mice model compared with human AD brain (Figs. 2 and 5 ), even though in these mice mutant tau is more aggregate-prone and is also overexpressed, suggesting that the time factor is important in tau aggregation and its change in solubility. Thus, in the mouse model, tau is not insoluble enough compared with human disease as previously reported (38). Considering that coaggregation of BRCA1 with PHF was demonstrated by EM, we speculate that a relatively short time (e.g., 12 mo) is sufficient for endogenous BRCA1 to become Sarkosyl-insoluble and dysfunctional, but insufficient to become a highly aggregated state like SDS-insoluble.

BRCA1 is a DNA repair protein, whose mutation drastically increases the risk of breast cancer (39). It is well studied in the field of cancer research, but its role in neurodegeneration remains unknown. Recent reports have indicated the importance of DNA damage in AD pathogenesis (40), with evidence of activity-induced DNA breaks governing the expression of neuronal early-response genes (41). Despite this potential relevant activity, its role remains elusive (32, 35, 40, 42). Since differentiated neurons do not proliferate, impaired DNA cannot be repaired by the precise process of homologous recombination, which uses homologous sister chromatids as templates during mitosis. Thus, nonhomologous end joining, including microhomology-mediated end joining, is the only alternative for neurons. Preserving the integrity of neuronal DNA at the best possible level over its life span would be a key to maintaining the cellular functions of neurons (43–45). Therefore, the recruitment of DNA repair proteins, such as BRCA1, would appear to be important for neurons. Interestingly, several epidemiological studies reported that the risk for developing AD is increased by diabetes mellitus or atherosclerosis, which are known to induce DNA damage by producing oxidative stress (46–48). Furthermore, an inverse association between cancer and AD has been reported (49–52). These observations suggest that deterioration in the maintenance of DNA integrity could play an important role in AD pathogenesis. We demonstrated that knockdown of BRCA1 in neuronal cells in the presence of Aβ burden in vitro and in vivo resulted in a decrease in the number of short neurites and dendritic spines. Considering the decreased spine density in brains of human AD patients and AD model mice (53, 54), these results suggest a potential mechanism underlying neuronal dysfunction in the AD brain. That is, Aβ may impair synaptic plasticity in neuronal cells by causing deterioration of DNA integrity, thereby resulting in memory impairment.

BRCA1 Dysfunction Leads to Deterioration of DNA Integrity.

Our results clearly demonstrate that Aβ-induced DSBs were efficiently repaired by BRCA1 in the absence of aggregated tau, and that repair by this mechanism was impaired in its presence due to mislocalization and insolubilization of BRCA1. Among many DNA repair markers, only γ-H2ax was positive in AD brain. Considering γ-H2ax is a scaffold for recruitment of DNA repair proteins at an early stage of DNA repair (55–57). BRCA1 is one of the early-phase proteins recruited to the γ-H2ax foci (58) and, therefore, could possibly fail to form the foci of DNA repair proteins due to its mislocalization and insolubilization impeding the subsequent DNA repair process. A recent study showed that BRCA1 nuclear export was related to susceptibility to DNA damage (59), which suggests that the mislocalization of BRCA1 to the cytoplasm and insolubility once there, interfered with DNA damage repair. Our data with cellular and animal model only with Aβ burden suggested that BRCA1 is up-regulated for a protective purpose. However, in human AD brain, BRCA1 is coaggregated with tau that makes BRCA1 dysfunctional. Despite its increased expression, coaggregation leads to insufficient DNA repair. In other words, although BRCA1 is up-regulated, it does not reach to the sufficient level necessary to counteract Aβ-induced DNA damage leading to its accumulation.

BRCA1 is involved primarily in homologous recombination, but also, to some extent, in nonhomologous end joining (60). This raises the question of why BRCA1, rather than other DNA repair proteins, is recruited in response to DNA damage in nondividing neuronal cells. Several reports have shown the importance of BRCA1 function in the developing brain (61, 62) and the accumulation of somatic mutations in neuronal cells during the developmental process in a transcriptional activity-dependent manner (63). However, activity-induced DNA DSBs have previously been shown to govern the expression of neuronal early-response genes (41). These reports collectively suggest that DNA DSBs occur frequently in neurons of highly active regions, such as the hippocampus and entorhinal cortex, and are efficiently repaired by BRCA1 in normal brains. However, in AD brains, this maintenance system seems to be impaired by BRCA1 dysfunction, resulting in neuronal deterioration in these regions. Indeed, comet assay on postmortem brains showed increased DNA fragmentation in AD samples.

In this study, we present findings of neuron-specific DNA methylation alterations in AD brains. Taken together with expression analyses, these results suggest that BRCA1 plays an important role in the pathogenesis of AD, since methylome, transcriptome, proteome, and neuropathology results are all consistent with its up-regulation. Aβ induces DNA DSBs, which should be efficiently repaired by BRCA1 when aggregated tau is absent. The balance of DNA damage and repair is compromised by dysfunction of BRCA1, due to cytoplasmic mislocalization and consequential insolubility as a result of aggregated tau accumulation (Fig. 6). In this study, methylome analysis has directly facilitated identification of biologically relevant molecules for a neurodegenerative disease. Future analyses of neuron-specific methylome may also prove useful in understanding the mechanisms of pathogenesis for sporadic neurodegenerative disorders.

Fig. 6.

Schematic illustration summarizing our current hypothesis regarding Aβ-induced DNA damage, tau, and epigenetic regulation of BRCA1 in AD. In normal brain without Aβ or tau accumulation, there is no need for BRCA1 up-regulation (Left). At an early stage of AD with no accumulated tau, BRCA1 efficiently repairs DNA DSBs induced by toxic Aβ (Middle). However, at an advanced stage of AD, cytoplasmic aggregated tau sequestrates BRCA1 to an insoluble fraction, resulting in its dysfunction (Right). While neurons try to cope with this situation by up-regulating expression of the BRCA1 gene through epigenetic mechanisms, they are eventually overwhelmed by the accumulation of DNA damage.

Materials and Methods

Sample Preparation and DNA Extraction.

This research received approval from the ethics committee of the University of Tokyo (approval 2183-15). We collected postmortem brains with written consent from patients’ families and maintained them at −80 °C until use. According to criteria established by Braak and McKeith (64–67), trained neuropathologists made a diagnosis of AD, DLB, or NC, using hematoxylin–eosin, Nissl and silver staining, and immunostaining (68). We diagnosed samples fulfilling Braak stage ≥4 and amyloid stage ≥C as AD (64, 65), and samples fulfilling Lewy body score ≥4, Braak stage ≤3, and amyloid stage ≤B as DLB (66, 67). Brain samples of 18 AD, 15 NC, and 21 DLB subjects were obtained from Tokyo Metropolitan Geriatric Hospital brain bank; 11 AD and 3 NC subjects were from the University of Tsukuba; and 1 AD and 12 NC subjects were from the University of Tokyo. Detailed clinical demographics of the postmortem brains, including sex, age, and APOE ε4 and TREM2 variants, are shown in Table 1 and . For the preparation of neuronal nuclei, we excised the gray matter from the inferior temporal gyrus and the cerebellum. Nuclear preparation and FACS nucleus sorting (FACSAria; BD Biosciences) were performed at the Support Unit for Bio-Material Analysis in RIKEN Brain Science Institute, according to a published protocol (16). Following the extraction of genomic DNA by the phenol–chloroform method, its concentration was measured using a Qubit dsDNA BR assay kit (Invitrogen).

Infinium 450k Methylation Microarray Analysis.

Genomic DNA (500 ng) was bisulfite-converted using an EZ DNA methylation kit (Zymo Research) and analyzed using an Illumina Infinium HumanMethylation450 BeadChip (Illumina) and iScan scanner (Illumina). We processed the methylation array data using the statistical computing package R 3.0.2 (www.r-project.org/) and the Bioconductor package ChAMP (69). Idat files were imported to calculate methylation levels (β values) for each probe. The β-mixture quantile (BMIQ) normalization method (70) was used to reduce probe-type bias, and batch correction was performed using the ComBat method (71). For each probe, we compared β values between AD and NC samples using two one-sided Student’s t tests: P (AD < NC) and P (AD > NC). We defined probes with P values < 0.05 and mean β differences between AD and NC > 0.05 as DMPs. Regions with no less than three consecutive DMPs were defined as DMRs.

Single-Cell Gel Electrophoresis (Comet Assay).

The comet assay was performed under neutral pH conditions, using a CometAssay kit (Trevigen), following the manufacturer’s instructions. DNA was visualized with SYBR Gold, and images were obtained using an Axioplan 2 fluorescent microscope (Carl Zeiss) and an Axiocam HRc CCD camera system (Carl Zeiss). The head and tail of each individual comet were detected semiautomatically to calculate Tail DNA% using CASPLab software (72).

Additional material and methods are included in .