Journal List > J Korean Orthop Assoc > v.45(5) > 1013021

Kim, Han, Kwak, Chun, Kim, and Lee: Association of GnRH1 Polymorphisms with Rheumatoid Arthritis in a Korean Female

Abstract

Purpose

Rheumatoid arthritis (RA) is a common, chronic inflammatory arthritis that develops most often in women. Gonadal hormones may account for the sexual dimorphism in the immune response and for the greater incidence of autoimmune disease in females. Gonadotrophin-releasing hormone (GnRH), one of the gonadal hormones, plays an important role in immune system modulation. This study examined the effects of single nucleotide polymorphisms (SNP) in GnRH on gender differences in the pathophysiology of RA.

Materials and Methods

The presence of SNPs rs2659590, rs2321248, rs6186, rs6185, and rs2321049 in the human GnRH1 gene was confirmed in Korean RA patients by Taqman® SNP genotyping assays. A total of 153 unrelated female, Korean RA patients and 96 female Korean controls participated.

Results

There were no significant associations between GnRH1 polymorphisms and RA. However, we found that the rs2659590, rs6185 and rs2321248 polymorphism might be associated with a susceptibility to aberrantly high erythrocyte sedimentation rates in female RA patients.

Conclusion

Additional studies, with a larger number of patients and in different populations will be required to assess whether GnRH1 polymorphisms and these haplotypes could be used as susceptibility or resistance markers in RA. To our knowledge, this study is the first to analyze associations between SNPs of GnRH1 and RA.

Introduction

Rheumatoid arthritis (RA) is a systemic inflammatory autoimmune disease associated with swelling and pain in multiple joints.1) Despite intensive research, the cause of RA remains obscure. Proposals of the initiating event include an infectious agent or other environmental exposure but genetic, hormonal, and reproductive factors may all contribute to developing RA. RA affects 1% to 2% of adults and, like most autoimmune diseases, develops most often in women, with a female/male ratio of 3:1.2) Previous report reviewed that hormonal involvement is the higher incidence in female RA patients, because disease symptoms lessen during pregnancy and exacerbate in the postpartal period.3) Evidence for endocrine involvement in RA tissue and cells include presence of androgen and estrogen receptors, high concentrations of biologically active steroids, key enzymes of steroid metabolism, and significant changes of estrogen to androgen ratio. These data suggest that individual immune cells, including synovial macrophages, may behave as steroid-sensitive cells.4) Previous studies showed that the immunomodulatory effects of androgens and estrogens may be mediated indirectly.5) One of the best-characterized targets for feedback actions of gonadal steroids is the hypothalamic hormone gonadotrophin-releasing hormone (GnRH).
GnRH is a peptide synthesized and released by the hypothalamus that is responsible for the release of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) from the anterior pituitary. GnRH1 gene located on chromosome 8 p21-11.2 synthesizes precursor GnRH1 protein of 92 amino acids, which is processed to GnRH1 decapeptide in mammals.6) The gonadotroph cells in the anterior pituitary are bathed with pulsatile GnRH which then stimulates its own receptor (GnRHR) to cause secretion of LH and FSH. In females, the acute rise of LH during the LH surge, brought about by the process of GnRH-self priming, triggers ovulation.7) GnRH itself has been shown to be immunostimulatory. Immune cells produce bioactive GnRH and express receptors for GnRH.8-11) Previous studies showed that GnRH exerts stimulatory influences on expression of the interleukin-2 receptor, on B and T lymphocyte proliferation, and on serum IgG levels.12,13) Based on those data, it was hypothesized that differences in responsiveness to GnRH might be due to gender differences in expression of GnRH receptor or due to gender differences in expression of G proteins.13)
We examined five single nucleotide polymorphisms and haplotypes in the GnRH1 regions and their association with susceptibility and severity with RA in female Korean patients.

Materials and Methods

1. Subjects

A total of 153 unrelated Korean female RA patients were participated (Table 1). Each patient was diagnosed by a rheumatologist and diagnosed according to diagnostic criteria, such as erythrocyte sedimentation rate (ESR), C-reactive protein (CRP) and bone erosion. The control subjects (96 female) were recruited after they had been designated as no clinical evidence of RA in a general health check-up program. Written informed consent was obtained from all subjects. This study was approved by the Ethics Committee of Kyung Hee University, Faculty of Medicine.

2. Genotyping

