BioMed Research International
Volume 2016 (2016), Article ID 1627143, 6 pages
http://dx.doi.org/10.1155/2016/1627143

Research Article

Matrix Metalloproteinase-9 −1562C/T Gene Polymorphism Is Associated with Diabetic Nephropathy

Shufen Feng,¹ Gang Ye,¹ Shi Bai,² Hongcheng Wei,¹ Xueling Liao,³ and Lu Li¹

¹Department of Gastroenterology, The First Affiliated Hospital of Jinan University, Guangzhou 510630, China
²Clinical Laboratory, People’s Hospital of New District Longhua, Shenzhen 518109, China
³Department of Nephrology, Affiliated Hospital of Guilin Medical College, Guangxi 541001, China

Received 17 May 2016; Accepted 25 July 2016

Academic Editor: Konstantinos Kantartzis

Copyright © 2016 Shufen Feng 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.

To investigate the association between the metalloproteinase-9 (MMP9) −1562C/T polymorphism and diabetic nephropathy (DN) in Han Chinese, the patients with type 2 diabetes were collected and divided into the non-DN (NDN) and DN groups; controls were recruited. Genotype and allele frequencies were assessed using polymerase chain reaction and restriction fragment length polymorphism. Results showed that SBP, DBP, HbA1c, UAER, Cr, BUN, TG, and TC were higher in the DN group compared with the control and NDN groups. SBP, HbA1c, and TC in DN patients with the TT and CT genotypes were lower than in those with CC. Compared with controls, the frequency of the T allele in the DN group was significantly lower. The MMP9 −1562C allele, SBP, Cr, BUN, TG, and TC were independent risk factors for DN. All of the above suggested that the MMP9 −1562C/T polymorphism was associated with DN in Han Chinese.

1. Introduction

Type 2 diabetes mellitus (T2DM) is a metabolic disorder characterized by hyperglycemia resulting from insulin resistance [1]. T2DM is a severe condition and is often associated with life-threatening conditions such as heart disease and stroke. The worldwide prevalence of T2DM is about 387 million. The incidence of T2DM is increasing worldwide, particularly in Asian countries. The prevalence of T2DM in China in 2014 was about 96.3 million [2]. Diabetic nephropathy (DN) occurs in almost 40% of patients with T2DM [3] and is one of the leading causes of T2DM-related death. DN is characterized by progressive kidney damage resulting from T2DM, but the risk factors for DN vary among populations. Therefore, a better understanding of the pathogenesis of DN is critical for early disease intervention.

DN development is related to inflammation and matrix metalloproteinases (MMP) play key roles in the progression of DN [4]. Several studies indicated that MMP9 is an important inflammatory marker involved in the pathophysiology of DN [5]. Proteins belonging to the MMP family are involved in the breakdown of extracellular matrix (ECM) in normal physiological processes such as tissue remodeling, but degradation of the ECM by MMP is also involved in the pathogenesis of a number of diseases [6]. Major pathological changes are associated with DN and result from the accumulation of ECM in the glomerular mesangium and the renal interstitium, leading to progressive glomerulosclerosis and renal interstitial fibrosis [7]. MMP9 plays a critical role in diabetic complications such as diabetic retinopathy and DN [8]. In addition, MMP9 is detectable in a number of other renal conditions such as immunoglobulin A nephropathy, membranoproliferative glomerulonephritis, and lupus nephritis [9]. Bernerine has renal protective effects against DN through the regulation of MMPs [10]. Therefore, taken together, there might be an association between MMP9 and DN, but whether MMP9 promotes or inhibits the development of DN is still controversial [11].

Several mechanisms may be responsible for regulating MMP9 expression and activity in DN, and single nucleotide polymorphisms (SNPs) within the MMP9 gene might influence its expression. The effects of a number of SNPs have been described for a number of MMPs [12]. A recent paper reported SNP in the MMP9 gene promoter that reduces the risk of diabetic microvascular complications [13], but there is no study specifically examining MMP9 SNPs in DN. A previous study showed that the MMP9 −1562C/T SNP was associated with diabetic macroangiopathy [14], but there is no study specifically examining MMP9 SNPs in DN.

Therefore, the present study examined SNPs in the MMP9 gene among Han Chinese individuals with DN, without DN (NDN), and without kidney disease (controls).

2. Materials and Methods

2.1. Subjects

