- Research article
- Open Access
- Open Peer Review
Association of rs1137101 polymorphism in LEPR and susceptibility to knee osteoarthritis in a Northwest Chinese Han population
BMC Musculoskeletal Disordersvolume 17, Article number: 311 (2016)
Osteoarthritis (OA) is a complex arthritic condition in which the genetic factor plays a crucial role. A single nucleotide polymorphism (SNP), rs1137101 (Gln223Arg) of leptin receptor (LEPR) gene has been demonstrated to be associated with susceptibility to knee OA. However, this association in Chinese Han population has never been examined. The present study aimed to investigate whether Gln223Arg was related to knee OA susceptibility in a Northwest Chinese population with Han ethnicity.
Gln223Arg polymorphisms were genotyped in 587 patients with confirmed knee OA and in 628 age- and sex-matched healthy controls using polymerase chain reaction-restriction fragment length polymorphism analysis. Besides, LEPR genotypes were verified by direct DNA sequencing analysis on PCR products.
The genotype and allele frequencies in LEPR SNP rs1137101 were significantly different between cases and control groups (chi-square = 6.52, P = 0.038 for genotype and chi-square = 5.06, P = 0.024 for allele frequencies; respectively). Rs1137101 was correlated with knee OA in the dominant genetic model (GG + GA versus AA) (P = 0.016) and a higher G allele frequency existed (P = 0.024) among all patients with knee OA and controls. On stratification analysis, the genotype GG and G allele were associated with susceptibility to knee OA in females, both young (≤65 years) and old groups (>65 years) and patients with mild knee OA.
Our finding suggested that the genetic variant of LEPR gene rs1137101 is independently related to knee OA susceptibility in Northwest Chinese population with Han ethnicity and may serve as a potential biomarker to determine risk of knee OA.
Osteoarthritis (OA) is a common form of progressive joint disorder characterized by age-related regressive changes in articular cartilage, with joint space narrowing, osteophyte formation and subcondral sclerosis [1, 2]. Common symptoms of OA consists of severe joint pain, stiffness, swelling, reduced motion, and a reduced quality of life. Among different joints part in the human body, the knee is the most common prevalent site that primary OA might be involved. A recent survey reported the morbidity of primary knee OA in people over 60 years is 19.4 % in China . Although the exact causes of this disease were still not well-defined till now, it has been demonstrated to be associated with several attributing factors including a variety of non-genetic factors (obesity, sex, age, and joint injury) and a combinatorial effect of genetic factors [4–6]. Therefore, there is an attractive prospect for us to find a susceptibility gene of knee OA, which could not only bring us a better understanding of the molecular mechanism of the disease, but also benefit for the development of OA targets drugs.
Leptin, a endogenous hormone expressed by the obese gene mainly in adipocytes, functions to control the metabolism of glucose and lipid through binding to the leptin receptor (LEPR, also known as Ob-R) . Dumond et al.  firstly reported the expression of leptin in articular cartilage and synovial fluid (SF), with only few leptin expression in normal cartilage chondrocytes, while in osteoarthritic cartilage and the areas of osteophytes formation, the expression level of leptin was remarkably increased and closely related to the severity of cartilage damage. Simopoulou et al.  demonstrated that leptin could induce the expression of metalloproteinases 9 and 13, 2 major enzymes that mediate the destructive process in OA and production of inflammatory cytokine IL-1β in chondrocytes, which contributes to the development and progression of OA. Ku et al.  reported that SF leptin concentrations in patients with knee OA was positively associated with Kellgren/Lawrence (K/L) grading and could be used as an effective marker for quantitative detection of OA severity. All these results intensively suggested that leptin mediated a detrimental effect on cartilage metabolism and may thus play potential roles in the pathogenesis of OA.
LEPR is a membrane-spanning glycoprotein belonging to the class I cytokine-receptor family which mediates the effects of over 20 kinds of cytokines. Among the widely expressed alternative spliced isforms (including Ob-Ra, Ob-Rb, Ob-Rc, Ob-Rd, and Ob-Re), the long isoform Ob-Rb is the only 1 indentified to be functional . Recent study showed that LEPR was also expressed within the native human cartilage chondrocytes, mediating the inflammatory and destructive responses of leptin in cartilage and other joint tissues . Simopoulou et al.  detected the expression of LEPR in articular cartilage of controls and patients with minimal OA and demonstrated the expressions of LEPR mRNA and protein were both significantly increased in advanced osteoarthritic cartilage compared to that of patients with minimal OA and normal controls. Since LEPR is of functional importance for the risk factor of OA-obesity, numerous investigations regarding the association between LEPR SNPs with obesity have been extensively carried out [13–16]. Among the various SNPs of LEPR, Gln223Arg (rs1137101) variant of LEPR exon 4 received more attention. By conducting a population-based case–control study, Portoles et al.  discovered the Gln223Arg polymorphism was significantly associated with obesity in Spain population. Other studies demonstrated that the Gln223Arg polymorphism is highly related to obesity in Pacific Islanders , body mass index (BMI), body fat amount, serum leptin concentrations in Caucasian woman  and antipsychotic-induced body weight gain .
Although there have been studies demonstrating the association of leptin gene SNPs with knee OA susceptibility  and it seems that these LEPR polymorphisms are associated with risk factors for knee OA, yet few studies have been conducted directly regarding the association of LEPR SNPs and susceptibility to knee OA. In 2003, Ma et al. . genotyped LEPR SNP A668G (Gln223Arg) in 148 patients with knee OA and 155 controls, concluding that Gln223Arg is associated with susceptibility to primary knee OA in Ningxia Hui population, but this association in Chinese Han population still remains unknown. Therefore, the present case − control study aims to investigate whether the genetic polymorphisms of Gln223Arg in LEPR gene is associated with OA susceptibility in Northwest Chinese population with Han ethnicity.
The procedure of the study was reviewed by the Ethical Review Committee of Health Science Center, First Affiliated Hospital of Xi’an Jiaotong University and conducted in line with the Declaration of Helsinki. Written informed consent was obtained from all participants prior to their participation in the study. A total of 1 215 subjects were recruited in this study, consisting of 587 patients (134 males and 453 females) with confirmed knee OA and 628 age − and sex − matched healthy controls (143 males and 485 females). The diagnosis of OA was based on the symptomatic criteria of American College of Rheumatology and radiographic signs (K/L grading ≥ 2) for OA  of at least 1 knee. The patients were recruited at the Department of Orthopaedics, First Affiliated Hospital of Health Science Center, Xi’an Jiaotong University between January 2010 and December 2014. Controls with normal knee radiographs (K/L grading < 2) undergoing physical examination were consecutively collected in a large member of individuals contemporaneously. Participants were excluded on the basis of having arthropathy due to gout, pseudogout, rheumatoid arthritis, with histories of corticosteroid medication, bilateral knee replacements, other forms of arthritis, cancer, or other chronic inflammation diseases. All individuals enrolled in the present study were drawn from a local Chinese Han ethnic population of Shaanxi Province in Northwest China. All subjects were stable residents in the areas and no 1 dropped out during the process of study.
Venous blood samples were collected from all subjects with informed consent after overnight fasting. Genomic DNA was all extracted from blood samples with the commercially available DNA extraction kit (TIANGEN Biotech, Bejing, China) according to the manufacturer’s instructions and stored at − 80 °C until utilized. The genotypes of Gln223Arg polymorphism in LEPR gene were determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods. Sequencing primers for PCR were designed with Primer Express 3.0 software (Applied Biosystems, Waltham, MA). The primers used for amplification were: Forward 5′-TCCTGCTTTAAAAGCCTAATCCAGTATTT-3′ and Reverse 5′-AGCTAGCAAATATTTTTGTAAGCAAT-3′ for Gln223Arg. The PCR reaction was performed on 25 μL reaction mixture containing 4 μL mixed DNA template, 1 μL of each primer, 12.5 μL 2 × Reaction Mix and 6.25 μL ddH2O. The PCR protocol was carried out as followed: initial denaturation at 95 °C for 4 min followed by 30 cycles of 94 °C for 30 s, annealing at 59.5 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 7 min. PCR products were digested over night at 37 °C with 5U of MspI restriction enzymes (Thermo Scientific, Waltham, MA, USA) for Gln223Arg polymorphisms. The digested mixtures were then separated on a 2.0 % agarose gel and visualized by ethidium bromide staining and ultraviolet light transillumination. The expected fragment of Gln223Arg was 242 bp and 125 bp in AA, 367 bp in GG and 367, 242 and 125 bp in AG genotypes, respectively. Fifty percent of the samples of both cases and control subjects were randomly selected for reproducibility tests at least twice in different days by using direct DNA sequencing with ABI3500DX Genetic Analyzer (Thermo Scientific, Waltham, MA, USA) to verify the accuracy of PCR-RFLP method.
Statistical analyses were performed using SPSS software, version 16.0 for windows (SPSS Inc.; Chicago, IL, USA). Baseline quantitative data are represented as mean ± SD. The Chi-square goodness-of-fit test was employed to assess the Hardy-Weinberg equilibrium (HWE) for genotype and allelic proportions in the cases and control groups. The associations between the Gln223Arg polymorphism and the risk of knee OA was evaluated with odds ratios (ORs) and 95 % confidence intervals (CI) after adjusting for age and BMI with multivariate logistic regression analysis. Further stratifications analysis for sex, age and the radiographic severity of OA were subsequently conducted. Kolmogorov-Smirnov test was performed to analyze the data normality and unpaired t test, Mann–Whitney U test, or chi-square test was used to assess significance in clinical characteristics between patients with knee OA and healthy controls, as appropriate. Clinical characteristics of OA patients with different genotype were compared using chi-square or Kruskal-Wallis tests. For all the tests, P < 0.05 was considered as statistically significant. The power of the study was estimated with the published data in the literature and then calculated with Power and Precision software (Biostat, Englewood, NJ, USA), with results reaching over 75 %.
Baseline clinical parameters of the study groups
A total of 1 215 Northwest Chinese subjects with Han ethnicity were recruited in this case–control study, comprising 587 patients with knee OA (134 males and 453 females; mean age 61.37 ± 10.64 years) and 628 healthy controls (143 males and 485 females; mean age 60.46 ± 9.85 years). The baseline characteristics of both groups were summarized in Table 1. No significant differences in age and gender were found between the 2 groups (P > 0.05). Similar with the previous studies that reported high BMI increased the risk of developing OA , BMI was found significantly higher in patients group compared to control group (P < 0.001). According to the radiographic KL grading system, approximately 50.3 % of the patients had a K/L grading of 3 or 4. There were no significant differences in the age, gender and BMI between different K/L grading subgroups (P > 0.05).
Genotypic distribution and allelic frequencies of Gln223Arg polymorphism
The results of the genotype distribution and allele frequencies of LEPR Gln223Arg for cases and control groups were presented in Table 2. LEPR Gln223Arg variant genotype frequencies in either the cases or the control group was found to be in Hardy–Weinberg equilibrium (P = 0.85 and 0.83, respectively). The genotype GG and allele-G were the most common genotype and allele in both cases and control groups. The allele frequencies of Gln223Arg in cases group (G, 68.57 %; A, 31.43 %) were remarkably different from those of control group (G, 64.25 %; A, 35.75 %, chi-square = 5.06, P = 0.024). Besides, the genotype distributions of Gln223Arg in cases group (GG, 45.49 %; GA, 46.17 %; AA, 8.35 %) were also significantly different from those in control group (GG, 41.08 %; GA, 46.34 %; AA, 12.58 %; chi-square = 6.52, P = 0.038).
Association of Gln223Arg polymorphism and susceptibility of knee OA
Table 3 illustrates the genetic association between the Gln223Arg polymorphism of LEPR gene and risk of knee OA. The data demonstrated that the Gln223Arg genetic polymorphism of LEPR gene was significantly related to the susceptibility of knee OA in patients as a whole and controls in the dominant model (GG + GA versus AA) (OR = 1.58; 95 % CI, 1.09–2.30; P = 0.016) and G allele was correlated with the higher risk of knee OA (OR = 1.21; 95 % CI, 1.03–1.44; P = 0.024). When further stratified by sex after adjusting for other variables such as age and BMI, multivariate logistic regression analysis showed the evident differences of genotype distribution in female subjects (OR = 1.89; 95 % CI, 1.25–2.87; P = 0.008). Additionally, G allele carriers were still related to the increased risk of knee OA in females (OR = 1.56; 95 % CI, 1.08–2.49; P = 0.012). However, both of the genotype distributions and allele frequencies in males were found no significant differences. When stratification by age, both young (≤65 years) (OR = 1.48; 95 % CI, 1.18–2.36; P = 0.028) and old (>65 years) (OR = 1.72; 95 % CI, 1.13–2.53; P = 0.020) patients showed significant differences in combined genotype (GG + GA) versus AA after adjustment for BMI and gender compared to control group. Besides, G allele carriers suffer higher risk of knee OA in both young (OR = 1.26; 95 % CI, 1.04–1.58; P = 0.045) and old patients (OR = 1.38; 95 % CI, 1.13–1.57; P = 0.037). After stratification according to the radiographic severity of knee OA, significant differences were still found to be existed in the comparison of combined genotype (GG + GA) versus AA after adjusting for variables such as age, BMI, and gender in cases with mild knee OA (OR = 1.53; 95 % CI, 1.05–2.47; P = 0.022). Additionally, G allele was moderately associated with an increased risk of knee OA in mild subgroup (OR = 1.30; 95 % CI, 1.05–1.62; P = 0.013). Unexpectedly, we didn’t detect any association between genotype distributions or allele frequencies of rs1137101 and the risk of knee OA in severe subgroup.
Identification of patients at risk for incident disease or disease progression in OA remains challenging, because the right now preferred radiography method is an insensitive reflection of genetic changes that presage cartilage and bone abnormalities. Although the etiology of OA is multifactorial, there have been compelling evidence suggesting the involvement of genetic factors in the pathogenesis of OA . Thus, there is a widely appreciated need for identification of susceptibility genes of OA. In recent years, many studies have been intensively carried out to determine the related genes that may cause knee OA [6, 22–24]. Among which, the impaired leptin signal pathway was demonstrated to be highly related to the pathophysiology of OA. Considering the crucial role of LEPR in regulating leptin activity, here we investigated the effects of the common exonic polymorphism of Gln223Arg in LEPR gene on the risk of knee OA in Northwest Chinese population with Han ethnicity.
The present case–control study firstly examined whether the LEPR SNP polymorphism was correlated with the presence of primary knee OA in Northwest Chinese population. Our results indicated that LEPR gene SNP rs1137101 (Gln223Arg) polymorphism was correlated with risk of knee OA. We have found GG genotype distribution and G allelic frequencies of Gln223Arg were both significantly higher in OA patients than that of control group, and G allele of rs1137101 is a common susceptibility allele for knee OA at least in the Han population in Northwest China. Multivariate logistic regression indicated that rs1137101 was associated with knee OA in the dominant genetic model (GG + GA versus AA) and G allele carriers of Gln223Arg had higher risks of knee OA, which is consistent with reports from Spanish Mediterranean population and Pacific Islanders [13, 14]. These 2 studies both detected the significant association of Q223R (Gln223Arg) polymorphism and prevalence of obesity in the dominant model for Q223R (QQ + QR versus RR) and revealed a significant association between the 223Q carrier genotype and obesity [13, 14]. However, the population we examined here had a BMI in a normal range (not obese) and significant associations of Gln223Arg and risk of knee OA were still detected, which indicates this association is not the only risk factor for knee OA in the Chinese Han population and implies the possibility that the effects of Gln223Arg appear in correlating with other genetic or environmental factors. Therefore, further studies are needed to investigate the potential role of Gln223Arg for predicting risk of knee OA in Chinese Han population in the absence of obesity.
On stratification analysis, the G allele carriers (GG + GA) and G allele were related with knee OA susceptibility in females, both young (≤65 years) and old groups (>65 years) and patients with mild knee OA after adjustment for other potential confounders such as age, gender and BMI. But the association between the risk of knee OA and rs1137101 G allele carriers (GG + GA) was stronger in female patients (1.72-fold) and those >65 years old (1.89-fold) compared with AA homozygote. However, it should be noted there were no associations between Gln223Arg polymorphism and knee OA susceptibility in male participants and patients with higher K/L grading. We assume that the disappearance of association in male participants may be ascribed to the relative small sample size and/or sex disparities in the genetics of OA. Although sex disparities in the genetics of OA have not been fully understood, in addition to the effect of estrogen, a higher incidence and heritability of knee OA has been reported for women than for men [5, 25–27], as well as sex differences in cartilage volume . The lack of difference in patients with higher K/L grading can be attributed to the possibility that rs1137101 might just play important roles in knee OA induction (with relative lower K/L grading), not in the radiographic progression of OA and that rs1137101 may be in linkage disequilibrium with some other more relevant alleles within the LEPR gene region. Previous studies of candidate genes for OA demonstrated different SNPs may function differently in the induction and progression of knee OA. ADAM12 gene was found to be involved in both the induction and radiographic progression of knee OA . NCOR2 and CD36 genes were associated with induction of radiographic knee OA measured as either a K/L grading of ≥2.0 or the presence of osteophytes , while the polymorphisms at CILP, TNA, IL1RN and variation at inflmmatory genes appeared to correlate with radiographic progression over time [29–31]. Besides, mitochondrial DNA haplogroups are also associated with the induction and progression of knee OA . Thus, polygenic effects on the pathogenesis of knee OA should not be ignored and there exists a necessity to explore other potentially contributing genes involved in the development of knee OA for a better comprehension of genetic regulation on knee OA.
Originally, the SNP rs1137101 designates a change of glutamine to an arginine at codon 223 (Gln223Arg) in the extracellular domain coding region of the LEPR protein . This transition results in an exchange of a neutral amino acid for 1 with a positive charge, altering receptor structure and function and affecting signaling capacity and circulating leptin levels . Previous studies have mostly reported the association of Gln223Arg polymorphism and risk of different kinds of cancers [34–36] and cardiovascular disease [37–39]. Gln223Arg was recently recognized as a risk predictor of knee OA in Ningxia Hui population, however, there have been no detailed data concerning the association LEPR SNP rs1137101 with risk of OA in Chinese Han population thus far. We were the first to show that Gln223Arg polymorphism was correlated with risk of knee OA in Northwest Chinese population with Han ethnicity. The correlation between the Gln223Arg polymorphism and knee OA susceptibility could be possible interpreted by physiological role of leptin. As described previously, leptin can also be produced by many other tissues and cells other than adipose tissue, such as joint tissues [8, 10]. Besides, LEPR was also expressed within the native human cartilage, mediating the inflammatory and destructive responses of leptin in cartilage and other joint tissues . Further mechanical studies indicated that the Gln223Arg polymorphism might affect leptin signaling through phosphorylation of signal transducers and activators of transcription 3 (STAT3)  and was in linkage disequilibrium with other functional polymorphisms that can impair the signaling capacity of the LEPR . The G to A mutation in rs1137101 led to reduced proliferation and synthesis of extracellular matrix and an amelioration of the anti-apoptotic effect of leptin on articular chondrocytes and eventually causes OA happen.
However, we have been aware of several potential limitations existing in this study that could affect our conclusions. First, as this was a study comprising a relatively small number of subjects, it is difficult to draw a strong conclusion. Further studies with a larger population sample are warranted to substantiate our results. Second, the P-value (P = 0.537) for age-matching between patients and healthy controls appeared to be marginal, which might affect our ability to detect the association between Gln223Arg polymorphism and the susceptibility of OA to some extent. Third, we have just focused on the association and missed the aspects of the functional study of LEPR gene polymorphism, which might potentially neglect other gene SNPs that may strongly correlate with OA. Lastly, we did not analyze the relationship between Gln223Arg polymorphism and joint pain experienced by patients. Further studies are warranted to detect the role of Gln223Arg polymorphism for predicting joint pain in patients with knee OA.
In conclusion, this study firstly demonstrated that the LEPR rs1137101 (Gln223Arg) polymorphism was associated with the susceptibility of knee OA in the dominant genetic model (GG + GA versus AA) in Northwest Chinese population with Han ethnicity. This association remained in patients with mild radiographic severity and female subjects. Further studies with larger-scale population and different ethnicities should be conducted to substantiate the current findings.
95%CI, 95 % conidence intervals; BMI, body mass index; K/L, Kellgren/Lawrence; LEPR, leptin receptor; OA, osteoarthritis; Ors, odds ratios; PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphism; SNP, single nucleotide polymorphism
Reginster JY. The prevalence and burden of arthritis. Rheumatology (Oxford). 2002;41(Supp 1):3–6.
Sun BH, Wu CW, Kalunian KC. New developments in osteoarthritis. Rheum Dis Clin North Am. 2007;33(1):135–48.
Xiang YJ, Dai SM. Prevalence of rheumatic diseases and disability in China. Rheumatol Int. 2009;29(5):481–90.
Spector TD, MacGregor AJ. Risk factors for osteoarthritis: genetics. Osteoarthritis Cartilage. 2004;12:S39–44.
Teichtahl AJ, Wluka AE, Proietto J, Cicuttini FM. Obesity and the female sex, risk factors for knee osteoarthritis that may be attributable to systemic or local leptin biosynthesis and its cellular effects. Med Hypotheses. 2005;65(2):312–5.
Valdes AM, Doherty S, Muir KR, Zhang W, Maciewicz RA, Wheeler M, Arden N, Cooper C, Doherty M. Genetic contribution to radiographic severity in osteoarthritis of the knee. Ann Rheum Dis. 2012;71(9):1537–40.
Woods AJ, Stock MJ. Leptin activation in hypothalamus. Nature. 1996;381(6585):745.
Dumond H, Presle N, Terlain B, Mainard D, Loeuille D, Netter P, Pottie P. Evidence for a key role of leptin in osteoarthritis. Arthritis Rheum. 2003;48(11):3118–29.
Simopoulou T, Malizos KN, Iliopoulos D, Stefanou N, Papatheodorou L, Ioannou M, Tsezou A. Differential expression of leptin and leptin’s receptor isoform (Ob-Rb) mRNA between advanced and minimally affected osteoarthritic cartilage; effect on cartilage metabolism. Osteoarthritis Cartilage. 2007;15(8):872–83.
Ku JH, Lee CK, Joo BS, An BM, Choi SH, Wang TH, Cho HL. Correlation of synovial fluid leptin concentrations with the severity of osteoarthritis. Clin Rheumatol. 2009;28(12):1431–5.
Baumann H, Morella KK, White DW, Dembski M, Bailon PS, Kim H, Lai CF, Tartaglia LA. The full-length leptin receptor has signaling capabilities of interleukin 6-type cytokine receptors. Proc Natl Acad Sci U S A. 1996;93(16):8374–8.
Figenschau Y, Knutsen G, Shahazeydi S, Johansen O, Sveinbjornsson B. Human articular chondrocytes express functional leptin receptors. Biochem Biophys Res Commun. 2001;287(1):190–7.
Portoles O, Sorli JV, Frances F, Coltell O, Gonzalez JI, Saiz C, Corella D. Effect of genetic variation in the leptin gene promoter and the leptin receptor gene on obesity risk in a population-based case–control study in Spain. Eur J Epidemiol. 2006;21(8):605–12.
Furusawa T, Naka I, Yamauchi T, Natsuhara K, Kimura R, Nakazawa M, Ishida T, Inaoka T, Matsumura Y, Ataka Y, et al. The Q223R polymorphism in LEPR is associated with obesity in Pacific Islanders. Hum Genet. 2010;127(3):287–94.
Quinton ND, Lee AJ, Ross RJ, Eastell R, Blakemore AI. A single nucleotide polymorphism (SNP) in the leptin receptor is associated with BMI, fat mass and leptin levels in postmenopausal Caucasian women. Hum Genet. 2001;108(3):233–6.
Brandl EJ, Frydrychowicz C, Tiwari AK, Lett TA, Kitzrow W, Buttner S, Ehrlich S, Meltzer HY, Lieberman JA, Kennedy JL, et al. Association study of polymorphisms in leptin and leptin receptor genes with antipsychotic-induced body weight gain. Prog Neuropsychopharmacol Biol Psychiatry. 2012;38(2):134–41.
Qin J, Shi D, Dai J, Zhu L, Tsezou A, Jiang Q. Association of the leptin gene with knee osteoarthritis susceptibility in a Han Chinese population: a case–control study. J Hum Genet. 2010;55(10):704–6.
Ma XJ, Guo HH, Hao SW, Sun SX, Yang XC, Yu B, Jin QH. Association of single nucleotide polymorphisms (SNPs) in leptin receptor gene with knee osteoarthritis in the Ningxia Hui population. Yi Chuan. 2013;35(3):359–64.
Kellgren JH, Lawrence JS. Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957;16(4):494–502.
Manek NJ, Hart D, Spector TD, MacGregor AJ. The association of body mass index and osteoarthritis of the knee joint: an examination of genetic and environmental influences. Arthritis Rheum. 2003;48(4):1024–9.
Yucesoy B, Charles LE, Baker B, Burchfiel CM. Occupational and genetic risk factors for osteoarthritis: a review. Work. 2015;50(2):261–73.
Wang L, Guo L, Tian F, Hao R, Yang T. Analysis of single nucleotide polymorphisms within ADAM12 and risk of knee osteoarthritis in a Chinese Han population. Biomed Res Int. 2015;2015:518643.
Zhang R, Yao J, Xu P, Ji B, Voegeli G, Hou W, Li H, Wang Y, Kelsoe JR, Ma J. Association between genetic variants of DVWA and osteoarthritis of the knee and hip: a comprehensive meta-analysis. Int J Clin Exp Med. 2015;8(6):9430–7.
Tawonsawatruk T, Changthong T, Pingsuthiwong S, Trachoo O, Sura T, Wajanavisit W. A genetic association study between growth differentiation factor 5 (GDF 5) polymorphism and knee osteoarthritis in Thai population. J Orthop Surg Res. 2011;6:47.
Srikanth VK, Fryer JL, Zhai G, Winzenberg TM, Hosmer D, Jones G. A meta-analysis of sex differences prevalence, incidence and severity of osteoarthritis. Osteoarthritis Cartilage. 2005;13(9):769–81.
Kaprio J, Kujala UM, Peltonen L, Koskenvuo M. Genetic liability to osteoarthritis may be greater in women than men. BMJ. 1996;313(7051):232.
Valdes AM, Loughlin J, Oene MV, Chapman K, Surdulescu GL, Doherty M, Spector TD. Sex and ethnic differences in the association of ASPN, CALM1, COL2A1, COMP, and FRZB with genetic susceptibility to osteoarthritis of the knee. Arthritis Rheum. 2007;56(1):137–46.
Faber SC, Eckstein F, Lukasz S, Muhlbauer R, Hohe J, Englmeier KH, Reiser M. Gender differences in knee joint cartilage thickness, volume and articular surface areas: assessment with quantitative three-dimensional MR imaging. Skeletal Radiol. 2001;30(3):144–50.
Valdes AM, Hart DJ, Jones KA, Surdulescu G, Swarbrick P, Doyle DV, Schafer AJ, Spector TD. Association study of candidate genes for the prevalence and progression of knee osteoarthritis. Arthritis Rheum. 2004;50(8):2497–507.
Wu X, Kondragunta V, Kornman KS, Wang HY, Duff GW, Renner JB, Jordan JM. IL-1 receptor antagonist gene as a predictive biomarker of progression of knee osteoarthritis in a population cohort. Osteoarthritis Cartilage. 2013;21(7):930–8.
Valdes AM, Hassett G, Hart DJ, Spector TD. Radiographic progression of lumbar spine disc degeneration is influenced by variation at inflammatory genes: a candidate SNP association study in the Chingford cohort. Spine (Phila Pa 1976). 2005;30(21):2445–51.
Rego-Perez I, Fernandez-Moreno M, Fernandez-Lopez C, Arenas J, Blanco FJ. Mitochondrial DNA haplogroups: role in the prevalence and severity of knee osteoarthritis. Arthritis Rheum. 2008;58(8):2387–96.
Li Y, Geng J, Wang Y, Lu Q, Du Y, Wang W, Li Z. The role of leptin receptor gene polymorphisms in determining the susceptibility and prognosis of NSCLC in Chinese patients. J Cancer Res Clin Oncol. 2012;138(2):311–6.
Mohammadzadeh G, Ghaffari MA, Bafandeh A, Hosseini SM. Effect of leptin receptor Q223R polymorphism on breast cancer risk. Iran J Basic Med Sci. 2014;17(8):588–94.
Rodrigues PR, Maia LL, Santos M, Peterle GT, Alves LU, Takamori JT, Souza RP, Barbosa WM, Mercante AM, Nunes FD, et al. Leptin receptor expression and Gln223Arg polymorphism as prognostic markers in oral and oropharyngeal cancer. Genet Mol Res. 2015;14(4):14979–88.
Karimi K, Arkani M, Safaei A, Pourhoseingholi MA, Mohebbi SR, Fatemi SR, Vafaei M. Association of leptin receptor gene Gln223Arg polymorphism with susceptibility to colorectal cancer. Gastroenterol Hepatol Bed Bench. 2011;4(4):192–8.
Gu P, Jiang W, Chen M, Lu B, Shao J, Du H, Jiang S. Association of leptin receptor gene polymorphisms and essential hypertension in a Chinese population. J Endocrinol Invest. 2012;35(9):859–65.
Liu Y, Lou YQ, Liu K, Liu JL, Wang ZG, Wen J, Zhao Q, Wen SJ, Xiao L. Role of leptin receptor gene polymorphisms in susceptibility to the development of essential hypertension: a case–control association study in a Northern Han Chinese population. J Hum Hypertens. 2014;28(9):551–6.
Zheng H, Xie N, Xu H, Huang J, Xie X, Luo M. Association of Gln223Arg polymorphism of the leptin receptor with hypertensive left ventricular hypertrophy. Folia Biol (Praha). 2013;59(6):246–52.
Mackey-Lawrence NM, Guo X, Sturdevant DE, Virtaneva K, Hernandez MM, Houpt E, Sher A, Porcella SF, Petri Jr WA. Effect of the leptin receptor Q223R polymorphism on the host transcriptome following infection with Entamoeba histolytica. Infect Immun. 2013;81(5):1460–70.
Riestra P, Garcia-Anguita A, Torres-Cantero A, Bayonas MJ, De Oya M, Garces C. Association of the Q223R polymorphism with age at menarche in the leptin receptor gene in humans. Biol Reprod. 2011;84(4):752–5.
We are grateful to the patients for donating their samples to our research studies, staff within the research and clinical teams at Rehabilitation Center and Department of Orthopedics, First Affiliated Hospital of Xi’an Jiaotong University Health Science Center for help in obtaining and processing samples for funding this research.
This work was funded by Nature Science Foundation of Shaanxi Province (2014TM4193).
Availability of data and materials
The data cannot be publicly shared because permission was not obtained from the participants but are available from the corresponding author on reasonable request.
YJH and DH performed the genetic studies. DH and LJG recruited patients and were involved in the clinical part of the investigation. ZLH performed the statistical analysis. YJH conceived the study, participated in its design and coordination and helped to draft the manuscript. All of the authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethical approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the Ethical Review Committee of the First Affiliated Hospital of Xi’an Jiaotong University Health Science Center and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Written informed consent was obtained from all subjects prior to their participation in the present study.