PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Circulation. Author manuscript; available in PMC 2010 February 24.
Published in final edited form as:
PMCID: PMC2782762
NIHMSID: NIHMS126507

Genetic Variations in NOS1AP are Associated with Sudden Cardiac Death in U.S. White Community Based Populations

WH Linda Kao, PhD, MHS,1,3,* Dan E Arking, PhD,2,3,* Wendy Post, MD, MS,1,3,* Thomas D Rea, MD, MPH,4 Nona Sotoodehnia, MD, MPH,4 Ronald J Prineas, MD, PhD,5 Bryan Bishe, BS,2 Betty Q Doan, PhD, MPH,2 Eric Boerwinkle, PhD,6 Bruce M Psaty, MD, PhD,4,7 Gordon Tomaselli, MD,3 Josef Coresh, MD, PhD, MHS,1,3,8 David S Siscovick, MD,4,7 Eduardo Marban, MD, PhD,8 Peter M Spooner, PhD,3 Greg L Burke, MD, MSc,5 and Aravinda Chakravarti, PhD2,3,9,10

Abstract

Background

The electrocardiographic QT interval is associated with risk of sudden cardiac death (SCD). A previous genome-wide association study demonstrated that allelic variants (rs10494366 and rs4657139) in NOS1AP, which encodes a carboxy-terminal PDZ ligand of neuronal nitric oxide synthase, are associated with the QT interval in white adults. The present analysis was conducted to validate the association between NOS1AP variants and the QT interval and to examine the association with SCD in a combined population of 19,295 black and white adults from the Atherosclerosis Risk in Communities (ARIC) study and the Cardiovascular Health Study (CHS).

Methods and Results

We examined 19 tagging SNPs in the genomic blocks containing rs10494366 and rs4657139 in NOS1AP. SCD was defined as a sudden pulseless condition of cardiac origin in a previously stable individual. General linear models and Cox proportional hazards regression models were used. Multiple SNPs in NOS1AP, including rs10494366, rs4657139, and rs16847548 were significantly associated with adjusted QT interval in whites (P<0.0001). In whites, after adjusting for age, sex, and study, the relative hazard (RH) of SCD associated with each C allele at rs16847548 was 1.31 (95% CI: 1.10 to 1.56, P=0.002), assuming an additive model. In addition, a downstream neighboring SNP, rs12567209, not correlated with rs16847548 or QT interval, was also independently associated with SCD in whites (RH = 0.57, 95%CI: 0.39, 0.83; P = 0.003). Adjustment for QT interval and CHD risk factors attenuated, but did not eliminate, the association between rs16847548 and SCD, and such adjustment had no effect on the association between rs12567209 and SCD. No significant associations between tagging SNPs in NOS1AP and either QT interval or SCD were observed in blacks.

Conclusions

In a combined analysis of two population-based prospective cohort studies, sequence variations in NOS1AP were associated with baseline QT interval and the risk of SCD in white U.S. adults.

Keywords: death, sudden, QT interval, genetics, epidemiology

Sudden cardiac death (SCD) and cardiac arrhythmias remain a daunting public health problem. It is estimated that there are between 250,000 and 400,000 sudden cardiac deaths in the United States each year1, 2. Nearly half of all coronary heart disease (CHD) deaths are sudden and approximately 1/3 of these deaths are the first clinical manifestation of disease3. Thus it is important to identify risk factors, both genetic and environmental, for SCD.

Previous studies have identified family history of SCD as a powerful risk factor for SCD that is independent of traditional risk factors for CHD or a family history of myocardial infarction46. Moreover, a number of genes have been linked to rare, heritable arrhythmias that predispose to SCD. However, the genetic factors underlying SCD in the general population are largely unknown.

SCD is a complex phenotype of heterogeneous etiology with multiple factors contributing to risk that can be broadly classified into three categories: 1) atherosclerosis and thrombosis; 2) electrogenesis and propagation; and 3) initiating influences and triggers7. Indeed, SCD is a multi-factorial disorder involving the interaction of multiple genes in conjunction with environmental influences. A number of pathways modulate the electrophysiology of the heart and have been associated with enhanced risk of SCD. Altered ventricular repolarization reflected in abnormalities of the electrocardiographic QT interval is an intermediate phenotype that is not only associated with an increased risk of SCD in both the presence and absence of structural heart disease8 but also heritable within families9 and in population-based studies10. Using a genome-wide association study (GWAS) we have previously demonstrated that allelic variants in NOS1AP (nitric oxide synthase 1 adaptor protein), which encodes a cytosolic ligand of neuronal nitric oxide synthase (nNOS), are associated with altered QT intervals in white adults11. This association has subsequently been replicated in additional white populations1217. The objectives of the present study were to validate the association of NOS1AP variants with QT interval prolongation in a large U.S. population-based cohort of white and black adults and, more importantly, to establish the association between NOS1AP variants and the risk of SCD in these community-based individuals. We hypothesized that at least one of the NOS1AP SNPs examined would be associated with the QT interval; moreover, at least one such SNP would also be associated with the risk of SCD.

Methods

The Atherosclerosis Risk in Communities (ARIC) study and the Cardiovascular Health Study (CHS)

The ARIC study and CHS are both population-based prospective cohort studies of cardiovascular disease. The ARIC Study includes 15,792 persons aged 45–64 years at baseline (1987–89), randomly chosen from four US communities18. ARIC cohort members completed four clinic examinations, conducted approximately three years apart between 1987 and 1998. CHS includes 5,888 participants > 65 years of age identified from four U.S. communities using Medicare eligibility lists. The original cohort included 5201 participants recruited in 1989–1990 and 687 additional subjects were recruited in 1992–1993 to enhance the racial/ethnic diversity of the cohort19.

