
==== Front
Commun Biol
Commun Biol
Communications Biology
2399-3642
Nature Publishing Group UK London

5725
10.1038/s42003-023-05725-x
Article
Differential methylation analysis in neuropathologically confirmed dementia with Lewy bodies
Reho Paolo 12
http://orcid.org/0000-0002-1524-9584
Saez-Atienzar Sara 3
Ruffo Paola 34
Solaiman Sultana 1
Shah Zalak 1
Chia Ruth 3
http://orcid.org/0000-0003-1218-4014
Kaivola Karri 1
http://orcid.org/0000-0003-0527-2446
Traynor Bryan J. 3
Tilley Bension S. 5
http://orcid.org/0000-0002-7131-4643
Gentleman Steve M. 5
http://orcid.org/0000-0002-5565-6678
Hodges Angela K. 6
http://orcid.org/0000-0001-6314-216X
Aarsland Dag 67
Monuki Edwin S. 8
Newell Kathy L. 9
Woltjer Randy 10
Albert Marilyn S. 11
Dawson Ted M. 11121314
Rosenthal Liana S. 11
Troncoso Juan C. 15
http://orcid.org/0000-0001-6516-4944
Pletnikova Olga 16
Serrano Geidy E. 17
Beach Thomas G. 17
Easwaran Hariharan P. 18
http://orcid.org/0000-0002-6623-0429
Scholz Sonja W. sonja.scholz@nih.gov

111
1 https://ror.org/01s5ya894 grid.416870.c 0000 0001 2177 357X Neurodegenerative Diseases Research Unit, National Institute of Neurological Disorders and Stroke, Bethesda, MD USA
2 https://ror.org/00hj8s172 grid.21729.3f 0000 0004 1936 8729 Laboratory of Precision Environmental Health, Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY USA
3 https://ror.org/049v75w11 grid.419475.a 0000 0000 9372 4913 Neuromuscular Diseases Research Section, National Institute on Aging, Bethesda, MD USA
4 https://ror.org/02rc97e94 grid.7778.f 0000 0004 1937 0319 Medical Genetics Laboratory, Department of Pharmacy, Health and Nutritional Sciences, University of Calabria, Rende, Italy
5 https://ror.org/041kmwe10 grid.7445.2 0000 0001 2113 8111 Neuropathology Unit, Department of Brain Sciences, Imperial College London, London, UK
6 https://ror.org/0220mzb33 grid.13097.3c 0000 0001 2322 6764 Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, UK
7 https://ror.org/04zn72g03 grid.412835.9 0000 0004 0627 2891 Centre for Age-Related Medicine, Stavanger University Hospital, Stavanger, Norway
8 grid.266093.8 0000 0001 0668 7243 Department of Pathology & Laboratory Medicine, School of Medicine, University of California Irvine, Irvine, CA USA
9 grid.257413.6 0000 0001 2287 3919 Department of Pathology and Laboratory Medicine, Indiana University School of Medicine, Indianapolis, IN USA
10 https://ror.org/009avj582 grid.5288.7 0000 0000 9758 5690 Department of Neurology, Oregon Health & Sciences University, Portland, OR USA
11 https://ror.org/00za53h95 grid.21107.35 0000 0001 2171 9311 Department of Neurology, Johns Hopkins University Medical Center, Baltimore, MD USA
12 grid.21107.35 0000 0001 2171 9311 Neuroregeneration and Stem Cell Programs, Institute of Cell Engineering, Johns Hopkins University School of Medicine, Baltimore, MD USA
13 grid.21107.35 0000 0001 2171 9311 Department of Pharmacology and Molecular Sciences, Johns Hopkins University School of Medicine, Baltimore, MD USA
14 grid.21107.35 0000 0001 2171 9311 Solomon H. Snyder Department of Neuroscience, Johns Hopkins University School of Medicine, Baltimore, MD USA
15 https://ror.org/00za53h95 grid.21107.35 0000 0001 2171 9311 Department of Pathology (Neuropathology), Johns Hopkins University Medical Center, Baltimore, MD USA
16 https://ror.org/01y64my43 grid.273335.3 0000 0004 1936 9887 Department of Pathology and Anatomical Sciences, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, Buffalo, NY USA
17 https://ror.org/04gjkkf30 grid.414208.b 0000 0004 0619 8759 Civin Laboratory for Neuropathology, Banner Sun Health Research Institute, Sun City, AZ USA
18 grid.280502.d 0000 0000 8741 3625 The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Baltimore, MD USA
5 1 2024
5 1 2024
2024
7 3513 1 2023
19 12 2023
© This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Dementia with Lewy bodies (DLB) is a common form of dementia in the elderly population. We performed genome-wide DNA methylation mapping of cerebellar tissue from pathologically confirmed DLB cases and controls to study the epigenetic profile of this understudied disease. After quality control filtering, 728,197 CpG-sites in 278 cases and 172 controls were available for the analysis. We undertook an epigenome-wide association study, which found a differential methylation signature in DLB cases. Our analysis identified seven differentially methylated probes and three regions associated with DLB. The most significant CpGs were located in ARSB (cg16086807), LINC00173 (cg18800161), and MGRN1 (cg16250093). Functional enrichment evaluations found widespread epigenetic dysregulation in genes associated with neuron-to-neuron synapse, postsynaptic specialization, postsynaptic density, and CTCF-mediated synaptic plasticity. In conclusion, our study highlights the potential importance of epigenetic alterations in the pathogenesis of DLB and provides insights into the modified genes, regions and pathways that may guide therapeutic developments.

An epigenome-wide association study of dementia with Lewy bodies discovers differentially methylated probes, genomic regions and pathways associated with this disease.

Subject terms

DNA methylation
Dementia
https://doi.org/10.13039/100000002 U.S. Department of Health & Human Services | National Institutes of Health (NIH) ZIANS003154 Scholz Sonja W. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Dementia with Lewy bodies (DLB) is a heterogeneous neurodegenerative disease characterized by parkinsonism, visual hallucinations, fluctuating mental status, and REM-sleep behavior disorder1. There are an estimated 1.4 million cases living in the United States2, and current therapy is limited to symptomatic and supportive care. Although genetic research studies have identified heritable factors that are important in the etiology of this understudied disease3,4, the molecular causes remain poorly understood, and little is known about non-genetic contributors to its pathogenesis.

