Skip to main content

Probability of Alzheimer’s disease based on common and rare genetic variants

Abstract

Background

Alzheimer’s disease, among other neurodegenerative disorders, spans decades in individuals’ life and exhibits complex progression, symptoms and pathophysiology. Early diagnosis is essential for disease prevention and therapeutic intervention. Genetics may help identify individuals at high risk. As thousands of genetic variants may contribute to the genetic risk of Alzheimer’s disease, the polygenic risk score (PRS) approach has been shown to be useful for disease risk prediction. The APOE-ε4 allele is a known common variant associated with high risk to AD, but also associated with earlier onset. Rare variants usually have higher effect sizes than common ones; their impact may not be well captured by the PRS. Instead of standardised PRS, we propose to calculate the disease probability as a measure of disease risk that allows comparison between individuals.

Methods

We estimate AD risk as a probability based on PRS and separately accounting for APOE, AD rare variants and the disease prevalence in age groups. The mathematical framework makes use of genetic variants effect sizes from summary statistics and AD disease prevalence in age groups.

Results

The AD probability varies with respect to age, APOE status and presence of rare variants. In age group 65+, the probability of AD grows from 0.03 to 0.18 (without APOE) and 0.07 to 0.7 (APOE e4e4 carriers) as PRS increases. In 85+, these values are 0.08–0.6 and 0.3–0.85. Presence of rare mutations, e.g. in TREM2, may increase the probability (in 65+) from 0.02 at the negative tail of the PRS to 0.3.

Conclusions

Our approach accounts for the varying disease prevalence in different genotype and age groups when modelling the APOE and rare genetic variants risk in addition to PRS. This approach has potential for use in a clinical setting and can easily be updated for novel rare variants and for other populations or confounding factors when appropriate genome-wide association data become available.

Introduction

Genome-wide association studies (GWAS) identified genetic risk variants of late onset “sporadic” disease beyond the APOE locus [1,2,3,4], followed by exome chip analyses identifying rare variants with moderate risk [5,6,7]. While causal fully penetrant mutations almost certainly lead to development of the disease [8], most of the identified singular nucleotide polymorphism (SNP) risk alleles have not been proven to be causal but replicated as carrying an increased disease risk.

The PRSs are designed to aggregate genome-wide genotype data into a single variable indicating genetic liability to a disorder or trait. PRS studies often reach sufficiently high statistical significance to suggest trait polygenicity and, although the prediction accuracy is usually insufficient for clinical utility [9], PRS has been suggested as a useful tool for the selection for clinical trials of individuals of European ancestry across different traits [10,11,12,13]. The PRS prediction accuracy of risk for Alzheimer’s disease (AD) is comparatively high, especially when the diagnosis is based upon pathology confirmed rather than clinical assessment (AUC up to 84%) [14].

Designed to capture the risk of common variants, the PRS aggregates the effects of known genome-wide associated loci [15] and of loci that do not reach genome-wide statistical significance. However, the PRS may not well reflect the effect of rare variants in TREM2, PLCG2, ABI3, SORL1 [5,6,7] and very rare highly penetrant mutations in APP, PSEN1 and PSEN2 [16], as the cumulation of many SNPs of small effect sizes tends to mask the strong effect of a single variant. As LD between rare and common variants tends to be small—for example, the maximum r2 between a rare variant with minor allele frequency (MAF) = 0.01 and a typical common SNP with MAF = 0.2 is r2 = 0.04 when the rare alleles of both variants appear on the same haplotype (i.e. D’ = 1) [17]—rare variants and PRS are likely to be independent.

Variants in the APOE gene highly affect the AD risk (OR = 3.2, MAF = 0.14) [2]. APOE is also associated with lower odds of reaching the over 90th percentile age [18] as it modifies the age at onset; for example, the age at onset of AD for ε4ε4 carriers is ~ 68 years [19]. Also, ε4ε4 carriers are more likely to develop other conditions associated with lower life expectancy such as cardiovascular disease and diabetes [20]. People with AD diagnosed in their late 60s live on average 7 years after the clinical diagnosis, whereas AD diagnosis after age 90 is associated with an expected survival of only 2.8 years [21, 22]. Since age is the major confounding factor to the AD risk, it is difficult to disentangle the ageing and disease pathogenic components.

