Polygenic Analysis of Late-Onset Alzheimer’s Disease from Mainland China

Abstract

Recently, a number of single nucleotide polymorphisms (SNPs) were identified to be associated with late-onset Alzheimer disease (LOAD) through genome-wide association study data. Identification of SNP-SNP interaction played an important role in better understanding genetic basis of LOAD. In this study, fifty-eight SNPs were screened in a cohort of 229 LOAD cases and 318 controls from mainland China, and their interaction was evaluated by a series of analysis methods. Seven risk SNPs and six protective SNPs were identified to be associated with LOAD. Risk SNPs included rs9331888 (CLU), rs6691117 (CR1), rs4938933 (MS4A), rs9349407 (CD2AP), rs1160985 (TOMM40), rs4945261 (GAB2) and rs5984894 (PCDH11X); Protective SNPs consisted of rs744373 (BIN1), rs1562990 (MS4A), rs597668 (EXOC3L2), rs9271192 (HLA-DRB5/DRB1), rs157581 and rs11556505 (TOMM40). Among positive SNPs presented above, we found the interaction between rs4938933 (risk) and rs1562990 (protective) in MS4A weakened their each effect for LOAD; for three significant SNPs in TOMM40, their cumulative interaction induced the two protective SNPs effects lost and made the risk SNP effect aggravate for LOAD. Finally, we found rs6656401-rs3865444 (CR1-CD33) pairs were significantly associated with decreasing LOAD risk, while rs28834970-rs6656401 (PTK2B-CR1), and rs28834970-rs6656401 (PTK2B-CD33) were associated with increasing LOAD risk. In a word, our study indicates that SNP-SNP interaction existed in the same gene or cross different genes, which could weaken or aggravate their initial single effects for LOAD.

About the Authors:

Bin Jiao

Contributed equally to this work with: Bin Jiao, Xiaoyan Liu

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, China

Xiaoyan Liu

Contributed equally to this work with: Bin Jiao, Xiaoyan Liu

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, China

Lin Zhou

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China

Maggie Haitian Wang

Affiliation: Division of Biostatistics, School of Public Health and Primary Care, the Chinese University of Hong Kong, Shatin, Hong Kong Special Administrative Region

Yafang Zhou

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China

Tingting Xiao

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, China

Weiwei Zhang

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, China

Rui Sun

Affiliation: Division of Biostatistics, School of Public Health and Primary Care, the Chinese University of Hong Kong, Shatin, Hong Kong Special Administrative Region

Mary Miu Yee Waye

Affiliation: School of Biomedical Sciences, the Chinese University of Hong Kong, Shatin, Hong Kong Special Administrative Region

Beisha Tang

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, China, State Key Laboratory of Medical Genetics, Changsha, China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China

Lu Shen

* E-mail: [email protected]

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, China, State Key Laboratory of Medical Genetics, Changsha, China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China

Introduction

Alzheimer’s disease (AD) is a clinically complex neurodegenerative disorder, affecting up to 81.1 million people worldwide [1]. It is characterized by memory and other cognitive decline, a variety of neuropsychiatric symptoms and restriction in the activities of daily living [2]. According to the age of onset, AD was classified as early-onset AD (EOAD, equal or less than 65 years) and late-onset AD (LOAD, more than 65 years), while the latter is the most common type of AD. The pathogenesis of AD is complicated, mainly caused by genetics, environment and normal aging [3]. To date, there are still only three causative genes reported in familial EOAD patients, including presenilin 1 (PSEN1), presenilin 2 (PSEN2) and amyloid precursor protein (APP) genes [4]. As the majority of patients were LOAD, many genetic association studies have been conducted in recent years to uncover the genetic contributions to LOAD, but the only gene variant considered to be an established LOAD risk factor was the APOEε4 allele [5].

