JCDR - Register at Journal of Clinical and Diagnostic Research
Journal of Clinical and Diagnostic Research, ISSN - 0973 - 709X
Genetics Section DOI : 10.7860/JCDR/2021/46429.14990
Year : 2021 | Month : Jun | Volume : 15 | Issue : 06 Full Version Page : GC01 - GC04

Association of Interleukin-6-174G/C Gene Promoter Polymorphism with Rheumatoid Arthritis in Adolescence: A Case-control Study

Sergei Egorovich Khalchitsky1, Marina Vanikovna Sogoyan2, Alexei Nicolaevich Kozhevnikov3, Sergei Valentinovich Vissarionov4, Alexei Georgievich Baindurashvili5

1 Head, Department of Genetics, H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, Saint-Petersburg, Russia.
2 Research Associate, Department of Genetics, H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, Saint-Petersburg, Russia.
3 Paediatric Rheumatologist, Department of Trauma Effects and Rheumatoid Arthritis, H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, Saint-Petersburg, Russia.
4 Deputy Director for Research and Academic Affairs, Head, Department of Spinal Pathology and Neurosurgery, H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, Saint-Petersburg, Russia.
5 President, Department of Administration, H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, Saint-Petersburg, Russia.


NAME, ADDRESS, E-MAIL ID OF THE CORRESPONDING AUTHOR: Sergei Egorovich Khalchitsky, Parkovaya 64-68, Pushkin, Saint-Petersburg-196603, Russia.
E-mail: s_khalchitski@mail.ru
Abstract

Introduction

Rheumatoid Arthritis (RA) is a chronic autoimmune disease with unknown pathogenesis. The disease is multifactorial, however, the exact causes of the occurrence, as well as the targets of the autoimmune process, are unknown. Genes candidate for a predisposition to RA are cytokine genes. Of the cytokines, IL-6 is considered a key mediator of systemic and localised inflammation in RA.

Aim

The comparative analysis of the frequency distribution of alleles and genotypes IL-6-174G/C polymorphism in patients with RA and in the control group to determine the genotype most characteristic of this disease.

Materials and Methods

This case-control study was conducted in a group of 136 children with RA aged 14 to 18 years. In the control group, there were 143 practically healthy children of similar age without RA and orthopaedic pathology. Both groups were tested for polymorphism IL-6-174G/C using real-time Polymerase Chain Reaction (real-time PCR). The results were statistically processed using the Pearson’s Chi-square test.

Results

The distribution of 174G/C genotypes in the IL-6 gene was significantly different in patients with RA when compared with the control group. In patients with RA, heterozygous carriers of 174G/C (52.94%) prevailed, while in the control group the most numerous group were homozygous carriers of -174G/G (40.56%).

Conclusion

As a result of this study, it was possible to identify significant differences in distribution of genotypes 174G/C polymorphism in IL-6 gene between patients with RA and subjects of the control group. In relation to the European population of the Russian Federation, this polymorphism can serve as a diagnostic marker in the study of the pathogenesis of RA.

Keywords

Introduction

Rheumatoid Arthritis (RA) is a multisystem autoimmune disease that leads to inflammatory processes in the connective tissue, mainly affecting the small joints. The incidence of this disease worldwide is approximately 1% [1]. One of the variants of the clinical course of RA is juvenile RA, which occurs in childhood and adolescence, i.e., up to 16 years of age and, unlike adult RA, does not always become chronic. Despite the fact that the aetiology of RA is multifactorial, there are genetic markers that determine the predisposition to this disease. The HLA-DR loci are the genetic factor most associated with susceptibility to RA. Recent Genome-Wide Association Studies (GWAS) have identified 101 Single Nucleotide Polymorphism (SNPs) in total, showing the highest contribution of HLA-DRB1 gene to the development of RA [2-5]. Cytokine genes are other candidate genes for susceptibility to RA. Cytokines are a broad and loose category of small proteins (~5–20 kDa) important in cell signaling. Cytokines have been shown to be involved in signaling pathways and numerous proinflammatory cytokines such as Tumour Necrosis Factor alpha (TNF-α), Interleukin-1 (IL-1) and IL-6 are mediators that are involved in the inflammatory response and play an essential role in the pathogenesis of RA.

