Investigation of MMP-1 genetic polymorphisms and protein expression and their effects on the risk of Kashin-Beck disease in the northwest Chinese Han population
© The Author(s). 2016
Received: 19 March 2016
Accepted: 19 May 2016
Published: 31 May 2016
The etiology of Kashin-Beck disease (KBD), an endemic osteochondropathy, is largely unknown. Matrix metalloproteinase-1 (MMP-1) plays a central role in the initiation and progression of cartilage destruction; however, no study has reported on the relationship between KBD and MMP-1. This study was to investigate the role of MMP-1 in the pathogenesis and progression of KBD.
Single nucleotide polymorphism (SNP) genotyping was conducted for 274 KBD cases and 248 healthy controls using the Sequenom MassARRAY system. Additionally, the expression of MMP-1 in the knee articular cartilage of 22 KBD patients and 21 controls was analyzed by immunohistochemistry, and the concentration of MMP-1 in their joint fluid was also measured by enzyme-linked immunosorbent assay (ELISA).
The results showed that two SNPs (rs470221 and rs1144396) had a weak association with increased KBD risk; however, the significance of these results did not survive Bonferroni’s correction. Moreover, the percentages of cells expressing MMP-1 in each layer of cartilage were significantly higher in the KBD group than in the controls (F = 11.41–28.31, P = 0.002–0.000). The concentration of MMP-1 in KBD joint fluid was significantly higher than that in the controls (t = 9.83, P < 0.0001).
The increased expression of MMP-1 has a potential effect on the risk of KBD in the northwest Chinese Han population. However, six selected SNPs in the MMP-1 gene might not be useful as significant markers for predicting KBD susceptibility in Chinese Han population. Therefore, future studies in the association of MMP-1 with KBD should focus on other candidate SNPs.
KeywordsKashin-Beck disease Matrix metalloproteinase-1 Single nucleotide polymorphisms Enzyme-linked immunosorbent assay Immunohistochemistry
Kashin-Beck disease (KBD) is a chronic osteochondropathy affecting the bones and joints that is endemic to certain geographical areas. According to the 2010 Health Statistical Year Book of China, the KBD-affected area in China covers 366 counties of 14 provinces or autonomous regions encompassing more than 690,000 KBD patients and a population of 10.584 million individuals at high risk of developing the disease. A key pathological feature of KBD is chondrocyte necrosis in the deep zone of the growth plate of cartilage and articular cartilage [1, 2]. Clinically, the disease usually presents in childhood, between 5 and 13 years of age, and mainly attacks the growth plate cartilage. KBD presents as dwarfism, very short upper limbs, and deformed and painful joints. The most frequently affected sites are the distal limb joints (e.g., ankles, knees, wrists, and elbows). Following cartilage damage, atrophied muscles occur due to joint pain and limited mobility; these deformed joints make manual farming work more difficult and painful.
The etiology of KBD remains unclear; over the past 150 years, three environmental hypotheses have been proposed: selenium deficiency, mycotoxins from contaminated storage grains, and organic matter (e.g., fulvic acid or FA) in drinking water [3–5]. However, these theories cannot account for the variable differences in the occurrence of KBD among residents exposed to the same environmental risk factors . In our previous study, we found that the incidence of KBD occurred in familial clusters and that 35.10 % of cases could be attributed to the polymorphic variations of KBD . Certain susceptibility genes may affect susceptibility to environmental factors, such as selenium deficiency or other biologic factors [8, 9]. The role of matrix metalloproteinases (MMPs) in the degradative events invoked in the cartilage and bone of arthritic joints has long been appreciated [10, 11]. In a previous association analysis, we indicated that gene markers near locus D11S4094 should be further studied . The gene coding for MMP-1 (11q22-q23) is near D11S4094 (11q23); studies on the MMP family have shown that MMP-1 may play a central role in the initiation and progression of cartilage destruction [12, 13]. However, to the best of our knowledge, no study has reported on the relationship between KBD and MMP-1. In this study, for the first time, we evaluated whether MMP-1 has a potential effect on the risk and progression of KBD.
In total, 522 unrelated Chinese Han individuals were included in the current study; these individuals were from KBD-endemic areas of the Linyou and Yongshou counties of Shaanxi province, in northwest China. This group consisted of 274 KBD patients and 248 healthy controls. Radiographs of the right hand were taken for both the KBD patients and the healthy controls and read by veteran orthopedists. KBD was diagnosed according to the national diagnostic criteria of China (WS/T 207-2010), which are based on radiographic examination (i.e., an X-ray of the right hand) and clinical diagnosis (from degree I to III). Patients with clinical symptoms or radiographic changes associated with other types of osteochondropathy were excluded. A healthy case was defined as no KBD and no primary or secondary osteoarthritis (OA). The controls were randomly selected and were frequency-matched by age and sex. Fresh blood (5 mL) was collected from the antecubital vein of all 522 subjects while in a fasting state.
Knee articular cartilage and synovial joint fluids were collected from 22 KBD patients and 21 healthy controls. The KBD patients, consisting of ten males and 12 females with an average age of 51.00 ± 8.30 (32–66) years, underwent knee debridement or arthroplasty at a hospital. The healthy control subjects, consisting of 11 males and ten females with an average age of 48.23 ± 7.65 (33–61) years, had no history of osteochondropathy, but they had undergone amputation because of traffic accidents. No significant differences were observed between KBD and control group in age (t = 1.13, P > 0.05) and sex (x 2 = 0.21, P > 0.05). None of the KBD patients or controls was diagnosed with bone or cartilage genetic diseases or rheumatoid arthritis (RA).
The study was performed in accordance with the Declaration of Helsinki and approved by the Human Ethics Committee of Xi’an Jiaotong University, People’s Republic of China. Written informed consent was also obtained from the subjects or their relatives.
SNPs selection and genotyping
The loci information of the six SNPs in MMP-1
Minor allele frequency (%)
HWE test (P)
The cartilage tissues were immediately fixed in 4 % (w/v) paraformaldehyde following collection, washed in phosphate-buffered saline (PBS), decalcified, embedded in paraffin, and cut into 5–8-μm-thick slices for immunohistochemistry and hematoxylin and eosin (HE) staining. Immunochemical identification was performed using the streptavidin-peroxidase (SP) method. Briefly, after deparaffinization, endogenous peroxidase was blocked with 3 % H2O2 for 15 min, after which the slides were washed with PBS. The slides were then predigested using a digestive complex followed by rinsing with PBS. After blocking using 10 % normal goat serum, the sections were incubated with a primary antibody recognizing MMP-1 at 1:100 dilution (polyclonal rabbit anti-MMP-1, Bioss Co, Beijing, China) or with PBS, serving as a negative control, at 4 °C overnight. Next, the sections were incubated with 1:200 biotinylated goat anti-rabbit IgG (ZSGB-Bio Co, Beijing ,China) at 37 °C for 20 min, followed by incubation with horseradish peroxidase-labeled streptavidin solution at 37 °C for 15 min. Color development was continued for 5 min at room temperature using diaminobenzidine followed by rinsing with distilled water. Counterstaining was performed with hematoxylin.
The percentage of positive cells was evaluated by counting within 10 high-magnification power fields (40×) in six consecutive tissue sections. The level of expression was determined using the following semi-quantitative criteria: 0 = negative; 1+ = less than 25 % positive staining; 2+ = 25–50 % positive staining; 3+ = more than 50 % positive staining .
Enzyme-linked immunosorbent assay measurement of MMP-1 in joint fluid
The joint fluid samples were dispensed into 1-ml aliquots that were treated with 50 ug/ml of hyaluronidase for 1 h at room temperature. The joint fluid was then diluted with diluent’s buffer to the appropriate detection range for enzyme-linked immunosorbent assay (ELISA). The MMP-1 level was measured using a commercial ELISA kit (R&D System, Inc 614 McKinley Place NE, Minneapolis, MN 55413, USA) in accordance with the manufacturer’s manual.
The Hardy-Weinberg equilibrium (HWE) of each SNP was tested by the goodness-of-fit x 2 test to compare the expected frequencies of genotypes in controls, SNPs with P > 0.05 were considered to be in HWE . Differences in genotypes and allele frequencies between the KBD cases and the controls were determined using CLUMP22 software with 10,000 Monte Carlo simulations. Unconditional logistic regression analysis adjusted for age and gender was used to estimate the strength of the association through calculation of odds ratios (ORs) with their 95 % confidence intervals (95 % CIs) . Dominant and recessive genetic models were also evaluated and expressed as ORs with 95 % CI . To account for multiple testing, Bonferroni’s correction was applied, consisting of dividing the usual significance level of 5 % by the number of tests . Differences in MMP-1 expression in the cartilage between the two groups were examined by analysis of variance (ANOVA) ; the independent-sample t test was used to compare the concentrations of MMP-1 in joint fluid of the two groups. P < 0.05 was considered to indicate statistically significant differences.
A total of 274 KBD patients and 248 age and sex matched controls were included in the genetic polymorphism study. No significant differences were observed between KBD and control group in age (53.37 ± 10.79 vs 51.71 ± 17.85, t = 1.29, P > 0.05) and sex (male/female, 125/149 vs 124/124, x 2 = 1, P > 0.05). In addition, the clinical stages of the 274 KBD patients were divided into I, II, and III stages, and the constituent ratio was 57.66 % (158/274), 31.02 % (85/274), and 11.32 % (31/274), respectively.
Association of polymorphisms of MMP-1 with KBD susceptibility
Analysis of association of the six SNPs gene polymorphism with the risk of KBD
KBD (N = 274)
Controls (N = 248)
OR (95 % CI)
OR (95 % CI)
Expression of MMP-1 in articular cartilage and levels of MMP-1 in joint fluids of two groups
Immunohistochemical expression of MMP-1 in different layers of the articular cartilage in the two groups
Level of expression
Positive rate (%)
Mean ± SD
Control (n = 21)
45.27 ± 20.48
KBD (n = 22)
67.02 ± 16.19*
Control (n = 21)
35.71 ± 24.80
KBD (n = 22)
58.19 ± 18.01*
Control (n = 21)
3.00 ± 32.44
KBD (n = 22)
34.87 ± 27.23*
Control (n = 21)
27.92 ± 13.59
KBD (n = 22)
53.32 ± 14.58*
KBD is an endemic OA with a specific geographical distribution. However, not all families or individuals in endemic areas suffer from KBD. Additionally, recent studies reported that susceptibility genes might contribute to the pathogenesis of KBD and account for an individual’s susceptibility to KBD . MMP-1 is the most highly expressed interstitial collagenase that degrades fibrillar collagens, which are major constituents of the extracellular matrix (ECM) . Many studies have provided evidence for an association between the MMP-1 gene and arthritis [22, 23]. KBD has been found to have overlapping phenotypes and pathologic changes similar to those of OA and RA . The polymorphism study of KBD is very limited; a similar study showed a significant association between rs6910140 of COL9A1 and KBD . Used genome-wide association study (GWAS), ITPR2, was identified as a susceptibility gene for KBD . However, none of the six selected SNPs in the MMP-1 gene showed an association with KBD in our study; the finding of this study was not consistent with results of previous studies, and we think some limitations should be considered in relation to present results: First, it is common for a genetic variant to be associated with certain populations but not with others , and we failed to investigate the role of MMP-1 in the pathogenesis and progression of KBD in different populations. Second, the number of SNPs was very limited in this study; the susceptibility of KBD may be associated with other SNPs in the MMP-1 gene. So we think it should not dismiss a possible association between MMP-1 gene polymorphism and KBD; the gene’s roles in determining the expression of MMP-1 should be further studied, and an approach examining the common genetic variation of the gene should be used.
MMPs are involved in joint destruction in arthritis and are strongly associated with the levels of inflammation. MMP-1 is more abundant and also degrades collagens effectively. In the present study, the expression of MMP-1 in articular cartilage of knee was analyzed by immunohistochemistry; the negative control showed negative staining, which indicated that the test was not confounded by other factors and was reasonable. Our results showed increased expression of MMP-1 should play a role in the pathogenesis of KBD. These results are similar to those found in other degenerative joint diseases, such as OA and RA. Regarding KBD, a similar study has been reported, in which human chondrocytes were isolated and cultured on bone matrix gelatin to form an artificial cartilage model in vitro with or without T-2 toxin and selenium. With exposure to the KBD risk factors (with T-2 toxin and without selenium), Western blotting and RT-PCR analyses showed level of MMP-1 was increased significantly after exposure of T-2 toxin; increasing doses of T-2 toxin resulted in an increased expression of MMP-1; selenium increased TIMP-1 expression and decreased MMP-1 expression and partly blocked the effects of T-2 toxin on the balance between TIMP-1 and MMPs . More specifically, studies on the expression of MMP-1 in patients with RA or OA found that the circulating levels of MMP-1 were elevated. In our study, the levels of MMP-1 in joint fluid of knee were also detected and the data showed higher level of MMP-1 in patients with KBD than in controls. A previous study in the cases of the same KBD area showed that the serum MMP-1 levels in KBD cases were higher than those of controls, even though no significant difference was detected . All these studies showed that elevated levels of MMP-1 in cartilage and synovium may be a cause of degrading cartilage in the pathogenesis of KBD in Chinese Han population.
Tight regulation of MMP activity is critical for maintaining a fine balance between aggrecan and collagen anabolism and catabolism. The MMP-1 enzyme is known as collagenase or collagenase-1, and it degrades a broad range of substrates (e.g., collagens I, II, III, VII, VIII, and X; gelatin; aggrecan; versican; and proteoglycan link protein), leading to collagen network edema that can be inhibited by TIMP-1. Articular cartilage contains a large number of matrix macromolecules, such as proteoglycan, aggrecan, and type II collagen [28, 29], which provide elastic and strong tensile properties that resist mechanical loading. Under normal conditions, there is a balance between catabolism and anabolism in the ECM in articular cartilage [30, 31]. However, certain conditions, such as inflamed synovium and pannus tissue, may disturb this balance and accelerate catabolism of the ECM, leading to loss of ECM components from the cartilage, thereby severely impairing ECM functionality and causing pain and disability . Studies have reported that IL-1β and INF-α expression in the synovium, and their levels in synovial fluid and serum are significantly higher in patients with KBD than in controls [1, 33]. We think that MMP-1 overexpression induced by IL-1β, INF-α, and other inflammation factors may play an important role in destroying the articular cartilage of KBD patients.
However, there are also some limitations to this study: the sample populations were a little small and only MMP-1 was included in this study. And it has been reported that selenium has a protective role while T-2 toxin induces cartilage matrix degradation by regulation of MMP-1 /TIMP-1 expression; however, environmental factors such as low Se and T-2 toxin were not analyzed in this study. Therefore, the effect of gene-environment interactions needs further investigation in future well-designed high-quality studies among well-subgrouped populations.
In conclusion, up-regulation of MMP-1 in the ECM accelerates ECM catabolism, leading to matrix degradation, cartilage degeneration, and subsequent pathogenesis and progression of KBD. However, the six selected SNPs in the MMP-1 gene might not be useful as significant markers for predicting KBD susceptibility in the northwest Chinese Han population. Therefore, analysis of the genetic component of KBD needs further investigation.
This research was supported by grants from the National Natural Science Foundation of China (No. 81273007) and The Foundation of Shaanxi Province of China (2011K15-06-05). For the sample collection, the authors acknowledge the help from the Centers for Disease Control and Prevention of Linyou and Yongshou counties.
XW, SHI, and XG designed the study. XW, SHI, ALL, J Ma,YW, ZT, and ZH performed the experiments. ALL and F Zh performed the statistical analysis. All authors analyzed and interpreted the results. XW, SHI, and F Zh drafted the manuscript. XG critically revised the manuscript for important intellectual content. All authors attest to the integrity of the original data and the analysis reported in this manuscript. The final manuscript was approved by all authors, who all qualify for authorship.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Cao J, Li S, Shi Z, Yue Y, Sun J, Chen J, et al. Articular cartilage metabolism in patients with Kashine-Beck disease: an endemic osteoarthropathy in China. Osteoarthr Cartil. 2008;16:680–8.View ArticlePubMedGoogle Scholar
- Mathieu F, Begaux F, Lan ZY, Suetens C, Hinsenkamp M. Clinical manifestations of Kashin-Beck disease in Nyemo Valley, Tibet. Int Orthop. 1997;21:151–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Moreno R. Kashin-Beck osteoarthropathy in rural Tibet in relation to selenium and iodine status. N Engl J Med. 1998;339:1112–20.View ArticleGoogle Scholar
- Chasseur C, Suetens V, Michel V, Mathieu F, Begaux F, Nolard N, et al. A 4-year study of the mycological aspects of Kashin-Beck disease in Tibet. Int Orthop. 2001;25:154–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Sudre P, Mathieu F. Kashin-Beck disease: from etiology to prevention or from prevention to etiology? Int Orthop. 2001;25:175–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Suete C, Moreno R, Chassear C, Mathieu F, Begaux F, Haubruge E, et al. Epidemiological support for a multifactorial aetiology of Kashin-Beck disease in Tibet. Int Orthop (SICOT). 2001;25:180–7.View ArticleGoogle Scholar
- Shi XW, Guo X, Lv AL, Kang L, Zhou YL, Zhang YZ, et al. Heritability estimates and linkage analysis of 23 short tandem repeat loci on chromosomes 2, 11, and 12 in an endemic osteochondropathy in China. Scand J Rheumatol. 2010;39:259–65.View ArticlePubMedGoogle Scholar
- Zhang F, Guo X, Wang WZ, Yan H, Li C. Genome-wide gene expression analysis suggests an important role of hypoxia in the pathogenesis of endemic osteochondropathy Kashin-Beck disease. PLoS One. 2011;6:e22983.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi XW, Guo X, Ren FL, Li J, Wu XM. The effect of short tandem repeat loci and low selenium levels on endemic osteoarthritis in China. J Bone Joint Surg Am. 2010;92:72–80.View ArticlePubMedGoogle Scholar
- Murphy G, Knäuper V, Atkinson S, Butler G, English W, Hutton M, et al. Matrix metalloproteinases in arthritic disease. Arthritis Res. 2002;4:S39–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Mengshol JA, Mix KS, Brinckerhoff CE. Matrix metalloproteinases as therapeutic targets in arthritic diseases: bull’s-eye or missing the mark? Arthritis Rheum. 2002;46:13–20.View ArticlePubMedGoogle Scholar
- Du F, Lü LJ, Teng JL, Shen N, Ye P, Bao CD. T-614 alters the production of matrix metalloproteinases (MMP-1 andMMP-3) and inhibits the migratory expansion of rheumatoid synovial fibroblasts, in vitro. Int Immunopharmacol. 2012;13:54–60.View ArticlePubMedGoogle Scholar
- Yoshihara Y, Nakamura H, Obata K, Yamada H, Hayakawa T, Fujikawa K, et al. Matrix metalloproteinases and tissue inhibitors of metalloproteinases in synovial fluids from patients with rheumatoid arthritis or osteoarthritis. Ann Rheum Dis. 2000;9:455–61.View ArticleGoogle Scholar
- Sun T, Gao Y, Tan W, Zhang X, Wang Y, Zhang Q, et al. Haplotypes in matrix metalloproteinase gene cluster on chromosome 11q22 contribute to the risk of lung cancer development and progression. Clin Cancer Res. 2006;12:7009–17.View ArticlePubMedGoogle Scholar
- Shi XW, Zhang F, Lv A, Wen Y, Guo X. COL9A1 Gene polymorphism is associated with Kashin-Beck disease in a northwest Chinese Han population. PLoS One. 2015;10:e01203650.Google Scholar
- Zhou J, Huang Y, Huang RS, Wang F, Xu L, Le Y, et al. A case-control study provides evidence of association for a common SNP rs974819 in PDGFD to coronary heart disease and suggests a sex-dependent effect. Thromb Res. 2012;130:602–6.View ArticlePubMedGoogle Scholar
- Rouault K, Scotet V, Autret S, Gaucher F, Dubrana F, Tanguy D, et al. Do HOXB9 and COL1A1 genes play a role in congenital dislocation of the hip? Study in a Caucasian population. Osteoarthritis Cartilage. 2009;17:1099–105.View ArticlePubMedGoogle Scholar
- Yuan LJ, Jin TB, Yin JK, Du XL, Wang Q, Dong R, et al. Polymorphisms of tumor-related genes IL-10, PSCA, MTRR and NOC3L are associated with the risk of gastric cancer in the Chinese Han population. Cancer Epidemiol. 2012;36:366–72.View ArticleGoogle Scholar
- Ren FL, Guo X, Zhang RJ, Wang SJ, Zuo H, Zhang ZT, et al. Effects of selenium and iodine deficiency on bone, cartilage growth plate and chondrocyte differentiation in two generations of rats. Osteoarthritis Cartilage. 2007;15:1171–7.View ArticlePubMedGoogle Scholar
- Xiong YM, Mo XY, Zou XZ, Song RX, Sun WY, Lu W, et al. Association study between polymorphisms in selenoprotein genes and susceptibility to Kashin-Beck disease. Osteoarthritis Cartilage. 2010;18:817–24.View ArticlePubMedGoogle Scholar
- Sauter W, Rosenberger A, Beckmann L, Kropp S, Mittelstrass K, Timofeeva M, et al. Matrix metalloproteinase 1 (MMP1) is associated with early-onset lung cancer. Cancer Epidemiol Biomarkers Prev. 2008;17:1127–35.View ArticlePubMedGoogle Scholar
- Dorr S, Lechtenbohmer N, Rau R, Herborn G, Wagner U, Müller-Myhsok B, et al. Association of a specific haplotype across the genes MMP1 and MMP3 with radiographic joint destruction in rheumatoid arthritis. Arthritis Res Ther. 2004;6:199–207.View ArticleGoogle Scholar
- Barlas IO, Sezgin M, Erdal ME, Sahin G, Ankarali HC, Altintas ZM, et al. Association of (-1,607) 1G/2G polymorphism of matrix metalloproteinase-1 gene with knee osteoarthritis in the Turkish population (knee osteoarthritis and MMPs gene polymorphisms). Rheumatol Int. 2009;29:383–8.View ArticlePubMedGoogle Scholar
- Shi Y, Lu F, Liu X, Wang Y, Huang L, Liu X, et al. Genetic variants in the HLA-DRB1 gene are associated with Kashin-Beck disease in the Tibetan population. Arthritis Rheum. 2011;63:3408–16.View ArticlePubMedGoogle Scholar
- Zhang F, Wen Y, Guo X, Zhang Y, Wang X. Genome-wide association study identifies ITPR2 as a susceptibility gene for Kashin-Beck disease in Han Chinese. Arthritis Rheum. 2015;67:176–81.View ArticleGoogle Scholar
- Seki S, Kawaguchi Y, Mori M, Mio F, Chiba K, Mikami Y, et al. Association study of COL9A2 with lumbar disc disease in the Japanese population. J Hum Genet. 2006;51:1063–7.View ArticlePubMedGoogle Scholar
- Chen J, Chu Y, Cao J, Wang W, Liu J, Wang J. Effects of T-2 toxin and selenium on chondrocyte expression of matrix metalloproteinases (MMP-1, MMP-13), α2- macroglobulin (α2M) and TIMPs. Toxicol In Vitro. 2011;25:492–9.View ArticlePubMedGoogle Scholar
- Giannelli G, Erriquez R, Iannone F, Marinosci F, Lapadula G, Antonaci S. MMP-2, MMP-9, TIMP-1 and TIMP-2 levels in patients with rheumatoid arthritis and psoriatic arthritis. Clin Exp Rheumatol. 2004;22:335–8.PubMedGoogle Scholar
- Sternlicht MD, Werb Z. How matrix metalloproteinase regulate cell behavior. Annu Rev Cell Dev Biol. 2001;17:463–516.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Mattey DL. Age at onset of rheumatoid arthritis: association with polymorphisms in the vascular endothelial growth factor A(VEGFA) gene and an intergenic locus between matrix metalloproteinase (MMP) 1 and 3 genes. Clin Exp Rheumatol. 2012;30:894–8.PubMedGoogle Scholar
- Zhang A, Cao JL, Yang B, Chen JH, Zhang ZT, Li SY, et al. Effects of moniliformin and selenium on human articular cartilage metabolism and their potential relationships to the pathogenesis of Kashin-Beck disease. J Zhejiang Univ Sci B. 2010;11:200–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Cunnane G, Fitzgerald O, Beeton C, Cawston TE, Bresnihan B. Early joint erosions and serum levels of matrix metalloproteinase 1, matrix metalloproteinase 3, and tissue inhibitor of metalloproteinases-1 in rheumatoid arthritis. Arthritis Rheum. 2001;44:2263–74.View ArticlePubMedGoogle Scholar
- Sun ZM, Ling M, Liu M, Zhang YG. Expression of interleukin-1β and tumor necrosis factor-α in the synovium and synovial fluid. J South Med Univ. 2009;29:5–8.Google Scholar