Clinic examinations for both ARIC and CHS participants included assessment of cardiovascular risk factors, self-reported medical family history, employment and educational status, diet, physical activity, comorbidities, and clinical and laboratory measurements. ARIC and CHS participants were contacted annually by telephone for identification of all hospitalizations and deaths, and lists of discharges from local hospitals were scanned for events. Deaths were identified from death certificates, and potential out-of-hospital fatal coronary heart disease events were investigated by interviewing one or more next of kin and by the completion of a questionnaire by the patient’s physician. ARIC and CHS staff abstract discharge diagnoses on all hospitalizations, as well as conduct standardized committee review of all CHD, stroke and cardiovascular death events18, 20, 21. Comprehensive data were gathered on cardiovascular events and deaths from hospital records, interviews with physicians, next of kin and/or witnesses, death certificates, and autopsy reports22. In addition to similar study protocols between ARIC and CHS, extensive review of the data definitions and study sources was carried out with review of files by an independent set of investigators so that only comparable clinical variables were included in the current analysis of the combined cohorts.

The following exclusion criteria, based on missing exposure or outcome data, were applied to obtain the final sample for the present analysis: 1) self-reported race/ethnicity other than black or white (48 in ARIC, 39 in CHS), 2) samples not genotyped due to lack of DNA or consent for genetic research (103 in ARIC; 432 in CHS), 3) samples with < 75% of genotypes called (1052 in ARIC, 305 in CHS), 4) missing electrocardiograms (not performed or not transmitted) or poor quality data (missing leads or artifacts) for either QT or heart rate (189 in ARIC, 124 in CHS), and 5) unconfirmed SCD (82 in ARIC, 2 in CHS). After these exclusions, 14,309 of 15,783 ARIC participants (91%) and 4,986 of 5,888 CHS participants (85%) were included in the present analysis for a combined sample size of 19,295 individuals.

Assessment of Sudden Cardiac Death

Each parent study classified all cases of fatal CHD according to standard protocols. To identify cases of SCD in ARIC and CHS for the present study, all cases of fatal CHD that occurred by July 31, 2002 in CHS and December 31, 2002 in ARIC were reviewed and adjudicated by a committee of physicians. SCD was operationally defined as a sudden pulseless condition from a cardiac origin in a previously stable individual. After review of data available from death certificates, informant interviews, physician questionnaires, coroner reports, and hospital discharge summaries, the reviewers classified each CHD death as definite sudden arrhythmic death, possible sudden arrhythmic death, definite non-sudden death, or unclassifiable. We a priori sought to exclude cases with non-arrhythmic characteristics including those with evidence of progressive hypotension or advanced congestive heart failure prior to death. We also excluded those cases with advanced dementia or terminal illness such as end stage cancer or liver disease. Each event was independently adjudicated by two investigators. If disagreement existed between the first two reviewers, a third investigator independently reviewed the event to provide final classification. As part of event review, information was systematically abstracted regarding duration of symptoms, whether the event was witnessed, other circumstances of the event, and medical co-morbidities of the patient in order to help classify whether the subject had experienced SCD. Those classified as “definite sudden arrhythmic death” were either confirmed by evidence of “instantaneous death” or in the case of unwitnessed deaths, there was descriptive information regarding the position of the body that indicated a sudden event had occurred. All suspected SCD, defined as a sudden pulseless condition from a cardiac origin in a previously stable individual, that we could not classify as “definite” were classified as “possible SCD”. Cases were identified as either in or out of hospital deaths. The primary outcome of SCD described in the present study combines both definite and possible sudden arrhythmic death. For the present analysis, participants were censored at time of loss to follow up or death if the cause of death was other than SCD. The administrative censoring date was July 31, 2002 for CHS and December 31, 2002 for ARIC, based on the study’s adjudication schedules.

Assessment of QT interval at the baseline examinations of each study

At the baseline visit of the ARIC study, participants were asked not to smoke or ingest caffeine for at least 1 hour prior to the electrocardiogram. After resting for 5–10 minutes while the electrodes were being placed, a standard supine 12-lead electrocardiogram and a 2-minute paper recording of a three-lead (leads V1, II, and V5) rhythm strip were made. The ECGs were digitally recorded, and identical methods (MAC personal computer, Marquette Electronics, Milwaukee, Wisconsin) were used in all clinical centers. A similar protocol was used at the baseline visit of CHS. MAC PC-DT ECG acquisition units (Marquette Electronics, Inc., Milwaukee, WI) were used to record a 10-second 12-lead simultaneous ECG at a sample rate of 250 per second per lead. All tracings from the baseline visits of both CHS and ARIC were transmitted over analogue phone lines to a central ECG Reading Center in Edmonton, Alberta, Canada for analysis. The QT interval from the digital 12-lead ECG was determined using the Novacode ECG measurement and classification program23.

Assessment of Covariates

Both the ARIC study and CHS have extensive data on behavioral, clinical, and serologic factors relevant to selected cardiovascular phenotypes and outcomes. At each visit, demographic, anthropometric, and cardiovascular risk factor data were collected. Data from the baseline visits of both ARIC and CHS were used for the present analyses. Participants described themselves as white or black in response to an interviewer-administered questionnaire, which also contained questions on highest education attained, smoking status, and marital status. Collection of fasting blood samples and processing for total cholesterol followed standard study protocols19, 24. Hypertension was defined as systolic blood pressure (SBP) ≥ 140 mmHg, diastolic blood pressure (DBP) ≥ 90 mmHg, or use of antihypertensive medications. Diabetes was defined as fasting glucose ≥ 126 mg/dL, non-fasting glucose ≥ 200 mg/dL, or history/treatment of diabetes. History of myocardial infarction (MI) at the baseline examination was defined by either a self-reported heart attack requiring hospitalization, self-reported history of physician-diagnosed MI, or a history of MI identified on the baseline electrocardiogram, which was verified in CHS25. A positive self-reported family history of CVD was defined as having at least one parent with CHD or stroke in ARIC and having at least one full sibling with heart attack or stroke in CHS.