Of all cytokines, IL-6 is considered a key mediator of systemic and localized inflammation in RA. IL-6 is a strong inducer which can cause fever, anaemia, secondary amyloidosis, elevations in acute-phase proteins {C-Reactive Protein (CRP)}. The capability of IL-6 to influence B-cell differentiation can leads to the formation of rheumatoid factor and other autoantibodies. IL-6 aids osteoclast activation in joints and induces the release of matrix metalloproteinases, thus contributing to joint damage [6-9]. Several studies have found high levels of IL-6 in the synovial fluid of inflamed joints [10,11]. Changes in IL-6 levels are regulated by several factors [12], including variants of polymorphisms in promoter regions of the IL-6 gene.

The human IL-6 gene is located on chromosome 7p21. The most frequently studied polymorphism is the SNP-174C and -174G in the promoter region, which has been associated with transcription rates of IL-6. This polymorphism is characterised by a single nucleotide substitution of guanine (G) for cytosine (C) at position -174 in the promoter region [13,14], which leads to three variants of possible genotypes- GG, GC, and CC. This polymorphism is functionally important as it affects the rate of gene transcription and the plasma concentration of the IL-6 protein.

According to number of studies, 174G/C gene polymorphism plays an important role in the pathogenesis of various diseases, such as idiopathic sick sinus syndrome, chronic hepatitis B and C, chronic HCV infection and others [15-17]. Nevertheless, each of the allelic variants of this polymorphism can play a positive, negative or neutral role in the pathogenesis of the diseases under consideration. According to studies carried out in different countries and ethnic populations, there is a wide scatter of data on the distribution of IL-6 gene genotypes both in the populations of RA patients and in controls [18-24]. Taking into account the conflicting data available in the literature, authors carried out a molecular genetic study of 174G/C polymorphism among RA patients and in control group in the Russian population (St. Petersburg) and carried out a comparative analysis of the distribution of this polymorphism among RA patients and healthy people of the same age and compared with other populations of Northern and Southern Europe, Asia and Latin America to clarify the role of this gene in the pathogenesis of RA.

Materials and Methods

This case-control study was conducted from September 2018 to September 2020 on patients who underwent treatment for juvenile RA at the H. Turner National Medical Research Center for Children’s Orthopaedics and Trauma Surgery, St-Petersburg, Russia of the Ministry of Health of the Russian Federation. It was in accordance with the Declaration of Helsinki after being approved by the Local Ethics Committee of the Institute. Parents of the all underage patients with RA, participating in the study, gave informed consent for participation of children in the molecular genetic study of the IL-6 gene by PCR diagnostics.

Inclusion criteria for cases: A total of 136 cases of RA were included ranging in age from 14 to 18 years (50 males and 86 females), with confirmed diagnosis of juvenile RA according to the International League of Rheumatology Associations (ILAR) International Standard Definition [25], which must have occurred at least two months before enrollment at age <18 years. Arthritis of one or more joints with fever lasting for at least two weeks, swelling or limited mobility of the joint, morning stiffness and arthralgia during the day, joint pain or abnormal use of the joints were included.

Exclusion criteria for cases: Patients with active or recurrent bacterial, fungal or viral infections, including those with signs of HIV infection, hepatitis B and hepatitis C infection, patients with chronic diseases (kidney and liver disease).

Inclusion criteria for controls: The control group consisted of 143 apparently healthy males (61) and females (82) adolescents of the same age group, who were students of educational institutions of St. Petersburg, Russia without signs of RA and orthopaedic pathology according to the results of medical examinations.

Exclusion criteria for controls: It was the presence of RA and any orthopaedic diseases.

Genomic DNA Extraction and SNP Selection

2 mL Peripheral blood was taken from all patients and individuals of the control group for molecular genetic analysis and was frozen at -70oC. Genomic DNA was isolated from 100 μL blood using «ExtractDNA Blood» kit from Evrogen Company, Russia, according to the manufacturer’s instructions. Analysis of the polymorphic region of the IL-6 gene was carried out using the SNP-Screen reagent kit for IL-6-174G/C polymorphism manufactured by Syntol Company, Russia, by real-time PCR on a CFX96 Touch Real-Time PCR Detection System (Bio-Rad, USA).

