Aberrant DNA methylation associated with Alzheimer's disease in the superior temporal gyrus.

Abnormal DNA methylation patterns have been demonstrated to be associated with the pathogenesis of Alzheimer's disease (AD). The present study aimed to identify differential methylation in the superior temporal gyrus (STG) of patients with late-onset AD based on epigenome-wide DNA methylation data by bioinformatics analysis. The genome-wide DNA methylation data in the STG region of 34 patients with late-onset AD and 34 controls without dementia were recruited from the Gene Expression Omnibus database. Through systemic quality control, differentially methylated CpG sites were determined by the Student's t-test and mean methylation value differences between the two conditions. Hierarchical clustering analysis was applied to assess the classification performance of differentially methylated CpGs. Functional analysis was performed to investigate the biological functions of the genes associated with differentially methylated CpGs. A total of 17,895 differentially methylated CpG sites were initially identified, including 11,822 hypermethylated CpGs and 6,073 hypomethylated CpGs. Further analysis examined 2,211 differentially methylated CpGs (covering 1,991 genes). AD subjects demonstrated distinctive DNA methylation patterns when compared with the controls, with a classification accuracy value of 1. Hypermethylation was mainly detected for genes regulating the cell cycle progression, whereas hypomethylation was observed in genes involved in transcription factor binding. The present study demonstrated widespread and distinctive DNA methylation alterations in late-onset AD. Identification of AD-associated epigenetic biomarkers may allow for the development of novel diagnostic and therapeutic targets.

Introduction

Alzheimer's disease (AD) is a chronic degenerative brain disorder and the leading cause of dementia in the elderly (1,2). AD is characterized by gradual neuronal death, memory loss and impaired cognitive ability, and contributes to 60–70% of dementia cases (2). Current advances are able to temporarily improve symptoms, while no effective treatment is available that can cure AD or reverse the progressive course (3). Understanding the risk factors for AD contributes to the investigations on reducing the disease occurrence. The vast majority of patients with AD have late-onset disease, and the susceptible risk is determined by complex mechanisms (2,4).

Previous studies have revealed the roles of genetic, environmental and epigenetic factors in AD risk (4–6). Epigenetics refers to DNA modifications that alter gene expression and phenotype without causing changes to the nucleotide sequences (7). As the best-studied epigenetic modification, DNA methylation is closely associated with several key cellular processes, including cell differentiation, gene regulation and genomic imprinting (8–10). Several lines of evidence has indicated the effect of DNA methylation on the pathogenesis of AD (11–13). Furthermore, the genome-wide methylation profiling technique has been applied for the characterization of methylation patterns across the genome, and has been considered as a reliable tool in identifying methylation differences associated with complex diseases (13,14).

Numerous studies have focused on the identification of specific brain regions that are particularly vulnerable throughout the progression of AD (15–17). Haroutunian et al (18) performed an extensive study on the transcriptional vulnerability involved in the progression of late-onset AD across 15 brain regions, and observed that the superior temporal gyrus (STG) region presented significant transcriptional abnormalities. Thus, the present study focused on the STG region of patients with late-onset AD. Recently, Watson et al (19) performed a genome-wide DNA methylation profiling in the STG region of patients with late-onset AD, and focused on the identification of differentially methylated regions. In the present study, the publicly available DNA methylation data of the study by Watson et al (19) were used to extract AD-associated methylated genes and identify the underlying biological processes impacted by hypomethylated and hypermethylated genes. Furthermore, epigenetic biomarkers in the STG region associated with AD risk were identified in the current study, and the findings may contribute to the development of novel diagnostic and therapeutic targets for late-onset AD.

Materials and methods

Genome-wide DNA methylation data

The genome-wide DNA methylation data of late-onset AD were obtained from the Gene Expression Omnibus database (GEO; https://www.ncbi.nlm.nih.gov/geo/), under the GSE76105 accession number (19). The current study was performed using the Illumina Infinium Human Methylation 450 array platform (Illumina, Inc., San Diego, CA, USA). The data analysed from the database included tissue samples obtained from the STG region of 34 patients with confirmed late-onset AD and 34 controls without dementia. Detailed information of these tissue samples is described in a previous study (19).

Prior to analysis, the raw DNA methylation data were subjected to a rigorous preprocessing procedure following the method developed by Huynh et al (20). Probes were removed from the data according to the following criteria: i) Single nucleotide polymorphisms within 5 bp upstream of the targeted CpG; ii) minimum allelic frequency of <0.05; iii) probes on X and Y chromosomes; and iv) cross-hybridising probes. Subsequent to applying the exclusion criteria, a total of 424,497 CpGs were selected and used for subsequent analysis. These remaining CpGs were quantile normalized using the lumi package in R software (21) and the beta-mixture quantile method (22).

Differential methylation analysis

