Galectin-9 gene (LGALS9) polymorphisms are associated with rheumatoid arthritis in Brazilian patients

Abstract

Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease characterized by synovial inflammation and hyperplasia, as well as cartilage and bone destruction. Several proteins are associated with the pathogenesis of the disease. Galectin-9 belongs to the family of lectins that are involved in various biological processes and have anti-inflammatory activity.

Objective

To investigate associations between the SNPs of the GAL-9 gene (LGALS9) and serum levels in rheumatoid arthritis patients. We extracted DNA from 356 subjects, 156 RA patients and 200 healthy controls from northeastern Brazil. Three polymorphisms (rs4795835, rs3763959, and rs4239242) in the LGALS9 gene were selected and genotyped using TaqMan SNP genotyping assay. Serum concentrations of galectin-9 were analyzed by ELISA.

Results

The rs4239242 TT genotype showed a positive association with RA (p = 0.0032, odds ratio = 0.28), and heterozygous TC were prevalent in the control group compared to RA patients (p = 0.0001, odds ratio = 7.99). Galectin-9 serum levels were significantly increased in RA patients compared to the control group (p<0.0001). Patients in remission had high levels of galectin compared to the moderate activity group (p<0.0001). Regarding the Clinical Disease Activity Index (CDAI), patients in remission or low activity presented high levels of galectin when compared to patients in severity (p<0.0001). Patients performing moderate activity had a significant value compared to patients who were in high disease severity (p = 0.0064). Interestingly, the AG genotype (rs3763959) has been associated with a higher presence of bone erosion in RA patients (p = 0.0436). The SNP rs4239242 TT genotype showed a positive association with RA in comparison to the control group. The AG genotype (rs3763959) has been associated with a higher presence of bone erosion in RA patients.

About the Authors:

Kamila de Melo Vilar

Roles Methodology, Writing – original draft

Affiliation: Laboratory of Immunomodulation and New Therapeutic Approaches (LINAT), Suely-Galdino Therapeutic Innovation Research Center (NUPIT-SG), Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil

Michelly Cristiny Pereira

Roles Investigation, Writing – review & editing

Affiliation: Laboratory of Immunomodulation and New Therapeutic Approaches (LINAT), Suely-Galdino Therapeutic Innovation Research Center (NUPIT-SG), Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil

ORCID logo http://orcid.org/0000-0002-1672-8202

Andrea Tavares Dantas

Roles Conceptualization, Formal analysis

Affiliations Laboratory of Immunomodulation and New Therapeutic Approaches (LINAT), Suely-Galdino Therapeutic Innovation Research Center (NUPIT-SG), Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil, Rheumatology Service, Clinical Hospital, UFPE, Recife, Pernambuco, Brazil

Moacyr Jesus Barreto de Melo Rêgo

Roles Funding acquisition, Supervision, Writing – review & editing

Affiliation: Laboratory of Immunomodulation and New Therapeutic Approaches (LINAT), Suely-Galdino Therapeutic Innovation Research Center (NUPIT-SG), Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil

Ivan da Rocha Pitta

Roles Funding acquisition, Supervision

Affiliation: Laboratory for Planning and Synthesis of Drugs (LPSF) of the Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil

Ângela Luzia Branco Pinto Duarte

Roles Formal analysis, Project administration

Affiliation: Rheumatology Service, Clinical Hospital, UFPE, Recife, Pernambuco, Brazil

Maira Galdino da Rocha Pitta

Roles Formal analysis, Funding acquisition, Writing – review & editing

* E-mail: [email protected]

Affiliation: Laboratory of Immunomodulation and New Therapeutic Approaches (LINAT), Suely-Galdino Therapeutic Innovation Research Center (NUPIT-SG), Federal University of Pernambuco (UFPE), Recife, Pernambuco, Brazil

Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease characterized by synovial inflammation and hyperplasia, as well as cartilage and bone destruction, and it may also cause extra-synovial symptoms [1]. The important characteristic of RA is the infiltration of multiple leukocytes into the joints, including B cells, T cells, macrophages, dendritic cells, and neutrophils. Numerous proteins have been shown to play a role in RA pathogenesis, such as autoantibodies used against citrullinated antigens (ACPA). TNFα, IL-1, IL-6, and Th17 cells have recently been identified in the promotion of autoimmune processes, joint destruction, and angiogenesis [2,3]. Its pathogenesis is related to a combination of genetic, environmental, hormonal, smoking, alcohol consumption, and infectious factors, and its progression requires activation of innate immunity [4].

