- Research article
- Open Access
Association of single-nucleotide polymorphisms in RHOB and TXNDC3 with knee osteoarthritis susceptibility: two case-control studies in East Asian populations and a meta-analysis
Arthritis Research & Therapy volume 10, Article number: R54 (2008)
Conflicting findings on the association of single nucleotide polymorphisms (SNPs) in RHOB and TXNDC3 with susceptibility to knee osteoarthritis (OA) have been reported in European Caucasians. To examine the associations of these SNPs with OA in East Asian populations and to evaluate their global significance, we conducted two case-control studies in 955 Chinese and 750 Japanese patients.
We genotyped the previously implicated SNPs rs585017 (in RHOB) and rs4720262 (in TXNDC3) in patients with primary symptomatic knee OA with radiographic confirmation and in matched control individuals, and analyzed their associations. We further conducted a meta-analysis of the study findings together with those of previously reported European studies using the DerSimonian-Laird procedure.
A significant association of RHOB with knee OA was observed in male Chinese patients (P = 0.02). No significant associations were found for RHOB in any other comparisons in the East Asian populations. The association of TXNDC3 was replicated in Chinese female (P = 0.04) and Japanese (P = 0.03) patients, although none of these associations persisted after Bonferroni correction. Significant association (P = 0.02 for the allelic frequency) with nonsignificant heterogeneity was found in the East Asian replication study. No significant association was found in any comparison in the meta-analysis for all studies.
Our study replicates the association, previously reported in European Caucasians, of TXNDC3 with knee OA susceptibility in an East Asian population.
Osteoarthritis (OA; OMIM [Online Mendelian Inheritance in Man] 165720) is a type of arthritis that is caused by breakdown and eventual loss of the cartilage of synovial joints. OA is the most common type of arthritis, with a high incidence in East Asian populations [1, 2]. Epidemiological studies have shown that OA has a strong genetic component, and several susceptibility genes for OA have been identified [3–6].
Marh and coworkers  examined regulatory polymorphisms in the 5' regions of the genes that potentially allow for differential expression in vivo, and they found positive association for RHOB and TXNDC3 in European Caucasians living in Germany. RHOB belongs to the family of small GTPases, and is constitutively expressed and essential in adult articular chondrocytes, but it is significantly downregulated in OA chondrocytes . RHOB (ras homolog gene family, member B), the protein encoded by RHOB, is important in the induction of apoptotic cell death that occurs in response to DNA damage , and chondrocytes are known to undergo significant DNA damage in OA .TXNDC3 encodes a thioredoxin protein, and its transcripts appear to be alternatively spliced in monocytes and chondrocytes, with the chondrocyte-specific transcripts lacking exon 2 . Thus, both are good candidate OA genes. However, a replication study in the UK found no association of RHOB and TXNDC3 with knee OA, even though the ethnicity and disease ascertainment were basically the same in the two studies .
To examine the replication of association between the RHOB and TXNDC3 SNPs and knee OA susceptibility, we conducted two case-control studies in Han Chinese and Japanese populations and a meta-analysis.
Materials and methods
A population of 955 Chinese individuals were studied; 469 (323 women and 146 men) were consecutively enrolled patients at the Center of Diagnosis and Treatment for Joint Disease (Drum Tower Hospital, affiliated to Medical School of Nanjing University), and 486 were healthy control individuals (316 women and 170 men), who were enrolled at the Center of Physical Examination. All individuals included in the study were Han Chinese living in and around Nanjing. None dropped out during the study. In addition, a population of 750 Japanese individuals were included; 376 patients (327 women and 49 men) and 374 healthy control individuals (116 women and 258 men) who were living in and around Tokyo were recruited by participating hospitals.
Included were patients with knee OA who not only had definite signs and symptoms of OA but also had radiographic evidence of OA. All patients had pain with rest and/or night pain over 5-month period. Patients with knee diseases such as inflammatory, post-traumatic and post-septic arthritis, and skeletal dysplasia were excluded. Radiographic OA was assessed using the Kellgren-Lawrence grading system . Only patients with Kellgren-Lawrence grades of 2 or higher were included. More than 70% of patients had a Kellgren-Lawrence score of 3 or 4. None of the control individuals had ever exhibited any signs or symptoms of arthritis or joint diseases.
The ages of the patients and the controls (mean ± standard deviation) in the Chinese study were 58 ± 13 (range 32 to 93) years and 57 ± 12 (range 40 to 97) years, respectively. We calculated the body mass index to assess the effect of obesity. The body mass index (mean ± standard deviation) was 25 ± 4 kg/m2 for patients and 24 ± 7 kg/m2 for control individuals. There was no statistically significant difference between the two groups. The ages of the patients and the control individuals (mean ± standard deviation) in the Japanese study were 73 ± 7 years and 63 ± 10 years, respectively.
The study protocol was approved by the ethnical committees of the Medical School of Nanjing University, and single nucleotide polymorphism (SNP) Research Center of RIKEN, and informed consent was obtained from patients and control individuals.
The RHOB SNP rs585017 and TXNDC3 SNP rs4720262  were genotyped by the 5-nuclease assay (Taqman) using the Mx3000P Real Time PCR System (Stratagene, La Jolla, CA, USA) or ABI 7700 (Applied Biosystems, Foster City, CA, USA). The primers, probes and reaction conditions are available upon request. Genotyping was conducted by laboratory personnel who were blinded to subject status. A random 5% of the samples were subjected to repeat analysis to validate the genotyping procedures. Two authors independently reviewed the genotyping results, data entry and statistical analyses. The distributions of genotypes in all knee OA and control groups conformed to Hardy-Weinberg equilibrium.
Standard χ2 analysis-of-contingency tables were used to compare genotype and allele distributions of RHOB and TXNDC3 in the case-control study. Differences in clinical characteristics, age, sex and body mass index between genotypes were tested by Mann-Whitney test, Kruskal-Wallis test, or χ2 test using SPSS 12.0 system software (SPSS Inc., Chicago, IL, USA). We assessed association of rs585017 and rs4720262 with stratifications and Hardy-Weinberg equilibrium using the χ2 test. Software R was used in the meta-analysis. Genotype data from previous reports [7, 10] and our own genotype data were analyzed using the DerSimonian-Laird procedure , based on the random effects model. Power estimates were calculated using the software R.
The allelic frequency of the susceptible G allele of rs585017 was similar in Chinese and Japanese individuals, and was low (< 4%), unlike that in European Caucasians (> 20%) [8, 10]. No GG genotype was detected in East Asian individuals (Table 1). Positive associations with OA were found in comparisons of genotypes and alleles in Chinese men (both P = 0.02; Table 2), but these associations did not persist after Bonferroni correction, and no such association was found in Japanese individuals (Table 2). No significant association was detected in any other comparison in East Asian populations (P = 0.28 and P = 0.27 in recessive and allele modes, respectively). We had 99% and 63% power to detect odds ratios of 2.1 and 1.5 in our study.
Frequencies of the susceptible T allele of rs4720262 were different between Chinese (9.7% in control individuals) and Japanese (14.2% in control individuals), which were significantly different from the frequency in the UK control individuals  (Table 3). A positive association was found only for Japanese patients versus control individuals in the comparisons of TT versus other genotypes combined (P = 0.03) and allelic frequency (P = 0.04; Table 4). No significant association was detected after stratification by sex. Because the results were mixed, we performed the meta-analysis for the four studies (first German study  and three replication studies) and for replication studies. No significant associations were observed in any comparison. Because marginally significant heterogeneity was detected in these analyses (Table 5), we further examined associations stratified by ethnicity. There was a significant association for TXNDC3 in East Asian populations with nonsignificant heterogeneity (P = 0.03 for CC versus CT+TT; P = 0.02 for C allele versus T allele; Table 5).
In RHOB, the minor allele frequency in East Asian individuals was below 0.05 and much lower than that in European Caucasians. We could not detect a definite association with OA in East Asian individuals. A weak association was found in male Chinese when the patients were stratified by sex, but the association did not persist after correction with multiple testing, and the trend of difference was reversed in female patients. Meta-analysis also yielded negative results. More than 80% power at a significance level of 5% is desirable for a 'negative' meta-analysis . Our study had adequate power if we had identified the same odds ratio as the original study , but we did not have sufficient power to detect an odds ratio as small as 1.5. This lack of association requires confirmation in further replication studies.
The minor allele frequencies of TXNDC3 in East Asian individuals were significantly different from those in UK control individuals. The corresponding German frequency was 13%, which was within the Asian frequency range. It appears that rs4720262 provides a good example for the variability of SNPs between different ethnicities or geographic locations. Weak associations were observed in the Chinese and Japanese populations in some comparisons. These associations did not persist after rigorous correction for multiple testing, and the meta-analysis found no associations in the analyses for all and replication studies; however, the analysis after stratification by ethnicity detected a significant association in the East Asian population. Had the meta-analysis of the East Asian study data found an absence of correlation, then it wold be possible to conclude with confidence that there is no association.
The association of TXNDC3 with knee OA was replicated in East Asian individuals, whereas a lack of association for RHOB was identified. Ethnic differences have been noted in OA susceptibility genes, especially in minor allele frequency between two genes. The association of LRCH1 in European  has not been replicated in Asian  patients, and the association for CALM1 identified in Japanese hip OA  has not been detected in UK Caucasian patients . In ASPN, frequencies of the susceptible D14 allele differ among European patients and between European and Asian patients . The association of ASPN with knee OA is found globally, but the effect of ASPN exhibits considerable ethnic differences . There is a clear, global link between GDF5 and knee OA, but the effect size differs considerably between Europeans and Asian patients [5, 20]. As noted in the LRCH1 study , the difference in the ascertainment criteria is unlikely to account for the discrepancy. A worldwide association study along with functional studies of the susceptibility SNP should be performed to clarify the significance of TXNDC3 as an OA susceptibility gene.
Our study replicates the association, previously reported in European Caucasians, of TXNDC3 with knee OA susceptibility in an East Asian population.
= single nucleotide polymorphism.
Zhang Y, Terkeltaub R, Nevitt M, Xu L, Neogi T, Aliabadi P, Niu J, Felson DT: Lower prevalence of chondrocalcinosis in Chinese subjects in Beijing than in white subjects in the United States: the Beijing Osteoarthritis Study. Arthritis Rheum. 2006, 54: 3508-3512. 10.1002/art.22189.
Kim HA, Kim S, Seo YI, Choi HJ, Seong SC, Song YW, Hunter D, Zhang Y: The epidemiology of total knee replacement in South Korea: national registry data. Rheumatology (Oxford). 2008, 47: 88-91. 10.1093/rheumatology/kem308.
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: 137-146. 10.1002/art.22301.
Ikegawa S: New gene associations in osteoarthritis: what do they provide, and where are we going?. Curr Opin Rheumatol. 2007, 19: 429-434. 10.1097/BOR.0b013e32825b079d.
Miyamoto Y, Mabuchi A, Shi D, Kubo T, Takatori Y, Saito S, Fujioka M, Sudo A, Uchida A, Yamamoto S, Ozaki K, Takigawa M, Tanaka T, Nakamura Y, Jiang Q, Ikegawa S: A functional polymorphism in the 5' UTR of GDF5 is associated with susceptibility to osteoarthritis. Nat Genet. 2007, 39: 529-533. 10.1038/2005.
Kizawa H, Kou I, Iida A, Sudo A, Miyamoto Y, Fukuda A, Mabuchi A, Kotani A, Kawakami A, Yamamoto S, Uchida A, Nakamura K, Notoya K, Nakamura Y, Ikegawa S: An aspartic acid repeat polymorphism in asporin inhibits chondrogenesis and increases susceptibility to osteoarthritis. Nat Genet. 2005, 37: 138-144. 10.1038/ng1496.
Mahr S, Burmester GR, Hilke D, Göbel U, Grützkau A, Häupl T, Hauschild M, Koczan D, Krenn V, Neidel J, Perka C, Radbruch A, Thiesen HJ, Müller B: Cis- and trans-acting gene regulation is associated with osteoarthritis. Am J Hum Genet. 2006, 78: 793-803. 10.1086/503849.
Gebhard PM, Soder S, Bau B, Aigner T: Down-regulation of the GTPase RhoB might be involved in the pre-apoptotic phenotype of osteoarthritic chondrocytes. Front Biosci. 2004, 9: 827-833. 10.2741/1254.
Davies CM, Guilak F, Weinberg JB, Fermor B: Reactive nitrogen and oxygen species in interleukin-1-mediated DNA damage associated with osteoarthritis. Osteoarthritis Cartilage. 2008, 16: 624-630. 10.1016/j.joca.2007.09.012.
Loughlin J, Meulenbelt I, Min J, Mustafa Z, Sinsheimer JS, Carr A, Slagboom PE: Genetic association analysis of RHOB and TXNDC3 in osteoarthritis. Am J Hum Genet. 2007, 80: 383-386. 10.1086/511443.
Kellgren JH, Lawrence JS: Radiological assessment of osteoarthrosis. Ann Rheum Dis. 1957, 16: 494-502.
DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7: 177-188. 10.1016/0197-2456(86)90046-2.
Bertram L, McQueen MB, Mullin K, Blacker D, Tanzi1 RE: Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet. 2007, 39: 17-23. 10.1038/ng1934.
Spector TD, Reneland RH, Mah S, Valdes AM, Hart DJ, Kammerer S, Langdown M, Hoyal CR, Atienza J, Doherty M, Rahman P, Nelson MR, Braun A: Association between a variation in LRCH1 and knee osteoarthritis: a genome-wide single-nucleotide polymorphism association study using DNA pooling. Arthritis Rheum. 2006, 54: 524-532. 10.1002/art.21624.
Jiang Q, Shi D, Nakajima M, Dai J, Wei J, Malizos KN, Qin J, Miyamoto Y, Kamatani N, Liu B, Tsezou A, Nakamura T, Ikegawa S: Lack of association of single nucleotide polymorphism in LRCH1 with knee osteoarthritis susceptibility. J Hum Genet. 2008, 53: 42-47. 10.1007/s10038-007-0216-4.
Mototani H, Mabuchi A, Saito S, Fujioka M, Iida A, Takatori Y, Kotani A, Kubo T, Nakamura K, Sekine A, Murakami Y, Tsunoda T, Notoya K, Nakamura Y, Ikegawa S: A functional single nucleotide polymorphism in the core promoter region of CALM1 is associated with hip osteoarthritis in Japanese. Hum Mol Genet. 2005, 14: 1009-1017. 10.1093/hmg/ddi093.
Loughlin J, Sinsheimer JS, Carr A, Chapman K: The CALM1 core promoter polymorphism is not associated with hip osteoarthritis in a United Kingdom Caucasian population. Osteoarthritis Cartilage. 2006, 14: 295-298. 10.1016/j.joca.2005.11.001.
Ikegawa S, Kawamura S, Takahashi A, Nakamura T, Kamatani N: Replication of association of the D-repeat polymorphism in asporin with osteoarthristis. Arthritis Res Ther. 2006, 8: 403-
Nakamura T, Shi D, Tzetis M, Rodriguez-Lopez J, Miyamoto Y, Tsezou A, Gonzalez A, Jiang Q, Kamatani N, Loughlin J, Ikegawa S: Meta-analysis of association between the ASPN D-repeat and osteoarthritis. Hum Mol Genet. 2007, 16: 1676-1681. 10.1093/hmg/ddm115.
Chapman K, Takahashi A, Meulenbelt I, Watson C, Rodriguez-Lopez J, Egli R, Tsezou A, Malizos KN, Kloppenburg M, Shi DQ, Southam L, Breggen van der R, Donn R, Qin J, Doherty M, Slagboom PE, Wallis G, Kamatani N, Jiang Q, Gonzalez A, Loughlin J, Ikegawa S: A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5' UTR of GDF5 with osteoarthritis susceptibility. Hum Mol Genet. 2008, 17: 1497-1504. 10.1093/hmg/ddn038.
This work was supported by National Nature Science Foundation of China (3057874; to DS and QJ), Programme of Technology Development of Nanjing (200603001; to DS and QJ) and Scientific Research Foundation for Young Doctors of Drum Tower Hospital (2007).
The authors declare that they have no competing interests.
All authors contributed to the final manuscript. In addition, DQS genotyped the samples and participated in the design and analysis of the study, and MN genotyped the Japanese samples. JD, JQ, NJ, YX, CY and JW evaluated the patients and genotyped Chinese samples, and TN helped with the meta-analysis. BL coordinated the study. QJ and SI supervised the whole study.
About this article
Cite this article
Shi, D., Nakamura, T., Nakajima, M. et al. Association of single-nucleotide polymorphisms in RHOB and TXNDC3 with knee osteoarthritis susceptibility: two case-control studies in East Asian populations and a meta-analysis. Arthritis Res Ther 10, R54 (2008). https://doi.org/10.1186/ar2423
- Control Individual
- Replication Study
- East Asian Population
- Healthy Control Individual
- Drum Tower Hospital