Published online before print August 29, 2005, doi: 10.1161/CIRCULATIONAHA.105.544619
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2005;112:1419-1427.)
© 2005 American Heart Association, Inc.
Epidemiology |
Common Genetic Variation at the Endothelial Nitric Oxide Synthase Locus and Relations to Brachial Artery Vasodilator Function in the Community
From the Framingham Heart Study of the National Heart, Lung, and Blood Institute (S.K., M.G.L., R.S.V., C.G., M.J.K., C.N.-C., D.L., C.J.O., E.J.B.), Framingham, Mass; Department of Mathematics and Statistics (M.G.L.), Evans Department of Medicine (R.S.V., J.A.V., E.J.B.), and Whitaker Cardiovascular Institute (R.S.V., J.A.V., E.J.B.), Boston University, Boston, Mass; Cardiovascular Engineering, Inc, Waltham, Mass (G.F.M.); Program in Medical and Population Genetics, Broad Institute of Massachusetts Institute of Technology and Harvard University, Cambridge, Mass (S.K., C.N.-C., S.L.M., A.L.L., J.A.D., C.J.O., J.N.H.); and Cardiology Division (S.K., C.N.-C., C.J.O.), Massachusetts General Hospital, Divisions of Genetics and Endocrinology, Children’s Hospital (J.A.D., J.N.H.), and Department of Genetics (J.N.H.), Harvard Medical School, Boston, Mass.
Correspondence to Emelia J. Benjamin, MD, ScM, Department of Medicine, Boston University School of Medicine, Framingham Heart Study, 73 Mt Wayte Ave, Suite 2, Framingham, MA 01702-5827. E-mail emelia{at}bu.edu
Received February 21, 2005; revision received May 25, 2005; accepted June 14, 2005.
 |
Abstract |
|---|
 Top Abstract IntroductionMethodsResultsDiscussionReferences |