Genomic DNA was extracted from whole blood cells using a Nucleo Spin® Blood kit (Macherey-Nagel, Easton, PA, USA). Five polymorphisms of a GnRH1 gene were selected (Fig. 1). Of the 5 SNPs known to be exist in GnRH1 (http://www.ensembl.org; www.ncbi.nlm.nih.gov/SNP), two SNPs (rs6186, rs6185) were existed in coding region. The corresponding probes were ordered at the company Applied Biosystems, https://products.appliedbiosystems.com (Table 2). Allelic discrimination using fluorogenic probes (5' nuclease assay, Taqman®) was chosen for genotyping on the ABI PRISM 7700 Sequence Detection Systems (Applied Biosystems, Foster City, CA, USA) and consisted of the use of allele-specific fluorogenic probes.14) Twenty nanograms of DNA were amplified by AmpliTaq Gold DNA polymerase which was included in Taqman® Universal Master Mix (Applied Biosystems). The PCR reactions were as follows: one step of 10 min at 95℃ followed by 40 cycles of two-step PCR with denaturation at 92℃ for 15 s and annealing and extension at 60℃ for 1 min. Five percent of all samples were genotyped again for quality control.

3. Statistics

For theses association study, Hardy-Weinberg equilibrium (HWE) for all SNPs was assessed using SNPstats (http://bioinfo.iconcologia.net/index.php).15) A linkage disequilibrium (LD) block of polymorphisms was tested using Haploview (version 4.0).16) The haplotypes and their frequencies were calculated by the EM algorithm.17) Multiple logistic regression models were calculated for the odds ratio (OR), 95% confidence interval (CI) and corresponding p-values, controlling for age as covariables. For all statistical tests, A p-value of p<0.05 was taken as significant.

Results

1. Clinical characteristics of RA in the female Korean

Clinical characteristics of RA patients and control subjects are shown in Table 2. One hundred fifty three female RA patients and 96 female healthy controls were recruited for this study. The mean age in RA patients was 49.3±5.4 years and control subjects 68.9±5.4. The mean disease duration was 4.68±3.53 years. ESR in RA patients was 37.77±26.52 (mm/hr) and CRP was 2.23±6.18 (mg/dL). Seventy three of a total of 153 RA patients exhibited bone erosion (47%).

2. Association between GnRH1 polymorphisms and RA

To assess whether polymorphisms of GnRH1 have any correlation with the risk of rheumatoid arthritis, genotype frequencies of 5 SNPs (rs2659590, rs2321248, rs6186, rs6185 and rs2321049) were examined in 153 rheumatoid arthritis and 96 control subjects, and multiple logistic regression analysis with adjustment for age was done.
Genotype distributions of 5 polymorphisms in this study were in Hardy-Weinberg equilibrium (p>0.05 data not shown). Genotype distributions and allele frequency of 5 SNPs in RA and control subjects are revealed in Table 2 (Table 3). Among 5 SNPs, rs6186 was monomorphic SNP. There was no significant association with RA in genotype distribution and allele frequency between RA female patients and controls (Table 3).

3. Associations between GnRH1 Polymorphisms and parameters of disease activity

As phenotyping stratification, association between polymorphisms and biomarker, ESR and bone erosion was compared as parameters of disease activity.
ESR between the ESR 30< and ≥30 was compared (Table 4). Of the five SNPs, 2 SNPs were significantly associated with the risk of RA: rs2659590, OR=4.37, 95% CI=1.21-15.72, p=0.011 in the recessive model; rs6185, OR=1.79, 95% CI=1.09-2.99, p=0.019 in the codominant model and OR=2.29, 95% CI=1.14-4.59, p=0.02 in the dominant model; respectively (Table 5). In the case of rs2321248, there was significant association with ESR level in the G/A genotype; OR=6.67, 95% CI= 1.47-30.16, p=0.0028. Existence of bone erosion was also compared in RA patients. There were no significance in the genotype distribution and allele frequency between bone erosion (Table 6).

4. Associations between GnRH1 Haplotypes and ESR parameter in disease activity

Three among 5 SNPs formed LD block. Haplotypes (frequency>0.01) in the LD block were analyzed using HapAnalyzer software 1.0.18)
Four haplotype were shown. On the haplotype association test, haplotypes, CAG and TGG, showed a significant association between the ESR<30 and ESR 30 RA patients, respectively (Table 7).

Discussion

RA is a common autoimmune disease with a complex etiology in which genetic and environmental factors contribute to disease. The genetic component of RA is largely undefined and, up until very recently, there were only two reproducible associations. The strongest of these associations is of genes within the HLA region, particularly the HLA-DRB1 gene.19) A second, more modest, association identified has been of the protein tyrosine phosphatase non-receptor 22 (PTPN22) gene.20,21) In this study, we demonstrated that whether GnRH1 gene polymorphisms are associated with RA. There was no association between RA susceptibility and GnRH1 in female population.
Altogether the genetics of RA susceptibility has made enormous progress during the past years. In contrast, the field of the genetics of RA severity is relatively unexplored. Radiological joint damage is conceived as the most objective measure to determine the severity of RA. The number of studies investigating the relation between joint destruction and genetic variants is limited. Some SNPs are observed to associate with the disease outcome, but none of these associations are convincingly replicated in independent cohorts. Moreover, although a recent twin study indicated that genetic factors play a role in determining the severity of RA,22) the heritability of the level of joint destruction (the variance of joint destruction that can be ascribed to genetic factors) is still unknown.
Joint damage in RA is highly variable between patients and the cumulative level is associated with the level of (persistent) inflammation. Nonetheless, several studies have provided evidence for incomplete coupling between inflammation and destruction, which suggests that other individual factors also play a role. This suggestion leads to the hypothesis that candidate genes for the severity of RA are, on the one hand, genes that regulate the level of (local) inflammation and, on the other, genes that mediate the process of bone/cartilage destruction or bone regeneration. Gómez-Vaquero et al.23) found significant difference in HAQ scores, ESR, current and accumulated dose of glucocorticoids, and number of disease modifying antirheumatic drugs taken, and lower serum hemoglobin and albumin according to the BsmI VDR gene polymorphisms. For the association analysis of RA-related phenotypes, we assessed whether GnRH1 polymorphisms are related to laboratory markers of disease activity such as ESR and bone erosion. Based on ESR (30 mm/hr) and the existence of bone erosion, we divided the RA patients into two clinical subgroups and analyzed the relationship between clinical data of RA and GnRH1 polymorphisms. Significant genetic associations were found between GnRH1 polymorphisms and clinical data of RA, such as ESR. RA patients were slightly less likely to have elevated ESR/low CRP.24) Infection, renal insufficiency, and low albumin were associated with having elevated ESR/low CRP; low albumin predicted elevated CRP/low ESR and elevated ESR/low CRP discordance. RA patients were less likely to have elevated ESR/depressed CRP. ESR as a measure inflammation in systemic rheumatic disease may be limited in settings of infection, renal insufficiency, and low albumin.
GnRH plays an important role in immune system modulation. GnRH and its receptor are produced locally by immune cells, suggesting an autocrine role for GnRH. Several studies showed that exogenous GnRH stimulates actions in immune response, such as increase the levels of interleukin (IL)-2 receptor, interferon-γ (IFN-γ) and T helper cells. Jacobson et al.13) suggested that GnRH might play a role in the exacerbation of autoimmune disorders. They assessed disease severity in intact and castrated, male and female, systemic lupus erythematosus (SLE)-prone (SWR×NZB) F1 hybrid mice during treatment with a GnRH agonist and antagonist. After 12 weeks of GnRH antagonist administration, total serum levels of IgG and anti-DNA antibodies were reduced, severity of renal disease decreased and survival improved significantly. On the other hand, they also demonstrated sexually dimorphic actions (disease severity was higher in females) of GnRH agonist. Based on their data, these results may partially explain gender differences in immune function.