There is little research on whether APOE and PRS can be modelled as independent variables. Leonenko et al. [23] show that ε4 frequency decreases with age in both cases and controls, whereas the PRS values are higher in older AD patients, indicating a negative correlation in cases, but apparently not in controls. These contravariant effects in cases cancel out when a PRS is formed including APOE alongside other SNPs, so subsequent adjustment for age is ineffective. The use of APOE genotypes and the PRS (calculated without APOE) as two separate predictors accounts for this effect and increases the case/control prediction accuracy but cannot be extended to disease prediction in age groups with different disease prevalence. The approach suggested in the present study accounts both for the age related APOE effect and different disease prevalence. We propose estimating the disease probability (between 0 and 1) based on the PRS while accounting separately for high effect size variants and rare highly penetrant mutations. We show the utility of our probability calculations in application to AD.

Our calculations only require the mean and variance of PRS in cases and in controls and the disease prevalence as reference data, thus avoiding the need to share background sensitive genetic data. While the PRS from different studies (with different SNP selection and/or standardisation) cannot be directly compared, disease probability, as a general quantity, can be used for comparative prioritisation of individuals.

Material and methods

PRS distribution

The PRS aggregates the effects of multiple genetic markers identified by GWAS. Generally, the PRS is expected to be higher in cases than in controls, indicating a higher genetic risk for the disorder, but the difference in mean PRS between case and control samples may be small. It is important to note that the PRS calculated for an individual does not provide an absolute measure of risk and is meaningless except in relation with the distribution of PRS in cases and non-cases in the underlying population.

The polygenic risk score for individual j {1, …, Nind} is \( {PRS}_j=\frac{1}{N_{snps}}{\sum}_{i=1}^{N_{snps}}{g}_{ij}{\beta}_i \), where Nind and Nsnps are the numbers of individuals and of SNPs contributing to the PRS, respectively, gij {0, 1, 2} is the genotype of SNP i for individual j, and βi is the effect size (logarithm of the odds ratio or logistic regression coefficient) of SNP i in an independent GWAS for the disease. The sample mean and variance are

$$ m(PRS)=\frac{1}{N_{ind}}{\sum}_{j=1}^{N_{ind}}{PRS}_j\kern1em \mathrm{and}\kern1em \operatorname{var}(PRS)=\frac{1}{N_{ind}}{\sum}_{j=1}^{N_{ind}}{\left({PRS}_j-m(PRS)\right)}^2. $$
(1)

Estimation of PRS distribution parameters for unscreened controls

Our calculations require the distribution parameters of the PRS in cases and non-cases (putative non-affected controls). If the mean m0 and variance \( {\sigma}_0^2 \) of the PRS distribution in non-cases are unknown as unscreened population controls are used, they can be inferred from the means m1, mp and variances \( {\sigma}_1^2 \), \( {\sigma}_p^2 \) of the PRS distributions in cases and in the population, respectively, and the disease prevalence K as

$$ {m}_0=\frac{m_p-K\ {m}_1}{1-K} $$

and

$$ {\sigma}_0^2=\frac{\sigma_p^2-K{\sigma}_1^2}{1-K}-\frac{K{\left({m}_p-{m}_1\right)}^2}{{\left(1-K\right)}^2} $$

(see Supplemental Note 1 for details).

Estimation of the probability of disease development

By Bayesian inversion, a raw probability \( \hat{P} \) to be affected by the disease can be inferred from an individual’s PRS value x and the distribution densities of PRS in cases, p1, and in controls, p0, as

$$ \hat{P}(x)=\frac{K\ {p}_1(x)}{K\ {p}_1(x)+\left(1-K\right){p}_0(x)} $$
(2)

However, \( \hat{P} \) cannot be directly interpreted as a probability of disease. Logistic regression from case/control samples gives the probability of disease in the logistic model

$$ P(x)=\frac{1}{1+{e}^{-\left(\alpha +\beta x\right)}} $$
(3)

with coefficients α, β arising as regression parameters from the maximum likelihood estimate. We use linear regression with the logit link function, taking as data the log odds ratio corresponding to (2),

$$ y=\log\ \frac{K\ {p}_1(x)}{\left(1-K\right){p}_0(x)}, $$
(4)

at every PRS value x and the joint probability density of PRS in the population, pp(x) = Kp1(x) + (1 − K)p0(x), as weight. For normal densities p1 and p0, the coefficients α, β can be expressed as

