p22phox C242T gene polymorphism and overt diabetic nephropathy: a meta-analysis of 1,452 participants

Article information

Korean J Intern Med. 2017;32(6):1045-1052

Publication date (electronic) : 2016 December 8

doi : https://doi.org/10.3904/kjim.2015.106

Yan-yan Li, Ge Gong, Hong-yu Geng, Yun Qian

Department of Geriatrics, First Affiliated Hospital of Nanjing Medical University, Nanjing, China

Correspondence to Yan-yan Li, Ph.D Department of Geriatrics, First Affiliated Hospital of Nanjing Medical University, No. 300, Guangzhou Rd, Nanjing 210029, China Tel: +86-25-83718836 Fax: +86-25-83780170 E-mail: lyynjmu123@126.com

Received 2015 April 18; Revised 2015 July 3; Accepted 2015 July 20.

Abstract

Background/Aims

The p22phox C242T gene polymorphism (rs4673) may be linked to an increased susceptibility for overt diabetic nephropathy (ODN), but the study results are still inconclusive.

Methods

To explore the relationship between p22phox C242T gene polymorphism and ODN, the current meta-analysis of 707 ODN patients and 745 controls from five individual studies was conducted. The pooled odds ratio (OR) and its corresponding 95% confidence interval (CI) were evaluated by either a random or fixed effect model.

Results

In our meta-analysis, a significant relationship between the p22phox C242T gene polymorphism and ODN was found under allelic (OR, 2.760; 95% CI, 1.400 to 5.450; p = 0.004), recessive (OR, 5.080; 95% CI, 1.020 to 25.430; p = 0.05), dominant (OR, 1.700; 95% CI, 1.167 to 2.477; p = 0.006), homozygous (OR, 3.900; 95% CI, 1.022 to 14.889; p = 0.046), heterozygous (OR, 1.523; 95% CI, 1.167 to 1.986; p = 0.002), and additive genetic models (OR, 2.019; 95% CI, 1.232 to 3.309; p = 0.005).

Conclusions

A positive correlation between p22phox C242T gene polymorphism and ODN risk was found. The T allele carriers of p22phox C242T gene polymorphism might be predisposed to ODN.

Keywords: p22-phox; C242T; Polymorphism; Overt diabetic nephropathy; Diabetes

INTRODUCTION

Diabetic nephropathy (DN) is a severe complication of diabetes mellitus and is characterized by persistent proteinuria, decreased glomerular filtration rate, and hypertension with high morbidity and mortality due to cardiovascular diseases. Although DN can occur in both types of diabetes mellitus, studies indicate that the incidence of DN is higher in patients with type 1 diabetes mellitus (T1DM) than in those with type 2 diabetes mellitus (T2DM) (T1DM, 30% to 40%; T2DM, 20%). In the United States and other developed countries, overt diabetic nephropathy (ODN) has become the primary cause of end-stage renal failure, and 60% of patients undergoing hemodialysis in the United States suffer from ODN [1-3].

The underlying mechanism of ODN is multifactorial, resulting from the interactions among abnormal metabolism, hemodynamic changes, and various growth and genetic factors [4]. The principal pathological characteristic of glomerulosclerosis is the thickening of the glomerular basement membrane and the accumulation of the extracellular matrix [5]. Recently, Najafi et al. [6] have discovered that nicotinamide adenine dinucleotide phosphate (NADPH) oxidase activity is closely associated with p22phox (CYBA) C242T polymorphism (rs4673) because p22phox is an essential subunit for NADPH.

The p22phox gene, located in 16q24, spans 8.5 kb and contains six exons and five introns. In the C242T polymorphism (rs4673), the 242nd base cytosine (C) is substituted by thymine (T), resulting in the wild-type histidine (His) being replaced by tyrosine (Tyr) at the 72nd codon of the fourth exon. Given that the 72nd codon is the heme-binding site, the His > Tyr variant reduces the heme-binding affinity. Consequently, NADPH oxidase activity is affected, and the production of reactive oxygen species is significantly altered, ultimately contributing to renal cellular or organ injury [7].

