- 1Biodemography of Aging Research Unit, Social Science Research Institute, Duke University, Durham, NC, United States
- 2Danish Aging Research Center, Department of Public Health, University of Southern Denmark, Odense, Denmark
- 3Division of Statistical Genomics, Department of Genetics, Washington University School of Medicine, St. Louis, MO, United States
- 4Division of Endocrinology, Diabetes and Nutrition and Program for Personalized and Genomic Medicine, University of Maryland School of Medicine, Baltimore, MD, United States
- 5Duke Center for the Study of Aging and Human Development, Duke University, Durham, NC, United States
- 6Durham VA Geriatrics Research Education and Clinical Center, Durham, NC, United States
Dysregulation of physiological processes may contribute to Alzheimer’s disease (AD) development. We previously found that an increase in the level of physiological dysregulation (PD) in the aging body is associated with declining resilience and robustness to major diseases. Also, our genome-wide association study found that genes associated with the age-related increase in PD frequently represented pathways implicated in axon guidance and synaptic function, which in turn were linked to AD and related traits (e.g., amyloid, tau, neurodegeneration) in the literature. Here, we tested the hypothesis that genes involved in PD and axon guidance/synapse function may jointly influence onset of AD. We assessed the impact of interactions between SNPs in such genes on AD onset in the Long Life Family Study and sought to replicate the findings in the Health and Retirement Study. We found significant interactions between SNPs in the UNC5C and CNTN6, and PLXNA4 and EPHB2 genes that influenced AD onset in both datasets. Associations with individual SNPs were not statistically significant. Our findings, thus, support a major role of genetic interactions in the heterogeneity of AD and suggest the joint contribution of genes involved in PD and axon guidance/synapse function (essential for the maintenance of complex neural networks) to AD development.
1 Introduction
Age is the strongest risk factor for Alzheimer’s disease (AD), suggesting that aging-related changes in the body may contribute to AD development. Indeed, human aging is characterized by an increasing level of dysregulation across body systems, leading to declines in resilience and robustness, with increased susceptibility to multiple health disorders and reduced ability to recover, as people age (Cohen et al., 2013; Arbeev et al., 2016; Whitson et al., 2016; Arbeev et al., 2019; Dansereau et al., 2019; Müller et al., 2019; Chanti-Ketterl et al., 2021; Ukraintseva et al., 2021; Lu et al., 2022). AD has also been linked to the dysregulation of metabolic processes and neural networks (An et al., 2018; de la Monte et al., 2019; Jeong et al., 2019; Yan et al., 2020; Fagiani et al., 2021; Jouini et al., 2021; Nanclares et al., 2021; Joshi et al., 2022; Martínez-Serra et al., 2022), which suggests that the aging-related increase in physiological dysregulation may contribute to the development of AD (Tan et al., 2011; von Bernhardi et al., 2015; Lutshumba et al., 2021; Piehl et al., 2022).
In a number of recent works, we explored a composite index of physiological dysregulation (PD) based on the statistical Mahalanobis distance that measures deviations of multiple biomarkers from their baseline/normal values, thereby reducing the high-dimensional biomarker space into a single PD estimate (Cohen et al., 2014; Arbeev et al., 2016). We found that the rate of increase in PD is a significant predictor of worsening health status manifested by onset of major complex diseases such as cancer, cardiovascular diseases or type 2 diabetes (Arbeev et al., 2019). We also found that genes associated with the increase in PD in our genome-wide association study (GWAS) of the Long Life Family Study participants are significantly enriched for involvement in axon guidance and synaptic function (Arbeev et al., 2020). Notably, pathways and molecules regulating axon guidance and synaptic function have also been broadly connected to AD and related traits (e.g., amyloid, tau, neurodegeneration) in the literature (Cissé et al., 2011; Jun et al., 2014 Wetzel-Smith et al., 2014; Jiang et al., 2015; Giri et al., 2016; Hashimoto et al., 2016; Kang et al., 2016; Zuko et al., 2016; Oguro-Ando et al., 2017; Han et al., 2018; Li et al., 2018; Liu et al., 2018; Lee et al., 2019; Whelan et al., 2019; Yang et al., 2019; Bamford et al., 2020; Cuestas Torres and Cardenas, 2020; Chen et al., 2021; Zhang et al., 2021). These findings indicate a possibility that genes involved in PD and axon guidance/synaptic function may jointly contribute to onset of AD via intertwined mechanisms. In this paper, we explore whether such genes influence AD onset in older adults, individually or as a result of epistatic interactions, using longitudinal human data from the Long Life Family Study and the Health and Retirement Study.
2 Materials and Methods
2.1 Data
The Long Life Family Study (LLFS) is a family-based, longitudinal study of healthy aging and longevity that enrolled participants at three United States field centers (Boston, New York, Pittsburgh) and in Denmark. From 2006 to 2009, more than 4,900 participants from 583 families selected for exceptional familial longevity based on the Family Longevity Selection Score (Sebastiani et al., 2009) were enrolled in the study. The second in-person visit of surviving participants from visit 1 and new enrollees was completed in 2014–2017. The third in-person visit started in 2020 and is ongoing (these data are not yet available for this study). At enrollment, the participants provided information on socio-demographic indicators, past and current medical conditions, medication use, and physical and cognitive functioning (Wojczynski et al., 2022). LLFS participants were followed up annually to track their vital and health status. The analyses reported in this paper used the July 2020 release of LLFS data, with the latest recorded follow-up date on 22 July 2019. Ages at the baseline visit were validated using dates of birth from official documents (such as a birth certificate or driver’s license) (Elo et al., 2013) in the United States and through the civil registration system in Denmark. Ages at censoring for those who did not die within the follow-up period were determined from dates of birth and the last follow-up. Information on past and present health conditions (including AD or dementia) was collected during the baseline and follow-up interviews from study participants or proxies (when the participant was unable to provide an answer). The question asked at the baseline visit was, “Please respond “yes” or “no” if you have EVER been told by a doctor that you had this condition.” Similar questions were asked during follow-up interviews (“Please respond “yes” or “no” if you have EVER been told by a doctor that you had this condition since we last interviewed you on … ”). Participants or their proxies were then asked the age at which they were diagnosed or their best guess. Also, during follow-up interviews participants who had cognitive concern were given a dementia questionnaire adopted by the study and were asked if they saw a neurologist and “If “Yes” was the diagnosis Alzheimer’s Disease?” Age at AD or dementia diagnosis (approximating the age of AD onset) was determined by the age they responded affirmatively and, if this information was unavailable, their age when the questionnaire was administered. Ages at censoring, computed as described above, were used in the time-to-event modeling for those who did not report the event during the follow-up.
Written informed consent was obtained from all subjects following protocols approved in the United States by the respective field center’s Institutional Review Boards (IRBs) and, in Denmark, by the Regional Committees on Health Research Ethics for Southern Denmark. In this paper, we performed secondary analyses of LLFS data collected at all field centers. This study was approved by the Duke University Health System IRB.
The University of Michigan’s Health and Retirement Study (HRS) is a longitudinal panel study that surveys a representative sample of approximately 20,000 people in America over the age of 50 every 2 years. The study launched data collection in 1992 and has re-interviewed the original sample of respondents every 2 years. The participants have provided information about income, work, assets, pension plans, health insurance, disability, physical health and functioning, cognitive functioning, and healthcare expenditures. Data were collected from different sources, including in-person, telephone, mail, and the Internet. The target population for the original HRS cohort included all adults in the contiguous United States born during the years 1931–1941 who reside in households, with a 2:1 oversample of African-American and Hispanic populations. The original sample is refreshed with new birth cohorts (51–56 years of age) every 6 years. The sample has been expanded over the years to include a broader range of birth cohorts and described in more detail by Sonnega et al. (2014). The RAND HRS Longitudinal File is an easy-to-use dataset based on the HRS core data. This file was developed at RAND with funding from the National Institute on Aging and the Social Security Administration. This file (version 2016 V1) was used in this study. Entry ages were determined from the first exam date (called “wave” in HRS). Exit ages were determined from ages at death and dates of birth for individuals who died during the follow-up. For individuals who survived the whole study, the lifespan was censored on 31 December 2017 (the latest available date in the sample used). Starting from wave 10 (year 2010), respondents (or proxies, if respondents were not able to provide answers) were asked if they have ever been told they have Alzheimer’s disease or dementia. If the condition was reported at a prior interview, then the question was skipped. Ages of AD onset were determined from the dates of the interview at the respective wave when the condition was first reported and dates of birth. For individuals who did not report the condition, ages at censoring, computed as described above, were used in the time-to-event modeling.
These two datasets are well-suited to serve as discovery and validation sets for two reasons. First, they are complementary because the LLFS is a family-based study in which participants were selected from longevity-enriched families presumably sharing genetic factors likely related to biological aging. Because the HRS is a nationally-representative population-based study, it can be used to determine whether relations detected in LLFS are replicated in the general United States population. Second, the studies use the same genotyping platform (see Section 2.3), thus simplifying the direct validation of discovered associations. As the LLFS sample has predominantly (>99%) White participants, we selected a subsample of White HRS participants for replication (see Section 2.3).
2.2 Candidate genes
Table 1 shows the list of candidate genes selected for this study, with the number of SNPs in each gene that remained after the genetic quality control (QC, see Section 2.3). This list includes six genes found in our earlier GWAS (TRIO, ALCAM, CADM1, CNTN6, RTN4, PLXNA4) that were significantly associated with the rate of increase in PD with age and enriched in biological processes of axon guidance/synaptic function (Arbeev et al., 2020). Since AD has been linked to physiological dysregulation and synaptic functioning in other research (Martínez-Serra et al., 2022; Ma et al., 2023; see also Introduction), the above genes were selected as candidates for the association analyses with AD risk in this study. We also included in the analysis six additional genes (PLXNB2, EPHA2, EPHB2, SEMA4D, UNC5B, UNC5C) for plexins, semaphorins, ephrins, and receptors to netrins that are major molecular players in axon guidance (during development) and synaptic function (in adults) and have been repeatedly linked to AD and related traits, such as amyloid beta, tau phosphorylation, neurodegeneration, cognitive function and neuroinflammation, in the literature. These additional genes were used as the literature-based candidate AD genes that are involved in axon guidance/synaptic function (Tang et al., 2008; Cissé et al., 2011; Jun et al., 2014; Wetzel-Smith et al., 2014; Jiang et al., 2015; Giri et al., 2016; Hashimoto et al., 2016; Zuko et al., 2016; Han et al., 2018; Li et al., 2018; Liu et al., 2018; McDermott et al., 2018; Lee et al., 2019; Yang et al., 2019; Chen et al., 2021; Zhang et al., 2021; Evans et al., 2022; Ma et al., 2022). The resulting set of twelve candidate genes was, therefore, highly enriched for putative roles in axon guidance, synaptic function, PD, and AD-related traits, and was based on our own recent findings, as well as on the literature (Table 1).
2.3 Genotyping, quality control, and selection of SNPs from the data
Blood assays in LLFS were collected via in-person visits and centrally processed at a Laboratory Core (University of Minnesota). Protocols were standardized, monitored, and coordinated through the LLFS Data Management and Coordinating Center. Genotyping in LLFS was performed by the Center for Inherited Disease Research (CIDR) using Illumina Human Omni 2.5 v1 BeadChip array [details on genotyping are provided in Lee et al. (2014)]. Genotyping in HRS was performed by the CIDR using a similar array. Details on genotyping in HRS are provided by Sonnega et al. (2014).
The genetic QC procedure was performed prior to selection of SNPs for analyses based on recommended procedures (Anderson et al., 2010; Marees et al., 2018). For LLFS, the original data contained 5,036 individuals, and among those, 4,693 (2,581 females, 2,112 males) were genotyped with 2,225,478 SNPs available. For HRS, the starting sample was 10,958 (6,285 females, 4,673 males) White participants (as reported in the data) with available genetic data on 2,443,179 SNPs. In both studies, the sample QC check removed individuals with a call rate below 95% and/or heterozygosity rate beyond ± 4 standard deviations (SD) from the mean as well as individuals of divergent ancestry, i.e., those for whom the first two principal components (PC) scores (computed as described below) were beyond 8 SD from the respective mean scores for the HapMap Phase III European reference populations. The SNPs QC check removed duplicated SNPs, variants with missing allele code information, indels, SNPs with a call rate below 95%, minor allele frequency below 1% (in the entire sample and in the sex-specific subsamples), and those with a significant deviation from Hardy-Weinberg equilibrium with p-value <10–10. In analyses of SNP interactions, an additional criterion was used to keep only SNP pairs that have the sum of the four genotype combinations containing minor alleles to be at least 20 subjects. For LLFS, the resulting sample after QC contained 4,635 individuals of European ancestry (2,548 females, 2,087 males) and 1,468,536 autosomal SNPs, which passed the QC procedure. For HRS, the respective sample was 10,715 participants of European ancestry (6,152 females, 4,563 males) and 1,267,413 post-QC autosomal SNPs. We used these SNPs to select SNPs in the candidate genes shown in Table 1. All available SNPs were first selected if their positions were within the boundary of the genes of interest using information from dbSNP build 151 GRCh38, available from ftp://ftp.ncbi.nlm.nih.gov/snp. SNPs from this initial list were retained for the analyses if they appeared in the post-QC genetic data of the respective dataset. The resulting numbers of SNPs in each candidate gene in LLFS and HRS are shown in Table 1.
The R-package GENESIS version 2.18 (Conomos et al., 2019) was used for the computation of PCs in LLFS using the PC-AiR method (Conomos et al., 2015) to take into account the relatedness among individuals in the LLFS sample. Linkage disequilibrium (LD) pruning was performed before running PC analysis to select a set of independent SNPs for computations of PCs (defined according to the LD threshold 0.2 and a 500 kb sliding window, as implemented in the SNPRelate package (Zheng et al., 2012)). The KING-robust kinship coefficient estimator (Manichaikul et al., 2010) was used to measure ancestry divergence to identify a mutually unrelated and ancestry representative subset of individuals as needed for the PC-AiR algorithm. In HRS, pruning was performed using PLINK (with window size 100 variants, shifted by 10 variants at the end of each step, and a pairwise squared correlation threshold of 0.8 at each step) and then PCs were computed using EIGENSOFT version 6.1.4 (Patterson et al., 2006; Price et al., 2006).
2.4 Association analyses of individual SNPs and SNP pairs
We first estimated associations of SNPs in the selected twelve candidate genes (Table 1; Section 2.2) with the onset of AD in LLFS data. Next, we validated associations of the top SNPs identified in LLFS using the HRS data. We estimated the associations of both individual SNPs and epistatic interactions between SNPs in these candidate genes.
LLFS is a family-based study, and as such, it has related individuals. In addition, our time-to-event outcome of interest is the age of AD onset. To account for relatedness between individuals in analyses of time-to-event outcomes, we used a recently developed R-package coxmeg (He and Kulminski, 2020) that implements a Cox mixed-effects model for time-to-event outcomes in family-based studies that is substantially faster than existing implementations (e.g., R-package coxme). We analyzed individual SNPs in the candidate genes (Table 1), using a block-diagonal family cluster matrix with blocks defined by LLFS family ID. Similarly, we used coxmeg to analyze interactions of pairs of SNPs (only pairs of SNPs taken from different candidate genes were analyzed, e.g., one SNP from ALCAM and one from CADM1 but not pairs of SNPs from ALCAM alone). For the HRS, we used the standard Cox regression model (R-package survival) as the HRS does not contain related individuals.
In addition to genetic variants, the models used in the association analyses included several covariates to adjust for potential confounders. For the LLFS, the following covariates were used: age at first in-person visit (when blood assays were collected), sex (1 - male, 0 - female), country (1 - Denmark, 0 - United States), education (1 - below high school, 0 - otherwise), smoking (smoked >100 cigarettes in a lifetime: yes [1]/no [0]; asked at the baseline visit), and medication use (anti-diabetic, lipid-lowering, anti-hypertensive, angina meds; 1 - used at visit 1 or visit 2, 0 - did not use). For the HRS, we used age at genotyping (which is different from age at first exam in HRS), sex (1 - male, 0 - female), smoking status (1 – ever smoked, 0 – never smoked), and education (1 – below high school, 0 – high school graduate and above). In addition, the top 5 PCs with the largest eigenvalues were included in the regression models to account for the population structure (the number of PCs was determined using scree plots).
Post-analysis QC was performed using easyQC (Winkler et al., 2014) to remove erroneous output, such as missing/excessive estimates or standard errors. After the post-analysis QC, we performed an LD-based clumping of results using the PLINK (Price et al., 2006) clumping procedure with an r-squared threshold of 0.1. Because the HRS does not have all covariates that were included in the LLFS analyses, we repeated the LLFS analyses using the shorter list of covariates (age, sex, smoking, education, PCs) from the HRS analyses. The LLFS results were similar to those using the original list of covariates and thus are not included in the paper.
3 Results
Characteristics of the LLFS and HRS samples are shown in Table 2. We analyzed 4,273 LLFS participants (2,372 females, 1,901 males). Among these, 176 participants (115 for females, 61 for males) reported a new diagnosis of AD at follow-up. For the HRS, the analytic sample included 10,624 participants (6,104 females, 4,520 males) with 336 new cases of AD during the follow-up period (199 among females, 137 among males). See Notes under Table 2 and Materials and Methods that describe how the analytic sample was constructed from the original sample.
TABLE 2. Characteristics of the Long Life Family Study and Health and Retirement Study participants used in the study of associations of SNPs in candidate genes with the onset of Alzheimer’s disease.
No SNPs in the candidate genes from Table 1 individually showed a significant association with AD in the LLFS after correction for multiple comparisons (Supplementary Table S1). However, we found statistically significant interactions between SNPs in the candidate genes. Table 3 shows the top results for associations of these interactions with onset of AD in the LLFS (in total, 29,635 SNP pairs remained after post-analysis QC and clumping). Twenty one SNP pairs had a false discovery rate (FDR) less than 0.05 (Table 3; Supplementary Table S2). The CNTN6, UNC5C, PLXNA4, and EPHB2 genes were broadly involved in the SNPxSNP interactions and appeared in most of the top significant interactions (see Supplementary Figure S1).
Table 4 shows the results of genetic interactions in the HRS data that support (with p-values <0.05) the top two SNP pairs in the CNTN6 – UNC5C and PLXNA4 – EPHB2 pairs found in the LLFS (Table 3). These SNPs are in LD with the top SNPs found in the LLFS (see Notes under Table 4). Note that the interactions between the CNTN6 – UNC5C and PLXNA4 – EPHB2 genes in the HRS also include other SNPs that are not in LD with the SNPs in the genes found in the LLFS (Supplementary Figure S1). This indicates a broad involvement of interactions between these genes in AD, which is not limited by the few specific SNPs. Indeed, different SNPs in the same gene may sometimes similarly influence its expression or protein yield and thus they can be potentially interchangeable in terms of their biological effects (Yashin et al., 2012). Such interactions should also be explored in future studies at the gene level aggregating many SNPs to better understand how multiple functionally similar SNPs may contribute to gene-by-gene interactions that influence AD onset.
TABLE 4. HRS results supporting the top significant SNP pairs in CNTN6 and UNC5C, and PLXNA4 and EPHB2 genes found in LLFS.
4 Discussion
We estimated associations of individual SNPs in candidate genes involved in PD and axon guidance, as well as SNPxSNP interactions, with AD onset among LLFS and HRS participants. Our analysis revealed significant interactions between the SNPs in the UNC5C, CNTN6, PLXNA4, and EPHB2 genes that influenced AD onset in both datasets. For the LLFS, the top significant interaction was observed between rs9307159 (UNC5C) and rs71309807 (CNTN6) (HR = 3.62; Bonferroni p-value: 0.026). The same gene pair was also among the top significant interactions associated with AD in the HRS. Results for individual SNPs were not significant. Our findings, thus, strongly support a major role of genetic interactions in the heterogeneity of AD, including the interactions among genes involved in aging-related increase in physiological dysregulation and changes in axon guidance/synaptic functioning.
Our results also show that selecting candidate genes based on their involvement in the same biological process relevant to AD is a viable strategy for identifying significant genetic interactions associated with AD. Indeed, molecular products of the UNC5C, PLXNA4, and EPHB2 genes are jointly involved in axon guidance (in development) and synaptic plasticity (in adults), essential for the maintenance of complex neural networks. These genes and their products have also been linked to AD and related traits in the literature (Cissé et al., 2011; Barthet et al., 2013; Jun et al., 2014 Wetzel-Smith et al., 2014; Giri et al., 2016; Hashimoto et al., 2016; Kang et al., 2016; Sun et al., 2016; van Dijken et al., 2017; Han et al., 2018; Li et al., 2018; Liu et al., 2018; Lee et al., 2019; Yang et al., 2019; Chen et al., 2021). Results of our analyses suggest that such genes are plausible candidates for exploring the role of genetic interactions in AD heterogeneity. One hypothetical mechanism connecting AD with the axon guidance genes could be reactivation in the aging body of the developmental synaptic pruning program involving such genes, which may lead to synapse loss and subsequent neurodegeneration (Vanderhaeghen and Cheng, 2010). As for CNTN6, it is a member of the contactins family of neuronal membrane proteins that function as cell adhesion molecules playing roles in the formation of axon connections, neurite outgrowth, synaptogenesis, cell survival, and neural circuit formation and regeneration (Oguro-Ando et al., 2017; Chatterjee et al., 2019). Various studies have implicated CNTN6, and more broadly contactins as a group, in metabolic, neurological, and mental health conditions, brain disorders, intellectual disability, depression, and AD-related traits (Jia et al., 2017; Oguro-Ando et al., 2017; Chatterjee et al., 2019; Bamford et al., 2020; Morris et al., 2020). All genes selected for the interaction analysis in this study are involved in the same biological processes relevant to AD, likely contributing to their propensity for interactions affecting AD risk.
Earlier we emphasized that the increase in physiological dysregulation is a fundamental feature of human aging contributing to declines in both resistance to disease occurrence (robustness), and ability to recover and survive after disease occurs (resilience) (Arbeev et al., 2019; Ukraintseva et al., 2021). If so, then the PD-related genes, and their interactions, may potentially influence not only AD onset, but also survival after it, as well as chances of extreme longevity. Therefore, a logical next step would be to investigate associations of the PD-related genes with survival outcomes at different age intervals. Given that the age-related increase in PD may affect the functioning of multiple organs, the PD-related genetic variants might also be associated with comorbidities, offering another avenue in this line of research.
We acknowledge that this study has some limitations. First, our outcome variable is an approximated age of onset derived from questionnaires. As these data are recorded during the interviews, they may not accurately reflect true ages of onset. This is a common limitation of the studies using such outcomes. Second, our studied samples were White, so that more analyses (and data) are necessary to investigate if the results hold for other populations. Third, our samples have relatively small numbers of AD cases. Therefore, additional validation studies are needed to determine if the interactions among UNC5C, CNTN6, PLXNA4, and EPHB2 genes can influence AD onset in other, larger-scale, datasets with time-to-event information on ages of AD onset.
In conclusion, interactions between SNPs (but not individual SNPs) in genes involved in increase in physiological dysregulation with age, a crucial feature of biological aging, and in axon guidance/synaptic functioning, were significantly associated with the age of AD onset in this study. Our findings support the idea that genetic interactions contribute to the heterogeneity of late-onset Alzheimer’s disease. Inherent differences in the rate of PD change with age appear to interact with one’s genetic propensity for synaptic dysfunction in determining AD risk. The results of our study also indicate that system-level changes in the aging body, such as the whole-body increase in physiological dysregulation, might contribute to the brain changes required for AD development.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: the database of Genotypes and Phenotypes (dbGaP) provides access to phenotypic and genetic LLFS data (dbGaP Study Accession: phs000397.v3.p3). The HRS data were provided by dbGaP (dbGaP Study Accession: phs000428.v2.p2) and the University of Michigan (https://hrs.isr.umich.edu/).
Ethics statement
The studies involving humans were approved by the Duke University Health System IRB. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in these studies.
Author contributions
KA conceived and designed the study, supervised statistical analyses and data preparations, and wrote the manuscript. SU contributed to development of research hypotheses, selection of candidate genes, and interpretation of study results, and wrote the manuscript. OB prepared data, performed statistical analyses, and contributed to writing Materials and Methods section. HD and DW contributed to data preparation. IA, AK, ES, KC, MF, JO’C, and AY contributed to writing the manuscript. DP and HW contributed to discussion of the manuscript content. All authors contributed to the article and approved the submitted version.
Funding
Research reported in this publication was supported by the National Institute on Aging of the National Institutes of Health (NIA/NIH) under Award Number. U19AG063893. The work of KA, SU, OB, HD, DW, and HW was partly supported by the NIA/NIH grant R01AG062623. The work of HW was further supported by the Duke/UNC Alzheimer’s Disease Research Center under Award Number. P30AG072958. The work of ES and AY was partly supported by the NIA/NIH grants RF1AG046860 and R01AG063971.
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.
Author disclaimer
This content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1236509/full#supplementary-material
References
An, Y., Varma, V. R., Varma, S., Casanova, R., Dammer, E., Pletnikova, O., et al. (2018). Evidence for brain glucose dysregulation in Alzheimer's disease. Alzheimers Dement. 14, 318–329. doi:10.1016/j.jalz.2017.09.011
Anderson, C. A., Pettersson, F. H., Clarke, G. M., Cardon, L. R., Morris, A. P., and Zondervan, K. T. (2010). Data quality control in genetic case-control association studies. Nat. Protoc. 5, 1564–1573. doi:10.1038/nprot.2010.116
Arbeev, K. G., Bagley, O., Ukraintseva, S. V., Wu, D., Duan, H., Kulminski, A. M., et al. (2020). Genetics of physiological dysregulation: findings from the long life family study using joint models. Aging (Albany NY) 12, 5920–5947. doi:10.18632/aging.102987
Arbeev, K. G., Cohen, A. A., Arbeeva, L. S., Milot, E., Stallard, E., Kulminski, A. M., et al. (2016). Optimal versus realized trajectories of physiological dysregulation in aging and their relation to sex-specific mortality risk. Front. Public Health 4, 3. doi:10.3389/fpubh.2016.00003
Arbeev, K. G., Ukraintseva, S. V., Bagley, O., Zhbannikov, I. Y., Cohen, A. A., Kulminski, A. M., et al. (2019). Physiological dysregulation" as a promising measure of robustness and resilience in studies of aging and a new indicator of preclinical disease. J. Gerontol. A Biol. Sci. Med. Sci. 74, 462–468. doi:10.1093/gerona/gly136
Bamford, R. A., Widagdo, J., Takamura, N., Eve, M., Anggono, V., and Oguro-Ando, A. (2020). The interaction between contactin and amyloid precursor protein and its role in Alzheimer's disease. Neuroscience 424, 184–202. doi:10.1016/j.neuroscience.2019.10.006
Barthet, G., Dunys, J., Shao, Z., Xuan, Z., Ren, Y., Xu, J., et al. (2013). Presenilin mediates neuroprotective functions of ephrinB and brain-derived neurotrophic factor and regulates ligand-induced internalization and metabolism of EphB2 and TrkB receptors. Neurobiol. Aging 34, 499–510. doi:10.1016/j.neurobiolaging.2012.02.024
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Braschi, B., Denny, P., Gray, K., Jones, T., Seal, R., Tweedie, S., et al. (2019). Genenames.org: the HGNC and VGNC resources in 2019. Nucleic Acids Res. 47, D786–d792. doi:10.1093/nar/gky930
Chanti-Ketterl, M., Stebbins, R. C., Obhi, H. K., Belsky, D. W., Plassman, B. L., and Yang, Y. C. (2021). Sex differences in the association between metabolic dysregulation and cognitive aging: the health and retirement study. J. Gerontol. A Biol. Sci. Med. Sci. 77, 1827–1835. doi:10.1093/gerona/glab285
Chatterjee, M., Schild, D., and Teunissen, C. E. (2019). Contactins in the central nervous system: role in health and disease. Neural Regen. Res. 14, 206–216. doi:10.4103/1673-5374.244776
Chen, G., Kang, S. S., Wang, Z., Ahn, E. H., Xia, Y., Liu, X., et al. (2021). Netrin-1 receptor UNC5C cleavage by active δ-secretase enhances neurodegeneration, promoting Alzheimer's disease pathologies. Sci. Adv. 7, eabe4499. doi:10.1126/sciadv.abe4499
Cissé, M., Halabisky, B., Harris, J., Devidze, N., Dubal, D. B., Sun, B., et al. (2011). Reversing EphB2 depletion rescues cognitive functions in Alzheimer model. Nature 469, 47–52. doi:10.1038/nature09635
Conomos, M. P., Miller, M. B., and Thornton, T. A. (2015). Robust inference of population structure for ancestry prediction and correction of stratification in the presence of relatedness. Genet. Epidemiol. 39, 276–293. doi:10.1002/gepi.21896
Cohen, A. A., Milot, E., Li, Q., Legault, V., Fried, L. P., and Ferrucci, L. (2014). Cross-population validation of statistical distance as a measure of physiological dysregulation during aging. Exp. Gerontol. 57, 203–210. doi:10.1016/j.exger.2014.04.016
Cohen, A. A., Milot, E., Yong, J., Seplaki, C. L., Fülöp, T., Bandeen-Roche, K., et al. (2013). A novel statistical approach shows evidence for multi-system physiological dysregulation during aging. Mech. Ageing Dev. 134, 110–117. doi:10.1016/j.mad.2013.01.004
Conomos, M. P., Gogarten, S. M., Brown, L., Chen, H., Rice, K., Sofer, T., et al. (2019). Genesis: GENetic EStimation and inference in structured samples (GENESIS): Statistical methods for analyzing genetic data from samples with population structure and/or relatedness. R package version 2.14.1, Available at: https://github.com/UW-GAC/GENESIS.
Cuestas Torres, D. M., and Cardenas, F. P. (2020). Synaptic plasticity in Alzheimer's disease and healthy aging. Rev. Neurosci. 31, 245–268. doi:10.1515/revneuro-2019-0058
Dansereau, G., Wey, T. W., Legault, V., Brunet, M. A., Kemnitz, J. W., Ferrucci, L., et al. (2019). Conservation of physiological dysregulation signatures of aging across primates. Aging Cell 18, e12925. doi:10.1111/acel.12925
De La Monte, S. M., Tong, M., Daiello, L. A., and Ott, B. R. (2019). Early-stage Alzheimer's disease is associated with simultaneous systemic and central nervous system dysregulation of insulin-linked metabolic pathways. J. Alzheimers Dis. 68, 657–668. doi:10.3233/JAD-180906
Elo, I. T., Mykyta, L., Sebastiani, P., Christensen, K., Glynn, N. W., and Perls, T. (2013). Age validation in the long life family study through a linkage to early-life census records. Journals Gerontology Ser. B Psychol. Sci. Soc. Sci. 68, 580–585. doi:10.1093/geronb/gbt033
Evans, E. E., Mishra, V., Mallow, C., Gersz, E. M., Balch, L., Howell, A., et al. (2022). Semaphorin 4D is upregulated in neurons of diseased brains and triggers astrocyte reactivity. J. Neuroinflammation 19, 200. doi:10.1186/s12974-022-02509-8
Fagiani, F., Lanni, C., Racchi, M., and Govoni, S. (2021). Dys)regulation of synaptic activity and neurotransmitter release by β-amyloid: A look beyond Alzheimer's disease pathogenesis. Front. Mol. Neurosci. 14, 635880. doi:10.3389/fnmol.2021.635880
Giri, M., Zhang, M., and Lü, Y. (2016). Genes associated with Alzheimer's disease: an overview and current status. Clin. Interv. Aging 11, 665–681. doi:10.2147/CIA.S105769
Han, Q., Sun, Y. A., Zong, Y., Chen, C., Wang, H. F., Tan, L., et al. (2018). Common variants in PLXNA4 and correlation to CSF-related Phenotypes in Alzheimer's disease. Front. Neurosci. 12, 946. doi:10.3389/fnins.2018.00946
Hashimoto, Y., Toyama, Y., Kusakari, S., Nawa, M., and Matsuoka, M. (2016). An alzheimer disease-linked rare mutation potentiates netrin receptor uncoordinated-5C-induced signaling that merges with amyloid β precursor protein signaling. J. Biol. Chem. 291, 12282–12293. doi:10.1074/jbc.M115.698092
He, L., and Kulminski, A. M. (2020). Fast algorithms for conducting large-scale GWAS of age-at-onset traits using Cox mixed-effects models. Genetics 215, 41–58. doi:10.1534/genetics.119.302940
Jeong, W., Lee, H., Cho, S., and Seo, J. (2019). ApoE4-Induced cholesterol dysregulation and its brain cell type-specific implications in the pathogenesis of Alzheimer's disease. Mol. Cells 42, 739–746. doi:10.14348/molcells.2019.0200
Jia, P., Zhao, Z., Hulgan, T., Bush, W. S., Samuels, D. C., Bloss, C. S., et al. (2017). Genome-wide association study of HIV-associated neurocognitive disorder (hand): A charter group study. Am. J. Med. Genet. B Neuropsychiatr. Genet. 174, 413–426. doi:10.1002/ajmg.b.32530
Jiang, J., Wang, Z. H., Qu, M., Gao, D., Liu, X. P., Zhu, L. Q., et al. (2015). Stimulation of EphB2 attenuates tau phosphorylation through PI3K/Akt-mediated inactivation of glycogen synthase kinase-3β. Sci. Rep. 5, 11765. doi:10.1038/srep11765
Joshi, M., Joshi, S., Khambete, M., and Degani, M. (2022). Role of calcium dysregulation in Alzheimer's disease and its therapeutic implications. Chem. Biol. Drug Des. 101, 453–468. doi:10.1111/cbdd.14175
Jouini, N., Saied, Z., Ben Sassi, S., Nebli, F., Messaoud, T., Hentati, F., et al. (2021). Impacts of iron metabolism dysregulation on Alzheimer's disease. J. Alzheimers Dis. 80, 1439–1450. doi:10.3233/JAD-201250
Jun, G., Asai, H., Zeldich, E., Drapeau, E., Chen, C., Chung, J., et al. (2014). PLXNA4 is associated with Alzheimer disease and modulates tau phosphorylation. Ann. Neurol. 76, 379–392. doi:10.1002/ana.24219
Kang, S. S., Kurti, A., Wojtas, A., Baker, K. E., Liu, C.-C., Kanekiyo, T., et al. (2016). Identification of plexin A4 as a novel clusterin receptor links two Alzheimer’s disease risk genes. Hum. Mol. Genet. 25, 3467–3475. doi:10.1093/hmg/ddw188
Lee, J., Cheng, R., Honig, L., Feitosa, M., Kammerer, C., Kang, M., et al. (2014). Genome wide association and linkage analyses identified three loci—4q25, 17q23.2, and 10q11.21—associated with variation in leukocyte telomere length: the long life family study. Front. Genet. 4, 310. doi:10.3389/fgene.2013.00310
Lee, W. S., Lee, W. H., Bae, Y. C., and Suk, K. (2019). Axon guidance molecules guiding neuroinflammation. Exp. Neurobiol. 28, 311–319. doi:10.5607/en.2019.28.3.311
Li, Q., Wang, B. L., Sun, F. R., Li, J. Q., Cao, X. P., and Tan, L. (2018). The role of UNC5C in Alzheimer's disease. Ann. Transl. Med. 6, 178. doi:10.21037/atm.2018.04.43
Liu, Z., Thakar, A., Santoro, S. W., and Pratt, K. G. (2018). Presenilin regulates retinotectal synapse formation through EphB2 receptor processing. Dev. Neurobiol. 78, 1171–1190. doi:10.1002/dneu.22638
Lu, Y., Pike, J. R., Kucharska-Newton, A., Palta, P., Whitsel, E. A., Bey, G., et al. (2022). Aging-related multisystem dysregulation over the adult lifespan and physical function in later life: the atherosclerosis risk in communities (ARIC) study. J. Gerontol. A Biol. Sci. Med. Sci. 78, 1497–1503. doi:10.1093/gerona/glac236
Lutshumba, J., Nikolajczyk, B. S., and Bachstetter, A. D. (2021). Dysregulation of systemic immunity in aging and dementia. Front. Cell Neurosci. 15, 652111. doi:10.3389/fncel.2021.652111
Ma, T., Chang, R. C.-C., and Macauley, S. L. (2023). Editorial: metabolic signaling dysregulation and cognitive impairments in aging and Alzheimer's disease, volume II. Front. Aging Neurosci. 15, 1150101. doi:10.3389/fnagi.2023.1150101
Ma, X., Zhang, Y., Gou, D., Ma, J., Du, J., Wang, C., et al. (2022). Metabolic reprogramming of microglia enhances proinflammatory cytokine release through EphA2/p38 MAPK pathway in Alzheimer's disease. J. Alzheimers Dis. 88, 771–785. doi:10.3233/JAD-220227
Manichaikul, A., Mychaleckyj, J. C., Daly, K., Sale, M., Rich, S. S., and Chen, W.-M. (2010). Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873. doi:10.1093/bioinformatics/btq559
Marees, A. T., De Kluiver, H., Stringer, S., Vorspan, F., Curis, E., Marie-Claire, C., et al. (2018). A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int. J. Methods Psychiatric Res. 27, e1608. doi:10.1002/mpr.1608
Martínez-Serra, R., Alonso-Nanclares, L., Cho, K., and Giese, K. P. (2022). Emerging insights into synapse dysregulation in Alzheimer’s disease. Brain Commun. 4, fcac083. doi:10.1093/braincomms/fcac083
Mcdermott, J. E., Goldblatt, D., and Paradis, S. (2018). Class 4 Semaphorins and Plexin-B receptors regulate GABAergic and glutamatergic synapse development in the mammalian hippocampus. Mol. Cell Neurosci. 92, 50–66. doi:10.1016/j.mcn.2018.06.008
Morris, J., Leung, S. S. Y., Bailey, M. E. S., Cullen, B., Ferguson, A., Graham, N., et al. (2020). Exploring the role of contactins across psychological, psychiatric and cardiometabolic traits within UK biobank. Genes (Basel) 11, 1326. doi:10.3390/genes11111326
Müller, L., Di Benedetto, S., and Pawelec, G. (2019). The immune system and its dysregulation with aging. Subcell. Biochem. 91, 21–43. doi:10.1007/978-981-13-3681-2_2
Nanclares, C., Baraibar, A. M., Araque, A., and Kofuji, P. (2021). Dysregulation of astrocyte-neuronal communication in Alzheimer's disease. Int. J. Mol. Sci. 22, 7887. doi:10.3390/ijms22157887
Oguro-Ando, A., Zuko, A., Kleijer, K. T. E., and Burbach, J. P. H. (2017). A current view on contactin-4, -5, and -6: implications in neurodevelopmental disorders. Mol. Cell Neurosci. 81, 72–83. doi:10.1016/j.mcn.2016.12.004
Patterson, N., Price, A. L., and Reich, D. (2006). Population structure and eigenanalysis. PLOS Genet. 2, e190. doi:10.1371/journal.pgen.0020190
Piehl, N., Van Olst, L., Ramakrishnan, A., Teregulova, V., Simonton, B., Zhang, Z., et al. (2022). Cerebrospinal fluid immune dysregulation during healthy brain aging and cognitive impairment. Cell 185, 5028–5039.e13. doi:10.1016/j.cell.2022.11.019
Price, A. L., Patterson, N. J., Plenge, R. M., Weinblatt, M. E., Shadick, N. A., and Reich, D. (2006). Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909. doi:10.1038/ng1847
Sebastiani, P., Hadley, E. C., Province, M., Christensen, K., Rossi, W., Perls, T. T., et al. (2009). A family longevity selection score: ranking sibships by their longevity, size, and availability for study. Am. J. Epidemiol. 170, 1555–1562. doi:10.1093/aje/kwp309
Sonnega, A., Faul, J. D., Ofstedal, M. B., Langa, K. M., Phillips, J. W. R., and Weir, D. R. (2014). Cohort profile: the health and retirement study (HRS). Int. J. Epidemiol. 43, 576–585. doi:10.1093/ije/dyu067
Sun, J. H., Wang, H. F., Zhu, X. C., Yu, W. J., Tan, C. C., Jiang, T., et al. (2016). The impact of UNC5C genetic variations on neuroimaging in Alzheimer's disease. Mol. Neurobiol. 53, 6759–6767. doi:10.1007/s12035-015-9589-0
Tan, Z. S., Beiser, A. S., Fox, C. S., Au, R., Himali, J. J., Debette, S., et al. (2011). Association of metabolic dysregulation with volumetric brain magnetic resonance imaging and cognitive markers of subclinical brain aging in middle-aged adults: the framingham offspring study. Diabetes Care 34, 1766–1770. doi:10.2337/dc11-0308
Tang, X., Jang, S. W., Okada, M., Chan, C. B., Feng, Y., Liu, Y., et al. (2008). Netrin-1 mediates neuronal survival through PIKE-L interaction with the dependence receptor UNC5B. Nat. Cell Biol. 10, 698–706. doi:10.1038/ncb1732
Ukraintseva, S., Arbeev, K., Duan, M., Akushevich, I., Kulminski, A., Stallard, E., et al. (2021). Decline in biological resilience as key manifestation of aging: potential mechanisms and role in health and longevity. Mech. Ageing Dev. 194, 111418. doi:10.1016/j.mad.2020.111418
Van Dijken, I., Van Der Vlag, M., Flores Hernández, R., and Ross, A. (2017). Perspectives on treatment of Alzheimer's disease: A closer look into EphB2 depletion. J. Neurosci. 37, 11296–11297. doi:10.1523/JNEUROSCI.0214-17.2017
Vanderhaeghen, P., and Cheng, H. J. (2010). Guidance molecules in axon pruning and cell death. Cold Spring Harb. Perspect. Biol. 2, a001859. doi:10.1101/cshperspect.a001859
Von Bernhardi, R., Eugenín-Von Bernhardi, L., and Eugenín, J. (2015). Microglial cell dysregulation in brain aging and neurodegeneration. Front. Aging Neurosci. 7, 124. doi:10.3389/fnagi.2015.00124
Wetzel-Smith, M. K., Hunkapiller, J., Bhangale, T. R., Srinivasan, K., Maloney, J. A., Atwal, J. K., et al. (2014). A rare mutation in UNC5C predisposes to late-onset Alzheimer's disease and increases neuronal cell death. Nat. Med. 20, 1452–1457. doi:10.1038/nm.3736
Whelan, C. D., Mattsson, N., Nagle, M. W., Vijayaraghavan, S., Hyde, C., Janelidze, S., et al. (2019). Multiplex proteomics identifies novel CSF and plasma biomarkers of early Alzheimer’s disease. Acta Neuropathol. Commun. 7, 169. doi:10.1186/s40478-019-0795-2
Whitson, H. E., Duan-Porter, W., Schmader, K. E., Morey, M. C., Cohen, H. J., and Colón-Emeric, C. S. (2016). Physical resilience in older adults: systematic Review and development of an emerging construct. J. Gerontol. A Biol. Sci. Med. Sci. 71, 489–495. doi:10.1093/gerona/glv202
Winkler, T. W., Day, F. R., Croteau-Chonka, D. C., Wood, A. R., Locke, A. E., Mägi, R., et al. (2014). Quality control and conduct of genome-wide association meta-analyses. Nat. Protoc. 9, 1192–1212. doi:10.1038/nprot.2014.071
Wojczynski, M. K., Jiuan Lin, S., Sebastiani, P., Perls, T. T., Lee, J., Kulminski, A., et al. (2022). NIA long life family study: objectives, design, and heritability of cross-sectional and longitudinal Phenotypes. Journals Gerontology. 77, 717–727. doi:10.1093/gerona/glab333
Yan, X., Hu, Y., Wang, B., Wang, S., and Zhang, X. (2020). Metabolic dysregulation contributes to the progression of Alzheimer's disease. Front. Neurosci. 14, 530219. doi:10.3389/fnins.2020.530219
Yang, H. S., Chhatwal, J. P., Xu, J., White, C. C., Hanseeuw, B., Rabin, J. S., et al. (2019). An UNC5C allele predicts cognitive decline and hippocampal atrophy in clinically normal older adults. J. Alzheimers Dis. 68, 1161–1170. doi:10.3233/JAD-180788
Yashin, A. I., Wu, D., Arbeev, K. G., and Ukraintseva, S. V. (2012). Polygenic effects of common single-nucleotide polymorphisms on life span: when association meets causality. Rejuvenation Res. 15, 381–394. doi:10.1089/rej.2011.1257
Zhang, L., Qi, Z., Li, J., Li, M., Du, X., Wang, S., et al. (2021). Roles and mechanisms of axon-guidance molecules in Alzheimer's disease. Mol. Neurobiol. 58, 3290–3307. doi:10.1007/s12035-021-02311-2
Zheng, X., Levine, D., Shen, J., Gogarten, S. M., Laurie, C., and Weir, B. S. (2012). A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 28, 3326–3328. doi:10.1093/bioinformatics/bts606
Keywords: aging, Alzheimer’s disease, physiological dysregulation, resilience, axon guidance, synaptic function, genetic interactions, candidate genes
Citation: Arbeev KG, Ukraintseva S, Bagley O, Duan H, Wu D, Akushevich I, Stallard E, Kulminski A, Christensen K, Feitosa MF, O’Connell JR, Parker D, Whitson H and Yashin AI (2023) Interactions between genes involved in physiological dysregulation and axon guidance: role in Alzheimer’s disease. Front. Genet. 14:1236509. doi: 10.3389/fgene.2023.1236509
Received: 07 June 2023; Accepted: 17 August 2023;
Published: 31 August 2023.
Edited by:
Xiaomin Zhang, University of Arkansas for Medical Sciences, United StatesReviewed by:
Dandan Chu, Nantong University, ChinaWinnie Liang, Translational Genomics Research Institute (TGen), United States
Copyright © 2023 Arbeev, Ukraintseva, Bagley, Duan, Wu, Akushevich, Stallard, Kulminski, Christensen, Feitosa, O’Connell, Parker, Whitson and Yashin. 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: Konstantin G. Arbeev, ka29@duke.edu; Svetlana Ukraintseva, svo@duke.edu
†These authors have contributed equally to this work and share first authorship