Summary

In summary, the present study suggests that GnRH1 polymorphisms might not be associated with susceptibility to RA in female population. However, this study suggests that GnR1 gene polymorphisms may play a relevant role in the susceptibility to severity of ESR level in female RA patients. Additional studies, with a larger number of patients and in different populations, may be required to assess whether GnRH1 polymorphisms and these haplotype could be used as susceptibility or resistance markers in RA disease. This study is the first to analyze the association between SNPs of GnRH1 and RA, and may contribute to future studies on GnRH1.

Figures and Tables

Figure 1
Schematic of the gonadotrophin-releasing hormone (GnRH1) gene. Shown are the 3 exons and 2 introns of the GnRH1 gene and the approximate location of each polymorphism identified in the present study.
jkoa-45-379-g001
Table 1
List of Selected SNPs and Context Sequence of Taqman® Probe in GnRH1 Gene
jkoa-45-379-i001
Table 2
Demographics and Clinical Variables in Controls and RA
jkoa-45-379-i002
Table 3
Genotype Distributions and Allele Frequencies of Each Polymorphism of the GnRH1 in the RA and Control Groups
jkoa-45-379-i003
Table 4
Genotype Distributions and Allele Frequencies of Each Polymorphism of the GnRH1 in the ESR<30 and ESR≥30 Groups of RA
jkoa-45-379-i004
Table 5
Logistic Regression Analysis of 3 SNPs in GnRH1 between the ESR<30 and ESR ≥30 Groups of RA
jkoa-45-379-i005
Table 6
Genotype Distributions and Allele Frequencies of Each Polymorphism of the GnRH1 in the Existence of Bone Erosion of RA
jkoa-45-379-i006
Table 7
Haplotype Analysis Consisting of Significant Polymorphisms in the Genes between the ESR<30 and ESR≥30 Groups of RA
jkoa-45-379-i007

Notes

*This work is supported by Clinical Medicine Research Grant (2007) funded by Kangwon National University Hospital.

References