$$ \alpha =\log \frac{K{\sigma}_0}{\left(1-K\right){\sigma}_1}+\frac{1}{2}\left(\left({r}_0-1\right)K+\left(1-{r}_1\right)\left(1-K\right)\right)-{m}_p\beta, \beta =\frac{m_1-{m}_0}{\sigma_p^2}\left(K\left(1-K\right)\left(\frac{r_0+{r}_1}{2}-1\right)+K\frac{\sigma_1^2}{\sigma_0^2}+\left(1-K\right)\frac{\sigma_0^2}{\sigma_1^2}\right) $$
(5)

where \( {r}_1=\frac{\sigma_0^2+{\left({m}_1-{m}_0\right)}^2}{\sigma_1^2} \) , \( {r}_0=\frac{\sigma_1^2+{\left({m}_1-{m}_0\right)}^2}{\sigma_0^2} \) , mp = Km1 + (1 − K)m0, and\( {\sigma}_p^2=K\ {\sigma}_1^2+\left(1-K\right){\sigma}_0^2+K\left(1-K\right){\left({m}_1-{m}_0\right)}^2 \) (see Supplemental Note 2 for details).

Formulae (5) determine the parameters of the logistic probability model (3) from the disease prevalence and the parameters of the distribution of PRS in cases and non-cases, dispensing with the need to obtain or simulate individual genotypes and perform logistic regression on the resulting PRS. They rely on the assumptions that the PRS distributions are normal and that the raw probability (2) represents well the fraction of cases for any value of PRS. For validation, we compared the outcome of (5) with the following three procedures of increasing abstraction, (a) simulation of genotypes in HWE with given MAF in cases and in non-cases and logistic regression of the resulting PRS, (b) sampling from normal distributions for PRS in cases and in non-cases with parameters \( {m}_0,{m}_1,{\sigma}_0^2,{\sigma}_1^2 \) and logistic regression, (c) sampling from the population distribution pp and linear regression of the raw log odds ratio (4).

Inclusion of rare variants in the probability

The effects of rare genetic variants with high (or medium) disease penetrance may be obscured if modelled as part of PRS including a large number of other SNPs, and the fraction of correctly identified cases carrying a rare mutation will be small in a sample and have little influence on the overall prediction accuracy. Therefore, it seems better to account for them at the level of the disease probability. Suppose we have the logistic regression model for the probability of disease PPRS in terms of the PRS by formulae (3) and (5), excluding the rare variant from the calculation of the PRS. An individual with PRS value x who carries a rare genetic variant with intrinsic probability prare to cause the disease has, assuming the effects of the rare variant and of the polygenic risk are independent, the probability of disease

$$ P(x)={P}_{PRS}(x)+{p}_{rare}\left(1-{P}_{PRS}(x)\right)=\frac{1+{p}_{rare}\ {e}^{-\left(\alpha +\beta x\right)}}{1+{e}^{-\left(\alpha +\beta x\right)}} $$

where x is the PRS for the individual. For very rare variant alleles that do not affect the disease prevalence K in the population, the intrinsic probability can be estimated as

$$ {p}_{rare}=\frac{K\ \left( OR-1\right)}{K\ \left( OR-1\right)+1}, $$

where OR is the odds ratio (see Supplemental Note 3). The probability P(x) takes values between prare and 1, reflecting the liability of the rare variant to cause the disease even in absence of polygenic risk. In case of several rare variants with mutually independent effect and intrinsic probabilities prare, 1, …, prare, ν, the above formula can be applied with \( {p}_{rare}=1-\prod \limits_{j=1}^{\nu}\left(1-{p}_{rare,j}\right) \). However, due to the assumption of very small allele frequencies, it is unlikely that an individual would carry more than one independent rare variant.

Inclusion of APOE

It may be advantageous to treat a high-effect common variant such as APOE separately from the PRS. The distributions in cases and non-cases of a PRS formed from SNPs excluding APOE can be assumed to be approximately equal for carriers and non-carriers of the APOE risk allele. Considering formulae (5), the probability of disease as a function of PRS will then differ between the groups only due to the higher disease prevalence in carriers of the risk allele. Applying (5) with the disease prevalence for the different APOE genotypes, separate probability curves are obtained. The prevalence in different genotype groups is not usually directly available but can be inferred as follows from the overall prevalence K, the overall allele frequency f and the odds ratio OR for the variant, under the assumption of HWE both in the general population and in the subpopulation of non-cases. This assumption is justified when the disease prevalence in the population is low (e.g. 2% for AD), but problematic when it is high [24] (e.g. major depression 30%). The prevalence K0, K1 and K2 for carriers of non-risk homozygotes, heterozygotes and risk homozygotes, respectively, can be calculated as