SNP Selection and Genotyping

Since over 80% of the present study subjects were whites, SNPs were selected to tag the linkage disequilibrium (LD) block containing rs10494366 and rs4657139 (the most significant SNPs from fine mapping in our initial study in whites) and its neighboring LD blocks in the thirty trio samples of U.S. residents with northern and western European ancestry (CEU population) used in the HapMap Project26, 27. Twenty-one SNPs were selected using the computer program Tagger with criteria of r2 >0.65 and minor allele frequency (MAF) >0.05 in CEU (data from 19 SNPs passed QC tests). SNPs from coding regions were not specifically interrogated.

Genotyping was performed using TaqMan assays (Applied Biosystems) in conjunction with the BioTrove OpenArray SNP genotyping platform28. Data completeness (reported as the percentage of samples that have a genotype call at a given SNP) across all SNPs and samples was 93.7%, and the accuracy of genotyping, determined by comparison to concordance calls generated for 58 samples genotyped multiple times (range 2–19 times, median 6, resulting in ~350 comparison per SNP), was 99.1%. Two SNPs were dropped from further analysis due to significantly lower accuracy (96.5% and 96.1%). Individual samples with <75% complete data were also removed from further analysis, as low data completeness was strongly associated with higher genotyping error rates. The overall accuracy and data completeness were 99.4% and 97.8%, respectively, after removing poor quality DNA samples and SNPs.

Statistical Analysis

Differences in baseline characteristics by subsequent SCD status were assessed using t-tests and χ2 tests within each self-reported race/ethnicity. All analyses were stratified by self-reported ethnicity. Because the study protocols were similar across ARIC and CHS, the number of SCD events was modest within each study, and initial genotype-phenotype analyses indicated that the associations did not differ between the two studies, all analyses were pooled across the two studies to increase statistical power. The rare allele of each SNP in whites was designated as the minor allele. Deviations from Hardy-Weinberg proportions were assessed using the chi-squared goodness of fit test within each ethnicity group.

To analyze the QT-interval, a Z score for QT interval was created for each individual by standardizing an individual’s QT-interval to the study-specific population mean and standard error. A generalized linear model was then used to assess the association between SNPs and these Z scores representing QT intervals assuming the additive genetic model while adjusting for age, sex, and heart rate. All results presented in the present paper were generated using the study-specific age-, sex-, and heart-rate-adjusted QT intervals. Parallel analyses were also performed using Bazett’s heart rate-corrected QT duration (QTc)29, and similar results and inferences were obtained (data not shown). Analysis of the QT interval included both individuals who later experienced SCD and those who did not. Since 19 SNPs were tested, we used a Bonferroni-corrected-alpha of 0.003 (0.05/19) as the threshold for statistical significance in the age, sex, and study-adjusted model.

For SCD risk, single SNP genotype--based analyses were performed. Cumulative incidence of SCD in the presence of competing events (deaths due to other causes) was estimated overall and by SNP genotype30. To estimate the relative hazards and the significance of the association between each SNP genotype and SCD risk while adjusting for covariates, Cox proportional hazards models were constructed, and a Bonferroni-corrected-alpha of 0.003 was used to declare statistical significance in the age, sex, and study adjusted analysis. An additive model was assumed for each SNP. For significant SNPs, a model assuming 3 genotypic risks was also constructed to confirm the use of the additive model. The proportional hazards assumption was checked with Schoenfeld’s residual31.

As exploratory analyses, the role of genotypic effects across various high-risk subgroups was examined both with stratified analyses and by fitting interaction terms into the regression models. All analyses were performed with either SAS (version 9.0) or STATA (version 9.2).

Statement of Responsibility

The authors had full access to and take full responsibility for the integrity of the data. All authors have read agree to the manuscript as written.

Results

Clinical characteristics of ARIC and CHS participants at the baseline examinations

Baseline demographic and cardiovascular risk factors are shown by SCD status and by self-reported race in Table 1. As expected, many of the well established cardiovascular risk factors were significantly associated with SCD risk. Among both whites and blacks, individuals who died from SCD were significantly older, had higher systolic blood pressure, and fibrinogen and lower HDL cholesterol. They were also more likely to be male, smokers, less well educated, and have a history of diabetes, hypertension, and myocardial infarction at the baseline examination. Individuals who experienced SCD had significantly longer mean QT intervals at baseline, prior to the event. In whites, the age-, sex-, heart-rate, and study-adjusted mean QT duration was 403 ms for those who did not ultimately have SCD and 411 for those who did (P<0.001); and the corresponding values were 402 and 410 ms, respectively, in blacks (P<0.001).

Table 1
Baseline characteristics of 14,737 white and 4,558 black ARIC and CHS participants by SCD status and race

While the two studies are largely comparable, there were modest differences in the significance and magnitude of associations between cardiovascular risk factors and SCD risk between the two studies. For example, smoking, BMI, and total cholesterol were not significantly associated with SCD in CHS but were in ARIC (Table 2).

Table 2
Baseline characteristics of 14,309 ARIC and 4,986 CHS participants by SCD status

Association between QT interval and SCD

Over a median follow up of 14.1 years in ARIC and 12.2 years in CHS, 334 whites and 164 blacks experienced sudden cardiac death (cumulative incidence rate per 1,000 person-years: 2.0 overall; 1.4 in ARIC and 4.5 in CHS; 1.8 in whites and 2.9 in blacks). Given the older age of the CHS cohort, 222 of the 498 SCD events (45%) occurred in CHS while 276 events occurred in ARIC. Of all CHD mortality adjudicated (N = 985), 40.2% were classified as definite SCD and 10.4% as possible SCD. The majority of cases of SCD occurred out of hospital (90%).

