Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Neurosci., 13 April 2023
Sec. Brain Disease Mechanisms

Brain region-specific genome-wide deoxyribonucleic acid methylation analysis in patients with Alzheimer’s disease

\r\nGang RenGang Ren1Shan Song,Shan Song2,3Sheng-Xiao Zhang,Sheng-Xiao Zhang2,3Yan Liu,Yan Liu2,3Yan LvYan Lv4Yan-Hong WangYan-Hong Wang1Rong Zhao,Rong Zhao2,3Xin-Yi Li*Xin-Yi Li1*
  • 1Department of Neurology, Third Hospital of Shanxi Medical University, Shanxi Bethune Hospital, Shanxi Academy of Medical Sciences, Tongji Shanxi Hospital, Taiyuan, China
  • 2Department of Rheumatology, The Second Hospital of Shanxi Medical University, Taiyuan, China
  • 3Key Laboratory of Cellular Physiology at Shanxi Medical University, Ministry of Education, Taiyuan, China
  • 4Department of Nephrology, Third Hospital of Shanxi Medical University, Shanxi Bethune Hospital, Shanxi Academy of Medical Sciences, Tongji Shanxi Hospital, Taiyuan, China

Objective: Alzheimer’s disease (AD) is a neurodegenerative disease characterized by neuropathology and cognitive decline and associated with age. The comprehensive deoxyribonucleic acid methylation (DNAm)-transcriptome profile association analysis conducted in this study aimed to establish whole-genome DNAm profiles and explore DNAm-related genes and their potential functions. More appropriate biomarkers were expected to be identified in terms of AD.

Materials and methods: Illumina 450KGSE59685 dataset AD (n = 54) and HC (n = 21) and ribonucleic-acid-sequencing data GSE118553 dataset AD patients (n = 21) and HCs (n = 13) were obtained from the gene expression omnibus database before a comprehensive DNAm-transcriptome profile association analysis, and we performed functional enrichment analysis by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses (KEGG). Three transgenic mice and three wild-type mice were used to validate the hub genes.

Results: A total of 18,104 DNAm sites in healthy controls (n = 21) and AD patients (n = 54) were surveyed across three brain regions (superior temporal gyrus, entorhinal cortex, and dorsolateral prefrontal cortex). With the addition of the transcriptome analysis, eight hypomethylated-related highly expressed genes and 61 hypermethylated-related lowly expressed genes were identified. Based on 69 shared differentially methylated genes (DMGs), the function enrichment analysis indicated Guanosine triphosphate enzymes (GTPase) regulator activity, a synaptic vesicle cycle, and tight junction functioning. Following this, mice-based models of AD were constructed, and five hub DMGs were verified, which represented a powerful, disease-specific DNAm signature for AD.

Conclusion: The results revealed that the cross-brain region DNAm was altered in those with AD. The alterations in DNAm affected the target gene expression and participated in the key biological processes of AD. The study provides a valuable epigenetic resource for identifying DNAm-based diagnostic biomarkers, developing effective drugs, and studying AD pathogenesis.

1. Introduction

Alzheimer’s disease (AD) is one of the most common age-related chronic neurodegenerative diseases that is characterized by progressive cognitive dysfunction and memory loss. The pathological features of AD involve the accumulation of amyloid-beta (Aβ) and neurofibrillary tangles. At present, there is no cure or treatment that can halt the progression of the disease. While the occurrence and development of AD is 70% dependent on genetic heterogeneity (Ballard et al., 2011), there are various risk factors, including cardiovascular risk factors (diabetes, hypertension, and obesity), age, and cognitive reserve (Daviglus et al., 2010; Mayeux and Stern, 2012). The involvement of specific epigenetic factors in the development of AD, including both pathological and physiological conditions, has become increasingly clear, (Berson et al., 2018; Perkovic et al., 2021) with deoxyribonucleic acid methylation (DNAm) being the most common epigenetic phenomenon. DNAm is crucial to regulating synaptic plasticity and maintaining the central nervous system (Perkovic et al., 2021). Previous research demonstrates that different regions of the brain exhibit distinct manifestations of AD (Wenk, 2003). Compared with the cerebellum, epigenetic perturbations in AD occur more frequently among multiple cortical regions, including the entorhinal cortex, prefrontal cortex, and temporal gyrus (Jager et al., 2014; Semick et al., 2019).

While previous studies successfully identify various DNAm differences associated with AD (Mastroeni et al., 2009; Chouliaras et al., 2013; Coppieters et al., 2014), existing research on methylation only identifies differentially methylated genes (DMGs) in a single brain region, while the relationship between DNAm and gene expression is not yet fully understood. Therefore, the present study aims to systematically assess multiple brain regions to identify DNAm alterations to drive the functionality of a transcriptome. A better understanding of the epigenetic mechanism of gene expression in AD could enable the discovery of more viable strategies for the prevention and treatment of the disease.