In the dataset, the methylation values for individual CpGs in each sample were expressed as β-values. β-value is a quantitative measure of methylation for each CpG site, ranging between 0 (completely unmethylated site) to 1 (fully methylated site). Initially, the β-values of individual CpG sites were obtained in AD subjects and controls, respectively. Next, the mean β-value difference of each CpG between two conditions was calculated, and CpG sites with an absolute mean β-value difference of >0.05 were included into the study. For each probe, β-values in AD subjects were compared against the controls using the Student's t-test. The CpGs with an absolute mean β-value differences of >0.05, as well as P<0.05, were considered to be differentially methylated sites.

Previous studies indicated that methylation alterations rarely exist in genomic regions at the extreme ends of methylation values (<0.2 or >0.8), but are preferentially detected in genomic locations at intermediate methylation levels, which is a range associated with active distal regulatory regions, such as enhancers (20,23). Thus, in order to reduce the number of non-variable CpGs and improve the statistical power of subsequent analyses, the CpG sites with β-values of >0.8 and <0.2 in all samples were removed from the current study. In addition, it is reported that numerous differentially methylated CpGs with a β-value difference of <0.2 may be difficult to reproduce by alternative methodologies or in replication studies (24). Therefore, only CpGs that had a mean β-value difference of ≥0.2 were retained in the current study.

Hierarchical clustering analysis

Hierarchical clustering is a method of cluster analysis that seeks to build a hierarchy of clusters in data mining and statistics. Theoretically, samples with similar methylation profiles can be clustered together. To assess the classification performance of differentially methylated CpGs, a hierarchical clustering analysis was applied using the Euclidian distance and average linkage criteria (25). Ideally, the samples should be classified into two distinct clusters, namely the AD subjects and controls. In order to assess the classification efficiency, the accuracy was measured, which is a fraction of correctly classified samples over all samples, according to the following formula: Accuracy=(TP+TN)/(TP+FP+TN+FN). TP represents the number of positive samples correctly predicted as positive, TN is the number of negative samples correctly predicted as negative, FP is the number of negative samples incorrectly predicted as positive, and FN represents the number of positive samples incorrectly predicted as negative.

Functional enrichment analysis

