Abstract
Background Dementia is a common disease influenced by both genetic and environmental factors. APOE ε4 is well-known to increase the risk of dementia, and it has been shown to attenuate the protective association of fish oil supplementation and the incidence of dementia. To identify more genetic factors with similar modifying effects, we performed a genome-wide scan.
Methods We first performed time-to-event genome-wide association study (GWAS) of all-cause dementia and two of its subtypes, Alzheimer’s disease (AD) and vascular dementia, in the UK Biobank. GWAS were performed in all participants (N = 357,631) and in two subgroups with or without fish oil supplementation (N = 113,267 and 244,364, respectively). Single nucleotide polymorphisms (SNPs) suggestively associated with dementia were then evaluated for their interactions with fish oil status in Cox-regression models. Furthermore, we conducted gene set enrichment analysis to identify the relevant cell types for these interaction signals.
Results Time-to-event GWAS identified 6, 5, and 2 genome-wide significant loci (p < 5e-8) for the incidence of all-cause dementia, AD, and vascular dementia, respectively. Most of them overlapped with previously known GWAS loci for AD and related dementia. A total of 178 suggestive GWAS loci (p < 1e-5) were passed onto interaction analysis, and 43 of them were found to significantly modify the association between fish oil supplementation and dementia incidence (p < 2.8e-4 with Bonferroni correction). One locus overlapped with a known AD GWAS locus (EED/PICALM) and two overlapped with GWAS loci for circulating omega-3 fatty acids (SRSF4, PSMG1). Gene set enrichment analysis found that candidate genes of interaction signals demonstrated tissue or cell-type specificity in the brain.
Conclusion We identified 43 genetic loci that modify the association between fish oil supplementation and dementia. These findings indicate a need for genome-informed personalized nutrition of fish oil supplementation for the purpose of dementia prevention.
Introduction
Dementia is a significant global health issue exacerbated by population aging and growth.1 This syndrome, characterized by the loss of cognitive functioning, is influenced by both genetic and environmental factors.2–4 Nutrition and dietary components are considered as the modifiable risk factors of dementia.2 Observational studies reported that omega-3 fatty acids were associated with a lower risk of cognitive decline.5 Habitual intake of fish oil supplements, a rich source of long-chain omega-3 fatty acids, was found to be inversely associated with the risk of all-cause dementia and several subtypes.6 The impacts of genetic factors have been examined by multiple genome-wide association studies (GWAS).4,7,8 A meta-study in 2023 revealed more than 75 loci associated with incident Alzheimer’s disease and related dementias (ADRD).4 One notable gene is APOE, with its risk allele APOE ε4 linked to vascular dysfunction, amyloid-β pathology, neurodegeneration, and ultimately, dementia.9–13
Gene-environment (G×E) interaction studies on human health are gaining widespread attention. These studies focus on how specific genetic and environmental factors synergistically increase disease risk, providing new insights into the etiology of complex diseases and promoting the development of preventative strategies and therapies.14 Recent studies showed that APOE ε4 modifies the protective associations of fish oil supplementation with all-cause dementia and vascular dementia, suggesting that only non-carriers of APOE ε4 may benefit from fish oil supplementation.15,16 Besides APOE, the potential modifying effects of other genetic factors remain unknown.
In this study, we sought to identify additional genetic variants that modify the association of fish oil supplementation with all-cause dementia and its two main subtypes, Alzheimer’s disease and vascular dementia. Time-to-event models used throughout the analysis allowed the detection of genetic variants associated with the onset of diseases.17 We adopted a two-step approach by first identifying suggestive loci associated with dementia phenotypes in the overall sample and the two fish oil supplementation subgroups. Then, time-to-event interaction analysis was used to identify significant interaction loci with Bonferroni correction. Finally, we assessed the expression signatures of candidate genes at these interaction loci and identified their biologically relevant tissue or cell types.
Methods
Study design and participants
We included participants of genetic European ancestry from UK Biobank, who completed the touch-screen questionnaire for the fish oil supplementation, did not have a diagnosis of dementia at baseline, and did not withdraw from the study by April 25th, 2023. The ancestry background was obtained from the Pan-UK Biobank Project by projecting the UKB individuals onto the principal components (PCs) inferred from the reference panels of the Human Genome Diversity Panel (HGDP) and 1000 Genomes Project.18
Outcomes
The dementia outcomes provided by UK Biobank were algorithmically defined using information from self-reported verbal interviews, linked hospital admissions, and death registries. Our study considered incident all-cause dementia, Alzheimer’s disease, and vascular dementia as the longitudinal outcomes. Incidence was defined as the first diagnosed date following the participant’s visit to the assessment center at baseline. The end of follow-up for each participant was recorded as the earliest of the following dates: dementia occurrence, death, or the end of follow-up records (Dec 13th, 2022).
Exposure
The touch-screen questionnaire completed at the assessment center provided information on fish oil supplementation. Participants were asked, “Do you regularly take any of the following?”. Participants who selected the option “Fish oil (including cod liver oil)” were considered fish oil users, and those who only selected other supplementation or “None of the above” were considered non-fish oil users. To evaluate the robustness of this questionnaire, we performed a sensitivity analysis on fish oil intake status based on the 24-hour dietary recall questionnaire. More details of sensitivity analysis were described in Supplementary Methods.
APOE ε4 dosage
APOE ε4 dosage was known to modify the relationship between fish oil supplementation and dementia onset.15 The APOE genotype was defined by two single nucleotide polymorphisms (SNPs), rs7412 and rs429358. The dosage of APOE ε4 was divided into 0 (ε1/ε1, ε1/ε2, ε2/ε2, ε2/ε3, ε3/ε3), 1 (ε1/ε4, ε3/ε4) and 2 (ε4/ε4), respectively. Since the determination of ε1/ε3 and ε2/ε4 from the genotype data was ambiguous, and ε2/ε4 may lead to conflated effect, we did not include them in the interaction analysis.
Genotype data processing
We used the imputed genotype dataset released by UK Biobank in 2017.19 About 92 million autosomal variants were imputed with the Haplotype Reference Consortium (HRC) panel and the UK10K + 1000 Genomes panel as references. Participants were removed if they had mismatches between self-reported sex and genetically inferred sex, exhibited sex chromosome aneuploidy, or were identified as outliers for heterozygosity or missing rate. We also excluded participants who had excessive third-degree relatives and randomly removed one individual among pairs of related individuals closer than the third degree.
For marker-based quality control, we first filtered out low-quality imputed variants with an info score below 0.3 and retained only biallelic variants. SNPs with a minor allele frequency (MAF) less than 0.1%, a missing calling rate larger than 5%, and a Hardy-Weinberg equilibrium P-value less than 10-8 were also removed. Finally, we ensured that the missing call rate per individual was below 2% after these processes. Genotype data was all processed by PLINK version 1.9 or version 2.0.20–22
Candidate SNPs selection
We performed time-to-event GWAS by a scalable and accurate method, saddlepoint approximation implementation based on the Cox proportional hazards regression model (SPACox), for all-cause dementia, Alzheimer’s disease, and vascular dementia, respectively.17 Age, sex, and the top 10 genetic PCs were included as covariates. Independent loci of genome-wide significant signals (p-value < 5e-8) for each outcome were defined based on clumping (p-value < 5e-8 and r² < 0.1 within a 250 kb physical distance) and then merging if two signal clusters are within 250 kb of each other. The novelty of these loci was defined with comparison to previous GWAS on the development of ADRD.4
Since we are looking for genetic factors that modify the association of fish oil supplementation and dementia outcomes, the associations of these genetic factors with the dementia outcomes are dependent on the fish oil supplementation status and may be masked or weakened in the GWAS in the entire dataset. Therefore, we also performed time-to-event GWAS of each outcome in the two subgroups of fish oil users and non-users, separately. Candidate variants for subsequent interaction analysis were selected using a lenient threshold of p-value < 1e-5 from GWAS in the entire dataset or in the two subgroups.
Time-to-event G×E interaction analysis
Interactions between candidate SNPs and fish oil supplementation on incident dementia were evaluated using the Cox proportional hazards models. Each model was adjusted by age, sex, the top 10 genetic PCs, fish oil supplementation, and the corresponding SNP. Interaction signals were identified with a Bonferroni-corrected threshold based on the number of independent candidate loci. These loci were obtained from the clumped regions of suggestive signals (r² < 0.1, p-value < 1e-5 within a 250 kb physical distance) with ±250 kb extension. The interaction analyses with APOE ε4 dosage were performed as positive controls. Furthermore, considering the potential effects of lifestyle, dietary patterns, and related medical history, we performed multiple sensitivity analyses by adjusting for additional covariates. Details were described in the Supplementary Methods.
Gene set enrichment analysis
Genome-wide significant interaction signals were annotated by the Variant Effect Predictor (VEP) toolset on the Ensembl platform.23 Each candidate variant was assigned location information for the genes, transcripts, and protein sequences that it might affect. Then, gene set enrichment analysis was conducted on the Web-based Cell-type-Specific Enrichment Analysis of Genes (WebCSEA) website for the list of candidate protein-coding genes included in the WebCSEA background gene list, to identify the tissue or cell-type-specific expression signatures of these interaction signals.24
Statistical analysis
Descriptive and inferential analyses were performed using R version 3.4.1.25 Imputation of missing covariate data was conducted by the mice package using predictive mean matching, logistic regression, and polytomous logistic regression methods, respectively, for continuous variables, binary variables, and unordered multiple categorial variables.26
Results
As shown in Figure 1, 425,157 participants from the UK Biobank are of genetic European ancestry, completed the touchscreen questionnaire on fish oil supplementation, and did not have a diagnosis of dementia at recruitment. After quality control with genotype data and removing related individuals, 357,631 participants with a mean age of 56.81±8.01 were included in our study. Around 31.7% of participants reported consuming fish oil supplements, as documented in the touchscreen questionnaire at baseline (Table S1-S2). Older individuals and women appeared more likely to take fish oil supplements.
Associations between fish oil supplementation and incident dementias
We evaluated the associations between fish oil supplementation and incident all-cause dementia and its subtypes using Cox proportional hazards models over a median follow-up period of 13.8 years. Regular intake of fish oil supplements demonstrated protective associations with incident all-cause dementia (HR, 0.92; 95 % CI, 0.87–0.96) and vascular dementia (HR, 0.80; 95 % CI, 0.72–0.89) after adjustment for age and sex (Table S3). However, fish oil supplementation showed no significant association with the onset of Alzheimer’s disease.
Time-to-event GWAS on dementias
We first performed time-to-event GWAS on three dementia outcomes separately, to identify disease-related loci and select candidate variants for further interaction analysis (Figure 2A). In the entire dataset of 357,631 participants and at the genome-wide significance (p-value <5e-8), we identified 6 loci for all-cause dementia, 5 for Alzheimer’s disease, and 2 for vascular dementia (Table S4). APOE was the locus shared among these three dementia outcomes. Bridging integrator 1 (BIN1), triggering receptor expressed on myeloid cells 2 (TREM2), and zinc finger CW-type (ZCWPW1)/PWWP domain containing 1 (NYAP1) were identified in both all-cause dementia and Alzheimer’s disease outcomes. Protein tyrosine kinase 2 beta (PTK2B)/clusterin (CLU) and membrane-spanning 4A (MS4A) were significantly associated with the onset of all-cause dementia, while aph-1 homolog B (APH1B) was specifically linked to Alzheimer’s disease. A novel locus for vascular dementia was discovered on chromosome 1, with castor zinc finger 1 (CASZ1) being the closest gene to the top signal.
Considering that genetic variants with interaction effects have exposure-dependent effects on the onset of dementia, we conducted stratified GWAS for the two subgroups of fish oil users and non-users. A total of 1,819 suggestive common SNPs (p-value<1e-5; MAF > 0.01) in 178 independent loci were then selected for interaction analysis (Figure S1-S3).
Time-to-event Gene × Fish Oil interaction analysis
We performed time-to-event gene × fish oil interaction analysis with the selected candidate SNPs for all-cause dementia, Alzheimer’s disease, and vascular dementia. Our primary interaction analyses used the Cox proportional hazards model, referred to as Model 1, which adjusted for age, sex, status of fish oil supplementation, and dosage of the corresponding SNP. A total of 81 SNPs in 43 loci were found to have significant interaction signals (p-value<2.8e-4, or 0.05/178 candidate loci) for one of the three dementia outcomes (Figure 2A, Table 1, Table S5-S8). A total of 22 protein-coding genes from these interaction loci were annotated through the VEP tool. One locus was annotated to the calcium/calmodulin-dependent protein kinase II delta (CAMK2D) gene, which is highly expressed in heart muscle tissue. The top variant at this locus (rs13123494) exhibited a pattern similar to that of APOE ε4, but with stronger effect sizes and statistical significance (Figure 2B and 2C). Among individuals with no or one C allele, taking fish oil supplements was associated with lower risk of all-cause dementia and vascular dementia, whereas a risk-increasing association was observed in those homozygous for the C allele. The association between AD and fish oil supplementation did not vary across different genotypes. Other interaction loci, such as rs2042528 and rs61804494, with different interaction patterns than that of APOE ε4, were also shown in Figure 2D and 2E.
We compared the GWAS loci and interaction loci. One locus on chromosome 11 overlapped with the reported ADRD locus embryonic ectoderm development (EED), and two loci on chromosomes 1 and 21 overlapped with known GWAS loci, serine and arginine rich splicing factor 4 (SRSF4) and proteasome assembly chaperone 1 (PSMG1), for polyunsaturated fatty acids (PUFAs) (Table S8). Of note, none of the interaction loci overlapped with the genome-wide significant loci in our time-to-event GWAS on dementias (Table S9).
While the total number of interaction loci across the three dementia outcomes was 43, the numbers were 17, 14, and 16, respectively, for all-cause dementia, Alzheimer’s disease, and vascular dementia. Among them, 15, 12, and 15 loci reached suggestive significance in GWAS in either fish oil subgroup. Interestingly, these interaction loci demonstrated a fish oil-specific genetic association pattern. 13 out of the 15 loci (87%, p-value = 0.0074) for all-cause dementia and 12 out of the 15 loci for vascular dementia (80%, p-value = 0.035) showing suggestive significance exclusively in the fish oil subgroup GWAS but not in the non-fish oil users. Alzheimer’s disease did not display such a preference, with 7 out of 12 loci (58%) exhibiting suggestive significance only in fish oil users (Table S9, Figure S4).
Sensitivity analyses
We performed a series of sensitivity analyses to evaluate the robustness of the identified interaction signals by adjusting for socioeconomic status, lifestyle, dietary patterns, and related medical history. Models adjusting for some or all of these additional covariates yielded similarly significant interaction signals (Figure S5). We also compared the interaction signals from Model 1 by using different questionnaires to define fish oil supplementation (Figure S6). Interaction signals on peroxisomal biogenesis factor 26 (PEX26) locus, EED locus, and CAMK2D locus were replicated in all-cause dementia, AD, and vascular dementia outcomes, respectively.
Gene set enrichment analysis
Gene set enrichment analysis was performed on 21 protein-coding genes included in the background gene list, to identify their expression signatures across 111 human tissue panels and 1,355 tissue-cell types spanning 12 human organ systems using the WebCSEA website. These 21 genes were highly expressed in various cell types of the nervous system (Figure 3). Notably, they were enriched in layer 6b excitatory neuron subtypes (Ex6b) in the frontal and visual cortex (p-value < 1.67e-4 and 4.64e-4, respectively), oligodendrocyte progenitor cells (OPCs) in the cerebellar hemisphere of adults (p-value < 5.39e-4), and oligodendrocytes in both the cerebrum and cerebellum (p-value < 1.49e-4 and 1.01e-3, respectively), as well as astrocytes in the cerebellum of fetuses (p-value < 2.56e-3) through permutation-based tests (Figure S7).
Discussion
In this study, we identified 6, 5, and 2 independent loci, respectively, associated with incident all-cause dementia, Alzheimer’s disease, and vascular dementia. Most of these loci were reported in the previous ADRD GWAS meta-analysis.4 Among 178 suggestive GWAS loci, 43 interaction signals significantly modified the association between fish oil supplementation and dementia development. Three loci were previously known to be associated with ADRD or PUFAs.4,27–31 However, none of these loci overlapped with the time-to-event GWAS loci we found, indicating that the interaction loci may have been missed in typical GWAS. Additionally, candidate protein-coding genes at interaction loci are highly expressed in the nervous system, especially in some specific neuron tissues and cell types. These results support that the new interaction signals may be involved in the pathway of brain functioning and PUFA metabolism, contributing to dementia development.
APOE were the strongest loci related to all-cause dementia and its two major subtypes in the time-to-event GWAS. It plays a key role in lipid metabolism and homeostasis. One of the major isoforms, APOE ε4 was known to be positively associated with amyloid plaques and neurofibrillary tangles (NFTs) formation, contributing to AD pathogenesis.32,33 Other known AD associated genes, BIN1, TREM2, ZCWPW1, and APH1B, were also replicated in our study.4 These genes were differentially expressed in specific cell types in the brains of AD patients compared to healthy individuals.34 Additionally, we identified a novel locus exclusively associated with vascular dementia. The top SNP within this locus is an intergenic variant likely involved in gene regulation. Castor zinc finger 1 (CASZ1), located in this region, plays roles in both neural and cardiovascular development.35 Other genes at this unknown locus, such as TAR DNA binding protein (TARDBP) and MBL associated serine protease 2 (MASP2), showed significant interactions with genes involved in tau-related pathologies.36
We identified 43 interaction loci modifying the associations between fish oil supplementation and dementia at statistical significance stronger than that of APOE ε4. These findings suggested a need and potential strategies for personalizing fish oil recommendations for high-risk populations. Among the interaction loci, EED, SRSF4, and PSMG1, were reported to be associated with incident ADD and PUFA traits, omega-3 PUFA level or percentage, and the ratio of omega-6 to omega-3 PUFA levels. These suggest that the interaction signals may be involved in dementia pathology and PUFA metabolism, although the specific mechanism is still unclear. Notably, none of the interaction loci overlapped with the time-to-event GWAS loci mentioned earlier. This disparity indicated that directly investigating G×E relationships based solely on GWAS hits may have missed these signals. Instead, we included candidate SNPs from subgroup-specific GWAS to consider the possibility of different or opposite effects of SNPs on outcomes between fish oil users and non-users, increasing the chances of catching the interaction signals. Previous interaction analyses between fish oil and polygenic risk score (PRS) of AD, excluding the effect of APOE ε4, showed no significant results in all-cause dementia and AD development.15 These also indicate the differences in the effects of GWAS loci and interaction loci in dementia and call for the consideration of G×E when applying PRS for disease prediction.
Protein-coding genes affected by interaction signals were enriched in high expression in oligodendrocytes, OPCs, and Ex6b cell types in adult brain tissues according to the gene set enrichment analysis. OPCs are known as oligodendrocyte precursor cells, which would divide and regenerate oligodendrocytes for myelin repair and myelination. This process slows down with aging and is associated with reduced myelination in AD, contributing to increased cognitive decline.37–39 The balance of excitatory and inhibitory inputs onto cortical neurons is strictly regulated by homeostatic mechanisms.40 A recent study revealed that the positive feedback loop between hyperexcitability and amyloidosis may lead to neuronal dysfunction and neurodegeneration.41 Additionally, omega-3 fatty acids have also been reported to protect against myelin sheath damage and affect synaptic plasticity.42,43 Therefore, we speculated that genetic signaling and fish oil supplementation may play some synergistic roles in myelin maintenance and synaptic homeostasis. Furthermore, we found that half of the interaction signals were mainly located in intronic or intergenic regions, suggesting gene regulation as their possible molecular mechanisms.
In this study, we mainly focused on all-cause dementia and its two subtypes, AD and vascular dementia, in the European population due to the limitation of sample and case sizes. Further studies on different populations and dementia subtypes, such as frontotemporal dementia, are needed to explore similarities and differences in dementia pathogenesis and across ancestral groups. Moreover, we applied a two-step approach to select a subset of SNPs for further time-to-event interaction analysis, reducing the computational cost of genome-wide interaction analysis of longitudinal outcomes. However, this approach was limited to categorial exposure and had the possibility of missing potential interaction signals during the candidate loci selection process. More effective algorithms need to be developed for broader and faster screening. Additionally, the binary fish oil exposure lacks information about dose or supplementation duration. Likely due to this imperfect exposure measurement, our results from sensitivity analyses with fish oil supplementation from different questionnaires only had limited replication. More reliable indicators of fish oil supplementation or dietary omega-3 intake, such as the circulating omega-3 fatty acid level, may be used for future studies as validation.
Conclusions
We identified 43 genomic loci that modify the association between fish oil supplementation and dementias. Gene set enrichment analysis showed that protein-coding genes around the interaction signals demonstrated high expression signatures in several neuron cell types in brain tissues. This study revealed potential interaction relationships between genetic factors and fish oil supplementation on the development of dementias. It provides new insight into the roles of fatty acids in disease pathogenesis and calls for genome-informed personalized dietary recommendations for dementia prevention.
Data Availability
The datasets analyzed during the current study are available from the UK Biobank through an application process.
Ethics Statement
The UK Biobank received ethical approval from the research ethics committee (reference ID: 11/ NW/0382). Written informed consent was obtained from all participants.
Data availability statement
The datasets analyzed during the current study are available from the UK Biobank through an application process (www.ukbiobank.ac.uk/).
Conflict of Interest
The authors declare no conflict of interest.
Funding information
Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institute of Health under the award number R35GM143060 (KY). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Supplementary materials
Table S1. Characteristics of the participants by fish oil supplementation
Table S2. Demographic and clinical characteristics of the participants by fish oil supplementation after imputation
Table S3. New locus and known loci related to all-cause dementia, Alzheimer’s disease, and vascular dementia
Table S4. The risk of 81 SNPs in 43 loci and APOE ε4 genotype interacted with fish oil supplementation in the development of all-cause dementia in the whole participants.
Table S5. The risk of 81 SNPs in 43 loci and APOE ε4 genotype interacted with fish oil supplementation in the development of Alzheimer’s disease in the whole participants.
Table S6. The risk of 81 SNPs in 43 loci and APOE ε4 genotype interacted with fish oil supplementation in the development of vascular dementia in the whole participants.
Table S7. 43 Loci interacted with fish oil supplementation in the development of dementias compared to previous discoveries.
Table S8. Significance of interaction loci among GWAS in the whole dataset and subgroup dataset.
Table S9. Association between interaction loci and dementia outcomes in the whole dataset and fish oil supplementation subgroups.
Figure S1. Manhattan plots of GWAS of all-cause dementia by fish oil supplementation
Figure S2. Manhattan plots of GWAS of Alzheimer’s disease by fish oil supplementation
Figure S3. Manhattan plots of GWAS of vascular dementia by fish oil supplementation
Figure S4. Comparison of the associations between interaction loci and dementia outcomes stratified by fish oil supplementation.
Figure S5. Interaction analysis of candidate SNPs with fish oil supplementation in the development of dementia using Model 1-4.
Figure S6. Interaction analysis of candidate SNPs with fish oil supplementation from the touchscreen questionnaire and 24-hour recall questionnaire in the development of dementia by model 1.
Figure S7. Protein-coding genes enrich in tissues by adult and fetal status.
Acknowledgements
This research has been conducted using the UK Biobank Resource under Application Number 48818. This work uses data provided by patients and collected by the NHS as part of their care and support. These data are copyrighted, 2022, NHS England. Reused with the permission of the NHS England and UK Biobank. All rights reserved. This research used data assets made available by National Safe Haven as part of the Data and Connectivity National Core Study, led by Health Data Research UK in partnership with the Office for National Statistics and funded by UK Research and Innovation (grant ref MC_PC_20058). We would like to express our sincere gratitude to the UK Biobank participants and administrative staff.