Combination analysis of NOS3, ABCB1 and IL23R polymorphisms with alcohol-induced osteonecrosis of the femoral head risk in Chinese males

Background Common variants of multiple genes played a crucial role in osteonecrosis of the femoral head (ONFH) onset which was proved by many previous reports. We hypothesized that polymorphisms in NOS3, ABCB1 and IL23R were related to individual differences in alcohol sensitivity and the development of alcohol-induced ONFH. Methods In this case-control study, we evaluated 8 SNPs in three genes in the Chinese Han population including 355 male cases and 355 healthy male controls. These SNPs were genotyped by Sequenom MassARRAY RS1000. To identify their relationship with alcohol-induced ONFH susceptibility using χ2 test and genetic model analysis. Results We found an association with alcohol-induced ONFH susceptibility for 4 SNPs (rs743506, rs3918184, rs13233308 and rs6693831) in three genes after adjusted by age. The genotype “G/A” of rs743506 in NOS3 gene acts as a risk factor in genotype (P = 0.003), dominant (P = 0.048), recessive (P = 0.005) and additive model(P = 0.006); The genotype “T/C” of rs3918184 in NOS3 gene acts as a risk factor in genotype (P = 0.012) and recessive model (P = 0.009); The genotype “T/C” of rs13233308 in ABCB1 gene acts as a risk factor in genotype (P = 0.038) and additive model(P = 0.041); The genotype “T/C” of rs6693831 in IL23R gene acts as a protective factor in genotype model (P = 0.046). Conclusions This study provides evidence for three alcohol-induced ONFH susceptibility genes (NOS3, ABCB1 and IL23R) in Chinese males and polymorphisms of them may be associated with alcohol-induced ONFH risk.


INTRODUCTION
Osteonecrosis of the femoral head (ONFH) indicates a debilitating disorder that the bone collapses, happening in a special anatomical site of femoral head and resulting in osteocyte death [1,2]. It mainly affects younger males who is between 30 and 50 years old. In China, there are seven million patients with osteonecrosis and be found annually new patients to increase by 150,000 to 200,000 patients [3]. The seventy-five to ninety percent cases of ONFH have been associated with various risk factors, involving alcohol Research Paper www.impactjournals.com/oncotarget abuse, corticosteroid use, hip trauma and smoking [4]. Increasing of adipose vesicles in the blood circulating is induced these risk factors, adding lipid deposition in osteocytes of the femoral head and resulting in embolism, which influence finally blood flow [5].
Current evidence suggests that being related to the risk of ONFH and eight SNPs in three genes of our research can be identified through the genome-wide association studies [1,6,7]. Nitric oxide is synthesised from L-arginine and has three isoforms of synthases (NOS): endothelial (eNOS) neuronal (nNOS) and inducible(iNOS). The eNOS is expressed in normal adult bone as a constitutive isoform [8]. A previous study in Korean patients indicated that idiopathic osteonecrosis of the femoral head was significantly associated with polymorphism of eNOS gene, but it also indicate that polymorphism in intron 4 of eNOS was not significant differences [9]. Glueck et al. found that the pathogenesis of idiopathic ONFH are related with the intron 4 of eNOS polymorphism and synthetic loss of nitric oxide outcome in African American and Caucasian patients. The ABCB1, is also called adenosine triphosphatebinding cassette B1, can encode the transport protein and made an significant influence on distribution and absorption in the mammalian body [10]. So far a few single nucleotide polymorphisms (SNPs) of ABCB1 gene have been identified, of which mutations in exon 21 (G2677T) and exon 26 (C3435T) are associated with alteration of P-gp expression or function, as recently reviewed by Sakaeda et al. and Fromm et al [11,12]. Interleukin 23 receptor (IL23) is a proinflammatory cytokine. According to previous reports indicated it inhibited osteoclastogenesis by Receptor Activator for Nuclear Factor κB Ligand with T cell affection [13]. Kim et al. discover polymorphisms in the Interleukin 23 receptor(IL23R) gene are associated with ONFH in Korean population [14].
A few studies investigated the association of NOS3, ABCB1 and IL23R polymorphism with the pathogenesis of ONFH. We hypothesized that NOS3, ABCB1 and IL23R genes have crucial roles in the development of alcoholinduced ONFH. To identify whether polymorphisms of these genes to be associated with alcohol-induced ONFH in Chinese males.

RESULTS
In the current study, a total of 355 male cases (median age at diagnosis 44.91 ± 9.85 years) and 355 healthy male controls (median age 46.02 ± 9.61 years) were included. The basic characteristics of cases and controls were showed in Table 1. The primers of the 8 candidate SNPs are revealed in Table 2. Eight SNP loci(rs6693831, rs790631, rs4148749, rs10808072, rs13233308, rs3918227, rs3918184 and rs743506) in the NOS3, ABCB1 and IL23R were evaluated in this study. SNP ID, HWE P value, allele A/B, MAF control/ case, odds ratio, 95 % confidence interval and P value were listed in Table 3. Using x 2 test, we found significant differences in frequency of alleles and rs13233308 in the ABCB1 gene and rs743506 in the NOS3 gene were associated with increased risk of alcohol-induced ONFH risk by allele model analysis (rs13233308, P =0.032, odds ratio [OR]: 1.26, 95% confidence interval [CI]: 1.02-1.55 and rs743506, P = 0.006, OR: 1.39, 95%CI: 1.09-1.78). All of the SNPs were in Hardy-Weinberg equilibrium (HWE) in the control population of this study.