Among whites, longer QT interval was associated with the development of SCD after adjust for age, sex, heart rate, and study. Compared to whites in the first quintile of QT interval, those in the 2nd, 3rd, 4th, and 5th quintiles were 1.43 (95% CI 0.97 to 2.13), 1.77 (95% CI 1.16 to 2.69), 2.33 (95% CI 1.49 to 3.64), and 3.58 (95% CI 2.20 to 5.81) times more likely to have suffered SCD, respectively (P for trend <0.0001). Similar dose response relationship was observed when the analysis was repeated using QT deciles (P for trend <0.0001). Among blacks, longer QT interval was also associated with SCD risk; however, the dose response relationship was less apparent, possibly due to the smaller number of cases. Compared to blacks in the first quintile of QT interval, those in the 2nd, 3rd, 4th, and 5th quintiles were about 1.75 (95% CI 1.07 to 2.88), 2.04 (95% CI 1.19 to 3.49), 1.62 (95% CI 0.86 to 3.05), and 2.54 (95% CI 1.34 to 4.82) times more likely to have suffered SCD (P for trend =0.02).

Association between 19 SNPs in NOS1AP and QT interval

Allele frequencies of the 19 SNPs examined are shown in Table 3 by ethnicity. The LD pattern of these 19 SNPs were quite different in blacks as compared with whites, with considerably less LD among SNPs in blacks (Figure 1). Eleven of 19 SNPs examined were significantly associated with age-, sex-, and heart rate-adjusted QT interval in whites with P≤0.003. However, no SNPs were significantly associated with QT interval in blacks.

Figure 1Figure 1
Plots showing the linkage disequilibrium (LD) pattern and association results for both QT interval and SCD in whites (A) and blacks (B) for 19 SNPs genotyped to tag the NOS1AP locus and surrounding region that exhibited the strongest association with ...
Table 3
Genotypic association between 19 SNPs in NOS1AP and QT interval and SCD by self-reported race

In whites, the most significant SNP in the present study was rs16847548 (P=2.2 × 10−18), which is in LD with rs4657139 but was not typed in previous studies1115, 17. The frequency of the C allele of rs16847548 was 0.22 in whites. After adjusting for age, sex, and heart rate, the mean QT interval of individuals with TT, TC, and CC genotypes at rs16847548 were 399, 401, and 403 ms respectively in ARIC; 411, 414, and 416 ms respectively in CHS; and 402, 404, and 407 ms respectively in the combined dataset. This effect size of approximately 5 ms difference between the two homozygous groups is consistent with our previous observations. In whites, the percent variation (R2) in the QT interval distribution (uncorrected QT interval) that was explained by rs16847548 was 0.2% in both the individual studies and the combined dataset. In comparison, the R2 associated with other variables was: 0.5% for age, 0.1% for sex, 0.2% for diabetes, 0.4% for history of MI at baseline, and 67% for heart rate in white ARIC participants. Among white CHS participants, the R2 associated with age, sex, diabetes, history of MI at baseline, and heart rate were 0.03%, 1.1%, 0.4%, 0.5%, and 63%, respectively.

Associations between NOS1AP genotypes and SCD

Consistent with the observation of longer mean QT interval associated with the C allele of rs16847548, this allele was also significantly associated with increased risk of SCD in whites. Indeed, only 179 of the 8,905 individuals (2%) carrying the TT genotype at rs16847548 suffered from SCD (Table 4), whereas, 25 of the 706 (3.5%) individuals with the CC genotype experienced SCD. In whites, the crude relative hazards that were estimated using a co-dominant model suggested a dose-response relationship between copies of the C allele at rs16847548 and SCD. Using an additive model, the age-, sex-, and study-adjusted RH for each C allele was 1.31, 95% CI 1.10 to 1.56; P=0.002 (Table 4).

Table 4
Unadjusted and adjusted relative hazard (RH) of SCD by rs16847548 and rs12567209 genotypes in whites from ARIC and CHS

In addition, a downstream neighboring SNP, rs12567209, not correlated with rs16847548 (r2 = 0.02), was also associated with SCD in whites (age-, sex-, and study-adjusted RH for each A allele =0.57 assuming an additive model, 95% CI 0.39 to 0.83; P=0.003). Due to the low frequency of the A allele (MAF=0.07), a dominant model was also used for the analysis of rs12567209. The age-, sex-, and study-adjusted relative hazard of SCD comparing those with at least one copy of the A allele to those with the GG genotype was 0.53 (95% CI 0.36 to 0.79; P=0.002), thus both the additive and the dominant models were consistent with the data. The present study is not able to distinguish whether one model was a better fit than the other (additive model shown in Table 4). Surprisingly, rs12567209 was only modestly associated with QT interval (Table 3), suggesting that the effect on risk for SCD was not necessarily conveyed through modulation of QT interval. The mean age-, sex-, heart-rate, and study-adjusted QT interval for GG, AG, and AA genotypes were 403, 403, and 401 ms respectively (P for additive =0.05; P for dominant model =0.08).

To demonstrate the independent effect on risk of SCD for rs16847548 and rs12567209, we included both SNPs in the same model and found that both SNPs remained associated with SCD. Moreover, there was no significant interaction between these two SNPs (P for interaction=0.68 adjusted for age, sex, and study). Assuming an additive model for both SNPs and after adjusting for age, sex, and study, the RH for each copy of the C allele of rs16847548 was 1.27 (95% CI 1.06 to 1.51; P=0.008) and 0.62 (95% CI 0.42 to 0.90; P=0.012) for each copy of the A allele of rs12567209 (Table 4). The associations between both SNPs and SCD risk were relatively consistent across the two studies although the study-specific p values did not reach the Bonferroni-corrected significance level due to their smaller sample sizes. The corresponding RH for rs16847548 was 1.32 (95% CI 1.01 to 1.70; P=0.04) for ARIC and 1.24 (95% CI 0.98 to 1.57; P=0.08) for CHS. The corresponding RH for rs12567209 was 0.84 (95% CI 0.52 to 1.39; P=0.51) for ARIC and 0.44 (95% CI 0.25 to 0.79; P=0.006).

