BioMed Research International
Volume 2019 (2019), Article ID 8298193, 10 pages
https://doi.org/10.1155/2019/8298193

Research Article

MMP2 and MMP10 Polymorphisms Are Related to Steroid-Induced Osteonecrosis of the Femoral Head among Chinese Han Population

Ye Tian,1,2 Feimeng An,1,2 Jiaqi Wang,1,2 Chang Liu,1,2 Huiqiang Wu,³ Yuju Cao,⁴ Jianzhong Wang,² and Guoqiang Wang²

¹Inner Mongolia Medical University, Hohhot, Inner Mongolia, China
²Department of Trauma Orthopedics, The Second Affiliated Hospital of Inner Mongolia Medical University, Hohhot, Inner Mongolia, China
³Inner Mongolia Autonomous Region Hospital of Traditional Chinese Medicine, Hohhot, Inner Mongolia, China
⁴Zhengzhou Traditional Chinese Medicine (TCM) Traumatology Hospital, Zhengzhou, Henan Province, China

Correspondence should be addressed to Jianzhong Wang; wangjianzhong0503@126.com and Guoqiang Wang; wangguoqianggq@163.com

Received 26 October 2018; Accepted 26 March 2019; Published 5 May 2019

Academic Editor: Yuan Yang

Copyright © 2019 Ye Tian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Steroid-induced osteonecrosis of the femoral head is a relatively serious condition which seriously reduces patient quality of life. However, the pathogenesis of steroid-induced ONFH is still unclear. In recent years, more scholars have found that the pathogenesis of steroid-induced ONFH is related to susceptibility factors such as MMPs/TIMPs system. The main purpose of this study is to investigate the correlation between MMP2 and MMP10 gene polymorphisms and steroid-induced ONFH in Chinese Han population. Methods. Six SNPs in MMP2 and two SNPs in MMP10 were genotyped using Agena MassARRAY RS1000 system from 286 patients of steroid-induced ONFH and in 309 healthy controls. The association between MMP2 and MMP10 polymorphisms and steroid-induced ONFH risk were estimated by the Chi-squared test, genetic model analysis, haplotype analysis, and stratification analysis. The relative risk was estimated by odd ratios (ORs) and 95% confidence intervals (CIs). Result. We found that the minor TG allele of rs470154 in MMP10 was associated with an increased risk of steroid-induced ONFH (OR = 1.45, 95% CI, 1.03 – 2.05, p = 0.032). In the genetic model analysis, we found that rs2241146 in MMP2 gene and rs470154 in MMP10 gene showed a statistically significant association with increased risk of steroid-induced ONFH. The six SNPs (rs470154, rs243866, rs243864, rs865094, rs11646643, and rs2241146) showed a statistically significant association with different clinical phenotypes. Conclusion. Our results verify that genetic polymorphisms of MMP2 and MMP10 contribute to steroid-induced ONFH susceptibility in the population of Chinese Han population, and our study provides new insights into the role that MMP2 and MMP10 plays in the mechanism of ONFH.

1. Introduction

Osteonecrosis of the femoral head (ONFH) is a debilitating bone disease, characterized by collapse of the femoral head and subsequent loss of hip joint function [1, 2]. ONFH frequently occurs in young and middle-aged people and seriously affects the quality of life for patients. At present, the number of new patients diagnosed with ONFH each year is 150,000 to 200,000 in China [3]. Widespread use of corticosteroids is widely believed to be one of the common causes of nontraumatic osteonecrosis. Steroid-induced ONFH accounts for 24.1% of all femoral head necrosis [4]. In clinical practice, we can often find that not all patients with a few months of steroid administration will develop steroid-induced ONFH, suggesting that genetic factors may confer susceptibility or resistance to steroid-induced ONFH. The certain genetic background has been shown related to ONFH in previous studies; single nucleotide polymorphism (SNP) is believed to be positively correlated with steroid-induced ONFH risk [5, 6].

Matrix metalloproteinases (MMPs), an enzyme family that is involved in the degradation of various proteins in the extracellular matrix (ECM), play an important part in the tissue remodeling process and physiological and pathological repair [7–9]. A complete set of MMPs produced by human cells is defined based on the complete human genome sequence. Recent genomic studies have shown that the MMP family has 24 distinct genetically encoded members, including MMP2 and MMP10 [10]. MMP2, encoded by MMP2 gene in chromosome 16q12.2, is a 72-kDa type IV collagenase. MMP2 is a rich enzyme secreted by osteoblasts and there is a significant negative correlation between MMP2 levels and bone mineral density [11]. Keiichi Inoue’s previous experiments demonstrated that MMP2 is important to the progress of forming and maintaining the osteocyte canalicular network and its formation is a determinant of bone remodeling and mineralization [12, 13]. In addition, studies have shown that osteolysis can be caused by abnormal expression of MMP2 [14]. In the course of ONFH, articular cartilage degenerates in the early stages of ONFH before significant radiological changes occur in the femoral head; ONFH articular cartilage degradation indicators include MMP2 [15]. MMP2 activates various cytokines and chemokines in inflammatory response [16]. Another gene MMP10, encoded by MMP10 gene in chromosome 11q22.2, is a secreted metalloproteinase involved in the physiological process of bone growth and wound healing [17]. In the occurrence of periodontitis, MMP10 may increase the risk of periodontitis; its mechanism is consistent with the local inflammatory response [18].