The family of galectins is involved in a range of biological processes. Galectin-9 belongs to the group containing two distinct carbohydrate recognition domains (CRDs) linked by a peptide linker and having anti-inflammatory activity. Galectin-9 is expressed by T cells, macrophages, fibroblasts, and endothelial cells, playing an important role in the regulation of inflammation and immune responses by down-regulating pro-inflammatory T cells [5,6]. An important cell surface receptor for galectin-9 is T cell immunoglobulin and mucin-domain-containing molecule-3 (Tim-3), with an expression on CD4+ Th1 cells, CD8+ cytotoxic T cells, and CD11b+ dendritic cells (DC), but not on Th2 cells or macrophages [7].

Galectin-9 is also related to the regulation of the differentiation of T cell subsets in vitro and in vivo. In cell culture, treatment with galectin-9 induced differentiation of T cells to regulatory T cells (Treg) and suppressed the differentiation of Th17 cells [8]. One study evaluated the expression levels of galectin-9 mRNA in peripheral blood mononuclear cells (PBMCs), where expression was significantly lower in patients suffering from RA with moderate to high disease activity when compared to patients with low disease activity [9].

To date, there are no studies correlating polymorphisms with the susceptibility for rheumatoid arthritis. In this context, our study was the first to analyze whether genetic polymorphisms in the LGALS9 gene were associated with susceptibility for RA development. Based on the broad functions of galectin-9, we hypothesized that functional genetic polymorphism in the gene coding for galectin-9 (LGALS9) would be associated with susceptibility for rheumatoid arthritis.

Materials and methods

Patients and controls

RA patient blood samples were acquired from the Rheumatology Outpatient Clinical Hospital Federal University of Pernambuco (UFPE). A total of 156 patients fulfilled four or more American College of Rheumatology (ACR) 1987 classification criteria [10]. Demographic, clinical, and laboratory data were collected from hospital records. Laboratory features of RA patients such as erythrocyte sedimentation rate (ESR) and rheumatoid factor positivity were recorded. Individual disease activity was quantified using the disease activity score (DAS28) [11] and the Clinical Activity Index (CDAI) [12]. The Health Assessment Questionnaire (HAQ) [13] score was also applied to patients to assess functional disability in arthritis. Radiographs of hands were obtained from patients with RA and evaluated for the presence of erosions by an experienced rheumatologist blinded to the clinical data. The control group consisted of 200 matched, unrelated healthy blood donors free of any rheumatologic disease. The clinical data pertaining to patients with RA are presented in Table 1. All subjects gave their written consent to participate. The study was approved by the research ethics committee (CCAE-63648616.0.0000.5208) of the UFPE.

[Figure omitted. See PDF.]

Table 1. Demographic and clinical parameters of RA patients and control group.

https://doi.org/10.1371/journal.pone.0223191.t001

SNP selection and genotyping

Single Nucleotide Polymorphisms (SNPs) for the galectin-9 gene (LGALS9) were identified by querying the HapMap. The LGALS9 gene is located on human chromosome 17 (Chr17:25,958,147–25,976,586), with a minor allele frequency (MAF) of at least 0.1% in Caucasians Genomic. The SNPs rs4795835, rs3763959, and rs4239242 were evaluated. Genomic DNA was extracted from the peripheral blood of RA patients and controls using the technical standard phenol-chloroform. Allelic discrimination with TaqMan real-time PCR was used to genotype for polymorphisms in the LGALS9. The PCR products were generated in a 12 μl reaction containing 20 ng DNA, and the cycling parameters were denaturation at 95°C for 10 minutes, followed by 45 cycles of denaturation at 95°C for 15 seconds and the annealing and extension at 60°C for 1 minute in a StepOnePlus Real-Time PCR System according to the conditions recommended by the manufacturer (Applied Biosystems, Foster City, CA, USA).

Serum Galectin-9 levels determination

Cytokine serum level was determined by Human galectin-9 ELISA Kit (R&D Systems, Ltd., Abingdon, UK) according to the manufacturer’s recommendation. The lower limit of detection for the ELISA galectin-9 kit was 15.6 pg/ml.

Statistical analysis