However, APOE was limited to account for approximately 50% of individuals with LOAD [6], indicating other genes could contribute to LOAD risk. Until now, at least two methods are suggested to identify the risk gene variants for LOAD, including Whole Exome Sequence (WES) and genome-wide association study (GWAS) [7,8]. The former is mainly to identify rare coding variants, such as rs75932628 in TREM2, rs145999145 in PLD3 and rs137875858 in UNC5C, which were recently recognized as risk variants for LOAD [7,9–11]. With regard to common variants for AD, since 2009, five large GWAS and one Meta-analysis have identified more than 20 loci significantly associated with LOAD. According to the possible role in the process of AD, these genes were classified as several groups: (1) Lipid metabolism: APOE CLU, ABCA7, SORL1; (2) Immune response: CR1, CD33, MS4A, EPHA1, CLU, ABCA7, HLA-DRB5/DRB1 and INPP5D; (3) Endocytosis: BIN1, PICALM, CD2AP, EPHA1, RIN3, SORL1, MEF2C and MADD [8,12–16].

However, in most cases, the identified single nucleotide polymorphisms (SNPs) have small or moderate effect sizes, and the proportion of heritability explained is quite modest. Like other complicated diseases, a polygenic analysis has been suggested to explain genetic contribution to the pathogenesis of the majority of LOAD cases. We hypothesized that SNPs might interact in subtle ways that led to substantially greater effects than the effect of any single SNP. Therefore, in this study, we adopted a series of statistical analysis methods to evaluate SNP-SNP interactions in a Chinese cohort consisting of 547 individuals. It was found that SNP-SNP interaction could weak or aggravate their single effect for LOAD. In addition, even if some variants had no effect on LOAD, their interactions could become significant effects associated with LOAD.

Materials and Methods

Sample subjects

A total of 547 subjects were recruited in this study, including 229 patients with LOAD (male 43.1%; age at onset: 75.2±5.0 years) and 318 controls (male 47.8%; age:71.6±2.5 years). All patients diagnosed as probable or definite AD met with the NINCDS-ADRDA criteria, and they were collected from the department of neurology, Xiangya Hospital. Analyses also included all unaffected individuals of matched geographical ancestry as healthy controls. The study was approved by the Ethics Committee of Xiangya Hospital, Central South University in China (equivalent to an Institutional Review Board). Written informed consent was obtained from all subjects (if the patient was no capacity to understand this study due to cognitive impairment, written informed consent was obtained from their legal guardians).

Genotyping methods

Genomic DNA was extracted from peripheral blood leukocytes of all subjects. We used a fluorometer to evaluate the quality and quantity of DNA, which were normalized to 50 ng/μl for sequencing. APOE genotype was identified through polymerase chain reaction (PCR), and the primer information was listed in Table A in S1 File. We sequenced all PCR products with BigDye terminator v3.1 sequencing chemistry and performed on an ABI 3730xl DNA analyzer (Applied Biosystems). Finally, the DNA sequences were read by Sequencher software.