Studies have shown that the lack of functional MMP10 may be involved in the calcification process in vivo and in vitro VSMC, leading to a decrease in the degree of calcification. This result was consistent with reports of impaired atherosclerotic calcification in MMP2-deficient mice [17]. In addition, Duits et al. demonstrated higher plasma MMP2 and CSF MMP10 levels in AD patients compared to VaD patients and controls [19]. In previous study, investigator found the association between genetic polymorphisms of MMP2 and the risk of alcohol-induced osteonecrosis of the femoral head in the population of China [20]. Although MMP2 has been studied more extensively than MMP10, CSF MMP10, also known as lysostromin-2, has been linked to inflammation. However, there are no reports on the association between genetic polymorphisms of MMP2 and MMP10 and the risk of steroid-induced ONFH. Given the crucial role of MMP2 and MMP10 in the inflammatory response, bone metabolism, and bone loss, we speculated that the polymorphism of MMP2 and MMP10 might be associated with the risk of steroid-induced ONFH. Thus, we undertook the case-control study included 286 patients with steroid-induced ONFH and 309 healthy controls to elucidate the role of MMP2 and MMP10 gene polymorphisms in the occurrence of steroid-induced ONFH among Chinese Han population.

2. Materials and Methods

2.1. Study Participants

From September 2014 to January 2016, a total of 595 subjects included 286 patients with steroid-induced ONFH and 309 healthy controls were consecutively recruited from Zhengzhou Traditional Chinese Medicine (TCM) Traumatology Hospital. The controls were enrolled based on a physical examination, with no steroid-induced ONFH or other related diseases. All subjects were Han people living in or near Zhengzhou City.

The inclusion criteria of cases were as follows: ONFH was diagnosed based on MRI and CT scanning. Steroid-induced osteonecrosis was diagnosed based on a history of taking a mean daily dose of 16.6 mg or an equivalent maximum daily dose of 80 mg of prednisolone within 1 year [21]. The exclusion criteria of cases were as follows: patients did not meet the diagnostic criteria of steroid-induced ONFH; patients with traumatic ONFH, dislocation of the hip joint, or other hip diseases were excluded; patients drank more than 400 ml per week or had significant familial hereditary disease.

This study complied with the principles of the Declaration of Helsinki. The Ethical Committee of Zhengzhou Traditional Chinese Medicine Traumatology Hospital agreed to the study and all participants provided written informed consent.

2.2. SNPs Selection and Genotyping

Based on data from the database (dbSNP and GWASs), eight SNPs of MMP2 and MMP10 were selected. All SNPs we selected had minor allele frequencies > 5% in the global population from the 1000 Genome HapMap database. We collected blood samples into EDTA tubes after centrifugation at 2000 rpm for 10 minutes. Blood samples kept at -80°C for future experiments. Genomic DNA was obtained from whole blood using GoldMag extraction method (GoldMag, China) and then stored at -20°C. We evaluated DNA quantity using spectrophotometry (DU530UV/VIS spectrophotometer, Beckman Instruments, Fullerton, CA, USA). Agena MassARRAY Assay Design 4.0 Software was implemented to develop the Multiplexed SNP Mass EXTEND assay. Agena MassARRAY RS1000 system was carried out to measure the SNP genotypes by the manufacturer's protocol. Data management and analyses were performed using the Agena Bioscience Typer 4.0 software.

2.3. Statistical Analyses

Statistical analyses were performed using the SPSS 20.0 software (SPSS, Chicago, IL). The Chi-square test was used to check whether the genotype distribution matched Hardy–Weinberg equilibrium (HWE) in the control group. In this study, all the tests were two-tailed, and p values less than 0.05 were regarded as statistically significant. We used Chi-squared test to calculate the distribution differences of allele, genotype, and haplotype. We analyzed associations between the genotypes of MMP2 and MMP10 polymorphisms and the risk of steroid-induced ONFH through different genetic models, codominant, dominant, recessive, and log-additive models by computing odds ratios (ORs) and 95% confidence intervals (CIs) from logistic regression adjusted for age and gender. The stratification analysis was performed with clinical stages, hip lesions, and course condition. Finally, haploview software (version 4.2) was used to estimate the linkage disequilibrium (LD) and haplotype construction.

3. Result

A total of 286 patients with steroid-induced ONFH (mean age: 41.83 ± 13.12) and 309 healthy controls (mean age: 48.78 ± 8.30) were recruited for our study. Basic characteristics, including gender, age, and the clinical phenotypes (unilateral or bilateral hip lesions, >12 months or ≤12 months course, and II, III, IV clinical stages of steroid-induced ONFH) were listed in Table 1. Steroid-induced ONFH cases and healthy controls were matched in sex (p = 0.412). However, we found the cases and healthy controls were not matched for age (p < 0.05), so we adjusted the age factor in the subsequent analysis.

Table 1: Characteristics of the individuals in controls and steroid-induced ONFH patients.

Variables
Cases (n=286)
Controls (n=309)
p value

