Skip to main content
  • Research article
  • Open access
  • Published:

Both common variations and rare non-synonymous substitutions and small insertion/deletions in CLU are associated with increased Alzheimer risk

Abstract

Background

We have followed-up on the recent genome-wide association (GWA) of the clusterin gene (CLU) with increased risk for Alzheimer disease (AD), by performing an unbiased resequencing of all CLU coding exons and regulatory regions in an extended Flanders-Belgian cohort of Caucasian AD patients and control individuals (n = 1930). Moreover, we have replicated genetic findings by targeted resequencing in independent Caucasian cohorts of French (n = 2182) and Canadian (n = 573) origin and by performing meta-analysis combining our data with previous genetic CLU screenings.

Results

In the Flanders-Belgian cohort, we identified significant clustering in exons 5-8 of rare genetic variations leading to non-synonymous substitutions and a 9-bp insertion/deletion affecting the CLU β-chain (p = 0.02). Replicating this observation by targeted resequencing of CLU exons 5-8 in 2 independent Caucasian cohorts of French and Canadian origin identified identical as well as novel non-synonymous substitutions and small insertion/deletions. A meta-analysis, combining the datasets of the 3 cohorts with published CLU sequencing data, confirmed that rare coding variations in the CLU β-chain were significantly enriched in AD patients (ORMH = 1.96 [95% CI = 1.18-3.25]; p = 0.009). Single nucleotide polymorphisms (SNPs) association analysis indicated the common AD risk association (GWA SNP rs11136000, p = 0.013) in the 3 combined datasets could not be explained by the presence of the rare coding variations we identified. Further, high-density SNP mapping in the CLU locus mapped the common association signal to a more 5' CLU region.

Conclusions

We identified a new genetic risk association of AD with rare coding CLU variations that is independent of the 5' common association signal identified in the GWA studies. At this stage the role of these coding variations and their likely effect on the β-chain domain and CLU protein functioning remains unclear and requires further studies.

Background

Genome-wide association (GWA) studies lead to long-awaited breakthroughs in the genetics of late-onset Alzheimer disease (AD) [MIM 104300] [1] by providing conclusive genetic association evidence for novel AD risk genes [25]. Notably, GWA significance with similar effect sizes was reached for the top single nucleotide polymorphism (SNP) rs11136000 in the clusterin gene (CLU) [MIM 185430] [1]. The CLU protein (also known as apolipoprotein J) is a multifunctional protein showing functional similarities with the major apolipoprotein of the brain, apolipoprotein E (APOE) [6]. In relation to AD, CLU expression is increased in pyramidal neurons and astrocytes of the hippocampus and entorhinal cortex, the most severely affected brain regions in AD [7]. CLU is present in senile plaques [8], binds to Aβ and is involved in Aβ42 clearance across the blood brain barrier [9]. Moreover, CLU enhances endocytosis of Aβ aggregates to brain phagocytes [10]. Taking its variety of physiological functions, CLU could be a guardian or enemy in AD [11].

The CLU transcriptional unit is located in the chromosomal region 8p21-p12 and comprises 9 exons in the longest transcript that translates in the main CLU protein isoform of a 449 amino-acid residues. The CLU precursor peptide is internally cleaved to produce an α- and β-subunit, held together by disulphide bridges and is subsequently secreted from the cell (Figure 1)[12].

Figure 1
figure 1

Schematic location of rare CLU coding variants identified in stage I, II and III resequencing. (A) Schematic presentation of CLU gene structure, CLU transcript 1 [NM_001831.2] and CLU protein [NP_001822.2]. Coding variants observed in AD patients only are indicated in red, variants observed in patients and controls in blue, variants detected in control individuals only in green. All predicted pathogenic variants are indicated in bold. After cleavage of the signal peptide, the secreted CLU form (449 AA) contains two coiled-coiled domains (pink), three amphipathic domains (blue) and a cysteine rich region (yellow) with 5 disulfide bridges (grey). Six N-glycosylation sites are marked in purple. For ease of interpretation, amino acids are given for specific CLU domains and for detected protein variants only. (B) Conservation alignment of amino acids of CLU beta-chain variants is shown in different species; Homo sapiens (ENSP00000315130), Gorilla gorilla (ENSGGOP00000016521), Pan troglodytes (ENSPTRP00000034423), Pongo abelii (ENSPPYP00000020696), Macaca mulata (ENSMMUP0000003216), Nomascus leucogenys (ENSNLEP00000020015), Tarsius syrichta (ENSTSYP00000001230), Mus musculus (ENSMUSP00000022616), Rattus norvegicus (ENSRNOP00000022095), Canis lupus familiaris (ENSCAFP00000012350) and Bos taurus (ENSBTAP00000007324). Similar to panel A, patient specific variants are marked in red, variants observed in patients and controls in blue and variants in control individuals in green. All predicted pathogenic variants are marked in bold.

In this follow-up study aiming at identifying the genetic variant underlying the CLU association with increased AD risk, we examined the CLU genetic variability using a resequencing approach including all coding exons and regulatory regions in a well-documented Flanders-Belgian patient/control cohort (n = 1930 subjects) [2, 3]. Significant genetic findings were verified by targeted resequencing in two independently ascertained French [3] and Canadian [13] AD cohorts and by meta-analyses of the genetic data sets generated in this study with published data sets obtained of previous genetic screenings of CLU in AD cohorts [14, 15].

Results

CLUresequencing

