The association between rs12901499 polymorphism in SMAD3 gene and risk of osteoarthritis: a meta-analysis

Abstract

Purpose

This study was conducted to assess and synthesize the current evidence on the association between rs12901499 polymorphism in SMAD3 gene and risk of osteoarthritis (OA).

Materials and methods

Four electronic databases, including PubMed, Embase, ISI Web of Science, and CENTRAL were systematically searched for potential studies. Summary odds ratio and corresponding 95% CI were calculated to evaluate the association. Risk of bias was assessed through the Newcastle–Ottawa Scale. Subgroups and sensitivity analyses were performed using the RevMan 5.3 software. Publication bias was evaluated by Egger’s and Begg’s tests. Power analysis was conducted using the Power and Sample Size Calculation program.

Results

Eight case–control studies containing 5,625 patients with OA and 5,600 healthy controls were obtained for the meta-analysis. After excluding cohorts with inadequate power, the pooled data supported that G allele carriers of rs12901499 had a significantly increased risk of OA (odds ratio 1.31, 95% CI: 1.21 to 1.43, P<0.00001). When stratified by OA site and ethnicity, the association remained statistically significant.

Conclusion

The combined results evidently supported that rs12901499 polymorphism in SMAD3 gene is significantly associated with OA vulnerability across both Caucasian and Asian populations.

Keywords:

• osteoarthritis
• polymorphism
• SMAD3
• rs12901499
• meta-analysis

Introduction

Osteoarthritis (OA), a disease of the joint as an organ, is one of the most common degenerative musculoskeletal disorders predominantly affecting synovial joints, such as knees, hips, spines, and hands.¹ An individual joint could be involved, but in more cases, several joints are affected concomitantly. Disease progression of OA is typically slow, and it might take years or even decades to develop. Therefore, the dangerousness of OA always eludes the attention of patients until an advanced stage of this disorder. It is reported that ~240,000,000 elderly population suffers from OA around the world.² With the prolongation of human lifespan, activity incapacity caused by knee OA alone is greater than any other medical conditions among population aged >65 years. The highly prevalent and attendant disability of OA has a formidable impact on both individuals and society.^3,4 While, up to date, no effective treatment strategy for this common disease is still available. There is a great necessity to look for novel strategies to decrease the influence of this disturbing condition.

The underlying etiopathogenesis of OA remains largely indeterminate. Epidemiological studies have demonstrated that OA is driven by a host of factors, including age,⁵ gender,⁶ body mass index,⁷ as well as genetics.⁸ Involvement of genetic components in OA has attracted significant interest during the past 2 decades. In reality, researchers have observed familial aggregation phenomenon of OA for a long time.⁹ Evidence based on twin studies suggested an estimated heritability of 37%–65% of OA in different skeletal sites.¹⁰ Considerable single-nucleotide polymorphisms (SNPs) in candidates such as the estrogen-receptor,^11,12 calmodulin,^13,14 vitamin D receptor,¹⁵ ADAM12,¹⁶ and growth differentiation factor 5 genes¹⁷ have been reported to be linked to OA.

SMAD3 belongs to SMAD family, which consists of at least 8 different SMADs in human. It is an intracellular mediator molecule regulating the extracellular transforming growth factor beta (TGF-β) signals to transduce into nucleus. TGF-β/SMAD3 signaling is crucial to maintain the integrity of articular cartilage.¹⁸ SMAD3 gene-deficient mice are characterized with hypertrophic chondrocytes, which induce the progressive loss of cartilage accompanied by formation of osteophytes in joints.¹⁸ A clinical study led by Yao et al found that the SMAD3 gene mutations were associated with the pathogenesis of human OA.¹⁹

Albeit several reports have been published in an effort to demonstrate the correlation between rs12901499 polymorphism in SMAD3 gene and OA predisposition among different populations, it is quite regrettable that no convincing evidence is yet available, which might be owing to reasons such as limited sample size, inadequate statistical power, ethnicity specificity, heterogeneity of participants. Therefore, we turn to exploit a meta-analysis to synthesize the existing evidence on the association between rs12901499 polymorphism in SMAD3 gene and risk of OA.