N (%)
N (%)

Age, years

Mean ± SD
41.83 ± 13.12
48.78 ± 8.30
< 0.001

≤48
203 (71.0)
170 (55.0)

>48
83 (29.0)
139 (45.0)

Gender

0.412

female
113 (39.5)
112 (36.2)

male
173 (60.5)
197 (63.8)

Clinical stages

Stage II
64 (22.4)

Stage III
126 (44.1)

Stage IV
96 (33.6)

Hip lesions

Unilateral
83 (29.0)

Bilateral
203 (71.0)

Course, months

>12
157 (45.1)

≤12
129 (54.9)

P < 0.05 indicates statistical significance.

Independent samples t test.

Two-sided Chi-squared test.

Table 2 summarizes basic information about the candidate in the eight SNPs. None of the SNPs that were excluded deviated from Hardy–Weinberg equilibrium (HWE). The χ² test was utilized to calculate the differences in frequency distribution of alleles between cases and controls. We found that rs470154 in MMP10 was linked to an increased risk of steroid-induced ONFH (OR = 1.45, 95% CI: 1.03 –2.05, p = 0.035).

Table 2: Basic information of candidate SNPs in this study.

SNP
Gene
Position
Alleles
MAF
p value
ORs
95% CI
p value

A/B
case
control
for HWE

rs243866
MMP2
16q12.2
A/G
0.112
0.173
0.329
0.79
0.56
1.13
0.194

rs243864
MMP2
16q12.2
G/T
0.112
0.140
0.342
0.78
0.55
1.10
0.151

rs865094
MMP2
16q12.2
G/A
0.290
0.266
1.000
1.13
0.87
1.45
0.357

rs11646643
MMP2
16q12.2
G/A
0.175
0.205
0.601
0.83
0.62
1.11
0.202

rs2241146
MMP2
16q12.2
A/G
0.234
0.213
0.611
1.13
0.86
1.48
0.388

rs9928731
MMP2
16q12.2
T/C
0.486
0.476
0.170
1.04
0.83
1.31
0.719

rs470154
MMP10
11q22.2
T/G
0.149
0.107
1.000
1.45
1.03
2.05
0.032*

rs17293607
MMP10
11q22.2
G/A
0.003
0.005
1.000
0.72
0.12
4.32
0.931

SNP: single nucleotide polymorphism, HWE: Hardy-Weinberg equilibrium, OR: odds ratio, 95% CI: 95% confidence interval, and MAF: minor allele frequency.

P values were calculated by exact test.

P values were calculated by Pearson Chi-squared test.

Subsequently, we evaluated the association between these SNPs and steroid-induced ONFH risk using four genetic models (codominant, dominant, recessive, and log-additive) by logistic-regression analysis. Table 3 showed that rs2241146 in the MMP2 gene was associated with an increased risk of steroid-induced ONFH in the recessive model before (OR = 2.15, 95% CI: 1.05–4.41, p = 0.032) and after (adjusted OR = 2.42, 95% CI: 1.14–5.14, p = 0.018) adjustment. Moreover, rs470154 in MMP10 polymorphism was correlated with an increased risk of steroid-induced ONFH in the log-additive model before adjustment for age and gender (OR = 1.44, 95% CI: 1.02–2.02, p = 0.035).

Table 3: Genotypic model analysis of relationship between SNPs and steroid-induced ONFH.

SNP ID
Model
Genotype
Controls
Cases
Without Adjustment
With Adjustment

OR (95% CI)
p value
OR (95% CI)
p value

rs470154
Codominant
G/G
245 (79.5%)
210 (73.4%)
1
0.067
1
0.120

G/T
60 (19.5%)
67 (23.4%)
1.30 (0.88-1.93)
1.27 (0.84-1.92)

T/T
3 (1%)
9 (3.1%)
3.50 (0.94-13.10)
3.37 (0.83-13.67)

Dominant
G/G
245 (79.5%)
210 (73.4%)
1
0.078
1
0.130

G/T-T/T
63 (20.4%)
76 (26.6%)
1.41 (0.96-2.06)
1.36 (0.91-2.04)

Recessive
G/G-G/T
305 (99%)
277 (96.8%)
1
0.055
1
0.083

T/T
3 (1%)
9 (3.1%)
3.30 (0.89-12.33)
3.20 (0.79-12.95)

Log-additive
- - -
- - -
- - -
1.44 (1.02-2.02)
0.035
1.40 (0.98-2.00)
0.066

rs2241146
Codominant
G/G
188 (61.2%)
175 (61.2%)
1
0.078
1
0.058

G/A
107 (34.9%)
88 (30.8%)
0.88 (0.62-1.25)
0.93 (0.65-1.34)

A/A
12 (3.9%)
23 (8%)
2.06 (0.99-4.26)
2.36 (1.10-5.07)

Dominant
G/G
188 (61.2%)
175 (61.2%)
1
0.990
1
0.710

G/A-A/A
119 (38.8%)
111 (38.8%)
1.00 (0.72-1.39)
1.07 (0.75-1.51)

Recessive
G/G-G/A
295 (96.1%)
263 (92%)
1
0.032
1
0.018