2. Materials and methods

2.1. Ribonucleic acid isolation, construction of nucleic acid library, and illumina sequencing

The feeding environment of SPF level mice should meet the following requirements: temperature 18∼22°C, Relative humidity: 40–70%, the humidity is 50∼60%. The experiment was carried out after 1 week of adaptive feeding. Mice were used three biological replicates, each of which comprised three frozen cortices from triple transgenic (3xTg) mice (females aged 5 months) that were processed for ribonucleic acid sequencing (RNA-seq). Meanwhile, the healthy control (HC) group contained three wild-type (WT) mice (females aged 5 months). The animals used in all experiments followed the National Institutes of Health guidelines, and all protocols were approved by the Institution Animal Care and Use Committee of Shanxi Medical University. The total RNA was extracted using TRIzol™ reagents (Invitrogen) in accordance with the manufacturer’s instructions, while the DNA was removed using DNase I reaction buffer (New England Biolabs). Prior to the RNA-seq, the quality of the RNA was assessed using the Agilent Bioanalyzer 2100 system, with 1.5 μg of the RNA sample subsequently taken for the sequencing procedure. The library preparations were sequenced on the Illumina system, with clean reads obtained by removing any reads containing adapter or ploy-N sequences and any low-quality reads containing > 50% bases with qualities of < 20 from the raw data. Meanwhile, reference genome and gene model annotation files were downloaded from the genome website directly, with paired-end clean reads aligned to the reference genome using Hisat2 v.2.0.5 software. FeatureCounts v.1.5.0-p3 was used to count the read numbers mapped to each gene.

2.2. Data collection

Illumina’s human methylation 450k array and specific gene expression profile datasets of AD were acquired from the gene expression omnibus database. In this search, the following three brain regions in a cohort of individuals were studied: the entorhinal cortex, the superior temporal gyrus, and the prefrontal cortex. The DNAm profiles of the AD (n = 54) and HC (n = 21) groups were obtained from the GSE59685 dataset (a cross-tissue analysis of methylomic variation in AD using samples from three independent human post-mortem brain cohorts), in addition to the transcriptome analysis data pertaining to AD patients (n = 21) and HCs (n = 13) from the GSE118553 individuals with intact cognition and neuropathology consistent with Alzheimer’s disease (AD) datasets. A flowchart of the study approach is shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. A flowchart of the study approach.

2.3. Data pre-processing

To estimate the methylated intensity at individual 5’-cytosine–phosphate–guanine-3’ (CpG) sites, raw intensity files were preprocessed using Illumina’s (Taub et al., 2013) internal normalization probes and algorithms package, lumiMethyR, a pre-processing package that includes color balance adjustment, background level correction, and data normalization. The relative DNAm levels, represented as beta values (β), were derived from the intensities. Comparisons of AD-affected and normal brain tissues, differentially methylated probes (DMPs), differentially methylated regions (DMRs), and differentially methylated blocks (DMBs) were carried out using ChAMP (Tian et al., 2017). By applying the default setting, the raw data were filtered for probes with a detection value of non-CpG sites, the multi-hit probe list, or X and Y chromosomes. All DMPs with an adjusted P-value of < 0.05 were selected for further analysis.

The expression profile was preprocessed using the oligo package (Carvalho and Irizarry, 2010). Here, probes with more than one gene were removed, along with empty probes and duplicate gene symbols, with the maximum values regarded as the gene manifestation. The differentially expressed genes (DEGs) among the AD and HCs groups were screened using the limma R package (Federico and Monti, 2020). P < 0.05 set as the cut-off criterion for the DEGs.

2.4. Integration of transcriptome and deoxyribonucleic acid methylation

The hypermethylation of the DNA promoter region has traditionally been associated with transcriptional silence. Here, AD-related DMGs were identified, including hypomethylated-mediated highly expressed genes (hypo-HGs) and hypermethylated-mediated lowly expressed genes (hyper-LGs).

The finite element algorithm method (FEM) (Jiao et al., 2014) is a supervised clustering method that uses protein–protein interaction networks to perform a correlation analysis between Illumina Infinium 450k data and matched gene expression data to identify specific epigenetic gene modules or signaling pathways. Functional modules with methylation and gene expression levels undergoing significant changes can be determined using a FEM analysis.

2.5. Identification of cerebral-cortex-shared differentially methylated genes