This case-control study included 310 patients diagnosed with T2DM at the Endocrinology Department of The First Affiliated Hospital of Jinan University (Guangzhou, China) between January 2013 and August 2013. Diagnosis of T2DM was in line with the diagnostic and classification criteria for diabetes established by the World Health Organization (WHO) in 1999 [15]. Patients with T2DM were categorized into the DN group (; 78 males and 90 females) and the NDN group (; 80 males and 62 females). Diagnostic criteria for DN included urinary albumin excretion rate (UAER) >30 mg/24 h, serum creatinine (Cr) >167 μmol/L, and blood urea nitrogen (BUN) >7.8 mmol/L, confirmed in at least two independent consecutive examinations [16]. Patients with other primary/secondary renal diseases or those who had used nephrotoxic drugs were excluded.

During the same period, 100 healthy volunteers were recruited as controls (55 males and 45 females). Individuals with cardiovascular diseases, a history of diabetes, or a family history of liver, kidney, endocrine, or metabolic disease were excluded.

All subjects were unrelated to one another. The present study was approved by the Ethics Committee of The First Affiliated Hospital of Jinan University. All subjects agreed to have blood drawn and provided signed informed consent.

2.2. Clinical Data

Whole blood (5 mL) was sampled from each fasting participant in EDTA-containing tubes; 2 mL was stored at −20°C for DNA extraction and 3 mL was used for plasma separation and detection of fasting glucose (FPG), Cr, BUN, triglycerides (TG), total cholesterol (TC), and other biochemical markers using an automatic biochemical analyzer. Glycated hemoglobin (HbA1c) was detected using an HbA1c analyzer. UAER was detected using a radioimmunoassay. Systolic blood pressure (SBP) and diastolic blood pressure (DBP) were measured after a 5 min rest in the sitting position.

2.3. Polymorphism Analysis

Genomic DNA was extracted from whole blood using a blood genomic DNA extraction kit (Shanghai Generay Biological Engineering Co., Ltd., Shanghai, China) and stored at −20°C. PCR primers were designed and synthesized by Shanghai Sangon Biological Engineering Technology Services Ltd. (Shanghai, China). The primers for the MMP9 promoter region were as follows: sense 5′-CTT CCT AGC CAG CCG GCA TC-3′ and antisense 5′-GCC TGG CCT ATA GTA GGC CC-3′. The target sequence was amplified by PCR, generating an amplification product of 435 bp. The PCR reaction volume was 25 μL and included 0.5 μL of each primer, 1.5 μL of template DNA, 12.5 μL of TaqDNA polymerase (1.25 U), and sterile double-distilled water. The reaction conditions were as follows: (1) 94°C for 5 min; (2) 35 cycles of amplification at 94°C for 35 s, 60°C for 30 s, and 72°C for 45 s; and (3) 72°C for 5 min. PCR was performed using an Applied Biosystems PCR System. PCR amplification product (3 μL) was digested with 2 U of the SphI restriction endonuclease at 37°C for 16 h to detect the MMP9 −1562C/T SNP. After digestion, the products were separated on 3% ethidium bromide-stained agarose gels (Biowest, Madrid, Spain). Images were acquired and analyzed using a gel imaging system.

2.4. Statistical Analysis

Data analysis was performed using SPSS 18.0 (IBM, Armonk, NY, USA). The Hardy-Weinberg equilibrium was used to verify that all samples were representative of the general population. Continuous data were expressed as mean ± standard deviation. Categorical data were expressed as frequencies. Intergroup comparisons were performed using analysis of variance (ANOVA) and the Bonferroni post hoc test. Genotype and allele frequencies were compared among groups with the chi-square test. Logistic regression analysis was conducted to explore the relationship between genotypes and relative risk of DN and was expressed as odds ratio (OR) and 95% confidence intervals (CI). Two-sided values < 0.05 were considered statistically significant.