A/A
12 (3.9%)
23 (8%)
2.15 (1.05-4.41)
2.42 (1.14-5.14)

Log-additive
- - -
- - -
- - -
1.12 (0.86-1.46)
0.400
1.19 (0.90-1.58)
0.230

P values were calculated by Wald test by unconditional logistic regression adjusted for age and gender.

P < 0.05 indicates statistical significance.

Next, cases were stratified in terms of unilateral or bilateral hip lesions, > 12 months or ≤ 12 months course, and clinical stages. We further evaluated the effects of alleles and genotypes of SNPs on the risk of steroid-induced ONFH among the different subgroups (Table 4). In the allele and genotype models, the locus associated with steroid-induced risk of ONFH in patients with unilateral femoral head necrosis was rs243866 (adjusted OR = 0.40, 95% CI: 0.20–0.80, p = 0.007; p = 0.013), rs243864 (adjusted OR = 0.40, 95% CI: 0.20–0.78, p = 0.006; p = 0.010), and rs865094 (adjusted OR = 1.56, 95% CI: 1.08–2.24, p = 0.016; p = 0.029). We found that rs11646643 and rs2241146 were only associated with steroid-induced ONFH risk in patients with unilateral femoral head necrosis in the allele model. Among them, rs2241146 was associated with steroid-induced femoral head necrosis in the genotype model (unilateral, p = 0.038; aged > 12 months, p = 0.013). Moreover, in the allele model, rs470154 was associated with steroid-induced ONFH risk in patients with bilateral femoral head necrosis (adjusted OR = 1.50, 95% CI: 1.04–2.18, p = 0.031), clinical stage III (adjusted OR = 1.76, 95% CI: 1.17–2.67; p = 0.007), ≤ 12 months course (adjusted OR = 1.71, 95% CI: 1.13–2.59, p = 0.010). In the genotype model, it was related to clinical stage III (p = 0.015). In general, the SNPs with protective factors were rs243866, rs243864, and rs11646643. Separately, the SNPs with increase risk were rs470154, 865094, and rs2241146. In the case-control study, we performed a power analysis of selected MMP2 and MMP10 SNPs (Table 5).

Table 4: Stratified analysis for association of MMP2 and MMP10 gene polymorphism with the clinical phenotypes of steroid-induced ONFH.

SNPs
Subgroups
Genotype(n)
Allele(%)

AA/AB/BB
p
OR
95% CI
p

rs470154(T/G)
Controls
3 / 60 / 245

Case-Unilateral
4 / 15 / 64
0.094
1.34
0.81
2.23
0.258

Case-Bilateral
5 / 52 / 146
0.250
1.50
1.04
2.18
0.031

Case-Stage II
3 / 11 / 50
0.130
1.28
0.72
2.26
0.401

Case-Stage III
5 / 34 / 87
0.015
1.76
1.17
2.67
0.007

Case-Stage IV
1 / 22 / 73
0.820
1.19
0.72
1.96
0.492

Case > 12 months
5 / 31 / 121
0.190
1.25
0.83
1.90
0.290

Case ≤ 12 months
4 / 36 / 89
0.074
1.71
1.13
2.59
0.010

rs243866 (A/G)
Controls
8 / 68 / 231

Case-Unilateral
0 / 10 / 73
0.013
0.40
0.20
0.80
0.007

Case-Bilateral
3 / 48 / 152
0.470
0.97
0.67
1.40
0.862

Case-Stage II
1 / 10 / 53
0.400
0.65
0.35
1.23
0.187

Case-Stage III
1 / 21 / 104
0.240
0.63
0.39
1.03
0.064

Case-Stage IV
1 / 27 / 68
0.240
1.12
0.71
1.77
0.620

Case > 12 months
2 / 32 / 123
0.470
0.82
0.54
1.24
0.341

Case ≤ 12 months
1 / 26 / 102
0.240
0.77
0.49
1.21
0.255

rs243864 (G/T)
Controls
8 / 70 / 230

Case-Unilateral
0 / 10 / 73
0.010
0.40
0.20
0.78
0.006

Case-Bilateral
3 / 48 / 152
0.500
0.95
0.66
1.36
0.764

Case-Stage II
1 / 10 / 53
0.360
0.64
0.34
1.20
0.163

Case-Stage III
1 / 21 / 104
0.200
0.62
0.38
1.01
0.051

Case-Stage IV
1 / 27 / 68
0.270
1.10
0.69
1.73
0.692

Case > 12 months
2 / 32 / 123
0.470
0.80
0.53
1.21
0.286

Case ≤ 12 months
1 / 26 / 102
0.230
0.75
0.48
1.18
0.213

rs865094 (G/A)
Controls
22 / 120 / 166

Case-Unilateral
8 / 44 / 31
0.029
1.56
1.08
2.24
0.016

Case-Bilateral
9 / 88 /106
0.430
0.97
0.73
1.29
0.855

Case-Stage II
4 / 30 / 30
0.500
1.16
0.77
1.77
0.478

Case-Stage III
9 / 58 / 59
0.780
1.19
0.86
1.64
0.290

Case-Stage IV
4 / 92 / 48
0.370
1.02
0.71
1.47
0.900