Statistical Analysis

Statistical data processing was carried out using online medical statistics calculators. We used the Pearson’s Chi-square (χ2) test for a comparative analysis of the distribution of allele and genotype frequencies in groups of patients with RA, as well as in the control group. Differences in indicators between the groups of RA patients and people from the control group were considered statistically significant at p<0.05. The Odds Ratio (OR) and Hardy-Weinberg Equilibrium (HWE) methods were also used.

Results

Mean age of the patients was 16 years. Study results of IL-6-174G/C polymorphism in patients with RA and in the control group is presented in the [Table-Fig-1].

Distribution of genotypes and alleles in study groups.

IL-6-174G/C genotypes, allelesRA patients (n=136)%Control group (n=143)%OR (95% CI)
G/G4230.885840.560.655 (0.400-1.073)
G/C7252.945437.761.854 (1.151-2.988)
C/C2216.183121.680.697 (0.381-1.277)
G15657.3517059.44χ2=0.250p=0.6
C11642.6511640.56

p-value for genotypes=0.04


During the statistical processing of the results of the distribution of genotypes between RA patients and the control group, the following indicators were obtained: the value of the χ2=6,488 with the number of degrees of freedom equal to 2. Significance level p=0.04, that is, the relationship between the factorial and effective traits was statistically significant with an acceptable value of p<0.05. As can be seen from the table, as a result of PCR genotyping, in patients with RA (n=136), heterozygous carriers of 174G/C (52.94%) prevailed, while in the control group (n=143) the group of homozygous carriers of -174 G/G (40.56%) was the most numerous.

At the same time, no significant difference was observed in the ratio of the number of G and C alleles between RA patients and the control group. Here, the value of the χ2=0.250 and the significance level was p=0.617, which is much higher than p<0.05.

When calculating the HWE index, it can be seen that in patients and in controls there were significant deviations from the expected equilibrium in opposite directions. Thus, in the case of heterozygous genotype 174 (G/C), the deviation from the expected value of HWE in RA patients occurred by 4.04% upwards, and in the control group by 10.44% downwards [Table/Fig-2].

Hardy-Weinberg test (HWE) for the RA group and control group.

IL-6-174G/C genotypesRA patients, % (n=136)HWEControl group, % (n=143)HWE
G/G30.8832.940.5635.3
G/C52.9448.937.7648.2
C/C16.1818.221.6816.5

Discussion

Our study of the IL-6-174G/C polymorphism in patients with juvenile RA showed the following regularities: in the group of RA patients, heterozygous (G/C) carriage of this polymorphism predominates (52.94%). At the same time, in the control group of healthy volunteers, the major genotype G/G was most represented (40.56%), as in many other populations.

This distribution of genotypes in patients with RA in North-West Russia is similar in nature to patients in the North European populations [21,23]. The most numerous group was the group of heterozygous GC carriers [Table-Fig-3].

Distribution of genotypes and alleles in the group of patients with RA and among residents of northern Europe.

IL-6-174G/C genotypesRA, % Great Britain [21]n=383Control group, % Great Britain [21] n=422RA, % Poland [23] n=98Control group, % Poland [23] n=105RA, % Russia, present studyControl group, % Russia, present study
G/G35.235.126.523.830.8840.56
G/C46.05054.155.252.9437.76
C/C18.814.919.421.016.1821.68
p=0.29p=0.9p=0.04

p<0.05 is considered to be significant


In the countries of southern Europe, the proportion of heterozygous GC carriers decreases and the GG genotype becomes dominant [Table-Fig-4] [20,22].

Distribution of genotypes and alleles in the group of patients with RA and among residents of southern Europe.

IL-6-174G/C genotypesRA, % North Macedonia, Trajkov D et al., 2009 [20] n=85Control group, % North Macedonia, Trajkov D et al, 2009 [20] n=301RA, % Spain, Pascual M et al., 2000 [22] n=163Control group, % Spain, Pascual M et al., 2000 [22] n=157RA, % Russia, present studyControl group, % Russia, present study
G/G45.247.846.046.530.8840.56
G/C35.743.944.242.052.9437.76
C/C19.18.39.811.516.1821.68
p=0.068p=0.864p=0.04