All genotype data were checked for derivation from Hardy-Weinberg equilibrium (HWE) using http://www.oege.org/software/hardy-weinberg.html [14]. The differences in allele/genotype frequencies between patients and controls were analyzed by Fisher’s exact test. The association of serum Galectin-9 levels with clinical and laboratory measures of patients with RA were analyzed by univariate comparisons using nonparametric tests (Mann-Whitney tests). A logistic regression analysis was performed to detect potential associations in our case-control sample and associations with clinical parameters. Statistical significance was assumed at p<0.05. Statistical analyses of the data were performed using the GraphPad Prism (version 6.0) statistical program.

Results

A total of 356 subjects were genotyped, among whom 156 RA patients and 200 controls were age and matched controls. The mean age for RA patients and controls was 52.15 ± 12.10 years and 48.01 ± 12.40 years, respectively. Information on the main clinical characteristics of the RA patients is summarized in Table 1.

We examined serum galectin-9 levels of RA patients and the control group as shown in Fig 1. Galectin-9 serum levels were significantly increased in RA patients compared to the control group (median 27.82 pg/ml; p<0.0001) (Fig 1A). Patients in remission (DAS28<3.2) had high levels of galectin compared to the moderate activity group (DAS28 3.2> <5.1) (p <0.0001). We evaluated the group in remission in relation to the group with high disease activity(DAS28>5.1) (p <0.0001) (Fig 1B). This finding may be related to mechanisms not yet known or to the use of drugs that modify the course of the disease, as shown in Table 1.

[Figure omitted. See PDF.]

Fig 1. Association of serum Galectin-9 levels in patients with rheumatoid arthritis.

(A) serum levels of galectin-9 in patients with RA and control group. (B) association of serum levels of galectin-9 with disease activity score (DAS28). (C) association of serum levels of galectin-9 with Clinical Disease Activity Index (CDAI).

https://doi.org/10.1371/journal.pone.0223191.g001

The genotype frequencies of rs4239242, rs3763959, and rs4795835 were distributed according to the Hardy–Weinberg equilibrium in normal controls (all p>0.05). The genotypes and alleles’ frequencies in RA patients and control individuals are shown in Table 2. The rs4239242 TT genotype showed a positive association with RA (15.59% in RA patients and 5.20% in controls, p = 0.0032, odds ratio = 0.28; 95%IC = 0.12–0.64). Heterozygous TC were prevalent in the control group (66.14%) compared to RA patients (19.26%; p = 0.0001, odds ratio = 7.99; 95%IC = 4.59–13.91). Comparison of the homozygous CC versus TC genotypes showed association between patients and controls (odds ratio = 3.36; p = 0.0054; 95%IC = 3.36–7.64). The rs3763959 AA genotype showed a negative association with RA (6.67% in RA patients and 15.10% in controls, p = 0.040, odds ratio = 2.49; 95%IC = 1.05–5.90). Homozygous GG were prevalent in RA patients (55.23%) compared to the control group (34.89%; p = 0.0009, odds ratio = 0.43; 95%IC = 0.26–0.70). Comparison of the homozygous GG versus AG genotypes showed association between patients and controls (odds ratio = 0.40; p = 0.040; 95%IC = 0.16–0.95). On the other hand, the allele and genotype frequencies of rs4795835 showed no association with RA (Table 2).

[Figure omitted. See PDF.]

Table 2. Distribution of the LGALS9 alleles and genotypes in patients with RA and healthy individuals.

https://doi.org/10.1371/journal.pone.0223191.t002

Regarding the Clinical Disease Activity Index (CDAI), patients in remission or low activity (CDA<10) presented high levels of galectin when compared to patients in severity (CDA>22) (p<0.0001). We also evaluated patients with CDAI>10<22 (moderate) and CDAI>22, where patients in moderate activity had a significant value compared to patients who were in high disease severity (p = 0.0064) (Fig 1C). The Health Assessment Questionnaire (HAQ) score to assess functional disability in rheumatoid arthritis and rheumatoid factor was also evaluated, but no significant differences were found based on galectin-9 levels.

To address the functional significance of the selected polymorphisms in the galectin-9 gene, we examined serum levels of galectin-9 with each genotype of RA patients, however, genotypic and allelic frequencies were not associated with serum levels. There were no statistical associations between galectin-9 polymorphisms with DAS28, CDAI, HAQ, ERS, or rheumatoid factor positivity. Interestingly, the AG genotype (rs3763959) has been associated with a higher presence of bone erosions in RA patients (p = 0.0436) (Fig 2).

