Skip to main content
  • Research Article
  • Open access
  • Published:

Genetic variant in SPAG16 is associated with the susceptibility of ACPA-positive rheumatoid arthritis possibly via regulation of MMP-3

Abstract

Objectives

In two previously published genome-wide association studies, a cluster of variants of sperm-associated antigen16 (SPAG16) were reported to be associated with the radiological progression rate of ACPA-positive rheumatoid arthritis (RA) patients from North American and Southern European ancestry. In this study, we aimed to investigate whether the reported RA-risk loci in SPAG16 are associated with the disease in the Chinese population and to further validate the functional role of the susceptible locus in RA tissues.

Methods

A total of 500 ACPA-positive RA patients and 1000 age-matched healthy subjects were recruited. Two SNPs of SPAG16, including rs7607479 (C/T) and rs6435818 (A/C), were genotyped, and the genotyping data were compared with chi-square test. Gene expression analysis was performed in synovial tissues obtained from 40 RA patients and 30 non-RA controls surgically treated for bone fracture. The tissue expression of SPAG16 and matrix metalloproteinase 3 (MMP-3) was compared between the two groups by the Student’s t test. The relationship between serum indexes and mRNA expression of SPAG16 and MMP-3 were evaluated by Spearman’s correlation analysis.

Result

For rs7607479, the frequency of genotype TT was significantly higher in RA patients than in the controls (49.0% vs. 40.4%, p = 0.002). The RA patients were found to have significantly lower frequency of allele C than the controls (30.9% vs. 36.8%, p = 0.001). As for rs6435818, there was no significant difference of genotype or allele frequency between the two groups. The mRNA expression of MMP-3 was 1.63-fold higher in the RA patients than in the controls (p < 0.001). The expression of SPAG16 was comparable between the two groups (p = 0.43). The mRNA expression of MMP-3 was 1.39-fold higher in patients with genotype TT than in the patients with genotype CC (p = 0.006). The mRNA expression level of MMP-3 was significantly correlated with serum rheumatoid factor (r = 0.498, p < 0.001) and C-reactive protein (r = 0.272, p = 0.01), weakly correlated with erythrocyte sedimentation rate (r = 0.236, p = 0.09).

Conclusions

We validated a common genetic risk factor in ACPA-positive patients with RA, which is associated with the tissue production of MMP-3 and disease progression. Further functional analysis into the role of rs7607479 in MMP-3 expression can shed new light on the genetic architecture of ACPA-positive RA.

Introduction

Rheumatoid arthritis (RA) is a disease characterized by inflammatory changes of the synovial tissues [1, 2]. RA arises based on both genetic and epigenetic components including direct cigarette smoking, dust exposure and other environmental factors [3,4,5]. Previous studies have well documented the role of genetic predisposition in the development of RA [6,7,8]. As reported in earlier studies, the rate of RA for the general population is approximately 0.2–1% [9, 10]. For identical twins, the concordance of RA was reported to be 15% [11], which was much higher than that in the general population and thus indicated the importance of genetic factors for RA. In the past decades, a large number of genetic risk loci associated with RA has been discovered through genome-wide association studies (GWASs) in various populations worldwide [12,13,14]. The contribution of the genetic variation to RA was estimated to be about 60% [12,13,14].

For patients with RA, the presence of rheumatoid factor (RF), elevated C-reactive protein (CRP) level and high erythrocyte sedimentation rate (ESR) were reported to indicate the severity of disease [15]. In addition to RF, CRP and ESR, antibodies against citrullinated proteins (ACPA) have been also reported as a risk factor in the development of joint damage in RA patients [16, 17]. Positive detection of ACPA at early stage of RA can be prognostic of an increased likelihood of disease severity [16, 17]. As specific subgroups of RA patients, ACPA-positive individuals were found with higher risk of severe erosive phenotype and mortality rate than ACPA-negative individuals [18,19,20]. This variability of the clinical presentation between these two subgroups might reflect the presence of different underlying genetic risk factors. Analyzing ACPA-positive RA patients is therefore a useful strategy to reduce heterogeneity and increase the power to identify novel genetic factors associated with the severity of RA.