p<0.05 is considered to be significant


In Latin American and Asian populations, the GG genotype is even more prevalent, and the C allele has an extremely low representation [Table-Fig-5] [19,24].

Distribution of genotypes and alleles in the group of patients with RA and among residents of Asia and Latin America.

IL-6-174G/C genotypesRA, % Mexico, Zavaleta-Muñiz SA et al., 2013 [24] n=137Control group, % Mexico, Zavaleta-Muñiz SA et al., 2013 [24] n=102RA, %, China, Li F et al., 2014 [19] n=256Control group, % China, Li F et al., 2014 [19]n=331RA, % Russia, present studyControl group, % Russia, present study
G/G77.478.496.599.430.8840.56
G/C21.919.62.70.352.9437.76
C/C0.720.80.316.1821.68
p=0,651p=0.049p=0.04

p<0.05 is considered to be significant.


At the same time, it should be emphasised that in the Russian Federation, due to the fact that it is a multinational country, there are significant differences in the distribution of genotypes and alleles in regions where the majority of the population are either Caucasian (Slavic) or Mongoloid (Asian) population. Nevertheless, a more important indicator, in our opinion, is the assessment of the significance of differences in the distribution of genotypes between patients and the control group. However, there is also mixed evidence in various studies. So, in RA groups in Egypt, China, Macedonia [18-20], these differences were significant, and in RA groups in UK, Spain, Poland [21-23] there were no significant differences.

There can be several explanations for these facts. First, it is logical to assume that with such a wide scatter of data, the 174G/C polymorphism does not play a significant role in the aetiology and pathogenesis of the RA, but there are other, more significant factors, possibly epigenetic methylation factors [26]. Second, where the difference between RA patients and the control group was significant, the situation changes, and obviously, this factor takes a more important place, while the importance of other genetic and epigenetic determinants decreases. Some of the increased risk is due to environmental factors such as acute viral or bacterial infections, prolonged hypothermia and other stochastic influences. Smoking is the best defined risk and interacts with HLA-DR to amplify the risk of developing RA. Although the search for genetic variance continues, we have likely reached the point of diminishing returns. For example, the disease concordance rate between identical twins is surprisingly low. If one twin develops RA the chance that the second twin will develop the disease is only 12% to 15% [27]. In this regard, in order to understand the determining genetic and epigenetic factors of the aetiology and pathogenesis of RA, it is necessary to search for factors that would have an equally high penetrance for worldwide populations. In this regard, we continue to look for such factors, expanding our research and investigating genes for other proinflammatory cytokines, as well as genetic and epigenetic mechanisms leading to increased expression of HLA-antigens.

Limitation(s)

Due to limited possibilities, the authors did not conduct a study of the IL-6 level in the blood and tissues, which could supplement the information on the ratio of IL-6 genotypes with the concentration distribution of IL-6 in tissues.

Conclusion(s)

This study revealed significant differences in the distribution of genotypes 174G/C of the IL-6 gene between patients with RA and a control group in the Russian population (St. Petersburg). This indicates that this polymorphism can serve as a diagnostic marker when studying the pathogenesis of RA in the European population of the Russian Federation. The authors are currently continuing the research on other proinflammatory cytokines along with investigating other possible genetic, pharmacogenetic and epigenetic markers of the aetiopathogenesis of RA. In particular, genes associated with apoptosis and genes of the folate cycle are being investigated.

The main recommendations for future studies of RA lie in the plane of epigenetic modifications in inflammatory processes that leads to the triggering of autoimmune switching in a predisposed person. The presented work makes a certain contribution to the accumulation of population information on the role of proinflammatory cytokines in the pathogenesis of various diseases.

p-value for genotypes=0.04p<0.05 is considered to be significantp<0.05 is considered to be significantp<0.05 is considered to be significant.

References