To investigate common AD-related DMGs in the three cerebral cortex regions, Venn diagrams were created to reveal the degree of overlap between the DMGs in each brain region. A principal component analysis and unsupervised hierarchical clustering were used to examine the effects of the DMGs among the AD and HC groups. The Metascape database1 was used to conduct pathway and process enrichment analyses based on the DMGs, including gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) biological functional enrichment analyses.

2.6. Identification of the hub genes

The Boruta algorithm was used as part of a feature selection method to identify the optimal diagnostic gene biomarkers for AD. The hub gene was evaluated in terms of the receiver operating characteristics (ROCs), with the diagnostic ability subsequently compared. Mice-based models of AD are widely used to study the pathogenesis of the disease and are critical for validating candidate biomarkers (Kursa and Rudnicki, 2010).

3. Results

3.1. Identification of differentially methylated 5’-cytosine–phosphate–guanine-3’ sites and differentially expressed genes in the three regions

In the entorhinal cortex, 26,350 DMPs were identified, which were annotated to 19,685 genes, including 7,978 (59.47%) with a higher methylation status, while 11,707 (40.53%) were hypomethylated. In the transcriptomic profiling dataset, 4,570 genes were up-regulated and 2,568 were down-regulated in the entorhinal cortex. In terms of the prefrontal cortex, 20,089 DMPs were identified and annotated to 15,835 genes, including 12,635 (79.79%) with higher methylation status, while 3,200 (20.21%) were hypomethylated, and 3,487 up-regulated DEGs and 2,382 down-regulated DEGs were identified. In terms of the superior temporal gyrus, 2,824 DMPs were identified, which corresponded to 2,082 genes, including 1,797 (86.31%) with a higher methylation status, while 285 (13.69%) were hypomethylated and 2,688 up-regulated DEGs and 1,929 down-regulated DEGs were identified. More CpG sites exhibited hypermethylation, while an analysis of the genomic locations of differentially methylated CpG sites revealed that place sites existed in the body for all brain regions, followed by TSS1500, 5’-UTR, and 1stExon (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Volcano map of 5’-cytosine–phosphate–guanine-3’ (CpG) sites and differentially expressed genes (DEGs) after filtering. The three panels show the distribution characteristics of differentially methylated CpG sites and differentially expressed genes, respectively, in the (A,B) entorhinal cortex, (C,D) prefrontal cortex, and (E,F) superior temporal gyrus, as well as the ratio of hypermethylated to hypomethylated CpG sites.

3.2. Identification of Alzheimer’s-related cerebral-cortex-shared differentially methylated genes

The DMGs were assessed by examining the DNAm and then any transcriptomic reversed alterations in the three brain regions. Using a Venn diagram, eight hypo-HGs and 61 hyper-LGs were revealed (Figure 3). Collectively, 69 aberrantly methylated DEGs were identified as AD-related cerebral cortex DMGs. An unsupervised hierarchical clustering heatmap of the shared DMG signatures revealed that the samples could be separated based on specific subject groups (Figure 4).

FIGURE 3
www.frontiersin.org

Figure 3. Identification of Alzheimer’s disease (AD)-related cerebral-cortex-shared differentially expressed genes (DMGs). The Venn diagram reveals eight hypomethylated-mediated highly expressed genes (hypo-HGs) (A) and 61 hypermethylated-mediated lowly expressed genes (hyper-LGs) (B).

FIGURE 4
www.frontiersin.org

Figure 4. Heatmap of the differentially expressed genes (DMGs) among the Alzheimer’s disease (AD) and healthy control (HC) groups in three brain regions. The heatmap reveals eight hypo-HG and 61 hyper-LG expressions in the (A,B) entorhinal cortex, (C,D) prefrontal cortex, and (E,F) superior temporal gyrus. Each column represents a sample, and each row represents the transcription expression level of all the samples involved in one gene.

3.3. Functional enrichment of the differentially methylated genes

To identify the functional characteristics of the DMGs, an enrichment analysis was performed in terms of the shared DMGs (Figure 5). The GO analysis revealed that the hypo-HGs were enriched in Rho protein signal transduction, the regulation of small Guanosine triphosphate enzymes (GTPase) mediated signal transduction, and GTPase regulator activity, while the hyper-LGs significantly increased the GTPase regulator activity, the synaptic vesicle cycle, protein domain-specific binding, the positive regulation of the protein catabolic process, and the response to the growth factor. Meanwhile, the KEGG analysis results indicated that hyper-LGs were mainly involved in the relaxin signaling pathway, as well as endocytosis, tight junction functioning, and aldosterone synthesis and secretion. These outcomes indicated that methylation-driven DEGs may play essential roles in AD development.

FIGURE 5
www.frontiersin.org

Figure 5. The gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) functional enrichment analysis results for differentially methylated genes in Alzheimer’s disease (AD). (A) The GO terms in the eight hypo-LGs of AD; (B,C) the GO and KEGG terms in the 61 hyper-LGs of AD.