Materials and methods

This systematic review and meta-analysis were conducted following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines.²⁰

Literature search strategy

Literature search of the study was conducted on PubMed, Embase, ISI Web of Science, and Cochrane Library (CENTRAL) 4 electronic databases, with all articles published previous to January 1, 2018. A combination of Medical Subject Headings (MeSH) together with free terms was applied to retrieve all the potentially eligible publications without language restriction. The following search string was employed: (Single Nucleotide Polymorphism or SNP or SNPs or “Polymorphism, Single Nucleotide” [Mesh]) and SMAD3 and (osteoarthritis or OA or osteoarthritides or osteoarthrosis or “Osteoarthritis” [Mesh]). Bibliographic lists of relevant reviews and full-text articles were also manually searched to capture more possible studies.

Inclusion and exclusion criteria

Studies meeting the following criteria were included for our meta-analysis: 1) studies examining the correlation between the rs12901499 polymorphism in SMAD3 and OA susceptibility that have been published; 2) case–control designed study; 3) studies with sufficient data to calculate the odds ratio (OR) with its 95% CI. Correspondingly, studies satisfying the following criteria were excluded: 1) animal study, review, expert opinion, conference abstract, case report, or case series; 2) data overlapping with previous publications; 3) studies failed to provide sufficient data to evaluate OR and the associated 95% CI.

Methodological quality assessment

Methodological quality of eligible studies was assessed independently by the first 2 reviewers using the Newcastle–Ottawa Scale (NOS) for observational studies.²¹ A “star system” has been developed in which a study is judged on 3 broad perspectives, including the selection of the study groups, the comparability of the groups, and the ascertainment of either the exposure or outcome of interest. Discrepancies between reviewers were settled through discussion until consensus was reached, a third reviewer could be consulted when necessary (W-BS).

Data extraction

Two independent investigators (S-TG and Z-TL) conducted the data extraction of all the included studies independently. A standardized data extraction list was utilized to collect information, including first author, year of the publication, country, sample size, ethnicity of enrolled subjects, site of OA, risk allele percentage of cases and controls, the calculated OR and 95% CI of individual cohort, and P-value for Hardy–Weinberg equilibrium (HWE) test of control subjects in each study. In the event of any discrepancy through the process, 2 reviewers reinvestigated the article together and discussed with each other until coming to an agreement.

Statistical analysis

To test the potential association between SMAD3 gene polymorphism (rs12901499) and risk of OA, the estimated allelic effect of G (the variant allele) versus A (the common allele) was indicated as OR and corresponding 95% CI using the allele frequency data extracted from included studies. Heterogeneity across included studies was estimated using a Q-test and the Higgins I² test, where P>0.1 and I²<50% indicated acceptable heterogeneity. I² values of 0%, 25%, 50%, and 75% were defined as no, low, moderate, and high heterogeneity, respectively. In our present study, we combined the OR of each study using a random-effect model because it accommodates diversity between studies and is preferable in the presence or anticipation of inter-study variances.²²

Prior to the meta-analysis, power analysis was performed using the Power and Sample Size Calculation program (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize) to evaluate whether the selected studies could offer adequate power to detect the association between SMAD3 polymorphism and OA at a level of significance of 0.05. Studies that failed to provide adequate <0.50) were not included in the meta-analysis. Regardless of the heterogeneity detected between included studies, we performed subgroup-analysis by ethnicity and site of OA. The leave-one-out sensitivity analysis was conducted by removing each included study at a time and reevaluating the resulting effect on pooled results to test the robustness of summary risk estimate. Begg’s rank correlation test and Egger’s linear regression test using Stata version 12.0 (Stata Corp LP, College Station, TX, USA) were combined to assess the publication bias, which was considered present with P<0.05.²³ Publication bias was also assessed by visually inspecting funnel plots, a visual asymmetry suggested possible publication bias. The forest and funnel plots were generated via RevMan 5.3 software (The Nordic Cochrane Centre, Copenhagen, Denmark).

