OBJECTIVE To determine whether glutathione S-transferase GSTM1, GSTM3, GSTT1, and GSTP1 genotypes influence susceptibility or outcome in rheumatoid arthritis (RA).
METHODS 277 RA patients were compared with 577 controls to examine any associations between GST genotypes and susceptibility to RA. The effect of genotypes on outcome (Larsen and functional scores) and time integrated acute phase responses (erythrocyte sedimentation rate and C reactive protein) was assessed in 122 patients with disease duration of 5–10 years. GST and HLA-DRB1 genotypes were determined using polymerase chain reaction based assays. Data were analysed using multiple regression analysis with correction for age, sex, disease duration, and the DRB1 associated shared epitope (SE) and rheumatoid factor (RF) positivity where appropriate.
RESULTS TheGSTM1*A/*B genotype was less common in RA cases (3 of 276) than in controls (22 of 591) (exact p=0.047), though significance was lost when adjustment was made for multiple comparisons. The Larsen score was higher (p=0.039) in the GSTM1 null patients (89.9) than those with other GSTM1 genotypes (74.7), and this was independent of the SE. Again, correction for multiple testing resulted in loss of significance. The difference in Larsen scores between patients homozygous or negative for the SE (87.9v 74.3) was similar to that between GSTM1 null and non-null patients. No associations between GSTM3 or GSTT1 genotypes and disease markers were identified although the association between GSTP1*B/*B and Larsen score approached significance (p=0.096).
CONCLUSION It is proposed that certain GSTs may influence susceptibility and radiological progression in RA and that this is independent of the effect of the HLA-DRB1 associated SE. The mechanism for this effect is presumed to be because of differences in the ability of various GST enzymes to utilise the cytotoxic products of oxidant stress. Although significance was lost after correction for multiple testing, the data indicate that further studies may be of value in RA to determine the influence of the GST and other genes involved in cellular protection against oxidative stress.
- rheumatoid arthritis
- glutathione S-transferase
- shared epitope
Statistics from Altmetric.com
Rheumatoid arthritis (RA) is characterised by chronic inflammation of synovial joints with progressive destruction of cartilage and bone. The genetic factors that determine susceptibility and outcome are unclear though the HLA system, particularly the HLA-DRB1 molecule, is important.1-3 This includes subtypes of DR4 (DRB1*0401, *0404, *0405, *0408), DR1 (DRB1*0101, *0102) and DRB1*1001 and *1402. All the HLA-DRB1 associated alleles encode a conserved amino acid sequence (QKRAA, QRRAA or RRRAA) known as the shared epitope SE.4 Several studies indicate that certain combinations of SE carrying alleles are associated with more severe disease as measured by radiological outcome, or extra-articular manifestations, or both.3 5-7 Other susceptibility or outcome candidate genes for RA have yet to be identified.
Inflammation is a central feature of RA with the resulting reactive oxygen species (ROS) causing oxidation of DNA and lipids giving rise to a variety of cytotoxic products including lipid and DNA hydroperoxides and alkenals.8 There is evidence implicating ROS and their products in the pathology of RA. For example, ROS is produced by phagocytes in the synovial fluid and pannus and by synovial endothelial cells during hypoxia-reperfusion events.9 This suggests that variations in host effectiveness in the detoxification of products of ROS activity is important. Accordingly, we propose that polymorphism in enzymes that detoxify ROS and their products may contribute to the wide variation seen in the amount of joint damage and functional impairment.
The widely expressed, glutathione S-transferase (GST) supergene family seems to be critical in cellular protection against ROS.8 10 These enzymes catalyse the conjugation of glutathione with electrophiles such as lipid hydroperoxides, 4-hydroxy non-2-enal, DNA hydroperoxides. Polymorphism has been identified in several GST genes. Three alleles are described for GSTM1;GSTM1*0 is deleted whileGSTM1*A andGSTM1*B encode monomers that form active enzymes.8 GSTM3 is biallelic withGSTM3*A and GSTM3*B,differing by a 3 base pair deletion in the latter that creates a motif for the YY1 transcription factor. Furthermore,GSTM3*B andGSTM1*A are in linkage disequilibrium. A null polymorphism has also been identified at GSTT1.8 In the case of GSTP1, three alleles have been identified apart from the wild type GSTP1*A; one contains an A-G transition at nucleotide +313 (GSTP1*B), one contains a C-T transition at +341 in addition to the A-G transition (GSTP1*C), while another contains only the C-T transition at +341 (GSTP1*D). The +313 transition has been shown to result in altered catalytic activity,11 12 though there is no evidence to date of a functional effect of the +341 transition. Therefore we have studied only the +313 variant alleles. We have designated homozygotes for this transition, GSTP1* B/*B, though it is possible that a proportion of these will have theGSTP1*C allele.
Support for the view that GST polymorphisms influence susceptibility/outcome to diseases with an inflammatory component comes from case control studies. Thus, GSTM1null is associated with increased production of anti-Ro antibodies in patients with systemic lupus erythematosus and GSTT1 null confers increased risk of ulcerative colitis.8 Pilot data also suggest that particular GSTM1/GSTM3 genotype combinations confer increased susceptibility to oxidant stress.10 There are few data on the recently identified GSTP1 polymorphism. On the basis of in vitro and molecular epidemiological studies, particular genotypes can be considered high risk in the context of detoxification of ROS products. These areGSTM1*0/*0,8 10 GSTT1*0/*0,8 10 GSTM3*A/*A,10 GSTP1*B/*B.11 We now describe studies to determine if these putatively high risk genotypes are associated with susceptibility or outcome in RA.
GST GENOTYPES AND SUSCEPTIBLITY TO RA
The association between GST genotypes and susceptibility was studied in 277 unrelated northern European white RA patients resident in north Staffordshire (40.4% male, median age 58.5, range 32.1–80 years) with a median disease duration of 11 years (range, 5–39). Of these, 60.5% were rheumatoid factor positive. They were recruited in a clinic established to examine the effects of slow acting anti-rheumatic drugs. Treatment was administered as clinically indicated. The ARA criteria of 1958 were recorded at the first presentation and have been used as the basic definition of the cohort. The 1987 ARA criteria were documented retrospectively from data in the case notes. The control group comprised 577 ethnically matched people (48% male, mean age 70 years) suffering varicose veins, hernias, haemorrhoids, or mild iron deficiency anaemia.10 Patients with inflammatory disorders such as ulcerative colitis, diabetes or asthma were excluded.
GST, DISEASE OUTCOME, AND ACUTE PHASE RESPONSE
The association between disease outcome, acute phase response, and GST was studied in a subset of 122 RA patients (38.5% male, median age 58.1 years, range 32.1–79.0) selected as having a disease duration of 5–10 years (median 8.3). This was to minimise the considerable between patient variation in disease duration (5–39 years) present in the total case group and to reduce the possible ceiling effect in radiological scores seen with increasing disease duration.13 These patients had been reviewed annually and their disease extensively characterised. Outcome measures were recorded at final review and consisted of assessment of functional status using the Health Assessment Questionnaire (HAQ)14 and radiographic outcome, obtained by scoring radiographs of the hands and feet using the standard radiographs of Larsen.15 The time integrated acute phase response was assessed by measuring erythrocyte sedimentation rate (ESR) and serum C reactive protein (CRP) concentrations at yearly intervals for at least five years and calculating the area under the curve (AUC) for these parameters. Division of the AUC by the number of years follow up gave the mean area under the curve (MAUC).16
GST AND HLA-DR GENOTYPING
Genotypes were identified from lymphocyte DNA. GSTM1 genotypes were defined using a PCR assay that identifiesGSTM1*0 homozygotes,GSTM1*A/*B heterozygotes and the GSTM1 A and GSTM1 B phenotypes. It does not distinguishGSTM1*0/*A andGSTM1*A/*A genotypes or the equivalent GSTM1 B genotypes.8 The GSTM3 *A/*A, *A/*B and *B/*B genotypes were identified by amplifying the exon 6/7 regions of GSTM3 and differentiating GSTM3*A fromGSTM3*B by digestion withMnlI.8 TheGSTP1 *A/*A, *A/*B and*B/*B genotypes were identified using primers to exon 5.10 The GSTT1*0 and expressing subjects were also identified using poymerase chain reaction.8 Presence of the HLA-DRB1 SE was identified in people using methods previously described.17 18
Genotype frequencies were compared between cases and controls using χ2 tests (StatXact Turbo). Disease severity and activity variables were first transformed to normality, where appropriate, and the effects of genotypes were assessed using multiple regression analysis with the addition of the independent variables, age, sex, disease duration, SE, and rheumatoid factor (RF) positivity where these were significant. Correction for potential multiple testing errors was performed using Holm’s procedure.19 This is less conservative than the widely used Bonferroni procedure and has been considered as being uniformly better.20 In each set of analyses six tests were performed (GSTM1, GSTT1, GSTM3, GSTP1, GSTM1+GSTM3, GSTM1 +GSTT1).
CASE CONTROL STUDY: GST AND SUSCEPTIBLITY
The clinical indices of the total group of 277 RA patients, and the 122 patients with disease duration 5–10 years were similar to other studies, and there was typical heterogeneity in terms of age of onset and disease duration. Table 1 summarises GST genotype frequencies. Differences were identified for GSTM1; the frequency ofGSTM1*A/*B in the cases was significantly lower than in the controls (χ2 1=4.7, exact p=0.047, odds ratio 0.28, 95% CI 0.054, 0.962). However, this was not significant after correction for potential multiple testing errors. The frequency of GSTM3*B/*B was lower in the cases though the difference did not achieve significance (exact p=0.12). No significant differences in the frequencies of GSTT1 or GSTP1 genotypes were identified.
ASSOCIATION OF GST GENOTYPES WITH MARKERS OF DISEASE OUTCOME AND ACUTE PHASE RESPONSE
Larsen scores together with HAQ, MAUCESR and MAUCCRP data from these patients are shown in table 2. Comparisons were made between individual GST genotypes, and also between putatively high risk genotypes and the corresponding complement (that is, GSTM1 null compared with non-null genotypes).
The highest mean Larsen score (89.9) was found in patients withGSTM1*0/*0, which was significantly different (p=0.039) to non-null people (74.7) after adjustment for age, sex, and disease duration. Correction for the presence of the SE did not change the association of Larsen score withGSTM1*0/*0 (p=0.045). Furthermore, this association was not confounded by RF positivity as addition of this variable to the multiple regression analysis had little impact (p=0.031). Adjustment for possible mutiple testing errors resulted in loss of significance. The lower Larsen score in non-null people seemed to be accounted for largely by people with the GSTM1 A phenotype (70.7) although the difference from GSTM1 B people (83.2) was not significant (p= 0.6) after correction for age, sex, and disease duration. No significant associations between GSTM1 genotypes and the age of onset, HAQ score or, measures of the acute phase response were identified.
GSTM3 and GSTT1
No significant associations between GSTM3 or GSTT1 genotypes and Larsen score, HAQ, MAUCESR or MAUCCRP were seen. Furthermore, no interactions between GSTM1 null and GSTT1 null or between GSTM1 null and GSTM3*A/*A were identified.
The differences in Larsen scores between individual GSTP1 genotypes were not significant, although an association between theGSTP1 *B/*B and the Larsen score approached significance (p=0.096; corrected for SE, p=0.064) when compared with non-*B/*B people. Additional correction for RF positivity as well as the SE achieved a significant difference (p=0.04), though significance was lost after adjusment for multiple testing. There were no associations between GSTP1 genotypes and age of onset, HAQ, MAUCCRP or MAUCESR, with or without correction for the SE and RF positivity.
COMPARISON WITH THE SE
As the presence of the SE has been associated with more severe disease we examined the association of the SE with markers of outcome and acute phase response, for comparison with the GST genotypes. The highest Larsen score (87.9) was found in SE homozygous patients, and the lowest (74.3) in patients negative for the SE (table3). Although not statistically significant (p=0.14, after correction for age, sex, and disease duration) the trend was consistent with previous studies. Interestingly the Larsen score for the SE homozygous patients (87.9) was similar to that for GSTM1 null patients (89.9), while the score for SE negative people (74.3) was essentially the same as that for patients with a positive GSTM1 allele (74.7). No significant differences were found between SE homozygous and SE negative patients for age of onset, HAQ, MAUCESR and MAUCCRP.
Correction for the presence of the SE in our multiple regression analysis indicated that it had little influence on the difference in Larsen scores between GSTM1 null and non-null patients. This is illustrated clearly in table 4 where GSTM1 null and non-null patients have been divided into SE+/+, SE+/−, and SE−/− groups. The frequency of SE−/− and SE+/+ people is not significantly different between GSTM1 null and non-null groups, and in GSTM1 null patients no statistical differences in Larsen score was found between patients with or without the SE. In non-null patients the difference between the groups was not significant although the trend suggested a lower Larsen score in SE−/− patients.
We have proposed that GSTM1, GSTM3, GSTT1, and GSTP1 are candidate genes for susceptibility and outcome in RA. We used a case-control approach to identify associations between GST genotypes and susceptibility to RA. Our purpose was exploratory and we realise there are dangers in making inferences on individual factors in multiple testing. Results are considered as hypothesis setting and needing confirmation in further studies. Our data suggest thatGSTM1 *A/*B, but not the GSTM1 A or GSTM1 B phenotypes, may be associated with reduced susceptibility. This indicates a gene dose effect; thus, GSTM1 *A/*B is protective because most people defined as GSTM1 A or GSTM1 B phenotypes are GSTM1*0/*A andM1*0/*B heterozygotes respectively.8 These data and others reporting that the genotype is protective against other oxidative stress related diseases such as basal cell carcinoma suggest that GSTM1 is important in cellular protection, though the mechanism is unclear.8Given the reported linkage disequilibrium with GSTM3 alleles, it is possible that the biological effect of GSTM1 is exerted through other mu class GST or another neighbouring gene. In our particular RA population the extent of linkage between GSTM1 and GSTM3 was no different to that in controls (data not shown).
We also assessed the influence of GST polymorphisms on outcome and the acute phase response in RA. A more severe radiological outcome, was associated with GSTM1 null, while the association withGSTP1 *B/*B approached significance. The association was independent of RF positivity and the presence of the SE, both of which have been associated with more severe disease.5 6 21-23 Interestingly the Larsen scores in GSTM1 null and GSTP1*B/*B patients were equivalent or greater than Larsen scores for patients with a double dose of the SE. A small but non-significant increase in Larsen score was seen in patients with both GSTM1 null and SE+/+. This preliminary study suggests that the GSTM1 null genotype may be an independent marker for development of more erosive disease. However, the association seems to be weak and did not achieve significance after correction for potential multiple testing errors. Thus its clinical usefulness may be limited and further studies are needed to determine more clearly the significance of this finding.
We did not identify any associations between GST genotypes and functional outcome assessed using the HAQ score. This may partly be explained by the fact that such scores include factors like pain perception, neuromuscular power and psychosocial factors that are difficult to assess precisely and may therefore, be difficult to correlate with single genotypes. Other factors may also cause discrepancy between radiologically assessed and functional outcomes. For example, Larsen score primarily assesses damage to cartilage and bone but does not necessarily reflect damage to other tissues and organs. Thus, in the joint, unknown levels of damage to tendons, ligaments and soft tissue, together with neurological changes and muscle wasting will all be important in determining clinical and functional outcome.
Previous studies have shown that measures of the acute phase response over 5–10 years are significantly correlated with radiological assessment of joint deterioration, and that persistently raised ESR and CRP values are associated with more radiological progression.16 However, there is also evidence that erosion and inflammatory response may follow diverging paths.24 In this study no associations between GST genotypes and MAUCCRP or MAUCESR were seen. We interpret this as indicating that the inflammatory response is not dependent on GST genotype, it being more probable that the GSTs serve to modify the effects of inflammation. Thus the radiological outcome in people with the same amount of inflammation over time may be different depending on GST status. Although one might expect to see a general association beween inflammatory markers and radiological outcome it would not necessarily be seen between these markers and the GST genotypes. The same applies to the DRB1 SE where the relation with the inflammatory response and radiological outcome is still far from clear.
Our results suggest that genes involved in protection from oxidant stress may influence the disease process in RA. The influence of these genotypes may contribute with other genes that exert a relatively weak influence.8 Though such weak effects may be real they tend to disappear when adjustment for multiple testing is carried out. It is unclear which GST substrates are important and why certain GST seem influential in some but not all diseases in which oxidant stress is a critical feature. This may reflect tissue specific patterns of gene expression. GSTP1 is of interest as it is widely expressed in human cells including synovial fibroblasts and lymphocytes.25 In contrast, expression of GSTM1 protein is more restricted, being found in lymphocytes but not in synovial fibroblasts.25
We gratefully acknowledge the generous support of the Arthritis Research Campaign, the Haywood Rheumatism Research and Development Foundation, the European Commission (contracts BMH4-CT-0396 and BMH4-CT96–0087) and Sheila Clarke and June Fisher, metrologists at the Staffordshire Rheumatology Centre.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.