Epigenetic changes are modifications to the DNA that regulate gene expression. These modifications are influenced by aging, the environment, lifestyle, disease state, and other factors. They allow cells to respond dynamically to the outside world and are considered the interface between genetic and environmental components. One type of epigenetic alteration is DNA methylation. This tissue-specific mechanism occurs when a methyl group is transferred onto a cytosine, mostly occurring in the context of a cytosine-phosphate-guanine dinucleotide sequence (CpG) in higher eukaryotes. These modifications change the DNA accessibility to the transcriptional machinery complex and fine-modulate gene expression. Importantly, epigenetic modifications are thought to play a prominent role in age-associated neurological diseases, such as Alzheimer’s disease and Parkinson’s disease5–7. Evidence is also emerging that DNA methylation changes influence the risk of developing DLB4,8–10.

In this study, we investigated the role of DNA methylation in the pathogenesis of DLB. We performed an epigenome-wide association study (EWAS) to characterize the differential methylation patterns in cerebellar tissue obtained from 298 pathologically confirmed DLB cases and 203 neurologically healthy controls. We further performed gene-region and pathway analyses, demonstrating widespread epigenetic changes associated with this fatal neurodegenerative disease.

Results

Epigenome-wide association study design

We performed an EWAS using cerebellar brain tissue obtained from patients diagnosed with DLB and healthy individuals. We profiled their DNA methylation status with the Illumina MethylationEPIC arrays. After quality control filters were applied, 728,197 sites were tested for association with DLB in 278 cases and 172 controls. We then designed the regression model adjusting the data by age, sex, experimental batch, five principal components from genotyping data, cell type proportion, and 44 surrogate variables from methylation data.

EWAS identifies differentially methylated probes

We identified seven differentially methylated probes (DMPs) that surpassed the genome-wide significance threshold (Bonferroni adjusted p value < 0.05) (Fig. 1, Table 1). Among these, four probes were hypomethylated (cg16086807, cg04866173, cg11099930, and cg24435966) and three were hypermethylated (cg18800161, cg16250093, and cg06951630) in DLB patients compared to healthy controls (Fig. 2). Overall, we detected mild epigenetic modulations, where the differences between β-values (Δβ) in cases and controls ranged from −0.10 to 0.091 (Fig. 1). Six out of the seven DMPs were mapped within gene regions, overlapping the gene body of ARSB (cg16086807, Bonferroni-adjusted p value = 1.19E-03, Δβ = −0.020), LINC00173 (cg18800161, adjusted p value = 5.26E-03, Δβ = 0.006), MGRN1 (cg16250093, adjusted p value = 1.02E-02, Δβ = 0.008), FHL2 (cg04866173, adjusted p value = 1.15E-02, Δβ = −0.023), IQSEC1 (cg06951630, adjusted p value = 3.62E-02, Δβ = 0.009), and the promoter region of NIPBL (cg11099930, adjusted p value = 1.50E-02, Δβ = −0.010) (Table 1).Fig. 1 Volcano plot, QQ-plot, and Manhattan plot of EWAS results.

Volcano plot (a) showing statistical significance (-log10 p value) and magnitude of change (delta beta, Δβ) of all CpG sites included in the DLB EWAS analysis. Red and blue dots indicate significantly hypomethylated DMPs and hypermethylated sites in DLB cases compared to controls, respectively. The Bonferroni adjusted p value < 0.05 threshold is shown as a red dashed line. QQ-plot (b) showing the p value distribution and inflation (lambda value). Density plot (c) illustrating the observed p value distribution. Manhattan plot (d) showing the p values of the tested probes across the genome. The genome-wide significance threshold (Bonferroni adjusted p value < 0.05) is shown as a red dashed line, while the orange dashed line represents the FDR threshold. Probes surpassing the genome-wide significance are shown as red dots.

Table 1 Significant differentially methylated probes in the DLB EWAS.

Probe ID	Chr	Position	Gene	Gene Region	Δβ	P value	Adj.P	Δβ Males vs Females	P value Males vs Females	Adj.P.Males vs Females	
cg16086807	5	78,271,020	ARSB	Body	−0.020	1.63E-09	1.19E-03	−0.037	1.02E-09	7.39E-04	
cg18800161	12	116,971,933	LINC00173	Body	0.006	7.22E-09	5.26E-03	0.007	2.432E-08	1.77E-02	
cg16250093	16	4,691,998	MGRN1	Body	0.008	1.41E-08	1.02E-02	0.008	2.34E-07	1.71E-01	
cg04866173	2	105,990,524	FHL2	Body	−0.023	1.58E-08	1.15E-02	−0.023	2.98E-07	2.17E-01	
cg11099930	5	36,876,680	NIPBL	TSS200	−0.010	2.06E-08	1.50E-02	−0.011	6.71E-08	4.88E-02	
cg06951630	3	13,249,045	IQSEC1	Body	0.009	4.97E-08	3.62E-02	0.009	3.94E-07	2.87E-01	
cg24435966	4	43,342,431	–	–	−0.029	5.33E-08	3.89E-02	−0.037	6.47E-10	4.71E-04	
Chr. chromosome, Adj.P adjusted p value

Differentially methylated probes in the DLB EWAS. Chromosome positions are shown relative to the human reference genome (hg19). Gene names are shown according to UCSC RefGen. The Δβ values refer to the difference between DNA methylation (β-values) in cases compared to controls (e.g., −0.020 indicates that DLB cases show a 2% decrease in DNA methylation compared to controls). Adjusted p value refer to Bonferroni corrected p values.

Fig. 2 Differentially methylated probes in DLB.

The violin plots show the DNA methylation (beta value) distribution (violin shape) in the seven differentially methylated probes (panels a–g). The vertical axis represents the range of values in the dataset, where 0 and 1 mean fully unmethylated and fully methylated respectively. The box plot represents the interquartile range of the dataset (25% bottom, 75% top), the middle line represents the median of the distribution, and the central line shows the value distribution. Black dots represent outliers. The overall difference between DLB cases and controls is shown as delta beta (Δβ); negative and positive values refer to hypomethylation and hypermethylation in DLB cases, respectively (e.g., −0.020 indicates that DLB cases show a 2% decrease in DNA methylation compared to controls). P values refer to Bonferroni corrected p values.

Among the DMP-associated genes, three have been previously implicated in neurodegenerative or neurodevelopmental disorders: ARSB, encoding for the Arylsulfatase B, has been associated with Alzheimer’s disease and Parkinson’s disease;11,12 MGRN1, encoding Mahogunin Ring Finger 1, has been implicated in late-onset spongiform neurodegeneration;13 and IQSEC1, encoding for the IQ Motif And Sec7 Domain ArfGEF 1, has been associated with a neurodevelopmental disorder14.