Results

Literature search

The comprehensive literature retrieval of 4 electronic databases totally yielded 48 potentially relevant records comprising 16 from PubMed, 19 from Embase, 13 from ISI Web of Science, but none from CENTRAL. After the first scan, 20 duplicate records were deleted. Of the remaining 28 records, a further 18 citations were eliminated after title and abstract screening. Among the rest 10 articles, 2 unrelated were removed after assessment in the light of predefined inclusion and exclusion criteria. Ultimately, 8 studies^24–31 were incorporated into the qualitative synthesis. The process of literature selection process is presented in Figure 1.

Figure 1 Flow chart of literature search and screen.

Main characteristics

The main characteristics of included studies are displayed in Table 1. A total of 8 case–control studies with 16 cohorts involving OA of the knee, hip, hand, spine, as well as temporomandibular joint (TMJ) were identified and included in the meta-analysis. All the included articles were published between 2010 and 2018. The sample size of each study ranged from 240 to 6,100. Six studies^{24,26,27,29–31} were performed among Asian population, and the other 2^25,28 across Caucasian population. The study by Valdes et al²⁸ consisted of 5 independent studies involving 8 cohorts with knee or hip OA patients in case groups, and the study by Jiang et al²⁴ comprised 2 cohorts with hand or knee OA patients. Liva et al study²⁵ was about spine OA, and Xiao et al study²⁹ was on TMJ OA. All these included cohorts conformed to HWE with the exception of Zhong et al.³¹ Of note, the statistical powers of these studies varied from 0.11 to 1.00, and several cohorts lacked adequate power (detailed data presented in Table 1). On the basis of the NOS, each study received no less than 5 stars for methodological quality assessment, as shown in Table 2.

Table 1 Main characteristics of included studies

Study	Country	Case/control	Ethnicity	Site of OA	G allele (%)		Association		HWE	Power
					Case	Control	OR	95% CI
Jiang et al,^Citation24 2013 (1)*	China	111/236	Asian	Knee	45.10	29.32	1.71	1.22 to 2.40	0.41	0.97
Jiang et al,^Citation24 2013 (2)*	China	121/236	Asian	Hand	44.59	29.32	1.68	1.21 to 2.33	0.41	0.97
Liva et al,^Citation25 2017	Greece	258/243	Caucasian	Spine	55.80	53.50	1.10	0.86 to 1.41	>0.05	0.11
Sharma et al,^Citation26 2017	India	450/458	Asian	Knee	48.78	43.89	1.22	1.01 to 1.46	>0.05	0.55
Su et al,^Citation27 2015	Taiwan	545/468	Asian	Knee	48.81	50.85	0.92	0.77 to 1.10	0.86	0.15
Valdes et al,^Citation28 2010 (1)**	UK	313/520	Caucasian	Knee	60.70	52.98	1.37	1.12 to 1.68	>0.40	0.87
Valdes et al,^Citation28 2010 (2)**	UK	214/520	Caucasian	Hip	59.58	52.98	1.31	1.12 to 1.68	>0.40	0.64
Valdes et al,^Citation28 2010 (3)**	UK	1,083/733	Caucasian	Knee	56.69	52.39	1.19	1.04 to 1.36	>0.40	0.72
Valdes et al,^Citation28 2010 (4)**	UK	979/733	Caucasian	Hip	56.59	52.39	1.18	1.03 to 1.36	>0.40	0.71
Valdes et al,^Citation28 2010 (5)**	UK	167/867	Caucasian	Knee	59.88	55.07	1.22	0.96 to 1.55	>0.40	0.37
Valdes et al,^Citation28 2010 (6)**	Estonia	68/449	Caucasian	Knee	58.09	46.44	1.60	1.11 to 2.30	>0.40	0.72
Valdes et al,^Citation28 2010 (7)**	UK	257/488	Caucasian	Knee	54.47	51.13	1.14	0.92 to 1.42	>0.40	0.23
Valdes et al,^Citation28 2010 (8)**	UK	95/488	Caucasian	Hip	56.84	51.13	1.26	0.92 to 1.72	>0.40	0.30
Xiao et al,^Citation29 2015	China	114/126	Asian	TMJ	49.56	38.49	1.57	1.09 to 2.26	0.17	0.69
Zhang et al,^Citation30 2018	China	350/400	Asian	Knee	50.71	53.00	0.90	0.74 to 1.11	0.82	0.14
Zhong et al,^Citation31 2018	China	500/1,080	Asian	Hip	78.00	69.91	1.53	1.28 to 1.82	<0.01	1.00