[Figure omitted. See PDF.]

Fig 2. Erosion stratified by genotypes of LGALS9 gene polymorphisms (rs3763959) in RA patients.

https://doi.org/10.1371/journal.pone.0223191.g002

Discussion

Many studies have evaluated the role of polymorphisms in galectin genes in rheumatoid arthritis [15,16,17,18]. However, to date there is no study on polymorphisms in the galectin-9 gene in RA. Therefore, this study was carried out to investigate the association of LGALS9 gene polymorphisms in Brazilian patients with rheumatoid arthritis. We did not find associations between the polymorphisms selected and the serum levels of galectin.

Our results show the first positive association of the rs4239242 TT genotype with RA. While genotype AA rs3763959 showed a negative association with RA, genotype frequencies of rs4795835 showed no association with RA. Rosen et al. [19] reported that LGALS9 polymorphisms are associated with the risk of developing alcoholic liver disease (ALD), where an association between the T rs4239242 allele and the protection against ALD and the presence of the C allele were associated with a 70% difference in the risk of developing ALD. This result differs from that found where the TT genotype showed a positive association with RA (p = 0.0032).

The rs3751093 was evaluated, and the proportion of G allele carriers was significantly higher among subjects who were at risk of developing alcoholic liver injury compared to those who were more likely to have the A allele bring on ALD. For different diseases, the polymorphism can generate diverse effects in the individual [19]. The AG genotype (rs3763959) has been associated with a higher presence of bone erosion in RA patients (Fig 2). Moriyama et al. [20] suggest that galectin-9 is indeed present in the bone marrow cavity in vivo. The Tim-3 surface receptor is crucial for the exogenous function of galectin-9, where this complex is related to the induction of apoptosis of CD4+ T-helper 1 (Th1), Th17, CD8+ cells, and cytotoxic T cells. Therefore, galactin-9 is an important molecule in the regulation of the immune system due to induction of apoptosis in Th1 and Th17 cells because the activation of stem cells controls the inflammatory process [21]. In order to clarify the Tim-3/gal-9 complex, a large amount of galectin-9 was administered to rats with arthritis, and bone destruction was consequently observed. A recombinant galectin-9 was administered, respectively, and it was found that recombinant galectin-9 lightly suppressed joint inflammation, however, the suppression did not reach significance according to clinical scores. Thus, the bone destruction that accompanies chronic inflammatory disease permeated by osteoclasts can be attenuated by the presence of the potential control system Tim-3/galectin-9, especially at the sites of destruction in RA, indicating that the presence of this complex Tim-3/galectin-9 can lead to increased inflammation and inflammatory bone destruction [20].

There have been no reports regarding an association between rheumatoid arthritis and LGALS9 polymorphism, although galectin-9 has been shown to play an important role in RA pathophysiology [22,23]. Patients with RA are associated with an inflammation in chronic synovium that leads to hyperplasia and activation of a set of fibroblasts similar to macrophages that express degradative enzymes, causing the destruction of joint structures. Several cell types are recruited into the inflammatory environment that develops within the synovium, where the fibroblasts, macrophages, and proinflammatory B and T cells contribute significantly to the pathology of the disease [24,25].

Most polymorphisms evaluated (genotypes and polymorphic alleles) had a protective effect. Pearson et al. [21] evaluated the difference in abundant galectin-9 expression in synovial fibroblasts compared to dermal fibroblasts. Thus, endogenous galectin-9 protects against apoptosis and increases the viability of synovial fibroblasts [9]. Other studies have found that galectin-9 is being used in patients with RA and may be beneficial in collagen-induced arthritis of mice when administered exogenously. In RA patients, a decrease in galectin-9-Tim-3 complex signalling was observed. Tim-3 expression levels in CD4+ T cells from RA patients were lower compared to those from healthy controls, resulting in a decrease in galectin-9 mediated apoptosis of CD4+ T cells [26]. Thus, it is suggested that this galectin has both pathogenic and pro-inflammatory roles.

One of the limitations of our study was that the vast majority of patients were in clinical remission, and perhaps because of this the frequency of protection alleles was higher and might not have shown the actual function of the polymorphism. Another limitation was the evaluation of the functional significance of the polymorphisms selected in the galectin-9 gene. When we examined serum levels of galectin-9 with each genotype of RA patients, however, genotypic and allelic frequencies were not associated with serum levels. Taking into account that AR is a microenvironment-dependent disease, polymorphisms may provide some function with galectin in the synovial fluid rather than at the systemic level. In addition, this polymorphism may be influenced by other proteins related to the synthesis, excretion, or degradation of galectin-9. Nevertheless, our study is the first to investigate the role of LGALS9 polymorphism in RA patients, contributing to the knowledge of the genetic background of RA.