3.4. Hub gene identification and validation of hub differentially expressed gene expression using the Alzheimer’s mice model

A total of 16 DMGs were identified as optimal biomarkers for AD via feature selection using the Boruta algorithm. Following this, a ROC analysis was performed to evaluate the diagnostic utility of the 16 shared DMGs. Ultimately, a 14-gene combination was left behind, the area under the curve of which was ≥ 0.7, including three up-regulated and 11 down-regulated genes. Finally, in the AD mice model, We tested all the candidate genes and found that a total of five genes were confirmed with consistent differences as follows: epidermal growth factor receptor pathway substrate 8 (EPS8), zinc finger MIZ-domain containing 1 (ZMIZ1), telomere maintenance 2 gene (TELO2), protein kinase C zeta (PRKCZ), and mannosidase alpha class 1B member 1 (Man1b1) (Figure 6 and Table 1).

FIGURE 6
www.frontiersin.org

Figure 6. Five significantly expressed genes in the Alzheimer’s disease (AD) model mice compared with the control group. Five of the 14 genes had a significant expression level among the AD and control groups (*P < 0.05 and **P < 0.01). The red and blue boxes represent the AD and control mice, respectively.

TABLE 1
www.frontiersin.org

Table 1. The five genes, which have been validated by Alzheimer’s mice model, show high ability to discriminate patients with AD and healthy controls.

4. Discussion

The DNAm in various brain regions is an epigenetic factor that is interrupted in AD. Over the past few years, epigenetic therapies for AD have achieved a certain amount of success (Berson et al., 2018; Chatterjee et al., 2018). To further validate the relationship between DNAm and AD, a cross-brain region analysis was conducted. In the analysis, key differential genes among the AD and HC groups were identified by integrating the AD methylation data and gene expression measures in multiple brain regions.

The GO and KEGG analyses indicated that these DMGs were closely related to the small GTPase mediated signal transduction. Small guanosine triphosphatases (small GTPases) comprises five families: Ras, Rho, Rab, Arf, and Ran. The various roles of small GTPases super familys in AD pathogenesis have been extensively reviewed (Zhu et al., 2000; Bolognin et al., 2014; Arrazola Sastre et al., 2020, 2021), several studies have examined the relationship between Rho GTPases, amyloid precursor protein (APP) synthesis, and β-amyloid (Aβ) production in various cell lines (Boo et al., 2008; Wang et al., 2009). And Rho GTPases appear to contribute to the increase in β-amyloid (Aβ) production and resulting neurotoxicity (Aguilar et al., 2017). RhoA co-localized with hyperphosphorylated tau in neurofibrillary tangles (Huesa et al., 2010). Rho protein signal transduction pathway was also enriched in our results. Rab GTPases and Arf GTPases control the processing and trafficking of toxic peptides such as the APP, the axonal transport of proteins such as membrane receptors and autophagy. Therefore, pathological conditions lead to membrane trafficking defects (Kiral et al., 2018; Arrazola Sastre et al., 2020, 2021; Jordan et al., 2022).

The functional significance of the EPS8 in the brain has only just begun to be elucidated. The EPS8 activates the GTPase, Rac, leading to actin cytoskeletal remodeling (Offenhäuser et al., 2004). This is abundantly expressed in many brain regions, (Disanza et al., 2004) mainly in neurons in the cerebral cortex and hippocampus, which are two areas typically linked to higher cognitive functions (Offenhäuser et al., 2006). In addition, the EPS8 is enriched in the growth cone and filopodia of developing hippocampal neurons, where it down-regulates axonal filopodia formation (Menna et al., 2009). In the hippocampus, the conditional deletion of the EPS8 reduces the synaptic plasticity and impairs cognitive functioning (Wang et al., 2017). The synaptic vesicle cycle pathway in our study was directly enriched. In summary, the lack of EPS8 affects the synapse number but does not significantly impact the structural organization of the presynaptic compartment (Menna et al., 2013).

The ZMIZ1-related DNAm status is regarded as a molecular marker in multiple cancer types, including astrocytoma, bladder cancer, and renal cell carcinoma (Mathios et al., 2019). While the methylation differences among ZMIZ1 genes have not been well-studied in relation to AD, there is some circumstantial evidence for the association between AD risk and the expression of ZMIZ1. ZMIZ1 expression involves the regulation of neuronal and glial cell development through the Notch signaling pathway (Angulo-Rojo et al., 2013). In addition, the possible neurologic developmental function of ZMIZ1 may serve as a transcriptional coregulator in chromatin remodeling complexes (Carapito et al., 2019).