$$ {K}_0=1-\frac{{\left(1-f-\nu \right)}^2}{\left(1-K\right){\left(1-f\right)}^2},{K}_1=1-\frac{\left(1-f-\nu \right)\left(\nu +f-K\right)}{\left(1-K\right)f\left(1-f\right)},{K}_2=1-\frac{{\left(\nu +f-K\right)}^2}{\left(1-K\right){f}^2}, $$

where

\( \nu =\frac{b-\sqrt{b^2-16\left(1- OR\right)\left(1-f\right)K}}{4\ \left(1- OR\right)} \) with b = 2 (1 + (1 − OR)(K − f)) (see Supplemental Note 4).

Standardisation of the probability curve

PRSs calculated from different sets of SNPs cannot be directly compared. We therefore standardise the PRS axis by expressing the PRS in terms of standard deviations difference from the population mean,

$$ {x}_{st}=\frac{x-{m}_p}{\sigma_p}, $$

where x is the PRS and xst is the standardised PRS variable.

Simulated and real data

Firstly, we simulated independent genotypes in a sample of 10,000 cases and 10,000 controls and used previously published effect sizes for genome-wide significant SNPs [2, 15]. We calculated an Oligogenic Risk Score (ORS) in the simulated sample using only 39 genome-wide significant SNPs (Supplementary Table 2, adopted from [15]), excluding the APOE proxy SNP (rs429358). The PRS was calculated for 10,039 SNPs, including the above 39 genome-wide significant SNPs and further 10,000 SNPs pruned for LD with r2 = 0.1 and allele frequencies and effect sizes taken from (2).

Secondly, to illustrate the probability of disease in the presence of rare variants, we used effect sizes for rare variants corresponding to the APP, SORL1, TREM2, ABI3 and PLCG2 genes [6, 7, 25]. We used the distribution parameters m0, m1 and \( {\sigma}_0^2 \),\( {\sigma}_1^2 \) for ORS and PRS as reported in [23] and calculated the disease probabilities with the suggested formulae. To demonstrate the APOE modelling with our approach, we also took the distribution parameters of APOE, ORS and PRS from the real case/control study [23] (Supplemental Table 1).

The simulations and probability calculations were implemented with R-statistical software. The codes (Simulations.R and Probability.R) can be downloaded from https://github.com/DRI-Cardiff/AD-probability/.

Results

As the validity of the formulae (5) was established by simulations (see Supplemental Figure 1), we used formulae (3) and (5) to calculate the probability of disease for an individual with PRS value x. This probability depends on the disease prevalence in the population of interest, e.g. the general population or a specific subpopulation. The prevalence of AD in the population depends strongly on age. Recent estimates show a 3%, 17% and 33% prevalence in the 65–74, 75–84 and 85+ age groups, respectively [26]. For illustration, we calculated the probability of AD for 2%, 10% and 30% prevalence during lifetime and in 65+ and 85+ age groups, respectively. The parameters of the PRS distributions were taken from a real case/control study [23].

Figure 1 shows the dependency of the AD probability (y-axis) on standardised PRS (x-axis). The solid thick line corresponds to PRS.AD, calculated as weighted sum of APOE and PRS.noAPOE with the relative weight of APOE not taken directly from combining the corresponding effect sizes as a part of PRS but from bivariate logistic regression using APOE and PRS.noAPOE as predictors. The dashed line shows the probability of AD with PRS calculated in the standard way including all SNPs weighted with their single-SNP effect sizes. In the black and blue scenarios (corresponding to the lifetime and 65+ prevalences), the standard PRS shows clearer discrimination between low and high probabilities than APOE alone, somewhat similar to ORS. When the disease prevalence is high (red scenario), then PRS.AD is considerably more discriminative than PRS. This dramatic difference between PRS.AD and PRS is due to the fact that the means and variances for the latter are calculated in cases and controls, ignoring the change in APOE-ε4 allele frequency due to age, whereas PRS.AD indirectly accounts for it via the interplay of APOE and PRS.noAPOE.

Fig. 1
figure1

Probability of Alzheimer’s disease, ignoring change in APOE-ε4 allele frequency due to age. ORS - Oligogenic Risk Score including SNPs with a p-value threshold pT ≤ 10-5. PRS - Polygenic Risk Score including SNPs with a p-value threshold pT ≤ 0.1. PRS.AD - PRS calculated as a weighted sum of PRS.no.APOE, including SNPs with a pT ≤ 0.1 excluding APOE region (CHR19:44.4-46.5), and APOE (ε2 + ε4), where APOE effects were weighted with effect sizes (B(ε2) = -0.47 and B(ε4) = 1.12) as in Kunkle et al. 2019 [2]