No SNPs were significantly associated with SCD in blacks at the α=0.003 level (Table 3). The RH for each C allele of rs16847548 was 1.07 (95% CI 0.81 to 1.41; P=0.64) after adjusting for age, sex, and study. On the other hand, the RH for each A allele of rs12567209 was 1.40 (opposite direction as the association in whites; 95% CI 0.97 to 2.03; P=0.07) after adjusting for age, sex, and study (Supplementary Table 1).

Multivariate genotype association analyses with SCD

To explore whether the effect of NOS1AP SNPs on the risk of SCD is entirely mediated through modulation of the QT interval, we added both QT interval and heart rate as variables into the Cox proportional hazards model in whites, after adjusting for age, sex, and study. The relative hazard of SCD associated with each additional copy of the C allele at rs16847548 decreased from 1.27 (model 2 in Table 4) to 1.22 (model 3 in Table 4), and the relative hazard for each additional A allele at rs12567209 changed from 0.62 to 0.63. Further adjustment for existing cardiovascular risk factors that were associated with SCD modestly attenuated the significance of the associations for rs16847548 (RH=1.17; 95% CI 0.97 to 1.42).

Additional analyses were performed to examine the impact of either co-morbidities or other SNPs on the robustness of the association between rs16847548 and rs12567209 and SCD risk in whites. First, exclusion of whites with a previous history of MI strengthened both associations (for each C allele of rs16847548 RH=1.39 for model 1, 95% CI 1.13 to 1.69; P=0.002; for each A allele of rs12567209 RH=0.49 for model 1; 95% CI 0.30 to 0.78; P=0.003). Second, exclusion of 623 individuals with electrocardiographic QRS complex >120 ms, which is indicative of a bundle branch block or other conduction defect, also resulted in a stronger association between rs16847548 and SCD risk (RH=1.40 for model 1; 95% CI 1.16 to 1.68; P<0.001). However, when QRS duration was included in the fully adjusted model (model 3), the relative hazard changed minimally from 1.23 to 1.25 (95% CI 1.04 to 1.51). The age-, sex-, and study-adjusted RH for rs12567209 changed from 0.57 to 0.62 (95% CI 0.42 to 0.92; P=0.02) upon exclusion of QRS complex >120 ms. Finally, among whites, none of the tests of interaction between either SNP, SCD, and known cardiovascular risk factors (study, history of MI, sex, diabetes, age at last follow up, diabetes, hypertension, family history of CVD, obesity, dyslipdemia, and smoking) was statistically significant (Supplemental Tables S2 and S3).

Lack of association between rs16847548 and rs12567209 and non-sudden cardiac mortality in whites

It is possible that the association between rs16847548 and/or rs12567209 with SCD is due to an association with overall CHD mortality. Therefore, survival analyses were also conducted for CHD mortality that was not coded as SCD (non-SCD CHD mortality) and all other mortality that were neither SCD nor CHD (non-SCD & non-CHD) mortality. Figure 2 shows the cumulative incidences of SCD, non-SCD CHD, and non-SCD & non-CHD mortality, while accounting for each other as competing case of death, by rs16847548 or rs12567209 in whites. The cumulative incidence of SCD per 1,000 person-years of whites with the TT, TC, and CC genotypes at rs16847548 were 1.5, 2.0, and 2.8, respectively (Figure 2A). On the other hand, rs16847548 was not associated with non- sudden CHD mortality (age-, sex-, and study-adjusted RH=0.98, 95% CI 0.83 to 1.17; P=0.86; Figure 2C) nor with non-SCD & non-CHD mortality (age-, sex-, and study-adjusted RH=1.00, 95% CI 0.94 to 1.07, P=0.94; Figure 2E). For rs12567209, the cumulative incidence of SCD per 1,000 person-years of whites with the GG, AG, and AA genotypes were 1.9, 0.9, and 2.2, respectively (Figure 2B). As for rs16847548, no association was observed with non- sudden CHD mortality (age-, sex-, and study-adjusted RH=0.83, 95% CI 0.62 to 1.11); P=0.21; Figure 2D) nor with non-SCD & non-CHD mortality (age-, sex-, and study-adjusted RH=99 95% CI 0.89 to 1.10, P=0.86; Figure 2F).

Figure 2
Cumulative incidence curves, accounting for competing cause of death, in whites. Kaplan-Meier survival curves in whites by rs16847548 (A,C,E) and by rs1267209 (B,D,F). (A,B) SCD, (C,D) non-SCD CHD mortality, and (E,F) non-SCD & non-CHD mortality. ...

Discussion

In this study, common sequence variations in NOS1AP are associated with both inter-individual variation in the QT interval and risk of SCD in whites from two large cohorts of adults in the U.S. The at-risk allele (C at rs16847548) is common, with 39% of the general white population carrying one copy of the C allele and 5% carrying two copies. More specifically, in the combined population, white U.S. adults carrying the CC genotype at rs16847548 of NOS1AP were about 72% more likely to die of SCD and had a mean QT interval that was approximately 5 ms longer than their counterparts with the TT genotype, even after accounting for age, sex, and heart rate. On the other hand, the less common A allele of rs12567209 (~13% of the general white population are carriers) was independently associated with a decreased risk of SCD (RH=0.57) in whites. Notably, both of the genetic effects were specific for SCD rather than other forms of death from CHD. Finally, in spite of demographic differences between ARIC and CHS whites, the genetic effect estimates were comparable in the two populations separately and there was no evidence of significant heterogeneity (Tables S2 and S3).