The TELO2 gene is essential to preserving genomic stability, and any disorder of the telomere maintenance mechanism results in telomere attrition (Kass-Eisler and Greider, 2000). Telomerase persists in the neurons of the adult human brain, where it may have a protective role against tau pathology (Boccardi et al., 2015). Forero et al. (2016) reported the presence of shorter telomeres in AD patients and highlighted the importance of the analysis of epigenomic markers associated with the risk of AD. Elsewhere, a multidomain lifestyle intervention facilitated telomere length maintenance among the subgroups of older people at risk of dementia and was linked to more pronounced cognitive intervention benefits (Sindi et al., 2021).

The internal promoter in the PRKCZ gene transcribes PKMζ, which is a constitutively active atypical kinase (Hernandez et al., 2003). Lee et al. (2013) targeted the exon 9 of the PRKCZ gene for developing experimental PRKZC–/– in mice and observed that an absence of PKMζ does not impair the learning and memory in mice. However, several studies reported that PKMζ plays a significant role in maintaining synaptic long-term potentiation and long-term memory (Serrano et al., 2008). Therefore, the role of PKMζ in maintaining memory remains debatable and requires further discussion.

There are few reports related to the Man1b1 gene in AD. However, it has been widely reported in relation to Parkinson’s disease (Balducci et al., 2007; Tasegian et al., 2017; Federico and Monti, 2020; Zhang et al., 2021), with the activities of alpha (α)-mannosidase found to be substantially decreased in the cerebrospinal fluid of Parkinson’s disease patients (Balducci et al., 2007). However, Tasegian et al. (2017) recently reported that the α-mannosidase activity in cerebrospinal fluid mirrors the pathological changes in neurodegenerative disorders. Therefore, the α-mannosidase activity in AD is worthy of further exploration.

Finally, we tested all the candidate genes in the AD mice model to verify the biomarkeEPS8, ZMIZ1, TELO2, PRKCZ, and Man1b1, although the number of these genes is significantly less than that of candidate genes (possibly related to species differences), this has partly suggested that the genes we screened may be universal in the biological world, so these five genes can be well-verified. As for the accurate biological mechanism, it may need to be verified by subsequent experiments.

In conclusion, a combined analysis involving transcriptome and methylation was conducted. Further study of these hub genes in relation to AD and the pathways they participate in may provide new ideas for AD treatment.

This study had several limitations. First, only the DEGs have between verified in the two groups, while the methylation level of these genes was not verified. Therefore, it is necessary to further verify gene expression via a polymerase chain reaction or immunohistochemistry analysis.

Data availability statement

The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent from the patients/participants or patients/participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author contributions

GR, SS, S-XZ, and X-YL: conception and design of the research. GR, SS, YLv, Y-HW, and S-XZ: acquisition of data. GR, SS, and RZ: analysis and interpretation of the data. SS, Y-HW, and RZ: statistical analysis. GR, SS, and YLiu: writing of the manuscript. YLv, Y-HW, S-XZ, and X-YL: critical revision of the manuscript for intellectual content. All authors read and approved the final draft.

Funding

This study was supported by the Science Foundation of Shanxi Health Commission No. 2023010.

Acknowledgments

We are particularly grateful to all the people who have given us help on our article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2023.971565/full#supplementary-material

Footnotes

  1. ^ www.metascape.org

References

Angulo-Rojo, C., Manning-Cela, R., Aguirre, A., Ortega, A., and López-Bayghen, E. (2013). Involvement of the notch pathway in terminal astrocytic differentiation: Role of PKA. ASN Neuro 5:e00130. doi: 10.1042/AN20130023

PubMed Abstract | CrossRef Full Text | Google Scholar

Aguilar, B. J., Zhu, Y., and Lu, Q. (2017). Rho GTPases as therapeutic targets in Alzheimer’s disease. Alzheimers Res. Ther. 9:97.

Google Scholar

Arrazola Sastre, A., Luque Montoro, M., Galvez-Martin, P., Lacerda, H. M., Lucia, A. M., Llavero, F., et al. (2020). Small GTPases of the ras and rho families switch on/off signaling pathways in neurodegenerative diseases. Int. J. Mol. Sci. 21:6312.

Google Scholar

Arrazola Sastre, A., Luque Montoro, M., Lacerda, H. M., Llavero, F., and Zugaza, J. L. (2021). Small GTPases of the Rab and Arf families: Key regulators of intracellular trafficking in neurodegeneration. Int. J. Mol. Sci. 22.