Figure 2 demonstrates the results of adding independent rare variant effects in SORL1 and TREM2 to the probability of the disease in both the general population (K = 0.02) and in the age group 65+ (K = 0.1). Both graphs show an elevated disease probability (solid lines), with lower values in the population (black) and higher values in the 65+ group. Dashed lines show the disease probability across the range of PRS if individuals have no rare risk variants. The results for rare variants such as in genes APP, PSEN1 and PSEN2 [27] are shown in Supplemental Figure 2.

Fig. 2
figure2

Combined probability of AD calculated with 2% lifetime prevalence of AD, 10% prevalence of AD in 65+ age group, including probability due to presence of a rare high-effect variant. The rare variant effect sizes and minor allele frequencies correspond to known variants in SORL1 (OR = 7.2) [25] and TREM2 (OR = 2.46) [7] genes. In age group 65+ (red) the presence of SORL1 mutations (left) increases the AD probability from ~ 0 to 0.4 when PRS is the lowest and from 0.6 to 0.76 when the PRS is highest (solid line vs dashed line). For TREM2 (right), these values are 0 to 0.13 (low PRS) and 0.6 to 0.66 (high PRS)

Finally, Figure 3 shows the probability of disease in early onset (left) and late onset (right) age groups. As expected, the late onset group shows elevated AD probability even if the PRS is low and APOE ε4 carriers show consistently higher probability than any other genotypes. The black dashed line shows the disease probability with PRS when the APOE region is excluded. It is slightly higher than for APOE-ε4 non-carriers (thin blue line) as excluding the APOE region removes the information whether the individuals have lower disease risk due to absence of ε4 or have protective ε2 alleles.

Fig. 3
figure3

Modelling APOE separately, assuming the same effect of APOE-ε4 (OR ~ 3) in all age groups, and accounting for age related differences in frequency of APOE-ε4 allele (MAF = 0.18 in 55+ and MAF = 0.05 in 85+). In age group 65+ (left) the presence of APOE-ɛ4 allele increases the AD probability from 0.01 to 0.07 when PRS is the lowest and from 0.2 to 0.66 when the PRS is highest (top vs bottom lines). For 85+ age group (right), these values are 0 .06 to 0.3 (low PRS) and 0.65 to 0.92 (high PRS)

Discussion

PRS do not directly indicate an individual’s liability to develop a disease, as they depend on a variety of study parameters such as the number and selection of SNPs included in their calculation and are therefore not comparable between different studies. For comparability, it is not sufficient to standardise the PRS against the population mean and variance; the difference in PRS means in cases and in non-cases is also essential for the interpretation of an individual’s PRS. Our proposed calculation of a probability of disease takes all of these properties of the PRS into account and provides a unified measure to assess a PRS value in view of the PRS distribution and disease prevalence in the population or subpopulation of interest. Of course, the choice of SNPs included in a PRS remains decisive for its indicative power.

While the probability curve as a function of standardised PRS can be calculated by logistic regression using genotyped case/control samples from the population of interest, we offer a method to achieve the same outcome using more easily available summary data. The theoretical formula (5) derives the parameters for the logistic probability function (3) from the disease prevalence and the PRS distribution parameters (mean and standard deviation) in cases and non-cases. This conveniently allows estimating an individual’s probability of disease from their PRS value using only a small set of parameters. We have shown that (5) gives a highly accurate proxy for case-control sampling of PRS and logistic regression if either the PRS aggregates a high number of SNPs or very highly associated SNPs (such as APOE variants for AD) are excluded from the PRS. If the PRS is calculated from a small number of SNPs including some that are highly associated, deviation from normality in the PRS may lead to some discrepancy, but the formula still gives a close approximation.

Moreover, the probability of disease framework allows for separate consideration of high-effect variants. It has been observed that inclusion of high-effect variants in the PRS calculation does not always give optimal results. The impact of common high-effect variants such as APOE for AD can vary due to confounders such as age stratification, so taking a summary effect size from a reference study may result in a suboptimal PRS; this becomes apparent when APOE and PRS are used as separate predictors in bivariate logistic regression. We propose including high-effect variants in two different ways in the calculation of disease probability.