Notes:

* This study comprised two cohorts: one with knee OA(1) and one with hand OA(2).

** This study comprised eight cohorts: discovery sets performed in Nottingham (1) and (2); replication sets performed in Nottingham (3) and (4); replication set performed in Hertfordshire (5); replication set performed in Estonia (6); and replication sets performed in Chingford (7) and (8).

Abbreviations: HWE, Hardy–Weinberg Equilibrium; OA, osteoarthritis; OR, odds ratio; TMJ, temporomandibular joint.

Table 2 Quality assessment of included studies

Item/study	Jiang et al,^Citation24 2013	Liva et al,^Citation25 2017	Sharma et al,^Citation26 2017	Su et al,^Citation27 2015	Valdes et al,^Citation28 2010	Xiao et al,^Citation29 2015	Zhang et al,^Citation30 2018	Zhong et al,^Citation31 2018
Adequate definition of cases	*	*	*	*	*	*	*	*
Representativeness of cases	–	–	*	–	–	–	–	–
Selection of control subjects	–	–	–	–	–	–	–	–
Definition of control subjects	*	*	*	*	*	*	*	*
Control for important factor	–	–	**	–	–	–	–	–
or additional factor
Exposure assessment	*	*	*	*	*	*	*	*
Same method of ascertainment	*	*	*	*	*	*	*	*
for all subjects
Non-response rate	*	*	*	*	*	*	*	*

Notes: A study could be awarded a maximum of 1 star for each item except for the item “Control for important factor or additional factor”. The definition/explanation of each column of the Newcastle–Ottawa Scale is available from http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp.

Meta-analysis and subgroup-analysis

After removing the cohorts with inadequate power, the pooled data suggested a significant correlation between rs12901499 polymorphism and risk of OA in the overall population with a modest heterogeneity (G versus A: OR 1.31, 95% CI: 1.21 to 1.43<0.00001; I²=29%; Figure 2). We performed subgroup-analysis by site of OA (knee, hip, hand, and TMJ) and ethnicity (Asian and Caucasian) to reappraise the association in different subgroups. Regarding OA site, subgroups of knee OA (OR 1.32, 95% CI: 1.17 to 1.48, P<0.00001; I²=35%), hip OA (OR 1.22, 95% CI: 1.08 to 1.37, P=0.001; I²=0), hand OA (OR 1.68, 95% CI: 1.21 to 2.33, P=0.002), and TMJ OA (OR 1.57, 95% CI: 1.09 to 2.26, P=0.01) all conformed to the association (Figure 3). In terms of ethnicity, there still existed a statistically significant association of rs12901499 polymorphism and OA in population of Asians (OR 1.47, 95% CI: 1.22 to 1.78, P<0.0001; I²=42%) and Caucasians (OR 1.24, 95% CI: 1.15 to 1.35, P<0.00001; I²=0), as shown in Figure 4.

Figure 2 Forest plot of rs12904199 in SMAD3 gene and risk of osteoarthritis using a random effect model.