Although many studies on the relationship between the p22phox C242T gene polymorphism and DN have been conducted, individual results are slightly contradictory. In 2003, Hodgkinson et al. [8] first reported that within the population of the United Kingdom, the frequency of the 242 TT genotype of the p22phox C242T gene polymorphism significantly increases in T1DM patients with ODN compared with those with retinopathy alone or uncomplicated diabetes (no microvascular disease after 20 years of diabetes duration), as well as normal healthy controls (33.3% vs. 6.5%, 5.7%, and 0.0%, respectively; p < 0.000001). Similarly, Liu et al. [9] found that the p22phox C242T gene polymorphism positively correlates with ODN and concluded that the T allele likely predisposes the Chinese population with T2DM to ODN. Yang et al. [10] drew the same conclusion in their study on the T2DM population in China. By contrast, Lim et al. [11] reported that the p22phox C242T gene polymorphism does not significantly increase the susceptibility to DN in another T2DM Chinese population.

The present meta-analysis involving 707 ODN cases and 745 controls was performed to explore the relationship between the p22phox C242T gene polymorphism and ODN (Supplementary Table 1) [12]. To avoid classifying subjects without significant renal impairment as cases and then erroneously enriching the cases with pathogenetic genotype, we only considered patients with ODN as cases in this meta-analysis.

METHODS

Publication search and inclusion criteria

Databases, including Web of Science, PubMed, Embase, China National Knowledge Infrastructure, and China Biological Medicine Database, were searched using the keywords ‘overt diabetic nephropathy,’ ‘p22phox,’ ‘C242T,’ and ‘polymorphism’ for our primary search of the literature. The retrieved studies were searched from 1980 (last research updated on July 1, 2015).

The inclusion criteria for our meta-analysis were as follows. (1) The studies must have evaluated the relationship between the p22phox C242T gene polymorphism and ODN. (2) ODN must be diagnosed by the presence of macroalbuminuria (albumin excretion rate ≥ 200 µg/min) or chronic renal failure treated by dialysis. Other renal diseases or other causes of proteinuria were excluded by renal biopsy. Furthermore, proteinuria must be measured more than once. In addition, the ODN patients must present a long duration of diabetes and be accompanied by diabetic retinopathy and other complications, or the ODN was confirmed by renal biopsy. (3) The studies must be published in official journals as case-control or cohort studies. (4) The Hardy-Weinberg equilibrium (HWE) must be followed by the studies.

Data extraction

Studies that failed to meet inclusion criteria or those with insufficient data were excluded. Data published multiple times was used in our analysis only once. The extracted data had to include the first author’s name, publication year, study region, number of genotypes, genotyping method, study design, matching criteria, and total number of ODN cases and controls.

Statistical analysis

Six genetic models, the allelic (T allele distribution frequency), recessive (TT vs. CC + CT), dominant (TT + CT vs. CC), homozygous (TT vs. CC), heterozygous (CT vs. CC), and additive genetic models (T vs. C), were adopted for our meta-analysis. The odds ratio (OR) corresponding to 95% confidence interval (CI) was used to compare the relationship between the p22phox C242T gene polymorphism and ODN. The heterogeneity between studies was calculated by chi-square-based q test with the threshold of significance set at p < 0.05 [13]. If heterogeneity were among the individual studies, the random-effects model (the DerSimonian and Laird method) would be used to estimate the pooled OR [14]. If not, a fixed-effects model (the Mantel-Haenszel method) would be adopted [15]. The Z test was used to determine the pooled OR with the threshold of significance was set at p ≤ 0.05 level.

The Fisher exact test was used to assess the HWE and significance was set at p < 0.05 level. The potential publication bias was estimated by funnel plot. The funnel plot asymmetry was assessed by Egger’s linear regression test on the natural logarithm scale of the OR with significance set as p < 0.05 [16]. The statistical analysis was performed using STATA version 11.0 software (Stata Corp., College Station, TX, USA).

RESULTS

Studies and populations

Our initial search of the literature retrieved 15 papers. Five of these papers met our inclusion criteria. Only one of these papers analyzed the relationship between the p22phox C242T gene polymorphism and ODN within the context of T1DM; the remaining four enrolled patients with T2DM. Among the 10 rejected papers, two papers were repeated publications, four papers were reviews, and four papers were unrelated to the subject at hand. In total, our data comprised 707 ODN patients and 745 controls (Table 1, Supplementary Fig. 1) [8-11,12,17]. Study regions comprised the United Kingdom, Brazil, and China. The populations were categorized into two ethnicities, namely, Caucasian and Chinese.