Our data showed modest inflation (lambda = 1.19). To further explore inflation, we also corrected the p values using the Bioconductor package bacon (lambda = 1.02). This approach replicated the results of the main study, where all seven DMPs reached the FDR threshold and the top two probes surpassed the genome-wide significance threshold based on Bonferroni correction (Supplementary Fig. 1).

Genes related to neurodegenerative diseases identified among the sub-significant signals

In addition to the top seven DMPs, our analysis identified 41 probes surpassing the less conservative false discovery rate (FDR) threshold of 0.05 (Supplementary Data 1). These sub-significant probes were localized within genes previously associated with neurological disorders: EPS8 (cg06658698, FDR adjusted p value = 1.49E-02, Δβ = 0.030), ARL6IP1 (cg20872370, FDR adjusted p value = 2.24E-02, Δβ = −0.007), FRMD4A (cg03775372, FDR adjusted p value = 2.82E-02, Δβ = −0.020), and PIAS1 (cg05751215, FDR adjusted p value = 3.40E-02, Δβ = −0.009). In particular, EPS8, encoding for the Epidermal Growth Factor Receptor Pathway Substrate 8, and FRMD4A, that encodes for FERM Domain Containing 4 A, have been recently associated with Alzheimer’s disease; 15,16 ARL6IP1, encoding for ADP Ribosylation Factor Like GTPase 6 Interacting Protein 1, has been associated with hereditary spastic paraplegia;17 and PIAS1, that encodes for Protein Inhibitor Of Activated STAT 1, has been implicated in Huntington’s disease18.

Identification of differentially methylated genomic regions

Recent research has shown that differentially methylated regions (DMRs) are more highly associated with diseases than differential methylation at individual CpG sites alone19. For this reason, we examined DMRs in our case-control dataset. Our analysis identified 32 CpG sites that clustered in three different DMRs characterized by a hypermethylation signature in DLB cases compared to controls (Table 2). For example, we identified a DMR overlapping the promoter region of DHRS4 and its antisense lncRNA DHRS4-AS1 (adjusted p value = 1.59E-10, mean Δβ = 0.022). DHRS4 has been recently suggested as a novel risk gene for inducing neurodegeneration in mouse models of amyotrophic lateral sclerosis20.Table 2 Differentially methylated regions in the DLB EWAS.

Chr.	Start	End	Width (bp)	No. CpGs	P value	Mean Δβ	Overlapping Genes	Gene Function	
6	29,648,161	29,649,024	864	20	5.85E-11	0.048	ZFP57	Regulation of gene expression controlling DNA methylation	
14	24,422,749	24,422,956	208	6	1.59E-10	0.022	DHRS4, DHRS4-AS1	Cell metabolism	
1	183,154,778	183,155,154	377	6	4.37E-10	0.018	LAMC2	Cell adhesion, motor neuron axon guidance, epidermis development	
Chr. chromosome, No.CpGs number of CpGs per region, Mean Δβ mean probe delta beta.

Chromosome positions are shown relative to the human reference genome (hg19). P value refers to smoothed FDR corrected p value.

Sex-specific EWAS delineates a male-driven effect

To explore differential methylation patterns among men and women, we performed an interaction model comparing male study participants (n = 165 DLB cases and n = 122 controls) and female subjects (n = 113 DLB cases and n = 50 controls). We identified five differentially methylated probes showing a statistically significant epigenetic modulation in males compared to females (Bonferroni adjusted p value < 0.05) (Fig. 3, Supplementary Data 2). All five DMPs were also identified in the overall cohort analysis. The most associated site in males was located 461 base pairs upstream of RP11-1E6.1 (cg24435966, Bonferroni-adjusted p value = 4.71E-04, Δβ = −0.037), a predicted lncRNA with unknown function. The four male-specific DMPs overlapped the following genes: ARSB (cg16086807, adjusted p value = 7.39E-04, Δβ = −0.022), FRMD4A (cg03775372, adjusted p value = 1.17E-02, Δβ = −0.026), LINC00173 (cg18800161, adjusted p value = 1.77E-02, Δβ = 0.007), and NIPBL (cg11099930, adjusted p value = 4.88E-02, Δβ = −0.011). Finally, our sex-specific EWAS did not identify any significant differentially methylated regions.Fig. 3 Volcano plot, QQ-plot, and Manhattan plot of sex-specific EWAS.

Volcano plot a) showing genome-wide significance (-log10 p value) and magnitude of change (Δβ) of all sites included in the analysis. Negative values (dots on the left side of the volcano plot) indicate hypomethylated DMPs, and hypermethylated sites are displayed as positive values (right-sided dots). The Bonferroni adjusted p value < 0.05 threshold is shown as a red dashed line. QQ-plot (b) showing p values distribution and inflation (lambda value). Density plot (c) illustrating the observed p value distribution. Manhattan plot (d) showing the p value of the probes across the genome. The Bonferroni adjusted p value < 0.05 threshold is shown as a red dashed line. Probes surpassing the genome-wide significance are shown as red dots.

Functional enrichment analysis determined pathways associated with DLB

We performed a pathway enrichment analysis of differentially methylated genes in our case-control cohort. Specifically, we investigated the 43 genes that overlapped FDR-significant differentially methylated probes. This analysis detected eight biological processes, cellular components, and molecular functions that were significantly associated with DLB. Among the Gene Ontology (GO) terms, we identified an association with the terms “regulation of response to stimulus” (Bonferroni-adjusted p value = 0.0237), “neuron to neuron synapse” (adjusted p value = 0.0494), “vesicle” (adjusted p value = 0.0489), “postsynaptic specialization” (adjusted p value = 0.0385), and “postsynaptic density” (adjusted p value = 0.0276) (Fig. 4). Of note, 29 out of 44 genes showed an interaction with the CTCF protein, a transcription factor that has been associated with Alzheimer’s disease and synaptic organization (adjusted p value = 0.0011) (Table 3)21,22.Fig. 4 Pathway enrichment analysis in pathologically confirmed DLB.

Functional enrichment of significant gene ontology pathways for biological processes (GO:BP, orange dot), cellular components (GO:CC, blue dots), and TRANSFAC database (TF, red dots) in pathologically confirmed DLB versus controls. The x-axis shows the p value associated with each pathway on a -log10 scale and the size of each dot indicates the number of genes involved.

Table 3 Biological pathways associated with differential methylation in DLB.