The present study identifies a novel gene, NOS1AP, along with a new set of cellular interactions, which can potentially affect SCD risk in the general population. The multivariate analyses shows that even after adjusting for QT interval and heart rate a significant association still remains between both rs16847548 and rs12567209 and SCD risk. This result was somewhat unanticipated given that our initial hypothesis evolved from a model in which SNPs influence QT interval, and that increasing QT interval would increase risk for SCD. However, given that adjusting for QT interval largely does not attenuate the risk for SCD associated with these SNPs, this suggests an alternate model, in which these SNPs modulate an unmeasured, or hidden, factor, which itself modulates both QT interval and risk for SCD, and need not do so equivalently. However, the possibility that the remaining significant association between these SNPs and SCD risk is due to the QT interval as assessed by ECG representing an imperfect measure of cardiac repolarization, the potential misclassification of the actual QT interval measurements and SCD definition, or the presence of additional genetic variation (i.e. not having identified the causal SNPs) cannot be excluded. Although a recent study reported no association between NOS1AP and SCD in the Rotterdam Study, it is important to note that rs16847548 and rs12567209 were not directly studied nor efficiently tagged and that the number of SCD events was small12. In addition, the positive association between rs16847548 of NOS1AP and SCD risk supports the approach of using either precursors or intermediate phenotypes in genetic studies of complex diseases7, 32, 33 as NOS1AP was first identified to be a candidate gene for SCD through a previous GWAS of QT interval11.

The results of the present study, together with previous reports of associations between NOS1AP and the QT interval in multiple populations of European descent11, 1317, suggest novel and potentially causal mechanisms linking NOS1AP and SCD risk. As an adapter protein, the gene product of NOS1AP (CAPON) serves to physically bridge neuronal nitric oxide synthase (nNOS) and its targets and modulator proteins. In guinea pig ventricular myocytes, CAPON is localized near ryanodine receptors, and the over expression of CAPON results in shortening of the cardiac action potential, a decrease in L-type Ca current and a smaller increase in the delayed rectifier potassium current, IKr, resulting in prolongation of the QT interval34.

Several limitations are warranted in the interpretation of these findings. First, although we have identified the association of sequence variation at the NOS1AP locus with QT interval and SCD risk, it is not known whether we have identified the functional variants. For example, it is likely that rs12567209 is only in linkage disequilibrium with the causal SNP since its association with SCD in whites was in the opposite direction as its association in blacks. Even though rs16847548 had the strongest association (judging by p-values of all 19 SNPs), with both QT interval and SCD in both the ARIC and CHS cohorts, simulation studies have shown that the causative SNP may not necessarily have the smallest P-value since P-values fluctuate by chance due to the nature of random sampling, dependent on the sample size and allele frequency35. Thus, it is possible that rs16847548 is also in linkage disequilibrium with another causal SNP. Second, no significant association between the 19 NOS1AP SNPs and QT interval or SCD was observed in the black participants at a conservative α=0.003. The discordance in the associations between blacks and whites may represent the result of lower statistical power (due to inappropriate tagging SNPs and smaller numbers of events) in the blacks or it may represent a real genetic difference. Given the observed carrier frequency for rs16847548 of 34% (based on allele frequency of 0.19 in blacks) and assuming an overall genotypic relative hazard of 1.37, as observed in whites, at least 649 SCD cases in the blacks would be necessary for the study to have 80% power with an alpha of 0.00331. If the correlation between our genotyped SNPs and the putative ungenotyped functional variant is lower in blacks, then we would have less power to detect an effect in blacks. On the other hand, it is possible that the causal allele in blacks is not rs16847548 or that there exists only one causal allele in NOS1AP but the pattern of linkage disequilibrium between the causal variant and rs16847548 differs between blacks and whites. Third, despite corroborating functional data from guinea pig cardiomyocytes34, it is still possible that the associated variants in the NOS1AP locus actually influence (or are in LD with) a distant gene rather than NOS1AP, as this “action at a distance” has on rare occasions been observed at other rare diseases36. Fourth, s with all genetic association studies of complex traits, an independent replication study of comparable size and phenotype is not only the best defense against possible false positive reporting in our study but is necessary before certainty of the observed associations can be established. Lastly, the present study was not ideal for assessing the utility of genetic risk prediction among those already at high risk as the number of high-risk individuals is relatively modest in these studies.

In summary, in this study, we report that sequence variations in NOS1AP, a novel candidate gene that was previously identified through GWAS of the QT interval11, are associated with both QT interval and the subsequent risk of SCD in a large cohort of 14,737 white U.S. adults. As expected, individuals carrying the at-risk allele at rs16847548 had a modest increased risk of SCD, with each allele increasing SCD risk by about 30% compared to those who did not carry the risk allele. On the other hand, the A minor allele of rs12567209 was associated with a reduced risk of SCD, and the associations of these two SNPs were independent of each of other. Although the genetic effects described here are modest, if replicated in other populations, this effort may be an important step towards the identification of a panel of susceptibility alleles that may potentially used for risk assessment in the general population. Future studies that explore the pathways mediating the association between variations of NOS1AP and SCD risk will also be crucial and may shed light on both targeted prevention strategies and novel therapeutic targets for abnormal cardiac repolarization and SCD risk.

Supplementary Material

Acknowledgments

The authors thank the staff and participants of the ARIC study and CHS for their important contributions. The authors would like to thank Ms. Ashley O’Connor for her effort in genotyping theses two large cohorts.

Funding Sources

This work was supported by the Donald W. Reynolds Foundation (WHLK, DEA, WP, BB, BD, GT, EM, PMS, AC), NHGRI HG02757 (AC), the Leducq Foundation (EM, PMS, AC), and K01DK067207 (WHLK).

