Skip to main content

Considering the APOE locus in Alzheimer’s disease polygenic scores in the Health and Retirement Study: a longitudinal panel study

Abstract

Background

Polygenic scores are a strategy to aggregate the small, additive effects of single nucleotide polymorphisms across the genome. With phenotypes like Alzheimer’s disease, which have a strong and well-established genomic locus (APOE), the cumulative effect of genetic variants outside of this area has not been well established in a population-representative sample.

Methods

Here we examine the association between polygenic scores for Alzheimer’s disease both with and without the APOE region (chr19: 45,384,477 to 45,432,606, build 37/hg 19) at different P value thresholds and dementia. We also investigate the addition of APOE-ε4 carrier status and its effect on the polygenic score—dementia association in the Health and Retirement Study using generalized linear models accounting for repeated measures by individual and use a binomial distribution, logit link, and unstructured correlation structure.

Results

In a large sample of European ancestry participants of the Health and Retirement Study (n = 9872) with an average of 5.2 (standard deviation 1.8) visit spaced two years apart, we found that including the APOE region through weighted variants in a polygenic score was insufficient to capture the large amount of risk attributed to this region. We also found that a polygenic score with a P value threshold of 0.01 had the strongest association with the odds of dementia in this sample (odds ratio = 1.10 95%CI 1.0 to 1.2).

Conclusion

We recommend removing the APOE region from polygenic score calculation and treating the APOE locus as an independent covariate when modeling dementia. We also recommend using a moderately conservative P value threshold (e.g. 0.01) when creating polygenic scores for Alzheimer’s disease on dementia. These recommendations may help elucidate relationships between polygenic scores and regions of strong significance for phenotypes similar to Alzheimer’s disease.

Peer Review reports

Background

The most common form of dementia is Alzheimer’s disease (AD), represented in roughly 65% of dementia cases [1]. Alzheimer’s disease is thought to arise from a combination of both genetics, environment, and lifestyle factors [2]. The estimated heritability of late onset Alzheimer’s disease is around 74% [3]. While large-scale genome-wide association studies (GWAS) have identified several genetic loci associated with Alzheimer’s disease [4,5,6,7,8,9,10], being a carrier of the Apolipoprotein E (APOE-ε4) allele remains the strongest genetic predictor of late-onset Alzheimer’s disease [11]. One copy of APOE-ε4 (inheriting a CC at these two locations from either parent) confers a threefold risk of Alzheimer’s disease while two copies (inheriting a CC at these two locations from both parents) a 15-fold increase in risk [12]. The effect of APOE-ε4 is all the more difficult to capture in a single variant GWAS as APOE-ε4 is a haplotype composed of two SNPs: rs7412 and rs429358—which will never be fully be captured in a traditional linear model GWAS framework. However, GWAS have identified many independent SNPs in and near the APOE gene locus. The APOE gene region contains many variants in high linkage disequilibrium within roughly 100 kilobases, including several additional high-risk sites in the translocase of outer mitochondrial membrane 40 (TOMM40) gene.

The largest Alzheimer’s disease GWAS meta-analysis to date (N = 94,437) is the from the International Genomics of Alzheimer’s Project (IGAP) [6]. This meta-analysis confirmed 20 previously identified Alzheimer’s disease risk loci [4] and identified five new genome-wide loci including (IQCK, ACE, ADAM10, ADAMTS1, and WWOX). The IGAP used a three-stage strategy where Stage 1 consisted of genotyped and imputed data on 9,456,058 common and 2,024,574 rare single nucleotide polymorphisms (SNPs) to meta-analyze GWAS from four cohorts (ncases = 21,982; ncontrols = 41,944). Stage 2 included replication with a custom I-select genotyping chip developed in Lambert et al. 2013 [4] and included 11,632 variants and 18,845 individuals with a meta-analysis of Stage 1 and Stage 2. Finally, Stage 3 replicated 44 variants and meta-analyzed Stages 1 and 2 and 3 for a total of 35,274 cases and 59,163 controls. The associations between millions of genetic loci and Alzheimer’s disease are documented in IGAP and available for testing in independent populations.