ID	Pathway	Source	No. Genes	P value	
GO:0048583	regulation of response to stimulus	GO:BP	18	4.57E-03	
GO:0051716	cellular response to stimulus	GO:BP	20	3.11E-02	
GO:0014069	postsynaptic density	GO:CC	5	2.74E-02	
GO:0031253	cell projection membrane	GO:CC	5	3.27E-02	
GO:0032279	asymmetric synapse	GO:CC	5	3.14E-02	
GO:0099572	postsynaptic specialization	GO:CC	5	3.82E-02	
GO:0031982	vesicle	GO:CC	6	4.85E-02	
GO:0098984	neuron to neuron synapse	GO:CC	5	4.90E-02	
TF:M12593	Factor: CTCF; motif: NNRSYGCCMCCTGSTGGCCN	TF	29	1.11E-03	
TF:M00918	Factor: E2F; motif: TTTSGCGSG	TF	14	3.49E-02	
No. Genes number of genes, GO Gene Ontology, BP Biological Process, CC Cellular Component, TF TRANSFACT.

P value refers to Bonferroni corrected p value.

Discussion

Our analyses illustrate the value of EWAS in unraveling the multiplex architectures of neurodegenerative conditions and highlight the potential contributions of methylation to the pathogenesis of DLB. We identified several probes that were differentially methylated in the DLB cases (Fig. 1). Interestingly, many of these genes are highly expressed in the brain (Fig. S2) and have been previously implicated in neurological disorders or central nervous system development.

Chief among the loci identified by our study was ARSB, encoding for Arylsulfatase B, a member of the sulfatase family that removes sulfate groups from chondroitin-4-sulfate, triggering its degradation. This lysosomal enzyme is involved in cell adhesion, migration, and invasion in colonic epithelium23. In cultured astrocytes, ARSB silencing increases chondroitin-4-sulfate and neurocan levels, and inhibits astrocyte-mediated neurite outgrowth, suggesting that ARSB may play an important role in neuronal plasticity in the central nervous system24. Homozygous or compound heterozygous mutations of the gene lead to Mucopolysaccharidosis type VI, a lysosomal storage disorder characterized by skeletal anomalies, short stature, and cardiac abnormalities25. Our data complement previous studies implicating ARSB variants in Alzheimer’s disease and Parkinson’s disease11,12, expanding the pathogenic role of the gene in neurodegeneration.

DLB affects males more than females26, but little is known about sex-specific contributions to the pathogenesis. Our sex-specific EWAS showed a male-driven effect, suggesting that epigenetic modifications could modulate the risk in males and females separately. However, our analysis did not include sex chromosomes, and future studies exploring epigenetic changes in these chromosomes are needed to draw conclusions.

Interestingly, our study also showed an enrichment of genes that regulate the neuron-to-neuron synapse, postsynaptic specialization, and postsynaptic density in DLB. Moreover, we highlighted the transcriptional repressor CTCF as a key factor able to orchestrate these biological processes. This protein is actively involved in maintaining the three-dimensional structure of the chromatin, creating topologically associated functional domains within the nucleus27. Our data suggest that the chromatin architecture could play a role in the pathogenesis of DLB, and support recent evidence implicating CTCF-mediated synaptic plasticity in Alzheimer’s disease22. These observations expand the role of CTCF and synaptic organization-related genes in neurodegeneration.

Our data corroborate emerging evidence implicating aberrant epigenetic modulation in neurodegeneration28,29. Only a few studies have investigated the epigenetic changes associated with DLB, and they differ in sample size, targeted tissue, and study design8,10,30. For example, Shao and colleagues performed an EWAS of the Brodman area 7 of the brain in a cohort consisting of fifteen pathologically confirmed DLB cases and sixteen neurologically healthy controls. Despite the limited sample size and the different brain regions investigated, their study design does represent the closest structure to our EWAS. In contrast, Nasamran and collaborators profiled blood epigenetic modulations comparing 42 DLB patients and 50 Parkinson’s disease dementia cases, while Pihlstrom and colleagues explored the epigenetic modulations associated with different Braak Lewy body disease stages in 322 Parkinson’s disease and DLB cases. Our study identified none of the DMPs or DMRs previously associated with DLB (Supplementary Data 3). However, the outlined differences in the design of these studies may account for these discrepancies.

Significant changes in methylation have also been observed in Alzheimer’s disease and Parkinson’s disease31–33. Sharma and collaborators showed that epigenetic modifications across single nucleotide polymorphisms, located within the first intron of the SNCA gene, modulate the susceptibility to Parkinson’s disease. A meta-analysis of 1,453 individuals with Alzheimer’s disease, investigating the epigenetic changes associated with Braak neurofibrillary tangle stage, identified differentially methylated sites and regions in the prefrontal cortex, temporal gyrus, and entorhinal cortex, but not in the cerebellum. Interestingly, our EWAS replicated some of the findings from Smith’s study, identifying epigenetic modification involving a common CpG site within FRMD4A (cg03775372) and a common DMP-associated gene, MKL2. This finding raises the scientific interest surrounding cg03775372 and the FRMD4A gene, since the same CpG reached the FDR significance in our main study and achieved genome-wide significance in the sex-specific EWAS.

We selected cerebellar tissue for our research as it is relatively spared in the terminal stages of DLB, unlike cortical tissue in which most cells of interest have been lost to neurodegeneration. Selecting a relatively spared tissue source provides a more accurate window into the epigenetic plasticity of a disease. This detail needs to be weighed against the fact that the disease-relevant changes are likely more prominent and representative in the regions primarily affected by the disease. Sampling multiple regions for comparison would have been ideal. Future efforts will likely increase our ability to identify disease-associated methylation patterns across the brain, and we have made the EWAS results from our study publicly available to facilitate this unfolding research.

Our EWAS has several limitations. The Illumina MethylationEPIC array contains only a fraction of the human genome’s CpGs. Furthermore, the bisulfite conversion chemistry does not distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), a brain-specific intermediate product of 5mC demethylation34. We also did not adjust the data for environmental factors that may impact the methylation status, such as vascular comorbidities, smoking, and alcohol use.