The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, and N01-HC-55022.

The Cardiovascular Heart Study is supported by contracts N01-HC-35129, N01-HC-45133, N01-HC-75150, N01-HC-85079 through N01-HC-85086, N01 HC-15103, N01 HC-55222, and U01 HL080295 from the National Heart, Lung, and Blood Institute, with additional contribution from the National Institute of Neurological Disorders and Stroke.

Footnotes

Short summary of potential clinical impact

Nearly half of all coronary heart disease (CHD) deaths are sudden. Family history of sudden cardiac death (SCD) is a powerful risk factor for SCD; however, the genetic factors underlying SCD in the general population are largely unknown. The electrocardiographic QT interval is associated with risk of SCD. A previous genome-wide association study reported that allelic variants in NOS1AP, which encodes a ligand of neuronal nitric oxide synthase, are associated with the QT interval in white adults. The present analysis was conducted to validate the association between NOS1AP variants and the QT interval and to further examine the association with SCD in a combined population of 19,295 black and white adults from two population-based cohort studies. Among whites, we found that multiple SNPs in NOS1AP were associated with adjusted QT interval in whites (P<0.0001), and two SNPs were independently associated with SCD. One SNP, with a minor allele frequency (MAF) of 22%, was associated with a 31% greater risk of SCD for each copy of the variant allele, while a neighboring SNP (MAF 7%) was associated with a 43% lower risk for SCD. No associations between SNPs in NOS1AP and either QT interval or SCD were observed in blacks. Although the genetic effects described here are modest, if replicated in other populations, this effort may represent one step towards using genetic risk markers, along with other risk factors, to help identify patients who warrant our most aggressive SCD preventive strategies.

Disclosure Statement

None

References

