Abstract
Background:
The E4 allele of the APOE gene is known to be associated with cognitive impairment. However, a limited number of studies have examined the association between the E2 allele and longitudinal changes of cognitive function.
Objective:
To determine whether rates of cognitive change differ in carriers of the APOE E2 allele compared to other genotypes.
Methods:
We conducted a secondary analysis of data from two ongoing longitudinal cohort studies, the Long Life Family Study (LLFS) and New England Centenarian Study (NECS). We included participants who had APOE genotyping data, data from longitudinal administrations of the Telephone Interview for Cognitive Status (TICS), and age, sex, and education available. We assessed whether cognitive change as measured by rate of decline in TICS score differed among people with different APOE genotypes. We used a hierarchical mixed effect model with APOE genotypes, their interactions with age, and potential confounders.
Results:
After adjusting for sex and education, in carriers of the common E3/E3 genotype, TICS score decreased by 0.15 points per year of age. In those with the E2/E2 genotype, TICS score decreased by 0.05 points per year of age, a significantly slower rate of decline (p = 0.017). We observed no protective effect of the E2/E3 genotype on cognitive decline.
Conclusion:
These results suggest a protective effect of the E2/E2 genotype on a measure of global cognitive function.
INTRODUCTION
Cognitive decline reduces the independence and quality of life of older individuals. As our population ages, it hinders productivity and adds significant costs to the health care system. Identification of genes associated with cognitive decline may improve screening and early detection, thereby providing the greatest benefit from treatment. Further, identification of protective loci may clarify the biological pathways causing cognitive decline, leading to more effective, personalized treatments [1].
An important genetic marker for cognitive decline is apolipoprotein E (APOE) [2]. The APOE protein is involved in lipid transport and may affect neurodegeneration through mediating amyloid-β clearance [3]. APOE gene has three alleles, of which the E4 allele is well established as a major risk factor for Alzheimer’s disease [3]. Less is known about E2, which is the rarest allele in whites, but there is evidence that it may have a protective effect against neurodegeneration [4]. APOE E2 is associated with extreme human longevity [5], as well as with decreased risk of Alzheimer’s disease [6, 7]. In the area of cognition, there is evidence that the E2 allele is protective against declines in episodic memory function [8, 9]. However, much remains unknown about the protective effect of APOE E2. For instance, it is unclear how age modifies the protective effect, as studies on middle-aged adults have frequently not found evidence of improved cognition [10] or differences in brain structure [4] among E2 carriers.
Another important knowledge gap is whether E2 homozygous and E2 heterozygous genotypes are associated with differential effects. There is evidence that this is the case for the E4 allele, with the homozygous genotype being more deleterious, and the heterozygous genotype possibly exerting an age-varying effect [11, 12]. Presumably because of rarity, longitudinal studies have included relatively few E2 carriers [8] and thus often grouped E2 homozygous and heterozygous genotypes for analysis. To address this knowledge gap, we analyzed data from two large longitudinal cohort studies of healthy aging populations with repeated cognitive testing and relatively large numbers of E2 homozygotes. We sought to separately assess whether the E2/E2, E2/E3, and E2/E4 genotypes are associated with differences in rate of cognitive decline.
MATERIALS AND METHODS
Study design and population
We conducted a longitudinal analysis of participants enrolled in two observational cohort studies of familial longevity. The Long Life Family Study (LLFS) is a multisite study of healthy aging and longevity with approximately 4,900 enrolled participants who belong to families with long-lived individuals [13, 14]. The New England Centenarian Study (NECS) is a study of longevity among centenarians and their family members, with approximately 4,000 enrolled participants [15]. For both studies, participants are followed up annually to assess vital and health status, including assessing cognitive function longitudinally. These studies seek to identify factors that contribute to longevity and healthy aging. The analysis presented here was not prespecified.
Genotype data
APOE alleles were defined from the SNPs rs7412 and rs429358 as E2: rs7412 = T; rs429358 = T, E3: rs7412 = C; rs429358 = T, E4: rs7412 = C; rs429358 = C that were directly genotyped in the LLFS participants and in 80% of the NECS participants and imputed using the Michigan Imputation Server to the HRC panel (Version r1.1 2016) of ∼65,000 haplotypes in 20% of NECS participants [5, 16].
Outcome of interest
Our outcome of interest was participants’ score on a modified version of the Telephone Interview for Cognitive Status (TICS) test [17, 18], and we estimated the average annual change in TICS score as a measure of cognitive decline. For NECS, the TICS test is administered every two years. For LLFS, the TICS test is administered every year for participants aged 70 and above, and every three years for participants aged less than 70. For each TICS administration, we calculated a score based on four tasks that cover the domains of episodic memory, working memory, and executive function as in [19]. Tasks included counting backwards, subtracting sevens, and immediate and delayed word list recall resulting in a score that ranged from 0 to 27. We chose these domains because previous studies have suggested that the E2 allele is protective against decline in episodic memory, while evidence for other cognitive domains is weaker [9]. We included in the analysis participants with at least one completed cognitive function assessment, APOE genotyping data, and age, sex, and education available. The same version of the TICS was administered in LLFS and NECS using the same protocol.
Statistical analysis
The data from the two studies were aggregated for the main analyses. For descriptive analysis, we calculated the proportion of participants with demographic and medical history characteristics of interest for each genotype. We described age and first TICS score in terms of mean and standard deviation, and number of TICS assessments as median and interquartile range. We tested whether these characteristics differed significantly across genotypes using a chi-squared test for proportions, ANOVA for age at baseline and first TICS score, and Kruskal-Wallis test for number of TICS assessments.
Our initial model to determine whether APOE genotypes modified the association between age and TICS score used a linear mixed effect model to predict TICS score based on age, APOE genotype, and interaction terms between age and APOE genotypes. Statistical significance in the interaction term was interpreted as evidence that genotypes modify the rate of change in TICS score with age. Age was centered at the baseline mean. We chose to use age as a timescale rather than years since enrollment (while adjusting for baseline age) because the range of ages represented by the study participants was far greater than the average length of time over which individuals have TICS measurements. Thus, we assume cross-sectional and longitudinal effects to be similar. We used a random intercept to account for repeated observations of the same individuals and assessed the need for a random slope for age using a likelihood ratio test. The E3/E3 genotype was the reference group for all analyses. The formula for the initial model is as follows:
Where
TICS ij : The j th score for the i th participant
cAge ij : The age, centered at the baseline mean, of the i th participant at the time of his/her j th score
E2E2i, E2E3i, E2E4i, E3E4i, E4E4i: Indicator variables for the APOE genotype of the i th participant
∈ ij ∼ N (0, σ2): Random within-participant error term
We then assessed potential confounders by separately adding to this initial model sex, education, self-reported history of medical conditions known to be associated with cognitive decline (diabetes, heart failure, heart attack, depression or anxiety, stroke, hypertension), history of smoking, and study cohort (LLFS or NECS). We identified covariates that were both significantly associated with TICS score when assessed individually in the context of the initial model, and that, when included changed the estimate for the E2/E2*age interaction term by more than 10%. Our final adjusted model was created by adding to the initial model any covariates that met these two criteria for potential confounders, after first assessing them for statistical interaction with APOE genotype.
Finally, we conducted subgroup analyses because our study population was drawn from two cohort studies with different inclusion criteria, and associations could differ between the two cohorts. The subgroup analyses applied the final adjusted model to participants from each study separately. For all analyses, p < 0.05 in a two-sided test was considered statistically significant; all analyses were conducted in SAS v9.4 (SAS Institute, Cary, NC).
RESULTS
We analyzed data from 2,933 people from the LLFS and 679 people from the NECS who had APOE genotyping data, TICS data, and age, sex, and education available. Of the combined analytic population, 2,019 (55.9%) were female and 3,575 (99.1%) were non-Hispanic white (Table 1; see Supplementary Tables 1 and 2 for participant characteristics from the two study cohorts). The mean age at first TICS administration was 72.1 (standard deviation 13.4) years. Subjects had a median of 3 TICS assessments (interquartile range [IQR] 2, 4) and a median follow-up time of 5.0 years (IQR 2.9, 7.0). The most common genotype in this population was E3/E3, comprising 66.3% of subjects (n = 2,394). The two rarest genotypes were E2/E2 (n = 23, 0.6%) and E4/E4 (n = 30, 0.8%).
Demographic characteristics of study population, by genotype (N = 3,612)
Percentages shown in column headers are the percent of the subjects with the indicated genotype. aParticipants missing data: Race (n = 5), Diabetes (n = 10), Heart Failure (n = 13), Heart Attack (n = 10), Depression or Anxiety (n = 12), Stroke (n = 13), Hypertension (n = 11), Ever Smoked (n = 15). Percentages are taken out of those with data.
In all genotype groups, TICS score decreased with age (Fig. 1). In the initial model using age, genotype, and an age–genotype interaction term to predict TICS score, the inclusion of a random slope for age did not improve the fit (p = 0.061), therefore we included only a random intercept. The E2/E2 genotype was associated with a significantly slower rate of decline with age in TICS score than the E3/E3 genotype (p = 0.049), but none of the other genotypes differed significantly from E3/E3 (Table 2). In this model the E3/E3 genotype was associated with an annual change in TICS score of –0.17 (95% confidence interval [CI] –0.18, –0.16) points per year of age (Table 3). In contrast, the E2/E2 genotype was associated with a change in TICS score of –0.08 points per year of age which was not significantly different from 0 (95% CI –0.17, 0.01). The E4/E4 genotype was associated with a change in TICS score of –0.23 (95% CI –0.32, –0.14) points per year of age, but this was not significantly different from E3/E3 since the interaction between the E4/E4 genotype and age was not significant (p = 0.170).