Conclusion

Our results suggest for the first time that genetic alterations in the LGALS9 gene are associated with rheumatoid arthritis in the Brazilian population. The SNP rs4239242 TT genotype showed a positive association with RA in comparison with the control group. As the AG genotype (rs3763959) has been associated with a higher presence of bone erosion in RA patients, these findings suggest that polymorphisms in LGALS9 may be an interesting aspect of rheumatoid arthritis for further exploration.

Supporting information

[Figure omitted. See PDF.]

S1 File. Supporting Information files.

https://doi.org/10.1371/journal.pone.0223191.s001

(XLSX)

Acknowledgments

We are grateful to patients with RA who agreed to participate. We also thank the Rheumatology Service, Clinical Hospital, for providing the collection of data and blood samples from the patients.

Citation: Vilar KdM, Pereira MC, Tavares Dantas A, de Melo Rêgo MJB, Pitta IdR, Pinto Duarte ÂLB, et al. (2019) Galectin-9 gene (LGALS9) polymorphisms are associated with rheumatoid arthritis in Brazilian patients. PLoS ONE 14(10): e0223191. https://doi.org/10.1371/journal.pone.0223191

References

1. Doody KM, Bottini N, Firestein GS. Epigenetic alterations in rheumatoid arthritis fibroblast-like synoviocytes. Epigenomics. 2017 Apr;9(4):479–492. pmid:28322585

2. Glant TT, Mikecz K, Rauch TA. Epigenetics in the pathogenesis of rheumatoid arthritis. BMC Med. 2014 Feb 26;12:35. pmid:24568138

3. Yap HY, Tee SZ, Wong MM, Chow SK, Peh SC, Teow SY. “Pathogenic Role of Immune Cells in Rheumatoid Arthritis: Implications in Clinical Treatment and Biomarker Development.” Cells vol. 7,10 161. 9 Oct. 2018.

4. Dhaouadi T, Chahbi M, Haouami Y, Sfar I, Abdelmoula L, Ben Abdallah T, et al. IL-17A, IL-17RC polymorphisms and IL17 plasma levels in Tunisian patients with rheumatoid arthritis. PLoS One. 2018 Mar 27;13(3):e0194883.

5. John S, Mishra R. Galectin-9: From cell biology to complex disease dynamics. J Biosci. 2016 Sep;41(3):507–34. pmid:27581941

6. Dai SY, Nakagawa R, Itoh A, Murakami H, Kashio Y, Abe H, et al. Galectin-9 induces maturation of human monocyte-derived dendritic cells. J Immunol. 2005 Sep 1;175(5):2974–81. pmid:16116184

7. Ocaña-Guzman R, Torre-Bouscoulet L, Sada-Ovalle I. TIM-3 Regulates Distinct Functions in Macrophages. Front Immunol. 2016 Jun 13;7:229. pmid:27379093

8. Seki M, Oomizu S, Sakata KM, Sakata A, Arikawa T, Watanabe K, et al. Galectin-9 suppresses the generation of Th17, promotes the induction of regulatory T cells, and regulates experimental autoimmune arthritis. Clin Immunol. 2008 Apr;127(1):78–88. pmid:18282810

9. Seki M, Sakata KM, Oomizu S, Arikawa T, Sakata A, Ueno M, et al. Beneficial effect of galectin 9 on rheumatoid arthritis by induction of apoptosis of synovial fibroblasts. Arthritis Rheum. 2007; 56:3968–3976. pmid:18050192

10. Arnett FC, Edworthy SM, Bioch DA, McShane DJ, Fries JF, Copper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31:315–24. pmid:3358796

11. Prevoo ML, van 't Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, van Riel PL. Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum. 1995 Jan; 38(1):44–8. pmid:7818570

12. Aletaha D, Nell VP, Stamm T, Uffmann M, Pflugbeil S, Machold K, et al. Acute phase reactants add little to composite disease activity indices for rheumatoid arthritis: validation of a clinical activity score. Arthritis Res Ther. 2005; 7(4):R796–806. pmid:15987481

