OBJECTIVE To study the association between rheumatoid arthritis (RA) and HLA and tumour necrosis factor (TNF) polymorphism in Peruvian mestizo patients in comparison with ethnically similar controls.
METHODS Seventy nine patients with RA and 65 ethnically matched healthy controls were genotyped for HLA-DRB1, HLA-DQA1, HLA-DQB1, and TNFα and TNFβ alleles using PCR amplification. Clinical severity was assessed as mild, moderate, or severe in 35 of the patients.
RESULTS TNFα6 showed the strongest association with disease susceptibility. The TNFα6 allele was more common in patients than in controls (p<0.0076) and the proportion of patients with at least one copy of this allele was greater (p<0.015, relative risk 2.35). Among the HLA-DRB1* alleles with the shared epitope sequence, only the DRB1*1402 allele was significantly increased in patients compared with controls (p<0.0311), as was the proportion of patients with at least one copy of this allele (p<0.0232, relative risk 2.74). In contrast, the overall frequency of alleles with the shared epitope was not different in patients and controls. The haplotype HLA-DRB1*1402-DQB1*0301-DQA1*0401 was significantly more common in patients. TNFα6 was more common in patients whether or not they had this haplotype. None of the 11 patients lacking the TNFα6 allele had severe disease.
CONCLUSIONS This study shows for the first time that TNF gene polymorphism is associated with susceptibility to RA in a non-white population. TNFα6 and HLA-DRB1*1402 independently conferred significantly increased risk in Peruvian mestizo patients.
- rheumatoid arthritis
- tumour necrosis factor
Statistics from Altmetric.com
Rheumatoid arthritis (RA) is a chronic inflammatory disease with worldwide distribution. Its manifestations and progression rate show pronounced individual variation both within and between cohorts of different geographical origin.1 ,2 Concordance in monozygotic twins is 10–30%, indicating that both genetic and environmental factors are important pathogenetically.3Several studies in different populations have shown that only a few of the 15 or more DRB1* alleles are associated with susceptibility to RA and these contain the amino acid sequence QRRAA, QKRAA, or RRRAA in position 70–74 on the third hypervariable region of the DRB1*β chain, the so called shared epitope.4 The presence of a shared epitope was found in some studies to be correlated with a more destructive form of RA, particularly if it was present on both alleles.5 ,14 More recently, polymorphism has been described in the genes for tumour necrosis factors (TNF) located on the short arm of chromosome 6 within the major histocompatibility complex.6 Disease correlations have been described with systemic lupus erythematosus (SLE), but these have been secondary mainly to linkage disequilibrium due to the presence of alleles with previously established disease associations.7 The best studied polymorphisms relate to the TNFα and TNFβ microsatellites located 3.5 kB upstream of the TNF gene, and 13 and eight alleles, respectively, have been identified.8 Interestingly, TNF polymorphism has been linked to the in vitro response of macrophages to lipopolysaccharide9 ,10 and therefore it could be of potential importance in the regulation of the host reaction in inflammation.
RA is an important cause of chronic and disabling polyarthritis in Peru, but very little is known about HLA and TNF polymorphisms in this population. We report such data in unselected Hispano-Indian (mestizo) patients with RA seen in two hospitals in Lima and in ethnically matched local controls.
Subjects and methods
PATIENTS AND CONTROLS
Seventy nine patients with RA and 65 ethnically matched healthy controls were studied (table 1). All patients fulfilled the ACR criteria for RA. All patient samples were available for DQA1* typing but, for technical reasons, accurate DQB1* typing was possible in only 76 patients and DRB1* typing in 78 patients. All 79 samples were tested for TNFα but only 78 for TNFβ alleles. In 35 patients the disease was graded as mild, moderate, or severe based on functional status, extent and severity of joint involvement, and presence or absence of extra-articular manifestations. Rheumatoid factor was determined in the routine laboratory by the sheep cell agglutination test.
Genomic DNA was extracted from whole blood using a standard phenol chloroform extraction procedure. One μg DNA was amplified by the polymerase chain reaction (PCR) using specific primers for DRB1* and DQB1* loci. The sequences of the primers and oligonucleotide probes used for the genetic typing have been described elsewhere, as well as the conditions for amplification and hybridisation.11 Five μl of the PCR product was denatured, spotted onto nylon membranes (Pall, CA, USA) using a dot blot apparatus (Biorad, CA, USA), and hybridised with either horseradish peroxidase (HRP) or biotin labelled oligonucleotide probes. Negative controls containing no DNA and positive controls with DNA from homozygous cell lines were included on each membrane. The DQA1* typing was performed with a commercially available reverse dot blot assay (Amplitype, Cetus, CA, USA).
A non-radioactive method was used for typing TNFα and TNFβ microsatellites.12 The method included two consecutive PCR amplifications, in the second of which biotinylated primers were used. In the first PCR we obtained a fragment of 200 bp. In the second amplification 1 μl of 1:100 dilution of the product obtained from the first PCR was used for the TNFα and TNFβ specific amplification. The PCR mixture contained 1 μM of each primer, 200 μM of each dNTP, 1.5 mM of MgCl2, and 2.5 units of Taq polymerase. The total volume was 50 μl. This mixture was subjected to 30 cycles of amplification (94°C for 1 minute, 68°C for 1 minute, and 72°C for 1 minute). To the final specific PCR product, one volume of loading buffer (98% formamide) was added. The samples were heated to 80°C for 2 minutes and loaded on a routine sequencing gel at 90 W for 3 hours in 1 × TBE (0.1 M Tris base, 0.009 M boric acid, 1 mM EDTA). After electrophoresis the PCR products were blotted onto a positively charged nylon membrane for 3 hours (Amersham-Pharmacia, Sweden). The membrane was air dried and saturated with 5% SDS and 0.9% NaCl and then left for 20 minutes in a solution containing 0.5% SDS, 0.9% NaCl, and 1 μg/ml streptavidin-HRP. The membrane was washed twice in 0.5% SDS, 0.9% NaCl for 5 minutes. Detection of TNFα and TNFβ microsatellites was performed by chemiluminescence (ECL, Amersham-Pharmacia, Sweden) according to the manufacturer's protocol.
For unambiguous typing, control DNA—that is, previously sequenced DNA—was included in all typing experiments.
Statistical analysis was performed using the χ2test. The level of significance was set at 0.05. Relative risk (RR) was estimated according to Wolf with Haldane's correction.13
Table 2 shows the results of DRB1* typing. DRB1*1402 was the only allele with the shared epitope sequence which was significantly more common in patients than in controls (17.9%v 8.5%, p<0.0311). The proportion of patients with at least one copy of this allele was similarly significantly increased compared with controls (33.3%v 15.4%, p<0.0232, relative risk 2.74). In contrast, the overall frequency of alleles with the shared epitope was not different between patients and controls. The proportion of subjects with none, one, or two alleles with the shared epitope was similar in patients and controls (fig 1). The DRB1*11 containing alleles were more common in controls than in patients (7.7% v0.6%, p<0.0045).
Data on clinical severity (assessed as mild, moderate, or severe) were available in 35 patients. Those with none, one, or two alleles with the shared epitope were equally represented in all three categories of disease severity (data not shown).
Tables 5 and 6 show the results of TNFα and TNFβ typing, respectively. The TNFα6 allele occurred more frequently in patients than in controls (36.7% v 21.5%, p<0.0076), and the proportion of patients with at least one copy of this allele was greater (59.5% v 38.5%, p<0.015, relative risk 2.35). There were more homozygous and heterozygous TNFα6 individuals in the patient group than in the control group (fig 2). TNF6 seems to be more closely related to RA than DRB1*1402. The number of patients with the TNFα10 allele was decreased compared with the controls (44.3%v 66.2%, p<0.014).
The overall TNFβ5 allele frequency was significantly increased among the patients (53.8% v 37.7%, p<0.009) and a similar trend was found when comparing phenotypic frequencies (75.6%v 61.5%, p<0.10). All but one of the TNFα6 positive patients were also positive for TNFβ5. Both TNFβ1 and TNFβ4 phenotype frequencies (10.3% v24.6%, p<0.019, and 59.0% v 78.5%, p<0.021) were decreased in patients.
There was no correlation between the occurrence of rheumatoid factor and TNF status. Furthermore, there was no significant correlation between disease severity and the presence of TNFα6. However, five of 24 (20.8%) patients with TNFα6 were classified as severe compared with none of 11 TNFα6 negative patients (NS).
The extended haplotype DRB1*1402-DQB1*0301-DQA1*0401 (based on known linkage disequilibrium) was found in 26 of 78 patients (33.3%) compared with 10 of 65 controls (15.4%; p<0.023). TNFα6 was more common among patients, whether or not they had this haplotype (16 of 26 (61.5%) v 30 of 52 (57.7%)). An association between DRB1*1402 and TNFα6 in this population is likely since about half of DRB1*1402 individuals also have TNFα6β5, with no differences between RA patients and controls. Moreover, an association between these two alleles is described in the cell line ASHI no EK 9054. However, since we did not perform a family study, we cannot formally confirm the linkage. To test for independent associations, TNFα6 frequencies were also analysed in DRB1*1402 negative subjects and a statistically significant increase was detected in patients with RA compared with controls (57.7%v 36.4%; χ2=4.065; p<0.044).
HLA-DR genetics has been studied in patients with RA for two decades, establishing susceptibility associated with a structure in the third hypervariable region of the DRB1 chain in position 70–74. Susceptibility is conferred by the sequence QRRAA, QKRAA, or RRRAA, called the shared epitope.4 Of particular interest is the fact that the correlation with DRB1* genes varies between ethnic groups.14 The most common association found in white populations and Eskimos is with DRB1*0401 and 0404, in Asians with DRB1*0405, in Mediterranean populations with DRB1*0101 and 1001, and in North American Indians with DRB1*1402.14-17 Our study shows that DRB1*1402 is also over-represented in Peruvian mestizos, a population of descendants of Hispanics and American Indians. Its presence conferred susceptibility with a relative risk of 2.74. This significant association was shown, despite the limited number of individuals studied. The modest level of the association could be explained in part by the relatively high prevalence of DRB1*1402 in the general mestizo population. Likewise, the overall prevalence of alleles with the shared epitope in the Peruvian population was as high as 58.5% compared with 66.7% in patients with RA. The presence of either a single or two shared epitopes is therefore not significantly predictive of RA in this population (fig 1). Some, but not all, studies have found a positive correlation between disease progression and the dose of the shared epitope.14 The limited clinical data in our study do not allow definite conclusions in this regard.
It has been proposed that only one third of the genetic predisposition for RA is related to HLA-DRB1*.14 ,18 The remaining candidate genes are still largely unknown. DQB1* has been suggested in some studies, but this was later found to be secondary to linkage disequilibrium with HLA-DRB1*.19 Our findings support the previous observations. Other gene polymorphisms that did not vary between patients with RA and controls independently of HLA-DRB1* include TAP1 and TAP2.20 ,21
The Gm system has been studied in white populations and, although anti-Gm antibodies showed strong correlation with destructive disease in one study,22 no association between the genes of the Gm system and RA was detected. Conflicting results have been reported regarding HLA-DM.23 ,24 The class III C4 null genes have been found to be associated with susceptibility to RA in Asians and with radiographic progression in a Finnish study.25 ,26
In our study the strongest association was with the TNFα6 allele, which was present in almost 60% of the patients compared with less than 40% of the controls. All but one of the TNFα6 positive patients were also TNFβ5 positive. Figure 2 shows that there is a significant difference in the distribution of TNFα6 alleles between patients and controls. TNF polymorphism has recently attracted increasing interest in relation to RA. Early work suggested no correlation with gene polymorphisms in TNFβ or TNFα,27 ,28 but an Irish study reported results from multicase families where a TNFα6β5 containing haplotype was present in 35.3% of affected individuals but in only 20.5% of non-affected individuals.29 In a British study the TNFα6 allele was found to be increased with an odds ratio of 2.5.19 Importantly, as in our study, these studies indicate an influence that is independent of HLA-DRB1* status. It has been reported that the association is more pronounced in female patients.30 ,31 A Spanish study found that the presence of TNFα6β5 increased susceptibility to RA, regardless of the presence or absence of the shared epitope.32 One recent study described an association with disease progression and guanine to adenine transition at position –238 of the promoter region of the TNFα gene.33 ,34 Robust evidence, gathered in several populations and looking at different regions, is therefore emerging which links TNF gene polymorphism with susceptibility to RA, and perhaps with disease progression. Our study shows for the first time an association in a Hispano-Indian population with an HLA polymorphism, HLA-DRB1*1402, that differs from that in white patients with RA. Our results also show that the TNFα6 allele is associated with RA independently of HLA-DRB1*1402.
As a result of the strong linkage disequilibrium between TNFα and TNFβ alleles, an increase in TNFβ5 (which is in linkage with TNFα6) and a decrease in TNFβ4 (which is in linkage with TNFα10) was an expected result in our study.
An intriguing finding was the low prevalence of the TNFα10 allele. Whereas no information is available regarding altered TNF production in individuals with TNFα6, high TNF production has been connected with TNFα10 haplotypes.10 Complex interactions between HLA and TNF gene polymorphism have been suggested in a recent report.35 Future studies will address the question if and how this polymorphism affects disease manifestations, and explore the correlation between TNFα6 and disease severity.
Drs Castro and Acevedo contributed equally to this work. The skilful assistance of Mrs Karin Sjögren is gratefully acknowledged.
↵* Present address: Istituto Nationale Neurologico “Carlo Besta”, Milano, Italy
This study was supported by The Swedish Medical Research Council (Project No. K00-06265-18C), Alfred Österlund's Foundation, and Greta and Johan Kock's Foundations.
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.