Table 1.

Characteristics of the investigated studies of the association of the p22phox C242T gene polymorphism and ODN

Pooled analyses

A significant relationship existed between the p22phox C242T gene polymorphism and ODN under allelic (OR, 2.760; 95% CI, 1.400 to 5.450; p = 0.004), recessive (OR, 5.080; 95% CI, 1.020 to 25.430; p = 0.05), dominant (OR, 1.700; 95% CI, 1.167 to 2.477; p = 0.006), homozygous (OR, 3.900; 95% CI, 1.022 to 14.889; p = 0.046), heterozygous (OR, 1.523; 95% CI, 1.167 to 1.986; p = 0.002), and additive genetic models (OR, 2.019; 95% CI, 1.232 to 3.309; p = 0.005) (Table 2, Figs. 1 and 2).

Table 2.

Summary of meta-analysis of association of p22phox C242T gene polymorphism and ODN

Figure 1.

Forest plot of overt diabetic nephropathy associated with p22phox C242T gene polymorphism under allelic genetic model (distribution of T allelic frequency of p22phox gene). DN, diabetic nephropathy; OR, odds ratio; CI, confidence interval.

Figure 2.

Forest plot of overt diabetic nephropathy (ODN) associated with p22phox C242T gene polymorphism under recessive genetic model (TT vs. CT + CC). DN, diabetic nephropathy; OR, odds ratio; CI, confidence interval.

After excluding the study on ODN from T1DM [8], the four remaining studies were adopted to perform the meta-analysis, and similar results were obtained. A significant relationship between the p22phox C242T gene polymorphism and ODN was detected under the allelic (OR, 2.320; 95% CI, 1.160 to 4.650; p = 0.02), dominant (OR, 2.317; 95% CI, 1.186 to 4.527; p = 0.014), homozygous (OR, 2.062; 95% CI, 1.128 to 3.767; p = 0.019), heterozygous (OR, 1.612; 95% CI, 1.210 to 2.165; p = 0.001), and additive genetic models (OR, 2.343; 95% CI, 1.152 to 4.764; p = 0.019). No significant relationship existed between the p22phox C242T gene polymorphism and ODN under the recessive genetic model (OR, 1.675; 95% CI, 0.952 to 2.947; p = 0.074) (Table 3, Figs. 3 and 4).

Table 3.

Summary of meta-analysis of association of p22phox C242T gene polymorphism and ODN confined to type 2 diabetes mellitus

Figure 3.

Forest plot of p22phox C242T gene polymorphism associated with overt diabetic nephropathy (ODN) confined to type 2 diabetes mellitus under allelic genetic model (distribution of T allelic frequency of p22phox gene). IV, inverse variance method; CI, confidence interval; DN, diabetic nephropathy.

Figure 4.

Forest plot of p22phox C242T gene polymorphism associated with overt diabetic nephropathy confined to type 2 diabetes mellitus under dominant genetic model (distribution of T allelic frequency of p22phox gene). OR, odds ratio; CI, confidence interval.

Bias diagnostics

The potential for publication bias of the individual studies was assessed using the funnel plot and Egger’s test. No visual publication bias was detected in the funnel plot (Fig. 5). Moreover, no significant difference was found in the Egger’s test, and this finding implied that no publication bias existed in the present meta-analysis, as shown by the allelic genetic model (T = −1.19, p = 0.321).

Figure 5.

Funnel plot for studies of the association of overt diabetic nephropathy associated and p22phox C242T gene polymorphism under an allelic genetic model (distribution of T allelic frequency of p22phox gene). The horizontal and vertical axis correspond to the odds ratio (OR) and confidence limits. SE, standard error.

DISCUSSION

In this meta-analysis, a significant positive correlation was found between the p22phox C242T gene polymorphism and ODN under all six genetic models, namely, allelic (OR, 2.760), recessive (OR, 5.080), dominant (OR, 1.700), homozygous (OR, 3.900), heterozygous (OR, 1.523), and additive genetic models (OR, 2.019). Even after the study on ODN from T1DM was excluded, a similar conclusion was drawn. Thus, we believe that our results are significant. Additionally, the T allele carriers of the p22phox C242T gene polymorphism might predispose to ODN.