1. van der Helm-van Mil AH, Breedveld FC, Huizinga TW. Aspects of early arthritis. Definition of disease states in early arthritis: remission versus minimal disease activity. Arthritis Res Ther. 2006. 8:216.
2. Gabriel SE, Michaud K. Epidemiological studies in incidence, prevalence, mortality, and comorbidity of the rheumatic diseases. Arthritis Res Ther. 2009. 11:229.
crossref
3. Gerosa M, De Angelis V, Riboldi P, Meroni PL. Rheumatoid arthritis: a female challenge. Womens Health (Lond Engl). 2008. 4:195–201.
crossref
4. Silverman MN, Sternberg EM. Neuroendocrine-immune interactions in rheumatoid arthritis: mechanisms of glucocorticoid resistance. Neuroimmunomodulation. 2008. 15:19–28.
crossref
5. Graff RJ, Lappé MA, Snell GD. The influence of the gonads and adrenal glands on the immune response to skin grafts. Transplantation. 1969. 7:105–111.
crossref
6. Millar RP, Lu ZL, Pawson AJ, Flanagan CA, Morgan K, Maudsley SR. Gonadotropin-releasing hormone receptors. Endocr Rev. 2004. 25:235–275.
crossref
7. Fink G. Oestrogen and progesterone interactions in the control of gonadotrophin and prolactin secretion. J Steroid Biochem. 1988. 30:169–178.
crossref
8. Chen ZG, Yu KL, Zheng HM, Dong KW. Estrogen receptor-mediated repression of gonadotropin-releasing hormone (gnRH) promoter activity in transfected CHO-K1 cells. Mol Cell Endocrinol. 1999. 158:131–142.
crossref
9. Emanuele NV, Emanuele MA, Tentler J, Kirsteins L, Azad N, Lawrence AM. Rat spleen lymphocytes contain an immunoactive and bioactive luteinizing hormone-releasing hormone. Endocrinology. 1990. 126:2482–2486.
crossref
10. Maier CC, Marchetti B, LeBoeuf RD, Blalock JE. Thymocytes express a mRNA that is identical to hypothalamic luteinizing hormone-releasing hormone mRNA. Cell Mol Neurobiol. 1992. 12:447–454.
crossref
11. Weesner GD, Becker BA, Matteri RL. Expression of luteinizing hormone-releasing hormone and its receptor in porcine immune tissues. Life Sci. 1997. 61:1643–1649.
crossref
12. Batticane N, Morale MC, Gallo F, Farinella Z, Marchetti B. Luteinizing hormone-releasing hormone signaling at the lymphocyte involves stimulation of interleukin-2 receptor expression. Endocrinology. 1991. 129:277–286.
crossref
13. Jacobson JD, Nisula BC, Steinberg AD. Modulation of the expression of murine lupus by gonadotropin-releasing hormone analogs. Endocrinology. 1994. 134:2516–2523.
crossref
14. Livak KJ. Allelic discrimination using fluorogenic probes and the 5' nuclease assay. Genet Anal. 1999. 14:143–149.
crossref
15. Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006. 22:1928–1929.
crossref
16. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005. 21:263–265.
crossref
17. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001. 68:978–989.
crossref
18. Inform G. HapAnalyzer: minimum haplotype analysis system for association studies. Genomics Inform. 2004. 2:107–109.
19. Newton JL, Harney SM, Timms AE, et al. Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis. Arthritis Rheum. 2004. 50:2122–2129.
crossref
20. Hinks A, Barton A, John S, et al. Association between the PTPN22 gene and rheumatoid arthritis and juvenile idiopathic arthritis in a UK population: further support that PTPN22 is an autoimmunity gene. Arthritis Rheum. 2005. 52:1694–1699.
21. Plenge RM, Padyukov L, Remmers EF, et al. Replication of putative candidate-gene associations with rheumatoid arthritis in >4,000 samples from North America and Sweden: association of susceptibility with PTPN22, CTLA4, and PADI4. Am J Hum Genet. 2005. 77:1044–1060.
22. van der Helm-van Mil AH, Kern M, Gregersen PK, Huizinga TW. Variation in radiologic joint destruction in rheumatoid arthritis differs between monozygotic and dizygotic twins and pairs of unrelated patients. Arthritis Rheum. 2006. 54:2028–2030.
crossref
23. Gómez-Vaquero C, Fiter J, Enjuanes A, Nogués X, Díez-Pérez A, Nolla JM. Influence of the BsmI polymorphism of the vitamin D receptor gene on rheumatoid arthritis clinical activity. J Rheumatol. 2007. 34:1823–1826.
24. Costenbader KH, Chibnik LB, Schur PH. Discordance between erythrocyte sedimentation rate and C-reactive protein measurements: clinical significance. Clin Exp Rheumatol. 2007. 25:746–749.
TOOLS
Similar articles