In a GWAS for ACPA-positive RA patients with radiographic joint damage, Knevel et al. [21] reported that a cluster of single nucleotide polymorphisms (SNPs), including rs7607479 (C/T) and rs6435818 (A/C), were associated with the radiological progression rate of RA in European and North American patients. Of note, rs7607479 reached genome-wide significance [21]. Interestingly, in a recently published case–control GWAS using ACPA-positive RA patients [22], rs6435818 was successfully located as a risk locus of RA in Caucasian and of Southern European ancestry. SNPs rs7607479 and rs6435818 were located in the intron of sperm-associated antigen16 (SPAG16), which is predicted to have a role in the basal structure of the primary cilium [23]. To date, there was a lack of knowledge concerning the function of SPAG16. It is mainly expressed in sperm, testis and other tissues such as bone marrow [23]. Knevel et al. [21] confirmed that SPAG16 was also expressed in synovium tissues of RA, while there was a lack of relationship between SPAG16 expression and genotypes of rs7607479. Interestingly, Knevel et al. [21] observed that the genotype of rs7607479 was correlated with the serum level of a RA-associated gene, matrix metalloproteinase 3 (MMP-3). The MMPs are composed of a family of zinc- and calcium-dependent enzymes that contribute to the destruction of articular cartilage [24,25,26,27]. Abundantly expressed in the synovium tissues, MMP-3 is considered to be the main MMP involved in cartilage degradation of RA patients [28]. As reported in previous studies [29, 30], MMP-3 may be associated with the progression of the RA. Elevation of serum MMP-3 can be detected at both early stage and advanced stage of RA patients. Significant association between the promoter polymorphism of MMP-3 and joint destruction was reported in patients with RA [31]. To our knowledge, there was a lack of replication of the association between SPAG16 variants and the risk of RA in other populations. Moreover, the relationship between the genotype of SPAG16 variants and the synovium mRNA expression of SPAG16/MMP-3 was worthy of further investigation. In this study, we aimed to investigate whether the reported RA-risk loci rs7607479 and rs6435818 are associated with the disease in the Chinese population and to further validate the functional role of RA-risk loci.

Methods

Subjects

In this multicenter retrospective study, all RA patients fulfilling ACR/EULAR 2010 RA classification who came to our clinics between 2015 and 2021 were screened. Patients with other connective tissue diseases, including systemic lupus erythematosus, dermatomyositis or systemic sclerosis, were excluded from the study. ACPA was analyzed by a third-generation anti-cyclic citrullinated peptides (anti-CCP) enzyme-linked immunosorbent assay (ELISA) kit (DLD Diagnostika, Hamburg, Germany). Under the approval of the Ethics Committee of our institutions, a total of 500 ACPA-positive RA patients were recruited. The controls were composed of 1000 age-matched healthy subjects who underwent routine health examination organized by their employees. All the controls were excluded to have RA or other autoimmune diseases by clinical examinations. The baseline characteristics of the patients were recorded at their first visit, including age, gender, body mass index (BMI), ESR, CRP, and the serum level of RF.

Blood collection

The ethylenediaminetetraacetic (EDTA)-anti-coagulated whole blood samples were collected from the patients and controls for DNA extraction. Genomic DNA was extracted from the leukocyte sediments by DNA extraction kit (QIAGEN, Tokyo, Japan) according to manufacturer’s guideline, which was then stored in − 20 °C for genotyping assays. Written informed consent forms were signed by each subject.

Genotyping assays of rs7607479 and rs6435818

Two SNPs, including rs7607479 and rs6435818, were investigated in this study. Approximately 20 ng of the DNA sample was used for genotyping with TaqMan SNP Genotyping Assay. A total of 30 μl reaction volume, containing 4 μl of distilled deionized water, 4 μl of TaqMan Genotyping assay mix, 12 μl of genomic DNA and 10 μl of the TaqMan Genotyping master mix, was used for polymerase chain reaction (PCR) amplification. The results of genotyping assay were then analyzed on Roche LightCycler 480 II thermocycler (Roche Diagnostic Ltd., Basel, Switzerland). Ten percent of the samples were randomly selected to reproduce the genotyping results. Reproducibility of 100% was confirmed.