NADPH oxidase is a specific electron-transfer compound composed of five elements: p40phox, p47phox, p67phox, and cytochrome b588 (p22phox, gp91phox) [18]. Cytochrome b588 is the main element of NADPH and includes two subunits as heavy and light chains. The heavy chain is a 91 kDa glycoprotein, and the light chain is a 22 kDa polypeptide. The former is the core section of NADPH oxidase, and the latter is the key position of NADPH-producing oxygen radicals [19]. NADPH oxidase is a key enzyme that produces superoxides. In particular, p22phox plays a crucial role in NADPH normal function and O²⁻ generation.

Whether the p22phox C242T gene variant contributes to vascular damage by increasing or decreasing NADPH oxidase activity was previously unclear. In 2000, Guzik et al. [20] found that the p22phox 242T allele is associated with reduced NADPH oxidase activity in human blood vessels and concluded that genetic variation in NADPH oxidase components might play a significant role in modulating superoxide production in human atherosclerosis. By contrast, Perianayagam et al. [7] found that patients with acute renal failure, who also exhibited the p22phox 242TT polymorphism, present higher levels of plasma nitrotyrosine than those with the wild-type genotype. This finding indicates an increase in NADPH activity induced by the T allele [7]. Enzyme activity is modulated by many factors in vivo, and the different experimental approaches of these two studies may contribute to the varying results [21]. NADPH activity is likely increased by the genetic polymorphism, considering that the latter experiment was performed in vivo.

The differential effects of the C242T variants on other biological processes, such as the upregulation of antioxidant defenses, inflammation, lipid peroxidation, gene expression, and apoptosis, can serve as potential explanations for their involvement in vascular damage [22,23]. In 2003, Hodgkinson et al. [8] found that NADPH oxidase with the polyol pathway may contribute to the pathogenesis of DN.

In the current meta-analysis, we selected the cases of ODN to avoid enrolling incipient renal impairment, but this selection can also enrich the group of cases. A significant heterogeneity among the individual studies might be associated with the gene variant distribution difference in diverse ethnicities. Additionally, DN may behave differently in T1DM and T2DM, but the present meta-analysis reached the same conclusion even after the study regarding ODN and T1DM was excluded. In 2003, Ji et al. [24] found that the p22phox C242T allele and genotype distribution in the Chinese population are similar to those in the Japanese population and much lower than those in Caucasians and Indians [23-26]. Given that only five manuscripts were included, meta-regression was not performed to explore the heterogeneity source in the current meta-analysis.

Despite the robust results of this analysis, which is the first to explore the association between the p22phox C242T gene polymorphism and ODN, some limitations still affect this meta-analysis. First, few individual studies were included in this meta-analysis, and large-scale studies on the relationship of ODN to the p22phox C242T gene polymorphism are lacking. Second, NADPH oxidase activity can also be influenced by other genetic polymorphisms, such as A640G (rs1049255), A8897T (rs4782390), and G383A (rs3794624).

In conclusion, we found a positive correlation between the p22phox C242T gene polymorphism and ODN susceptibility. The T allele carriers may be ascribed to the increased susceptibility. This conclusion may facilitate individual diagnosis of ODN and guide therapeutic strategies. Nevertheless, considering the above limitations, further studies are necessary to elucidate the relationship between the p22phox C242T gene polymorphism and ODN susceptibility.

KEY MESSAGE

1. The essential subunit for nicotinamide adenine dinucleotide phosphate, p22phox C242T gene polymorphism (rs4673) was positively correlated with overt diabetic nephropathy (ODN).

2. The T allele carriers of p22phox C242T gene polymorphism might be predisposed to ODN.

Notes

No potential conflict of interest relevant to this article was reported.

Acknowledgements

This work was funded by the National Natural Science Foundation of China (NSFC 81100073 to Dr Yan-yan Li), Excellent Young and Middle-Aged Teachers Assistance Program of Nanjing Medical University for Dr Yan-yan Li (2013-2015, JX2161015034), Jiangsu Overseas Research & Training Program for University Prominent Young & Middle-aged Teachers and Presidents (2014), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). Thank all our colleagues working in the Department of Geriatrics, the First Affiliated Hospital of Nanjing Medical University.

Supplementary Materials

Supplementary Table 1.