[1]Gibofsky A, Overview of epidemiology, pathophysiology, and diagnosis of rheumatoid arthritis Am J Manag Care 2012 18(13 Suppl):S295-302.  [Google Scholar]

[2]Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis Nat Genet 2012 44:1336-40.10.1038/ng.246223143596  [Google Scholar]  [CrossRef]  [PubMed]

[3]Okada Y, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population Nat Genet 2012 44:511-16.10.1038/ng.223122446963  [Google Scholar]  [CrossRef]  [PubMed]

[4]Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, Genetics of rheumatoid arthritis contributes to biology and drug discovery Nature 2014 506:376-81.10.1038/nature1287324390342  [Google Scholar]  [CrossRef]  [PubMed]

[5]Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci Nat Genet 2010 42:508-14.10.1038/ng.58220453842  [Google Scholar]  [CrossRef]  [PubMed]

[6]Srirangan S, Choy EH, The Role of Interleukin 6 in the Pathophysiology of Rheumatoid Arthritis Ther Adv Musculoskelet Dis 2010 2(5):247-56.10.1177/1759720X1037837222870451  [Google Scholar]  [CrossRef]  [PubMed]

[7]Alunno A, Carubbi F, Giacomelli R, Gerli R, Cytokines in the pathogenesis of rheumatoid arthritis: New players and therapeutic targets BMC Rheumatol 2017 1:310.1186/s41927-017-0001-830886947  [Google Scholar]  [CrossRef]  [PubMed]

[8]Narazaki M, Kishimoto T, The two-faced cytokine IL-6 in host defense and diseases Int J Mol Sci 2018 19(11):352810.3390/ijms1911352830423923  [Google Scholar]  [CrossRef]  [PubMed]

[9]Kany S, Vollrath JT, Relja B, Cytokines in inflammatory disease Int J Mol Sci 2019 20(23):600810.3390/ijms2023600831795299  [Google Scholar]  [CrossRef]  [PubMed]

[10]Hirano T, Matsuda T, Turner M, Miyasaka N, Buchan G, Tang B, Excessive production of interleukin 6/B cell stimulatory factor-2 in rheumatoid arthritis European Journal of Immunology 1988 18(11):1797-801.10.1002/eji.18301811222462501  [Google Scholar]  [CrossRef]  [PubMed]

[11]Nishimoto N, Kishimoto T, Interleukin 6: From bench to bedside Nature Clinical Practice Rheumatology 2006 2(11):619-26.10.1038/ncprheum033817075601  [Google Scholar]  [CrossRef]  [PubMed]

[12]Luo Y, Zheng SG, Hall of fame among pro-inflammatory cytokines: Interleukin-6 gene and its transcriptional regulation mechanisms Front Immunol 2016 7:60410.3389/fimmu.2016.0060428066415  [Google Scholar]  [CrossRef]  [PubMed]

[13]Fishman D, Faulds G, Jeffey R, Mohamed-Ali V, Yudkin JS, Humphries S, The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with systemic-onset juvenile chronic arthritis Journal of Clinical Investigation 1998 102(7):1369-76.10.1172/JCI26299769329  [Google Scholar]  [CrossRef]  [PubMed]

[14]Terry CF, Loukaci V, Green FR, Cooperative influence of genetic polymorphisms on interleukin 6 transcriptional regulation Journal of Biological Chemistry 2000 275(24):18138-44.10.1074/jbc.M00037920010747905  [Google Scholar]  [CrossRef]  [PubMed]

[15]Nikulina SYu, Marilovtseva OV, Chernova AA, Tretyakova SS, Nikulin DA, Maksimov VN, The role of interleukin-6 gene in development of idiopathic sick sinus syndrome Russian Journal of Cardiology 2016 138(10):32-36.(In Russ.)10.15829/1560-4071-2016-10-32-36  [Google Scholar]  [CrossRef]

[16]Romanova SV, Vidmanova TA, Zhukova EA, Ermolina EV, Tolkacheva NI, Mayanskaya IV, The role of genetic polymorphisms of gene IL-1β (-31 T/C) and IL-6 (-174 G/C) in the course of chronic hepatitis B and C in children Medical Immunology 2014 15(6):535-42.(In Russ.)10.15789/1563-0625-2013-6-535-542  [Google Scholar]  [CrossRef]