Many complex diseases may result from the consideration of small individual effects across the genome. Polygenic scores (PGS)—also known as polygenic risk scores/PRS, though they do not always model a risky phenotype (e.g. “risk” of high education, “risk” of increased height)—are generally derived from the sum of weighted variants across an individual [13, 14]. PGS incorporate genome-wide genetic variation into a single, quantitative measure that can be used in modeling as a tool to assess susceptibility. Though conceptually simple, many analytic decisions contribute to different qualities of PGSs including different coefficient of variation (R2), correlations between scores, and areas under the curve. An important consideration for diseases such as Alzheimer’s disease—which have a genetic locus like APOE/TOMM40 conferring much of the genetic risk to the disease—is to determine how the remaining variants in the genome contribute to the disease. Removing a region with many risk variants and deciding which and how many variants to include in a PGS can offer substantively different conclusions. For instance, one study reported a PGS area under the curve of 0.57 for Alzheimer’s disease (parental proxy) using 21 SNPs and excluding the APOE region [7], while another study reported using more than 200,000 variants (including APOE) and a PGS area under the curve of 0.84 for Alzheimer’s disease [15]. While Alzheimer’s disease has a strong genetic locus, there is no benchmark across metrics of PGS construction for the rest of the genome’s polygenic contribution.

Population-based studies often assess dementia status, rather that AD, as a trade-off between feasibility for longitudinal, larger samples of more diverse participants versus more specific and intensive clinical assessments. An Alzheimer’s disease PGS may be informative for dementia more broadly, and Alzheimer’s disease PGS has not been assessed in population-based studies of dementia, other than by proxy in the UK Biobank cohort [5]. Further PGS construction metrics regarding the APOE region and additional SNPs have not been compared. The goals of this manuscript are three-fold. The first aim is to assess the utility of using Alzheimer’s disease PGSs in population-based analyses of dementia. The second aim is to evaluate the inclusion of the APOE region in these PGSs with and without a covariate modeling risk directly from the APOE-ε4 allele. The third aim is to test SNP inclusion thresholds in PGS on dementia. We conduct this analysis using the Health and Retirement Study (HRS) in the European ancestries (intentionally plural, as there is no single “European” ancestry) participants.

Methods

Health and Retirement Study

The Health and Retirement Study (HRS) is a nationally representative panel study featuring a biennial survey of adults over age 50 and their spouses in the United States [16]. The HRS is sponsored by the National Institute on Aging (NIA U01AG009740) and is conducted by the University of Michigan. The HRS began in 1992 as a means to provide a national resource for data on changing health and economic circumstances associated with ageing at both the individual and population levels. These changes are focused on four broads topics: income and wealth; health, cognition, and use of healthcare services; work and retirement; and family connections [17].

The HRS pre-selected a random one-half of the sample to receive an enhanced face-to-face interview in 2006, which included physical performance tests, anthropometric measurements, blood and saliva samples, and a psychosocial self-administered questionnaire in addition to the HRS core interview. The remaining one-half sample received the same enhanced face-to-face protocol in 2008. The HRS randomly assigned the new 2010 cohort to receive an enhanced face-to-face interview in either 2010 or 2012. Those participants who were not interviewed or did not consent to saliva in 2006 were asked again in 2010.