PRISMA 2009 checklist. Adapted from Moher et al. [12].

Supplementary Figure 1.

PRISMA 2009 flow diagram. Adapted from Moher et al. [12]. HWE, Hardy-Weinberg equilibrium; ODN, overt diabetic nephropathy

References

1. Adarkwah CC, Gandjour A, Akkerman M, Evers SM. Cost-effectiveness of angiotensin-converting enzyme inhibitors for the prevention of diabetic nephropathy in the Netherlands: a Markov model. PLoS One 2011;6e26139.

2. Foley RN, Collins AJ. The growing economic burden of diabetic kidney disease. Curr Diab Rep 2009;9:460–465.

3. Tzamaloukas AH, Murata GH. Prevention of nephropathy in patients with type 2 diabetes mellitus. Int Urol Nephrol 2005;37:655–663.

4. Conserva F, Pontrelli P, Accetturo M, Gesualdo L. The pathogenesis of diabetic nephropathy: focus on microRNAs and proteomics. J Nephrol 2013;26:811–20.

5. Bhatt AP, Gupta A, Vibha D, Sharma A, Mahajan S. Florid diabetic complications in impaired glucose tolerance. Saudi J Kidney Dis Transpl 2013;24:86–88.

6. Najafi M, Alipoor B, Shabani M, Amirfarhangi A, Ghasemi H. Association between rs4673 (C/T) and rs13306294 (A/G) haplotypes of NAD(P)H oxidase p22phox gene and severity of stenosis in coronary arteries. Gene 2012;499:213–217.

7. Perianayagam MC, Liangos O, Kolyada AY, et al. NADPH oxidase p22phox and catalase gene variants are associated with biomarkers of oxidative stress and adverse outcomes in acute renal failure. J Am Soc Nephrol 2007;18:255–263.

8. Hodgkinson AD, Millward BA, Demaine AG. Association of the p22phox component of NAD(P)H oxidase with susceptibility to diabetic nephropathy in patients with type 1 diabetes. Diabetes Care 2003;26:3111–5.

9. Liu LY, Zou DJ, Zhang J. Association between NADPH oxidase subunit p22phox gene polymorphism and type 2 diabetic nephropathy patients in Shanghai. Chin J Nephrol 2006;22:84–87.

10. Yang ST, Liu H, Wang YM, Song DP. The polymorphism of C242T of p22phox gene related to diabetic nephropathy. J Kunming Med Univ 2006;27:111–115.

11. Lim SC, Goh SK, Lai YR, et al. Relationship between common functional polymorphisms of the p22phox gene (-930A > G and +242C > T) and nephropathy as a result of type 2 diabetes in a Chinese population. Diabet Med 2006;23:1037–1041.

12. Moher D, Liberati A, Tetzlaff J, Altman DG; PRISMA Group. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med 2009;6e1000097.

13. Cochran WG. The effectiveness of adjustment by subclassification in removing bias in observational studies. Biometrics 1968;24:295–313.

14. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials 1986;7:177–188.

15. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst 1959;22:719–748.

16. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ 1997;315:629–634.

17. Santos KG, Canani LH, Gross JL, Tschiedel B, Souto KE, Roisenberg I. Relationship of p22phox C242T polymorphism with nephropathy in type 2 diabetic patients. J Nephrol 2005;18:733–738.

18. Hong YR, Han L, Huang J, Kleinberg ME. Elimination of false negative results in the two-hybrid system in the phagocyte NADPH oxidase. Kaohsiung J Med Sci 1996;12:301–305.

19. Inoue N, Kawashima S, Kanazawa K, Yamada S, Akita H, Yokoyama M. Polymorphism of the NADH/NADPH oxidase p22 phox gene in patients with coronary artery disease. Circulation 1998;97:135–137.

20. Guzik TJ, West NE, Black E, et al. Functional effect of the C242T polymorphism in the NAD(P)H oxidase p22phox gene on vascular superoxide production in atherosclerosis. Circulation 2000;102:1744–1747.

21. Liu CC, Karimi Galougahi K, Weisbrod RM, et al. Oxidative inhibition of the vascular Na+-K+ pump via NADPH oxidase-dependent β1-subunit glutathionylation: implications for angiotensin II-induced vascular dysfunction. Free Radic Biol Med 2013;65:563–572.