A post hoc power analysis was performed using the QUANTO software (http://biostats.usc.edu/software). Considering that the prevalence of DN among diabetic Chinese subjects is about 20% () and considering that the frequency of the null/null genotype among study subjects is 53%, . The hazard ratio of the null/null genotype patients compared with healthy subjects is therefore 2.27 (). Consequently, using a power of 0.8 and two-sided α of 0.05, the minimal sample size was 99.

3. Results

3.1. Characteristics of the Subjects

There was no statistical difference in gender, age, height, or BMI (all ) among the three groups. SBP, DBP, HbA1c, UAER, Cr, BUN, TG, and TC were all significantly higher in the DN group compared with the control group (all ). SBP, DBP, HbA1c, Cr, BUN, TG, and TC were also significantly higher in the NDN group compared with the control group (all ). Finally, SBP, HbA1c, UAER, Cr, BUN, TG, and TC were all significantly higher in the DN group compared with the NDN group (all ) (Table 1).

Table 1: Characteristics of the subjects.

NC
NDN
DN

Cases
100
142
168

Gender (M/F)
55/45
80/62
78/90

Age (years)
57.07 10.82
56.48 9.10
54.67 11.81

Height (cm)
164.52 6.53
165.37 7.12
166.17 7.08

BMI (kg/m²)
21.47 1.08
24.58 2.07
25.09 1.92

SBP (mmHg)
121.68 7.13
138.41 8.43
149.51 8.65

DBP (mmHg)
84.78 9.41
106.83 10.26
111.42 9.30

HbA1c (%)
4.98 0.85
6.51 0.81
8.01 0.96

UAER (mg/24 h)
17.14 2.51
20.78 1.59
301.93 51.32

Cr (mol/L)
60.51 7.21
97.92 15.23
129.71 23.53

BUN (mmol/L)
4.04 0.51
5.47 0.67
9.43 0.82

TG (mmol/L)
1.06 0.08
2.13 0.24
2.25 0.32

TC (mmol/L)
4.12 0.26
4.58 0.49
4.96 0.54

versus the NC group.

versus the NDN group.

BMI: body mass index; SBP: systolic blood pressure; DBP: diastolic blood pressure; UAER: urinary albumin excretion rate; TC: total cholesterol; TG: triglycerides; Cr: creatinine; BUN: blood urea nitrogen.

3.2. MMP9 −1562C/T Allele SNPs

The CC genotype was represented by a single 435-bp band. The TT genotype was identified by two bands (247 and 188 bp). Finally, the CT genotype was represented by three bands (435, 247, and 188 bp (Figure 1)).

Figure 1: Electrophoresis of amplification fragments of the MMP9 promoter region for the detection of the −1562C/T polymorphism. Lanes 1, 2, and 5 are the CC genotype; lanes 3, 4, and 7 are the TT genotype; lanes 6 and 8 are the CT genotype. M: DNA molecular weight marker.

3.3. Genotype Distribution and Allele Frequencies of the MMP9 −1562C/T Allele

The Hardy-Weinberg equilibrium test showed that the distribution of the MMP9 −1562C/T allele in the NC, NDN, and DN groups was in line with the law of genetic equilibrium () (Table 2). There were significant differences in allele and genotype frequencies between the DN and NC groups (). The difference in genotype between the DN and NC groups had a recessive inheritance pattern (CC versus CT + TT) and was statistically significant (). The frequency of the T allele in the DN group was significantly lower than in the NC group (). The relative risk of suffering from DN in subjects carrying the T allele was 0.478 compared with individuals with the C allele (OR = 0.478, 95% CI: 0.288–0.791). There was no difference in allele and genotype frequencies between the NC and NDN groups (). The genotype frequencies of the DN group were different compared with the NDN group (), but there was no significant difference in allele frequency () (Table 3).

Table 2: Genetic equilibrium test for MMP9 gene −1562C/T mutation.

NC
NDN
DN

Cases

100
142
168

Genotype frequency
CC
45/48.05
64/65.86
101/98.52

CT
13/11.25
30/31.26
37/35.01

TT
42/40.70
48/44.74
30/34.47

0.79
0.86
0.62

Table 3: Distribution frequency of MMP9 −1562C/T alleles [case (%)].

NC
NDN
DN

100
142
168

Genotype frequency
CC
45 (45.0)
64 (45.1)
101 (60.1)

CT
13 (13.0)
30 (21.1)
37 (22.0)

TT
42 (42.0)
48 (33.8)
30 (17.9)

Allele frequency
C alleles
53 (53.0)
75 (52.8)
119 (52.8)

T alleles
47 (47.0)
67 (47.2)
49 (29.2)

versus the NC group.

versus the NDN group.

3.4. Comparison of Blood Glucose, Blood Pressure, and Blood Lipids in Patients with DN

SBP, HbA1c, and TC levels in patients with DN and the MMP9 −1562 TT and CT genotypes were significantly lower than in those with the CC genotype (all ) (Table 4).

Table 4: Comparison of blood pressure, blood glucose, and blood lipid levels between the DN patients carrying MMP9 −1562T and MMP9 −1562C genotypes.

CT + TT
CC

Cases
67
101

SBP (mmHg)
138.41 8.43
149.51 8.65
0.032

DBP (mmHg)
106.83 10.26
111.42 9.30
0.612

HbA1c (%)
6.51 0.81
8.01 0.96
0.004

TC (mmol/L)
4.79 1.08
5.48 1.25
0.017

TG (mmol/L)
1.85 0.84
1.86 0.91
0.057

SBP: systolic blood pressure; DBP: diastolic blood pressure; TC: total cholesterol; TG: triglycerides.

3.5. Multivariate Analysis

To predict the probability of DN in the NDN and DN groups, MMP9 −1562C genotypes (dummy variables were created and the TT genotype was used as control), DBP, SBP, HbA1c, Cr, BUN, TG, TC, and UAER were entered as independent variables, and the binary logistic regression analysis was performed using the forward selection procedure. MMP9 −1562C genotypes (CT: OR = 1.531, 95% CI: 1.186–1.975, ; CC: OR = 7.402, 95% CI: 4.534–12.083, ), SBP (OR = 1.6, 95% CI: 1.243–2.081, ), Cr (OR = 6.798, 95% CI: 4.357–12.479, ), BUN (OR = 4.423, 95% CI: 2.765–9.831, ), TG (OR = 6.466, 95% CI: 3.047–11.926, ), and TC (OR = 3.522, 95% CI: 2.105–6.269, ) were independent risk factors for DN (Table 5).

Table 5: Multivariate logistic regression analysis.

Variables

OR
95% CI

MMP9 −1562C
TT
Reference
—
—

CT
1.531
1.186–1.975
0.001

CC
7.402
4.534–12.083
<0.001

SBP

1.600
1.243–2.081
<0.001

TC

3.522
2.105–6.269
0.001

TG

6.466
3.047–11.926
0.039

Cr

6.798
4.357–12.479
0.042

BUN

4.423
2.765–9.831
0.006

OR: odds ratio; 95% CI: 95% confidence interval; SBP: systolic blood pressure; TC: total cholesterol; TG: triglycerides; Cr: creatinine; BUN: blood urea nitrogen.

4. Discussion

In this study, the relationship between MMP9 SNPs and DN was examined in the Han Chinese population in Guangzhou (China). SBP, DBP, HbA1c, UAER, Cr, BUN, TG, and TC were higher in the DN group compared with the control and NDN groups. SBP, HbA1c, and TC in DN patients with the TT and CT genotypes were lower than in those with the CC genotype. The relative risk of DN in individuals carrying the T allele was 0.478 compared with individuals carrying the C allele. The MMP9 −1562C genotype, SBP, Cr, BUN, TG, and TC were independent risk factors for DN.

The MMP9 gene is located on the chromosome 20q11.2–13.1 in humans and contains four promoter regions. MMP9 −1562C/T is located upstream of the promoter and contains binding sites for transcriptional repressors [17, 18]. In the presence of the T allele, transcriptional repression is decreased or eliminated, probably leading to elevated expression of MMP9, increased degradation of ECM, and delayed interstitial fibrosis and renal glomerular sclerosis; as a result, DN development and progression are delayed. Therefore, the C allele probably alters MMP9 expression to promote the occurrence and development of DN. A previous study has shown that the T allele of the MMP9 −1562C/T SNP was associated with diabetic macroangiopathy [14]. Two studies showed that the T allele was associated with diabetic foot ulcers and nonhealing wounds in diabetic subjects [19, 20]. This discrepancy could be related to a number of factors, including the subjects themselves, the choice of the controls, and the nature and pathogenesis of the lesions (macroangiopathy versus DN). Additional studies are necessary to solve this issue.

Compared with the control and NDN groups, SBP, HbA1c, UAER, Cr, BUN, TG, and TC were all elevated in the DN group, suggesting that blood pressure, blood glucose level, and kidney function are important factors involved in DN. SBP, HbA1c, and TC levels in patients with DN and the TT/CT genotypes were significantly lower than in patients with the CC genotype. In addition, the recessive inheritance pattern (i.e., CC versus CT + TT) between the DN and control groups was statistically different. Compared with the control group, the T allele frequency in the DN group was significantly lower, and the relative risk of DN in individuals carrying the T allele was only 0.478 compared with the individuals carrying the C allele. In the present study, the MMP9 −1562C genotype, SBP, Cr, BUN, TG, and TC were all independent risk factors for the occurrence of DN. These findings strongly suggest that the T allele is probably a protective factor against the development of DN, while the C allele probably promotes DN development.

A number of groups have previously performed genome-wide association studies to identify genetic variants associated with the risk of developing T2DM [21–23]. To date, there are approximately 70 loci that are significantly associated with T2DM, but it is important to note that many of these SNPs are in noncoding regions. Thus, the exact mechanisms by which they influence the risk of disease are not very well understood. Along these same lines, many studies have identified genomic regions associated with susceptibility to DN, but the identification of a single causal variant has not been reported. The present study is unique in that the −1562C/T variant occurs within the MMP9 promoter and within a region that is known to influence the expression of MMP9. This provides a mechanistic basis explaining why the C allele contributes to DN and why the T allele is protective. It would be interesting to explore if the −1562C/T variant identified in this study adds to the predictive power of any of the previously identified T2DM susceptibility alleles, but additional comprehensive and in-depth studies are necessary to address this issue.

This study is not without limitations. The sample size was relatively small and from a single geographic region. Indeed, there are differences in genotypes and allele frequencies between different ethnic groups and/or regions [24]. Additional studies are necessary to examine the contribution of the MMP9 SNPs to DN.

5. Conclusions

The MMP9 −1562C/T SNP was associated with DN in a Han Chinese population. The T allele was a protective factor for DN, while the C allele contributed to the disease.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Shufen Feng and Gang Ye contributed equally and share the first authorship.

Acknowledgments

The authors thank Professor Junning Chen of the Nephrology Department of Affiliated Hospital of Guilin Medical College and Lu Luo of the Endocrinology Department of The First Affiliated Hospital of Jinan University for providing research ideas and clinical data; Fonan Li of the Hemodialysis Center of The First Affiliated Hospital of Jinan University for helping with the blood collection; and Professor Chunyan Feng of Clinical Laboratory of People’s Hospital of Shenzhen Longhua District for allowing them access to instruments and equipment. The authors also thank Dr. Ji Fang for critically reading and reviewing the paper. In particular, they sincerely thank Professor Jianquan Yang of the Gastroenterology Department of The First Affiliated Hospital of Jinan University and Professor Ronghua Zhang of Pharmacy College of Jinan University for providing financial support. This work was supported by the Science and Technology Planning Project of Guangdong Province of China (no. 2013B022000075) and the Administration of Traditional Chinese Medicine of Guangdong Province of China (no. 20131152).

References

R. Taylor, “Type 2 diabetes: etiology and reversibility,” Diabetes Care, vol. 36, no. 4, pp. 1047–1055, 2013.

International Diabetes Federation, IDF Diabetes Atlas, International Diabetes Federation, Brussels, Belgium, 2014, http://www.idf.org/diabetesatlas.

S. Rizvi, S. T. Raza, and F. Mahdi, “Association of genetic variants with diabetic nephropathy,” World Journal of Diabetes, vol. 5, no. 6, pp. 809–816, 2014.

X. Xu, L. Xiao, P. Xiao et al., “A glimpse of matrix metalloproteinases in diabetic nephropathy,” Current Medicinal Chemistry, vol. 21, no. 28, pp. 3244–3260, 2014.

N. J. van der Zijl, R. Hanemaaijer, M. E. Tushuizen et al., “Urinary matrix metalloproteinase-8 and -9 activities in type 2 diabetic subjects: a marker of incipient diabetic nephropathy?” Clinical Biochemistry, vol. 43, no. 7-8, pp. 635–639, 2010.

C. J. Malemud, “Matrix metalloproteinases (MMPs) in health and disease: an overview,” Frontiers in Bioscience, vol. 11, no. 2, pp. 1696–1701, 2006.

S. O. Kolset, F. P. Reinholt, and T. Jenssen, “Diabetic nephropathy and extracellular matrix,” Journal of Histochemistry and Cytochemistry, vol. 60, no. 12, pp. 976–986, 2012.

R. A. Kowluru, Q. Zhong, and J. M. Santos, “Matrix metalloproteinases in diabetic retinopathy: potential role of MMP-9,” Expert Opinion on Investigational Drugs, vol. 21, no. 6, pp. 797–805, 2012.

M. Sekiuchi, A. Kudo, K. Nakabayashi et al., “Expression of matrix metalloproteinases 2 and 9 and tissue inhibitors of matrix metalloproteinases 2 and 1 in the glomeruli of human glomerular diseases: the results of studies using immunofluorescence, in situ hybridization, and immunoelectron microscopy,” Clinical and Experimental Nephrology, vol. 16, no. 6, pp. 863–874, 2012.

W.-J. Ni, H.-H. Ding, H. Zhou, Y.-Y. Qiu, and L.-Q. Tang, “Renoprotective effects of berberine through regulation of the MMPs/TIMPs system in streptozocin-induced diabetic nephropathy in rats,” European Journal of Pharmacology, vol. 764, pp. 448–456, 2015.

K. M. Thrailkill, R. Clay Bunn, and J. L. Fowlkes, “Matrix metalloproteinases: their potential role in the pathogenesis of diabetic nephropathy,” Endocrine, vol. 35, no. 1, pp. 1–10, 2009.

G. Martin, V. Asensi, A. H. Montes et al., “Role of plasma matrix-metalloproteases (MMPs) and their polymorphisms (SNPs) in sepsis development and outcome in ICU patients,” Scientific Reports, vol. 4, article 5002, 2014.

Z. Zhang, X. Wu, T. Cai et al., “Matrix metalloproteinase 9 gene promoter (rs 3918242) mutation reduces the risk of diabetic microvascular complications,” International Journal of Environmental Research and Public Health, vol. 12, no. 7, pp. 8023–8033, 2015.

Y. Wang, Y. Su, Y. Xu, S.-H. Pan, and G.-D. Liu, “Genetic polymorphism c.1562C>T of the MMP-9 is associated with macroangiopathy in type 2 diabetes mellitus,” Biochemical and Biophysical Research Communications, vol. 391, no. 1, pp. 113–117, 2010.

K. G. Alberti and P. Z. Zimmet, “Definition, diagnosis and classification of diabetes mellitus and its complications—part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation,” Diabetic Medicine, vol. 15, no. 7, pp. 539–553, 1998.

J. L. Gross, M. J. de Azevedo, S. P. Silveiro, L. H. Canani, M. L. Caramori, and T. Zelmanovitz, “Diabetic nephropathy: diagnosis, prevention, and treatment,” Diabetes Care, vol. 28, no. 1, pp. 164–176, 2005.

C. Yan and D. D. Boyd, “Regulation of matrix metalloproteinase expression expression,” Journal of Cellular Physiology, vol. 211, no. 1, pp. 19–26, 2007.

I. M. Clark, T. E. Swingler, C. L. Sampieri, and D. R. Edwards, “The regulation of matrix metalloproteinases and their inhibitors,” International Journal of Biochemistry and Cell Biology, vol. 40, no. 6-7, pp. 1362–1378, 2008.

K. Singh, N. K. Agrawal, S. K. Gupta, G. Mohan, S. Chaturvedi, and K. Singh, “Differential expression of matrix metalloproteinase-9 gene in wounds of type 2 diabetes mellitus cases with susceptible -1562c>T genotypes and wound severity,” The International Journal of Lower Extremity Wounds, vol. 13, no. 2, pp. 94–102, 2014.

K. Singh, N. K. Agrawal, S. K. Gupta, and K. Singh, “A functional single nucleotide polymorphism −1562C>T in the matrix metalloproteinase-9 promoter is associated with type 2 diabetes and diabetic foot ulcers,” International Journal of Lower Extremity Wounds, vol. 12, no. 3, pp. 199–204, 2013.

E. Zeggini, L. J. Scott, R. Saxena, and B. F. Voight, “Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes,” Nature Genetics, vol. 40, no. 5, pp. 638–645, 2008.

B. F. Voight, L. J. Scott, V. Steinthorsdottir et al., “Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis,” Nature Genetics, vol. 42, no. 7, pp. 579–589, 2010.

L. K. Billings and J. C. Florez, “The genetics of type 2 diabetes: what have we learned from GWAS?” Annals of the New York Academy of Sciences, vol. 1212, pp. 59–77, 2010.

L. Chen, T. Wang, L. Liu, Y. Shen, C. Wan, and F. Wen, “Matrix metalloproteinase-9 -1562C/T promoter polymorphism confers risk for COPD: a meta-analysis,” PLoS ONE, vol. 8, no. 3, Article ID e60523, 2013.