Observed and predicted TICS scores by age and genotype. Connected points represent observed TICS scores for an individual. The solid line shows the predicted TICS score based on a model using age, genotype, and an age–genotype interaction term as predictors.
Initial model for TICS score based on age and genotype
aAge was centered at the baseline mean of 72.
Estimated change in TICS score per year of age, by genotype, based on initial and final adjusted models
Estimates are the sum of the main effect of age and the corresponding age*APOE genotype interaction effect (See Table 2 and Table 5 for the individual parameter estimates). aAdjusted for sex and education.
In our assessment of potential confounders, sex, education, diabetes, heart failure, heart attack, and stroke were significantly associated with TICS score when separately added to the individual model (Table 4). However, of these, only sex and education also changed the estimate for the E2/E2*age interaction term by more than 10% (Table 5). No significant statistical interaction was detected between genotype and either sex or education. When we assessed interaction terms between each potential covariate and age, none changed the E2/E2*age interaction estimate by more than 10%.
Bivariate associations between potential confounders and TICS score when added to initial model containing age and genotype
Race was not assessed because of lack of non-white patients.
Final adjusted model for TICS score based on age and genotype, adjusted for sex and education
aAge was centered at the baseline mean of 72.
In the final adjusted model, E2/E2 was again the only genotype associated with a difference in rate of TICS score decline with age compared to E3/E3 (p = 0.017, Table 5). In this model, the E3/E3 genotype was associated with an annual change in TICS score of –0.15 (95% CI –0.16, –0.15) points per year of age (Table 3). The E2/E2 genotype was associated with a change in TICS score of –0.05 that was not significantly different from 0 (95% CI –0.14, 0.04). The E4/E4 genotype was associated with a change in TICS score of –0.20 (95% CI –0.29, –0.12) points per year of age, but this was not significantly different from E3/E3 (p = 0.252).
In subgroup analyses, the magnitude of the interaction term between E2/E2 and age was comparable in the LLFS and NECS study cohorts, although the interaction did not achieve statistical significance in the NECS cohort (Supplementary Table 3).
DISCUSSION
We analyzed the association between APOE genotypes and rate of decline of TICS scores with age in a large set of older individuals enrolled in two longitudinal studies of healthy aging and longevity. The analysis showed that TICS score was positively correlated with education and female sex and, in carriers of the common E3/E3 genotype, TICS score was estimated to decrease on average by 0.15 points per year of age, which represents a change of 0.5% of the maximum score. The E2/E2 genotype was associated with a significantly slower rate of TICS decline with age while, in contrast, we observed no protective effect of the E2/E3 genotype that had previously been found to be associated with decreased risk of Alzheimer’s disease [2, 6]. Moreover, our analysis confirmed a significant negative effect of the E3/E4 and E4/E4 genotypes on TICS score, but no significantly increased rate of decline. Taken together, these results are consistent with a protective effect of the E2/E2 genotype on episodic memory, working memory, and executive function domains of cognitive function.
Our findings add to the growing evidence for a protective effect of the E2 allele, while potentially adding information to distinguish homozygous from heterozygous effects. Previous longitudinal studies of cognitive decline have shown that people with either the E2/E2 or E2/E3 genotype experience less decline in tests of episodic memory over time than people with E3 and E4 alleles [9, 20–22]. However, these studies grouped the E2/E2 homozygous and E2/E3 “carrier” genotypes and excluded the E2/E4 genotype from analysis. Our study separately assessed all the possible genotypes and found that only the E2/E2 homozygous genotype was associated with slower decline with age of TICS scores. Despite having over 500 people with the E2/E3 genotype in our analysis, we found no evidence for a protective effect of this genotype, as our model predicted indistinguishable rates of cognitive decline for the E2/E3 and E3/E3 genotypes. While previous studies have found that the E2/E3 is protective against Alzheimer’s disease [6], it is possible that the outcome of TICS score decline measures cognitive processes whose mechanistic relationship with APOE is distinct from that by which APOE mediates Alzheimer’s disease. Our analysis suggests that molecular mechanisms that are associated with homozygosity of the E2 allele may provide better targets for the discovery of therapeutics that can preserve cognitive functions. We recently discovered a serum protein signature that is associated with APOE genotypes and indeed showed that the distribution of some of these proteins such as APOB and BIRC2 is drastically different in carriers of the E2/E2 genotype [23].
It is notable that we did not see any association between the E4 allele and rate of cognitive decline with age. In previous longitudinal studies, the E4/E4 genotype has been found to be associated with faster cognitive decline by various measures [8, 24], including as measured by TICS [25]. While the estimated rate of cognitive decline in our study was faster for the E4/E4 genotype in our study population, it is possible that the small number of people with this genotype (n = 30) may have prevented us from detecting a statistically significant association. Our analysis confirmed a significant effect of sex on TICS, with females on average scoring 1.37 points higher than males. The sex difference did not change with age, or APOE genotype. A strength of our analysis is that we had repeated assessments of TICS over time, thus making the results more robust. Our study has limitations related to the population of the LLFS and NECS studies. First, the populations of participants enrolled in the LLFS and NECS studies differ from the general population in ways that could affect the association between genetics and cognition. In particular, both participant populations were selected because they are likely to have genetic factors promoting longevity, and both cohorts have high levels of education compared to the national age cohort [13]. Given that some of the genetic factors associated with longevity may also protect against cognitive decline, and that education is known to have a protective association with delayed onset of cognitive decline [26, 27], the results should be replicated in less selected populations. Moreover, the numbers of participants with the E2/E2, E2/E4, and E4/E4 genotypes were relatively small, and some of these groups differed in key demographic characteristics from the study population as a whole. We attempted to control for these differences by adjusting for age, sex, and education, and we assessed various comorbidities as well as study cohort as potential confounders. However, we lacked information about other factors that have been shown to be associated with cognitive decline, such as exercise and social engagement [28].
Our study also has limitations related to our ability to assess cognitive decline. While changes in TICS score have been used as a measure of cognitive decline by other studies [28], TICS is known to have reduced sensitivity for detecting mild cognitive impairment as opposed to dementia [29]. Because the length of follow-up time for individual participants was limited, we may have had insufficient outcome data to accurately detect cognitive decline in many cases. Moreover, because we used age rather than time since enrollment as a timescale, we were unable to distinguish individual-level effects (i.e., longitudinal associations with the TICS scores of individual participants) from age cohort effects (i.e., associations with population-level differences in TICS scores among participants of different ages) in our analysis.
In conclusion, our analysis confirms the protective effect of the E2 allele of APOE on cognitive decline and the deleterious effect of the E4 allele of APOE on level of cognitive function. The contribution of our work is in the detection of a protective effect of the homozygous genotype E2/E2 that could be helpful to target more specific molecular biomarkers of cognitive function and to develop therapeutic interventions.
Footnotes
ACKNOWLEDGMENTS
This work is submitted on behalf of the Long Life Family Study. This work was supported by the NIH: T32 GM74905 to BS; K01 AG057798 to SA; R01 AG061844 to PS, TP; U19 AG063893 to PS, TP, NS, SC; and the Boston University Digital Health Initiative Research Incubation Award to PS, TP, SA.