Tissue expression of SPAG16 and MMP-3

Synovial tissues were obtained from the knee joint of 40 RA patients with Parker-Pearson needle biopsy. Normal synovial tissues of 30 patients undergoing tibial plateau reconstruction due to traumatic fracture were used as the negative control. The total RNA was extracted from synovial tissues using TRIzol reagent (QIAGEN, Tokyo, Japan), which were then reversely transcribed with the PrimeScript RT Master Mix kit (TaKaRa, Tokyo, Japan). Tissue expression of SPAG16 and MMP-3 was quantified by Quantitative real-time PCR (qPCR) with beta-actin used as endogenous control. The specific primers were listed below, forward 5′- ATGTTCCAGATGTCTACACCCA -3′, reverse 5′- TGTAACTTCAACCCTTTGAGGTC -3′ for SPAG16, forward 5′- AGCAAGGACCTCGTTTTCATT -3′, reverse 5′- GTCAATCCCTGGAAAGTCTTCA -3′ for MMP-3, and forward 5′- CCTCGCCTTTGCCGATCC-3′, reverse 5′- GGATCTTCATGAGGTAGTCAGTC -3′ for beta-actin. All amplification reactions were carried out in triplicate.

Statistical analysis

The SPSS software (version 22.0, Chicago, USA) was applied to statistical analysis. The Hardy–Weinberg equilibrium (HWE) test was performed to detect whether there could be selection bias of the subjects. The frequency of genotype and minor allele was calculated and compared between the cases and the controls by the chi-square test. The effect of 3 genetic models was examined, including the additive model, dominant model and recessive model. The odds ratio (OR) and 95% confidential intervals (CIs) were then calculated for the two SNPs. The tissue expression of SPAG16 and MMP-3 was compared between the two groups by the Student’s t test. One-way ANOVA was performed to compare the gene expression among different genotypes of rs7607479, and the Tukey test was applied to the post hoc pairwise analysis. The relationship between serum indexes and mRNA expression of SPAG16 and MMP-3 was evaluated by Spearman’s correlation analysis. Descriptive data were presented in the form of mean value ± standard deviation (SD). p value less than 0.05 was considered to indicate statistical significance.

Results

Baseline characteristics of the subjects

For association analysis, the patients and the controls were matched in terms of age (42.7 ± 11.9 yrs for cases vs. 43.3 ± 12.6 yrs for controls, p = 0.34), the ratio of male/female (144/356 for cases vs. 323/677 for controls, p = 0.18) and BMI (23.1 ± 4.3 kg/m2 for cases vs. 23.4 ± 5.7 kg/m2 for controls, p = 0.26). The mean value of RF was 92.1 ± 7.9 (IU)/ml. The mean ESR level and CRP level were 34.9 ± 10.3 mm/h and 58.5 ± 21.3 mg/L, respectively. The clinical characteristics of the subjects are presented in Table 1.

Table 1 Demographic data of the subjects

Association of rs7607479 and rs6435818 with RA risk

HWE test showed normal distribution of the genotype frequency of the subjects (p > 0.05). For rs7607479 of SPAG16, the distribution of genotype frequency of RA patients was significantly different from that of the controls (p = 0.005). As shown by the dominant model (Additional file 1: Supplementary table 1), the frequency of genotype TT was significantly different between the two groups (49.0% for cases vs. 40.4% for controls, p = 0.002). RA patients were found to have significantly lower frequency of allele C than the controls (30.9% for cases vs. 36.8% for controls, p = 0.001). The OR was 0.77 with 95% CI ranging from 0.65 to 0.90. As for rs6435818, there was no significant difference regarding genotype or allele frequency between the two groups (Table 2).

