Polymorphism of Nitric Oxide Synthase 1 Affects the Clinical Phenotypes of Ischemic Stroke in Korean Population.

OBJECTIVE: To investigate whether four single nucleotide polymorphisms (SNPs) rs2293054 [Ile734Ile], rs1047735 [His902His], rs2293044 [Val1353Val], rs2682826 (3'UTR) of nitric oxide synthase 1 (NOS1) are associated with the development and clinical phenotypes of ischemic stroke.
METHODS: We enrolled 120 ischemic stroke patients and 314 control subjects. Ischemic stroke patients were divided into subgroups according to the scores of the National Institutes of Health Stroke Survey (NIHSS, <6 and ≥6) and Modified Barthel Index (MBI, <60 and ≥60). SNPStats, SNPAnalyzer, and HelixTree programs were used to calculate odds ratios (ORs), 95% confidence intervals (CIs), and p-values. Multiple logistic regression models were performed to analyze genetic data.
RESULTS: No SNPs of the NOS1 gene were found to be associated with ischemic stroke. However, in an analysis of clinical phenotypes, we found that rs2293054 was associated with the NIHSS scores of ischemic stroke patients in codominant (p=0.019), dominant (p=0.007), overdominant (p=0.033), and log-additive (p=0.0048) models. Also, rs2682826 revealed a significant association in the recessive model (p=0.034). In allele frequency analysis, we also found that the T alleles of rs2293054 were associated with lower NIHSS scores (p=0.007). Respectively, rs2293054 had a significant association in the MBI scores of ischemic stroke in codominant (p=0.038), dominant (p=0.031), overdominant (p=0.045), and log-additive (p=0.04) models.
CONCLUSION: These results suggest that NOS1 may be related to the clinical phenotypes of ischemic stroke in Korean population.

INTRODUCTION

Stroke is a serious neurological disease which is the leading cause of severe disability and the second most common cause of death [12]. Stroke can be subdivided into ischemic stroke (IS) and hemorrhagic stroke. IS accounts for approximately 80% of all strokes and results from vascular occlusion [3]. The IS may be a highly multifactorial disease consisting of a group of heterogeneous disorders with multiple genetic and environmental risk factors [4]. Several lines of evidence have consistently suggested a genetic influence on the risk and prognosis of IS [5]. Although genetic polymorphisms of several candidate genes for IS have been studied by case control analysis, few associations have been consistently replicated [5].

The role of nitric oxide (NO) was investigated in in vivo models of cerebral ischemia [6]. NO is produced using L-arginine as the substrate and nitric oxide synthase (NOS) is responsible for synthesis. In human beings, three isoforms of NOS, which are NOS1, NOS2, NOS3 are identified. NOS1 is the so-called neuronal NOS and is mainly expressed in the brain. It is reported to play a role in central nervous system synaptic plasticity, reduced sympathetic activity and reduced tone in various smooth muscles. NOS2 is known as inducible NOS and is mainly induced in macrophages, which represent a cytotoxic principle. NOS3 is so-called endothelial NOS and is expressed in endothelial cells. It reportedly dilates all types of blood vessels [7]. In the early phase after ischemia, NO is produced by the constitutive endothelial and neuronal isoforms of nitric oxide synthase (NOS3 and NOS1) [6]. Although NOS3 appears beneficial via vasodilatation and inhibition of leukocyte adhesion and platelet aggregation, NOS1 was demonstrated to be deleterious in cerebral ischemia. Moreover it is now thought that NO induced neuronal death is mainly mediated through the formation of very reactive peroxynitrite, resulting from the reaction between NO and superoxide anion [6]. In terms of variation in NO availability, a risk for IS could be explained in part by NOS1 genetic polymorphism [8910]. Although previous studies reported a relationship between NOS1 genetic polymorphism and the susceptibility of IS, important clinical phenotypes were not included in their analysis [8910]. Therefore, we evaluated not only the genetic association between NOS 1 polymorphisms and the susceptibility of IS, but also clinical phenotypes of IS patients in a Korean population.

MATERIALS AND METHODS

Subjects and clinical phenotypes

IS patients were recruited at the Stroke Center of Kyung Hee University Hospital and Emergency Department in Kyung Hee University Medical Center from October 2007 to December 2009. One hundred twenty-one IS patients (67 males, 53 females; 65.8±12.2 years) 1 week to 1 month from onset were enrolled (Table 1). Subjects with transient ischemic attack, traumatic brain injury, brain tumors, congenital brain disorders, cardiovascular malformation, chronic inflammatory diseases, auto-immune diseases and neuropsychiatric disorders were excluded. Each stroke patient was diagnosed through neurological examinations, computed tomography (CT), magnetic resonance imaging (MRI), angiography, and duplex sonography. All stroke patients were classified into two clinical phenotypes using the National Institutes of Health Stroke Survey (NIHSS) score and the Modified Barthel Index (MBI) score. For the prognosis of stroke patients, the severity of 13 neurological symptoms were assessed using the NIHSS score. For hospitalization, the quality of 10 activities of general life was estimated by the MBI score. Age-matched control subjects (n=314; 176 males, 138 females; 63.1±9.4 years) were selected through a general health check-up program (Table 1). Participants with symptoms or history of neurological diseases, ischemic heart diseases, or any severe diseases were excluded. This study was approved by the ethics review committee of Medical Research Institute, Kyung Hee University College of Medicine, and written informed consent were obtained from both the IS patients and control subjects.

Table 1

Patients' characteristics

Values are presented as number or mean±standard deviation.

IS, ischemic stroke; NIHSS, National Institute of Health Stroke Scale; MBI, Modified Barthel Index.

SNP selection and genotyping

We selected four, single nucleotide polymorphisms (SNPs) with greater than 0.1 heterozygosity among SNPs located in the coding regions or 3'-untranslated region (UTR) of the NOS1 gene which were searched at the dbSNP database (http://www.ncbi.nlm.nih.gov/SNP; dbSNP Build 131) [1112]. SNPs with unknown heterozygosity, minor allele frequency below 4% were excluded. Finally, four SNPs were selected and analyzed: rs2293054 [Ile734Ile], coding-synonymous; rs1047735 [His902His], coding-synonymous; rs2293044 [Val1353Val], codingsynonymous; rs2682826, 3'-UTR. Peripheral blood samples were collected from all participants and then preserved in a –70℃ refrigerator for further laboratory procedures. The polymerase chain reaction (PCR) comprised 35 cycles at 94℃ for 30 seconds, 58℃ for 30 seconds, 72℃ for 1 minute, and 1 cycle at 72℃ for 5 minutes for the final reaction. PCR products were sequenced using an ABI PRISM 3730XL analyzer (PE Applied Biosystems, Foster City, CA, USA). Sequence data were analyzed using the SeqMan II software (DNASTAR Inc., Madison, WI, USA).

Statistical analysis

Hardy–Weinberg equilibrium (HWE) was assessed using a chi-squared test. The HWE is a principle stating that the genetic variation in a population will remain constant from one generation to the next in the absence of disturbing factors. The HWE can be disturbed by several factors including mutations, natural selection, nonrandom mating, and genetic drift. A significant result (p<0.05) means that controls are not in HWE [13]. SNPAnalyzer-Pro (ISTECH, Goyang, Korea), HelixTree (Golden Helix, Bozeman, United States), and SNPStats (http://bioinfo.iconcologia.net/index.php?module=Snpstats) were used to analyze genetic data. Multiple logistic regression models (codominant 1, codominant 2, dominant, recessive, overdominant, and log-additive) were performed to determine the odds ratios (ORs), 95% confidence intervals (CIs) and p-values, (age and gender were controlled as co-variables). For the statistical power, the required case size in each SNP was calculated using the genetic power calculator (http://pngu.mgh.harvard.edu/~purcell/gpc/cc2.html). The statistical significance level was considered as a p-value below 0.05.

Genotyping of each SNP

For the analysis of the genetic association between NOS1 polymorphisms and susceptibility to IS, multiple logistic regression analysis, adjusted for age and gender, was performed using the following: codominant 1 (major allele homozygotes vs. heterozygotes), codominant 2 (major allele homozygotes vs. minor allele homozygotes), dominant (major allele homozygotes vs. heterozygotes/minor allele homozygotes), recessive (major allele homozygotes/heterozygotes vs. minor allele homozygotes), overdominant (major allele homozygotes/minor allele homozygotes vs. heterozygotes) and log-additive (major allele homozygotes vs. heterozygotes vs. minor allele homozygotes).

In the codominant model, every genotype gives a different and non additive risk. Therefore, we compared the major allele homozygotes that is the most frequent allele, to heterozygotes in codominant 1 and minor allele homozygotes in codominant 2. In the dominant model a single copy of an allele is enough to modify the risk and heterozygous and homozygous genotypes have the same risk. In the recessive model two copies of the allele are necessary to change the risk. In the overdominant model heterozygotes are compared to both allele homozygotes. In the log-additive model each copy of the allele modifies the risk in an additive form.

RESULTS

Clinical characteristics

The clinical characteristics of patients with IS and control subjects are shown in Table 1. The age of IS and control subjects (mean±standard deviation) were 65.8±12.2 and 63.1±9.4 years, respectively. Stroke patients were divided according to clinical phenotypes which are NIHSS (scores of 13 neurological symptoms; <6 or ≥6), MBI (scores of 10 activities of general life; <60 or ≥60), hypertension (present or absent) and diabetes mellitus (present or absent). The number of IS patients with NIHSS scores <6 and ≥6 were 56 and 56, respectively. The number of IS patients with MBI scores <60 and ≥60 were 70 and 25, respectively. The number of IS patients with hypertension present and absent were 80 and 35, respectively. The number of IS patients with diabetes mellitus present and absent were 45 and 71, respectively. No control subjects had hypertension or diabetes mellitus (Table 1).

Analysis of IS susceptibility

The genotype and allele distributions of the four selected SNPs are presented in Tables 2, 3, 4. The observed genotype distributions of the SNPs were in HWE (p>0.05). The genotype and allele distributions of the 4 SNPs among the IS group and control subjects are presented in Table 2. All 4 SNPs tested showed no statistically significant association between the IS and control subjects (Table 2).

Table 2

Genotype frequencies of SNPs of nitric oxide synthase 1 (NOS1) genes in ischemic stroke patients and control subjects

Values are presented as number (%).

SNP, single nucleotide polymorphism; IS, ischemic stroke; OR, odds ratio; CI, confidence interval; cds-synon, codon synonymous; UTR, untranslated region; NA, not applicable.

The p-values were calculated using logistic regression analyses adjusting for sex and age.

Table 3

Multiple logistic regression analysis of nitric oxide synthase 1 (NOS1) gene polymorphisms in ischemic stroke patients using NIHSS scores of <6 and ≥6

Values are presented as number (%).

SNP, single nucleotide polymorphism; NIHSS, National Institute of Health Stroke Scale; OR, odds ratio; CI, confidence interval; cds-synon, codon synonymous; UTR, untranslated region.

The p-values were calculated using logistic regression analyses adjusting for sex and age.

*p<0.05 statistically significant.

Table 4

Multiple logistic regression analysis of nitric oxide synthase 1 (NOS1) gene polymorphisms in ischemic stroke patients with MBI scores of <60 and ≥60

Values are presented as number (%).

SNP, single nucleotide polymorphism; MBI, Modified Barthel Index; OR, odds ratio; CI, confidence interval; cds-synon, codon synonymous; UTR, untranslated region.

The p-values were calculated using logistic regression analyses adjusting for sex and age.

*p<0.05 statistically significant.

Analysis of clinical phenotypes

All stroke patients were classified into clinical phenotypes according to the NIHSS score (<6 or ≥6) and the MBI score (<60 or ≥60). At first, we assessed the relationship between NOS1 SNPs and NIHSS scores of <6 or ≥6. In the analysis of NIHSS, we observed a difference between NIHSS <6 group and NIHSS ≥6 group (Table 3). The genotype frequencies of the CC, CT, and TT for rs2293054 were 0.48, 0.45, and 0.07 in the NIHSS <6 group, while they were 0.73, 0.25, and 0.02 in the NIHSS ≥6 group. The CT genotype frequency of rs2293054 differed between the NIHSS <6 group and NIHSS ≥6 group (44.6% vs. 25%). The rs2293054 was associated with IS in the codominant model (OR, 0.38; 95% CI, 0.17–0.85; p=0.019), dominant model (OR, 0.34; 95% CI, 0.16–0.76; p=0.007), overdominant (OR, 0.42; 95% CI, 0.19–0.94; p=0.033), and log-additive model (OR, 0.38; 95% CI, 0.19–0.77; p=0.0048). In the allele frequency analysis, we also found that the T alleles of rs2293054 were associated with lower NIHSS scores (OR, 0.40; 95% CI, 0.21–0.78; p=0.007) (Table 3). Also, the genotype frequencies of the CC, CT, and TT for rs2682826 were 0.45, 0.52, and 0.04 in the NIHSS <6 group, while they were 0.48, 0.38, and 0.14 in the NIHSS ≥6 group. The CT genotype frequency of rs2682826 was not different between the NIHSS <6 group and NIHSS ≥6 group (51.8% vs. 37.5%). The rs2682826 was not associated with NIHSS scores in the codominant model (OR, 0.67; 95% CI, 0.31–1.47; p=0.32), dominant model (OR, 0.87; 95% CI, 0.41–1.84; p=0.72), overdominant model (OR, 0.55; 95% CI, 0.26–1.18; p=0.12), and log-additive model (OR, 1.20; 95% CI, 0.67–2.15; p=0.53). But, rs2682826 revealed a significant association in recessive model (OR, 4.76; 95% CI, 0.95–23.77; p=0.034). Secondly, in multiple logistic regression analysis of NOS1 polymorphisms in IS patients with MBI scores of <60 or ≥60, rs2293054 had significant association within the codominant model (OR, 2.98; 95% CI, 1.06–8.32; p=0.038), dominant model (OR, 2.99; 95% CI, 1.09–8.21; p=0.031), overdominant model (OR, 2.77; 95% CI, 1.02–7.57; p=0.045), and log-additive model (OR, 2.47; 95% CI, 1.04–5.88; p=0.04), respectively. The rs1047735 and rs2293044 were not statistically associated with NIHSS scores (<6 or ≥6), and MBI scores (<60 or ≥60).

DISCUSSION

Several polymorphisms of the neuronal NOS1 gene have been reported to be associated with impulsivity, aggression, and suicide attempts [14], high-grade renal cell carcinoma [15], end stage renal disease susceptibility [16], asthma and cystic fibrosis [17181920], infantile hypertrophic pyloric stenosis [21], and systolic heart failure [9]. In the viewpoint of polymorphisms related to IS, NOS3 appears beneficial via vasodilatation and inhibition of leukocyte adhesion and platelet aggregation. By contrast, NOS1 and NOS2 were demonstrated to be deleterious in cases of cerebral ischemia [6]. Recent studies suggested that NOS1 may have a crucial vasculoprotective role in atherosclerosis and is involved in blood pressure control [22]. NOS1 influences stroke susceptibility through very different mechanisms from NOS3. A significant association of the NOS1 gene, with stroke susceptibility, have been reported [10]. Manso et al. [10] showed that associated SNPs and haplotypes of NOS1 were located in the critical gene regions, namely the NOS domain and the PDZ domain, required to anchor NOS1 to the cellular membrane. Four SNPs (rs2293050, rs2139733, rs7308402 and rs1483757) were found to have a significant relationship with the development of IS. However, no relationship was found with the clinical severity of IS. As mentioned above, the aim of our study was to investigate different SNPs from previous studies, and evaluate not only the genetic association between NOS1 polymorphisms and the susceptibility of IS, but also the clinical phenotypes of IS patients in a Korean population.

We classified stroke patients into clinical phenotypes according to the NIHSS score (<6 or ≥6) and the MBI score (<60 or ≥60). In a previous study, patients with MBI score higher than 60 required shorter home stays. In this study, an MBI score of 60 appeared to be a reference point to whether patients are dependent or independent or requires only partial assistance [23]. In terms of NIHSS score, score below 6 predicted a good outcome in a study consisting of 1,281 stroke patients [24].

The NOS1 (also known as nNOS and IHPS1), which is located on chromosome 12q24.2-q24.31, produce NO in a reaction that converts arginine and oxygen into citrulline and NO [25]. NO is a reactive free radical, which acts as a biologic mediator in several processes, including neurotransmission, antimicrobial, and antitumoral activities. In the brain and peripheral nervous system, NO displays many properties of a neurotransmitter and has been implicated in neurotoxicity associated with stroke and neurodegenerative diseases, neural regulation of smooth muscle including peristalsis, and penile erection. NOS1 protein (P29475, UniProtKB/Swiss-Prot) consists of 1,434 amino acids, and has a molecular mass of 160,970 Da. Mammalian NOS consist of an oxygenase domain (heme and H4B binding) and a reductase domain related by sequence, structure and function to NADPH-dependent microsomal cytochrome P450 reductase. The NOS reductase domain can be further divided into FMN and FAD binding domains, linked together by a single polypeptide chain. The FMN domain is connected to the oxygenase domain, via a segment of polypeptide about 30 amino acids in length, which forms a binding site for calmodulin (amino acids from 730 to 750) [8]. Our results suggest that the coding SNP (rs2293054, exon 16, 734th amino acid of NOS1 protein) (http://www.uniprot.org/uniprot/P29475) in the NOS1 gene contributes to the susceptibility to the clinical phenotype such as NIHSS and MBI in IS patients (Tables 3, 4). In terms of location on NOS1 gene, the rs2293054 SNPs may be related to the function of the calmodulin domain of NOS1.

Hancock et al. [26] found that significant associations with Parkinson's disease were detected for the NOS1 SNPs especially rs2682826 (3'UTR, untranslated region). In our study, we found an association between rs2682826 and NIHSS scores in the recessive model (OR, 4.76; 95% CI, 0.95–23.77; p=0.034). The genotype frequencies and T-allele of the rs2293054 SNP for the NOS1 in Parkinson's disease or stroke patients was not associated with disease susceptibility or was not reported [1026]. In another study, the relationship between rs2293054 SNP and cognitive aging was tested but no relationship could be found [27]. To our knowledge these two studies are the only ones to investigate the relationship of rs2293054 SNP and clinical diseases [2627]. As mentioned above, none of these studies found significant results of rs2293054 SNP. However, we found the association between rs2293054 and NIHSS or MBI scores in the codominant model, dominant model, overdominant model and log-additive model. Therefore, the T-allele of the rs2293054 may contribute to the pathogenesis of clinical phenotypes in patients with IS.

In previous studies conducted in the Chinese population, rs1047735 SNP showed no relationship with glioma and infertility [2829]. It also showed no relationship with Parkinson's disease [26]. Furthermore, there are no studies investigating the clinical role of rs2293044 SNP. Like these previous studies, we also found no clinical importance for these two SNPs. Further study is required to determine whether rs1047735 SNP and rs2293044 SNP play a clinical role.

In our study, 121 IS patients were enrolled. Compared to previous studies evaluating SNPs, the sample size is small and therefore our results need to be replicated in other populations. Another limitation of this study is that although NIHSS and MBI are one of the most widely used scales to measure stroke outcome, stroke outcome is influenced by multiple factors including rehabilitation treatment and pharmacotherapy, which were not assessed in our study.

In conclusion, we observed the association between clinical phenotype and coding polymorphisms (rs2293054, Ile734Ile) in the NOS1 gene in the Korean population. These results suggest that NOS1 gene may affect susceptibility to clinical phenotypes of IS.