The effect of rare, highly penetrant genetic variants tends to be masked by more common variants in the PRS. However, if they do not affect the disease prevalence in the population and act independently of the PRS, we can account for them directly as adding a certain intrinsic probability for carriers of the risk allele. The probability curve approaches this intrinsic probability for highly negative values of standardised PRS.

Common high-effect variants such as APOE for AD strongly affect the disease prevalence in the population and cannot be dealt with as above. Although APOE may not act independently from other genetic causal variants, a recent study [28] suggests that it is reasonable to assume that the distribution parameters for the PRS calculated without APOE are independent of the APOE genotype. We propose calculating separate probability curves for each APOE status, based on the disease prevalence calculated for each APOE status from the disease prevalence in the population, the risk allele frequency and the odds ratio. These data are available and can easily be further stratified into relevant subpopulations, e.g. by age for AD.

The disease probability allows comparison of PRS calculated from summary data of different reference studies and thus has the potential to be used in a clinical context to prioritise individuals for diagnostics and preventative intervention based on assessed risk of developing the disease.

There are other factors that can influence the disease development probability. For example, it has been reported that AD is more prevalent in women and PRS effects may also depend on gender [29]. Due to the lack of information on gender-interactions for genetic variants, incorporating gender information in the probability calculations may not be straightforward at present. However, in view of the emerging literature it is likely that this information will be available and reliable in future. It can then be incorporated in our calculations, e.g. by selecting the SNPs for the PRS in males and females separately and/or changing the disease prevalence not by age only, but by sex as well. It is possible to include other, non-genetic predictors. This can be achieved by adjusting the disease prevalence in different e.g. educational attainment groups.

Our approach can be used for other complex genetic disorders. For example, schizophrenia is a highly polygenic disorder [30] and has an increased burden of rare variants and CNVs [31]. It is, however, a neurodevelopmental disorder, and the disease prevalence does not depend on age. It also does not have strong genetic risk factors like APOE for AD. Therefore, the most relevant probability calculation approach for diseases like schizophrenia and depression is a combination of common and rare variants, while the method shown above for the inclusion of APOE-like variants is not required.

The proposed method relies on the availability of allele frequency and genetic effect size estimates derived from a representative reference population. Expanding this approach to other populations will be possible when the effect sizes of SNPs in other populations will be reliably identified and reported. They can then be used to estimate the disease probability in the relevant population. As in other complex genetic disorders, the disease risk estimates rely heavily on the individual SNP risk estimates and disease prevalence, which differ depending on the demographics, ethnicity and age groups.

Limitation

A limitation of the present study is that it is based on sound, but theoretical principles, uses SNP and PRS characteristics from the literature and employs simulated data for validation. It remains to test and validate the theory directly in real datasets. Since we are dealing with rare variants, validation of this approach in real data with a sufficient level of confidence will require large population datasets like e.g. the UK BioBank [32]. However, the UK Biobank is not directly suitable to study neurodegenerative disorders, as the cohort is relatively young and only a small proportion of individuals manifest the disease. In addition, it does not provide phenotypic variables which are used to assess cognitive decline in dementia in clinical settings. Publicly available AD-specific datasets are typically small and not suitable to extract a reliable number of people carrying rare mutations (e.g. Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (www.loni.ucla.edu/ADNI)). Nevertheless, the approach presented here can be used in small studies focusing on rare and common genetic variants, for example to identify individuals most at risk of developing the disease. It can easily and flexibly be updated as novel rare variants are discovered and as appropriate GWAS data become available for specific populations.

Conclusions

The proposed method gives an estimate of the probability of developing AD based on an individual’s PRS, APOE genotype and the presence or absence of rare genetic variants associated with AD. The computational framework uses as reference data the means and standard deviations of the PRS employed in cases and in (either screened or population) controls and the disease prevalence. The disease prevalence varies considerably in different age and APOE genotype groups, and the present method allows for taking these differences into account in a natural and transparent way.

Availability of data and materials

The relevant data and materials are available in the supplemental material.