Table 2 Case–control comparison of the frequency of the genotype and allele for rs7607479 and rs6435818

Expression of SPAG16 and MMP-3 in joint

For the 40 RA patients and 30 controls included in expression analysis, there was no significant difference in terms of age (42.3 ± 12.4 yrs for cases vs. 43.7 ± 11.2 yrs for controls, p = 0.62), the ratio of male/female (12/28 for cases vs. 14/16 for controls, p = 0.23), or BMI (23.7 ± 3.2 for cases vs. 24.1 ± 5.5 for controls, p = 0.70) between the two groups. As shown in Fig. 1a, the mRNA expression of MMP-3 was 1.63-fold higher in the RA patients than in the controls (0.00725 ± 0.00314 for cases vs. 0.00443 ± 0.00132 for controls, p < 0.001). The expression of SPAG16 was not statistically different between the two groups (0.0000625 ± 0.0000232 for cases vs. 0.0000585 ± 0.0000151 for controls, p = 0.43) (Fig. 1b).

Fig. 1
figure 1

Tissue expression of MMP-3 and SPAG16 in RA patients. a, b Significantly higher expression of MMP-3 was observed in RA tissues as compared with the non-RA controls. There was no significant difference regarding the tissue expression of SPAG16 between the two groups; c, d There was significant difference regarding MMP-3 expression among the three genotypes (p = 0.01). The mRNA expression of MMP-3 was 1.39-fold higher in patients with genotype TT than in the patients with genotype CC (post hoc p = 0.006). There was significant difference regarding MMP-3 expression between genotype CT and genotype TT (post hoc p = 0.002). There was no significant difference between genotype CT and genotype CC (post hoc p = 0.16). There was no significant difference regarding the tissue expression of SPAG16 among the three genotypes

For RA patients, there were 10 cases with genotype CC of rs7607479, 5 cases with genotype CT and 25 cases with genotype TT. There was significant difference regarding MMP-3 expression among the three genotypes (p = 0.01). As shown in Fig. 1c, the mRNA expression of MMP-3 was 1.39-fold higher in patients with genotype TT than in the patients with genotype CC (0.00833 ± 0.00148 for genotype TT vs. 0.00597 ± 0.00214 for genotype CC, post hoc p = 0.006). Also, there was significant difference between genotype CT and genotype TT (0.00457 ± 0.00241 for genotype CT vs. 0.00833 ± 0.00148 for genotype TT, post hoc p = 0.002). There was no significant difference between genotype CT and genotype CC (post hoc p = 0.16). As for the expression of SPAG16, there was no significant difference among the three genotypes (p = 0.95) (Fig. 1d). The tissue expression of MMP-3 and SPAG16 for different genotypes of rs7607479 is summarized in Additional file 1: Table S2.

As shown in Fig. 2, the mRNA expression level of MMP-3 was significantly correlated with serum RF (r = 0.498, p < 0.001) and CRP (r = 0.272, p = 0.01), partially correlated with ESR (r = 0.236, p = 0.09). However, no significant correlation between SPAG16 expression and these 3 clinical variables was found (r = 0.177 p = 0.21 for RF; r = 0.205, p = 0.15 for CRP; r = -0.166, p = 0.24 for ESR).

Fig. 2
figure 2

The relationship between MMP-3 expression and clinical parameters of RA. The mRNA expression level of MMP-3 was significantly correlated with the serum RF (r = 0.498, p < 0.001) and CRP (r = 0.272, p = 0.01). It was partially correlated with ESR (r = 0.236, p = 0.09)

Discussion