Google Scholar

Balducci, C., Pierguidi, L., Persichetti, E., Parnetti, L., Sbaragli, M., Tassi, C., et al. (2007). Lysosomal hydrolases in cerebrospinal fluid from subjects with Parkinson’s disease. Mov. Disord. 22, 1481–1484. doi: 10.1002/mds.21399

PubMed Abstract | CrossRef Full Text | Google Scholar

Ballard, C., Gauthier, S., Corbett, A., Brayne, C., Aarsland, D., and Jones, E. (2011). Alzheimer’s disease. Lancet 377, 1019–1031. doi: 10.1016/S0140-6736(10)61349-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Berson, A., Nativio, R., Berger, S., and Bonini, N. (2018). Epigenetic regulation in neurodegenerative diseases. Trends Neurosci. 41, 587–598. doi: 10.1016/j.tins.2018.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Boccardi, V., Pelini, L., Ercolani, S., Ruggiero, C., and Mecocci, P. (2015). From cellular senescence to Alzheimer’s disease: The role of telomere shortening. Ageing Res. Rev. 22, 1–8. doi: 10.1016/j.arr.2015.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolognin, S., Lorenzetto, E., Diana, G., and Buffelli, M. (2014). The potential role of rho GTPases in Alzheimer’s disease pathogenesis. Mol. Neurobiol. 50, 406–422.

Google Scholar

Boo, J. H., Sohn, J. H., Kim, J. E., Song, H., and Mook-Jung, I. (2008). Rac1 changes the substrate specificity of gamma-secretase between amyloid precursor protein and Notch1. Biochem. Biophys. Res. Commun. 372, 913–917.

Google Scholar

Carapito, R., Ivanova, E., Morlon, A., Meng, L., Molitor, A., Erdmann, E., et al. (2019). ZMIZ1 variants cause a syndromic neurodevelopmental disorder. Am. J. Hum. Genet. 104, 319–330. doi: 10.1016/j.ajhg.2018.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, B. S., and Irizarry, R. A. (2010). A framework for oligonucleotide microarray preprocessing. Bioinformatics 26, 2363–2367. doi: 10.1093/bioinformatics/btq431

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatterjee, P., Roy, D., and Rathi, N. (2018). Epigenetic drug repositioning for Alzheimer’s disease based on epigenetic targets in human interactome. J. Alzheimers Dis. 61, 53–65. doi: 10.3233/JAD-161104

PubMed Abstract | CrossRef Full Text | Google Scholar

Chouliaras, L., Mastroeni, D., Delvaux, E., Grover, A., Kenis, G., Hof, P. R., et al. (2013). Consistent decrease in global DNA methylation and hydroxymethylation in the hippocampus of Alzheimer’s disease patients. Neurobiol. Aging 34, 2091–2099. doi: 10.1016/j.neurobiolaging.2013.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Coppieters, N., Dieriks, B., Lill, C., Faull, R., Curtis, M., and Dragunow, M. (2014). Global changes in DNA methylation and hydroxymethylation in Alzheimer’s disease human brain. Neurobiol. Aging 35, 1334–1344. doi: 10.1016/j.neurobiolaging.2013.11.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Daviglus, M. L., Bell, C. C., Berrettini, W., Bowen, P. E., Connolly, E. S. Jr., Cox, N. J., et al. (2010). National institutes of health state-of-the-science conference statement: Preventing Alzheimer disease and cognitive decline. Ann. Intern. Med. 153, 176–181. doi: 10.7326/0003-4819-153-3-201008030-00260

PubMed Abstract | CrossRef Full Text | Google Scholar

Disanza, A., Carlier, M., Stradal, T. E., Didry, D., Frittoli, E., Confalonieri, S., et al. (2004). Eps8 controls actin-based motility by capping the barbed ends of actin filaments. Nat. Cell Biol. 6, 1180–1188. doi: 10.1038/ncb1199

PubMed Abstract | CrossRef Full Text | Google Scholar

Federico, A., and Monti, S. (2020). hypeR: An R package for geneset enrichment workflows. Bioinformatics 36, 1307–1308. doi: 10.1093/bioinformatics/btz700

PubMed Abstract | CrossRef Full Text | Google Scholar

Forero, D., González-Giraldo, Y., López-Quintero, C., Castro-Vega, L., Barreto, G., and Perry, G. (2016). Meta-analysis of telomere length in Alzheimer’s disease. J. Gerontol. A Biol. Sci. Med. Sci. 71, 1069–1073. doi: 10.1093/gerona/glw053

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, A. I., Blace, N., Crary, J. F., Serrano, P. A., Leitges, M., Libien, J. M., et al. (2003). Protein kinase M zeta synthesis from a brain mRNA encoding an independent protein kinase C zeta catalytic domain. Implications for the molecular mechanism of memory. J. Biol. Chem. 278, 40305–40316. doi: 10.1074/jbc.M307065200