22. Cahilly C, Ballantyne CM, Lim DS, Gotto A, Marian AJ. A variant of p22(phox), involved in generation of reactive oxygen species in the vessel wall, is associated with progression of coronary atherosclerosis. Circ Res 2000;86:391–395.

23. Ito D, Murata M, Watanabe K, et al. C242T polymorphism of NADPH oxidase p22 PHOX gene and ischemic cerebrovascular disease in the Japanese population. Stroke 2000;31:936–939.

24. Ji Z, Lu BX, Wang XY, Pan SY. C242T p22phox gene polymorphism in the population of Han nationality in Guangdong province. Di Yi Jun Yi Da Xue Xue Bao 2003;23:387–388. 390.

25. Li A, Prasad A, Mincemoyer R, et al. Relationship of the C242T p22phox gene polymorphism to angiographic coronary artery disease and endothelial function. Am J Med Genet 1999;86:57–61.

26. Saha N, Sanghera DK, Kamboh MI. The p22 phox polymorphism C242T is not associated with CHD risk in Asian Indians and Chinese. Eur J Clin Invest 1999;29:999–1002.

Article information Continued

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Table 1.

Characteristics of the investigated studies of the association of the p22phox C242T gene polymorphism and ODN

Study	Region	DM type	Ethnicity	ODN			Control			Matching criteria	Sample size, ODN/control
Study	Region	DM type	Ethnicity	CC	CT	TT	CC	CT	TT	Matching criteria	Sample size, ODN/control
Hodgkinson et al. (2003) [8]	UK	T1DM	Caucasian	22	56	39	41	25	4	Age, sex, duration	117/70
Santos et al. (2005) [16]	Brazil	T2DM	Caucasian	54	69	25	92	80	24	Age, BMI, HbA1c	148/201
Lim et al. (2006) [11]	China	T2DM	Chinese	259	46	1	268	36	2	Age, BMI, HbA1c, duration	306/306
Liu et al. (2006) [9]	China	T2DM	Chinese	52	6	13	101	2	0	Age, sex, BMI	71/103
Yang et al. (2006) [10]	China	T2DM	Chinese	30	32	3	49	21	0	BMI, HbA1c	65/70

The polymerase chain reaction-restriction fragment length polymorphism genotyping method and case-control study design were adopted in all of the above studies.

ODN, overt diabetic nephropathy; DM, diabetes mellitus; T1DM, type 1 diabetes mellitus; T2DM, type 2 diabetes mellitus; BMI, body mass index; HbA1c, glycated hemoglobin.

Genetic model	Pooled OR (95% CI)	p value	Literature no.	No. of ODN case	No. of control group	P_{heterogeneity} (I²%)
Allelic genetic model	2.760 (1.400–5.450)	0.004^a	5	707	745	<0.00001^a (88.5)
Recessive genetic model	5.080 (1.020–25.430)	0.05^a	5	707	745	0.0003^a (80.9)
Dominant genetic model	1.700 (1.167–2.477)	0.006^a	5	707	745	<0.00001^a (80.7)
Homozygous genetic model	3.900 (1.022–14.889)	0.046^a	5	707	745	0.001^a (77.2)
Heterozygous genetic model	1.523 (1.167–1.986)	0.002^a	5	707	745	0.07 (53.8)
Additive genetic model	2.019 (1.232–3.309)	0.005^a	5	707	745	<0.00001^a (87.2)

Genetic model	Pooled OR (95% CI)	p value	Literature no.	No. of ODN case	No. of control group	P_{heterogeneity} (I²%)
Allelic genetic model	2.320 (1.160–4.650)	0.02^a	4	590	675	0.0002^a (85.0)
Recessive genetic model	1.675 (0.952–2.947)	0.074^a	4	590	675	0.056^a (60.3)
Dominant genetic model	2.317 (1.186–4.527)	0.014^a	4	590	675	0.004^a (77.2)
Homozygous genetic model	2.062 (1.128–3.767)	0.019^a	4	590	675	0.053 (60.9)
Heterozygous genetic model	1.612 (1.210–2.165)	0.001^a	4	590	675	0.200 (35.3)
Additive genetic model	2.343 (1.152–4.764)	0.019^a	4	590	675	<0.00001^a (85.3)