Case > 12 months
9 / 70 / 78
0.750
1.07
0.79
1.45
0.649

Case ≤ 12 months
8 / 62 / 59
0.240
1.19
0.87
1.65
0.277

rs11646643 (G/A)
Controls
11 / 104 / 193

Case-Unilateral
1 / 18 / 64
0.062
0.53
0.32
0.88
0.014

Case-Bilateral
5 / 70 /127
0.320
0.96
0.70
1.31
0.800

Case-Stage II
2 / 22 / 40
0.640
0.99
0.62
1.59
0.971

Case-Stage III
2 / 38 / 86
0.360
0.78
0.53
1.14
0.200

Case-Stage IV
2 / 28 / 65
0.430
0.79
0.51
1.21
0.273

Case > 12 months
4 / 44 / 108
0.370
0.78
0.54
1.11
0.166

Case ≤ 12 months
2 / 44 / 83
0.320
0.89
0.61
1.29
0.532

rs2241146 (A/G)
Controls
12 / 107 / 188

Case-Unilateral
10 / 30 / 43
0.038
1.59
1.08
2.33
0.017

Case-Bilateral
13 / 58 / 132
0.230
0.96
0.71
1.31
0.804

Case-Stage II
6 / 14 / 44
0.076
0.94
0.59
1.51
0.797

Case-Stage III
7 / 44 / 75
0.650
1.10
0.78
1.57
0.587

Case-Stage IV
10 / 30 / 56
0.078
1.30
0.89
1.89
0.173

Case > 12 months
17 / 51 / 89
0.013
1.37
1.00
1.88
0.050

Case ≤ 12 months
6 / 37 / 86
0.710
0.86
0.60
1.25
0.435

rs9928731 (T/C)
Controls
63 / 166 / 78

Case-Unilateral
17 / 43 / 23
0.970
0.95
0.68
1.35
0.789

Case-Bilateral
47 / 107 / 49
0.860
1.08
0.84
1.39
0.542

Case-Stage II
15 / 29 / 20
0.550
0.94
0.64
1.38
0.763

Case-Stage III
28 / 72 / 26
0.600
1.14
0.85
1.53
0.387

Case-Stage IV
21 / 49 / 26
0.980
0.99
0.72
1.37
0.969

Case > 12 months
27 / 88 / 42
0.670
0.91
0.69
1.20
0.500

Case ≤ 12 months
37 / 62 / 30
0.270
1.23
0.92
1.64
0.164

P<0.05 indicates statistical significance.

P values were calculated by Wald test adjusted for age and gender.

P values were calculated by Pearson Chi-squared test.

Table 5: A power analysis of selected SNPs for MMP2 and MMP10 in case-control studies.

SNP_ID
Case n₁
Control n₂
Case A
Control A
Case p₁
Control p₂
p
zβ
power

rs243866
572
614
64
84
0.112
0.137
0.125
-0.665
0.253

rs243864
572
616
64
86
0.112
0.140
0.126
-0.525
0.300

rs865094
572
616
166
164
0.290
0.266
0.278
-1.037
0.150

rs11646643
570
616
100
126
0.175
0.205
0.191
-0.687
0.246

rs2241146
572
614
134
131
0.234
0.213
0.223
-1.095
0.137

rs9928731
572
614
278
292
0.486
0.476
0.481
-1.600
0.055

rs470154
572
616
85
66
0.149
0.107
0.127
0.183
0.573

rs17293607
572
618
2
3
0.003
0.005
0.004
-1.608
0.054

p<0.05 indicates statistical significance.

Finally, the parameter D′ was used to analyze the extent of linkage disequilibrium between SNPs. Then the haplotype LD block was determined according to the control group data (Figure 1). Two blocks were detected by haplotype analysis in the MMP2 gene (Figure 1), one block comprising rs243866 and rs243864 and the other one comprising rs11646643, rs2241146, and rs9928731. However, we did not find any association between the haplotype in the MMP2 and steroid-induced ONFH risk (Tables 6 and 7).

Table 6: The haplotype frequencies of MMP2 polymorphisms and their association with steroid-induced ONFH risk.

Haplotype
Freq
OR(95% CI)
p

rs11646643
rs2241146
rs9928731

A
G
T
0.481
1
-

A
A
C
0.220
1.13 (0.84 - 1.54)
0.420

G
G
C
0.188
0.84 (0.60 - 1.19)
0.330

A
G
C
0.108
0.86 (0.57 - 1.29)
0.470

p < 0.05 indicates statistical significance.

a= adjusted by gender and age.

Table 7: The haplotype frequencies of MMP2 polymorphisms and their association with steroid-induced ONFH risk.

Haplotype
Freq
OR(95% CI)
p

rs243866
rs243864

G
T
0.874
1
-

A
G
0.125
0.77 (0.53 - 1.10)
0.150

p < 0.05 indicates statistical significance.

a= adjusted by gender and age.

Figure 1: Linkage disequilibrium (LD) plots containing six SNPs from MMP2. Red squares display statistically significant associations between a pair of SNPs, as measured by D’; darker shades of red indicate higher D’.

4. Discussion