Previous GWASs have reported the association of different variants in SPAG16 with the risk or severity of RA in Caucasian, North American and Southern European cohorts [21, 22]. To further validate the role of SPAG16 in the susceptibility of RA, we replicated two previously reported variants and confirmed that rs7607479 of SPAG16 was associated with RA in the Chinese population. We found that allele C of rs7607479 played a protective role in the risk of RA, which can significantly decrease the risk of RA by 0.77-fold. Consistent with our findings, Knevel et al. [21] reported a decreased RA risk of 0.71-fold in subjects with allele C from the North American population. As for the other reported susceptible locus of RA, rs6435818, we found no association with RA in the Chinese population. Julia et al. [22] discovered that rs6435818 was associated with RA in the European population through GWAS, although the combination of data in discovery stage and replication stages did not reach genome-wide significance. Herein, the lack of replication could be probably attributed to the ethnic difference or relatively small sample size of the current study. However, the association of rs6435818 with RA was worthy of further replication in other populations.

To determine the potential functional role of rs7607479 in RA, for the first time, we analyzed the tissue expression of rs7607479-related genes including MMP-3 and SPAG16 in the synovial tissues of both RA patients and non-RA controls. Significantly higher tissue expression of MMP-3 was observed in RA patients than in the controls. Moreover, patients with genotype TT of rs7607479 were found to have significantly increased expression of MMP-3 than those with genotype CC. Interestingly, the expression of SPAG16 was not statistically different between the patients and the controls. Moreover, there was no significant difference regarding the expression of SPAG16 between patients with different genotypes of rs7607479. In line with our findings, Knevel et al. [21] observed that SPAG16 expression did not correlate with rs7607479 genotypes in fibroblast-like synoviocytes, while patients with the allele C of rs7607479 had lower MMP-3 expression in serum levels. Taken together, it was possible that rs7607479 might be involved in the development of RA via regulation of MMP-3, although it was located in the intronic region of SPAG16. Interestingly, in this study, we confirmed that MMP-3 is highly expressed in synovium tissues of ACPA-positive RA patients, and we observed significant correlation between MMP-3 expression and serum level of RF and ESR. Collectively, increased MMP-3 expression may indicate poor prognosis of the disease.

To date, the exact regulatory effect of rs7607479 on MMP-3 expression remains obscure. Since we found no significant difference regarding the expression of SPAG16 between patients and controls, it was not likely that the relationship between rs7607479 and MMP-3 is mediated through interaction between MMP-3 and SPAG16. Other unknown regulatory elements overlapped with the same genomic region of rs7607479 may play a role in this relationship. As speculated by Knevel et al. [21], miRNA encoded within the same intron of rs7607479 is worthy of further investigation to uncover the molecular mechanism underlying the regulation of MMP-3 production by rs7607479.

The limitations of our study should be addressed here. First, we did not perform in vivo experiments to validate the role of rs7607479 in MMP-3 expression. Further functional experiments are warranted to clarify the underlying regulatory mechanism. Second, the sample size of our study was smaller than previous GWAS study. Besides, no negative controls were included in the qPCR experiments, which could potentially lead to bias of the results. Although the power analysis confirmed that the sample size was sufficient, we believed that in future study larger sample size with negative controls may be included to further validate the relationship between SPAG16 and RA. Third, it cannot be ruled out that the expression of MMP-3 and SPAG16 may be changed due to injury. In future study, tissues from normal controls need to be collected to better investigate if the gene expression of MMP-3 and SPAG16 was changed in RA patients.

Conclusions

We validated a common genetic risk factor in ACPA-positive patients with RA, which is potentially associated with the tissue expression of MMP-3 and disease progression. Further functional analysis into the role of rs7607479 in MMP-3 expression can shed new light on the genetic architecture of ACPA-positive RA.

Availability of data and materials

All data used in this study are available at the request of editors, reviewers and the research community.

Abbreviations

RA:

Rheumatoid arthritis

GWASs:

Genome-wide association studies

BMI:

Body mass index

RF:

Rheumatoid factor

ESR:

Erythrocyte sedimentation rate

CRP:

C-reactive protein

HWE:

Hardy–Weinberg equilibrium

OR:

Odds ratio

Cis:

Confidential intervals

ACPA:

Antibodies against citrullinated proteins

MMP:

Matrix metalloproteinases