Notes: *This study comprised two cohorts: one with knee OA(1) and one with hand OA(2). **This study comprised eight cohorts: discovery sets performed in Nottingham (1) and (2); replication sets performed in Nottingham (3); replication set performed in Hertfordshire (5); replication set performed in Estonia (6); and replication sets performed in Chingford (7) and (8); not all cohorts are shown in this figure.
Abbreviations: M–H, Mantel–Haenszel; OA, osteoarthritis; OR, odds ratio.

Figure 3 Forest plot of rs12904199 in SMAD3 gene and risk of osteoarthritis using a random effect model: subgroup-analysis by osteoarthritis site.

Notes: *This study comprised two cohorts: one with knee OA(1) and one with hand OA(2). **This study comprised eight cohorts: discovery sets performed in Nottingham (1) and (2); replication sets performed in Nottingham (3); replication set performed in Hertfordshire (5); replication set performed in Estonia (6); and replication sets performed in Chingford (7) and (8); not all cohorts are shown in this figure.
Abbreviations: M–H, Mantel–Haenszel; OA, osteoarthritis; OR, odds ratio; TMJ, temporomandibular joint.

Figure 4 Forest plot of rs12904199 in SMAD3 gene and risk of osteoarthritis using a random effect model: subgroup-analysis by ethnicity.

Notes: *This study comprised two cohorts: one with knee OA(1) and one with hand OA(2). **This study comprised eight cohorts: discovery sets performed in Nottingham (1) and (2); replication sets performed in Nottingham (3); replication set performed in Hertfordshire (5); replication set performed in Estonia (6); and replication sets performed in Chingford (7) and (8); not all cohorts are shown in this figure.
Abbreviations: M–H, Mantel–Haenszel; OA, osteoarthritis; OR, odds ratio.

Sensitive analysis and publication bias

Sensitivity analyses were tested by removing an individual dataset sequentially. But the pooled results did not change remarkably after the removal of any dataset, which suggested the stability and credibility of the overall estimate of risk (detailed data not shown). For publication bias investigation, there was an obvious symmetry presented in the funnel plot (Figure 5). Begg’s rank correlation test and Egger’s linear regression test also confirmed significant publication bias (Begg’s test: z=2.40, P=0.016; Egger’s test: t=8.02, P<0.001).

Figure 5 Funnel plot of rs12901499 in SMAD3 gene and risk of osteoarthritis.

Abbreviations: OR, odds ratio; SE, standard error.

Discussion

Despite identification of the accurate risk factors for OA might greatly benefit earlier treatment in addition to preventive strategy, the precise etiology of OA remains unanswered. Researches of OA have focused on multiple areas and acknowledged the complexity of this disorder. Intrinsic factors in synergy with extrinsic factors are considered to be the potential etiology. As an intrinsic factor, SMAD3 gene has attracted a great deal of attention over the recent years. The association between rs12901499 polymorphism in SMAD3 gene and OA vulnerability has been documented in several studies with conflicting outcomes. To overcome the drawbacks of limited sample size and ethnicity specificity of an individual study, we designed the present meta-analysis to investigate the underlying correlation of rs12901499 and OA in populations from different ethnicities. The pooled data supported G allele carriers of rs12901499 had a significantly increased risk of OA (OR 1.31, 95% CI: 1.21 to 1.43, P<0.00001). Likewise, the association remained significant regardless of the site of OA or any ethnicity of the included populations, which suggested the convincing evidence that rs12901499 polymorphism in SMAD3 gene contributed to OA predisposition.