In this study, we found rs2241146 and rs470154 were significantly associated with an increased risk of steroid-induced ONFH. Six SNPs showed a statistically significant association with different clinical phenotypes (rs470154, rs243866, rs243864, rs865094, rs11646643, and rs2241146). To our knowledge, this is the first report to investigate an association between these loci and steroid-induced ONFH susceptibility.

MMP2 has been taken into account as a target candidate gene of genetic association studies for numerous human diseases. Robert Wojciechowski et al. [22] study showed that rs9928731 in MMP2 was significantly associated with ocular refraction in the AMISH families, which implicated MMP2 as regulator of nonpathologic refractive error. Ou Liu et al. [23] found that the MMP2 haplotype rs2241145-rs9928731 showed the significant association with Thoracic Aortic dissection (TAD), which may indicate the genetic predisposition to TAD in Chinese Han population. Rs9928731, rs11646643, and rs2241146 have a strong linkage between each pair of the three SNPs. In our study, however, there was no association between rs9928731 and steroid-induced ONFH. Inconsistent data observed between single-label analysis and haplotype analysis may indicate that rs9928731 was not independently associated with ONFH risk, but rather the SNPs in combination. Rs9928731 is located on chr16 of the MMP2 intron region between the sixth and seventh exons in MMP2, it does not represent a coding variant and located near the region where acetylation of histone H3 lysine 27 is significantly enriched, which is a reliable chromatin marker of active enhancers. Noncoding variants in these regions can result in damage to regulatory functions, and in some cases are supposed to cause disease phenotypes. In addition, epigenetic regulation of gene expression by histone modifications is cell-type specific. In the relevant region, H3 lysine 27 acetylation is highly enriched in human fibroblasts (MMP2 region) and human keratinocytes (MMP1 through MMP10 intergenic region); SNPs are also possible [22, 24].

Juiz et al. [25] indicated that rs243866 may be involved in the ECM remodeling. Additionally, previous studies showed that A allele carriers may be related to decreased risk of systolic heart failure in a Han Chinese population [26]. MMP2 rs243866 is located in 5 to a half palindromic potential estrogen receptor binding site. The 1575G allele acts as an enhancer in estrogen receptor-positive MCF-7 cells, while the 1575A allele reduces transcriptional activity. Gel migration assays confirmed that the binding of the estrogen receptor to this region was affected by differences in allelic expression [27]. Located in the GATA-1 site (CTATCT) of the promoter region of MMP2, the rs243866 polymorphism has binding sites for transcription factors AP-2, p53, Sp1, and Sp3. Protease expression is regulated by these transcription factors by the transcription rate of the gene [28]. For rs243864 polymorphism, a paper reported that the rs243864, rs243865, 243866 haplotype was statistically significant related to childhood obesity [29]. In addition, rs243864 is present in the gene promoter, which is a tag SNP of the rs243865 and rs243866 polymorphisms [30]. It has been reported that genetic associations in the MMP2 promoter may be caused by changes in MMP expression, because the three SNPs studied are functional SNPs that affect transcriptional activity.

Our study also found that MMP2 gene rs243866 and rs243864 may be the two protective SNPs of steroid-induced ONFH. The strong linkage between each pair of two SNPs (rs243866 and rs243864) confirmed that rs243866 and rs243864 could act not only alone but also coinitiate disease. In the Mexican mestizo population, rs243864 and rs11646643 of the MMP2 gene were associated with the high risk of COPD [28]. Our conclusion was that rs11646643 may be a low-risk factor of unilateral steroid-induced ONFH. And there was a strong linkage between each pair of the three SNPs (rs11646643, rs2241146, and rs9928731). But we did not find an association between steroid-induced ONFH and the above haplotype. In the study of Alicia Beeghly-Fadiel [31], rs2241146 may be unrelated to breast cancer risk. In this study, we found rs2241146 may be a high-risk factor of steroid-induced ONFH (especially unilateral ONFH).

Similarly, in our study, we found that rs470154 SNP in MMP10 gene had an increased risk steroid-induced of ONFH (especially bilateral ONFH, clinical stage III, and ≤ 12 months course). The other study showed that there is no association between rs470154 and the susceptibility to disc degeneration [32]. Information on specific genetic variants is inconsistent across different diseases, which can be explained by genetic differences between these diseases. It also illustrates the complexity of the mechanism involved in promoting SNP in ONFH. The differences observed here do not necessarily contradict previous reports. Further research is needed in populations with different ethnic and genetic backgrounds.

There were some important discoveries revealed in our study, but some limitations of this study should be considered when interpreting these results. First of all, our study did not include an analysis of biological functions, which will be crucial for elucidating the role of the gene in steroid-induced ONFH. Secondly, the sample size was relatively small, which might convert the positive findings into negative results. A larger case-control study is expected to circumvent those problems, which could make our conclusions more powerful. Although we obtained evidence for the association between MMP2 and MMP10 gene polymorphisms and steroid-induced ONFH risk, the pathogenesis of steroid-induced ONFH was still not completely understood. Further studies should be carried out to verify the present study in other ethnicity.

5. Conclusions