|---|
Background— Sequence variants at the endothelial nitric oxide synthase (NOS3) locus have been associated with endothelial function measures, but replication has been limited.
Methods and Results— In reference pedigrees, we characterized linkage disequilibrium structure at the NOS3 locus using 33 common single nucleotide polymorphisms (SNPs). Eighteen SNPs that capture underlying common variation were genotyped in unrelated Framingham Heart Study participants (49.5% women; mean age, 62 years) with measured brachial artery flow-mediated dilation (n=1446) or hyperemic flow velocity (n=1043). Within 3 defined blocks of strong linkage disequilibrium that spanned NOS3, 11 SNPs captured >80% of common haplotypic variation. Among men, there were nominally significant associations between 8 NOS3 SNPs (minimum P=0.002) and between haplotypes (minimum P=0.002) and either flow-mediated dilation or hyperemic flow velocity. In women, we did not observe significant associations between NOS3 SNPs or haplotypes and endothelial function measures. To correct for multiple testing, we constructed 1000 bootstrapped null data sets and found that empirical probability values exceeded 0.05 for both phenotypes.
Conclusions— A parsimonious set of SNPs captures common genetic variation at the NOS3 locus. A conservative interpretation of our results is that, accounting for multiple testing, we did not observe statistically significant relations between NOS3 sequence variants and endothelial function measures in either sex. The nominal associations of select NOS3 variants with endothelial function in men (unadjusted for multiple testing) should be viewed as hypothesis-generating observations and may merit testing in other cohorts and experimental designs.
Key Words: endothelium • epidemiology • genetics • nitric oxide synthase
|
Introduction |
|---|
|
TopAbstractIntroductionMethodsResultsDiscussionReferences |
|---|
The endothelial isoform of nitric oxide synthase (NOS) (encoded by the NOS3 gene) converts L-arginine to nitric oxide (NO), a potent vasodilator and inhibitor of inflammation and platelet activity.1 In humans, flow-mediated dilation (FMD) and reactive hyperemia have been shown to depend on the bioavailability of endothelium-derived NO.2,3 Impaired endothelium-dependent vasodilation predicts future cardiovascular events (reviewed in detail elsewhere4),5 suggesting that loss of endothelium-derived NO contributes to the pathogenesis and clinical expression of atherosclerosis.6 For these reasons, there is intense interest in defining the mechanisms of endothelial dysfunction, including the possible contributory role of altered expression and/or activity of endothelial NOS.
Experimental and clinical studies suggest that genetic variation in NOS may influence the bioavailability of NO. In mice, targeted deletion of NOS3 leads to endothelial dysfunction and hypertension.7 Young individuals with a parental history of premature cardiovascular disease have impaired brachial artery FMD,8,9 and we recently reported that brachial FMD is modestly heritable in the community.10 Other investigators have evaluated the relations of specific NOS3 genetic variants to endothelial function in small- to modest-sized referral samples, but these studies have yielded inconsistent results.11–16
It is now possible to comprehensively study single nucleotide polymorphisms (SNPs) and haplotypes (the combination of SNP alleles observed on a chromosome) at a genetic locus.17 SNPs at a locus may be correlated (termed linkage disequilibrium [LD]),18 and within regions of high LD, there are a limited number of common haplotype patterns.19 A parsimonious subset of SNPs (termed tag SNPs) can be chosen that reliably predict these common haplotypes and therefore capture the underlying common genetic variation within regions of strong LD.20,21 By testing an efficient set of tag SNPs and (within regions of strong LD) their common haplotypes, potential causal variants may be evaluated for an association with phenotype, either directly, or indirectly via LD.22
Thus, we sought to (1) define the LD pattern for common SNP variants at the NOS3 locus in reference pedigrees and identify a set of SNPs that capture the underlying common variation and (2) genotype these SNPs in a large community-based sample and test the hypothesis that common variation in the NOS3 gene is associated with brachial artery FMD and reactive hyperemia.
|
Methods |
|---|
|
TopAbstractIntroductionMethodsResultsDiscussionReferences |
|---|
Study Participants
The design of the Framingham Heart Study (FHS) Offspring cohort has been described.23
Among the 3539 participants who attended the seventh examination of the Framingham Offspring Study, 2883 participants had adequate FMD measurements10 and 2031 had adequate hyperemic flow velocity flow measurements24 as described previously. Separately, during the sixth examination (1995–1998), 1809 unrelated individuals provided blood samples for DNA extraction.25 FHS randomly selected these unrelated individuals for DNA collection without regard to any phenotype feature. In the present study, everyone who had FMD or hyperemic flow velocity measured was eligible provided they had DNA available. This resulted in 1446 eligible for genotype-FMD and 1043 for genotype–reactive hyperemia association analyses. The Boston Medical Center institutional review board approved the study; participants gave written informed consent.
Noninvasive Measures of Endothelial Function
We a priori selected FMD and hyperemic flow velocity as phenotypes to be consistent with prior publications.10,26 We selected FMD to assess conduit artery endothelial function. Because hyperemic flow is the principal stimulus for FMD and is an indicator of microvascular function, we investigated hyperemic flow velocity. The methodology and reproducibility for measuring brachial artery FMD (expressed as percent increase in diameter) and mean hyperemic flow velocity (cm/s) have been described.10,24 Briefly, using a Toshiba SSH-140A ultrasound system and commercially available software (Brachial Analyzer version 3.2.3, Medical Imaging Applications), investigators blinded to participant clinical and genetic status determined brachial artery diameter at baseline and 1 minute after reactive hyperemia induced by 5-minute forearm cuff occlusion. Mean hyperemic flow velocity was analyzed with the use of semiautomated signal averaging (Cardiovascular Engineering).
Genotyping Methods
Genotyping was performed at the Broad Institute of the Massachusetts Institute of Technology/Harvard University by using matrix-assisted laser desorption ionization–time of flight mass spectrometry (Sequenom) to resolve allele-specific single-base extension products. The genotyping protocol has been described.27
SNP Selection and Genotyping in Reference Pedigrees
We surveyed common genetic variation in the region spanning NOS3 (chromosome 7; accession number NM_000603). We successfully genotyped common SNPs (minor allele frequency 
5%) with an average density of 1 SNP every 1.2 kilobases (kb), such that there were a minimum of 6 common SNPs per block of strong LD and gaps between blocks were <2 kb. We selected 139 markers from the public National Center for Biotechnology Information SNP map (http://www.ncbi.nlm.nih.gov/SNP) and the Celera variation database (http://www.celera.com). We genotyped these SNPs in a panel of 12 multigenerational pedigrees containing 93 individuals of western and northern European ancestry from Utah–Centre d’Etude du Polymorphisme Humain (CEPH) panel (Coriell Institute for Medical Research, Camden, NJ).27
We developed working assays for 33 polymorphic common SNPs (Table I in the online-only Data Supplement for CEPH minor allele frequencies) as defined by the following criteria: (1) at least 75% genotyping call success; (2) Hardy-Weinberg equilibrium P>0.01; and (3) mendelian transmission errors in CEPH pedigrees 
1.
LD Structure in Reference Pedigrees and Identification of Haplotype Blocks
For each pair of markers, D' (LD strength estimate)28 and logarithm of the odds score (LD significance estimate)29 were calculated. With the use of the "spine of LD" setting in Haploview software, haplotype blocks were defined on the basis of each end marker of a block having a D' >0.8, with all intervening pairwise marker comparisons, allowing for 1 comparison exception (version 2.03; http://www.broad.mit.edu/personal/jcbarret/haploview/).30
Tag SNP Selection
In the reference population, we used tag SNPs selection software (http://www-rcf.usc.edu/
stram), which calculates the squared correlation, Rh2, between haplotypes within LD blocks and haplotypes predicted with the use of only a subset of SNPs.21 We selected 11 of the 33 SNPs (3, 6, and 2 SNPs for blocks 1, 2, and 3, respectively) as tag SNPs; they represented the minimum subset of SNPs that was required to predict all common haplotypes (5% frequency) within each block, with Rh2 0.90. Using customized software (http://www.broad.mit.edu/mpg/tagger/) that assesses the ability of tag SNPs or their multimarker combinations to capture other SNPs, we separately evaluated the ability of tag SNPs to capture common variation at the locus, independent of block designation.
Genotyping in the FHS
In the Framingham sample we genotyped the 11 tag SNPs, the previously reported common missense variant, Glu298Asp (rs1799983),11,13–16 the previously studied T–786
C promoter variant (rs2070744), and 5 SNPs that were redundant in CEPH pedigrees (to help assess LD block structure similarity between CEPH and FHS samples; data not shown). Thus, 18 NOS3 SNPs were genotyped in FHS (Table I in the online-only Data Supplement). All SNPs were in Hardy-Weinberg equilibrium in the FHS sample (
2 test P>0.05).
Statistical Analysis
Using a general model of inheritance, we conducted multivariable linear regression analyses (SAS Proc GLM31) to test the null hypothesis that FMD and hyperemic flow velocity means did not differ by SNP genotype. In exploratory analyses, because of prior reports,11,12 we evaluated specific interaction terms in the multivariable-adjusted model, including sex for all 18 SNPs and smoking for the Glu298Asp variant. Because of previously reported sex differences11 and statistically significant sex-phenotype-genotype interactions observed in our data, we report sex-specific analyses.
Haplotype-based association analyses were conducted with the use of weighted-regression as implemented in the haplo.stats program.32,33 Haplotype frequencies in Framingham participants were estimated with the use of the expectation-maximization algorithm.34 A global score statistic tested whether trait differences existed among all haplotypes simultaneously; a haplotype-specific score statistic tested whether trait differences existed between a single haplotype versus all other haplotypes combined.
In both individual SNP and haplotype analyses, we adjusted for covariates previously associated with FMD or hyperemic flow response in our cohort,10,24 including mean arterial pressure, brachial artery pulse pressure, heart rate, body mass index, total/HDL cholesterol ratio, fasting glucose, diabetes, smoking within past 6 hours, prevalent cardiovascular disease, hormone replacement therapy, hypertension (systolic blood or diastolic pressure 140/90 mm Hg or antihypertensive medication use), lipid-lowering medication, and walk test (before or after FMD determination).
To account for multiple statistical testing, we constructed null data sets through bootstrap resampling; genotypes and phenotypes were sampled randomly with replacement. For each of the 18 SNPs, we ran 6 regression models (male, female, and sex-pooled, age-adjusted, and multivariable-adjusted models) 1000 times and evaluated the distribution of minimum probability values. We obtained an empirical probability value for each phenotype by comparing the minimum nominal probability value with the distribution of probability values from the null data sets.
|
Results |
|---|
|
TopAbstractIntroductionMethodsResultsDiscussionReferences |
|---|
LD Pattern at the NOS3 Locus in Reference Pedigrees
The LD structure at the NOS3 locus was defined by genotyping 33 common SNPs (minor allele frequency
5%) in CEPH pedigrees (Table I in the online-only Data Supplement). The mean interval between SNPs was 1.2 kb, and the SNP map spanned 39 kb in genomic distance (7 kb upstream, 23 kb NOS3 coding region, and 9 kb downstream). Three blocks of strong LD were evident: The blocks encompassed 14.2, 10.1, and 3.2 kb for blocks 1, 2, and 3, respectively (Figure). Both interblock distances were 0.5 kb. Within each block, there was limited haplotypic diversity. In blocks 1, 2, and 3, there were 4, 6, and 3 common haplotypes (frequency 5%) capturing 97%, 83%, and 98% of chromosomes, respectively (Figure).
|
Adequacy of SNP Selection
To examine the extent to which SNPs genotyped in our association study predicted unmeasured SNPs without regard to block structure, we calculated pairwise correlation (r2).35 In pairwise correlation between the 18 NOS3 SNPs genotyped in FHS and the remaining 15 unmeasured SNPs, we observed that 30 of the 33 common SNPs (91%) were predicted with an r2 0.70 (mean r2=0.90).
We compared the correlation between SNPs genotyped in FHS and NOS3 SNPs from the SeattleSNPs Program for Genomics Applications resequencing effort.36 SeattleSNPs identified 37 common NOS3 SNPs over a 25-kb genomic span in 23 individuals of European ancestry. Thirteen SNPs overlapped between the 37 identified by SeattleSNPs resequencing and the 18 SNPs genotyped in FHS. These 13 SNPs genotyped had strong correlation (r2 >0.70) with the majority (76% [28 of 37 SNPs]) of SeattleSNPs variants (mean r2=0.72). Hence, we infer that the 18 SNPs genotyped in our association study captured most of the common genetic variation at the NOS3 gene.
FHS Participants
We examined 1446 (49.5% women) middle-aged to elderly FHS participants (mean age, 62±9 years; Table 1) with NOS3 genotype and FMD (1043 had hyperemic flow velocity measurements). The mean FMD and hyperemic flow velocity were 3.36±3.08% and 53.8±21.9 cm/s in women and 2.32±2.38% and 47.1±20.4 cm/s in men, respectively. The age- and sex-adjusted correlation between FMD and hyperemic flow velocity was 0.41 (P<0.0001).
|
Associations of Single SNPs With FMD or Hyperemic Flow Velocity
In men, 4 sequence variants showed modest nominally significant (ie, P<0.05, not corrected for multiple testing) associations with FMD (Table 2). For example, individuals homozygous or heterozygous for the minor allele of SNP 10 (rs1800781) had higher FMD versus major allele homozygotes (3.37±0.72%, 2.66±0.18%, 2.21±0.10% for AA, AG, and GG genotypes, respectively; P=0.027).
|
In men, 5 sequence variants showed modest nominally significant associations with hyperemic flow velocity (Table 3). For example, individuals heterozygous for the minor allele of SNP 16 (rs3918174) had higher hyperemic flow velocity compared with major allele homozygotes (51.9±1.54 cm/s versus 45.9±0.88 cm/s for AG versus AA; general model P=0.002). SNP 16 explains 
2% of the variability in hyperemic flow velocity in men.
|
Notably, we did not observe any significant association between the previously reported T–786C promoter variant (SNP 9: rs2070744) or Glu298Asp (SNP 15: rs1799983) coding variant and either FMD or hyperemic flow velocity. In women, we did not observe significant association between any of the 18 NOS3 SNPs and either phenotype.
Effect Modification by Gender or Smoking
The difference in association between hyperemic flow velocity and genotype in the 2 sexes was nominally significant in an interaction model (interaction term P<0.05 for 4 SNPs: SNPs 10, 11, 16, 29). There was not a significant sexxSNP interaction for FMD for the 18 SNPs tested. Given a previous report,11 we tested but did not observe any significant interaction between Glu298Asp (SNP 15: rs1799983) genotype and smoking for either phenotype (P>0.05).
Accounting for Multiple Testing
For association analyses involving FMD, 18 SNPs, and 6 models for each SNP, our minimum nominal P was 0.008 (in men, age-adjusted model for SNP 10). With 1000 resampling procedure runs, this corresponds to an empirical P=0.34 for FMD. For hyperemic flow velocity, our minimum nominal P was 0.0021 (in men, multivariable-adjusted model for SNP 16), corresponding to an empirical P=0.14.
Associations of Haplotypes With FMD or Hyperemic Flow Velocity
In women, we did not observe any significant association between NOS3 haplotypes and FMD or hyperemic flow velocity. Among men, global haplotype association tests were not significant in either block 1 or block 2 (Tables 4 and 5). However, given that global haplotype tests may miss more subtle associations, we compared individual haplotypes versus all other haplotypes combined. Among men, haplotype 1D in block 1 (12% frequency) showed a nominally significant association with higher FMD and hyperemic flow velocity (minimum haplotype P=0.008; Table 4). The minor A allele of SNP 10 (rs1800781), present exclusively on haplotype 1D (Figure), was nominally significantly associated with both FMD and hyperemic flow velocity in men (Tables 2 and 3
); thus, the single SNP and haplotype findings represent the same result.
|
Among men, haplotype 2C in block 2 (8% frequency) showed nominal significance for association with both phenotypes (haplotype P=0.04 and 0.03, respectively; Table 5). Haplotype 2F (5% frequency) showed a nominally significant association with higher FMD (haplotype P=0.006). The minor G allele of SNP 16 (rs3918174) was present on haplotypes 2C and 2F.
|
Among men, the global haplotype association test was nominally significant in block 3 (global P=0.006 and 0.02 for FMD and hyperemic flow velocity, respectively; Table 6). Men with haplotype 3B (33% frequency) had higher FMD and greater hyperemic flow velocity compared with all other haplotypes combined (haplotype P=0.007 and 0.01, respectively). The minor T allele of SNP 29 (rs1065299) was present on haplotype 3B; hence, SNP 29 represents the same association as haplotype 3B.
|
|
Discussion |
|---|
|
TopAbstractIntroductionMethodsResultsDiscussionReferences |
|---|
Principal Findings
In reference pedigrees, we characterized LD structure across a 39-kb genomic segment spanning the NOS3 gene using a dense SNP map and found 3 blocks of SNPs in strong LD. A subset of 11 tag SNPs captured a majority of common variation at the locus. With 33 common SNPs in reference pedigrees and the genotyping of 18 NOS3 SNPs in our association study, we have extended prior work by comprehensively studying the sex-specific association of common NOS3 genetic variation with endothelial function measures in a large community-based sample.
In women, we found no evidence of an association between endothelial function and genetic variation at the NOS3 locus. In men, however, we identified modest associations with both specific haplotypes and SNPs residing on these haplotypes. The first involves variants in the promoter region (SNPs 6 and 7). The second involves 1 of 3 intronic SNPs that are near-perfect proxies for each other: SNPs 10, 11, 16. The third possible site is SNP 29, which is located in intron 24. Although these previously unreported associations were nominally significant, the empirical probability value exceeded 0.05 when the multiple statistical tests performed for this locus were taken into account. Therefore, we view our findings as hypothesis generating, underscoring the need for additional studies to test our observations.
Comparison With Prior Studies
Prior association studies examining the relation between endothelial function and NOS3 variants have focused on Glu298Asp, T–786C (SNP 9: rs2070744 in our study; Table I in the online-only Data Supplement), and a variable number tandem repeat (VNTR) in intron 4.11–16 In general, those studies were modest in size (n <250) and yielded mixed results.11–16 In particular, the Glu298Asp variant has been a focus because of initial suggestions that the mutation led to differential susceptibility to intracellular cleavage37; however, a more recent in vitro study has challenged this assertion.38 Leeson et al11 evaluated 248 young participants and reported that the Glu298Asp genotype did not relate to FMD in the overall sample, a finding that is consistent with our results.
One study observed a relation between the T–786C promoter variant and forearm blood flow responses during acetylcholine infusion in hypertensive patients,13 but other studies found no relation between this polymorphism and FMD.14,15 Cattaruzza et al39 evaluated the response of endothelial cells (isolated from umbilical cord vein samples) to 24 to 36 hours of laminar shear stress and reported that in samples with the CC genotype, NOS3 mRNA and protein levels did not increase in response to shear stress. We observed that the T–786C variant (rs20720744, SNP 9) was not associated with FMD or hyperemic flow velocity, although the nominal P was of borderline significance (P=0.06). The intron 4 VNTR polymorphism has been related to lower FMD in diabetic patients12 but not in a cohort of black subjects.14 We did not directly test the intron 4 VNTR; however, this polymorphism has been previously shown to be in LD with the T–786C variant.40
In summary, differences in phenotype and variants studied limit direct comparison of our findings to prior work. However, the cumulative evidence argues against an effect of Glu298Asp coding variant on FMD or hyperemic flow. Similar to prior work, our study leaves open the possibility that endothelial function might relate to T–786C, intron 4 VNTR, or a SNP in LD.
SNPs and Haplotypes
In addition to the much larger sample size, the present study differed from prior work by including haplotype analyses. Testing haplotypes in association studies may offer several advantages.41 For example, specific multimarker combinations may be better predictors of unmeasured variants than an individual SNP. In addition, SNP alleles interacting in cis may produce a greater effect than any individual variant. The associations of 5-lipoxygenase activating protein haplotypes and myocardial infarction42 or ß2-adrenergic receptor haplotypes and bronchodilator responsiveness43 provide examples of haplotypes apparently providing information beyond testing specific SNPs. However, in the present study the haplotypes and SNP analyses provided similar insights.
Role of Multiple Testing
To facilitate discovery, it has been suggested that candidate gene association studies transition from testing single variants to gene-based approaches in which all common variation within a gene is considered jointly.17 We undertook such an approach and identified several modest nominal associations. However, considering multiple variants at a locus may increase type I error.44 Methods to determine whether putative associations identified are bona fide or statistical fluctuations include stringent statistical criteria to correct for multiple testing and replication.45
We constructed bootstrapped null-distribution data sets to guide the interpretation of our results. After consideration of the number of SNPs and analytical models studied, our strongest result was observed by chance >5% of the time for both phenotypes. Thus, despite observing a minimum nominal P=0.002, declaration of a bona fide association would be premature. Instead, we regard our putative associations as intriguing hypotheses (specifically SNPs 6, 7, 10, 11, 16, and 29) for independent replication studies. Replication would require determination of vascular function and NOS3 genotype determination, a similar-sized or larger sample, and significant association findings with the same alleles with the same direction of effect.
Study Limitations and Strengths
Our study has several limitations. There remains the possibility that multiple rare NOS3 mutations influence endothelial function measures. The allelic spectrum for complex traits such as FMD and hyperemic flow velocity is unknown. Because most human heterozygosity is accounted for by common variants and susceptibility alleles for common disease (and related subclinical disease measures) may not be under strong negative selection, common variants may contribute to common disease traits.46 Our approach is most suitable for the discovery of such variants.
Second, the optimal methodology to select a nonredundant tag SNP set at a locus is still evolving.21,47,48 We tested the performance of our set using several approaches; the tag SNPs we selected captured the majority of common variation at the NOS3 locus. Third, our cohort was middle-aged to elderly and white, potentially limiting generalizability to individuals who are younger or of other ethnicities/races.
Fourth, insufficient sample size or random measurement error may have limited our power to detect phenotype-genotype association. Although accepted as an in vivo measure of endothelial function, measurement of brachial artery FMD is subject to a number of well-recognized pitfalls and limitations.49 In particular, using below-elbow blood pressure cuff position may have led to low FMD values and less biological variability.50 To address this issue, we performed the following power calculations on the basis of our data. We set type I error at 
=0.0028 (ie, =0.05/18 due to 18 SNPs), fixed sex-specific sample sizes at 700 for FMD (600 for hyperemic flow velocity), and assumed that covariates explained 15% of interindividual variation in FMD (35% for hyperemic flow velocity). Given these conditions, we had 80% power to detect any SNP that is expected to account for 2.1% of interindividual variability in FMD (for hyperemic flow velocity, 1.9%). To place our power to detect phenotype-genotype associations in context, systolic blood pressure and age explained 3.2% to 12.4% of variance in FMD and 6.0% to 26.3% for hyperemic flow velocity, whereas cardiovascular disease explained 1.8% of variance in hyperemic flow velocity among men. If the true effects of SNPs on phenotypes were smaller, we would have had lower power and we may have failed to detect true associations.
Finally, by controlling the false-positive rate through empirical methods, we may have increased our false-negative rate.44 Rigid adherence to an empirical P<0.05 significance threshold across a study could be overly conservative and obscure some true positive associations.51
Study strengths distinguishing the present investigation include the characterization of NOS3 genetic variation in reference pedigrees, genotyping of a comprehensive set of common variants in a large community-based cohort with routinely ascertained phenotypes and covariates, sex-specific analyses, the inclusion of both single variant and haplotype association, and multiple statistical testing assessments.
Conclusion
We created a dense SNP map at the NOS3 locus and defined a parsimonious set of SNPs that captured common genetic variation. Overall, we did not observe statistically significant relations of NOS3 sequence variants with endothelial function measures in men or women in our sample after accounting for multiple testing. The nominal associations of select NOS3 variants with endothelial function in men may provide hypotheses for further testing in other samples or in different experimental designs. In addition, our results provide the tools necessary to comprehensively test common genetic variation at NOS3 for association with other phenotypes potentially influenced by NO such as hypertension and incident coronary heart disease.
|
Acknowledgments |
|---|
This study was supported by National Heart, Lung, and Blood Institute grant N01-HC-25195; CardioGenomics Programs for Genomic Applications (HL66582); HL060040; HL070100; HL071039, K24HL004334 (Dr Vasan); and the Donald W. Reynolds Foundation. Dr Kathiresan is supported by ACC Foundation/Merck Adult Cardiology Research Fellowship Award and GlaxoSmithKline Research and Education Foundation for Cardiovascular Disease Young Investigator Award.
|
Footnotes |
|---|
Guest Editor for this article was Mark A. Creager, MD.
The online-only Data Supplement, which contains 3 supplemental tables, can be found at http://circ.ahajournals.org/cgi/content/full/CIRCULATIONAHA.105.544619/DC1.
|
References |
|---|
|
TopAbstractIntroductionMethodsResultsDiscussionReferences |
|---|
1. Palmer RM, Ashton DS, Moncada S. Vascular endothelial cells synthesize nitric oxide from L-arginine. Nature. 1988; 333: 664–666.[CrossRef][Medline] [Order article via Infotrieve]
2. Corretti MC, Anderson TJ, Benjamin EJ, Celermajer D, Charbonneau F, Creager MA, Deanfield J, Drexler H, Gerhard-Herman M, Herrington D, Vallance P, Vita J, Vogel R. Guidelines for the ultrasound assessment of endothelial-dependent flow-mediated vasodilation of the brachial artery: a report of the International Brachial Artery Reactivity Task Force. J Am Coll Cardiol. 2002; 39: 257–265.
3. Dakak N, Husain S, Mulcahy D, Andrews NP, Panza JA, Waclawiw M, Schenke W, Quyyumi AA. Contribution of nitric oxide to reactive hyperemia: impact of endothelial dysfunction. Hypertension. 1998; 32: 9–15.
4. Widlansky ME, Gokce N, Keaney JF Jr, Vita JA. The clinical implications of endothelial dysfunction. J Am Coll Cardiol. 2003; 42: 1149–1160.
5. Targonski PV, Bonetti PO, Pumper GM, Higano ST, Holmes DR Jr, Lerman A. Coronary endothelial dysfunction is associated with an increased risk of cerebrovascular events. Circulation. 2003; 107: 2805–2809.
6. Ganz P, Vita JA. Testing endothelial vasomotor function: nitric oxide, a multipotent molecule. Circulation. 2003; 108: 2049–2053.
7. Huang PL, Huang Z, Mashimo H, Bloch KD, Moskowitz MA, Bevan JA, Fishman MC. Hypertension in mice lacking the gene for endothelial nitric oxide synthase. Nature. 1995; 377: 239–242.[CrossRef][Medline] [Order article via Infotrieve]
8. Celermajer DS, Sorensen KE, Bull C, Robinson J, Deanfield JE. Endothelium-dependent dilation in the systemic arteries of asymptomatic subjects relates to coronary risk factors and their interaction. J Am Coll Cardiol. 1994; 24: 1468–1474.[Abstract]
9. Gaeta G, De Michele M, Cuomo S, Guarini P, Foglia MC, Bond MG, Trevisan M. Arterial abnormalities in the offspring of patients with premature myocardial infarction. N Engl J Med. 2000; 343: 840–846.
10. Benjamin EJ, Larson MG, Kupka MJ, Mitchell GF, Vasan RS, Keaney JF Jr, Lehman BT, Fan S, Osypiuk E, Vita JA. Clinical correlates and heritability of flow-mediated dilation in the community: the Framingham Heart Study. Circulation. 2004; 110: 3604–3609.
11. Leeson CP, Hingorani AD, Mullen MJ, Jeerooburkhan N, Kattenhorn M, Cole TJ, Muller DP, Lucas A, Humphries SE, Deanfield JE. Glu298Asp endothelial nitric oxide synthase gene polymorphism interacts with environmental and dietary factors to influence endothelial function. Circ Res. 2002; 90: 1153–1158.
12. Komatsu M, Kawagishi T, Emoto M, Shoji T, Yamada A, Sato K, Hosoi M, Nishizawa Y. ecNOS gene polymorphism is associated with endothelium-dependent vasodilation in type 2 diabetes. Am J Physiol. 2002; 283: 557–561.
13. Rossi GP, Taddei S, Virdis A, Cavallin M, Ghiadoni L, Favilla S, Versari D, Sudano I, Pessina AC, Salvetti A. The T-786C and Glu298Asp polymorphisms of the endothelial nitric oxide gene affect the forearm blood flow responses of Caucasian hypertensive patients. J Am Coll Cardiol. 2003; 41: 938–945.
14. Li R, Lyn D, Lapu-Bula R, Oduwole A, Igho-Pemu P, Lankford B, Morgan J, Nkemdechi S, Liu G, Pack C, Silvestrov N, von Deutsch DA, Song Q, Abukhalaf IK, Ofili E. Relation of endothelial nitric oxide synthase gene to plasma nitric oxide level, endothelial function, and blood pressure in African Americans. Am J Hypertens. 2004; 17: 560–567.[CrossRef][Medline] [Order article via Infotrieve]
15. Paradossi U, Ciofini E, Clerico A, Botto N, Biagini A, Colombo MG. Endothelial function and carotid intima-media thickness in young healthy subjects among endothelial nitric oxide synthase Glu298Asp and T-786C polymorphisms. Stroke. 2004; 35: 1305–1309.
16. Naber CK, Baumgart D, Altmann C, Siffert W, Erbel R, Heusch G. eNOS 894T allele and coronary blood flow at rest and during adenosine-induced hyperemia. Am J Physiol. 2001; 281: H1908–H1912.
17. Neale BM, Sham PC. The future of association studies: gene-based analysis and replication. Am J Hum Genet. 2004; 75: 353–362.[CrossRef][Medline] [Order article via Infotrieve]
18. Kruglyak L. Prospects for whole-genome linkage disequilibrium mapping of common disease genes. Nat Genet. 1999; 22: 139–144.[CrossRef][Medline] [Order article via Infotrieve]
19. Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES. High-resolution haplotype structure in the human genome. Nat Genet. 2001; 29: 229–232.[CrossRef][Medline] [Order article via Infotrieve]
20. Johnson GC, Esposito L, Barratt BJ, Smith AN, Heward J, Di Genova G, Ueda H, Cordell HJ, Eaves IA, Dudbridge F, Twells RC, Payne F, Hughes W, Nutland S, Stevens H, Carr P, Tuomilehto-Wolf E, Tuomilehto J, Gough SC, Clayton DG, Todd JA. Haplotype tagging for the identification of common disease genes. Nat Genet. 2001; 29: 233–237.[CrossRef][Medline] [Order article via Infotrieve]
21. Stram DO, Haiman CA, Hirschhorn JN, Altshuler D, Kolonel LN, Henderson BE, Pike MC. Choosing haplotype-tagging SNPS based on unphased genotype data using a preliminary sample of unrelated subjects with an example from the Multiethnic Cohort Study. Hum Hered. 2003; 55: 27–36.[CrossRef][Medline] [Order article via Infotrieve]
22. The International HapMap Project. Nature. 2003; 426: 789–796.[CrossRef][Medline] [Order article via Infotrieve]
23. Kannel WB, Feinleib M, McNamara PM, Garrison RJ, Castelli WP. An investigation of coronary heart disease in families: the Framingham Offspring Study. Am J Epidemiol. 1979; 110: 281–290.
24. Mitchell GF, Parise H, Vita JA, Larson MG, Warner E, Keaney JF Jr, Keyes MJ, Levy D, Vasan RS, Benjamin EJ. Local shear stress and brachial artery flow-mediated dilation: the Framingham Heart Study. Hypertension. 2004; 44: 134–139.
25. Shearman AM, Cupples LA, Demissie S, Peter I, Schmid CH, Karas RH, Mendelsohn ME, Housman DE, Levy D. Association between estrogen receptor alpha gene variation and cardiovascular disease. JAMA. 2003; 290: 2263–2270.
26. Vita JA, Keaney JF Jr, Larson MG, Keyes MJ, Massaro JM, Lipinska I, Lehman BT, Fan S, Osypiuk E, Wilson PW, Vasan RS, Mitchell GF, Benjamin EJ. Brachial artery vasodilator function and systemic inflammation in the Framingham Offspring Study. Circulation. 2004; 110: 3604–3609.
27. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D. The structure of haplotype blocks in the human genome. Science. 2002; 296: 2225–2229.
28. Lewontin RC. The interaction of selection and linkage, II: optimum models. Genetics. 1964; 50: 757–782.
29. Pritchard JK, Przeworski M. Linkage disequilibrium in humans: models and data. Am J Hum Genet. 2001; 69: 1–14.[CrossRef][Medline] [Order article via Infotrieve]
30. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21: 263–265.
31. SAS/STAT User’s Guide, Version 8.1. Cary, NC: SAS Institute, Inc; 2000.
32. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA. Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002; 70: 425–434.[CrossRef][Medline] [Order article via Infotrieve]
33. Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, Ehm MG. Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered. 2002; 53: 79–91.[Medline] [Order article via Infotrieve]
34. Excoffier L, Slatkin M. Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995; 12: 921–927.[Abstract]
35. Devlin B, Risch N. A comparison of linkage disequilibrium measures for fine-scale mapping. Genomics. 1995; 29: 311–322.[CrossRef][Medline] [Order article via Infotrieve]
36. SeattleSNPs. NHLBI Program for Genomics Applications, UW-FHCRC. Available at: http://pga.gs.washington.edu. Accessed December 23, 2004.
37. Tesauro M, Thompson WC, Rogliani P, Qi L, Chaudhary PP, Moss J. Intracellular processing of endothelial nitric oxide synthase isoforms associated with differences in severity of cardiopulmonary diseases: cleavage of proteins with aspartate vs. glutamate at position 298. Proc Natl Acad Sci U S A. 2000; 97: 2832–2835.
38. Fairchild TA, Fulton D, Fontana JT, Gratton JP, McCabe TJ, Sessa WC. Acidic hydrolysis as a mechanism for the cleavage of the Glu(298)Asp variant of human endothelial nitric-oxide synthase. J Biol Chem. 2001; 276: 26674–26679.
39. Cattaruzza M, Guzik TJ, Slodowski W, Pelvan A, Becker J, Halle M, Buchwald AB, Channon KM, Hecker M. Shear stress insensitivity of endothelial nitric oxide synthase expression as a genetic risk factor for coronary heart disease. Circ Res. 2004; 95: 841–847.
40. Tanus-Santos JE, Desai M, Flockhart DA. Effects of ethnicity on the distribution of clinically relevant endothelial nitric oxide variants. Pharmacogenetics. 2001; 11: 719–725.[CrossRef][Medline] [Order article via Infotrieve]
41. Hoehe MR. Haplotypes and the systematic analysis of genetic variation in genes and genomes. Pharmacogenomics. 2003; 4: 547–570.[CrossRef][Medline] [Order article via Infotrieve]
42. Helgadottir A, Manolescu A, Thorleifsson G, Gretarsdottir S, Jonsdottir H, Thorsteinsdottir U, Samani NJ, Gudmundsson G, Grant SF, Thorgeirsson G, Sveinbjornsdottir S, Valdimarsson EM, Matthiasson SE, Johannsson H, Gudmundsdottir O, Gurney ME, Sainz J, Thorhallsdottir M, Andresdottir M, Frigge ML, Topol EJ, Kong A, Gudnason V, Hakonarson H, Gulcher JR, Stefansson K. The gene encoding 5-lipoxygenase activating protein confers risk of myocardial infarction and stroke. Nat Genet. 2004; 36: 233–239.[CrossRef][Medline] [Order article via Infotrieve]
43. Drysdale CM, McGraw DW, Stack CB, Stephens JC, Judson RS, Nandabalan K, Arnold K, Ruano G, Liggett SB. Complex promoter and coding region beta 2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc Natl Acad Sci U S A. 2000; 97: 10483–10488.
44. Wacholder S, Chanock S, Garcia-Closas M, El Ghormli L, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst. 2004; 96: 434–442.
45. Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN. Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003; 33: 177–182.[CrossRef][Medline] [Order article via Infotrieve]
46. Reich DE, Lander ES. On the allelic spectrum of human disease. Trends Genet. 2001; 17: 502–510.[CrossRef][Medline] [Order article via Infotrieve]
47. Byng MC, Whittaker JC, Cuthbert AP, Mathew CG, Lewis CM. SNP subset selection for genetic association studies. Ann Hum Genet. 2003; 67: 543–556.[CrossRef][Medline] [Order article via Infotrieve]
48. Chapman JM, Cooper JD, Todd JA, Clayton DG. Detecting disease associations due to linkage disequilibrium using haplotype tags: a class of tests and the determinants of statistical power. Hum Hered. 2003; 56: 18–31.[CrossRef][Medline] [Order article via Infotrieve]
49. Deanfield J, Donald A, Ferri C, Giannattasio C, Halcox J, Halligan S, Lerman A, Mancia G, Oliver JJ, Pessina AC, Rizzoni D, Rossi GP, Salvetti A, Schiffrin EL, Taddei S, Webb DJ. Endothelial function and dysfunction, part I: methodological issues for assessment in the different vascular beds: a statement by the Working Group on Endothelin and Endothelial Factors of the European Society of Hypertension. J Hypertens. 2005; 23: 7–17.[CrossRef][Medline] [Order article via Infotrieve]
50. Frick M, Alber HF, Weidinger F. Endothelial function in a large community. Circulation. 2004; 110: e24.
51. Thomas DC, Clayton DG. Betting odds and genetic associations. J Natl Cancer Inst. 2004; 96: 421–423.
This article has been cited by other articles:
 |
|
 |
|
  L. Moro, C. Pedone, S. Scarlata, V. Malafarina, F. Fimognari, and R. Antonelli-Incalzi Endothelial Dysfunction in Chronic Obstructive Pulmonary Disease Angiology, July 1, 2008; 59(3): 357 - 364. [Abstract] [PDF]
|
|
|
|
 |
|
 M.-X. Zhang, C. Zhang, Y. H. Shen, J. Wang, X. N. Li, Y. Zhang, J. Coselli, and X. L. Wang Biogenesis of Short Intronic Repeat 27-Nucleotide Small RNA from Endothelial Nitric-oxide Synthase Gene J. Biol. Chem., May 23, 2008; 283(21): 14685 - 14693. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 G. F. Mitchell, C.-Y. Guo, S. Kathiresan, R. S. Vasan, M. G. Larson, J. A. Vita, M. J. Keyes, M. Vyas, C. Newton-Cheh, S. L. Musone, et al. Vascular Stiffness and Genetic Variation at the Endothelial Nitric Oxide Synthase Locus: The Framingham Heart Study Hypertension, June 1, 2007; 49(6): 1285 - 1290. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||