Background Heat shock proteins (Hsps) play a role in the delivery and presentation of antigenic peptides and are thought to be involved in the pathogenesis of multifactorial diseases.
Objective To investigate genes encoding cytosolic Hsp70 proteins for associations of allelic variants with systemic lupus erythematosus (SLE).
Methods Case–control studies of two independent Caucasian SLE cohorts were performed. In a haplotype-tagging single-nucleotide polymorphism approach, common variants of HspA1L, HspA1A and HspA1B were genotyped and principal component analyses were performed for the cohort from the Oklahoma Medical Research Foundation (OMRF). Relative quantification of mRNA was carried out for each Hsp70 gene in healthy controls. Conditional regression analysis was performed to determine if allelic variants in Hsp70 act independently of HLA-DR3.
Results On analysis of common genetic variants of HspA1L, HspA1A and HspA1B, a haplotype significantly associated with SLE in the Erlangen-SLE cohort was identified, which was confirmed in the OMRF cohort. Depending on the cohorts, OR ranging from 1.43 to 1.88 and 2.64 to 3.16 was observed for individuals heterozygous and homozygous for the associated haplotype, respectively. Patients carrying the risk haplotype or the risk allele more often displayed autoantibodies to Ro and La in both cohorts. In healthy controls bearing this haplotype, the amount of HspA1A mRNA was significantly increased, whereas total Hsp70 protein concentration was not altered.
Conclusions Allelic variants of the Hsp70 genes are significantly associated with SLE in Caucasians, independently of HLA-DR3, and correlate with the presence of autoantibodies to Ro and La. Hence, the Hsp70 gene locus appears to be involved in SLE pathogenesis.
This paper is freely available online under the BMJ Journals unlocked scheme, see http://ard.bmj.com/info/unlocked.dtl
Statistics from Altmetric.com
Heat shock proteins (Hsps) have been discovered as a result of their inducible expression in response to endogenous and exogenous stimuli such as elevated temperature, osmotic shock and the presence of cytotoxic agents.1 Functionally, most Hsps act as molecular chaperones by selectively recognising and binding non-native proteins, thereby preventing irreversible aggregation under physiological and stress conditions.2 By virtue of their peptide binding ability, this protein family modulates antigen processing and presentation and also contributes to immune responses against pathogens.3,–,5 Hsps are fundamental regulatory elements of cellular networks as they interact with a large number of proteins such as kinases, transcription factors and several other proteins influencing key steps in protein homoeostasis, cell growth, division, apoptosis and development.6 7 It has been shown that mutations and polymorphisms in chaperones, particularly Hsp70, are associated with several human, mostly immune-mediated, diseases such as multiple sclerosis,8 Crohn's disease,9 Alzheimer's disease,10 Grave's disease,11 insulin-dependent diabetes mellitus12 and systemic lupus erythematosus (SLE).13 14 The chronic inflammatory autoimmune disease SLE involves multiple organs, including skin, joints, kidneys, brains, the cardiovascular system and serosal membranes. Immunologically, SLE is characterised by the presence of antinuclear autoantibodies, especially against double-stranded DNA and nucleosomes, activation of the complement system during flares, and type I interferon secretion.15 16 The aetiology of SLE remains elusive, but an interplay of genetic and environmental factors is considered to ultimately cause immune dysregulation, resulting in clinical symptoms. Many genetic variants on different chromosomes together with more than 31 different candidate genes have been reported and were confirmed as susceptibility loci for SLE by replication studies.17,–,20 A well-known susceptibility locus for SLE is the major histocompatibility complex (MHC) region on chromosome 6p21, which is linked to multiple autoimmune diseases.21 Mainly alleles of the MHC class II region such as human leucocyte antigen (HLA)-DR2 and HLA-DR3 were found to be associated with SLE, especially in Caucasians22 23; however, the exact mechanism by which these HLA molecules contribute to the immunopathogenesis is unknown. Owing to the extent of local linkage disequilibrium (LD), it is not yet clear whether the susceptibility to the disease is directly and exclusively associated with the products of specific MHC alleles or of a combination of certain MHC alleles with alleles of other linked genes, or just with certain polymorphisms of linked genes. Therefore, other genes mapping to the MHC loci have been considered to be additional candidates for disease susceptibility such as complement C2 and C4, tumour necrosis factor α,24 21-hydroxylase and the superkiller viralicidic activity 2-like (Saccharomyces cerevisiae) gene (SKIVL2).25 In humans, three members of the Hsp70 gene family, referred to as Hsp70-1 (HspA1A), Hsp70-2 (HspA1B) and Hsp70-Hom (HspA1L), have been mapped to this MHC class III locus.26 Proteins of the Hsp70 family are commonly among the most highly conserved proteins known. Subtle genetic variations can result in diversified thermotolerance between species27 and may also be important determinants for the specific propagation of diseases such as SLE.
In this study, we examined a potential genetic association of the Hsp70 locus on chromosome 6p21.33a with SLE. In two independent Caucasian case–control studies, we found one haplotype to be significantly associated with SLE.
Patients and methods
Study population and clinical evaluation
A cohort of 212 unrelated well-documented German patients with SLE (184 women, 28 men) from the Department of Internal Medicine 3, University of Erlangen–Nuremberg, were randomly selected irrespective of stage or severity of the disease. All patients fulfilled the criteria of the American College of Rheumatology for the diagnosis of SLE.15 16 A median of 13 visits was recorded for each patient, varying from 1 to 54, and the follow-up period extended from 1985 to 2006. Autoantibodies against double-stranded DNA were present in 99% of our SLE cohort, at least at some time during the course of the disease. All abnormalities associated with SLE were recorded from the first visit either in Erlangen or from a well-documented history and at each follow-up visit. A group of 223 healthy unrelated donors volunteered as controls. The study was approved by the ethics committee of the University of Erlangen–Nuremberg, and written informed consent was obtained from all participants.
A second cohort consisting of 934 unrelated American patients with SLE and normal healthy donors of European descent from the Lupus Family Registry and Repository (http://lupus.omrf.org) was kindly provided by the Oklahoma Medical Research Foundation (OMRF). Data on autoantibody titres were derived from both medical records and Clinical Laboratory Improvement Act-approved laboratory data.
Genomic DNA was isolated from 0.5 ml peripheral blood using standard procedures. Identification of haplotypes was performed with the Haploview 3.2 CI method based on the HapMap Data release 21 (http://www.hapmap.org/). Five different tagging single-nucleotide polymorphisms (SNPs) (rs2075800, rs2227956, rs1043618, rs2763979, rs3115673) were selected on the basis of LD to discriminate the major haplotypes of the Hsp70 gene locus that are present with a frequency above 1% in the Centre d'Etude du Polymorphisme Humaine (CEPH) population. Genotyping was carried out in an iCycler real-time PCR setup (Bio-Rad, Munich, Germany) using TaqMan SNP Genotyping Assays and TaqMan Universal PCR Master Mix, No AmpErase UNG (Applied Biosystems, Darmstadt, Germany).
Venous blood was collected from healthy donors and patients with SLE in 20 ml heparin-containing tubes for isolation of cells, or without anticoagulant for serum (Monovette; Sarstedt, Nuembrecht, Germany). Serum samples were centrifuged at 2000 g for 30 min at 4°C and stored at −20°C.
The human cell line HEK293T (CRL-11268; ATCC, Wesel, Germany) and peripheral blood mononuclear cells (PBMCs) were grown in Dulbecco's modified Eagle medium containing 10% fetal bovine serum, 2 mM L-glutamine, 100 U/ml penicillin, 100 µg/ml streptomycin (all from Invitrogen Life Technologies, Karlsruhe, Germany) at 5% CO2 and 37°C. PBMCs were isolated as described previously.28 Heat shock treatment was carried out for 2 h at 42°C followed by a recovery period (37°C) either for 3 h for subsequent RNA isolation or overnight for western blotting.
Sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and western blotting
Isolated PBMCs lysed in 50 µl SDS-PAGE sample buffer were subjected to SDS-PAGE as described.29 Proteins were transferred on to a nitrocellulose membrane (Millipore, Bedford, UK) using semidry apparatus (Serva, Heidelberg, Germany) according to the manufacturer's recommendations. After blocking (5% (w/v) milk powder, 10 mM Tris/HCl, 100 mM NaCl, 0.05% (v/v) Tween-20, pH 7.4), primary antibodies against Hsp72 (1:30 000) (SPA-812; Biomol, Hamburg, Germany) and β-actin (1:1000) (Sigma-Aldrich, Hamburg, Germany) were applied. After washing, blots were incubated with secondary IRDye 800CW goat anti-rabbit antibody (1:10 000) (LI-COR Bioscience, Bad Homburg, Germany). Integrated intensity of the bands was quantified using the Odyssey infrared imaging system (LI-COR Bioscience).
RNA isolation and quantitative PCR
Total RNA was extracted from 2 × 106 PBMCs using the RNeasy kit (Qiagen, Hilden, Germany). cDNA was prepared with Moloney murine leukaemia virus reverse transcriptase (Promega, Mannheim, Germany) and oligo(dT)20 primers according to the manufacturer's instructions. Quantification of cDNA was carried out using ABsolute QPCR SYBR Green ROX Mix (Thermo Fisher Scientific, Schwerte, Germany) in a StepOnePlus Real-Time PCR System (Applied Biosystems) with the following primers: HspA1A, 5′-GCTGCTGCGACAGTCCACTAC-3′ and 5′-CCCTGCTCTCTGTCGGCTCGGC-3′; HspA1B, 5′-AGCTGCTGCGAGGGTCCGCTTC-3′ and 5′-GTGCCCTGCTCTGTGGGCTCC-3′; HspA1L, 5′-ACAGGGTATGTGCCTGGAAG-3′ and 5′-GGTGGTGCATGAGACTCCTT-3′; GAPDH, 5′-CATGAGAAGTATGACAACAGCCT-3′ and 5′-AGTCCTTCCACGATACCAAAGT-3′. Relative quantification of cDNA was based on the comparative Ct method using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as endogenous control. ΔCt levels were calculated in triplicate for four different donors of each genotype and normalised to one randomly selected donor homozygous for cytosine.
Principal component analysis (PCA) was performed using Eigensoft30 on 905 of the 934 Oklahoma samples using case–control data from 17 875 SNPs collected from a separate candidate gene study. After removal of genetic outliers (12 cases and 7 controls), 493 cases and 393 controls remained for analysis. We then performed genomic control analysis to calculate the inflation factor (λ), which produced a value of 1.21 for these samples (the value for the pretrimmed sample was 1.24). After trimming of the sample, we performed a logistic regression analysis using the PCA values as covariates, using PLINK31 to remove the final source of confounding via admixture, which produced a corrected p=0.0145 (OR=1.292) at rs2763979 and p=0.4762 (OR=0.903) at rs3115673 in these samples.
Conditional regression analysis
For the OMRF samples, a surrogate SNP (rs3129860) was chosen from the OMRF Lupus Large Association Study 1 (LLAS1) data that is in strong LD (r2=0.95) with the HLA-DRB1 tagging SNP (rs701831). Conditional regression analysis was performed with the surrogate SNP and SNP Hsp70-4 on the OMRF samples using PLINK31.
Statistical analysis for the case–control study was carried out with the Haploview V3.32 software package.32 Haplotype-based association testing was carried out using PLINK.31 The Fisher exact test and Mann–Whitney U test were performed with Prism 4 software. p Values <0.05 were considered significant.
Genetic association of Hsp70 haplotypes with SLE
To investigate a possible contribution of Hsp7 polymorphisms to SLE susceptibility, we performed a case–control association study followed by an independent replication study. As shown in figure 1A, the three intronless Hsp70 genes, HspA1L, HspA1A and HspA1B, are located on chromosome 6p21.33a. On the basis of information on the CEPH population provided by the HapMap database, we selected five haplotype-tagging SNPs to define haplotypes with more than 1% frequency (figure 1B) and with high LD depicted in the LD plot (figure 1C). All observed SNP allele frequencies were concordant with those reported in the NCBI dbSNP database, and the genotype frequencies were in Hardy–Weinberg equilibrium. In the Erlangen cohort, polymorphism Hsp70-4 was identified to be significantly associated with SLE (p=0.016). The major allele of this SNP, cytosine, was present in 61% of patients with SLE compared with 69% in controls, suggesting either cytosine as a protective allele or thymine as part of a risk allele for SLE. Two reconstructed haplotypes were also strongly associated with SLE, with C–A–C–T–G being more frequently (p=0.0022) present in patients with SLE (26%) than in controls (17%). Conversely, the haplotype C–A–C–C–G (p=0.0444) was found more frequently in healthy individuals (table 1). In addition, we analysed whether haplotypes have independent or interactive effects. By testing each haplotype against all others, we found that only the risk-conferring (p=0.0018) and protective (p=0.045) haplotypes are significant in the haplotype-specific test. Even when all haplotype effects were tested together in an interaction model, only this combined haplotype was found to be significant.
Genotyping of the second cohort, OMRF, was performed only for SNPs Hsp70-4 and Hsp70-5, because the two significantly associated haplotypes C–A–C–T–G and C–A–C–C–G can be identified by these SNPs. In this independent replication cohort, the genetic associations of both SNP Hsp70-4 and the representative haplotypes T–G and C–G, respectively, were confirmed with even higher significance (table 2). Combined analysis of both cohorts reached a p value of 3.64×10−7 for the haplotype T–G and 0.0002 for SNP Hsp70-4.
To determine the effect size of SNP Hsp70-4 and the Hsp70 haplotypes for SLE, ORs were calculated using unstratified genotyping data. For our newly identified SNP, Hsp70-4, the OR for donors bearing a T allele is 1.43 (95% CI 1.08 to 1.89) and 1.3 (95% CI 1.07 to 1.58) for the Erlangen and OMRF cohort, respectively. When the analysis was performed for haplotypes, the OR for patients bearing one T–G haplotype compared with those homozygous for C–G is 1.88 (95% CI 1.18 to 3.02) in the Erlangen cohort and 1.43 (95% CI 1.04 to 1.96) in the OMRF cohort. When homozygous donors for T–G and C–G are compared, the OR is 2.64 (95% CI 1.02 to 6.85) and 3.16 (95% CI 1.67 to 5.97) in the Erlangen and OMRF cohort, respectively (table 3).
Conditional regression analysis
The distance between the three Hsp70 genes, located in the MHC class III region, and HLA-DRB1 (MHC class II) is ∼1 Mb. The LD plot for European-derived samples genotyped in the International HapMap Project (online supplementary figure 1) shows incomplete linkage in this region. However, owing to long-range linkage in this entire locus, we performed conditional regression analysis to investigate whether the genetic association of SNP Hsp70-4 is independent of HLA-DRB1. One surrogate SNP (rs2763979) that shows strong LD with SNP rs701831 in HLA-DRB1 was chosen from the OMRF LLAS1 genotyping data, as the HLA-DRB1 tagging SNP itself was not included in LLAS1 dataset. We found that SNP Hsp70-4 (rs2763929) acts independently of the HLA-DRB1 surrogate (p=0.019) and vice versa (p=0.045).
Correlation of the haplotype T–G with clinical variables in patients with SLE
To examine a potential effect of the risk haplotype T–G on disease manifestation, we correlated clinical variables of the Erlangen cohort with the associated haplotypes in a cumulative analysis. Patients with haplotype T–G were significantly more frequently positive for anti-Ro and anti-La autoantibody responses, lupus anticoagulant and oronasal ulcers than those with haplotype C–G (online supplementary figure 2A). In contrast, when SLE patients bearing HLA-DR3 alleles were compared with those lacking an HLA-DR3 allele (HLA-DRx/x), no statistically significant correlations with clinical or immunological disease variables could be detected (online supplementary figure 2B). These data also indicate that the Hsp70 haplotype T–G contributes independently of HLA-DR3 to the manifestation of SLE. As depicted in table 4, patients either heterozygous or homozygous for allele T or for the haplotype T–G were more frequently positive for anti-Ro and anti-La autoantibodies in both SLE cohorts.
Increased mRNA expression of HspA1A in PBMCs from healthy donors homozygous for thymine at SNP Hsp70-4
SNP Hsp70-4 is located between the genes for HspA1A and HspA1B, thus allelic variations at this position may influence the expression of HspA1A and HspA1B. To prevent the influence of medication, we used PBMCs from 12 healthy donors who were either homozygous for thymine or cytosine or were heterozygous at SNP Hsp70-4. Figure 2A,B shows a significantly increased expression of HspA1A mRNA in healthy donors homozygous or heterozygous for thymine compared with homozygous carriers of cytosine, both under physiological conditions (figure 2A) and after heat stress (figure 2B). In contrast, we detected no significant differences in mRNA levels for either HspA1B or HspA1L under all conditions tested. In accordance with other reports,33 HspA1A and HspA1B were found to be strongly heat inducible, whereas HspA1L was weakly induced after heat stress (data not shown). To further investigate the expression of the three Hsp70 genes (HspA1L, HspA1A and HspA1B), the total amount of all three Hsp72 proteins was determined by western blot (figure 2C), as gene products are 99% identical and currently indistinguishable by available antibodies. The overall Hsp72 protein level was not altered in PBMCs from donors with thymine or cytosine at SNP Hsp70-4, when kept at 37°C or after heat shock. To examine a potential direct influence of genotypes of SNP Hsp70-4, located upstream of the HspA1B promoter, on the expression of HspA1B, the promoter activity was analysed by reporter assays (online supplementary methods and supplementary figure 3). The promoter of HspA1B was highly inducible by heat stress and to a lower extent by osmotic shock. In agreement with mRNA quantification in PBMCs of healthy donors, similar luciferase activities for both promoter constructs were reproducibly observed for all conditions tested, arguing against a direct effect of SNP Hsp70-4 on the promoter activity of HspA1B.
In our study, we identified one haplotype (T–G) in the Hsp70 locus as a potential genetic risk factor for SLE using a haplotype-tagging SNP approach and confirmed the genetic association in an independent case–control cohort. In a cumulative analysis, this haplotype correlated in SLE patients with the presence of clinical variables such as autoantibodies to Ro and La, as well as lupus anticoagulant. To our knowledge this is the first genetic association study on autoimmune diseases such as SLE that takes into account the entire Hsp70 locus on 6p21.33 with respect to the underlying LD. Previous studies aimed to identify individual SNPs, such as the 8.5 kb PstI polymorphism at position 1267 within the HspA1B coding region representing a silent mutation. This SNP has been found to be associated with SLE in African–Americans and in a Spanish, but not a Mexican Mestizo, population.13 14 34 Polymorphisms in HspA1L (2437, coding), HspA1A (−110, 5′ flanking; and 190, 5′ untranslated region) and HspA1B (1267, coding) have been discovered to correlate with increased longevity in Danish populations.35 36 Most Hsps can positively influence cellular maintenance and repair mechanisms through their chaperone activity.37,–,39 Therefore Hsps are generally recognised as protective elements increasingly gaining attention in relation to diseases and immunological processes.6 Associations of polymorphisms within the three Hsp70 genes with Ménière's disease,40 noise-induced hearing loss,41 high-altitude pulmonary oedema42 and coronary heart disease43 have recently been discovered. In these reports, major importance was determined for two SNPs in HspA1A that are located at positions −110 (rs1008438) and +190 (rs1043618), the latter also being part of our study (Hsp70-3). As we found a genetic association for the haplotypes T–G and C–G (Hsp70-4 and Hsp70-5) with SLE, we stratified the samples according to SNP Hsp70-4 for functional analyses and discovered differences in mRNA regulation for HspA1A (figure 2A,B). In spite of divergent mRNA concentrations for HspA1A, the protein expression level of all three Hsp70 genes was not altered in PBMCs from different donors either under steady-state conditions or after heat shock (figure 2C), which might be a result of negative transcriptional feedback regulation. Whether increased transcription of HspA1A mRNA is functionally relevant for the pathogenesis of SLE remains to be elucidated. Because of the extent of LD, we are currently unable to rule out indirect effects that might arise from potentially causal variants of Hsp70-3 or other unknown factors. As we have analysed the expression of Hsp70 genes in PBMCs, a tissue-specific expression of Hsp70 is possible and may be relevant for SLE. Because of long-range linkage in the lupus-associated DR3 ‘autoimmune’ extended haplotype, we performed conditional regression analysis to investigate if the associated SNP Hsp70-4 acts independently of HLA-DR3. Tracking HLA-DR3 using a surrogate SNP, we found SNP Hsp70-4 to be independent, and vice versa. These results are in accordance with recent investigations carried out by high-density screening of the MHC genes, producing strong evidence for independent susceptibility regions in the MHC locus.44 In addition, we showed that the Hsp70 haplotype T–G, as well as the SNP Hsp70-4, correlated in both cohorts with the presence of autoantibodies to Ro and La in patients with SLE.
In summary, our data indicate an important role for Hsp70 molecules in the aetiology and pathogenesis of SLE, consistent with the involvement of Hsp70 molecules in antigen delivery and presentation as well as protein folding and homoeostasis.
We thank Gerhard Groer for critical reading of the manuscript, Rüdiger Müller for help with statistical analyses, and Birgit Lauer for HLA genotyping. Portions of the samples used in this study were obtained from the Lupus Family Registry and Repository (AR62277) (see http//lupus.omrf.org).
Funding This study was supported by the Interdisciplinary Center of Clinical Research (IZKF) of the University Hospital Erlangen, project N2. This work has been supported by the National Institutes of Health (grants AI24717, AR62277, AR42460, AR49084, RR20143, AR053483, N01-AI-50026, P20RR015577), the Mary Kirkland Scholarship, the Alliance for Lupus Research and the US Department of Veterans Affairs.
Competing interests None.
Patient consent Obtained.
Ethics approval This study was conducted with the approval of the ethics committee of the Medical Faculty of the University of Erlangen-Nuremberg.
Provenance and peer review Not commissioned; externally peer reviewed.