References

  1. 1.

    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. https://doi.org/10.1038/s41588-018-0311-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    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. https://doi.org/10.1038/s41588-019-0358-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    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. https://doi.org/10.1038/ng.2802.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Marioni RE, Harris SE, Zhang Q, McRae AF, Hagenaars SP, Hill WD, et al. GWAS on family history of Alzheimer’s disease. Transl Psychiatry. 2018;8(1):99. https://doi.org/10.1038/s41398-018-0150-6.

    Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. TREM2 Variants in Alzheimer’s Disease. N Engl J Med. 2013;368(2):117–27.

  6. 6.

    Jonsson T, Stefansson H, Steinberg S, Jonsdottir I, Jonsson PV, Snaedal J, et al. Variant of TREM2 associated with the risk of Alzheimer’s disease. N Engl J Med. 2013;368(2):107–16. https://doi.org/10.1056/NEJMoa1211103.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Sims R, van der Lee SJ, Naj AC, Bellenguez C, Badarinarayan N, Jakobsdottir J, et al. Rare coding variants in PLCG2, ABI3, and TREM2 implicate microglial-mediated innate immunity in Alzheimer’s disease. Nat Genet. 2017;49(9):1373–84. https://doi.org/10.1038/ng.3916.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Tanzi RE, Bertram L. Twenty years of the Alzheimer’s disease amyloid hypothesis: a genetic perspective. Cell. 2005;120(4):545–55. https://doi.org/10.1016/j.cell.2005.02.008.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Lewis CM, Vassos E. Prospects for using risk scores in polygenic medicine. Genome Med. 2017;9(1):96. https://doi.org/10.1186/s13073-017-0489-y.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Agerbo E, Sullivan PF, Vilhjalmsson BJ, Pedersen CB, Mors O, Borglum AD, et al. Polygenic Risk Score, Parental Socioeconomic Status, Family History of Psychiatric Disorders, and the Risk for Schizophrenia: A Danish Population-Based Study and Meta-analysis. JAMA Psychiatry. 2015;72(7):635–41. https://doi.org/10.1001/jamapsychiatry.2015.0346.

    Article  PubMed  Google Scholar 

  11. 11.

    Escott-Price V, Shoai M, Pither R, Williams J, Hardy J. Polygenic score prediction captures nearly all common genetic risk for Alzheimer’s disease. Neurobiol Aging. 2017;49:214 e7- e11.

    Article  Google Scholar 

  12. 12.

    Lall K, Magi R, Morris A, Metspalu A, Fischer K. Personalized risk prediction for type 2 diabetes: the potential of genetic risk scores. Genet Med. 2017;19(3):322–9. https://doi.org/10.1038/gim.2016.103.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Natarajan P, Young R, Stitziel NO, Padmanabhan S, Baber U, Mehran R, et al. Polygenic Risk Score Identifies Subgroup With Higher Burden of Atherosclerosis and Greater Relative Benefit From Statin Therapy in the Primary Prevention Setting. Circulation. 2017;135(22):2091–101. https://doi.org/10.1161/CIRCULATIONAHA.116.024436.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    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. https://doi.org/10.1002/ana.24999.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Andrews SJ, Fulton-Howard B, Goate A. Interpretation of risk loci from genome-wide association studies of Alzheimer’s disease. Lancet Neurol. 2020;19(4):326–35. https://doi.org/10.1016/S1474-4422(19)30435-1.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Tandon A, Fraser P. The presenilins. Genome Biol. 2002;3(11):reviews3014.

    Article  Google Scholar 

  17. 17.

    Slatkin M. Linkage disequilibrium--understanding the evolutionary past and mapping the medical future. Nat Rev Genet. 2008;9(6):477–85. https://doi.org/10.1038/nrg2361.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Deelen J, Evans DS, Arking DE, Tesi N, Nygaard M, Liu X, et al. A meta-analysis of genome-wide association studies identifies multiple longevity genes. Nat Commun. 2019;10(1):3669. https://doi.org/10.1038/s41467-019-11558-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    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. https://doi.org/10.1038/nrneurol.2012.263.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    El-Lebedy D, Raslan HM, Mohammed AM. Apolipoprotein E gene polymorphism and risk of type 2 diabetes and cardiovascular disease. Cardiovasc Diabetol. 2016;15(1):12. https://doi.org/10.1186/s12933-016-0329-1.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Garre-Olmo J, Ponjoan A, Inoriza JM, Blanch J, Sanchez-Perez I, Cubi R, et al. Survival, effect measures, and impact numbers after dementia diagnosis: a matched cohort study. Clin Epidemiol. 2019;11:525–42. https://doi.org/10.2147/CLEP.S213228.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Tom SE, Hubbard RA, Crane PK, Haneuse SJ, Bowen J, McCormick WC, et al. Characterization of dementia and Alzheimer’s disease in an older population: updated incidence and life expectancy with and without dementia. Am J Public Health. 2015;105(2):408–13. https://doi.org/10.2105/AJPH.2014.301935.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Leonenko G, Baker E, Stevenson-Hoare J, Sierksma A, Fiers M, Williams J, et al. Identifying individuals with high risk of Alzheimer’s disease using polygenic risk scores. Nat Commun. 2021;12:4506.

  24. 24.

    Wang J, Shete S. Using both cases and controls for testing hardy-weinberg proportions in a genetic association study. Hum Hered. 2010;69(3):212–8. https://doi.org/10.1159/000289597.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Bellenguez C, Charbonnier C, Grenier-Boley B, Quenez O, Le Guennec K, Nicolas G, et al. Contribution to Alzheimer’s disease risk of rare variants in TREM2, SORL1, and ABCA7 in 1779 cases and 1273 controls. Neurobiol Aging. 2017;59:220 e1- e9.

    Article  Google Scholar 

  26. 26.

    Hebert LE, Weuve J, Scherr PA, Evans DA. Alzheimer disease in the United States (2010-2050) estimated using the 2010 census. Neurology. 2013;80(19):1778–83. https://doi.org/10.1212/WNL.0b013e31828726f5.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    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. https://doi.org/10.1038/gim.2015.117.

    Article  PubMed  Google Scholar 

  28. 28.

    Bakulski KM, Vadari HS, Faul JD, Heeringa SG, Kardia SLR, Langa KM, et al. Cumulative Genetic Risk and APOE epsilon4 Are Independently Associated With Dementia Status in a Multiethnic, Population-Based Cohort. Neurol Genet. 2021;7(2):e576. https://doi.org/10.1212/NXG.0000000000000576.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Fan CC, Banks SJ, Thompson WK, Chen CH, McEvoy LK, Tan CH, et al. Sex-dependent autosomal effects on clinical progression of Alzheimer's disease. Brain. 2020;143(7):2272–80. https://doi.org/10.1093/brain/awaa164.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Pardinas AF, Holmans P, Pocklington AJ, Escott-Price V, Ripke S, Carrera N, et al. Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat Genet. 2018;50(3):381–9. https://doi.org/10.1038/s41588-018-0059-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Rees E, O'Donovan MC, Owen MJ. Genetics of schizophrenia. Curr Opin Behav Sci. 2015;2:8–14. https://doi.org/10.1016/j.cobeha.2014.07.001.

    Article  Google Scholar 

  32. 32.

    Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779. https://doi.org/10.1371/journal.pmed.1001779.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