SMADs are a cluster of proteins sharing similar structure but having distinct roles in signaling.³² Of them, SMAD3 has been popularly investigated. A large body of evidence has established that changes in the normal structure of SMAD3 molecule implicated in numerous physiological and disease processes.³³ SMAD3 gene is located in chromosome 15q22.33, a region, which has been reported to be linked with OA.³⁴ The expression of SMAD3 gene plays a key role in the TGF-β pathway which is critical to TGF-β signal transduction.^32,35 Reports have indicated that the TGF-β signaling promotes the metabolism of chondrocytes and assists in the maintenance of the normal phenotype of chondrocytes by regulating phosphorylation of SMAD3. In articular cartilage tissues, phosphorylated SMAD3 can protect chondrocyte phenotype from deleterious hypertrophy and production of catabolic enzymes.¹⁸ Li et al found that loss of SMAD3 could give rise to accelerated maturation of chondrocytes, the TGF-β signaling effects were as well decreased in SMAD3^−/− chondrocytes.³⁶ Van de Laar et al reported that mutations in SMAD3 gene could lead to serious aortic aneurysms and dissections accompanied by early-onset OA.³⁷ All these investigations supported the importance of SMAD3 gene to OA, and targeting of SMAD3 may provide a novel strategy for the treatment of OA.

It should be pointed out that the study led by Valdes et al²⁸ also involved a meta-analysis suggesting a significant association between rs12901499 polymorphism and OA susceptibility across Caucasians. Though the outcome of the present study was quite in line with the former, there still existed several notable differences. First of all, we attempted to conduct a more comprehensive research to find the underlying association between rs12901499 and OA. Hence, we incorporated several latest published studies with populations beyond Caucasians and conducted subgroup analysis by ethnicity to investigate whether there existed an ethnicity-specific effect of the association. Besides, for each of the included studies, we performed quality assessment with the NOS, which enabled us to judge potential risk of bias as well as grade the level of evidence. After that, considering genetic association studies might be overrepresented by underpowered studies, which indicated a significant association, the power calculation was employed to minimize the possibility of false negatives in the present study, and studies with low power were excluded from the final data synthesis. Last but not least, we carried out sensitivity analysis to validate stability and robustness of the outcome. As a result, all these advantages strongly support us to guarantee a more accurate and reliable conclusion.

Although the present meta-analysis has collected currently available data and provided a more comprehensive evaluation and precise evidence on the association between rs12901499 polymorphism and OA susceptibility, there were several limitations that could not be neglected. First, we only looked into the role of a single loci polymorphism in SMAD3 gene. But as with other complex disorders, OA is modulated by multiple loci polymorphisms and a range of genes with synergetic effects. Therefore, the effects of loci–loci and gene–gene interactions required to be addressed. Second, only 8 articles in English from 4 databanks were retrieved for the meta-analysis and the funnel plot suggested a significant publication bias, potential relevant articles published in other languages and databases might have been missed, which could lead to inclusion criteria bias.³⁸ Third, there was only 1 cohort focused on OA of the hand, spine, as well as TMJ in our study. Therefore, more investigations concerning these sites of OA are necessary to validate the association. Finally, despite the meta-analysis only incorporated case–control designed studies, no study had a satisfying match on age, sex, body mass index, or other confounding factors, which might have caused serious confounding bias in the present study.

Like dozens of other diseases, investigations focused on association of SNPs and OA have been going on for over 2 decades. A growing number of novel loci have been determined, and more are being discovered, which tend to provide new hints and shed light on the underlying mechanisms that affect OA susceptibility. Nevertheless, deviating from original expectation of researchers, most SNPs appear to have only a slight effect on OA, and almost any SNP on hand only seem to increase the risk of OA. In reality, present findings and advances in the genetic basic of OA are nothing more than initial insights into the uncharted territory, and are still far from application in clinical practice of early diagnosis and personalized therapy of OA. As the reach of science always turns ahead in the twists, through unremitting endeavor of researchers, we believe current and future findings will contribute to drive the prevention and treatment of OA through screening individuals if a first-degree relative has a history of OA, or by developing new drugs that target the specific loci.

Conclusion

Collectively, the combined result evidently supported rs12901499 polymorphism in SMAD3 gene was significantly associated with OA vulnerability across both Caucasian and Asian populations. However, the gene polymorphisms identified hitherto were associated with only a fraction of increase in OA risk, and therefore, the determinate etiology of OA may continue to elude researchers. Further work is still encouraged to reveal the etiopathogenesis of OA.

Disclosure

The authors report no conflicts of interest in this work.