Research
Research
Research
22 4519–4529
doi: 10.1093/hmg/ddx329
Advance Access Publication Date: 26 August 2017
Association Studies Article
Abstract
Genetic variants at PTK2B–CLU locus pose as high-risk factors for many age-related disorders. However, the role of these variants
in disease progression is less characterized. In this study, we aimed to investigate the functional significance of a clusterin
intronic SNP, rs2279590, that has been associated with pseudoexfoliation, Alzheimer’s disease (AD) and diabetes. We have previ-
ously shown that the alleles at rs2279590 differentially regulate clusterin (CLU) gene expression in lens capsule tissues. This
polymorphism resides in an active regulatory region marked by H3K27Ac and DNase I hypersensitive site and is an eQTL for CLU
expression. Here, we report the presence of an enhancer element in surrounding region of rs2279590. Deletion of a 115 bp
intronic region flanking the rs2279590 variant through CRISPR-Cas9 genome editing in HEK293 cells demonstrated a decreased
clusterin gene expression. Electrophoretic mobility shift and chromatin immunoprecipitation assays show that rs2279590 with
allele ‘A’ constitutes a transcription factor binding site for heat shock factor-1 (HSF1) but not with allele ‘G’. By binding to allele
‘A’, HSF1 abrogates the enhancer effect of the locus as validated by reporter assays. Interestingly, rs2279590 locus has a wide-
spread enhancer effect on two nearby genes, protein tyrosine kinase 2 beta (PTK2B) and epoxide hydrolase-2 (EPHX2); both of
which have been previously associated with AD as risk factors. To summarize, our study unveils a mechanistic role of the
common variant rs2279590 that can affect a variety of aging disorders by regulating the expression of a specific set of genes.
Introduction
The former two conditions, PEX (OMIM: 177650) and AD (OMIM:
Clusterin or Apolipoprotein-J (8p21.1) is a multifunctional pro- 104300), share similar pathological alterations like characteristic
tein with diverse roles. Stabilization of cell–cell and cell–matrix deposition of fibrillar protein aggregates and gradual deteriora-
interactions, lipid transport, apoptosis, inhibition of comple- tion of optic and brain nerves, respectively. Studies have shown
ment activation and clearance of extracellular deposits are that Clusterin accumulation differentiates the severe advanced
some of the crucial functions attributed to clusterin (CLU) (1,2). stage of PEX called pseudoexfoliation glaucoma (PEXG) from the
CLU has been associated at both genomic and proteomic levels less severe syndrome stage, pseudoexfoliation syndrome (PEXS)
with age-related disorders, such as, pseudoexfoliation (PEX), (3). Additionally, brain function in AD-affected individuals with
Alzheimer’s disease (AD) and Type 2 Diabetes mellitus (3–7). over accumulated Clusterin has been reported to deteriorate
Received: May 16, 2017. Revised: August 7, 2017. Accepted: August 22, 2017
C The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com
V
4519
4520 | Human Molecular Genetics, 2017, Vol. 26, No. 22
faster than normal individuals (8–10). It is however unclear, neighbouring AD associated candidate genes, protein tyrosine
whether CLU accumulation in such disorders is the cause or a kinase 2 beta (PTK2B) and epoxide hydrolase-2 (EPHX2); suggest-
consequence. Although studies suggest that CLU plays a crucial ing a distal enhancer effect contributed by this locus harbouring
role in preventing extracellular deposits; it also has been shown rs2279590. In conclusion, our study suggests that a common
that knockout of clusterin in PDAPP mice (AD mice model) risk variant at the CLU locus has a widespread enhancer effect
have reduced fibrillar plaque and neuritic dystrophy (11). in regulating multiple candidate genes associated with PEX/AD
Furthermore, an alternative truncated form of Clusterin called and assists in the pathogenesis of such aging disorders.
nuclear Clusterin (nCLU), is known to induce apoptosis by act-
ing as a pro-death protein (12). Although, nCLU level is in low
abundance in a healthy cell, stress induced increase in its Results
expression can lead to apoptosis of affected neurons; a patho-
Locus containing rs2279590 acts as an enhancer and
genic feature for both PEX and AD.
regulates clusterin gene expression
At the genomic level, common variants within CLU locus also
have been associated with PEX and AD (3,4,6). A variant, Using the data from ENCODE (Encyclopedia of DNA Elements)
rs2279590 within the clusterin gene is of particular importance as project in UCSC genome browser, we observed that the region
it has been picked up as a genetic risk factor for PEX, AD and dia- surrounding the polymorphic site, rs2279590 located in 7th
betes (3,4,6). Earlier, we have shown that individuals homozygous intron of clusterin gene, is within a DNase I hypersensitive site
for the risk allele ‘G’ at rs2279590 have a 2-fold higher clusterin and has modified H3K27Ac mark, suggesting the presence of
expression in lens capsule tissues than those with allele ‘A’ (3); an active regulatory region (Fig. 1). In addition, eQTL data from
suggesting a regulatory role of rs2279590 over CLU expression. GTEx project also indicate a regulatory role of rs2279590 over
Association of rs2279590 with such diverse disorders and its loca- CLU expression in both muscle-skeletal and skin-sun exposed
tion within a high-risk locus demands further investigation into tissue samples examined (Supplementary Material, Fig. S1).
its functional implications in disease progression. We, therefore, cloned a 201 bp region surrounding rs2279590
The current study focuses on finding the functional role, if into a luciferase reporter vector to check its regulatory effect
any, of rs2279590 in clusterin gene regulation. We observed that in HEK293 cells. As shown in Figure 2A, the luciferase activity
the genomic region surrounding rs2279590 is located in an was significantly higher (2-fold, P ¼ 0.03) in constructs con-
active locus marked by H3K27Ac and DNase I hypersensitivity taining the rs2279590 locus as compared with the cells with
site (ENCODE project) and is a quantitative trait loci (eQTL) for empty luciferase vector; implying a regulatory effect of this
CLU expression as evident through GTEx project. Luciferase locus.
reporter assays validated the rs2279590 region as a regulatory In addition, to check the in vivo effect of the locus on clus-
element. Deletion of a 115 bp genomic region flanking rs2279590 terin gene expression, we generated HEK293115/ cells in
polymorphic site by CRISPR/Cas9 system, leads to downregula- which 115 bp genomic region around rs2279590 was deleted
tion of clusterin gene expression in HEK293 cells; suggesting using CRISPR/Cas9 system (Fig. 2B and C). Two independent
that the variant resides within an enhancer element. Transition single-cell derived homozygous clones were subsequently
from risk allele ‘G’ to ‘A’ at rs2279590 was found to create a tran- used to check for CLU gene expression compared with that
scription factor binding site (TFBS) for heat shock factor 1 (HSF1) of control non-deleted cells. Quantitative real-time PCR (qRT-
that abrogated the enhancer activity of the region. The locus PCR) assays (Fig. 2D) confirmed a significant downregulation
surrounding this variant was also found to regulate two (P ¼ 0.01) of transcripts for secretory form of clusterin (sCLU)
Human Molecular Genetics, 2017, Vol. 26, No. 22 | 4521
was observed that a specific protein complex bound to the do not show a prominent shift as seen with allele ‘A’
labelled oligos (Lane 2, Fig. 3A) marked by a shift (arrowhead). (Supplementary Material, Fig. S2A).
With addition of unlabelled oligos comprising allele ‘A’ (Lane 3) A three-tier experimental validation was carried out to vali-
the shift disappears but not with unlabelled oligos with allele ‘G’ date the binding of HSFs at rs2279590. Figure 3B shows a compa-
(Lane 4); suggesting the binding complex is specific for allele ‘A’. rative shift in lanes 2 and 5 when a labelled consensus heat
Similar binding assays with labelled oligos comprising ‘G’ allele shock element (HSE) and labelled rs2279590 oligo with allele ‘A’
Human Molecular Genetics, 2017, Vol. 26, No. 22 | 4523
were used, respectively. A competitive EMSA on labelled oligo performed in HEK293 cells. As shown in Figure 4A and B the
containing rs2279590/A allele, when challenged with increased genomic region surrounding rs2279590 was enriched by qRT-PCR
concentrations (1- and 100-fold excess) of consensus unlabelled after immunoprecipitation with anti-HSF1 antibody but not by
HSE oligo (lanes 3 and 4, Fig. 3C) or unlabelled rs2279590/A oligo rabbit IgG (negative control). This confirms in vivo binding of HSF1
(1-, 10-, 50-, 100-, 200- and 400-fold excess) (lanes 5–10, Fig. 3C) to the sequences comprising of rs2279590.
the shift disappears; indicating that the 29 bp rs2279590/A
sequence binds to protein complexes similar to those binding to
a heat shock element. As shown in Figure 3C, 1-fold addition of Binding of HSF1 to allele ‘A’ at rs2279590 abrogates the
either unlabelled HSE or rs2279590/A (lane 3 and 5, respectively) enhancer effect of the locus
oligos was sufficient to compete for the binding complex, thereby In order to check the allele specific effect on the reporter activity, a
drastically reducing the intensity of the shift. 201 bp intronic region harbouring rs2279590 (with either ‘A’ or ‘G’
Later, when the binding complex was challenged with HSF1- allele) was cloned into pGL4.23 luciferase vector and transfected
specific antibody, disappearance of the shift was noted (lane 4, Fig. into HEK293 cells and the reporter activity was checked. As repre-
3D) indicating that the protein binding complex bound to allele ‘A’ sented in Figure 4C, changing the allele from ‘G’ to ‘A’ in HEK293
at rs2279590 comprises of HSF1 protein. However, absence of any transfected cells, under heat shock or MG-132 treatment drasti-
super-shift could be because the antibodies are bound to the epit- cally reduced the reporter activity; suggesting a loss of enhancer
ope of proteins necessary for DNA-protein complex formation. effect with rs2279590/A allele. Although a similar effect was also
Super-shift assays with other HSFs (HSF2 and HSF4) did not dis- found in non-treated cells growing at normal conditions but the
solve the shift completely, suggesting the binding complex con- difference was not found to be significant (Supplementary
sists of HSF1 but not HSF2 and HSF4 (Supplementary Material, Fig. Material, Fig. S3). To understand allele specific enhancer activity in
S2B). Additionally, to validate the binding of HSF1 to rs2279590 absence of HSF1, luciferase activity in HEK293 transfected cells
in vivo, chromatin immunoprecipitation (ChIP) assays were was checked after knocking down HSF1 expression by a pool of
4524 | Human Molecular Genetics, 2017, Vol. 26, No. 22
HSF1 specific siRNA (Supplementary Material, Fig. S4). As overexpression has a protective role in a variety of neurodege-
expected, no significant difference (P ¼ 0.31) was observed in the nerative disorders caused due to accumulation of misfolded
reporter activity between constructs containing either allele ‘G’ or proteins (15–17). Further, aberrant expression of HSF1 has also
‘A’ in HSF1 knocked down cells (Fig. 4C). This suggests binding of been implicated in neurodegeneration (18). Assuming an
HSF1 to allele ‘A’ abrogates the enhancer effect of the locus. involvement of proteotoxic stress in anterior eye tissues of PEX-
affected subjects, we checked the expression of HSF1.
Interestingly, HSF1 level was found to be significantly upregu-
Deletion of rs2279590 locus alters the expression of PTK2B lated in both lens capsule and conjunctiva of PEXS-affected
and EPHX2, known modulators in AD pathogenesis individuals (Fig. 6). However, we didn’t find any upregulation of
HSF1 in PEXG-affected individuals; a later stage of PEX with
GWAS studies have previously shown that the region containing
degenerated optic nerve head (ONH) cells.
CLU and PTK2B is situated within a high-risk locus for AD patho-
genesis (14). As rs2279590 is also within the same high-risk locus,
we wanted to investigate the effect of this variant in regulating Discussion
the neighbouring genes around CLU. For this, expression of nine
Clusterin, a multifunctional protein has divergent roles from
crucial genes clustered around 5’ and 3’ ends of CLU gene, within
being cytoprotective to cytotoxic. Extensive studies have reported
a stretch of around 558 kb was analysed in HEK293115/ cells
the association of CLU variants with the risk of developing vari-
lacking the regulatory locus (Fig. 5A). One of the genes, CHRNA2
ous diseases like PEX, diabetes and AD (3,4,6,7). This indicates
did not show any expression in control HEK293 cells and is there-
that a complex mechanism employed by CLU is responsible for
fore omitted from analysis. Interestingly, a significant downregu-
the progression of these diseases. Understanding the role of risk
lation in the expression of two upstream genes, PTK2B and EPHX2
variants within CLU can help in better diagnosis and treatment of
was observed in HEK293115/ cells as compared with control
such disorders. In this study, we aimed to characterize the func-
HEK293 cells (Fig. 5B). This indicates that the genomic region
tional significance of a risk variant, rs2279590, housed in 7th
around rs2279590 has an extended regulatory role with a wide-
intron of CLU gene. This particular variant has been reported to
spread enhancer effect on PTK2B and EPHX2.
be a risk factor for both AD and PEX (3,4,6). We, therefore,
intended to understand the functional role of this SNP.
rs2279590 enhancer element also regulates two distal HSF1 to the protective allele ‘A’ abolishes the enhancer effect of
genes: PTK2B and EPHX2 the locus, which consequent to decreased expression of the target
genes; thereby implying a lowered risk of developing PEX or AD.
GWAS studies have previously shown an association of PTK2B-
CLU loci with AD (14). Protein tyrosine kinase 2 beta or PTK2B,
belongs to a non-receptor protein kinase family involved in cal- Summary and Conclusion
cium induced regulation of ion channels and MAPK pathway
activation. It plays a role in inducing phosphorylation of GSK3 With this study, we confirm that a PEX/AD associated risk var-
(Glycogen synthase kinase 3) which then promotes Tau fibrillar iant, rs2279590, resides within an enhancer element and regu-
pathology. In hippocampus of AD-affected individuals, height- lates the expression of three candidate genes CLU, PTK2B and
ened phosphotyrosine immunoreactivity was found in the EPHX2, which were previously known to be modulators in the
neuritic plaques, tangle-bearing neurons and microglia, which progression of AD. Increase in CLU expression by the risk allele
are characteristic features of increased PTK2B activity (29). at rs2279590 provides a mechanistic insight into the cytotoxic
Accumulation of PTK2B was also reported in the early event of role of CLU in PEXG individuals. Future studies are needed to
AD pathogenesis with progressive Tau pathology. However, it is explore the role of PTK2B and EPHX2 genes in PEX pathogenesis.
unclear whether accumulation of PTK2B in aiding the Tau path-
ology, is a consequence or a cause (30). In this study, we found
decreased levels of PTK2B mRNA in cells with the deleted Materials and Methods
enhancer element containing rs2279590; suggesting a distal Study subjects recruitment
enhancer effect of the locus over PTK2B expression. Although,
regulatory effect of rs2279590 on PTK2B was not established in This study was approved by the ethics review boards of the
eQTL data from GTEx project, it may depend on various con- National Institute of Science Education and Research and JPM
founding factors including the type of tissue analysed. Further, Rotary Club of Cuttack Eye Hospital and Research Institute,
earlier reports have shown an increased PTK2B expression with India and adhered to the tenets of the Declaration of Helsinki.
AD associated risk alleles at rs28834970 and rs2718058 (a cis- All participants underwent a detailed ocular examination,
pQTL within PTK2B locus and a trans-pQTL within NME8 locus) including slit lamp, ocular biometry, Goldman applanation
(31). Cumulatively, these results suggest that AD risk variants tonometry, þ90 D biomicroscopic fundus evaluation and 4 mir-
are capable of regulating PTK2B expression to promote AD ror gonioscopy. Inclusion and exclusion criteria for the grouping
pathogenesis. of control and PEX-affected individuals were followed as
Deletion of rs2279590 locus also leads to downregulation of reported previously (3). Anterior eye tissues (lens capsules and
EPHX2, another candidate gene for AD pathogenesis. EPHX2 or conjunctiva) from PEX-affected study subjects and age matched
epoxide hydrolase-2 metabolises epoxyeicosatrienoic acids which controls were collected in RNAlater stabilisation solution
are proven to be neuroprotective. Upregulation of hydrolase activ- (Invitrogen, USA) during cataract surgery and stored at 80 C
ity of soluble EPHX2 leads to an increased OGD-induced (oxygen- until further use.
glucose deprived) neuronal cell death (32). Recently, EPHX2 too has
been associated with AD as a risk factor through GWAS (33).
According to all the cumulative data known till date, a model Cell culture
is being proposed as shown in Figure 7 that rs2279590 resides The human cell line, HEK293 was cultured in HiGlutaXL
within a regulatory genomic region and with risk allele ‘G’ shows Dulbecco’s Modified Eagle Medium, High Glucose (AL007G) with
a widespread enhancer effect on three PEX/AD risk associated 10% fetal bovine serum (RM9952) and 1% penicillin (100 U/ml)
candidate genes; CLU, PTK2B and EPHX2. However, binding of and streptomycin (0.1 mg/ml) (A001), maintained at 37 C and 5%
Human Molecular Genetics, 2017, Vol. 26, No. 22 | 4527
CO2. All cell culture chemicals were procured from HiMedia, biotin 5’-end labelling [Integrated DNA Technologies (IDT), Iowa,
Mumbai, India. USA]. Sequences of the oligomers are listed in Supplementary
Material, Table S1. Annealing of complementary oligos was done
by incubating them at 95 C for 5 min, followed by step-cooling to
Luciferase reporter assays room temperature.
For reporter assays, pGL4.23 luciferase reporter vector with min- Nuclear extract from heat shocked (1 h at 42 C) HEK293 was
imal promoter and pGL4.74 renilla vector were used (Promega, prepared by using NE-PER kit (Thermo Fisher Scientific) and pro-
Madison, Wisconsin). Genomic DNA was extracted from tein concentrations were estimated by Bradford’s assay. EMSA
peripheral blood leucocytes of the study subjects by phenol- was performed using LightShift Chemiluminescent EMSA Kit
chloroform extraction method. An intronic region of 201 bp (Thermo Fisher Scientific) following the manufacturer’s instruc-
surrounding rs2279590 variant (harbouring either ‘AA’ or ‘GG’ tions. DNA-protein binding assays were carried out using 3 mg of
genotype at the polymorphic site) was PCR-amplified using a total nuclear extract and 100 fmol of biotinylated annealed oli-
specific primer pair (Supplementary Material, Table S1) from the gonucleotides for each 20 ml total reaction volume. Competitive
described previously (34). A pair of sgRNA was designed (http:// digested chromatin was incubated overnight at 4 C with 5 mg of
crispr.mit.edu/; date last accessed September 7, 2016) with a lit- ChIP validated HSF1 antibody (Santa Cruz Biotechnology, sc-
tle off-target specificity to delete a 115 bp region around the 9144X). As negative control, normal Rabbit IgG (1 ml) provided in
SNP, rs2279590 (Supplementary Material, Table S1). Annealed the ChIP kit was used. Later, 20 ml of ChIP grade protein A/G plus
oligonucleotides were phosphorylated and ligated into BbsI agarose beads was added to each IP and incubated for 2 h at 4 C
digested PX459 (Plasmid #62988; Addgene, Cambridge, with continuous rocking. IP complexes were then washed twice
Massachusetts). with each wash buffer supplied in the kit. Elution of the IP com-
At 50% confluency, the HEK293 cells were transfected with plex was done by incubating in IP elution buffer at 65 C for 1 hr.
2.5 mg of each CRISPR construct (with sgRNA1 and sgRNA2) Reverse crosslinking and protein digestion were done by NaCl
simultaneously using lipofectamine 2000 (Invitrogen, Carlsbad, and Proteinase K, respectively. Subsequently, DNA was purified
California). Transfected cells were selected after 24-h post- by using DNA Clean-Up column supplied in the kit. Fold enrich-
transfection in complete media supplemented with 2.5 mg/ml ment of the target region in IP compared with that of input was
puromycin. Single cell clones were then isolated and cultured assayed by a specific primer set (Supplementary Material, Table
by dilution cloning and subsequently genomic DNA was iso- S1) through quantitative real time PCR (qRT-PCR).
lated using DNeasy Blood & Tissue Kit (Qiagen, Hilden,
Germany) following manufacturer’s instructions. Screening was
Knockdown assays
done to find out homozygous deletion of genomic region con-
taining rs2279590 by using a specific set of primers siRNA pool for targeting HSF1 was procured from Santa Cruz
(Supplementary Material, Table S1) outside the targeted sites Biotechnology (sc-35611) and transfection was done as per man-
and subsequently confirmed by sequencing. Positive clones ufacturer’s protocol. For luciferase assays, reporter vector and
were used for subsequent experiments. siRNA pool were transfected together, cell lysates were pre-
pared after 36 h and subsequently checked for reporter activity.
16. Fujimoto, M., Takaki, E., Hayashi, T., Kitaura, Y., Tanaka, Y., 26. Jiang, Y.Q., Wang, X.L., Cao, X.H., Ye, Z.Y., Li, L. and Cai, W.Q.
Inouye, S. and Nakai, A. (2005) Active HSF1 significantly sup- (2013) Increased heat shock transcription factor 1 in the
presses polyglutamine aggregate formation in cellular and cerebellum reverses the deficiency of Purkinje cells in
mouse models. J. Biol. Chem., 280, 34908–34916. Alzheimer’s disease. Brain Res., 1519, 105–111.
17. Neef, D.W., Jaeger, A.M. and Thiele, D.J. (2011) Heat shock 27. Gouras, G.K. (2013) mTOR: at the crossroads of aging, chaper-
transcription factor 1 as a therapeutic target in neurodege- ones, and Alzheimer’s disease. J. Neurochem., 124, 747–748.
nerative diseases. Nat. Rev. Drug Discov., 10, 930–944. 28. Goetzl, E.J., Boxer, A., Schwartz, J.B., Abner, E.L., Petersen,
18. Kim, E., Wang, B., Sastry, N., Masliah, E., Nelson, P.T., Cai, H. R.C., Miller, B.L., Carlson, O.D., Mustapic, M. and
and Liao, F.F. (2016) NEDD4-mediated HSF1 degradation Kapogiannis, D. (2015) Low neural exosomal levels of cellular
underlies alpha-synucleinopathy. Hum. Mol. Genet., 25, survival factors in Alzheimer’s disease. Ann. Clin. Transl.
211–222. Neurol., 2, 769–773.
19. Szymanski, M., Wang, R., Bassett, S.S. and Avramopoulos, D. 29. Kohler, C., Dinekov, M. and Gotz, J. (2013) Active glycogen
(2011) Alzheimer’s risk variants in the clusterin gene are synthase kinase-3 and tau pathology-related tyrosine phos-