References

  1. Prieto-Potin I, Largo R, Roman-Blas JA, Herrero-Beaumont G, Walsh DA. Characterization of multinucleated giant cells in synovium and subchondral bone in knee osteoarthritis and rheumatoid arthritis. BMC Musculoskelet Disord. 2015;16:226.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Senolt L, Housa D, Vernerova Z, Jirasek T, Svobodova R, Veigl D, Anderlova K, Muller-Ladner U, Pavelka K, Haluzik M. Resistin in rheumatoid arthritis synovial tissue, synovial fluid and serum. Ann Rheum Dis. 2007;66(4):458–63.

    Article  CAS  PubMed  Google Scholar 

  3. Barik RR, Bhatt LK. Emerging epigenetic targets in rheumatoid arthritis. Rheumatol Int. 2021;41(12):2047–67.

    Article  PubMed  Google Scholar 

  4. Frank-Bertoncelj M, Klein K, Gay S. Interplay between genetic and epigenetic mechanisms in rheumatoid arthritis. Epigenomics. 2017;9(4):493–504.

    Article  CAS  PubMed  Google Scholar 

  5. Nemtsova MV, Zaletaev DV, Bure IV, Mikhaylenko DS, Kuznetsova EB, Alekseeva EA, Beloukhova MI, Deviatkin AA, Lukashev AN, Zamyatnin AA Jr. Epigenetic changes in the pathogenesis of rheumatoid arthritis. Front Genet. 2019;10:570.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Kawai VK, Shi M, Feng Q, Chung CP, Liu G, Cox NJ, Jarvik GP, Lee MTM, Hebbring SJ, Harley JB, Kaufman KM, Namjou B, Larson E, Gordon AS, Roden DM, Stein CM, Mosley JD, e MI. Pleiotropy in the genetic predisposition to rheumatoid arthritis: a phenome-wide association study and inverse variance-weighted meta-analysis. Arthritis Rheumatol. 2020,72(9):1483–1492.

  7. Knevel R, Grondal G, Huizinga TW, Visser AW, Jonsson H, Vikingsson A, Geirsson AJ, Steinsson K, van der Helm-van Mil AH. Genetic predisposition of the severity of joint destruction in rheumatoid arthritis: a population-based study. Ann Rheum Dis. 2012;71(5):707–9.

    Article  PubMed  Google Scholar 

  8. Vetchinkina EA, Mikhaylenko DS, Kuznetsova EB, Deryagina TA, Alekseeva EA, Bure IV, Zamyatnin AA, Jr., Nemtsova MV. Genetic factors of predisposition and clinical characteristics of rheumatoid arthritis in Russian patients. J Pers Med,2021,11(6).

  9. Myasoedova E, Davis J, Matteson EL, Crowson CS. Is the epidemiology of rheumatoid arthritis changing? Results from a population-based incidence study, 1985–2014. Ann Rheum Dis. 2020;79(4):440–4.

    Article  PubMed  Google Scholar 

  10. van der Woude D, van der Helm-van Mil AHM. Update on the epidemiology, risk factors, and disease outcomes of rheumatoid arthritis. Best Pract Res Clin Rheumatol. 2018;32(2):174–87.

    Article  PubMed  Google Scholar 

  11. Prahalad S, Ryan MH, Shear ES, Thompson SD, Glass DN, Giannini EH. Twins concordant for juvenile rheumatoid arthritis. Arthritis Rheum. 2000;43(11):2611–2.

    Article  CAS  PubMed  Google Scholar 

  12. Hayashi S, Matsubara T, Fukuda K, Maeda T, Funahashi K, Hashimoto M, Kamenaga T, Takashima Y, Kuroda R. A genome-wide association study identifying the SNPs predictive of rapid joint destruction in patients with rheumatoid arthritis. Biomed Rep. 2021;14(3):31.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Kwon YC, Lim J, Bang SY, Ha E, Hwang MY, Yoon K, Choe JY, Yoo DH, Lee SS, Lee J, Chung WT, Kim TH, Sung YK, Shim SC, Choi CB, Jun JB, Kang YM, Shin JM, Lee YK, Cho SK, Kim BJ, Lee HS, Kim K, Bae SC. Genome-wide association study in a Korean population identifies six novel susceptibility loci for rheumatoid arthritis. Ann Rheum Dis. 2020;79(11):1438–45.

    Article  CAS  PubMed  Google Scholar 

  14. Sundbaum JK, Baecklund E, Eriksson N, Kohnke H, Wallenberg M, Cavalli M, Wadelius C, Wadelius M, Hallberg P. Genome-wide association study of liver enzyme elevation in rheumatoid arthritis patients starting methotrexate. Pharmacogenomics. 2021;22(15):973–82.

    Article  CAS  PubMed  Google Scholar 

  15. Wasserman AM. Diagnosis and management of rheumatoid arthritis. Am Fam Physician. 2011;84(11):1245–52.

    PubMed  Google Scholar 

  16. Navarro-Millan I, Darrah E, Westfall AO, Mikuls TR, Reynolds RJ, Danila MI, Curtis JR, Investigators C, Rosen A, Bridges SL Jr. Association of anti-peptidyl arginine deiminase antibodies with radiographic severity of rheumatoid arthritis in African Americans. Arthritis Res Ther. 2016;18(1):241.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Nowakowska-Plaza A, Potaczek DP, Gluszko P, Undas A. Antibodies to N-homocysteinylated albumin and haemoglobin in patients with rheumatoid arthritis: a potential new marker of disease severity. Scand J Rheumatol. 2014;43(1):17–21.

    Article  CAS  PubMed  Google Scholar 

  18. Daha NA, Toes RE. Rheumatoid arthritis: are ACPA-positive and ACPA-negative RA the same disease? Nat Rev Rheumatol. 2011;7(4):202–3.

    Article  PubMed  Google Scholar 

  19. Grosse J, Allado E, Roux C, Pierreisnard A, Couderc M, Clerc-Urmes I, Remen T, Albuisson E, De Carvalho-Bittencourt M, Chary-Valckenaere I, Loeuille D. ACPA-positive versus ACPA-negative rheumatoid arthritis: two distinct erosive disease entities on radiography and ultrasonography. Rheumatol Int. 2020;40(4):615–24.

    Article  PubMed  Google Scholar 

  20. Burgers LE, van Steenbergen HW, Ten Brinck RM, Huizinga TW, van der Helm-van Mil AH. Differences in the symptomatic phase preceding ACPA-positive and ACPA-negative RA: a longitudinal study in arthralgia during progression to clinical arthritis. Ann Rheum Dis. 2017;76(10):1751–4.

    Article  CAS  PubMed  Google Scholar 

  21. Knevel R, Klein K, Somers K, Ospelt C, Houwing-Duistermaat JJ, van Nies JA, de Rooy DP, de Bock L, Kurreeman FA, Schonkeren J, Stoeken-Rijsbergen G, Helmer Q, van der Linden MP, Kern M, Manjarrez-Orduno N, Rodriguez-Rodriquez L, Stinissen P, Huizinga TW, Toes RE, Gay S, Gregersen PK, Somers V, van der Helm-van Mil AH. Identification of a genetic variant for joint damage progression in autoantibody-positive rheumatoid arthritis. Ann Rheum Dis. 2014;73(11):2038–46.

    Article  CAS  PubMed  Google Scholar 

  22. Julia A, Gonzalez I, Fernandez-Nebro A, Blanco F, Rodriguez L, Gonzalez A, Canete JD, Maymo J, Alperi-Lopez M, Olive A, Corominas H, Martinez-Taboada V, Erra A, Sanchez-Fernandez S, Alonso A, Lopez-Lasanta M, Tortosa R, Codo L, Gelpi JL, Garcia-Montero AC, Bertranpetit J, Absher D, Bridges SL Jr, Myers RM, Tornero J, Marsal S. A genome-wide association study identifies SLC8A3 as a susceptibility locus for ACPA-positive rheumatoid arthritis. Rheumatology (Oxford). 2016;55(6):1106–11.

    Article  CAS  PubMed  Google Scholar 

  23. Alciaturi J, Anesetti G, Irigoin F, Skowronek F, Sapiro R. Distribution of sperm antigen 6 (SPAG6) and 16 (SPAG16) in mouse ciliated and non-ciliated tissues. J Mol Histol. 2019;50(3):189–202.

    Article  CAS  PubMed  Google Scholar 

  24. Castagna A, Cesari E, Garofalo R, Gigante A, Conti M, Markopoulos N, Maffulli N. Matrix metalloproteases and their inhibitors are altered in torn rotator cuff tendons, but also in the macroscopically and histologically intact portion of those tendons. Muscles Ligaments Tendons J. 2013;3(3):132–8.

    PubMed Central  PubMed  Google Scholar 

  25. Del Buono A, Oliva F, Osti L, Maffulli N. Metalloproteases and tendinopathy. Muscles Ligaments Tendons J. 2013;3(1):51–7.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Diniz-Fernandes T, Godoy-Santos AL, Santos MC, Pontin P, Pereira CAA, Jardim YJ, Velosa APP, Maffulli N, Teodoro WR, Capelozzi VL. Matrix metalloproteinase-1 (MMP-1) and (MMP-8) gene polymorphisms promote increase and remodeling of the collagen III and V in posterior tibial tendinopathy. Histol Histopathol. 2018;33(9):929–36.

    CAS  PubMed  Google Scholar 

  27. Magra M, Maffulli N. Matrix metalloproteases: a role in overuse tendinopathies. Br J Sports Med. 2005;39(11):789–91.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Yong W, Hongbin L, Jing W, Jing Z, Ning T, Lijie B. Associations of changes in serum inflammatory factors, MMP-3, 25(OH)D and intestinal flora with osteoporosis and disease activity in rheumatoid arthritis patients. Clin Lab. 2020,66(12).

  29. Huang J, Xie B, Li Q, Xie X, Zhu S, Wang M, Peng W, Gu J. Infliximab reduces CD147, MMP-3, and MMP-9 expression in peripheral blood monocytes in patients with active rheumatoid arthritis. Eur J Pharmacol. 2013;698(1–3):429–34.

    Article  CAS  PubMed  Google Scholar 

  30. Siebuhr AS, Kjelgaard-Petersen CF, Sun S, Byrjalsen I, Christiansen C, Karsdal MA, Bay-Jensen AC. Suppression of active, but not total MMP-3, is associated with treatment response in a phase III clinical study of rheumatoid arthritis. Clin Exp Rheumatol. 2018;36(1):94–101.

    PubMed  Google Scholar 

  31. Nemec P, Pavkova-Goldbergova M, Gatterova J, Vasku A, Soucek M. Association of the 5A/6A promoter polymorphism of the MMP-3 gene with the radiographic progression of rheumatoid arthritis. Ann N Y Acad Sci. 2007;1110:166–76.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all the subjects for their consent to be included in this research.

Funding

No funding was used in this research.

Author information

Authors and Affiliations

Authors

Contributions

LQ and ZB collected the data and drafted the manuscript. SX, YW, LQ and ST performed the laboratory experiments. CC, LY and WX reviewed and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yetian Li or Xing Wei.

Ethics declarations

Ethical approval and consent to participate

Written consent was obtained from each subject. The histological study of surgical samples was approved by the Institutional Ethics Committee. All experimental protocols and methods were carried out in accordance with relevant guidelines and regulations and complied with the Helsinki Declaration.

Consent for publication

All the authors agreed to publish this paper.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Supplementary table 1. Post-hoc analysis for the comparison of genotype frequency of rs7607479 between cases and controls.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lin, Q., Zhou, B., Song, X. et al. Genetic variant in SPAG16 is associated with the susceptibility of ACPA-positive rheumatoid arthritis possibly via regulation of MMP-3. J Orthop Surg Res 17, 511 (2022). https://doi.org/10.1186/s13018-022-03405-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-022-03405-w

Keywords