13. Ferraz MB, Oliveira LM, Araujo PM, Atra E, Tugwell P. Crosscultural reliability of the physical ability dimension of the health assessment questionnaire. J Rheumatol. 1990 Jun; 17(6):813–7. pmid:2388204

14. Rodriguez Santiago,Gaunt Tom R. and Day Ian N. M. Hardy-Weinberg Equilibrium Testing of Biological Ascertainment for Mendelian Randomization Studies. American Journal of Epidemiology Advance Access published on January 6,2009.

15. Pál Z, Antal P, Srivastava SK, Hullám G, Semsei AF, Gál J, et al. Non-synonymous single nucleotide polymorphisms in genes for immunoregulatory galectins: association of galectin-8 (F19Y) occurrence with autoimmune diseases in a Caucasian population. Biochim Biophys Acta. 2012 Oct;1820(10):1512–8. pmid:22683700

16. Hu CY, Chang SK, Wu CS, Tsai WI, Hsu PN. Galectin-3 gene (LGALS3) +292C allele is a genetic predisposition factor for rheumatoid arthritis in Taiwan. Clin Rheumatol. 2011 Sep;30(9):1227–33. pmid:21475983

17. Panoulas VF, Douglas KM, Smith JP, Metsios GS, Elisaf MS, Nightingale P, et al. Galectin-2 (LGALS2) 3279C/T polymorphism may be independently associated with diastolic blood pressure in patients with rheumatoid arthritis. Clin Exp Hypertens. 2009 Apr;31(2):93–104. pmid:19330599

18. Panoulas VF, Nikas SN, Smith JP, Douglas KM, Nightingale P, Milionis HJ, et al. Lymphotoxin 252A>G polymorphism is common and associates with myocardial infarction in patients with rheumatoid arthritis. Ann Rheum Dis. 2008 Nov;67(11):1550–6. pmid:18230628

19. Rosen HR, Golden-Mason L, Daly AK, Yang I, Day CP. Variants in the LGALS9 Gene Are Associated With Development of Liver Disease in Heavy Consumers of Alcohol. Clin Gastroenterol Hepatol. 2016 May;14(5):762–8.e1. pmid:26598225

20. Moriyama K, Kukita A, Li YJ, Uehara N, Zhang JQ, Takahashi I, et al. Regulation of osteoclastogenesis through Tim-3: possible involvement of the Tim-3/galectin-9 system in the modulation of inflammatory bone destruction. Lab Invest. 2014; 1200–11. pmid:25264706

21. Pearson MJ, Bik MA, Ospelt C, Naylor AJ, Wehmeyer C, Jones SW, et al. Endogenous Galectin-9 Suppresses Apoptosis in Human Rheumatoid Arthritis Synovial Fibroblasts. Sci Rep. 2018 Aug 27;8(1):12887. pmid:30150656

22. Lee J, Park EJ, Noh JW, Hwang JW, Bae EK, et al. Underexpression of TIM-3 and blunted galectin-9-induced apoptosis of CD4+ T cells in rheumatoid arthritis. Inflammation. 2012; 35:633–637. pmid:21717191

23. O'Brien MJ, Shu Q, Stinson WA, Tsou PS, Ruth JH, Isozaki T, et al. A unique role for galectin-9 in angiogenesis and inflammatory arthritis. Arthritis Res Ther. 2018 Feb 12;20(1):31. pmid:29433546

24. Bradfield PF, Amft N, Vernon-Wilson E, Exley AE, Parsonage G, Rainger GE, et al. Rheumatoid fibroblast-like synoviocytes overexpress the chemokine stromal cell-derived factor 1 (CXCL12),which supports distinct patterns and rates of CD4+ and CD8+ T cell migration within synovial tissue. Arthritis Rheum. 2003 Sep;48(9):2472–82. pmid:13130466

25. Yeo L, Toellner KM, Salmon M, Filer A, Buckley CD, Raza K, et al. Cytokine mRNA profiling identifies B cells as a major source of RANKL in rheumatoid arthritis. Ann Rheum Dis. 2011 Nov;70(11):2022–8. pmid:21742639

26. Arikawa T, Watanabe K, Seki M, Matsukawa A, Oomizu S, Sakata KM, et al. Galectin-9 ameliorates immune complex-induced arthritis by regulating Fc gamma R expression on macrophages. Clin Immunol. 2009; 133:382–392. pmid:19800850