We then selected a list of AD candidate genes based on published GWAS paper and Alzgene.org (http://www.alzgene.org/) [8,12–16]. Totally, 56 SNPs in 28 candidate genes were selected in this study, and the genotypes of 547 subjects were screened using MALDI-TOF mass array method. To identify the accuracy of this method, APOE genotypes (rs7412 and rs429358) were also screened, and the result of this method was consistent with that of Sanger sequencing, The 56 SNPs involved BIN1 (rs744373, rs7561528, rs6733839), CLU (rs2279590, rs11136000, rs1532278, rs9331888, rs9331896), ABCA7 (rs3764650, rs4147929), CR1 (rs6656401, rs3818361, rs6701713, rs6691117), PICALM (rs3851179, rs592297, rs541458, rs561655, rs10792832), MS4A (rs4938933, rs670139, rs1562990, rs983392), CD33 (rs3865444, rs3826656, rs12459419), CD2AP (rs9349407, rs10948363), EPHA1 (rs11767557, rs11771145), TOMM40 (rs2075650, rs157581, rs157580, rs8106922, rs11556505, rs1160985), LRAT (rs727153), TNK1 (rs1554948), ARID5B (rs2588969), GAB2 (rs10793294, rs4945261), PCDH11X (rs5984894) SORL1(rs11218343, rs668387) PTK2B (rs28834970), ATP7B (rs1801243), EXOC3L2 (rs597668), SLC24A4-RIN3 (rs10498633), ZCWPW1 (rs1476679), CELF1 (rs10838725), MEF2C (rs190982), NME8 (rs2718058), HLA-DRB5/DRB1 (rs9271192), FERMT2 (rs17125944), ECE1 (rs213045) and CASS4 (rs7274581) (Table A in S2 File).

Statistical methods

Hardy-Weinberg equilibrium was analyzed using SHEsis software [17]. Logistic regression analysis was used to test for associations between each SNP allele and LOAD risk after adjusting age and gender. We then put all positive SNPs into multivariable logistic regression model to evaluate the association between each SNP and LOAD susceptibility. All above statistics were analyzed using SPSS 17.0 version. The p value<0.05 was defined as statistical significance.

Finally, Lasso-multiple regression (LMR) method was used to identify SNP-SNP pairs’ interaction. The calculation involved two steps: In step one, the Lasso was used to select candidate SNPs with all pair-wise interaction terms in the regression model, which was conducted under the R package glnetmi. Ten-fold cross validation was performed in this step. Interaction terms were ranked by the absolute value of the coefficients from high to low. Thus, it could achieve the aim of feature selection. In the second step, multiple linear regression was performed on the candidate SNPs selected in the first step to evaluate the pure interaction effect. The variable significance was measured by the p-value of cross-product term. Stepwise AIC selection was applied on the multiple linear regression to obtain an optimal regression model [18–20].

Results

The frequencies of APOE alleles and genotypes of all subjects were listed in Table 1. The distribution of ε4 allele in LOAD cases was significantly higher than that of controls (p = 0.002), and lower than that of controls for ε2 allele (p = 0.009). With regard to genotype, there was a significant difference in ε3/4 and ε4/4 genotypes between them (p = 0.004, p = 0.047) (Table 1).

[Figure omitted. See PDF.]

Table 1. The distribution of APOE genotypes and alleles in cases and controls.

https://doi.org/10.1371/journal.pone.0144898.t001

All candidate SNPs were in Hardy-Weinberg equilibrium except rs2075650 and rs7274581, therefore, we deleted both of them in the further analysis. Totally, 13 significantly different allelic frequencies were identified between patients and controls after adjusting age and gender, including seven risk SNPs and six protective SNPs for LOAD. The former consisted of rs9331888-C in CLU, rs6691117-A in CR1, rs4938933-C in MS4A, rs9349407-C in CD2AP, rs1160985-C in TOMM40, rs4945261-A in GAB2, rs5984894-A in PCDH11X, and the protective SNPs involved rs744373-T in BIN1, rs1562990-C in MS4A, rs597668-T in EXOC3L2, rs9271192-C in HLA-DRB5/DRB1, rs157581-G and rs11556505-T in TOMM40 (Table 2). To explore whether the identified seven risk SNPs or six protective SNPs had an interaction among them, we further analyzed them using multivariable logistic regression model. Finally, we found eleven SNPs had an independent effect for LOAD, except rs1160985 and rs157581 after adjusting age, gender and APOE ε4. (Tables A and B in S3 File).

[Figure omitted. See PDF.]

Table 2. Independent association of the 54 SNPs with LOAD.

https://doi.org/10.1371/journal.pone.0144898.t002

It was noted that both risk and protective SNPs were identified in MS4A (rs4938933-risk and rs1562990-protective) and TOMM40 (rs1160985-risk, rs157581-protective, and rs11556505-protective). With regard to two SNPs in MS4A, a total of four subhaplotypes were grouped, while no significant difference was identified between cases and controls (Table 3). In addition, eight subhaplotypes were combined within three SNPs in TOMM40, and the only subhaplotype G-T-C (rs157581-rs11556505-rs1160985) was significantly susceptible to LOAD risk. (Table 4)

[Figure omitted. See PDF.]

Table 3. Frequencies of subhaplotypes (rs4938933-rs1562990^*) in MS4A.

https://doi.org/10.1371/journal.pone.0144898.t003

[Figure omitted. See PDF.]

Table 4. Frequencies of subhaplotypes in TOMM40^*.

https://doi.org/10.1371/journal.pone.0144898.t004

We further used LMR method to analyze SNP-SNP pairs’ interaction: in the first step, 20 SNPs were selected as candidates from 56 SNPs (Tables A and B in S4 File); in the second step, three SNP-SNP pairs were found to be statistically significant after adjusting multiple testing by Bonferroni correction (Table 5). Three SNP pairs were identified to be significantly associated with LOAD (adjusted p-value to be 0.017, 0.018 and 0.032), including rs6656401-rs3865444 (CR1-CD33), rs28834970-rs6656401 (PTK2B-CR1), rs28834970- rs6656401 (PTK2B-CD33). The pair (rs6656401-rs3865444) with combination (AG, GT) reduced the LOAD risk by 0.10 compared to (GG, GG), while the other two SNP pairs had an increase risk for LOAD. The pair (rs28834970-rs6656401) with genotype (CT, AG) increased the LOAD risk by 10.90 compared to (TT, GG) and genotype (CT, GT) on (rs28834970-rs3865444) had an increased risk by 1.54 compared to (TT, GG).

[Figure omitted. See PDF.]

Table 5. The significant pairs were detected using LMR method.

https://doi.org/10.1371/journal.pone.0144898.t005

Discussion

LOAD was one of common diseases with strong genetic component [3]. Although APOE ε4 explained a portion of LOAD, other loci, especially identified through GWAS of LOAD, probably also participated in the process of LOAD. [21,22]. Researching these variants will lead to better understanding its biological function that can help in LOAD risk assessment, diagnosis and development of new therapies for LOAD. In this study, we conducted a comprehensive analysis of 56 SNPs using univariate analysis, multiple logistic regression and LMR methods. To our knowledge, it was the most comprehensive genetic analysis for LOAD patients from mainland China.

In current study, we found seven risk SNPs in CLU, CR1, MS4A, CD2AP, TOMM40, GAB2 and PCDH11X genes conferred to LOAD risk. CLU, known as apolipoprotein J probably increased AD risk through interacting with APOE [23]. A Meta-analysis showed rs9331888 was significantly associated with AD risk in Caucasian population [24], and we successfully replicated this risk loci. The role of CR1 in AD development has been highlighted due to involving in erythrocyte amyloid β42 sequestration and clearance of Aβ from whole blood [25], and the missense variant rs6691117 (Ile→Val) may change the folding of CR1 protein and affect structural stability through affecting Hsp70 binding, then cause functional changes [26]. The third identified risk gene was MS4A, which may be associated with the control of intracellular free Ca2+ concentration, resulting in neuronal death and decline in cognition [27]. The SNP rs4938933-C in MS4A was previously reported to be associated with decreasing LOAD risk in white population [15], while our result displayed rs4938933-C was a susceptible loci for LOAD risk. CD2AP was probably linked to modulating amyloid β clearance and tau neurotoxicity [28]. In this study, we first reported the SNP rs9349407 in CD2AP was significantly associated with LOAD in Chinese Han population, which was first identified to be significantly associated with AD in European ancestry, while the recent studies failed to replicate this result in Japanese, African-American and Canadian [29]. The TOMM40 gene is located adjacent to APOE. These two proteins may interact with each other to affect mitochondrial dynamics, although the precise mechanism underlying this is unclear. An interesting study showed the SNP rs1160985 was found to be associated with serum triglyceride concentration [30]. GAB2 is well characterized as a risk gene for the development of AD, which probably interacts with APOE ε4 to further modify risk [31]. With regard to PCDH11X, the detail of its biological mechanism was unclear, which needs us do more research to know about their relation with LOAD.

Until now, only two variants were identified to protect individuals away from AD: rs63750847-A in APP and APOEε2 genotype [32]. However, these two genes were causative or risk for AD. Therefore, we hypothesized some risk genes might also carry some protective variants. In current study, six SNPs were found to be associated with decreasing LOAD risk in BIN1, TOMM40, MS4A, EXOC3L2 and HLA-DRB5/DRB1. However, among them, five SNPs were previously identified risk loci except rs1562990-C in MS4A. Our result indicated that the impact of susceptible genes varied in different ethnicities, which could help us better understand the contributions of genetics to LOAD from Chinese Han population. In addition, more samples should be confirmed these protective variants.

No matter risk SNPs or protective SNPs, further analysis using multivariable logistic regression indicated most of significant loci might have an independent effect on LOAD, which confirmed AD was one of multifactorial diseases with great complex genetic background. Meanwhile, due to our positive result including different genes, it is worth to be focused on exploring their relation between amyloid-β pathology or tauopathy and these genes functions (lipid metabolism, immune response and endocytosis) in the future.

Another goal for current study was to identify whether SNP-SNP interactions existed in the pathogenesis of LOAD. We first analyzed the interactions between positive SNPs in the same gene. For rs4938933 and rs1562990 in MS4A, although they had independent reverse effect on LOAD, both of their respective effects were loss when being considered together. With regard to rs157581, rs11556505 and rs1160985 in TOMM40, although they were involved in decline or increasing risk for LOAD, the only significant subhaplotype was G-T-C that could increase LOAD risk, which indicates the SNPs effect on LOAD could be influenced by other SNPs in the same gene.

In addition, we further used LMR analysis method to identify three SNP pairs significantly associated with LOAD risk. However, the significant effect was not to be found in any of them during the previous univariate analysis. Among them, one pair rs28834970-rs6656401 from PTK2B-CR1 gene should be concerned due to the high Odds Ratio score (10.90) for LOAD risk, suggesting these two genes may have an strong interaction effect in the pathogenesis pathway. Previous study has identified two significant SNP-pairs for LOAD risk: rs386544-rs670139 from CLU-MS4A4E and rs11136000-rs670139 from CD33-MS4A4E. Taken these results together presented above, we offered evidence that SNP-SNP interaction played a pivotal role in LOAD susceptibility. As we know, AD is a complicated disease with gene-environment interaction. Therefore, although three pairs PTK2B-CR1, PTK2B-CD33, CD33-CR1 interaction appeared to have strong effects on LOAD, there may be unmeasured environment factors participating in these interaction. In current study, three new SNP-SNP pairs left their inter-molecular mechanism unsolved. CD33 is a transmembrane protein which has been implicated as a negative regulator of myeloid cells. CR1 is found on myeloid cells as well. Both of them are key molecules in inflammatory cells, such as microglia [33], therefore, we speculate they may share the similar pathway in the progression of LOAD. One AD research team from the United States now is beginning to map the molecular consequences of CR1 and CD33 variants to uncover their functionally link susceptibility loci [33]. With regard to PTK2B, which was one of recent identified new susceptible gene in 74046 individuals from diverse ethnicities [16], it was a key component of signaling pathways involved in neurite growth and synapse formation [34]. Although PTK2B was located on chromosome 8 near the CLU gene, we did not find a SNP-SNP interaction in these two genes. However, we found PTK2B had a strong interaction with NEDD9, which was important for lymphocyte signaling and migration and played a role in T cell-mediated inflammation [34]. Therefore, we speculate PTK2B might participate in immune inflammation through regulating NEDD9 protein function, which was involved in CR1 and CD33.

In summary, we examined 56 candidate SNPs in a cohort of Chinese subjects using a series of analysis methods. We identified seven risk and six protective variants for LOAD. With regard to SNP-SNP interaction, firstly, we found out some SNPs in the same gene could be influenced by each other to weaken or aggravate their single effect. Secondly, although some variants had a weak effect, a strong interaction effect was found after interacted with each other.

Supporting Information

[Figure omitted. See PDF.]

S1 File. Table A, The PCR condition of APOE genotyping.

https://doi.org/10.1371/journal.pone.0144898.s001

(DOCX)

S2 File. Table A, The details of each SNP genotyping in all subjects.

https://doi.org/10.1371/journal.pone.0144898.s002

(TXT)

S3 File. Table A, The multivariable logistic regression analysis in identified SNPs.

A multivariable logistic analysis in seven risk SNPs. Table B, A multivariable logistic regression analysis in the six protective SNPs.

https://doi.org/10.1371/journal.pone.0144898.s003

(DOCX)

S4 File. Table A, Identification of SNP-SNP interaction among all SNPs.

20 Candidate SNPs were selected in the first step using LMR method. Table B, Multiple linear regression model after stepwise selection.

https://doi.org/10.1371/journal.pone.0144898.s004

(DOCX)

Acknowledgments

We are grateful to all subjects for participation in our study.

Author Contributions

Conceived and designed the experiments: LS BT. Performed the experiments: BJ XL WZ TX. Analyzed the data: MHW YZ RS MMYW. Contributed reagents/materials/analysis tools: BJ LZ. Wrote the paper: BJ LS.

Citation: Jiao B, Liu X, Zhou L, Wang MH, Zhou Y, Xiao T, et al. (2015) Polygenic Analysis of Late-Onset Alzheimer’s Disease from Mainland China. PLoS ONE 10(12): e0144898. https://doi.org/10.1371/journal.pone.0144898

References

1. Ballard C, Gauthier S, Corbett A, Brayne C, Aarsland D, Jones E (2011) Alzheimer's disease. Lancet 377: 1019–1031. pmid:21371747

2. Balasa M, Gelpi E, Antonell A, Rey MJ, Sanchez-Valle R, Molinuevo JL, et al. (2011) Clinical features and APOE genotype of pathologically proven early-onset Alzheimer disease. Neurology 76: 1720–1725. pmid:21576687

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

4. Jiao B, Tang B, Liu X, Xu J, Wang Y, Zhou L, et al. (2014) Mutational analysis in early-onset familial Alzheimer's disease in Mainland China. Neurobiol Aging 35: 1957 e1951–1956.

5. Farrer LA, Cupples LA, Haines JL, Hyman B, Kukull WA, Mayeux R, et al. (1997) Effects of age, sex, and ethnicity on the association between apolipoprotein E genotype and Alzheimer disease. A meta-analysis. APOE and Alzheimer Disease Meta Analysis Consortium. JAMA 278: 1349–1356. pmid:9343467

6. Tanzi RE (2012) The genetics of Alzheimer disease. Cold Spring Harb Perspect Med 2: a006296. pmid:23028126

7. Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. (2013) TREM2 variants in Alzheimer's disease. N Engl J Med 368: 117–127. pmid:23150934

8. Lambert JC, Heath S, Even G, Campion D, Sleegers K, Hiltunen M, et al. (2009) Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet 41: 1094–1099. pmid:19734903

9. Cruchaga C, Karch CM, Jin SC, Benitez BA, Cai Y, Guerreiro R, et al. (2014) Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease. Nature 505: 550–554. pmid:24336208

10. Jiao B, Tang B, Liu X, Yan X, Zhou L, Yang Y, et al. (2014) Identification of C9orf72 repeat expansions in patients with amyotrophic lateral sclerosis and frontotemporal dementia in mainland China. Neurobiol Aging 35: 936 e919–922.

11. Wetzel-Smith MK, Hunkapiller J, Bhangale TR, Srinivasan K, Maloney JA, Atwal JK, et al. (2014) A rare mutation in UNC5C predisposes to late-onset Alzheimer's disease and increases neuronal cell death. Nat Med 20: 1452–1457. pmid:25419706

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

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

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

15. Hollingworth P, Harold D, Sims R, Gerrish A, Lambert JC, Carrasquillo MM, et al. (2011) Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease. Nat Genet 43: 429–435. pmid:21460840

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

17. Shi YY, He L (2005) SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res 15: 97–98. pmid:15740637

18. Wang MH, Li J, Yeung VS, Zee BC, Yu RH, Ho S, et al. (2014) Four pairs of gene-gene interactions associated with increased risk for type 2 diabetes (CDKN2BAS-KCNJ11), obesity (SLC2A9-IGF2BP2, FTO-APOA5), and hypertension (MC4R-IGF2BP2) in Chinese women. Meta Gene 2: 384–391. pmid:25606423

19. Lockhart R, Taylor J, Tibshirani RJ, Tibshirani R (2014) A SIGNIFICANCE TEST FOR THE LASSO. Ann Stat 42: 413–468. pmid:25574062

20. Friedman J, Hastie T, Tibshirani R (2010) Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw 33: 1–22. pmid:20808728

21. Jun G, Vardarajan BN, Buros J, Yu CE, Hawk MV, Dombroski BA, et al. (2012) Comprehensive search for Alzheimer disease susceptibility loci in the APOE region. Arch Neurol 69: 1270–1279. pmid:22869155

22. Jun G, Ibrahim-Verbaas CA, Vronskaya M, Lambert JC, Chung J, Naj AC, et al. (2015) A novel Alzheimer disease locus located near the gene encoding tau protein. Mol Psychiatry. [Epub ahead of print].

23. Bertrand P, Poirier J, Oda T, Finch CE, Pasinetti GM (1995) Association of apolipoprotein E genotype with brain levels of apolipoprotein E and apolipoprotein J (clusterin) in Alzheimer disease. Brain Res Mol Brain Res 33: 174–178. pmid:8774959

24. Shuai P, Liu Y, Lu W, Liu Q, Li T, Gong B, et al. (2015) Genetic associations of CLU rs9331888 polymorphism with Alzheimer's disease: A meta-analysis. Neurosci Lett 591: 160–165. pmid:25703218

25. Crehan H, Holton P, Wray S, Pocock J, Guerreiro R, Hardy J, et al. (2012) Complement receptor 1 (CR1) and Alzheimer's disease. Immunobiology 217: 244–250. pmid:21840620

26. Kullo IJ, Ding K, Shameer K, McCarty CA, Jarvik GP, Denny JC, et al. (2011) Complement receptor 1 gene variants are associated with erythrocyte sedimentation rate. Am J Hum Genet 89: 131–138. pmid:21700265

27. Ma J, Yu JT, Tan L (2014) MS4A Cluster in Alzheimer's Disease. Mol Neurobiol.

28. Shulman JM, Chen K, Keenan BT, Chibnik LB, Fleisher A, Thiyyagura P, et al. (2013) Genetic susceptibility for Alzheimer disease neuritic plaque pathology. JAMA Neurol 70: 1150–1157. pmid:23836404

29. Chen H, Wu G, Jiang Y, Feng R, Liao M, Zhang L, et al. (2015) Analyzing 54,936 Samples Supports the Association Between CD2AP rs9349407 Polymorphism and Alzheimer's Disease Susceptibility. Mol Neurobiol. [Epub ahead of print].

30. Salakhov RR, Goncharova IA, Makeeva OA, Golubenko MV, Kulish EV, Kashtalap VV, et al. (2014) TOMM40 gene polymorphism association with lipid profile. Genetika 50: 222–229. pmid:25711031

31. Reiman EM, Webster JA, Myers AJ, Hardy J, Dunckley T, Zismann VL, et al. (2007) GAB2 alleles modify Alzheimer's risk in APOE epsilon4 carriers. Neuron 54: 713–720. pmid:17553421

32. Jonsson T, Atwal JK, Steinberg S, Snaedal J, Jonsson PV, Bjornsson S, et al. (2012) A mutation in APP protects against Alzheimer's disease and age-related cognitive decline. Nature 488: 96–99. pmid:22801501

33. De Jager P, Bradshaw E, Chibnik L, Keenan B, Shulman J, Schneider J, et al. (2012) Neuropathologic and immunologic traits integrate the effect of CR1 and CD33 variants on Alzheimer's-related amyloid pathology Alzheimer's & Dementia: The Journal of the Alzheimer's Association 8: 232.

34. Beck TN, Nicolas E, Kopp MC, Golemis EA (2014) Adaptors for disorders of the brain? The cancer signaling proteins NEDD9, CASS4, and PTK2B in Alzheimer's disease. Oncoscience 1: 486–503. pmid:25594051