Salivary DNA was collected using Oragene-250 saliva kits and protocol. DNA extracted from the saliva and was genotyped at the Center for Inherited Disease Research (CIDR) using the Illumina HumanOmni2.5 array (8v1 and 4v1). The Genetics Coordinating Center at the University of Washington, Seattle, WA performed Genotyping Quality Control. SNP annotation aligned to genome build 37/hg 19. The Genetics Coordinating Center calculated genetic principal components (PC) with HapMap controls [18, 19]. In addition to selecting independent SNPs with missing call rates < 5% and minor allele frequencies > 5%, the 2q21 (LCT), HLA, 8p23, and 17q21.31 regions were excluded from the initial pool [20]. The final European ancestries sample includes all self-reported non-Hispanic White persons that had PC loadings within ± one standard deviations of the mean for eigenvectors one and two in the PC analysis of all unrelated study subjects. The HRS re-calculated genetic PCs within the group of non-Hispanic Whites of European ancestries to further account for population stratification. These “ancestry-specific principal components” were used in subsequent analyses. Imputation was performed using IMPUTE2 on HRS data phased using SHAPEIT2. Data were imputed to the 1000 Genomes Project (1000GP) cosmopolitan reference panel phase 3 version 5 (initial release on May 2013, haplotypes released Oct 2014) and are available on the database of genotypes and phenotypes (dbGaP, https://www.ncbi.nlm.nih.gov/gap/,phs000428.v2.p2).

APOE-ε4

The apolipoprotein E (APOE) gene codes for a protein that binds and transports low-density lipids and is responsible, in part, for removing cholesterol from the bloodstream [21, 22]. Variations in this gene affect cholesterol metabolism and may lead to increases in the risk for stroke, heart disease, and may alter the odds of having Alzheimer’s disease. Two variants (rs7412 and rs429358) define the APOE genotype resulting in three common isoforms of protein apoE: E2, E3, and E4 encoded by ε2, ε3, and ε4. Of note, the genotyped rs7412 and rs429358 variants failed genotyping quality control in the HRS pipeline and are therefore not included as individual variants in any PGS using genotyped data alone. Using the imputed rs7412 (IMPUTE2 INFO score = 0.99) and rs429358 (IMPUTE2 INFO score = 0.99) variants, we categorized HRS participants as ε2/ε2, ε2/ε3, ε2/ε4, ε3/ε3, ε3/ε4, and ε4/ε4. We analyze two indicator variables for presence of one ε4 allele (1 = yes, 0 = no) or two ε4 allele (1 = yes, 0 = no).

Polygenic score for Alzheimer’s disease

We investigate using genome-wide raw genotyped variants and the two imputed APOE variants (rs7412 and rs429358) in the creation of our PGSs. We include SNPs in our PGS analysis at six AD-SNP association P-value thresholds (pT) from the Kunkle et al. [6] summary statistics: pT = (0.001, 0.01, 0.05, 0.1, 0.3, 1.0). For example, a PGS at pT = 0.01 includes only those variants for which the association p-value in the Kunkle et al. [6] GAP meta-analysis on Alzheimer’s disease was less than 0.01. We obtained summary statistics from National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site https://www.niagads.org/datasets/ng00075. Importantly, our study sample was not included in the Kunkle study of clinical Alzheimer’s disease. Thus, the weights are independent of our study sample. We do not include any linkage disequilibrium filtering in our scores (i.e. no clumping or pruning algorithms). That is, we include any directly genotyped variants that pass quality control from the HRS and overlap with those variants from the Alzheimer’s disease GWAS summary statistics in our scores, unless otherwise noted. We did not perform clumping or pruning because we are using the set of tag-SNPs from the genotyping chip, which do not contain imputed variants and represent, in essence, an already filtered set of variants.

For our scores with the APOE gene region removed, we removed all variants from the summary statistics on chromosome 19 near APOE (45,384,477 to 45,432,606, build 37/hg 19). This represents the start position of TOMM40 (45,394,477) − 10 kilobases and the stop position of APOC1 (45,422,606) + 10 kilobases. This region was removed in its entirety due to the dense linkage disequilibrium block in European ancestries overlapping these three genes (TOMM40, APOE, APOC1). We compare two PGSs: genotyped PGS including the APOE gene region, genotyped PGS with APOE region removed (Fig. 1), at six p-values thresholds from the Kunkle et al. [6] GWAS [(pT = (0.001, 0.01, 0.05, 0.1, 0.3, 1.0)].

Fig. 1
figure1

The APOE/TOMM40 genomic locus on chromosome 19. The y-axis corresponds to − log10(p-values) of association with Alzheimer’s disease in Kunkle et al. [6]. Single nucleotide polymorphisms within the bracketed genomic region were removed from consideration in polygenic score development for scores designed to exclude the APOE region. The variants that make the APOE isoforms are highlighted in purple diamond (rs7412, rs429358)

Cognition status

The HRS used a multidimensional measure of cognitive functioning, based on a telephone screening instrument: Telephone Interview for Cognitive Status [23]. Domains assessed using this measure include: memory, mental status, abstract reasoning, fluid reasoning, vocabulary, dementia, and numeracy. In 2009, Langa, Kabeto, and Weir developed an approach for defining dementia and cognitively impaired non-dementia (CIND) in the HRS. A team of dementia experts clinically validated this method using equipercentile equating in the HRS against the Aging, Demographics, and Memory Study (ADAMS). The ADAMS study is a sub-sample of the HRS who received a more extensive neurological battery [24, 25]. For self-respondents, the score consists of overall cognitive test performance while the proxy respondents’ scores are composed of proxy-rated memory, interviewer-perceived cognition, and impaired activities of daily living. The cut points for this method reflect the prevalence of dementia or cognitive impairment to the expected population prevalence from the ADAMS study. We used repeated measures of the classification of cognitive function data contributed for assessment years 2000−2014 from the HRS imputed cognition researcher contribution data set [26]. For self-respondents, a score from 0 to 6 is categorized as dementia, 7 to 11 is categorized as cognitive impaired not dementia, and 12 to 27 is categorized as normal cognition. For proxy respondents starting in 2000, a score of 6 or higher out of 11 is classified as having dementia, a score of 3 to 5 indicates cognitive impaired not dementia, while 0 to 2 indicates normal cognition [24]. In this analysis, we are only testing the odds of dementia versus normal cognition (1 = dementia, 0 = normal cognition).

Covariates

Educational attainment (years of school), birth cohort [(AHEAD: Asset and Health Dynamics Among the Oldest Old (b. < 1924); CODA: Children of the Depression (b. 1924–1930); HRS: Health and Retirement Study—original cohort (b. 1931–1941); WB: War babies (b. 1942–1947); EBB: early baby boomers (b. 1948–1953); MBB: mid-baby boomers (b. 1954–1959)], and sex (0 = female, 1 = male) are measured at a participant’s initial HRS exam. Age (years) and a self-report of doctor diagnosed stroke (0 = none, 1 = stroke, possible stroke/TIA/mini-strokes, respondent disputes previous waves that indicate condition) are assessed at the same wave as the corresponding cognition visit.

Statistical analysis

Due to the repeated measures in this analysis, we use generalized linear models and the GENMOD procedure in SAS 9.4. We specify repeated measures on the individual, a binomial distribution, and a logit link with an unstructured correlation structure. We include fixed effects covariates in every model: sex, years of education, and five ancestry-specific principal components. The time varying covariates chronological age, year, and stroke history at each interview wave are also included in every model. To evaluate the utility of the PGS in individuals with and without APOE-ε4, we also investigate an interaction effect between APOE-ε4 status and PGS. We use an α of 0.05 as a threshold for significance.

Results

There are 9,872 individuals in the HRS non-Hispanic White, European ancestries analytic sample collected between 2006 and 2010. We removed observations with missing cognition (m = 10,958), observations where the cognitive status was classified as CIND (m = 6905), and observations with a missing stroke status (m = 16). This removed 55 individuals from the analysis. The final analysis included n = 9817 h respondents of European ancestries with a total of m = 51,225 cognitive observations (Fig. 2).

Fig. 2
figure2

Sample inclusion flow chart for the Health and Retirement Study. HRS: Health and Retirement Study; n: number of individuals; m: number of observations; CIND: cognitively impaired, non-dementia

Our analytic sample was 57.8% female with an average age of 63.3 (SD = 10.3) at their first visit during the observation period. On average, participants had 13.2 (2.5) years of education. The average number of cognitive assessments per participant was 5.2 (1.8). A small portion of the sample had a history of stroke at their first visit during the observation period (n = 417, 4.3%), with a higher proportion of males reporting a history of stroke at their first visit (n = 201, 4.9%) than females (n = 216, 3.8%). In unadjusted analyses, the Alzheimer’s disease PGS at any pT was not significantly different between males and females and broadly centered at zero with a standard deviation of one. There was no difference in the distribution of number of copies of APOE-ε4 by sex, where the overall proportion with one copy of APOE-ε4 was 24.3% (n = 2382), and two copies of APOE-ε4 was 2.1% (n = 210). Across all observations, accounting for repeated measures within individuals, age and dementia at first visit during the observation period were not significantly different by sex. The proportion of stroke; however, was significantly different by sex (P = 0.0004) with a higher proportion in males than in females (Table 1). APOE-ε4 status is associated (all Ps < 0.05) with PGS (Fig. 3). Additional file 1: AF Table 1 contains a table of correlations between each Alzheimer’s disease PGS.

Table 1 Individual and observation-level descriptive statistics in the Health and Retirement Study, n = 9871, m = 51225
Fig. 3
figure3

Distribution of Alzheimer’s disease polygenic score, by number of copies of APOE-ε4 in the Health and Retirement Study, n = 9871. AD: Alzheimer’s disease; pT: p-value threshold for SNP-outcome association from the Alzheimer’s disease meta-analysis for inclusion into the polygenic score. APOE region defined as chromosome 19 (45,384,477 to 45,432,606, build 37/hg 19). This represents the start position of TOMM40 (45,394,477) − 10 KB and the stop position of APOC1 (45,422,606) + 10 KB

In adjusted, repeated measures regression models, increased age, later year of observation, history of stroke, and lower education were associated with increased odds of dementia in all models (all P < 0.0001; Table 2). In the models using an Alzheimer’s disease PGS with the APOE region included, having one copy of APOE-ε4 increased the odds of dementia relative to normal cognition by roughly twofold, while two copies of APOE-ε4 increased the odds of dementia by over fourfold (Table 2), holding all other variables constant. The Alzheimer’s disease PGS with the APOE region included was marginally associated with the odds of dementia relative to normal cognition only for the pT cutoff of 0.01 (OR = 1.1 95%CI 1.0 to 1.2), while the Alzheimer’s disease PGS at all other P-value thresholds was not significantly associated with the odds of dementia. In the models that included the Alzheimer’s disease PGS with the APOE region removed, we observed similar estimates for one and two copies of APOE-ε4 as before, and a similar effect size of ~ 8% increase in the odds of dementia for every one standard deviation increase in the Alzheimer’s disease PGSs. Both the Alzheimer’s disease PGS (without the APOE region) at pT = 0.001 and 0.01 were significantly associated with the odds of dementia (pT = 0.001 OR = 1.1 95%CI 1.0 to 1.2 and pT = 0.01 OR = 1.1 95%CI 1.0 to 1.2). The Alzheimer’s disease PGSs at pT = (0.05, 0.1, 0.3, and 1.0) were not associated with the odds of dementia, relative to normal cognition in this sample.

Table 2 Odds ratios from repeated measures generalized linear models regression analyses of dementia in the Health and Retirement Study, n = 9871, m = 51,225

Sensitivity analysis

Researchers have demonstrated some amount of mortality selection in the oldest individuals in the HRS genetic sample [27], and so we removed the two oldest cohorts (AHEAD and CODA) as a sensitivity analysis (Additional file 1: AF Table 2). We removed a total of 2020 individuals and 8957 observations. The effect of APOE-ε4 remained highly associated with the odds of dementia compared to normal cognition. Though the effect size for the Alzheimer’s disease PGSs at each pT and whether or not the APOE region was included were not significantly different from those in Table 2, the P-values associated with these effects were non-significant for all Alzheimer’s disease PGS. The slightly attenuated effects are not surprising as the younger cohorts are just now entering ages at which dementia becomes more prevalent. We additionally regressed out the effect of the APOE-ε4 from each PGS and performed the same analyses as in Table 2 (Additional file 1: AF Table 3). The PGSs with APOE-ε4 regressed out correlated with their counterparts at a Pearson’s correlation coefficient r > 0.97 (Additional file 1: AF Table 1). As expected, the results from Table 2 and these new analyses in Additional file 1: AF Table 3 are virtually identical. Associations between each PGS and dementia without adjustment for APOE-ε4, and also the association between the three level APOE-ε4 variable and dementia without adjustment for PGS are found in Additional file 1: AF Tables 4 and 5, respectively. Each of these models are adjusted for age, sex, education, year of visit, and stroke history. We found no statistical interaction between APOE-ε4 and a PGS at a P-value threshold of 0.01, with [(ORPGS*APOE(1) = 0.99 95%CI 0.83 to1.17; ORPGS*APOE(2) = 1.05 95%CI 0.69 to 1.59] or without [(ORPGS*APOE(1) = 1.00 95%CI 0.85 to 1.19]; [ORPGS*APOE(2) = 1.12 95%CI 0.74 to 1.68] the APOE region (Additional file 1: AF Tables 6).

Discussion

In a large, population-based cohort of older, European ancestries Americans, cumulative genetic risk summarized as a PGS is informative of longitudinal dementia odds. We observed that the APOE region requires handling with care in the development of PGS. Specifically, including the APOE region as weighted SNPs in a PGS was insufficient to account for the large risk attributed to the APOE region. We recommend removing the region in linkage disequilibrium around the APOE locus from the PGS and treating the APOE locus as an independent covariate. In addition, we observed greater performance from PGS developed at a P-value threshold of 0.01 for SNP inclusion, with greater noise from a PGS informed by the full genome in association with this dementia phenotype. Optimized measures of the polygenic nature of dementia allow for more powerful interrogations of genetic and environmental risk for dementia.

We observed the APOE-ε4 allele was longitudinally associated with higher risk of dementia, in a dose dependent manner. This observation is consistent with extensive prior research [28, 29]. The APOE-ε4 allele is neither necessary, nor sufficient to cause dementia, but the magnitude of increased risk attributed to each copy of the allele is relatively high. The APOE-ε4 allele is in linkage disequilibrium with a ~ 100 kilobase region involving the APOE, APOEC, and TOMM40 genes. Thus, an APOE independent PGS would need to remove the SNPs from the entire APOE region. In excess of the association between the APOE locus and dementia, we observed a small, but significant association between Alzheimer’s disease PGS and dementia. We also found that the effect of the PGS on dementia was not significantly different by APOE status. These findings are similar to those observed in clinical populations investigating APOE independent PGS risk of Alzheimer’s disease specifically [15, 30,31,32,33]. When building PGS, it is important to have independent study samples between the discovery GWAS and the application PGS. Notably, our study sample was not part of the Kunkle GWAS that generated the weights for the PGS. Our findings show the APOE independent Alzheimer’s disease PGS can be successfully implemented in population-based research of a broad dementia phenotype.

Dementia is a disorder with a strong genetic locus of effect (APOE) and substantially weaker effects are scattered throughout the genome. Including the APOE region in PGS without specific measurement of APOE-ε4 is insufficient, and overestimates the polygenic nature of dementia. Similarly, in Amyotrophic Lateral Sclerosis there is a strong main effect locus (C9orf72). Further, there is a significant, albeit modest, proportion of the phenotypic variance explained by a polygenic risk score over and above the C9orf72 region [34]. In contrast, other chronic disease traits, such as obesity, lack a dominant genetic locus and polygenic score development is successful across the entire genome at a higher p-value threshold [35]. Together, these results suggest that in traits with a strong genetic locus, polygenic scores should exclude the primary regions and seek to aggregate the remaining genetic risk as a separate predictor.

We acknowledge several limitations in this research. The first is that our study relied on imputed APOE variant calls. The APOE region is notoriously difficult and labor intensive to measure genotypes [36, 37]. Indeed, the two primary APOE SNPs of interest failed quality control metrics on the genotyping array in the HRS. We used the correlation structure of the genome from the 1000GP reference to impute these SNPs with ~ 99% confidence. Second, our study may be subject to mortality selection; however, dementia is primarily a disease of older age and requires survival long enough to manifest symptoms. Mortality selection related to the APOE genotype would only serve to make our observations more conservative. Third, we developed our PGS using weights from a GWAS of primarily European ancestry participants, thereby limiting generalizability to other ancestries. Last, our population-based study assessed a broad phenotype of dementia. There are many types of dementia including Alzheimer’s, vascular, and frontotemporal lobe, which have varying genetic architectures, to which we applied a PGS specific for Alzheimer’s. As future GWAS become available for dementia subtypes in a clinical population, investigators may be able to classify the utility of PGS in dementia subtypes. Future family-based studies may also consider SNPs strongly associated with early-onset AD, including APP, PSEN1, or PSEN2 mutations.

Conclusion

Dementia has considerable risk attributed to genetic factors. The APOE region is the strongest locus associated with disease and many additional sites confer small risk effects. Incorporating genetic risk from many sites in a polygenic risk score is a useful metric for risk prediction and etiologic testing in epidemiologic research of complex traits [38]. Our findings demonstrate the APOE region should be removed prior to polygenic risk score development and treated as an independent factor in dementia analyses. More work is needed to assess polygenic scores for Alzheimer’s disease for clinical utility and prediction and in diverse ancestries.

Availability of data and materials

Phenotype and covariate data are publicly available through the Health and Retirement Study, public use dataset. Produced and distributed by the University of Michigan with funding from the National Institute on Aging (Grant Number NIA U01AG009740) [39]. https://hrs.isr.umich.edu/data-products. The data/analyses presented in the current publication have been deposited in and are available from the dbGaP database under dbGaP accession phs000428.v2.p2 [40]. Summary statistics for Kunkle et al. [6] are available through The National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (https://www.niagads.org/datasets/ng00075). The list of specific genetic variants included in each polygenic score during the current study are available from the corresponding author on reasonable request.

Abbreviations

AD:

Alzheimer’s disease

APOE:

Apolipoprotein E

CI:

Confidence interval

GWAS:

Genome-wide association study

HRS:

Health and Retirement Study

IGAP:

International Genomics of Alzheimer’s Project

OR:

Odds ratio

PGS:

Polygenic risk score

PC:

Principal component

pT:

P-value threshold

SNP:

Single nucleotide polymorphism

1000GP:

1000 Genomes Project

References

  1. 1.

    Van Cauwenberghe C, Van Broeckhoven C, Sleegers K. The genetic landscape of Alzheimer disease: clinical implications and perspectives. Genet Med. 2016;18(5):421–30.

    Article  PubMed  Google Scholar 

  2. 2.

    Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA, Berg S, et al. Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry. 2006;63(2):168–74.

    Article  PubMed  Google Scholar 

  3. 3.

    Gatz M, Pedersen NL, Berg S, Johansson B, Johansson K, Mortimer JA, et al. Heritability for Alzheimer’s disease: the study of dementia in Swedish twins. J Gerontol A Biol Sci Med Sci. 1997;52(2):M117–25.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45(12):1452–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S, et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet. 2019;51(3):404–13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Abeta, tau, immunity and lipid processing. Nat Genet. 2019;51(3):414–30.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Tosto G, Bird TD, Tsuang D, Bennett DA, Boeve BF, Cruchaga C, et al. Polygenic risk scores in familial Alzheimer disease. Neurology. 2017;88(12):1180–6.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Seshadri S, Fitzpatrick AL, Ikram MA, DeStefano AL, Gudnason V, Boada M, et al. Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA. 2010;303(18):1832–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Naj AC, Jun G, Beecham GW, Wang LS, Vardarajan BN, Buros J, et al. Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer’s disease. Nat Genet. 2011;43(5):436–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, Hamshere ML, et al. Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer’s disease. Nat Genet. 2009;41(10):1088–93.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Ward A, Crean S, Mercaldi CJ, Collins JM, Boyd D, Cook MN, et al. Prevalence of apolipoprotein E4 genotype and homozygotes (APOE e4/4) among patients diagnosed with Alzheimer’s disease: a systematic review and meta-analysis. Neuroepidemiology. 2012;38(1):1–17.

    Article  PubMed  Google Scholar 

  12. 12.

    Liu CC, Liu CC, Kanekiyo T, Xu H, Bu G. Apolipoprotein E and Alzheimer disease: risk, mechanisms and therapy. Nat Rev Neurol. 2013;9(2):106–18.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    International Schizophrenia C, Purcell SM, Wray NR, Stone JL, Visscher PM, O'Donovan MC, et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 2009;460(7256):748–52.

  14. 14.

    Martin AR, Daly MJ, Robinson EB, Hyman SE, Neale BM. Predicting Polygenic Risk of Psychiatric Disorders. Biol Psychiatry. 2019;86(2):97–109.

    Article  PubMed  Google Scholar 

  15. 15.

    Escott-Price V, Myers AJ, Huentelman M, Hardy J. Polygenic risk score analysis of pathologically confirmed Alzheimer disease. Ann Neurol. 2017;82(2):311–4.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Health and Retirement Study. Health and Retirement Study, core public use dataset. In: U01AG009740) UoMwfftNIoAgnN, editor. Ann Arbor, MI2019.

  17. 17.

    Sonnega A, Faul J, Ofstedal M, Langa K, Phillips J, Weir D. Cohort Profile: the Health and Retirement Study (HRS). Int J Epidemiol. 2014;43(2):576–85.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS Genet. 2006;2(12):e190.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Center for Inherited Disease Research. Quality Control Report for Genotypic Data - Health and Retirement Study Phase 1–3. University of Washington; 2013.

  21. 21.

    Mahley RW. Apolipoprotein E: cholesterol transport protein with expanding role in cell biology. Science. 1988;240(4852):622–30.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Huang Y, Mahley RW. Apolipoprotein E: structure and function in lipid metabolism, neurobiology, and Alzheimer's diseases. Neurobiol Dis. 2014;72 Pt A:3–12.

  23. 23.

    Brandt J, Spencer M, Folstein M. The telephone interview for cognitive status. Neuropsychiatry Neuropsychol Behav Neurol. 1998;1:111–7.

    Google Scholar 

  24. 24.

    Langa KM, Kabeto M, Weir D. Report on race and cognitive impairment using HRS in 2010. Alzheimer's disease facts and figures; 2009.

  25. 25.

    Langa KM, Plassman BL, Wallace RB, Herzog AR, Heeringa SG, Ofstedal MB, et al. The Aging, Demographics, and Memory Study: study design and methods. Neuroepidemiology. 2005;25(4):181–91.

    Article  PubMed  Google Scholar 

  26. 26.

    Langa KM, Weir D, Kabeto M, Sonnega A. Researcher Contribution: Langa-Weir Classification of Cognitive Function. In: University of Michigan IfSR, editor. Survey Research CenterNovember 2018.

  27. 27.

    Domingue BW, Belsky DW, Harrati A, Conley D, Weir DR, Boardman JD. Mortality selection in a genetic sample and implications for association studies. Int J Epidemiol. 2017;46(4):1285–94.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Tanzi RE, Bertram L. New Frontiers in Alzheimer’s Disease Genetics. Neuron. 2001;32(2):181–4.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Logue MW, Panizzon MS, Elman JA, Gillespie NA, Hatton SN, Gustavson DE, et al. Use of an Alzheimer’s disease polygenic risk score to identify mild cognitive impairment in adults in their 50s. Molecular psychiatry. 2019;24(3):421.

    Article  PubMed  Google Scholar 

  30. 30.

    Desikan RS, Fan CC, Wang Y, Schork AJ, Cabral HJ, Cupples LA, et al. Genetic assessment of age-associated Alzheimer disease risk: Development and validation of a polygenic hazard score. PLoS Med. 2017;14(3):e1002258.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Tan CH, Fan CC, Mormino EC, Sugrue LP, Broce IJ, Hess CP, et al. Polygenic hazard score: an enrichment marker for Alzheimer’s associated amyloid and tau deposition. Acta Neuropathol. 2018;135(1):85–93.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Cruchaga C, Del-Aguila JL, Saef B, Black K, Fernandez MV, Budde J, et al. Polygenic risk score of sporadic late-onset Alzheimer’s disease reveals a shared architecture with the familial and early-onset forms. Alzheimers Dement. 2018;14(2):205–14.

    Article  PubMed  Google Scholar 

  33. 33.

    Lupton MK, Strike L, Hansell NK, Wen W, Mather KA, Armstrong NJ, et al. The effect of increased genetic risk for Alzheimer’s disease on hippocampal and amygdala volume. Neurobiol Aging. 2016;40:68–77.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    van Rheenen W, Shatunov A, Dekker AM, McLaughlin RL, Diekstra FP, Pulit SL, et al. Genome-wide association analyses identify new risk variants and the genetic architecture of amyotrophic lateral sclerosis. Nat Genet. 2016;48(9):1043–8.

    Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Ware EB, Schmitz LL, Faul JD, Gard A, Mitchell C, Smith JA, et al. Heterogeneity in polygenic scores for common human traits. bioRxiv. 2017.

  36. 36.

    Radmanesh F, Devan WJ, Anderson CD, Rosand J, Falcone GJ, Alzheimer's Disease Neuroimaging I. Accuracy of imputation to infer unobserved APOE epsilon alleles in genome-wide genotyping data. Eur J Hum Genet. 2014;22(10):1239–42.

  37. 37.

    Zhong L, Xie YZ, Cao TT, Wang Z, Wang T, Li X, et al. A rapid and cost-effective method for genotyping apolipoprotein E gene polymorphism. Mol Neurodegener. 2016;11:2.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Maher BS. Polygenic Scores in Epidemiology: Risk Prediction, Etiology, and Clinical Utility. Curr Epidemiol Rep. 2015;2(4):239–44.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Health and Retirement Study, public use dataset. Produced and distributed by the University of Michigan with funding from the National Institute on Aging (grant number NIA U01AG009740). Ann Arbor, MI, (2019). https://hrs.isr.umich.edu/data-products

  40. 40.

    The NCBI dbGaP Database of Genotypes and Phenotypes. dbGaP accession phs000428.v2.p2. (2019) https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000428.v2.p2

Download references

Acknowledgements

The authors would like to thank the participants of the Health and Retirement Study. We would additionally like to thank the scientific community for helpful comments on our medRxiv posting.

Funding

Research reported in this publication was supported by the National Institute on Aging and the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Numbers (NIA R01 AG055406, NIA R01 AG067592 (EBW, KMB, JDF), NIA RF1 AG055654 (EBW, JDF), NIMHD R01 MD011716 and R01 MD011716-S1 (EBW, CMM)). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The funders did not have any role in the design, analysis, interpretation, or writing of the manuscript.

Author information

Affiliations

Authors

Contributions

EBW analyzed the data. EBW and KMB interpreted the data analysis and were major contributors in writing the manuscript. EBW, KMB, JDF, and CMM were involved in critical editing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Erin B. Ware.

Ethics declarations

Ethics approval and consent to participate

The Health and Retirement Study is sponsored by the National Institute on Aging (NIA U01AG009740) and is conducted by the University of Michigan, where written informed consent was approved by the Institutional Review Board. This analysis was exempt and not regulated as determined by the Institutional Review Board at the University of Michigan (HUM00128220).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Supplemental tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ware, E.B., Faul, J.D., Mitchell, C.M. et al. Considering the APOE locus in Alzheimer’s disease polygenic scores in the Health and Retirement Study: a longitudinal panel study. BMC Med Genomics 13, 164 (2020). https://doi.org/10.1186/s12920-020-00815-9

Download citation

Keywords

  • Polygenic score
  • Alzheimer’s disease
  • Dementia
  • Apolipoprotein E
  • P-value
  • Thresholding