The minfi-approach that we used to estimate the cell type proportion was able to discriminate between neuronal and non-neuronal cells, using the frontal cortex region as reference. This method is based on the neuronal-specific protein NeuN, expressed in the vast majority of neurons, though Purkinje cells represent an exception. Our approach has already been successfully applied in EWAS based on DNA from cerebellar tissue35,36. Moreover, exploring the surrogate variables, we identified a mild negative correlation between NeuN-negative cell proportion and the surrogate variables 1 and 3 (Pearson correlation = −0.64 and −0.61, respectively), and a positive correlation between NeuN-positive cell proportion and surrogate variable 3 (Pearson correlation = 0.83). These data suggest that the surrogate variable analysis accounts for at least some of the variability due to the different cell types within a tissue. As such, it represents a valuable approach that could be employed in similar instances where there is no tool to estimate the cell proportion in a tissue (Supplementary Fig. 3).

Genome-wide association studies are affected by inflation and EWAS are not an exception, so exploring inflation is crucial to reduce the number of unreliable results. Our data showed moderate inflation (lambda = 1.19) (Fig. 1b). We further reduced genomic inflation using the R package bacon (lambda = 1.02, Supplementary Fig. 1), and we were able to replicate the results of the main EWAS. To consolidate our main results, we also performed the EWAS using the MLM-based omic association (MOA) tool from OSCA37. This tool represents a stringent approach to processing DNA methylation data. Not surprisingly, therefore, inflation was drastically reduced when this tool was applied (lambda = 1.005), but none of the probes surpassed the genome-wide significance threshold (Supplementary Fig. 4). However, the top 48 CpG sites identified in our study showed a consistent pattern when comparing the two approaches (Pearson correlation = 0.87789), and the directions of their modulations were coherent (Supplementary Fig. 5).

Sample size differences between the male- and the female-specific cohorts may represent a bias of the study, particularly for the sex-specific EWAS. Furthermore, the absence of matching for age and sex between the patients and controls might have led to an overestimation of the contribution of the identified loci in the pathogenesis of DLB. However, this effect was probably mitigated by including age and sex as covariates in the association model38. The use of post-mortem tissues cannot discriminate between causal effects and the downstream consequences of the DNA methylation changes observed. Finally, we detected only modest effects on DNA methylation, though this is consistent with previous EWAS efforts investigating neurodegenerative disorders8,10,30.

In summary, we investigated the differential methylation signature of a large cohort of patients with pathologically confirmed DLB. We delineated clear epigenetic modulation associated with this common form of neurodegeneration, defined differentially methylated probes and regions, and highlighted new loci and biological pathways affected by these changes. In particular, we provide evidence implicating the ARSB gene and the CTCF-mediated synaptic plasticity to DLB. Our study underlines the potential role that epigenetic modulation plays in the pathogenesis of DLB and represents an opportunity for the future identification of biomarkers and new therapeutic targets.

Methods

Study samples

Frozen cerebellar tissue from 298 DLB patients and 203 neurologically healthy controls were obtained from brain donation programs. The demographic and clinical characteristics of the study participants are summarized in Supplementary Table 1. The DLB patients were diagnosed with pathologically definite disease (limbic or neocortical subtype) according to the McKeith consensus criteria1. Neurologically healthy controls were selected based on the absence of neurological disease in their clinical history and the absence of neurodegenerative disease on pathological examination. All participants were of European ancestry. Informed consent for post mortem brain tissue donation was obtained from all subjects or their surrogate decision makers according to the Declaration of Helsinki. Each brain donation program was approved by its own institutional ethics committee. These convenience control samples were obtained from the same brain banks as the DLB cases and were of European ancestry. The controls were not specifically matched for age or sex, however, age and sex distributions among cases and controls were comparable (Supplementary Fig. 6, Supplementary Table 1).

SNP genotyping, quality control, and data processing

Genomic DNA was extracted from cerebellar tissue samples using the Maxwell RSC Tissue kit following the manufacturer’s instructions (Promega, Madison, WI, USA). Each DNA sample was processed on the Infinium Global Diversity Array + Neuro Booster chip (v.1.0, Illumina, San Diego, CA, USA). This array genotypes over 1.8 million single-nucleotide polymorphisms (SNPs) across the genome, including 75,000 genetic variants previously implicated in common neurodegenerative diseases. Quality control checks were performed using PLINK (v.1.9.0-beta 4.4)39. Samples were excluded from the analysis based on the following criteria: (1) genome-wide call rate < 0.98, (2) heterozygosity outliers (> ± 0.15 F-statistic), (3) mismatch between reported sex and genotypic sex, (4) non-European ancestry (based on principal components analysis when compared to the HapMap 3 Genome Reference Panel), and (5) duplicates and related samples (pi-hat > 0.125). For variant-level quality control, we excluded variants for the following reasons: (1) non-random missingness between cases and controls (p value ≤ 1×10-4), (2) haplotype-based non-random missingness (p value ≤ 1×10-6), (3) deviation from Hardy-Weinberg equilibrium test in controls (mid-p value ≤ 1×10-6), and (4) overall missingness rate of ≥ 5%.

DNA methylation profiling, quality control, and data processing