The MMP10 may be related to increased risk of steroid-induced ONFH. It is the first time that our results reveal an association between a MMP2 SNP at the rs2241146 locus and an increased risk of steroid-induced ONFH in Chinese Han population. And we found statistically significant differences in clinical phenotypes among the six SNPs (rs470154, rs243866, rs243864, rs865094, rs11646643, and rs2241146). From the genetic point of view to explore the pathogenesis of steroid-induced ONFH, early detection of susceptible population, timely intervention, and active prevention, the current treatment model can be transferred from postoperative surgery to predisease prevention. This study provides an important data source for the establishment of a bone necrosis gene bank in China. It can be used as a molecular marker for the diagnosis of osteonecrosis in the future, which not only benefits the majority of patients, but also has significant economic and social significance.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

Ethics committee approval was received for this study from Zhengzhou Traditional Chinese Medicine Traumatology Hospital.

Consent

All participants provided written informed consent.

Disclosure

Ye Tian and Feimeng An are co-first authors.

Conflicts of Interest

The authors declare there are no conflicts of interest.

Authors’ Contributions

Jiaqi Wang and Chang Liu conceived and designed the experiments. Ye Tian and Feimeng An wrote the paper. Huiqiang Wu contributed reagents/materials/analysis tools. Yuju Cao recruited samples. All authors reviewed the manuscript.

Acknowledgments

We gratefully acknowledge all of the participants, the clinicians, and all the hospital staff for their contributions to this study. We are also thankful to Zhengzhou TCM Traumatology Hospital, Zhengzhou, Henan Province, China, for their support. This work received the support of the National Natural Science Foundation of China (nos. 81160228, 81260284, and 81660378).

References

H. Kim, C. Cho, Y. Cho, S. Cho, K. Yoon, and K. Kim, “Significant associations of PAI-1 genetic polymorphisms with osteonecrosis of the femoral head,” BMC Musculoskeletal Disorders, vol. 12, article 160, 2011.

T.-H. Kim, J. M. Hong, H.-J. Kim, E. K. Park, and S.-Y. Kim, “Lack of association of MTHFR gene polymorphisms with the risk of osteonecrosis of the femoral head in a Korean population,” Molecules and Cells, vol. 29, no. 4, pp. 343–348, 2010.

Y. Song, Z. Du, M. Ren et al., “Association of gene variants of transcription factors PPARγ, RUNX2, Osterix genes and COL2A1, IGFBP3 genes with the development of osteonecrosis of the femoral head in Chinese population,” Bone, vol. 101, pp. 104–112, 2017.

L. Cui, Q. Zhuang, J. Lin et al., “Multicentric epidemiologic study on six thousand three hundred and ninety five cases of femoral head osteonecrosis in China,” International Orthopaedics, vol. 40, no. 2, pp. 267–276, 2016.

G. Hadjigeorgiou, E. Dardiotis, M. Dardioti, A. Karantanas, A. Dimitroulias, and K. Malizos, “Genetic association studies in osteonecrosis of the femoral head: Mini review of the literature,” Skeletal Radiology, vol. 37, no. 1, pp. 1–7, 2008.

J. Du, W. Liu, T. Jin et al., “A single-nucleotide polymorphism in MMP9 is associated with decreased risk of steroid-induced osteonecrosis of the femoral head,” Oncotarget , vol. 7, no. 42, pp. 68434–68441, 2016.

G. Murphy and A. J. Docherty, “The matrix metalloproteinases and their inhibitors.,” American Journal of Respiratory Cell and Molecular Biology, vol. 7, no. 2, pp. 120–125, 1992.

S. L. Parsons, S. A. Watson, P. D. Brown, H. M. Collins, and R. J. C. Steele, “Matrix metalloproteinases,” British Journal of Surgery, vol. 84, no. 2, pp. 160–166, 1997.

K. Białkowska, W. Marciniak, M. Muszyńska et al., “Association of zinc level and polymorphism in MMP-7 gene with prostate cancer in Polish population,” PLoS ONE, vol. 13, no. 7, Article ID e0201065, 2018.

C. Kapoor, S. Vaidya, V. Wadhwan, Hitesh, G. Kaur, and A. Pathak, “Seesaw of matrix metalloproteinases (MMPs),” Journal of Cancer Research and Therapeutics, vol. 12, no. 1, pp. 28–35, 2016.

J. S. Nyman, C. C. Lynch, D. S. Perrien et al., “Differential effects between the loss of MMP-2 and MMP-9 on structural and tissue-level properties of bone,” Journal of Bone and Mineral Research, vol. 26, no. 6, pp. 1252–1260, 2011.

K. Inoue, Y. Mikuni-Takagaki, K. Oikawa et al., “A crucial role for matrix metalloproteinase 2 in osteocytic canalicular formation and bone metabolism,” The Journal of Biological Chemistry, vol. 281, no. 44, pp. 33814–33824, 2006.

S. Grässel, J. Beckmann, B. Rath, M. Vogel, J. Grifka, and M. Tingart, “Expression profile of matrix metalloproteinase-2 and -9 and their endogenous tissue inhibitors in osteonecrotic femoral heads,” International Journal of Molecular Medicine, vol. 26, no. 1, pp. 127–133, 2010.