The stage I resequencing of the coding exons and the regulatory regions of CLU in patients and control individuals of the Flanders-Belgian AD cohort (n = 1930) (Table 1), identified in total 19 rare to intermediate rare non-synonymous single nucleotide variations predicting an amino acid substitution in the CLU protein of which only 5 had been reported earlier [14, 15] (Table 2). Further, we detected an in-frame 9-bp deletion predicting a 3 amino acid deletion p.T445_D447del. Fourteen of the 19 non-synonymous substitutions occurred in 31 AD patients (n = 849, 3.6%), of which 8 appeared only in patients (n = 11), and 6 in patients (n = 20) and control individuals (n = 20). One AD patient carried 2 non-synonymous substitutions p.R338W and p.T345M, and one healthy individual carried both p.S16R and p.R234H. All 3 AD patients with p.T445_D447del carried also p.A309T. The remaining 5 non-synonymous substitutions occurred only in control individuals (n = 5), adding up to 25 control carriers (n = 659, 3.8%) (Table 2 Figure 1).

Table 1 Description of cohorts.
Table 2 Rare non-synonymous CLU variants identified in Stage I Flanders-Belgium AD cohort.

Seven of the 14 non-synonymous substitutions observed in patients were predicted to have a possible or probable harmful effect on CLU protein structure and function using physical and comparative considerations embedded in the PolyPhen software [16](Table 2). SIFT software [17, 18] assessing whether the amino acid substitutions had an effect on CLU protein function based on evolutionary conservation identified 3 non-tolerated substitutions (p.R338W, p.S16R and p.H37Q) (Table 2) (Figure 1). For 4 patients carrying a rare non-synonymous CLU variant, we obtained positive disease history data. In their pedigree, we identified only one other affected first degree relative which suggested that these rare variants have an intermediate disease penetrance.

In contrast, all 5 non-synonymous substitutions detected only in control individuals were labeled by PolyPhen software as benign and were tolerated according to SIFT software (Table 2 Figure 1).

Genetic variations in regulatory elements were also detected in AD patients but none of these genetic variants were located within conserved transcription factor binding sites (See table part D in Additional file 1). In addition, no CLU variants were found affecting splicing using FSPLICE, Netgene2 and SPL prediction programs (See table part E in Additional file 1).

Association analyses of rare CLUvariants

Though overall the percentages of CLU rare variant carriers among AD patients and control individuals was similar, visible inspection of the distribution of the non-synonymous substitutions in AD patients suggested clustering in exons 5-8 which for the most part are coding for the CLU β-chain domain (Table 2) (Figure 1). Statistical analysis confirmed significant clustering in exons 5-8 with 3.8 times more carriers of non-synonymous substitutions or a deletion mutation among patients than control individuals (p = 0.02). To confirm this observation, we used targeted resequencing of exons 5-8 in two independent replication AD cohorts of French and Canadian origin (n = 2755). Again we identified non-synonymous CLU substitutions of which p.R338W, p.N369H and insertion/deletion p.T445_D447del that were already observed in Flanders-Belgian patients. Also, in the Lille AD cohort one 1 bp-insertion predicted a frameshift and premature stop codon - p.I303NfsX13 was found in one AD patient (Table 3) (Figure 1).

Table 3 Rare non-synonymous CLU variants identified in Stage II Lille and Toronto AD cohorts.

Allele sharing among p.T445_D447del carriers i.e. 3 Flanders-Belgian, 3 French, 1 Canadian patients and offspring of two Belgian patients was examined by genotyping short tandem repeat markers located in a 1.3 Mb region of the CLU locus (See table in Additional file 2). DNA available of two children (non-affected, inclusion ages < 50 years) allowed to reconstruct phased haplotypes. No other family members were available. The 3 Flanders-Belgian patients shared alleles at 10 neighboring markers constituting a haplotype of 1.3 Mb around CLU. Seven of these consecutive markers were shared by 2 French patients. A third French patient and a Canadian patient carrier shared 3 markers of this haplotype near CLU. Of note, only Flanders-Belgian patients also carried non-synonymous variant p.A309T.

All non-synonymous substitutions and the deletion/insertion mutations observed in the 3 AD cohorts affecting the CLU β-chain domain were carried forward in a meta-analysis (Table 4). Similarly, we included in the meta-analysis previously reported non-synonymous substitutions that were observed in Portuguese, UK and US Caucasian AD cohorts [14, 15]. Combining the information of in total 6724 patient and 4820 control chromosomes, the meta-analysis confirmed that non-synonymous substitutions and insertion/deletions in exons 5-8 were significantly more enriched in AD patients compared to healthy control individuals (ORMH 1.96 [95% CI 1.18-3.25]; p MH = 0.009, pwoolf = 0.551) (Table 4).

Table 4 Meta-analysis of rare CLU genetic variants.

To assess whether the mutation frequencies of rare coding variants differed by ethnicity, we also resequenced exons 5-8 in a Caribbean Hispanic AD patient/control cohort (Table 5). This identified 2 novel benign non-synonymous substitutions affecting the CLU β-chain domain - p.E431Q and p.D331N. In addition, 3 non-synonymous substitutions found in Belgian, French and Canadian individuals were more frequently observed in Caribbean Hispanics - p.N369H, p.D380N and p.S448L (Table 5) (Figure 1).

Table 5 Rare non-synonymous CLU variants identified in Stage III Caribbean Hispanic AD cohort.

Common CLUassociation analysis