For each sample, 450 ng of genomic DNA was bisulfite-converted (Zymo Research, Irvine, CA, USA) and processed on Infinium HumanMethylation EPIC BeadChips (v.1.0, Illumina), according to the manufacturer’s protocol. This array measures DNA methylation signals of 863,904 CpGs across the genome. Raw intensity data were processed in R (v.4.0.5) using the MethylAid (v.1.24.0)40, minfi (v.1.36.0)41, wateRmelon (v.1.34.0)42, and maxprobes (v.0.0.2, https://github.com/markgene/maxprobes) packages.

Quality checks and filtering were performed to exclude low-quality samples and probes (summarized in Fig. S7). Samples meeting the following criteria were excluded: (1) low overall quality based on sample-dependent and sample-independent control probes and methylated/unmethylated probes ratio (using default settings in MethylAid), (2) bisulfite conversion rate < 80.0%, (3) mismatch between reported sex and genotypic sex, (4) mean detection p value > 0.01, (5) samples that were flagged as outliers using the default setting in wateRmelon, and (6) samples that had > 1.0% of probes with a detection p value > 0.05.

We then performed a principal component analysis of the Infinium MethylationEPIC array control probes. Data normalization was done using the R package minfi (preprocessFunnorm function), with 24 principal components (explaining 99.0% of variance)43. Following normalization, we excluded probes that had: (1) detection p value < 0.01, (2) CpG sites containing SNPs of any minor allele frequency, (3) probes located on sex chromosomes, and (4) any Illumina methylation array 450 K and EPIC850K cross-reactive probes44.

We generated beta- and M-values using the minfi functions getBeta and getM, respectively. Significant surrogate variables have been generated from M-values through the sva package (v.3.46.0) in R45. We estimated the proportion of neuronal and non-neuronal cell types using the minfi package (estimateCellCounts function)35,36.

EWAS analysis in DLB

The DNA methylation status of each CpG site (as measured by ß-values, Fig. S8) was tested for association with DLB using the Bioconductor package limma (v.3.46.0)42 in R (v.4.0.5). In contrast to GWAS, there are no common guidelines established in the field for conducting EWAS analyses. We performed our study using beta values since they directly represent the proportion of methylated CpG sites, making it easier to relate the values to biological processes. Supporting this approach, beta values and M-values showed highly correlated outputs in our analysis (Pearson correlation of p values = 0.904), and the top 48 probes showed a consistent direction of effect between the two approaches (Fig. S9). Age, sex, experimental batch, the first five principal components (generated from the Infinium Global Diversity Array + Neuro Booster genotyping data to account for population stratification), NeuN-positive/NeuN-negative cell type proportion (minfi), and all significant surrogate variables (n = 44) (Bioconductor package sva, v.4.3) were included as covariates in the linear regression model. A two-sided p value and a Δβ-value was calculated for each CpG site. We used the Illumina EPIC annotation R package (IlluminaHumanMethylationEPICanno.ilm10b4.hg19) to define the overlapping genes. Four additional genes were manually annotated since the probes overlapped a gene (IQSEC1: cg06951630, LINC01158:cg07171538, and LINC00856:cg05757757) or mapped within 10 base-pairs upstream of a gene (EPHA8:cg25394625). P values were adjusted by Bonferroni correction. The Bonferroni threshold for declaring a differentially methylated probe (DMP) to be genome-wide significant was 6.87 × 10-8 ( = 0.05/728,197 markers). We identified sub-significant probes surpassing the False Discovery Rate threshold of 0.05, and we included them in the functional enrichment analysis. To further reduce genomic inflation removing unknown bias, we corrected the p values using the Bioconductor package bacon (v.1.26.0). EWAS evaluations were also performed comparing male (n = 165 cases and 122 controls) and female participants (n = 113 cases and 50 controls) through an interaction model, to assess possible sex-specific epigenetic modulation. Sex was not included as a covariate in those analyses.

The Bioconductor package DMRcate (v.2.4.1) was used with the recommended default settings (lambda = 1000 and C = 2, corresponding to 1 standard deviation of Gaussian kernel each 500 base pairs) to identify and evaluate regions in the DLB data for evidence of differential methylation43. This software calculated p values based on smoothed FDR of CpGs within the region.

Pathway enrichment analysis

Functional enrichment analysis was performed using the g:Profiler toolkit (v.0.7.0)46. We investigated the following pathways: 1) Gene Ontology (GO) biological processes (n = 15,808), 2) GO cellular components (n = 1973), 3) GO molecular functions (n = 5015)47,48, 4) pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG; n = 563)49–51, 5) pathways described in Reactome (n = 2532)52, 6) WikiPathways (n = 790)53, 7) transcription factors in the TRANSFAC database (n = 11,647)54, 8) microRNAs in miRTarBase (n = 2658)55, 9) proteins in the Human Protein Atlas (n = 830)56, 10) protein complexes in CORUM (n = 2885)57, and 11) traits in the Human Phenotype Ontology database (n = 10,668)58. Pathways with fewer than five genes and a p value > 0.05 were excluded from the analysis. Bonferroni corrections were applied to the p values in each pathway to correct for multiple testing.

Statistics and reproducibility