PubMed Abstract | CrossRef Full Text | Google Scholar

Huesa, G., Baltrons, M. A., Gomez-Ramos, P., Moran, A., Garcia, A., Hidalgo, J., et al. (2010). Altered distribution of RhoA in Alzheimer’s disease and AbetaPP overexpressing mice. J. Alzheimers Dis. 19, 37–56.

Google Scholar

Jager, P. L., Srivastava, G., Lunnon, K., Burgess, J., Schalkwyk, L. C., Yu, L., et al. (2014). Alzheimer’s disease: Early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci. Nat. Neurosci. 17, 1156–1163. doi: 10.1038/nn.3786

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, Y., Widschwendter, M., and Teschendorff, A. (2014). A systems-level integrative framework for genome-wide DNA methylation and gene expression data identifies differential gene expression modules under epigenetic control. Bioinformatics 30, 2360–2366. doi: 10.1093/bioinformatics/btu316

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordan, K. L., Koss, D. J., Outeiro, T. F., and Giorgini, F. (2022). Therapeutic targeting of Rab GTPases: Relevance for Alzheimer’s disease. Biomedicines 10:1141.

Google Scholar

Kass-Eisler, A., and Greider, C. (2000). Recombination in telomere-length maintenance. Trends Biochem. Sci. 25, 200–204. doi: 10.1016/S0968-0004(00)01557-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiral, F. R., Kohrs, F. E., Jin, E. J., and Hiesinger, P. R. (2018). Rab GTPases and membrane trafficking in neurodegeneration. Curr. Biol. 28, R471–R486.

Google Scholar

Kursa, M. B., and Rudnicki, W. R. (2010). “Feature selection with boruta package.”. J. Stat. Softw. 36, 1–13. doi: 10.18637/jss.v036.i11

CrossRef Full Text | Google Scholar

Lee, A. M., Kanter, B. R., Wang, D., Lim, J. P., Zou, M. E., Qiu, C., et al. (2013). Prkcz null mice show normal learning and memory. Nature 493, 416–419. doi: 10.1038/nature11803

PubMed Abstract | CrossRef Full Text | Google Scholar

Mastroeni, D., McKee, A., Grover, A., Rogers, J., and Coleman, P. (2009). Epigenetic differences in cortical neurons from a pair of monozygotic twins discordant for Alzheimer’s disease. PLoS One 4:e6617. doi: 10.1371/journal.pone.0006617

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathios, D., Hwang, T., Xia, Y., Phallen, J., Rui, Y., See, A. P., et al. (2019). Genome-wide investigation of intragenic DNA methylation identifies ZMIZ1 gene as a prognostic marker in glioblastoma and multiple cancer types. Int. J. Cancer 145, 3425–3435. doi: 10.1002/ijc.32587

PubMed Abstract | CrossRef Full Text | Google Scholar

Mayeux, R., and Stern, Y. (2012). Epidemiology of Alzheimer disease. Cold Spring Harb. Perspect. Med. 2:a006239. doi: 10.1101/cshperspect.a006239

PubMed Abstract | CrossRef Full Text | Google Scholar

Menna, E., Disanza, A., Cagnoli, C., Schenk, U., Gelsomino, G., Frittoli, E., et al. (2009). Eps8 regulates axonal filopodia in hippocampal neurons in response to brain-derived neurotrophic factor (BDNF). PLoS Biol 7:e1000138. doi: 10.1371/journal.pbio.1000138

PubMed Abstract | CrossRef Full Text | Google Scholar

Menna, E., Zambetti, S., Morini, R., Donzelli, A., Disanza, A., Calvigioni, D., et al. (2013). Eps8 controls dendritic spine density and synaptic plasticity through its actin-capping activity. EMBO J. 32, 1730–1744. doi: 10.1038/emboj.2013.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Offenhäuser, N., Borgonovo, A., Disanza, A., Romano, P., Ponzanelli, I., Iannolo, G., et al. (2004). The eps8 family of proteins links growth factor stimulation to actin reorganization generating functional redundancy in the Ras/Rac pathway. Mol. Biol. Cell 15, 91–98. doi: 10.1091/mbc.e03-06-0427

PubMed Abstract | CrossRef Full Text | Google Scholar

