Abstract
Recently, several genome-wide association studies (GWASs) have led to the discovery of nine new loci of genetic susceptibility in Alzheimer's disease (AD). However, the landscape of the AD genetic susceptibility is far away to be complete and in addition to single-SNP (single-nucleotide polymorphism) analyses as performed in conventional GWAS, complementary strategies need to be applied to overcome limitations inherent to this type of approaches. We performed a genome-wide haplotype association (GWHA) study in the EADI1 study (n=2025 AD cases and 5328 controls) by applying a sliding-windows approach. After exclusion of loci already known to be involved in AD (APOE, BIN1 and CR1), 91 regions with suggestive haplotype effects were identified. In a second step, we attempted to replicate the best suggestive haplotype associations in the GERAD1 consortium (2820 AD cases and 6356 controls) and observed that 9 of them showed nominal association. In a third step, we tested relevant haplotype associations in a combined analysis of five additional case–control studies (5093 AD cases and 4061 controls). We consistently replicated the association of a haplotype within FRMD4A on Chr.10p13 in all the data set analyzed (OR: 1.68; 95% CI: (1.43–1.96); P=1.1 × 10−10). We finally searched for association between SNPs within the FRMD4A locus and Aβ plasma concentrations in three independent non-demented populations (n=2579). We reported that polymorphisms were associated with plasma Aβ42/Aβ40 ratio (best signal, P=5.4 × 10−7). In conclusion, combining both GWHA study and a conservative three-stage replication approach, we characterised FRMD4A as a new genetic risk factor of AD.
Similar content being viewed by others
Introduction
The identification of genes involved in monogenic forms of Alzheimer's disease (AD) has significantly contributed to our knowledge of the disease mechanisms. The causal links between mutations, the functions of the mutated genes (APP, PS1 and PS2) and disease development prompted a pathophysiological hypothesis, which radically changed our understanding of AD: the amyloid cascade hypothesis.1 The systematic association of pathogenic mutations with changes in APP metabolism and, more particularly, a relative overproduction of Aβ42 peptides indicates that this metabolism is at the heart of the disease process (at least in the monogenic forms of the disease). The overproduction of these neurotoxic peptides is supposed to lead to or accentuate neuron-to-neuron propagation of the τ pathology (leading to neuronal death) by an unknown mechanism.2
By analogy, it was expected that the characterization of genetic factors involved in the common forms of AD (that is, lacking classical Mendelian inheritance), the most frequent form of the disease, should also help to better understand the AD physiopathological process. However, the characterization of these genetic factors has encountered significant difficulties. Until 2009, the APOE (apolipoprotein E) gene was the only globally valid genetic determinant of AD to have been unambiguously identified in 15 years of intensive research.3, 4
As with other multifactorial diseases, this systematic inability to detect new genetic determinants has prompted more comprehensive investigations using genome-wide association studies (GWASs). We and others performed five large GWASs in this field and reported that the CLU (clusterin), PICALM (phosphatidylinositol-binding clathrin assembly protein), CR1 (complement component (3b/4b) receptor 1), BIN1 (bridging integrator 1), ABCA7 (ATP-binding cassette, sub-family A, member 7), MS4A (membrane spanning 4A) cluster, EPHA1 (ephrin type-A receptor 1), CD33 (differentiation antigen 33) and CD2AP (CD2-associated protein) genes were associated with the AD risk.5, 6, 7, 8, 9 Most of these susceptibility genes have been already systematically replicated in Caucasians in large case–control studies and in families.10
However, our understanding of the AD genetics is far away to be complete and strong efforts have still to be done. At this level, classical GWAS approaches present an important limitation with systematic application of a conventional, highly conservative Bonferroni correction leading to select only the most statistically significant associations (commonly, P<1 × 10−8). This involves the risk of rejecting biologically valid hypotheses on purely statistical grounds, that is, false negatives.
To partly handle some of this limitation, other complementary approach consists in extracting pertinent information from single-nucleotide polymorphisms (SNPs) nominally associated with the risk of developing AD in GWAS by using complex statistical and bioinformatics multiple-SNP analyses such as genome-wide haplotype association (GWHA) study.11 We adopted this GWHA strategy to the AD genetic susceptibility through a three-step approach.
Materials and methods
Population description
The main characteristics of the populations used for the GWHA study are described in the Supplementary Notes and Supplementary Table 1. All AD cases met the criteria for either probable AD (NINCDS-ADRDA, DSM-IV)12 or definite AD (CERAD).13 All elderly controls were screened for dementia using the MMSE or the AD Assessment Scale-cognitive subscale and were determined to be free from dementia at neuropathological examination or had a Braak score of 2.5 or below. All subjects or, in those with substantial cognitive impairment, a caregiver, legal guardian, or other proxy gave written informed consent for participation in this study. The study protocols for all populations were reviewed and approved by the appropriate Institutional review boards of each country.
The main characteristics of the populations (nondemented individuals) used for the Aβ plasma study (3C study, Rotterdam and CHS) are described in the Supplementary Notes and Supplementary Table 2.14, 15, 16
Genotyping
Participants in the French GWA study (including the 3C study) were genotyped using an Illumina 610–quad array (Illumina, San Diego, CA, USA). Quality control and analytical parameters have been described in detail elsewhere.5 The GERAD participants were genotyped using an Illumina 610–quad array (Illumina), a HumanHap550 array or a HumanHap300 array (Illumina). Again, QC and analytical parameters have been described in detail elsewhere.6 The Rotterdam and CHS studies were genotyped using the Affymetrix 500K array (Affymetrix, Santa Clara, CA, USA). Again, QC and analytical parameters have been described in detail elsewhere.7
In the European populations (stage 3), genotyping was performed using Sequenom assays at the exception of the German population genotyped with a 610K ILLUMINA chip (Illumina). The primer and probe sequences used in the genotyping assays are available upon request. In order to avoid any genotyping bias, cases and controls were randomly mixed while genotyping and laboratory personnel were blinded to case/control status. The genotyping success rate was at least 95%. Departure from Hardy–Weinberg equilibrium was observed for rs2446581 in the Swedish control and case samples (P=9.1 × 10−7 for the whole population). The Swedish sample was consequently excluded from further analysis, as haplotype analyses that do not comply with Hardy–Weinberg equilibrium are likely to bias observations.
Statistical analyses
Missing age or gender data
Any individuals with missing age or gender data were excluded. This gave a maximum of 2025 AD cases and 5328 controls in step 1, 2820 AD cases and 6356 controls in step 2 and 5093 AD cases and 4061 controls in step 3.
Detection of haplotype effects using a sliding-window approach
This approach has been fully described elsewhere.11, 17 Briefly, the search for haplotype effects was carried out by applying a sliding-window approach18, 19 to the French GWA data set for each chromosome. After excluding SNPs not in Hardy–Weinberg equilibrium or with a minor allele frequency (MAF)<0.02, the first step of the strategy was to eliminate part of the redundancy between SNPs by using haplotype-tagging SNPs (htSNPs). For this, the same binning procedure as in11 was used: within each bin of 10 adjacent SNPs, we identified a minimal set of htSNPs that were able to characterize more than 95% of the inferred haplotypes with estimated frequency greater than 0.02. Once a bin had been characterized by a set of htSNPs, the same strategy was applied to the bin composed of the next 10 adjacent SNPs. The final set of htSNPs (n=287 956) was then fed into the sliding-windows approach.
Given a window of 10 htSNPs, the search for the most informative and parsimonious haplotype configuration in terms of disease prediction was performed for all possible 1 to 4 loci combinations of not necessarily adjacent SNPs. We used a strategy based on Akaike's information criterion, which has been previously described for candidate gene haplotype analysis.19, 20 It relies on the stochastic expectation–maximisation module21 in THESIAS software.22 If required, missing genotypes were inferred by applying,23 multiple imputation.21 In all, 37 330 050 combinations were investigated in our genome scan and this investigation was conducted thanks to the use of the grid technology developed by the European Grid Infrastructure (http://www.egi.eu).20 This technology enables several thousand computations to be run in parallel on a large number of different CPUs. The sliding-window haplotype approach was developed into a GridHaplo grid package for the EGI grid (http://genecanvas.ecgene.net).11
Replication of haplotype effects
Regions with ‘window P-values’ (i) below 10−5 and (ii) 100 times smaller than the smallest single-locus P-value (including all SNPs and not only htSNPs) were analyzed in terms of replication in the GERAD1 data set by using THESIAS software, with systematic adjustment for age and gender.
Association of FRMD4A SNPs with Aβ plasma concentrations
In each center, Aβ plasma variables are normally distributed. We excluded prior analyses all samples with a +/− two s.d. values in order to avoid potential associations driven by extreme observations. Finally, each quantitative variable was transformed into a z-score (equal to (observed value minus the sample mean), divided by the sample s.d.). The association between the Aβ1−40, Aβ1−42 and Aβ1−42/Aβ1−40 z-scores on one hand and imputed FRMD4A imputed SNPs on the other (see below) were assessed using a general linear model under an additive model adjusted for age, center and gender.
We used inverse-variance weighting (also known as fixed-effects meta-analysis) to investigate the homogeneity of haplotype effects from one study to another and to provide meta-analysed, age- and gender-adjusted ORs for haplotype effect estimates in the seven studies. A similar strategy was used to provide meta-analysed, age- and gender-adjusted association levels between Aβ plasma z-scores and FRMD4A SNPs.
Imputation analyses
We imputed SNPs by using MaCH (http://www.sph.umich.edu/csg/abecasis/mach/index.html) and minimac software (http://genome.sph.umich.edu/wiki/minimac). The reference haplotype data are provided by the MaCH website, which was built for the combined Caucasian populations as part of the 1000 genomes project. In our data set, all individuals were genotyped on the same platform (the Illumina Human660W-Quad Beadchip) and we used 492 941 observed SNP genotypes that passed quality filters as follows: genotyping call rate ⩾98%, Hardy–Weinberg Equilibrium Test P value ⩾1 × 10−6 and MAF ⩾1%. We first inferred haplotype combinations of each individual using the ‘phase’ option in the MaCH program and then imputed them with minimac. As minimac is a newly developed software tool, we compared the correlation between the imputed genetic dosage from minimac and those from the standard MaCH program for SNPs in chromosome 22. The results were very similar (data not shown). Doses for 7 704 555 million SNPs with a MAF>0.01 were available from the French GWA data set using the 1000 genomes data set. We selected 2538 SNPs within the Chromosome 10p13 locus of interest (chromosome10:13 655 705–14 402 866) and evaluated their associations with AD risk in an additive logistic regression model adjusted for age, gender and disease status. A graphic representation was then generated with Locuszoom software (http://csg.sph.umich.edu/locuszoom/).
In the 3C, Rotterdam and CHS cohorts, we used the genotype data to impute to the 2.5 million non-monomorphic, autosomal SNPs described in HapMap II (CEU population) as described elsewhere.7 We selected 1486 SNPs (MAF>0.01) within the Chromosome 10p13 locus of interest (chr10:13 655 705–14 402 866) and evaluated their associations with Aβ plasma concentrations as described above. Again, a graphic representation was then generated with Locuszoom software (http://csg.sph.umich.edu/locuszoom/).
Results
We developed a three-step approach. In the first step, the French GWA study (EADI1 for European Alzheimer's initiative 1),5 including 2025 AD cases and 5328 controls was used to select regions with potential haplotype associations with AD. Following a sliding-windows approach (see Material and methods), we applied two a priori criteria to select loci of interest as previously described:11 level of association with a P-value (i) below 10−5 and (ii) at least 100 times smaller than the smallest single-SNP P value observed in the corresponding locus. We were able to detect loci already known to be involved in AD (the APOE, BIN1 and CR1 locus)5, 6, 7, 8, 9 from previous GWASs. The obtained signal was systematically highly stronger than the one observed for each SNP taken separately (Table 1). After exclusion of these loci, we retained 91 regions of interest.
In the second step, we replicated these haplotype associations in the GWA database from another consortium involved in the study of AD genetic susceptibility, the GERAD1 consortium,6 including 2820 AD cases and 6356 controls. All 91 regions were available for investigation and 9 of them showed nominal association (P<0.05; Table 2). We decided to further investigate two regions showing the same best haplotypes associated with AD risk in both EADI1 and GERAD1, with similar magnitude and direction of association (Table 2). These two loci, not previously detected in single-SNP GWAS analyses, were located on chromosomes 6p21 and 10p13. The identified AD-associated haplotypes at the 6p21 region were tagged by rs2395760, rs991762 and rs4711652. The rs7081208, rs2446581 and rs17314229 tagged the 10p13 haplotypes. For the Chromosome 6p21 locus, the best association was attributable to the GGT haplotype in both GWASs (OR (Odds Ratio): 1.53; 95% CI: (1.31–1.79); P=8.1 × 10−8 after adjustment for age and gender when both EADI1 and GERAD1 studies were combined). For the Chromosome 10p13 locus, the highest level of association was attributable to the AAC haplotype (OR: 1.76; 95% CI: (1.44–2.15); P=2.3 × 10−8 after adjustment for age and gender when both EADI1 and GERAD1 studies were combined. Additional adjustment for the four main principal components did not modify the results, data not shown) (Table 2).
In the third step, the six tagging SNPs were further genotyped in five additional AD case–control studies from Flanders-Belgium (842 cases and 489 controls), Finland (560 cases and 623 controls), Germany (728 cases and 961 controls), Italy (1846 cases and 904 controls) and Spain (1117 cases and 1084 controls) (Supplementary Tables 3 and 4). The same common haplotypes were inferred from the six SNPs in the two loci, with the exception of the AGT haplotype in Chromosome 10p13 in the Italian population (frequency<0.01 in controls) (Supplementary Table 5 and 6).
In a combined analysis of the five replication data sets, the overall difference in haplotype distribution for the 6p21 locus was not significant between AD cases and controls (P=0.13) and the GGT haplotype was not associated with AD risk (OR: 0.89; 95% CI: (0.74–1.08); P=0.23 after adjustment for age and gender). We therefore considered that the haplotype association in this locus was not confirmed (see Supplementary Table 5).
Conversely, an overall significant difference in haplotype distribution between AD cases and controls was observed at the 10p13 locus in the replication sample (P=4.3 × 10−2 after adjustment for age, gender and country) and the AAC haplotype was associated with increased AD risk in the five data sets meta-analysis (OR: 1.55; 95% CI: (1.19–2.00); P=9.2 × 10−4 after adjustment for age and gender). When the seven data sets (EADI1, GERAD1 and the follow-up studies) were analyzed together, the AAC haplotype had a meta-analysed OR of 1.68 (95% CI: (1.43–1.96); P=1.1 × 10−10 adjusted for age and gender) when compared with the most frequent GGC haplotype, with no evidence of heterogeneity across the seven countries (P=0.92) (Figure 1). Among the whole sample, nine cases were homozygous for the AAC haplotype versus only five controls (OR: 2.85; 95% CI: (0.88–9.76); P=0.09 with Yates correction), which is consistent with a dose-dependent effect of the AAC haplotype.
Importantly, none of the individual SNPs were associated with AD risk in the total sample with a P value lower than 10−5. The best meta-analysed SNP was rs2446581 (OR: 1.15; 95% CI: (1.08–1.24); P=2.1 × 10−5 after adjustment for age and gender; Supplementary Table 7). The hypothesis that the rs2446581 could solely explain the association observed with the 10p13 haplotypes was rejected by use of the likelihood ratio test (P=9.1 × 10−7).
We also tested whether this haplotype association might be explained by one or more untyped SNPs located nearby the genotyped SNPs. Genotypes for 7 704 555 million SNPs with a MAF⩾0.01 were imputed from the French GWA data by using the 1000 genome dataset (http://www.1000genomes.org/). None of the single SNPs imputed at this locus (n=2538) showed stronger evidence of association with the AD risk than the haplotypes initially identified (Supplementary Figure 1). The 10p13 haplotype region appeared to be fully included within the FERM domain containing 4A (FRMD4A) gene, as indicated by the linkage disequilibrium map of this locus of interest (Figure 2).
We finally explored how FRMD4A might be involved in the AD process. According to a recent report, the FRMD4A gene was described to interact with Arf6.24 As this latter protein was reported to control APP processing,25 we postulated that the FRMD4A locus could be associated with endophenotypes susceptible to reflect modulation of the APP metabolism. We accordingly analyzed association of the FRMD4A locus with Aβ peptide plasma concentrations in three independent populations of nondemented individuals (n=2579) (Supplementary Table 2) for which FRMD4A imputed SNPs, Aβ1−40 and Aβ1−42 plasma concentrations were available. Strong associations were observed between several FRMD4A SNPs and Aβ1−42/Aβ1−40 (nine SNPs reaching a significant level after Bonferroni correction, P=3.4 × 10−5 for 1486 SNPs, Figure 3). The best signal was obtained for rs7921545 (meta-analysed z-score β coefficient: 0.12, CI: 95% (0.07–0.017); P=5.4 × 10−7) and was homogeneous between the three data sets (P for heterogeneity=3.7 × 10−1, Supplementary Figure 2). Of note, rs2446581 showed nominal association with Aβ1−42/Aβ1−40 (P=3.1 × 10−2, see Supplementary Figure 2). Furthermore, only nominal associations between FRMD4A SNPs, Aβ1−40 or Aβ1−42 were detected (see Supplementary Figure 3).
Discussion
From a specific haplotype-based GWAS approach, we were able to detect a new genetic susceptibility factor for AD that could not be identified through usual GWAS analyses. Owing to the high number of association tests performed in this GWHA study (37 330 050, not all of which were independent), a robust replication strategy was necessary. Our three-step approach was thus particularly conservative and we cannot rule out the possibility that we failed to identify other haplotype-based loci associated with AD risk. Nonetheless, our GWHA immediately identified two potential AD susceptibility loci among which the 10p13 locus showed an AAC haplotype that was strongly and consistently associated with AD risk in seven independent European populations. Interestingly, this locus is included in the large AD linkage region regularly identified on chromosome 10.26
Additional work will be necessary to determine whether the AD risk-haplotype association indicates an interaction between SNPs or whether they tag non-genotyped, functional variants. However, our GWHA study could have handled several of the inherent limitations of GWAS. As the AAC haplotype is rare (with a mean frequency of 2% in our Caucasian populations), our work suggests the possibility that rare variants may be responsible for the signal detected within the FRMD4A gene. This might explain why the locus was not detected (i) in our previous GWA studies based on single-SNPs analyses5, 6, 7, 8, 9 and (ii) through imputation, as SNPs with low frequency and/or SNPs within specific regions with low LD are poorly imputed even when using the 1000 genome data set (Figure 2). Furthermore, we cannot rule out the possibility that the AAC haplotype tags insertion–deletion variants or copy-number variations—neither of which is captured by imputation. Interestingly, several copy-number variations within the FRMD4A locus have been described (http://projects.tcag.ca/variation/).
Little is known about the function of the protein encoded by FRMD4A in animal cells. It belongs to the FERM super family, which includes ubiquitous components of the cytocortex involved in cell structure, transport and signalling functions.27 According to a recent report, the FRMD4A gene product may regulate epithelial polarity by interacting with Arf6 and the PAR complex.24 Interestingly, several other lines of evidence suggest that Arf6 modulates cell polarity in various systems—including neurons. Arf6 reportedly regulates dendritic branching in hippocampal neurons and neurite outgrowth in PC12 cells.28, 29 Finally, Arf6 was recently reported to control APP processing,27 suggesting that FRMD4A could also be implied in this metabolism. This hypothesis is sustained by the observation of an association of the FRMD4A locus with plasma Aβ1−42/Aβ1−40, at once reinforcing the plausibility of the association of this gene with AD risk and its potential implication in a subtle control of the APP metabolism. Unfortunately, as the SNPs defining the AAC haplotype were only available by imputation in the Rotterdam and CHS study, search for association of the AAC haplotype with plasma Aβ1−42/Aβ1−40 was not possible. Furthermore, it is important to keep in mind that this result is difficult to interpret in terms of AD pathophysiological process. First, it is not known whether plasma Aβ peptides reflect a dynamic equilibrium between the brain, CSF and plasma compartments.30, 31, 32, 33, 34 Second, the source of plasma Aβ species is not known and the Aβ peptides' physiological functions are still not fully understood.35, 36
However, taken as a whole, these data suggest that FRMD4A could be a relevant candidate gene for AD risk and the basis for another possible pathophysiological pathway for AD.
References
Hardy J, Selkoe D . The amyloid hypothesis of Alzheimer's disease: progress and problems on the road to therapeutics. Science 2002; 297: 353–356.
Small SA, Duff K . Linking Abeta and tau in late-onset Alzheimer's disease: a dual pathway hypothesis. Neuron 2008; 60: 534–542.
Lambert JC, Amouyel P . Genetic heterogeneity of Alzheimer's disease: complexity and advances. Psychoneuroendocrinology 2007; 32: S62–S70.
Genin E, Hannequin D, Wallon D, Sleegers K, Hiltunen M, Combarros O et al. APOE and Alzheimer disease: a major gene with semi-dominant inheritance. Mol Psychiatry 2011; 16: 903–907.
Lambert JC, Heath S, Even G, Campion D, Sleegers K, Hiltunen M et al. Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet 2009; 41: 1094–1099.
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: 1088–1093.
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: 1832–1840.
Hollingworth P, Harold D, Sims R, Gerrish A, Lambert JC, Carrasquillo MM et al. Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease. Nat Genet 2011; 43: 429–435.
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: 436–441.
Lambert JC, Amouyel P . Genetics of Alzheimer's disease: new evidences for an old hypothesis? Curr Opin Genet Dev 2011; 21: 295–301.
Trégouët DA, König IR, Erdmann J, Munteanu A, Braund PS, Hall AS et al. Genome-wide haplotype association study identifies the SLC22A3-LPAL2-LPA gene cluster as a risk locus for coronary artery disease. Nat Genet 2009; 41: 283–285.
McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan EM . Clinical diagnosis of Alzheimer's disease: report of the NINCDS-ADRDA work group under the auspices of department of Health and Human services task force on Alzheimer's disease. Neurology 1984; 34: 939–944.
Mirra SS, Heyman A, McKeel D, Sumi SM, Crain BJ, Brownlee LM et al. The consortium to establish a registry for Alzheiemr's disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer's disease. Neurology 1991; 41: 479–486.
Van Oijen M, Hofman A, Soares HD, Koudstaal PJ, Breteler MM . Plasma Abeta(1-40) and Abeta(1-42) and the risk of dementia: a prospective case-cohort study. Lancet Neurol 2006; 5: 655–660.
Lambert JC, Schraen-Maschke S, Richard F, Fievet N, Rouaud O, Berr C et al. Association of plasma amyloid beta with risk of dementia: the prospective three-city study. Neurology 2009; 73: 847–853.
Lopez OL, Kuller LH, Mehta PD, Becker JT, Gach HM, Sweet RA et al. Plasma amyloid levels and the risk of AD in normal subjects in the cardiovascular health study. Neurology 2008; 70: 1664–1671.
Mathias RA, Gao P, Goldstein JL, Wilson AF . A graphical assessment of P-values from sliding window haplotype tests of association to identify asthma susceptibility loci on chromosome 11q. BMC Genet 2006; 7: 38.
Yu Z, Schaid DJ . Sequential haplotype scan methods for association analysis. Genet Epidemiol 2007; 31: 553–564.
Huang BE, Amos CI, Lin DY . Detecting haplotype effects in genome wide association studies. Genet Epidemiol 2007; 31: 803–812.
Tregouet DA, Ricard S, Nicaud V, Arnould I, Soubigou S, Rosier M et al. In-depth haplotype analysis of ABCA1 gene polymorphisms in relation to plasma ApoA1 levels and myocardial infarction. Arterioscler Thromb Vasc Biol 2004; 24: 775–781.
Tregouet DA, Escolano S, Tiret L, Mallet A, Golmard JL . A new algorithm for haplotype-based association analysis: the stochastic-EM algorithm. Ann Hum Genet 2004; 68: 165–177.
Tregouet DA, Garelle V . A new JAVA interface implementation of THESIAS: testing haplotype effects in association studies. Bioinformatics 2007; 23: 1038–1039.
Gagliardi F, Jones B, Grey F, Bégin ME, Heikkurinen M . Building an infrastructure for scientific grid computing: status and goals of the EGEE project. Philos Transact A Math Phys Eng Sci 2005; 363: 1729–1742.
Ikenouchi J, Umeda M . FRMD4A regulates epithelial polarity by connecting Arf6 activation with the PAR complex. Proc Natl Acad Sci USA 2010; 107: 748–753.
Sannerud R, Declerck I, Peric A, Raemaekers T, Menendez G, Zhou L et al. ADP ribosylation factor 6 (ARF6) controls amyloid precursor protein (APP) processing by mediating the endosomal sorting of BACE1. Proc Natl Acad Sci USA 2011; 108: E559–E568.
Hamshere ML, Holmans PA, Avramopoulos D, Bassett SS, Blacker D, Bertram L et al. Genome-wide linkage analysis of 723 affected relative pairs with late-onset Alzheimer's disease. Hum Mol Genet 2007; 16: 2703–2712.
Tepass U . FERM proteins in animal morphogenesis. Curr Opin Genet Dev 2009; 19: 357–367.
Hernández-Deviez DJ, Casanova JE, Wilson JM . Regulation of dendritic development by the ARF exchange factor ARNO. Nat Neurosci 2002; 5: 623–624.
Albertinazzi C, Za L, Paris S, de Curtis I . ADP-ribosylation factor 6 and a functional PIX/p95-APP1 complex are required for Rac1B-mediated neurite outgrowth. Mol Biol Cell 2003; 14: 1295–1307.
Mehta PD, Pirttila T, Patrick BA, Barshatzky M, Mehta SP . Amyloid beta protein 1-40 and 1-42 levels in matched cerebrospinal fluid and plasma from patients with Alzheimer disease. Neurosci Lett 2001; 304: 102–106.
Vanderstichele H, Van Kerschaver E, Hesse C, Davidsson P, Buyse MA, Andreasen N et al. Standardization of measurement of beta-amyloid(1-42) in cerebrospinal fluid and plasma. Amyloid 2000; 7: 245–258.
Kawarabayashi T, Younkin LH, Saido TC, Shoji M, Ashe KH, Younkin SG . Age-dependent changes in brain, CSF, and plasma amyloid (beta) protein in the Tg2576 transgenic mouse model of Alzheimer's disease. J Neurosci 2001; 21: 372–381.
Ghersi-Egea JF, Gorevic PD, Ghiso J, Frangione B, Patlak CS, Fenstermacher JD . Fate of cerebrospinal fluid-borne amyloid beta-peptide: rapid clearance into blood and appreciable accumulation by cerebral arteries. J Neurochem 1996; 67: 880–883.
Giedraitis V, Sundelöf J, Irizarry MC, Gårevik N, Hyman BT, Wahlund LO et al. The normal equilibrium between CSF and plasma amyloid beta levels is disrupted in Alzheimer's disease. Neurosci Lett 2007; 427: 127–131.
Soscia SJ, Kirby JE, Washicosky KJ, Tucker SM, Ingelsson M, Hyman B et al. The Alzheimer's disease-associated amyloid beta-protein is an antimicrobial peptide. PloS One 2010; 5: e9505.
Lambert JC, Dallongeville J, Ellis KA, Schraen-Maschke S, Lui J, Laws S et al. Association of plasma Aβ peptides with blood pressure in the elderly. PloS One 2011; 6: e18536.
Acknowledgements
The work was made possible by the generous participation of the control subjects, the patients and their families. This work was supported by the National Foundation for Alzheimer's disease and related disorders, the Institut Pasteur de Lille, the Centre National de Génotypage, Inserm, FRC (fondation sur la recherche sur le cerveau) and Rotary. The Three-City Study was performed as part of a collaboration between the Institut National de la Santé et de la Recherche Médicale (Inserm), the Victor Segalen Bordeaux II University and Sanofi-Synthélabo. The Fondation pour la Recherche Médicale funded the preparation and initiation of the study. The 3C Study was also funded by the Caisse Nationale Maladie des Travailleurs Salariés, Direction Générale de la Santé, MGEN, Institut de la Longévité, Agence Française de Sécurité Sanitaire des Produits de Santé, the Aquitaine and Bourgogne Regional Councils, Fondation de France and the joint French Ministry of Research/INSERM ‘Cohortes et collections de données biologiques’ programme. Lille Génopôle received an unconditional grant from Eisai. Experiments presented in this paper were carried out using the EGI European grid infrastructure, supported by EGI, NGI France and Université Lille 1. We highly appreciate and thank the technical staff of the CRI-Lille 1 center for their strong and helpful support. Belgium sample collection: the research at the Antwerp site was in part supported by the Interuniversity Attraction Poles program P6/43 of the Belgian Science Policy Office, the Foundation for Alzheimer Research (SAO-FRMA), a Methusalem Excellence Grant of the Flemish Government, the Research Foundation Flanders (FWO), the University of Antwerp and the Special Research Fund of the University of Antwerp, the Antwerp Medical Research Foundation and Neurosearch, Belgium. KS and KB are postdoctoral fellows of the FWO. The Antwerp site authors thank the personnel of the VIB Genetic Service Facility, the Biobank of the Institute Born-Bunge and the Departments of Neurology and Memory Clinics at the ZNA Middelheim, ZNA Hoge Beuken and the University Hospitals Leuven. Italian sample collections: the Florence site was supported by a grant from the Italian Ministry of Health: Progetto Strategico 2007- convenzione 39 (sottoprogetto 2). The Milan site was supported by grants from Italian Ministry of Health (PS39) and ‘Fondazione Monzino’. We thank the expert contribution of Mr. Carmelo Romano. Spanish sample collection: the Madrid site (MB) was supported by grants of the Ministerio de Educación y Ciencia and the Ministerio de Sanidad y Consumo (Instituto de Salud Carlos III), and an institutional grant of the Fundación Ramón Areces to the CBMSO. We thank I Sastre and Dr A Martínez-García for the preparation and control of the DNA collection, and Drs P Gil and P Coria for their cooperation in the cases/controls recruitment. We are grateful to the Asociación de Familiares de Alzheimer de Madrid (AFAL) for continuous encouragement and help. The Gran Canaria site was supported by grants of the Ministerio de Sanidad y Consumo (Instituto de Salud Carlos III), and the Fundacion Canaria de Investigacion en Salud (FUNCIS). German sample collection: the Saarland site at the UdS was supported by the German National Genome Research Network (NGFN), grant 01GS08125 to MR, the Helmholtz Alliance for Mental Health in Ageing Societies (HelMA), and the Deutsche Forschungsgemeinschaft (DFG) grant INST 256/317-1 FUGG. Swedish sample collection: the department of medical Epidemiology and Biostatistics at the Karolinska Institute is funded by the US National Institutes of Health (AG028555, AG08724, AG 04563, AG10175, AG08861. The sample collection contributed by LF and CG were financially supported in part by the Swedish Brain Power network, the Marianne and Marcus Wallenberg Foundation, the Swedish Research Council (521-2010-3134), the King Gustaf V and Queen Victoria's Foundation of Freemasons, the Regional Agreement on Medical Training and Clinical Research (ALF) between Stockholm County Council and the Karolinska Institutet. We also wish to acknowledge Anne Kinhult-Ståhlbom, Charlotte Forsell, Håkan Thonberg and Lena Lilius for preparation of the DNA and samples.
Author information
Authors and Affiliations
Consortia
Corresponding author
Ethics declarations
Competing interests
The authors declare no conflict of interest.
Additional information
Supplementary Information accompanies the paper on the Molecular Psychiatry website
Supplementary information
Rights and permissions
This work is licensed under the Creative Commons Attribution-NonCommercial-Share Alike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/
About this article
Cite this article
Lambert, JC., Grenier-Boley, B., Harold, D. et al. Genome-wide haplotype association study identifies the FRMD4A gene as a risk locus for Alzheimer's disease. Mol Psychiatry 18, 461–470 (2013). https://doi.org/10.1038/mp.2012.14
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/mp.2012.14