1. Chugh SS, Jui J, Gunson K, Stecker EC, John BT, Thompson B, Ilias N, Vickers C, Dogra V, Daya M, Kron J, Zheng ZJ, Mensah G, McAnulty J. Current burden of sudden cardiac death: multiple source surveillance versus retrospective death certificate-based review in a large U.S. community. J Am Coll Cardiol. 2004;44:1268–1275. [PubMed]
2. Zheng ZJ, Croft JB, Giles WH, Mensah GA. Sudden cardiac death in the United States, 1989 to 1998. Circulation. 2001;104:2158–2163. [PubMed]
3. Myerburg RJ, Castellanos A. Emerging paradigms of the epidemiology and demographics of sudden cardiac arrest. Heart Rhythm. 2006;3:235–239. [PubMed]
4. Dekker LR, Bezzina CR, Henriques JP, Tanck MW, Koch KT, Alings MW, Arnold AE, de Boer MJ, Gorgels AP, Michels HR, Verkerk A, Verheugt FW, Zijlstra F, Wilde AA. Familial sudden death is an important risk factor for primary ventricular fibrillation: a case-control study in acute myocardial infarction patients. Circulation. 2006;114:1140–1145. [PubMed]
5. Friedlander Y, Siscovick DS, Weinmann S, Austin MA, Psaty BM, Lemaitre RN, Arbogast P, Raghunathan TE, Cobb LA. Family history as a risk factor for primary cardiac arrest. Circulation. 1998;97:155–160. [PubMed]
6. Jouven X, Desnos M, Guerot C, Ducimetiere P. Predicting sudden death in the population: the Paris Prospective Study I. Circulation. 1999;99:1978–1983. [PubMed]
7. Arking DE, Chugh SS, Chakravarti A, Spooner PM. Genomics in sudden cardiac death. Circ Res. 2004;94:712–723. [PubMed]
8. Vrtovec B, Delgado R, Zewail A, Thomas CD, Richartz BM, Radovancevic B. Prolonged QTc interval and high B-type natriuretic peptide levels together predict mortality in patients with advanced heart failure. Circulation. 2003;107:1764–1769. [PubMed]
9. Newton-Cheh C, Larson MG, Corey DC, Benjamin EJ, Herbert AG, Levy D, D’Agostino RB, O’Donnell CJ. QT interval is a heritable quantitative trait with evidence of linkage to chromosome 3 in a genome-wide linkage analysis: The Framingham Heart Study. Heart Rhythm. 2005;2:277–284. [PubMed]
10. Straus SM, Kors JA, De Bruin ML, van der Hooft CS, Hofman A, Heeringa J, Deckers JW, Kingma JH, Sturkenboom MC, Stricker BH, Witteman JC. Prolonged QTc interval and risk of sudden cardiac death in a population of older adults. J Am Coll Cardiol. 2006;47:362–367. [PubMed]
11. Arking DE, Pfeufer A, Post W, Kao WH, Newton-Cheh C, Ikeda M, West K, Kashuk C, Akyol M, Perz S, Jalilzadeh S, Illig T, Gieger C, Guo CY, Larson MG, Wichmann HE, Marban E, O’Donnell CJ, Hirschhorn JN, Kaab S, Spooner PM, Meitinger T, Chakravarti A. A common genetic variant in the NOS1 regulator NOS1AP modulates cardiac repolarization. Nat Genet. 2006;38:644–651. [PubMed]
12. Aarnoudse AJ, Newton-Cheh C, de Bakker PI, Straus SM, Kors JA, Hofman A, Uitterlinden AG, Witteman JC, Stricker BH. Common NOS1AP variants are associated with a prolonged QTc interval in the Rotterdam Study. Circulation. 2007;116:10–16. [PubMed]
13. Eijgelsheim M, Aarnoudse AL, Rivadeneira F, Kors JA, Witteman JC, Hofman A, van Duijn CM, Uitterlinden AG, Stricker BH. Identification of a Common Variant at the NOS1AP Locus Strongly Associated to QT - Interval Duration. Hum Mol Genet. 2008 [PubMed]
14. Lehtinen AB, Newton-Cheh C, Ziegler JT, Langefeld CD, Freedman BI, Daniel KR, Herrington DM, Bowden DW. Association of NOS1AP genetic variants with QT interval duration in families from the Diabetes Heart Study. Diabetes. 2008;57:1108–1114. [PubMed]
15. Post W, Shen H, Damcott C, Arking DE, Kao WH, Sack PA, Ryan KA, Chakravarti A, Mitchell BD, Shuldiner AR. Associations between genetic variants in the NOS1AP (CAPON) gene and cardiac repolarization in the old order Amish. Hum Hered. 2007;64:214–219. [PMC free article] [PubMed]
16. Raitakari OT, Blom-Nyholm J, Koskinen TA, Kahonen M, Viikari JS, Lehtimaki T. Common variation in NOS1AP and KCNH2 genes and QT interval duration in young adults. The Cardiovascular Risk in Young Finns Study. Ann Med. 2008:1–8. [PubMed]
17. Tobin MD, Kahonen M, Braund P, Nieminen T, Hajat C, Tomaszewski M, Viik J, Lehtinen R, Ng GA, Macfarlane PW, Burton PR, Lehtimaki T, Samani NJ. Gender and effects of a common genetic variant in the NOS1 regulator NOS1AP on cardiac repolarization in 3761 individuals from two independent populations. Int J Epidemiol. 2008;37:1132–1141. [PubMed]
18. The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. The ARIC investigators. Am J Epidemiol. 1989;129:687–702. [PubMed]
19. Fried LP, Borhani NO, Enright P, Furberg CD, Gardin JM, Kronmal RA, Kuller LH, Manolio TA, Mittelmark MB, Newman A. The Cardiovascular Health Study: design and rationale. Ann Epidemiol. 1991;1:263–276. [PubMed]
20. Ives DG, Fitzpatrick AL, Bild DE, Psaty BM, Kuller LH, Crowley PM, Cruise RG, Theroux S. Surveillance and ascertainment of cardiovascular events. The Cardiovascular Health Study. Ann Epidemiol. 1995;5:278–285. [PubMed]
21. Price TR, Psaty B, O’Leary D, Burke G, Gardin J. Assessment of cerebrovascular disease in the Cardiovascular Health Study. Ann Epidemiol. 1993;3:504–507. [PubMed]
22. Sotoodehnia N, Siscovick DS, Vatta M, Psaty BM, Tracy RP, Towbin JA, Lemaitre RN, Rea TD, Durda JP, Chang JM, Lumley TS, Kuller LH, Burke GL, Heckbert SR. Beta2-adrenergic receptor genetic variants and risk of sudden cardiac death. Circulation. 2006;113:1842–1848. [PubMed]
23. Rautaharju PM, Park LP, Chaitman BR, Rautaharju F, Zhang ZM. The Novacode criteria for classification of ECG abnormalities and their clinically significant progression and regression. J Electrocardiol. 1998;31:157–187. [PubMed]
24. General Description and Study Management. ARIC Protocol. 1987
25. Psaty BM, Kuller LH, Bild D, Burke GL, Kittner SJ, Mittelmark M, Price TR, Rautaharju PM, Robbins J. Methods of assessing prevalent cardiovascular disease in the Cardiovascular Health Study. Ann Epidemiol. 1995;5:270–277. [PubMed]
26. A haplotype map of the human genome. Nature. 2005;437:1299–1320. [PMC free article] [PubMed]
27. de Bakker PI, Yelensky R, Pe’er I, Gabriel SB, Daly MJ, Altshuler D. Efficiency and power in genetic association studies. Nat Genet. 2005;37:1217–1223. [PubMed]
28. www.biotrove.com. 2007.
29. Bazett H. An analysis of the time-relations of electrocardiograms. Heart-A Journal for the Study of the Circulation. 1920;7:353–370.
30. Gooley TA, Leisenring W, Crowley J, Storer BE. Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. Stat Med. 1999;18:695–706. [PubMed]
31. Schoenfeld DA. Sample-size formula for the proportional-hazards regression model. Biometrics. 1983;39:499–503. [PubMed]
32. Cohen JC, Boerwinkle E, Mosley TH, Jr, Hobbs HH. Sequence variations in PCSK9, low LDL, and protection against coronary heart disease. N Engl J Med. 2006;354:1264–1272. [PubMed]
33. Lange LA, Carlson CS, Hindorff LA, Lange EM, Walston J, Durda JP, Cushman M, Bis JC, Zeng D, Lin D, Kuller LH, Nickerson DA, Psaty BM, Tracy RP, Reiner AP. Association of polymorphisms in the CRP gene with circulating C-reactive protein levels and cardiovascular events. JAMA. 2006;296:2703–2711. [PubMed]
34. Chang KC, Barth AS, Sasano T, Kizana E, Kashiwakura Y, Zhang Y, Foster DB, Marban E. CAPON modulates cardiac repolarization via neuronal nitric oxide synthase signaling in the heart. Proc Natl Acad Sci U S A. 2008;105:4477–4482. [PubMed]
35. Huang Q, Morrison AC, Boerwinkle E. Linkage disequilibrium structure and its impact on the localization of a candidate functional mutation. Genet Epidemiol. 2001;21 (Suppl 1):S620–625. [PubMed]
36. Lettice LA, Horikoshi T, Heaney SJ, van Baren MJ, van der Linde HC, Breedveld GJ, Joosse M, Akarsu N, Oostra BA, Endo N, Shibata M, Suzuki M, Takahashi E, Shinka T, Nakahori Y, Ayusawa D, Nakabayashi K, Scherer SW, Heutink P, Hill RE, Noji S. Disruption of a long-range cis-acting regulator for Shh causes preaxial polydactyly. Proc Natl Acad Sci U S A. 2002;99:7548–7553. [PubMed]