To further investigate the biological functions of the genes associated with differentially methylated CpGs, functional enrichment analysis was performed based on Gene Ontology (GO; http://www.geneontology.org) and Kyoto Encyclopedia of Genes and Genomes database (KEGG; http://www.genome.jp/kegg/pathway.html). In addition, enrichment analysis was conducted using the online software Database for Annotation, Visualization and Integrated Discovery (DAVID 6.8; https://david.ncifcrf.gov). Fisher's test was utilized to measure the significance of GO terms and biological pathways. The P-values were adjusted using Benjamini-Hochberg false discovery rate (FDR). Statistically significantly different terms were determined under the criterion of P<0.01.

Results

The present study focused on comparing the DNA methylation in the STG region of 34 late-onset AD subjects with that in 34 controls without dementia. Following quality control processing of the DNA methylation data, probes were removed under the filtering criteria described earlier, and methylation data for 424,497 autosomal CpGs in each of the 68 individuals were used in subsequent differential methylation analysis. Under the criteria of absolute mean β-value differences between AD subjects and controls of >0.05, as well as P<0.05, a total of 17,895 differentially methylated CpG sites covering 8,678 genes were initially identified. Among these differentially methylated CpGs, a total of 11,822 CpGs were hypermethylated and 6,073 CpGs were hypomethylated. A volcano plot exhibiting the distribution of differentially methylated CpGs is shown in Fig. 1.

Figure 1.

Volcano plot exhibiting the distribution of differentially methylated CpGs in patients with AD compared with those in control individuals. The x-axis represents the mean methylation differences between AD patients and controls, while the y-axis represents the log transformed P-values. AD, Alzheimer's disease.

Since methylation changes were preferentially detected in genomic regions at intermediate methylation levels, further selection was performed to reduce non-variable CpGs and improve the statistical power. By removing CpGs at the extreme ends of methylation values (β-values of >0.8 and <0.2), 2,225 differentially methylated CpGs covering 2,001 genes were identified. Furthermore, differentially methylated CpGs with a β-value difference of <0.2 were considered to be difficult to reproduce and excluded. Thus, only differentially methylated CpGs with a mean β-value difference of ≥0.2 were extracted, and 2,211 differentially methylated CpGs (covering 1,991 genes) were finally identified.

Hierarchical clustering was implemented to determine whether the identified differentially methylated CpGs can be applied to distinguish AD subjects from controls. The results illustrated that these two conditions in the AD and control groups exhibited distinctive DNA methylation patterns (Fig. 2). In addition, it was observed that the samples were classified into two distinct clusters by these differentially methylated CpGs with an accuracy of 1.

Figure 2.

Hierarchical clustering analysis of differentially methylated CpGs between Alzheimer's disease samples and controls. The differentially methylated CpGs were able to distinctly distinguish between the two types of samples.

GO terms and KEGG pathway enrichment analysis for genes associated with the 2,211 differentially methylated CpGs were performed using the DAVID online tool. Under the criterion of P<0.01, hypermethylated genes were significantly enriched in 10 GO terms, with the top three (ranked according to their P-value) involved in the regulation of the mitotic cell cycle, regulation of the mitotic cell cycle phase transition, and the Wnt signaling pathway (Fig. 3). The hypomethylated genes were significantly enriched in 5 GO terms, with the top one associated with transcription factor binding (Fig. 4). However, KEGG pathway analysis demonstrated that no statistically significant pathways were enriched by both the hypermethylated and hypomethylated genes under an FDR-adjusted P-value of <0.01.

Figure 3.

Gene oncology analysis of genes associated with hypermethylated CpGs in Alzheimer's disease subjects.

Figure 4.

Gene oncology analysis of genes associated with hypomethylated CpGs in Alzheimer's disease subjects.

Discussion

Recently, epigenetic modification in diseases has attracted increasing attention from researchers. Late-onset AD results from the interactions of multiple genetic and environmental factors, along with the disruption of epigenetic mechanisms controlling gene expression (26). Abnormal DNA methylation patterns have been observed to be associated with the pathogenesis of AD. In the present study, the epigenome-wide DNA methylation patterns in the STG of patients with late-onset AD were assessed, and differentially methylated genes and biological functions were identified under the AD condition, which may contribute to the development of novel targets for late-onset AD diagnosis and therapy.

The present study revealed significant alterations to the DNA methylation profiles in AD, with 2,211 CpGs displaying significantly differential methylation levels. The hierarchical clustering results illustrated that AD subjects exhibited distinctive DNA methylation patterns compared with the controls, which was consistent with the observations of a previous study (19). Initially, the differential methylation identification demonstrated a strong bias to hypermethylated alterations in late-onset AD (11,822 hypermethylated and 6,073 hypomethylated CpGs). Given that the greatest risk factor for late-onset AD is older age (27), it may be hypothesized that CpG methylation increased with age. Previous studies also indicated that CpG methylation in AD was significantly associated with age (13,28).

To further examine the functions of genes associated with differentially methylated CpGs, GO and KEGG analysis were performed in the present study. The results revealed that genes associated with hypermethylated CpGs were enriched in 10 GO terms, 6 of which were cell cycle-associated biological process. The cell cycle involves a series of events that lead to cell division and DNA replication to produce two daughter cells. Cell cycle progression ensures a correct inheritance of epigenetic modifications. A previous study has indicated that epigenetic modifications arising within the cells of an individual are transmitted through multiple cell cycles, which is crucial in maintaining a given chromatin state (29). In the current study, it was identified that genes associated with hypermethylated CpGs participated in two Wnt-associated biological processes, including the Wnt signaling pathway and cell-cell signaling by Wnt. The Wnt signaling pathway is a vital signal transduction pathway involved in various biological activities. Increasing evidence suggested the protective effects of the Wnt pathway on synaptic modulation and cognitive processes in neurodegenerative diseases, implying the potential role of aberrant Wnt signaling in cognitive impairment (30,31). Riise et al (32) also confirmed an aberrant Wnt signaling pathway in medial temporal lobe structures of neurodegenerative AD samples. Recent studies have further affirmed the potential role of the Wnt signaling pathway as a therapeutic target for the treatment of AD (33,34).

In the present study, hypomethylation was detected in genes involved in transcription factor binding and core promoter binding. Numerous studies have revealed the effects of DNA methylation on the transcription factor binding (35–37). Transcription factors serve important roles in the epigenetic regulation of gene transcription, while DNA methylation prohibits the recruitment of transcription factors, resulting in transcription suppression. A protein microarray-based approach observed that methylation-dependent transcription factor binding was a widespread phenomenon in biological processes (38). Methylated CpG sites are typically considered to reduce gene expression by preventing transcription factors from binding to promoter regions (38). Promoter site-specific CpG methylation is associated with various biological processes. Cyclical methylation alterations in promoter CpGs represent a critical event in achieving transcription (39). Pieper et al (40) also identified a tendency for hypomethylation in the tumor necrosis factor α (TNF-α) core promoter region in neurodegenerative disorders, resulting in reduced binding of the transcription factors and suppressed TNF-α promoter activity.

The focus of the present study was the STG region of patients with late-onset AD. Increasing studies have already identified unique epigenetic signatures involved in the progression of AD in various brain regions. For instance, human AD case-cognitively normal control studies observed global hypomethylation in the entorhinal cortex (41) and in the temporal neocortex neuronal nuclei (42). Additionally, Bakulski et al (13) demonstrated widespread discordant DNA methylation in AD in the human frontal cortex. Chouliaras et al (43) identified a consistent decrease in global DNA methylation and hydroxymethylation in the hippocampus of AD subjects. Furthermore, the present study demonstrated a strong bias to hypermethylated alterations in the STG region in AD patients, implying different methylation changes across various brain regions.

In conclusion, the present study demonstrated significant and distinctive DNA methylation patterns in the STG region of patients with late-onset AD. Hypermethylation was mainly detected for genes regulating the cell cycle progression, while hypomethylation was identified in genes involved in transcription factor binding. These results suggested AD-associated epigenetic marks, potentially contributing to the understanding of underlying mechanisms involved in the epigenetic regulation of AD and the development of novel therapeutic targets for AD.