Offenhäuser, N., Castelletti, D., Mapelli, L., Soppo, B. E., Regondi, M. C., Rossi, P., et al. (2006). Increased ethanol resistance and consumption in Eps8 knockout mice correlates with altered actin dynamics. Cell 127, 213–226. doi: 10.1016/j.cell.2006.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Perkovic, M. N., Paska, A. V., Konjevod, M., Kouter, K., Strac, D. S., Erjavec, G. N., et al. (2021). Epigenetics of Alzheimer’s disease. Biomolecules 11:195. doi: 10.3390/biom11020195

PubMed Abstract | CrossRef Full Text | Google Scholar

Semick, S. A., Bharadwaj, R. A., Collado-Torres, L., Tao, R., Shin, J. H., Deep-Soboslay, A., et al. (2019). Integrated DNA methylation and gene expression profiling across multiple brain regions implicate novel genes in Alzheimer’s disease. Acta Neuropathol. 137, 557–569. doi: 10.1007/s00401-019-01966-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Serrano, P., Friedman, E. L., Kenney, J., Taubenfeld, S. M., Zimmerman, J. M., Hanna, J., et al. (2008). PKMzeta maintains spatial, instrumental, and classically conditioned long-term memories. PLoS Biol 6:2698–2706. doi: 10.1371/journal.pbio.0060318

PubMed Abstract | CrossRef Full Text | Google Scholar

Sindi, S., Solomon, A., Kåreholt, I., Hovatta, I., Antikainen, R., Hänninen, T., et al. (2021). Telomere length change in a multidomain lifestyle intervention to prevent cognitive decline: A randomized clinical trial. J. Gerontol. A Biol. Sci. Med. Sci. 76, 491–498. doi: 10.1093/gerona/glaa279

PubMed Abstract | CrossRef Full Text | Google Scholar

Tasegian, A., Paciotti, S., Ceccarini, M. R., Codini, M., Moors, T., Chiasserini, D., et al. (2017). Origin of α-mannosidase activity in CSF. Int. J. Biochem. Cell Biol. 87, 34–37. doi: 10.1016/j.biocel.2017.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Taub, P., Higginbotham, E., and Henry, R. (2013). Beneficial and detrimental effects of glycemic control on cardiovascular disease in type 2 diabetes. Curr. Cardiol. Rep. 15:332. doi: 10.1007/s11886-012-0332-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, Y., Morris, T. J., Webster, A. P., Yang, Z., Beck, S., Feber, A., et al. (2017). ChAMP: Updated methylation analysis pipeline for illumina beadchips. Bioinformatics 33, 3982–3984. doi: 10.1093/bioinformatics/btx513

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Huang, C., Lin, Y., Huang, W., Yang, C., Lee, C., et al. (2017). Conditional deletion of Eps8 reduces hippocampal synaptic plasticity and impairs cognitive function. Neuropharmacology 112, 113–123. doi: 10.1016/j.neuropharm.2016.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P. L., Niidome, T., Akaike, A., Kihara, T., and Sugimoto, H. (2009). Rac1 inhibition negatively regulates transcriptional activity of the amyloid precursor protein gene. J. Neurosci. Res. 87, 2105–2114.

Google Scholar

Wenk, G. (2003). Neuropathologic changes in Alzheimer’s disease. J. Clin. Psychiatry 64, 7–10.

Google Scholar

Zhang, Y., Hu, X., Li, Y., Zhao, C., Abudushalamu Jin, X., et al. (2021). International clinical practice guideline of chinese medicine alzheimer. World J. Tradit. Chin. Med. 7, 265–275.

Google Scholar

Zhu, X., Raina, A. K., Boux, H., Simmons, Z. L., Takeda, A., and Smith, M. A. (2000). Activation of oncogenic pathways in degenerating neurons in Alzheimer disease. Int. J. Dev. Neurosci. 18, 433–437.

Google Scholar

Keywords: Alzheimer’s disease, DNA methylation, CpG island, transcriptome, epigenetic age

Citation: Ren G, Song S, Zhang S-X, Liu Y, Lv Y, Wang Y-H, Zhao R and Li X-Y (2023) Brain region-specific genome-wide deoxyribonucleic acid methylation analysis in patients with Alzheimer’s disease. Front. Mol. Neurosci. 16:971565. doi: 10.3389/fnmol.2023.971565

Received: 17 June 2022; Accepted: 16 January 2023;
Published: 13 April 2023.

Edited by:

Michel J. Simonneau, École Normale Supérieure Paris-Saclay, France

Reviewed by:

Xiao Wang, Konge Larsen ApS, Denmark
Idil Yet, Hacettepe University, Türkiye

Copyright © 2023 Ren, Song, Zhang, Liu, Lv, Wang, Zhao and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xin-Yi Li, bGl4aW55aWx4eTkyQG91dGxvb2suY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.