We analyzed association in the Flanders-Belgian AD cohort with the GWA top SNP rs11136000 and 14 other SNPs in the CLU locus, aiming at replicating the common AD risk association with CLU as well as fine-mapping the common association signal. We observed significant allelic association with rs11136000 (p = 0.004) and 4 other SNPs (p < 0.05) (Table 6a). Stratification by APOE ε4 genotype indicated that allelic associations were only significant in the APOE ε4 + stratum (See table in Additional file 3). All five SNPs are positioned within the same LD block (D' between 0.12 and 1.00, r2 between 0.00 and 0.96) (See figure in Additional file 4 and table in Additional file 5). Conditional logistic regression analysis showed that the 4 associated SNPs were not independent from GWA top SNP rs11136000 (data not shown).

Table 6 Common CLU associations with AD.

Next, we genotyped the 5 associated SNPs in the Lille and Toronto AD cohorts and calculated allelic association with AD risk (Table 6a). Nominal significance was observed in the Lille AD cohort with 4 out of 5 SNPs and in the APOE ε4+ stratum of the Lille AD cohort with rs9331908 (p = 0.024) (See table in Additional file 6). In the Toronto cohort association was only significant in the APOE ε4+ stratum with rs867230 (p = 0.024) and rs7982 (p = 0.023) (See table in Additional file 6).

Meta-analysis combining data from all 3 AD cohorts, confirmed allelic association with GWA SNP rs11136000 (p = 0.013), but showed the strongest association after Bonferroni correction with rs1532278 (p = 0.005) (Table 6b, figure in Additional file 7).

Next, we investigated whether the significantly increased presence of rare non-synonymous substitutions and deletion/insertions in AD patients could have driven the common association of CLU with AD risk. Hereto, we recalculated the allelic association in the meta-analysis after excluding the 58 carriers from the Flanders-Belgian (26), the Lille (24) and the Toronto (8) AD cohorts. Allelic associations with the 5 common SNPs remained significant (p-values < 0.05), indicating that the associations with the common SNPs and the rare coding variations represent two independent observations (data not shown).

Discussion

We used an extensive resequencing approach to follow-up on the significant association of common SNPs in the CLU locus with increased risk for AD in 2 independent GWA studies. In the discovery stage, we used unbiased resequencing of all coding and regulatory regions of CLU in a Flanders-Belgian AD cohort followed by targeted sequencing of 2 independent replication cohorts of French and Canadian origin. In the Flanders-Belgian AD cohort, we obtained significant evidence for clustering of rare coding variants in exons 5 to 8, predicting 3 possible (p.T345M, p.N369H, p.T440M) and one probable (p.R338W) damaging non-synonymous substitutions and a 9-bp insertion/deletion p.T445_D447del affecting the CLU β-chain domain. Interestingly, 4 AD patients carried multiple CLU variations: a patient with a positive family history and onset age of 68 years harbored 2 predicted pathogenic substitutions (p.R338W, p.T345M) and three patients carried the deletion p.T445_D447del together with the benign variant p.A309T. Targeted resequencing in the Lille and Toronto AD cohorts identified predicted pathogenic substitutions (p.R338W and p.N369H) in French AD patients and the p.T445_D447del deletion in French and Canadian AD patients. Noteworthy, we observed haplotype sharing with 7 markers around the insertion/deletion site in 3 Belgian and 2 French AD patient carriers which is suggestive for a common ancestor. Partial haplotype sharing was further observed with another French and Canadian patient.

Moreover, we identified a 1-bp insertion/deletion variation p.I303NfsX13, predicting a premature stop codon after 315 amino acid residues (24 of which in β-chain), in a French patient (onset age 68 years). This frame-shift mutation would most likely produce an unstable transcript that is degraded by the nonsense mediated mRNA decay control system or an unstable C-truncated protein that is degraded within the cell, however, no patient material was available for further testing. Substitutions p.R338W and p.T345M are positioned inside the disulphide rich region which is involved in the folding of the α and β-domains into a heterodimer complex. Variant p.N369H deletes a CLU N-glycosylation signal [19], and since CLU is highly glycosylated (20-25% of its total mass comprises carbohydrates), this substitution might affect CLU processing or functioning in ligand interaction. Although absent from all Belgian and French control individuals (> 2500 control chromosomes), p.N369H was present in one Canadian control person aged 58 years, which is suggestive of an intermediate penetrance for this substitution. We also observed 3 non-synonymous substitutions with possible damaging effects (α-chain and β-chain) in both Flanders-Belgian patients and control individuals. As these variants were found in 1 to 4 control individuals only and the control group had on average a younger inclusion age than the patient group, these substitutions might represent risk factors of intermediate disease penetrance and might still contribute to AD risk. In contrast, all non-synonymous substitutions observed in control individuals were predicted to be benign and to unlikely affect CLU protein functioning, strengthening our observation that rare non-synonymous substitutions might contribute to AD risk.

Meta-analysis including all rare coding variants affecting the CLU β-chain domain observed in the discovery and replication cohorts plus those reported in 3 additional cohorts [14, 15], totaling 5772 patients and controls, lent further support to an increased occurrence of rare coding variants predicted to affect the CLU β-chain in AD patients compared to healthy individuals (ORMH 1.96 [95% CI 1.18-3.25]; p = 0.009). Assessing the prevalence of variations affecting the β-chain domain by targeted resequencing of CLU exons 5 to 8 in a Caribbean Hispanic AD cohort, detected novel benign variants while predicted pathogenic variants found in stage I and II AD cohorts were absent. Of note, predicted pathogenic variant p.N369H was more frequently observed in Hispanic patient and control individuals then in stage I and II Caucasian cohorts, suggesting this variant might be a polymorphism which is in agreement with previous reports in African-Americans [15]. Excluding p.N369H from our meta-analysis on rare variants did not alter the significance of our observation that rare coding variants affecting the CLU β-chain are more frequently observed in AD patients than control individuals (p = 0.001). It is possible that rare CLU variants contribute in a different way to disease risk in diverse populations and alternatively, that part of the genetic CLU variability is population-specific (e.g. Hispanic origin). Screening other and larger AD cohorts of different ethnicity and conducting functional analyses on putative pathogenic variants will help unraveling the role of these rare coding variants to AD risk. In this context it is relevant to mention that we have primarily used protein prediction algorithms that have their limitations as well. While our meta-analysis suggested that the CLU β-chain is involved in AD risk, to adequately discriminate between rare or low frequent pathogenic variants and benign polymorphisms, their pathogenicity needs to be experimentally evaluated. These experiments will have to be performed along with comparable efforts aiming at understanding the role of CLU in AD pathogenesis.

Although we identified an association of AD risk with clustering of rare coding variants affecting the CLU β chain, this association was independent from the association with rs11136000, the top GWA SNP [2, 3]. Genotyping rs11136000 in the stage I and II cohorts followed by a meta- analysis showed significant association (p = 0.013). This association remained significant after exclusion of carriers with predicted pathogenic β-chain variants. This is in line with recent LRRK2 observations in which independent common risk associations were found after exclusion of carriers with known pathogenic mutations [20]. Further, high-density association fine-mapping by genotyping common SNPs throughout the CLU locus in the stage I cohort, identified 4 additional SNP associations, that all except one (rs9331908) were in high LD with rs11136000 (r2 > 0.89) and located in the 5' CLU region. This association region ranges from CLU intron 1 to exon 5, which is only marginally overlapping with the rare variant association region in the β-chain (exon 5 to 8). Combined with the two independent association signals for common and rare CLU variants, this could imply two distinct mechanisms underlying AD risk.

Further genotyping of the 4 significant SNPs and combined meta-analyses on all 3 AD cohorts showed, after Bonferroni correction, the strongest evidence of association with rs1532278 (p = 0.005). This variant is located in an intron 3 sequence that strongly resembles a regulatory element based on sequence alignment of seven species. Also, association evidence was found for rs867230 in CLU intron 1, which potentially affects a MEF2 transcription factor binding site. These two novel associated common variants may represent variants of functional relevance underlying the common association signal with rs11136000 in various populations [2123]. In agreement with other reports [14, 21], we did not observe any association with common coding variants. Though we did not find association with the recently described splice site variant rs9331888 [24] in our Flanders-Belgian population (OR 1.11 [95% CI 0.94-1.31]; p = 0.2), our association findings do point to common variants with possible regulatory effects (rs1532278, rs867230). It is likely that a number of common and rare CLU variants contribute to an increased risk for AD. In addition, different mechanisms may contribute to disease: our common CLU associations might regulate CLU transcription while our rare CLU association signal results from non-synonymous variants affecting protein functioning. Whether these rare CLU variants represent loss or gain of function variants remains at present unclear and requires further studies.

Conclusions

In conclusion, in-depth CLU resequencing showed significant clustering and association of rare coding variants with AD risk in a combined meta-analysis of 3 Caucasian AD cohorts analyzed in this study and cohorts of previous published studies. The rare coding variants are non-synonymous substitutions and insertion/deletion mutations that affect the CLU β-chain domain. While the physiological properties of the β-chain domain remain unclear, our data suggests that this protein subunit may be of interest in AD pathogenesis, and merits follow-up with detailed functional analyses. Also, in our combined AD data set, the rare coding variant association was independent of the common AD risk association signal that we fine-mapped to a more 5' region of CLU. The strongest association was not obtained with GWA SNP rs11136000 but with an intron 3 SNP (rs1532278) located in a sequence with high regulatory potential. Altogether, our data suggest that CLU may be a risk gene in which multiple rare and common variants have independent effects on AD disease susceptibility.

Methods

AD cohorts

Flanders-Belgian AD cohort

In the stage I analysis, we used as discovery sample an extended Flanders-Belgian AD cohort of 1930 subjects. A fraction of this AD cohort (n = 1750, 91%) contributed to the replication phase of the 2 AD GWA studies confirming AD risk association at the CLU locus [2, 3]. Since then we extended the AD cohort dataset with newly ascertained subjects and additional genotypes. The AD patients have been ascertained at memory clinics and neurology divisions in a prospective study of neurodegenerative and vascular dementia in Flanders, the Dutch-speaking region of Belgium [25, 26], and in a comparable prospective study on molecular genetics of cognitive impairment [27] using the same clinical assessments and biosampling schemes (Table 1). Each AD patient underwent a neuropsychological examination, including Mini-Mental State Examination (MMSE)[28] and structural and/or functional neuroimaging [27]. Consensus diagnosis of possible or probable AD was obtained by minimal two neurologists based on the National Institute of Neurological and Communication Disorders and Stroke-Alzheimer's Disease and Related Disorders Association (NINCDS-ADRDA) criteria [29]. Unrelated control individuals had no neurological or psychiatric antecedents or had neurological complaints or organic disease involving the central nervous system. Additional community control individuals were included after interview concerning medical and family history and a Mini Mental State Examination (MMSE) > 24 [28].

French and Canadian AD cohorts

In the stage II analysis, we used two independent Caucasian AD patient/control cohorts for replication of rare coding variants in CLU exons 5-8 by targeted resequencing. A first French AD cohort consisted of patients ascertained in the North of France (Lille AD cohort) that were diagnosed with probable AD according to DSM-III-R and NINCDS/ADRDA criteria. Healthy control individuals were from the same geographical area and were cognitively intact, had no family history of AD or DSM-III-R criteria for dementia and a MMSE ≥ 25 (Table 1) [30, 31]. In a second Canadian AD cohort (Toronto AD cohort [13]), diagnoses of probable or possible AD were defined according to the NINCDS-ADRDA criteria at clinics specialized in memory disorders [29]. Subjects were classified as control individuals when they were without cognitive impairment or dementia at last visit (Table 1) [32].

Caribbean Hispanic AD cohort

In the stage III analysis, we initiated an assessment of the contribution of rare coding CLU variants in populations of different ethnicity by analyzing an AD cohort of Caribbean Hispanics ancestry for rare coding mutations by resequencing CLU exons 5-8 (Table 1). The diagnosis of AD was based on the NINCDS/ADRDA criteria [29, 33].

Ethical assurances

All clinical and genetic studies described in this manuscript were approved by the medical ethical committees of the respective hospital divisions and university genetic laboratories at the respective cohort sampling sites in Flanders-Belgium, France, Canada and USA. Informed consent was obtained from all participants using procedures approved by institutional review boards at each of the clinical research centers enrolling subjects.

CLUresequencing

Unbiased resequencing was performed in the Flanders-Belgian AD cohort by PCR sequencing of all 9 CLU exons and their flanking intron-exon boundaries [NM_001831.2], the 5' and 3'UTR regions of longest transcript 1 [NM_001831.2], the 5'UTR region of transcript 2 [NM_20339.1] and regulatory elements (Table 2 and see table in Additional file 1). Primers were designed using ExonPrimer software and primer3. In stages II and III, we performed targeted resequencing of CLU exons 5-8 in the Lille, Toronto and Caribbean Hispanics AD cohorts (Table 3 Table 5). All sequences were analyzed by two independent researchers using Seqman, the NovoSNP software package [34] or Vector NTI software (Invitrogen).

Statistical genetic analyses

Coding variant association

In stage I, an exon-by-exon approach was used to examine significance of putative clustering in AD patients of rare CLU non-synonymous variants with a MAF of < 0.01-0.05 (Table 4). Alleles were collapsed and overall frequency differences between AD patients and control individuals were compared using χ2 statistics. Additionally, we performed a literature search identifying in PubMed all published studies that were reporting a genetic screening of CLU in AD patients and healthy individuals of Caucasian origin, yielding two publications describing data in three AD cohorts from UK, Portugal and US [14, 15]. For the meta-analysis, we included all rare coding variants predicting non-synonymous and insertion/deletion mutations in the CLU β-chain domain. Rare variant counts in patients and control individuals were obtained for each dataset. To obtain summary ORs and 95% CI, a fixed-effect (Mantel-Haenszel) meta-analysis was performed in R using the library rmeta-version 2.16.

Common SNP associations

The stage I Flanders-Belgian AD cohort (n = 1930, Table 1) had an estimated power of 97% to detect a modest association with an odds ratio (OR) of ~1.5 and a common genetic variation with a minor allele frequency (MAF) of 0.2 [35]. Tagging SNPs were selected throughout the CLU locus (chr8:27.538.244 - 27.500.368; including 10 kb up- and downstream regions), using the CEPH (Centre d'Etude du Polymorphism Humain) population of HapMap (Data Rel24/phase II Nov08). SNPs with a MAF < 0.05, with a Hardy-Weinberg equilibrium (HWE) p-value < 0.05 or within repeat sequences (> 50%) were removed from the selection. A total of 15 SNPs including the GWA top SNP rs11136000 [2, 3], were used in the stage I Flanders-Belgian AD cohort for association analysis and for fine-mapping of the association signal in the CLU locus (See figure in Additional file 4).

Four SNPs were genotyped by sequencing (rs7982, rs867230, rs3216167, rs11136000), the others SNPs by Sequenom MassArray® assay, followed by MALDI-TOF mass spectrometry (Sequenom, Inc., Hamburg, Germany). PCR and extension primers were designed using Assay Design 3.1 Software. Genotypes were scored both automatically (MassArray Typer version 4.0) and by two researchers blind to disease status. Genotyping success rate for SNPs genotyped with Sequenom was > 95%. Inter-plate controls showed 100% concordance for all genotyped SNPs. Deviations from HWE were determined using an exact HWE test [36].

For common SNPs, differences in allele frequencies between AD patients and control individuals were tested using χ2 statistics. P-values and odds ratios (OR) with 95% confidence intervals for the minor allele were calculated relative to the common allele and corrected for gender, onset age or age at inclusion using binary logistic regression models (Table 6a). Allelic associations were further stratified for APOE status i.e. presence or absence of APOE ε4 alleles (See table in additional file 3 and figure in Additional file 4). Statistical analyses were performed using SPSS 16.0 (SPSS, Inc., Chicago, IL). SNPs showing significant association at stage I, were genotyped in the Lille and Toronto AD cohorts and associations calculated (See table in Additional file 6). Fixed-effect (Mantel-Haenszel) meta-analysis was performed based on effect estimates of the separate cohorts adjusted for age, gender and APOE ε4 status in R using the library rmeta-version 2.16. Woolf's test for heterogeneity was calculated. A Bonferroni p-value < 0.05, corrected for the number of SNPs included in the meta-analysis (n = 5), was considered significant (Table 6b and see forest plots figure in Additional file 7).

Abbreviations

AA:

Amino Acids

AD:

Alzheimer's disease

AAO:

age at onset

AAI:

age at inclusion

CI:

confidence interval

CLU :

clusterin

GWA:

genome-wide association

HWE:

Hardy-Weinberg equilibrium

MAF:

minor allele frequency

nt:

nucleotides

OR:

odds ratio

SD:

standard deviation

SNP:

single nucleotide polymorphism

UTR:

untranslated region.

References

  1. Online Mendelian Inheritance in Man (OMIMI) database. [http://www.ncbi.nlm.nih.gov/Omim]

  2. Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, Hamshere ML, Pahwa JS, Moskvina V, Dowzell K, Williams A, et al: Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat Genet. 2009, 41: 1088-1093. 10.1038/ng.440.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Lambert JC, Heath S, Even G, Campion D, Sleegers K, Hiltunen M, Combarros O, Zelenika D, Bullido MJ, Tavernier B, et al: Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet. 2009, 41: 1094-1099. 10.1038/ng.439.

    Article  CAS  PubMed  Google Scholar 

  4. Naj AC, Jun G, Beecham GW, Wang LS, Vardarajan BN, Buros J, Gallins PJ, Buxbaum JD, Jarvik GP, Crane PK, 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. 10.1038/ng.801.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Seshadri S, Fitzpatrick AL, Ikram MA, DeStefano AL, Gudnason V, Boada M, Bis JC, Smith AV, Carassquillo MM, Lambert JC, et al: Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA. 2010, 303: 1832-1840. 10.1001/jama.2010.574.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Sleegers K, Lambert JC, Bertram L, Cruts M, Amouyel P, Van BC: The pursuit of susceptibility genes for Alzheimer's disease: progress and prospects. Trends Genet. 2010, 26: 84-93. 10.1016/j.tig.2009.12.004.

    Article  CAS  PubMed  Google Scholar 

  7. May PC, Lampert-Etchells M, Johnson SA, Poirier J, Masters JN, Finch CE: Dynamics of gene expression for a hippocampal glycoprotein elevated in Alzheimer's disease and in response to experimental lesions in rat. Neuron. 1990, 5: 831-839. 10.1016/0896-6273(90)90342-D.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  9. Bell RD, Sagare AP, Friedman AE, Bedi GS, Holtzman DM, Deane R, Zlokovic BV: Transport pathways for clearance of human Alzheimer's amyloid beta-peptide and apolipoproteins E and J in the mouse central nervous system. J Cereb Blood Flow Metab. 2007, 27: 909-918.

    PubMed Central  CAS  PubMed  Google Scholar 

  10. Bartl MM, Luckenbach T, Bergner O, Ullrich O, Koch-Brandt C: Multiple receptors mediate apoJ-dependent clearance of cellular debris into nonprofessional phagocytes. Exp Cell Res. 2001, 271: 130-141. 10.1006/excr.2001.5358.

    Article  CAS  PubMed  Google Scholar 

  11. Nuutinen T, Suuronen T, Kauppinen A, Salminen A: Clusterin: A forgotten player in Alzheimer's disease. Brain Res Rev. 2009, 61: 89-104. 10.1016/j.brainresrev.2009.05.007.

    Article  CAS  PubMed  Google Scholar 

  12. Wilson MR, Easterbrook-Smith SB: Clusterin is a secreted mammalian chaperone. Trends Biochem Sci. 2000, 25: 95-98. 10.1016/S0968-0004(99)01534-0.

    Article  CAS  PubMed  Google Scholar 

  13. Rogaeva EA, Premkumar S, Grubber J, Serneels L, Scott WK, Kawarai T, Song Y, Hill DL, Abou-Donia SM, Martin ER, et al: An alpha-2-macroglobulin insertion-deletion polymorphism in Alzheimer disease. Nat Genet. 1999, 22: 19-22. 10.1038/8729.

    Article  CAS  PubMed  Google Scholar 

  14. Guerreiro RJ, Beck J, Gibbs JR, Santana I, Rossor MN, Schott JM, Nalls MA, Ribeiro H, Santiago B, Fox NC, et al: Genetic Variability in CLU and Its Association with Alzheimer's Disease. PLoS One. 2010, 5: e9510-10.1371/journal.pone.0009510.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Tycko B, Feng L, Nguyen L, Francis A, Hays A, Chung WY, Tang MX, Stern Y, Sahota A, Hendrie H, et al: Polymorphisms in the human apolipoprotein-J/clusterin gene: ethnic variation and distribution in Alzheimer's disease. Hum Genet. 1996, 98: 430-436. 10.1007/s004390050234.

    Article  CAS  PubMed  Google Scholar 

  16. Polyphen. [http://genetics.bwh.harvard.edu/pph/]

  17. Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31: 3812-3814. 10.1093/nar/gkg509.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. SIFT. [http://blocks.fhcrc.org/sift/SIFT.html]

  19. Kapron JT, Hilliard GM, Lakins JN, Tenniswood MP, West KA, Carr SA, Crabb JW: Identification and characterization of glycosylation sites in human serum clusterin. Protein Sci. 1997, 6: 2120-2133.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Ross Owen, Soto-Ortolaza Alexandra, Heckman Michael, Aasly Jan, Nadine Abahuni, Grazia Annesi, Bacon Justin, Soraya Bardien, Maria Bozi, Alexis Brice, et al: LRRK2 exonic variants and susceptibility to Parkinson's disease. Lancet Neurology. 2011

    Google Scholar 

  21. Carrasquillo MM, Belbin O, Hunter TA, Ma L, Bisceglio GD, Zou F, Crook JE, Pankratz VS, Dickson DW, Graff-Radford NR, et al: Replication of CLU, CR1, and PICALM Associations With Alzheimer Disease. Arch Neurol. 2010

    Google Scholar 

  22. Corneveaux JJ, Myers AJ, Allen AN, Pruzin JJ, Ramirez M, Engel A, Nalls MA, Chen K, Lee W, Chewning K, et al: Association of CR1, CLU and PICALM with Alzheimer's disease in a cohort of clinically characterized and neuropathologically verified individuals. Hum Mol Genet. 2010

    Google Scholar 

  23. Jun G, Naj AC, Beecham GW, Wang LS, Buros J, Gallins PJ, Buxbaum JD, Ertekin-Taner N, Fallin MD, Friedland R, et al: Meta-analysis Confirms CR1, CLU, and PICALM as Alzheimer Disease Risk Loci and Reveals Interactions With APOE Genotypes. Arch Neurol. 2010

    Google Scholar 

  24. Szymanski M, Wang R, Bassett SS, Avramopoulos D: Alzheimer's risk variants in the Clusterin gene are associated with alternative splicing. Transl Psychiatr. 2011, 1:

    Google Scholar 

  25. Engelborghs S, Dermaut B, Goeman J, Saerens J, Marien P, Pickut BA, van den Broeck M, Serneels S, Cruts M, Van Broeckhoven C, et al: Prospective Belgian study of neurodegenerative and vascular dementia: APOE genotype effects. J Neurol Neurosurg Psychiatry. 2003, 74: 1148-1151. 10.1136/jnnp.74.8.1148.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Engelborghs S, Dermaut B, Marien P, Symons A, Vloeberghs E, Maertens K, Somers N, Goeman J, Rademakers R, van den Broeck M, et al: Dose dependent effect of APOE epsilon4 on behavioral symptoms in frontal lobe dementia. Neurobiol Aging. 2006, 27: 285-292. 10.1016/j.neurobiolaging.2005.02.005.

    Article  CAS  PubMed  Google Scholar 

  27. Bettens K, Brouwers N, Van Miegroet H, Gil A, Engelborghs S, De Deyn PP, Vandenberghe R, Van Broeckhoven C, Sleegers K: Follow-Up Study of Susceptibility Loci for Alzheimer's Disease and Onset Age Identified by Genome-Wide Association. J Alzheimers Dis. 2009, 19: 1169-1175.

    Google Scholar 

  28. Folstein MF, Folstein SE, McHugh PR: "Mini-mental state". A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975, 12: 189-198. 10.1016/0022-3956(75)90026-6.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  Google Scholar 

  30. Chapuis J, Hot D, Hansmannel F, Kerdraon O, Ferreira S, Hubans C, Maurage CA, Huot L, Bensemain F, Laumet G, et al: Transcriptomic and genetic studies identify IL-33 as a candidate gene for Alzheimer's disease. Mol Psychiatry. 2009, 14: 1004-1016. 10.1038/mp.2009.10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Dreses-Werringloer U, Lambert JC, Vingtdeux V, Zhao H, Vais H, Siebert A, Jain A, Koppel J, Rovelet-Lecrux A, Hannequin D, et al: A polymorphism in CALHM1 influences Ca2+ homeostasis, Abeta levels, and Alzheimer's disease risk. Cell. 2008, 133: 1149-1161. 10.1016/j.cell.2008.05.048.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Reitz C, Tokuhiro S, Clark LN, Conrad C, Vonsattel JP, Hazrati LN, Palotas A, Lantigua R, Medrano M, Jimenez-Velazquez Z, et al: SORCS1 alters amyloid precursor protein processing and variants may increase Alzheimer's disease risk. Ann Neurol. 2011, 69: 47-64. 10.1002/ana.22308.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Lee JH, Cheng R, Barral S, Reitz C, Medrano M, Lantigua R, Jimenez-Velazquez IZ, Rogaeva E, St George-Hyslop PH, Mayeux R: Identification of novel loci for Alzheimer disease and replication of CLU, PICALM, and BIN1 in Caribbean Hispanic individuals. Arch Neurol. 2011, 68: 320-328. 10.1001/archneurol.2010.292.

    PubMed Central  PubMed  Google Scholar 

  34. Weckx S, Del Favero J, Rademakers R, Claes L, Cruts M, De Jonghe P, Van Broeckhoven C, De Rijk P: novoSNP, a novel computational tool for sequence variation discovery. Genome Res. 2005, 15: 436-442. 10.1101/gr.2754005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Purcell S, Cherny SS, Sham PC: Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19: 149-150. 10.1093/bioinformatics/19.1.149.

    Article  CAS  PubMed  Google Scholar 

  36. Exact HWE. [http://www.pharmgat.org/IIPGA2/Bioinformatics/exacthweform]

  37. Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36: D149-D153.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.

    Article  CAS  PubMed  Google Scholar 

  39. Grun D, Wang YL, Langenberger D, Gunsalus KC, Rajewsky N: microRNA target predictions across seven Drosophila species and comparison to mammalian targets. PLoS Comput Biol. 2005, 1: e13-10.1371/journal.pcbi.0010013.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Georges M, Coppieters W, Charlier C: Polymorphic miRNA-mediated gene regulation: contribution to phenotypic variation and disease. Curr Opin Genet Dev. 2007, 17: 166-176. 10.1016/j.gde.2007.04.005.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors acknowledge the important participation of patients and their relatives as well as control individuals. The authors are grateful to the personnel of the VIB Genetic Service Facility, the Biobank of the Institute Born-Bunge, the Departments of Neurology and Memory Clinics Hospital Network Antwerp Middelheim and Hoge Beuken and the University Hospitals Leuven. The research at the Antwerp site was in part supported by the Interuniversity Attraction Poles program P6/43 of the Belgian Science Policy Office (BELSPO, http://www.belspo.be/), the Foundation for Alzheimer Research (SAO/FRMA, http://alzh.org/), a Methusalem Excellence Grant of the Flemish Government (EWI, http://www.ewi-vlaanderen.be/), the Research Foundation Flanders (FWO, http://www.fwo.be/), the Special Research Fund of the University of Antwerp (UA, http://www.ua.ac.be/), the Antwerp Medical Research Foundation and Neurosearch, Belgium. KB, NB and KS are postdoctoral fellows, and RV a senior clinical investigator of the FWO. The Lille site was funded in part by the National Foundation for Alzheimer's disease and related disorders, the Institute Pasteur de Lille and INSERM. The Toronto site was supported by grants from the Canadian Institutes of Health Research, Ontario Research Fund, Weston foundation (PSH), the Howard Hughes Medical Institute, The Wellcome Trust, the Alzheimer Society of Ontario (PSH).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Kristel Sleegers or Christine Van Broeckhoven.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KB, KS, NB and CVB were involved in the conception and design of the study, and in the analysis and interpretation of the data. SE, RV, NLB, MM, KP, PDD, J-CL, ER, RM, PSH collected the patient materials included in this study. KB, SV and JVD participated in the sequencing studies, KB performed the statistical analyses, KB and KS drafted the manuscript. J-CL, ER, RV, FP, RM, PSH, PA, SE, RV, NLB, PDD have helped to revise the manuscript for intellectual content. CVB is the principal investigator. All authors have read and approved the final manuscript.

Electronic supplementary material

13024_2011_322_MOESM1_ESM.DOC

Additional file 1: Synonymous CLU variants and CLU variants in 3' UTR, 5' UTR, regulatory regions, and splice sites in Flanders-Belgian AD cohort. aGene location position according to the longest CLU transcript with 9 coding exons [NM_001831.2], bNumbering according to build GRCh37/hg19 (Feb.2009), nucleotide changes given for complementary negative strand; MAF = Minor Allele Frequency, calculated upon the minimum number of successful sequences (1698 AD alleles, 1318 control alleles) (B) The miRanda algorithm (microRNA resource) was used to predict whether 3' UTR variants are involved in miRNA-binding [37]. TargetScan [38] and Pictar [39] were applied for prediction of miRNA binding site positions. Polymorphisms in microRNA Target Site (PolymiRTS) and Patrocles databases were used [40] to estimate miRNA binding effects of variants. (D) Regulatory elements (proximal promoter, OREG0018671) predicted using FirstEF and OregAnno from UCSC browser Human Mar. 2006) were sequenced, encompassing 458 nt before the translation initiation codon (transcript 1) and 457 nt downstream of coding exon 1 [NM_001831.2]. No variants were positions in conserved transcription factor binding sites (according to UCSC). (E) Variants near splice sites. FSPLICE, Netgene2 and SPL were used to predict possible splicing effects. (DOC 176 KB)

13024_2011_322_MOESM2_ESM.DOC

Additional file 2: Haplotype sharing of p.T445_D447del carriers. Allele sharing of p.T445_D447del carriers (3 Belgian patients, 3 French and 1 Canadian AD patient) was examined by genotyping 10 short tandem repeat (STR) markers located in a 1.3 Mb region covering the CLU locus. aPhysical location of STR markers is relative to the NCBI genome build 36, DNA of two unaffected children (with inclusion age < 50 years) of Belgian Flanders-Belgian carriers allowed reconstructing haplotypes; Flanders-Belgian carriers of the insertion/deletion shared a 10-marker haplotype covering 1.3 Mb around CLU (shared alleles are indicated in bold), two French patient carriers shared alleles at 7 consecutive markers of this haplotype, while French patient 2 and the Canadian patient 1 only shared 3 consecutive markers with the Belgian patients. Of note, Belgian insertion/deletion carriers also carried p.A309T (Table 2), which was not detected in the French or Canadian individuals. (DOC 84 KB)

13024_2011_322_MOESM3_ESM.DOC

Additional file 3: Common CLU allelic associations in Flanders-Belgian APOE ε4 strata. aGene location position according to the longest CLU transcript with 9 coding exons [NM_001831.2], bNumbering according to build GRCh37/hg19 (Feb.2009), minor alleles given for complementary negative strand. Allele frequencies are shown with absolute numbers in brackets. Calculations of odds ratios, presented with 95% confidence intervals (CI), were performed using the common allele as reference allele. Nominal p-values were adjusted for age (onset age for patients, inclusion age for control individuals) and gender in the APOE subgroups. Nominally significant p-values are marked in bold. (DOC 84 KB)

13024_2011_322_MOESM4_ESM.JPEG

Additional file 4: Schematic overview of CLU allelic associations and LD pattern in Flanders-Belgian AD cohort. -log10 p-values for allelic associations of 15 common SNPs (MAF > 0.05) encompassing the CLU locus are given for the total cohort adjusted for age, gender and APOE, and for APOE ε4 genotype strata adjusted for age and gender (Additional file 3). Based upon 15 common SNP genotypes of stage I, the overall linkage disequilibrium (LD) plot was reconstructed with D' as LD measure drawn using the LDheatmap v0.2-8 package. The LD pattern consisted of a major LD block (12 consecutive SNPs starting from intron 3 to 3' intergenic region) and a minor LD block of 3 SNPs upstream from CLU. (JPEG 4 MB)

13024_2011_322_MOESM5_ESM.DOC

Additional file 5: Linkage disequilibrium measures in stage I and II AD cohorts. Pairwise linkage measures (D' and r2) are given for the associated SNPs in stage I and stage II cohorts. (DOC 84 KB)

13024_2011_322_MOESM6_ESM.DOC

Additional file 6: Common CLU allelic associations in replication AD cohorts APOE ε4 strata. Allele frequencies are shown with absolute numbers in brackets, minor alleles given for complementary negative strand. Calculations of odds ratios, presented with 95% confidence intervals (CI), were performed using the common allele as reference allele. Nominal p-values were adjusted for age (onset age for patients, inclusion age for control individuals) and gender in the APOE subgroups. Nominally significant p-values are marked in bold. (DOC 73 KB)

13024_2011_322_MOESM7_ESM.JPEG

Additional file 7: Forest plots of common CLU association in stage I and II AD cohorts. Odds ratio's and 95% confidence intervals are given for stage I (Flanders-Belgian) and stage II cohorts (Lille, Toronto) separately as well as overall combining stage I and II. (JPEG 2 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bettens, K., Brouwers, N., Engelborghs, S. et al. Both common variations and rare non-synonymous substitutions and small insertion/deletions in CLU are associated with increased Alzheimer risk. Mol Neurodegeneration 7, 3 (2012). https://doi.org/10.1186/1750-1326-7-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1750-1326-7-3

Keywords