DISCUSSION
Non-traumatic ONFH has been diffusely awaked as a pathological define with numerous etiologies, the precise    Freq, frequency; OR, odd ratio; CI, confidence interval. *P < 0.05, statistical significance: †Adjusted by age pathogenesis of osteonecrosis leaves to be evaluated. Alcoholism and chronic steroid use were considered to be the essential risk factors induced osteonecrosis. Genetic mutation of alcohol-metabolizing enzyme genes has been associated with alcohol-induced osteonecrosis and polymorphism of NOS3, ABCB1 and IL23R gene are related with developing osteonecrosis [6,[15][16][17]. The NOS3 gene located in chromosome 7q36.1 involves two SNPs ("G/A" in rs743506 and "T/C" in rs3918184) related with increased alcohol-induced ONFH risk. Endothelial nitric oxide synthase is the remarkable NOS isoform expressed in normal bone and eNOS adjusts bone resorption and formation and oestrogen are dependent on NO production via eNOS in the skeletal system [18]. Excessive NO production was found in various rheumatic diseases, including rheumatoid arthritis, systemic lupus erythematosus, osteoarthritis and vasculitis [8]. Güler S et al. found the genotype "GA" of rs743506 showed risk effect was associated with aura in migraine patients [19]. It proved that mutation of rs743506 would affect function of NO in cerebral blood flow regulation and include the activation of nociceptors in the vascular system and the release of vasoactive neuropeptides in the neurogenic inflammatory response [20]. Meanwhile, Juan P et al. found genotype "T/C" of rs3918184 be related with increased risk of hypertension, preeclampsia and stroke [21]. Combination with our research, rs743506 and rs3918184 associated with increased with alcoholinduced ONFH risk and it may also be effect NO in the femoral head blood flow regulation. The NOS is fundamentally expressed in vascular endothelium and had been proved to shorten NO norms in human plasma [22]. Our study suggests that the genotype "GA" of rs743506 and genotype "T/C" of rs3918184 polymorphisms may be risk factors for alcohol-induced ONFH, and NO brought in fundamentally decided eNOS may make a significant effect on the pathogenesis of alcohol-induced ONFH.
The ABCB1 gene, also known as Multidrugresistant transporter-1(MDR1), was located in chromosome 7q21.12 and the genotype "T/C" of rs13233308 related with increased alcohol-induced ONFH risk. This gene spans 4.5 kb encoding a 1280 amino acid (aa) polypeptide. After glycosylation of the amino acid polypeptide, a 170-kDa of protein named P-glycoprotein (P-gp) can be formed. P-gp is occurred in various tissue cells and is required for pumping lipophilic drugs (including glucocorticoids) out of the cells. This causes cytotoxic damage and the intracellular drug concentration to decrease is subsequently reduced. In addition, P-gp, encoded by the human ABCB1 gene, is an important determinant in absorption, elimination of drugs and tissue targeting [23]. Xue et al. found polymorphisms of ABCB1 gene be related with decreased risk of steroid-induced avascular necrosis of the femoral head in a Chinese population [16]. Zhang et al. also found ABCB1 gene polymorphisms associated with decrease risk for steroid-induced osteonecrosis of the femoral head in Chinese population [1]. To date, 28 SNPs in the ABCB1 gene have been reported at 27 positions in Caucasians and Africans [24,25]. Apparently, these results was inconsistent with ours. Meanwhile, we must notice that ethnic differences may be considered to take an important effect to ABCB1 gene and ONFH. The homozygosity of both the wildtype and variant alleles of MDR1 (C1236 T, G2677T/A, and C3435T) in healthy control subjects was similar to that reported by other Asian studies, which were performed with healthy Chinese, Japanese, and Korean subjects. However, the allele frequencies were significantly higher than those observed in European ethnicities, such as German, Russian, and Serbian [16]. Thus, further studies related with in other ONFH populations are required to confirm and realize ABCB1 gene different mechanisms in steroidinduced and alcohol-induced ONFH.
The IL23R gene was located in chromosome 1p31.3 and the genotype "T/C" of rs6693831 related with decreased alcohol-induced ONFH risk. Interleukin 23 receptor, a heterodimer protein, was composed of two subunits: IL-23R and IL-12Rb1. IL-23R, expressed on activated T cells, dendritic cells, macrophages and natural killer (NK) cells, binds particularly to IL-23, produced primarily by dendritic cells and activated macrophages from peripheral tissues [26]. Kim et al. found polymorphisms of the IL23R gene were expressively related with ONFH among the Korean population. Among others, the polymorphism of the IL23R gene were particularly related with idiopathic ONFH risk. However, there was no significant associated with alcohol-induced and steroid-induced ONFH [14]. Apparently, this results were inconsistent compared with our results. In our study, rs6693831 of IL23R gene were related with decrease alcohol-induced ONFH risk in Chinese males. Alcoholism, which can influence fat metabolism and backlogging cell pressure, which considers absorption of the blood flood, were considerable risk factors with osteonecrosis [27]. We may thought whether the effect of alcohol abuse was stronger than genetic factors in regards to alcohol-induced ONFH patients. Meanwhile, the difference of results must be considered the effect related ethnic. In further studies, we need larger different ethnic populations to verify our results.
In this case-control study, some limitations were intrinsic and must be marked. First, to avoid selection bias, alcohol-induced ONFH cases and healthy controls came from the same hospital. However, this bias not necessarily to be of consequence because of demographic variables. Secondly, the sample size (355 male cases and 355 male controls) is not large enough in our work. We performed a power analysis and just found the power of rs743506 in NOS3 gene was 0.91. The power of seven surplus SNPs was < 0.75. We considered that the number of samples was small in association studies and confirming our finding required a larger sample size [28].
Our research provides a new evidence for a relationship between NOS3, ABCB1 and IL23R gene and alcohol-induced ONFH in Chinese males onset, which may shed light on the etiology of alcohol-induced ONFH. Functional studies are further required to evaluate the correlation of genotype and phenotype in a large cohort of various ethnicities.

Ethics statement
The use of the protocol in this study was strictly affirmed to the principles expressed in the Declaration of Helsinki and were approved by the Ethical Committee of Zhengzhou Traditional Chinese Medicine Traumatology Hospital. Informed consent was signed from all of the participants.

Study participants
Three hundred and fifty male patients with alcoholinduced ONFH were consecutively enrolled in this study at the Zhengzhou Traditional Chinese Medicine Traumatology Hospital from June 2014 to January 2016. Patients with a history of ethanol consumption of at least 400 ml per week were categorized under alcohol-induced osteonecrosis [29]. The development of alcohol-induced ONFH was diagnosed according to assessment by X-rays, magnetic resonance imaging (MRI), and bone scans [2]. The ONFH was present in one hip in 135 patients and in both hips in 220 patients (440 hips). Patients with a demonstrable history of direct trauma or with the possibility of a combination of causes were excluded.
We also enrolled 355 healthy male controls between July 2014 and January 2016 based on medical examination at Zhengzhou Traditional Chinese Medicine Traumatology Hospital. The controls had a history of ethanol consumption of at least 400 ml per week, however, they had no alcohol-induced ONFH and other related diseases, no history of thromboembolic events and no symptoms of hip disease. All participants were restricted to Chinese Han population who lived in Zhengzhou city and its surrounding areas.

Genotyping
Genomic DNA was extracted and purified from the whole blood of all the participants using a kit (GoldMag, China) and their concentration was measured by spectrometry (DU530UV/VIS spectrophotometer; Beckman Instruments, Fullerton, CA). All PCR reactions were operated in a 25ul volume containing 10pmol of each primer, 1X PCR buffer, 10mM dNTP, 25 ng of genomic DNA, and 1U of Taq DNA polymerase (Solgent, Daejeon, Korea). The conditions of PCR amplification were as follows: initial denaturation at 95°C for 15 min followed by 35 cycles of denaturation at 95°C for 20 sec, annealing at 55°C for 40 sec, extension at 72°C for 1 min, and a final extension at 72°C for 5 min. SNP genotyping was performed on the SEQUENOM MassARRAY_Analyzer 4 (Sequenom, Inc., San Diego, CA, USA) using genomic DNA in a single multiplex reaction. According to the manufacturer's instructions, the primers for polymerase chain reaction amplification and single base extension were designed by Sequenom Assay Design 3.0 software (Sequenom, San Diego, CA, USA) [30]. SNP genotyping using the standard protocol recommended by the manufacturer was performed by Sequenom MassARRAY RS1000. Data analyses and management were conducted by Sequenom Typer 4.0 Software [30,31].

Statistical analysis
Microsoft Excel and SPSS 16.0 statistical package (SPSS, Chicago, IL) were used to perform statistical analyses. All P values in this study were two-sided, and P = 0.05 was considered the threshold of whether statistical significance was achieved or not. A significant departure of genotype frequency from Hardy-Weinberg equilibrium (HWE) for each SNP was estimated using SNPStats (http://bioinfo.iconcologia.net/snpstats/start.htm). The genotype frequencies of cases and controls were calculated using χ2 test [32]. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were determined using unconditional logistic regression analysis with adjustment by age [33]. The four genetic models (Genotype, Dominant, Recessive and Additive) were applied using PLINK software (http://pngu.mgh.harvard.edu/purcell/plink/) to assess the association of SNPs with the risk of alcohol-induced ONFH.