Dementia Research Institute (UKDRI supported by the Medical Research Council (UKDRI-3003), Alzheimer’s Research UK, and Alzheimer’s Society), Joint Programming for Neurodegeneration (JPND, MRC: MR/T04604X/1), DPUK (MR/L023784/2), MRC Centre for Neuropsychiatric Genetics and Genomics (MR/L010305/1).

Author information

Affiliations

Authors

Contributions

KMS and VEP wrote the manuscript and prepared figures and tables. All authors reviewed and approved the final version of the manuscript.

Corresponding author

Correspondence to Valentina Escott-Price.

Ethics declarations

Ethical approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have 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 Note 1 Estimate of distribution parameters with unscreened controls. Supplemental Note 2 Estimate of probability of disease by linear regression with logistic link function. Supplemental Note 3 Inclusion of rare variants in the probability. Supplemental Note 4 Inclusion of common variants of common with high effect. Supplemental Table 1 Mean and variance for the Alzheimer’s disease genetic risk scores. Supplementary Table 2 AD genome-wide significant SNPs (adopted from Andrewes et al 2020). Supplemental Figure 1 The comparison of the theoretical probability of disease using formulae (3), (5) with the probability of disease estimated by logistic regression of simulated genotypes. Supplemental Figure 2 Combined probability of AD calculated with 2% lifetime prevalence of AD, 10% prevalence of AD in 65+ age group, including probability due to presence of a rare variant. The effect size is set to OR=500 reflecting almost fully penetrant mutations in APP, PSEN1, and PSEN2 genes.

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

Escott-Price, V., Schmidt, K.M. Probability of Alzheimer’s disease based on common and rare genetic variants. Alz Res Therapy 13, 140 (2021). https://doi.org/10.1186/s13195-021-00884-7

Download citation

Keywords

  • Alzheimer’s disease
  • Disease risk
  • Polygenic risk score
  • Rare variants
  • APOE