J. A. Martignetti, A. A. Aqeel, W. A. Sewairi et al., “Mutation of the matrix metalloproteinase 2 gene (MMP2) causes a multicentric osteolysis and arthritis syndrome,” Nature Genetics, vol. 28, no. 3, pp. 261–265, 2001.

Y. Han, M. Si, Y. Zhao et al., “Progranulin protects against osteonecrosis of the femoral head by activating ERK1/2 pathway,” Inflammation, vol. 40, no. 3, pp. 946–955, 2017.

M. Onishi, T. Kobayashi, C. N. D'Alessandro-Gabazza et al., “Mice overexpressing latent matrix metalloproteinase-2 develop lung emphysema after short-term exposure to cigarette smoke extract,” Biochemical and Biophysical Research Communications, vol. 497, no. 1, pp. 332–338, 2018.

A. Purroy, C. Roncal, J. Orbe et al., “Matrix metalloproteinase-10 deficiency delays atherosclerosis progression and plaque calcification,” Atherosclerosis, vol. 278, pp. 124–134, 2018.

S. N. Pandruvada, O. A. Gonzalez, S. Kirakodu et al., “Bone biology-related gingival transcriptome in ageing and periodontitis in non-human primates,” Journal of Clinical Periodontology, vol. 43, no. 5, pp. 408–417, 2016.

F. H. Duits, M. Hernandez-Guillamon, J. Montaner et al., “Matrix metalloproteinases in alzheimer's disease and concurrent cerebral microbleeds,” Journal of Alzheimer's Disease, vol. 48, no. 3, pp. 711–720, 2015.

Y. Yu, Z. Xie, J. Wang et al., “Single-nucleotide polymorphisms of MMP2 in MMP/TIMP pathways associated with the risk of alcohol-induced osteonecrosis of the femoral head in Chinese males: A case-control study,” Medicine (United States), vol. 95, no. 49, p. e5407, 2016.

Q. Wang, Q. Yang, G. Chen et al., “LncRNA expression profiling of BMSCs in osteonecrosis of the femoral head associated with increased adipogenic and decreased osteogenic differentiation,” Scientific Reports, vol. 8, no. 1, 2018.

R. Wojciechowski, J. E. Bailey-Wilson, and D. Stambolian, “Association of matrix metalloproteinase gene polymorphisms with refractive error in amish and ashkenazi families,” Investigative Ophthalmology & Visual Science, vol. 51, no. 10, pp. 4989–4995, 2010.

O. Liu, J. Li, Y. Xin et al., “Association of MMP-2 gene haplotypes with thoracic aortic dissection in chinese han population,” BMC Cardiovascular Disorders, vol. 16, no. 1, p. 11, 2016.

R. Wojciechowski, S. S. Yee, C. L. Simpson, J. E. Bailey-Wilson, and D. Stambolian, “Matrix metalloproteinases and educational attainment in refractive error: Evidence of gene-environment interactions in the age-related eye disease study,” Ophthalmology, vol. 120, no. 2, pp. 298–305, 2013.

N. A. Juiz, N. M. Cayo, M. Burgos et al., “Human polymorphisms in placentally expressed genes and their association with susceptibility to congenital trypanosoma cruzi infection,” The Journal of Infectious Diseases, vol. 213, no. 8, pp. 1299–1306, 2016.

Y. Hua, N. Wu, Y. Yang et al., “Polymorphisms of MMP-2 gene are associated with systolic heart failure risk in han chinese,” The American Journal of the Medical Sciences, vol. 337, no. 5, pp. 344–348, 2009.

S. Harendza, D. H. Lovett, U. Panzer, Z. Lukacs, P. Kühnl, and R. A. K. Stahl, “Linked common polymorphisms in the gelatinase a promoter are associated with diminished transcriptional response to estrogen and genetic fitness,” The Journal of Biological Chemistry, vol. 278, no. 23, pp. 20490–20499, 2003.

J. Hernández-Montoya, J. Pérez-Ramos, M. Montaño et al., “Genetic polymorphisms of matrix metalloproteinases and protein levels in chronic obstructive pulmonary disease in a Mexican population,” Biomarkers in Medicine, vol. 9, no. 10, pp. 979–988, 2015.

A. R. Morgan, D. Y. Han, J. M. Thompson, E. A. Mitchell, and L. R. Ferguson, “Analysis of MMP2 promoter polymorphisms in childhood obesity,” BMC Research Notes, vol. 4, no. 1, 2011.

I. Ito, S. Nagai, T. Handa et al., “Matrix metalloproteinase-9 promoter polymorphism associated with upper lung dominant emphysema,” American Journal of Respiratory and Critical Care Medicine, vol. 172, no. 11, pp. 1378–1382, 2005.

A. Beeghly-Fadiel, W. Lu, J.-R. Long et al., “Matrix metalloproteinase-2 polymorphisms and breast cancer susceptibility,” Cancer Epidemiol Biomarkers Prev, vol. 18, no. 6, pp. 1770–1776, 2009.

A. Kelempisioti, P. J. Eskola, A. Okuloff et al., “Genetic susceptibility of intervertebral disc degeneration among young Finnish adults,” BMC Medical Genetics, vol. 12, no. 1, 2011.