We performed a case-control association study by fitting a linear regression model for each marker using the Bioconductor package limma (v.3.46.0). The topTable function was used to calculate the statistics of differentially methylated probes comparing DLB cases to healthy control subjects, adjusting the p values for multiple testing. Bonferroni-corrected genome-wide significance threshold was set to p < 6.87×10-8 ( = 0.05/728,197 sites tested). We applied a False Discovery Rate p value correction to declare sub-significant markers. To facilitate reproducible results, we made the analysis code publicly available on Github (https://github.com/pireho/EWAS-Lewy_body_dementia) and https://zenodo.org (DOI: 10.5281/zenodo.10365334).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Peer review file

Supplementary Information

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Reporting summary

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-023-05725-x.

Acknowledgements

We thank the patients and families whose help and participation made this work possible. We are grateful to the Banner Sun Health Research Institute Brain and Body Donation Program of Sun City, Arizona, for the provision human biological materials. The Brain and Body Donation Program is supported by the National Institute of Neurological Disorders and Stroke (U24 NS072026 National Brain and Tissue Resource for Parkinson’s Disease and Related Disorders), the National Institute on Aging (P30 AG19610 and P30 AG072980, Arizona Alzheimer’s Disease Center), the Arizona Department of Health Services (contract 211002, Arizona Alzheimer’s Research Center), the Arizona Biomedical Research Commission (contracts 4001, 0011, 05-901, and 1001 to the Arizona Parkinson’s Disease Consortium) and the Michael J. Fox Foundation for Parkinson’s Research. The study used tissue samples and data from the Johns Hopkins Morris K. Udall Center of Excellence for Parkinson’s Disease Research (NIH P50 NS38377). The Parkinson’s UK Brain Bank at Imperial College London is funded by Parkinson’s UK, a charity registered in England and Wales (948776) and in Scotland (SC037554). We thank the members of the Laboratory of Neurogenetics (NIH) for their collegial support and technical assistance. We thank the NABEC Consortium, the Virginia Commonwealth University Brain Bank, and the NIH NeuroBioBank for providing tissue samples. Tissue/data used in this research was obtained from the Human Brain Collection Core, Intramural Research Program, NIMH (http://www.nimh.nih.gov/hbcc). This work utilized the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov). This study was supported by the Intramural Research Program of the National Institutes of Health (National Institute of Neurological Disorders and Stroke; project number: ZIANS003154). K.K. was funded by the Finnish Cultural Foundation, The Finnish Parkinson Foundation, The Päivikki and Sakari Sohlberg Foundation, and The Finnish Brain Foundation. This study was supported by the Intramural Research Program of the National Institutes of Health (National Institute of Neurological Disorders and Stroke; project number: ZIANS003154). K.K. was funded by the Finnish Cultural Foundation, The Finnish Parkinson Foundation, The Päivikki and Sakari Sohlberg Foundation, and The Finnish Brain Foundation.

Author contributions

Conceptualization: P.R., S.S.-A., B.J.T., H.P.E., S.W.S.; Formal analysis: P.R., S.S.-A., P.Ru., Z.H., R.C., K.K.; Investigation: P.R., S.S.-A., P.Ru., S.S., Z.S., R.C., K.K.; Resources: B.S.T., S.M.G, A.K.H., D.A., E.S.M., K.L.N., R.W., M.S.A., T.M.D., L.S.R., J.C.T., O.P., G.E.S., T.G.B.; Data curation: P.R., S.S.-A., P.Ru., S.S., Z.S., R.C., K.K.; Writing – original draft: P.R.; Writing – review & editing: All authors; Visualization: P.R., S.S.-A., P.Ru.; Supervision, project administration and funding: S.W.S.

Peer review

Peer review information

Communications Biology thanks Moshe Szyf, Lasse Pihlstrom and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Melanie Bahlo, Karli Montague-Cardoso and George Inglis. A peer review file is available.

Data availability

The EWAS summary statistics have been deposited into the EWAS catalog (www.ewascatalog.org). Individual-level methylation array data are available on dbGaP (accession #: phs001963). Source data underlying Figs. 2 and 4 are provided in Supplementary Data 4.

Code availability

Analyses were performed using open-source tools and the code is available on GitHub at https://github.com/pireho/EWAS-Lewy_body_dementia and https://zenodo.org (10.5281/zenodo.10365334)59.

Competing interests

The authors declare the following competing interests: S.W.S. serves on the Scientific Advisory Council of the Lewy Body Dementia Association and the Multiple System Atrophy Coalition. S.W.S. and B.J.T. receive research support from Cerevel Therapeutics. S.W.S. serves on the editorial board of the journal JAMA Neurology. B.J.T. is associate editor for the journal Brain. B.J.T. holds patents on the clinical testing and therapeutic implication of the C9orf72 repeat expansion. All other authors have no conflicts of interest to declare that are relevant to the content of this article.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. McKeith IG Diagnosis and management of dementia with Lewy bodies: Fourth consensus report of the DLB Consortium Neurology 2017 89 88 100 10.1212/WNL.0000000000004058 28592453
2. Association L. B. D. 10 things you should know about LBD.).
3. Chia R Genome sequencing analysis identifies new loci associated with Lewy body dementia and provides insights into its genetic architecture Nat. Genet 2021 53 294 303 10.1038/s41588-021-00785-3 33589841
4. Kaivola K Genome-wide structural variant analysis identifies risk loci for non-Alzheimer’s dementias Cell Genom. 2023 3 100316 10.1016/j.xgen.2023.100316 37388914
5. Lardenoije R The epigenetics of aging and neurodegeneration Prog. Neurobiol. 2015 131 21 64 10.1016/j.pneurobio.2015.05.002 26072273
6. Guhathakurta S Bok E Evangelista BA Kim YS Deregulation of alpha-synuclein in Parkinson’s disease: Insight from epigenetic structure and transcriptional regulation of SNCA Prog. Neurobiol. 2017 154 21 36 10.1016/j.pneurobio.2017.04.004 28445713
7. Gasparoni G DNA methylation analysis on purified neurons and glia dissects age and Alzheimer’s disease-specific changes in the human cortex Epigenetics Chromatin 2018 11 41 10.1186/s13072-018-0211-3 30045751
8. Pihlstrom L Epigenome-wide association study of human frontal cortex identifies differential methylation in Lewy body pathology Nat. Commun. 2022 13 4932 10.1038/s41467-022-32619-z 35995800
9. Beyer K Domingo-Sabat M Ariza A Molecular pathology of Lewy body diseases Int J. Mol. Sci. 2009 10 724 745 10.3390/ijms10030724 19399218
10. Nasamran CA Differential blood DNA methylation across Lewy body dementias Alzheimers Dement (Amst.) 2021 13 e12156 10.1002/dad2.12156 33665346
11. Potkin SG Hippocampal atrophy as a quantitative trait in a genome-wide association study identifying novel susceptibility genes for Alzheimer’s disease PLoS One 2009 4 e6501 10.1371/journal.pone.0006501 19668339
12. Jansen IE Discovery and functional prioritization of Parkinson’s disease candidate genes from large-scale whole exome sequencing Genome Biol. 2017 18 22 10.1186/s13059-017-1147-9 28137300
13. Benvegnù S Aging triggers cytoplasmic depletion and nuclear translocation of the E3 ligase mahogunin: a function for ubiquitin in neuronal survival Mol. Cell 2017 66 358 372 10.1016/j.molcel.2017.04.005 28475871
14. Ansar M Bi-allelic variants in IQSEC1 cause intellectual disability, developmental delay, and short stature Am. J. Hum. Genet 2019 105 907 920 10.1016/j.ajhg.2019.09.013 31607425
15. Lancour D Analysis of brain region-specific co-expression networks reveals clustering of established and novel genes associated with Alzheimer disease Alzheimers Res. Ther. 2020 12 103 10.1186/s13195-020-00674-7 32878640
16. Lambert JC Genome-wide haplotype association study identifies the FRMD4A gene as a risk locus for Alzheimer’s disease Mol. Psychiatry 2012 18 461 470 10.1038/mp.2012.14 22430674
17. Novarino G Exome sequencing links corticospinal motor neuron disease to common neurodegenerative disorders Science 2014 343 506 511 10.1126/science.1247363 24482476
18. Ochaba J PIAS1 regulates mutant huntingtin accumulation and huntington’s disease-associated phenotypes in vivo Neuron 2016 90 507 520 10.1016/j.neuron.2016.03.016 27146268
19. Eckhardt F DNA methylation profiling of human chromosomes 6, 20 and 22 Nat. Genet 2006 38 1378 1385 10.1038/ng1909 17072317
20. Li S Identification of molecular correlations between DHRS4 and progressive neurodegeneration in amyotrophic lateral sclerosis by gene co-expression network analysis Front Immunol. 2022 13 874978 10.3389/fimmu.2022.874978 35479082
21. Vostrov AA Quitschke WW The zinc finger protein CTCF binds to the APBbeta domain of the amyloid beta-protein precursor promoter. Evidence for a role in transcriptional activation J. Biol. Chem. 1997 272 33353 33359 10.1074/jbc.272.52.33353 9407128
22. Patel PJ Ren Y Yan Z Epigenomic analysis of Alzheimer’s disease brains reveals diminished CTCF binding on genes involved in synaptic organization Neurobiol. Dis. 2023 184 106192 10.1016/j.nbd.2023.106192 37302762
23. Bhattacharyya S Tobacman JK Arylsulfatase B regulates colonic epithelial cell migration by effects on MMP9 expression and RhoA activation Clin. Exp. Metastasis 2009 26 535 545 10.1007/s10585-009-9253-z 19306108
24. Zhang X Arylsulfatase B modulates neurite outgrowth via astrocyte chondroitin-4-sulfate: dysregulation by ethanol Glia 2014 62 259 271 10.1002/glia.22604 24311516
25. Azevedo AC Clinical and biochemical study of 28 patients with mucopolysaccharidosis type VI Clin. Genet 2004 66 208 213 10.1111/j.1399-0004.2004.00277.x 15324318
26. Chiu SY Sex differences in dementia with Lewy bodies: Focused review of available evidence and future directions Parkinsonism Relat. Disord. 2023 107 105285 10.1016/j.parkreldis.2023.105285 36682958
27. Ong CT Corces VG CTCF: an architectural protein bridging genome topology and function Nat. Rev. Genet 2014 15 234 246 10.1038/nrg3663 24614316
28. Murthy M Cheng YY Holton JL Bettencourt C Neurodegenerative movement disorders: An epigenetics perspective and promise for the future Neuropathol. Appl Neurobiol. 2021 47 897 909 10.1111/nan.12757 34318515
29. Maloney B Lahiri DK Epigenetics of dementia: understanding the disease as a transformation rather than a state Lancet Neurol. 2016 15 760 774 10.1016/S1474-4422(16)00065-X 27302240
30. Shao X Dementia with Lewy bodies post-mortem brains reveal differentially methylated CpG sites with biomarker potential Commun. Biol. 2022 5 1279 10.1038/s42003-022-03965-x 36418427
31. Smith RG A meta-analysis of epigenome-wide association studies in Alzheimer’s disease highlights novel differentially methylated loci across cortex Nat. Commun. 2021 12 3517 10.1038/s41467-021-23243-4 34112773
32. Sharma A Common genetic variants associated with Parkinson’s disease display widespread signature of epigenetic plasticity Sci. Rep. 2019 9 18464 10.1038/s41598-019-54865-w 31804560
33. Moore K McKnight AJ Craig D O’Neill F Epigenome-wide association study for Parkinson’s disease Neuromol. Med. 2014 16 845 855 10.1007/s12017-014-8332-8
34. Huang Y The behaviour of 5-hydroxymethylcytosine in bisulfite sequencing PLoS One 2010 5 e8888 10.1371/journal.pone.0008888 20126651
35. Guintivano J Aryee MJ Kaminsky ZA A cell epigenotype specific model for the correction of brain cellular heterogeneity bias and its application to age, brain region and major depression Epigenetics 2013 8 290 302 10.4161/epi.23924 23426267
36. Pellegrini C A meta-analysis of brain DNA methylation across sex, age, and alzheimer’s disease points for accelerated epigenetic aging in neurodegeneration Front Aging Neurosci. 2021 13 639428 10.3389/fnagi.2021.639428 33790779
37. Zhang F OSCA: a tool for omic-data-based complex trait analysis Genome Biol. 2019 28 20
38. Luca D On the use of general control samples for genome-wide association studies genetic matching highlights causal variants Am. J. Hum. Genet 2008 82 453 463 10.1016/j.ajhg.2007.11.003 18252225
39. Purcell S PLINK: a tool set for whole-genome association and population-based linkage analyses Am. J. Hum. Genet 2007 81 559 575 10.1086/519795 17701901
40. van Iterson M MethylAid: visual and interactive quality control of large Illumina 450k datasets Bioinformatics 2014 30 3435 3437 10.1093/bioinformatics/btu566 25147358
41. Aryee MJ Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays Bioinformatics 2014 30 1363 1369 10.1093/bioinformatics/btu049 24478339
42. Pidsley R A data-driven approach to preprocessing Illumina 450K methylation array data BMC Genomics 2013 14 293 10.1186/1471-2164-14-293 23631413
43. Lehne B A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies Genome Biol. 2015 16 37 10.1186/s13059-015-0600-x 25853392
44. Chen YA R. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray Epigenetics 2013 8 203 209 10.4161/epi.23470 23314698
45. Leek JT Johnson WE Parker HS Jaffe AE Storey JD The sva package for removing batch effects and other unwanted variation in high-throughput experiments Bioinformatics 2012 28 882 823 10.1093/bioinformatics/bts034 22257669
46. Raudvere U g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update) Nucleic Acids Res 2019 47 191 198 10.1093/nar/gkz369
47. Ashburner M Gene Ontology: tool for the unification of biology Nat. Genet 2000 25 25 29 10.1038/75556 10802651
48. The Gene Ontology Consortium. The Gene Ontology knowledgebase in 2023 Genetics 2023 224 iyad031 10.1093/genetics/iyad031 36866529
49. Kanehisa M Goto S KEGG: kyoto encyclopedia of genes and genomes Nucleic Acids Res 2000 28 27 30 10.1093/nar/28.1.27 10592173
50. Kanehisa M Toward understanding the origin and evolution of cellular organisms Protein Sci. 2019 28 1947 1951 10.1002/pro.3715 31441146
51. Kanehisa M Furumichi M Sato Y Kawashima M Ishiguro-Watanabe M KEGG for taxonomy-based analysis of pathways and genomes Nucleic Acids Res. 2023 51 D587 D592 10.1093/nar/gkac963 36300620
52. Fabregat A Reactome pathway analysis: a high-performance in-memory approach BMC Bioinform. 2017 18 142 10.1186/s12859-017-1559-2
53. Pico AR WikiPathways: pathway editing for the people PLoS Biol. 2008 6 e184 10.1371/journal.pbio.0060184 18651794
54. Matys V TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes Nuc Ac Res 2006 34 D108 D110 10.1093/nar/gkj143
55. Hsi-Yuan H miRTarBase update 2022: an informative resource for experimentally validated miRNA–target interactions Nucleic Acids Res. 2022 50 D222 D230 10.1093/nar/gkab1079 34850920
56. Uhlén M Proteomics. Tissue-based map of the human proteome Science 2015 23 347
57. Ruepp A CORUM: the comprehensive resource of mammalian protein complexes Nucleic Acids Res. 2008 36 D646 D650 10.1093/nar/gkm936 17965090
58. Köhler S The Human Phenotype Ontology in 2021 Nucleic Acids Res. 2021 49 D1207 D1217 10.1093/nar/gkaa1043 33264411
59. Paolo, R. Zenodo repository “pireho/EWAS-Lewy_body_dementia: EWAS in Lewy body dementia”. https://zenodo.org/records/10365335 (2023).