[17]Semyonova NA, Ryazantseva NV, Novitsky VV, Bychkov VA, Chechina OY, Role of IL6 -174C/G gene polymorphism in development of chronic HCV infection Bulletin of Siberian Medicine 2010 5:93-97.(In Russ.)10.20538/1682-0363-2010-5-93-97  [Google Scholar]  [CrossRef]

[18]Amr K, El-Awady R, Raslan H, Assessment of the -174 G/C (rs1800795) and -572 G/C (rs1800796) Interleukin 6 gene polymorphisms in Egyptian patients with rheumatoid arthritis Open Access Maced J Med Sci 2016 4(4):574-77.10.3889/oamjms.2016.11028028393  [Google Scholar]  [CrossRef]  [PubMed]

[19]Li F, Xu J, Zheng J, Sokolove J, Zhu K, Zhang Y, Association between interleukin-6 gene polymorphisms and rheumatoid arthritis in Chinese Han population: A case-control study and a meta-analysis Sci Rep 2014 4:571410.1038/srep0571425030201  [Google Scholar]  [CrossRef]  [PubMed]

[20]Trajkov D, Mishevska-Perchinkova S, Karadzova-Stojanoska A, Petlichkovski A, Strezova A, Spiroski M, Association of 22 cytokine gene polymorphisms with rheumatoid arthritis in population of ethnic Macedonians Clin Rheumatol 2009 28:1291-300.10.1007/s10067-009-1238-419662331  [Google Scholar]  [CrossRef]  [PubMed]

[21]Panoulas VF, Stavropoulos-Kalinoglou A, Metsios GS, Smith JP, Milionis HJ, Douglas KM, Association of interleukin-6 (IL-6) -174G/C gene polymorphism with cardiovascular disease in patients with rheumatoid arthritis: The role of obesity and smoking Atherosclerosis 2009 204:178-83.10.1016/j.atherosclerosis.2008.08.03618848327  [Google Scholar]  [CrossRef]  [PubMed]

[22]Pascual M, Nieto A, Matarán L, Balsa A, Pascual-Salcedo D, Martín J, IL-6 promoter polymorphisms in rheumatoid arthritis Genes Immun 2000 1(5):338-40.10.1038/sj.gene.636367711196696  [Google Scholar]  [CrossRef]  [PubMed]

[23]Pawlik A, Wrzesniewska J, Florczak M, Gawronska-Szklarz B, Herczynska M, IL-6 promoter polymorphism in patients with rheumatoid arthritis Scand J Rheumatol 2005 34:109-13.10.1080/0300974051002637316095005  [Google Scholar]  [CrossRef]  [PubMed]

[24]Zavaleta-Muñiz SA, Martín-Márquez BT, Gonzalez-Lopez L, Gonzalez-Montoya NG, Díaz-Toscano ML, Ponce-Guarneros JM, The -174 G/C and -572 G/C Interleukin 6 promoter gene polymorphisms in mexican patients with rheumatoid arthritis: A case-control study Clin Dev Immunol 2013 2013:95908410.1155/2013/95908424223608  [Google Scholar]  [CrossRef]  [PubMed]

[25]Petty RE, Southwood TR, Manners P, Baum J, Glass DN, Goldenberg J, International League of Associations for RheumatologyInternational League of Associations for Rheumatology classification of juvenile idiopathic arthritis: Second revision, Edmonton, 2001 The Journal of Rheumatology 2004 31(2):390-92.  [Google Scholar]

[26]Zhu H, Wu L, Mo X, Lu X, Tang H, Zhu XW, Rheumatoid arthritis–associated DNA methylation sites in peripheral blood mononuclear cells Annals of the Rheumatic Diseases 2019 78:36-42.10.1136/annrheumdis-2018-21397030297333  [Google Scholar]  [CrossRef]  [PubMed]

[27]Aho K, Koskenvuo M, Tuominen J, Kaprio J, Occurrence of rheumatoid arthritis in a nationwide series of twins J Rheumatol 1986 13:899-902.  [Google Scholar]