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(22):2703-2711. doi:10.1001/jama.296.22.2703
Author Affiliations: Departments of Genetics (Drs L. Lange and E. Lange) and Biostatistics (Drs E. Lange, Zeng, and Lin), University of North Carolina, Chapel Hill; Fred Hutchinson Cancer Research Center, Seattle, Wash (Dr Carlson); Departments of Epidemiology (Drs Hindorff, Psaty, and Reiner, and Mr Bis), Medicine (Dr Psaty), Genome Sciences (Dr Nickerson), and Laboratory Medicine (Dr Reiner), University of Washington, Seattle; Center on Aging and Health, Johns Hopkins University School of Medicine, Baltimore, Md (Dr Walston); Departments of Pathology (Drs Cushman and Tracy, and Mr Durda), Medicine (Drs Cushman and Tracy), and Biochemistry (Dr Tracy), University of Vermont College of Medicine, Burlington; and Department of Epidemiology, University of Pittsburgh, Pittsburgh, Pa (Dr Kuller).
Context C-reactive protein (CRP) is an inflammation protein that may play a role in the pathogenesis of cardiovascular disease (CVD).
Objective To assess whether polymorphisms in the CRP gene are associated with plasma CRP, carotid intima-media thickness (CIMT), and CVD events.
Design, Setting, and Participants In the prospective, population-based Cardiovascular Health Study, 4 tag single-nucleotide polymorphisms (SNPs) (1919A/T, 2667G/C, 3872G/A, 5237A/G) were genotyped in 3941 white (European American) participants and 5 tag SNPs (addition of 790A/T) were genotyped in 700 black (African American) participants, aged 65 years or older, all of whom were without myocardial infarction (MI) or stroke before study entry. Median follow-up was 13 years (1989-2003).
Main Outcome Measures Baseline CIMT; occurrence of MI, stroke, and CVD mortality during follow-up.
Results In white participants, 461 incident MIs, 491 incident strokes, and 490 CVD-related deaths occurred; in black participants, 67 incident MIs, 78 incident strokes, and 75 CVD-related deaths occurred. The 1919T and 790T alleles were associated with higher CRP levels in white and black participants, respectively. The 3872A allele was associated with lower CRP levels in both populations, and the 2667C allele was associated with lower CRP levels in white participants only. There was no association between CIMT and any CRP gene polymorphism in either population. In white participants, the 1919T allele was associated with increased risk of stroke for TT vs AA (hazard ratio [HR], 1.40; 95% confidence interval [CI], 1.06-1.87) and for CVD mortality (HR, 1.40; 95% CI, 1.10-1.90). In black participants, homozygosity for the 790T allele was associated with a 4-fold increased risk of MI compared with homozygosity for the 790A allele (95% CI, 1.58-10.53). The minor alleles of the 2 SNPs associated with lower plasma CRP concentration in white participants (2667C and 3872A) were associated with decreased risk of CVD mortality.
Conclusions Genetic variation in the CRP gene is associated with plasma CRP levels and CVD risk in older adults.
C-reactive protein (CRP), an acute phase reactant, plays an important role in acute and chronic inflammation. Inflammation contributes to all phases of atherosclerosis, from fatty streak initiation to cardiovascular disease (CVD) events.1 Plasma CRP levels measured at a single point predict future incident CVD events, such as myocardial infarction (MI) and stroke,2- 7 leading to interest in the use of CRP as a biomarker of CVD risk.8 Although some data from in vitro observations and experimental models support a role of CRP in atherogenesis,9 the direct involvement of CRP remains in question.
Plasma CRP levels are under genetic influence.10,11 In healthy, young white (European American) and black (African American) adults, CRP gene haplotypes have been shown to be strongly associated with plasma concentration.12- 14 This led us to hypothesize that CRP is causally involved in atherosclerosis, and that alleles associated with higher CRP would confer an increased risk of CVD. Like the variant in PCSK9 (proprotein convertase subtilisin/kexin type 9), which was associated not only with lifelong lowering of low-density lipoprotein cholesterol but also with a reduced risk of CVD, CRP gene variants that affect CRP concentration levels may reflect lifetime exposure to CRP more accurately than CRP serum concentrations measured at a single point in time.15
The Cardiovascular Health Study (CHS) design, recruitment, and data collection have been previously described.16 There were 5888 men and women (≥65 years) recruited from 4 US field centers: Forsyth County, North Carolina; Sacramento County, California; Washington County, Maryland; and Pittsburgh, Pennsylvania. Race/ethnicity was defined by self-report. The original cohort (n = 5201) was recruited from 1989 to 1990. A second cohort (n = 687) of black participants was recruited from 1992 to 1993. All participants were followed up through 2003, unless lost to follow-up, dropped out of study, or had an event before that time. Only individuals who reported race as white or black were included in our study; we excluded participants with prior MI or stroke (n = 765), or for whom DNA was not collected or who did not consent to the use of their DNA (n = 482). All procedures were conducted under institutionally approved protocols for use of human subjects, and all participants provided written informed consent.
The CHS baseline evaluation included demographic, lifestyle, and medical histories, physical examination, and fasting blood collection.16,17 Common carotid artery wall intima-media thickness (CIMT) was determined at the baseline examination by high-resolution B-mode ultrasonography.18 The mean maximal intima-media thickness of the common or internal carotid wall was calculated from 4 scans on each carotid artery (right and left side) for each participant. C-reactive protein was measured on stored EDTA plasma from the baseline examination by using a high-sensitivity enzyme-linked immunosorbent assay (coefficient of variation, 6.2%).19 Plasma CRP was also measured 3 years after the baseline examination on a subset of 3256 white participants by immunonephelometry (coefficient of variation, 3.0%; BN-II instrument, Dade-Behring, Deerfield, Ill). Both the enzyme-linked immunosorbent assay and nephelometric assays are standardized according to the World Health Organization reference plasma and yield highly correlated values.
Details of event ascertainment in CHS have been previously published.20 Participants, family members, or other previously identified informants reported new cardiovascular events during semiannual contacts. Medical records were reviewed, and events were adjudicated by a physician review panel. Criteria for MI included history of chest pain, cardiac enzyme levels, and characteristic changes on serial electrocardiograms. Criteria for stroke included onset of symptoms, duration of deficits, and findings on computed tomography or magnetic resonance imaging. Adjudicated events occurring through June 30, 2003, were available. Participants were followed up for up to 14 years (median: 13 years for white cohort and 10 years for black cohort).
Primary clinical end points for this analysis were incident non-procedure–related fatal or nonfatal MI, stroke, and CVD mortality, defined as fatal event in which death was adjudicated due to atherosclerotic coronary heart disease or cerebrovascular disease, including definite fatal MI, definite fatal stroke, and definite or probable fatal coronary heart disease.20
Based on linkage disequilibrium and haplotype patterns across the CRP gene locus,12 we selected tag single-nucleotide polymorphisms (SNPs) for genotyping. The SNP numbering is relative to GenBank accession AF449713. Four tag SNPs (1919 [rs1417938], located in intron 1; 2667 [rs1800947], a synonymous SNP located in exon 2; and 3872 [rs1205] and 5237 [rs2808630], both located in the 3′ flanking region) were typed in both the white and black CHS samples. A fifth tag SNP, 790 (rs3093058), located in the 5′ promoter region is prevalent only among black participants, and therefore was genotyped only in the black CHS samples. Tag SNPs were genotyped using TaqMan Assays By Design (Applied Biosystems, Foster City, Calif) under standard conditions. Probe and primer sequences are available from the authors by request.
Each of the CRP polymorphisms was assessed to determine if the observed genotype frequencies were consistent with expected Hardy-Weinberg proportions by using Pearson χ2 tests. Marker-marker linkage disequilibrium was assessed using Lewontin’s D′ statistic21 and the squared correlation statistic Δ2.22 We considered statistical significance at the α=.05 level, before and after adjustment for multiple tests.
Genotype Association. All analyses were stratified by ethnicity to minimize potential confounding due to population stratification. Non-normal distributions were log-transformed as necessary. To assess SNP associations with plasma CRP levels and CIMT, we used analysis of covariance models (SAS version 8.0; SAS Institute, Cary, NC). Covariate adjustment was made for age, sex, clinic site, body mass index (calculated as weight in kilograms divided by height in meters squared), smoking status, triglycerides, and clinical or subclinical (for CRP only) disease at baseline.23 Genotypes were tested for general association (no mode-of-inheritance assumption) using a 2-df F test. Covariate-adjusted mean values (least square means) were calculated from the regression coefficients.
We performed Cox proportional hazards regression models to test for association between SNPs and risk of incident MI, stroke, or CVD mortality. Two df likelihood ratio tests were conducted to test for a general association between genotype and time-to-CVD event. Cox proportional hazards regression models were performed with 3 levels of covariate adjustment: (1) adjustment for baseline age, sex, clinic site, body mass index, systolic blood pressure, diabetes mellitus, hypertension, and smoking status; (2) level 1 plus additional adjustment for baseline CRP concentration; and (3) level 2 plus additional adjustment for 3-year CRP concentration. To control for population stratification in black participants, we also adjusted black regression models for a variable that reflects the estimated proportion of African ancestry for each individual.24
Haplotype Association. Associations between haplotypes and plasma CRP and CIMT were assessed using an expectation-maximization–derived score test,25 implemented in HAPLO.STAT (http://mayoresearch.mayo.edu/mayo/research/schaid_lab/software.cfm). For time-to-event data, we applied a method that formulates the effects of haplotypes and environmental variables on the time-to-disease occurrence in a Cox proportional hazards regression model by using an expectation-maximization–based algorithm to maximize the likelihood for the relative risks and other parameters.26 Both methods provide a global test of haplotype association and tests of association for individual haplotypes. An additive genetic mechanism was assumed. Only haplotypes with estimated frequencies of more than 5% were included. HAPLO.STAT performs individual haplotype tests that compare each haplotype to all other haplotypes pooled together, whereas the Cox proportional hazards regression model approach compares each haplotype to the most common haplotype.
Sample descriptions, including clinical end points, are shown in Table 1. Black participants had a higher prevalence of diabetes mellitus, current smoking, and higher median baseline plasma CRP compared with white participants.
The observed genotype distributions for the 5 CRP tag SNPs (Table 2) were consistent with Hardy-Weinberg equilibrium. There was little evidence of historical recombination between SNPs (estimated D′ range, 0.98-1.00 in white participants and 0.85-1.00 in black participants). Because the SNPs were selected to maximize linkage disequilibrium coverage across the CRP gene, there was only a modest degree of correlation in genotype data between SNPs (estimated Δ2 range, 0.03-0.22 in white participants and 0.00-0.06 in black participants).
Three of 4 SNPs in the white sample and 2 of 5 SNPs in the black sample were strongly associated with baseline CRP concentration (Table 2). The estimated proportion of variation in plasma CRP level explained by significant single SNP genotypes ranged from 0.7% to 1.8% in white participants and 1.2% to 4.4% in black participants. The estimated proportion of variation explained when including all CRP SNPs was 2.6% in white participants and 6.4% in black participants. Similar results were obtained using plasma CRP measurements obtained after 3 years of follow-up (data available from authors). Haplotype-based analyses of plasma CRP levels yielded results that were consistent with single SNP results (Table 3). There was no association between CIMT and any CRP SNP genotype or haplotype in either the white or black samples.
Among white participants, SNPs 1919, 2667, and 3872 were associated with CVD mortality, and SNP 1919 was associated with increased risk for stroke (Table 4). In black participants, SNPs 790 and 1919 were associated with incident MI and SNP 790 was associated with incident stroke (Table 5). To assess the overall significance of the association between CRP SNPs and the risk of clinical CVD events, we performed joint association analyses that included all SNPs simultaneously. Specifically, we conducted likelihood-ratio tests for Cox proportional hazards regression models, in each case comparing the model with all genotyped SNPs to the null hypothesis model that contained no SNPs. The association between all SNPs and CVD mortality in white participants (P = .004) and MI in black participants (P = .005) were statistically significant at an overall α = .05 level, even after accounting for the 3 events (Bonferroni significance threshold P = .05 divided by 3: P = .02). The association between all SNPs and stroke in black participants (P = .04) was statistically significant before, but not after, accounting for 3 events.
In this study, CRP level predicted CVD risk.7 Because CRP genotype likely influences risk through CRP synthesis, we additionally adjusted for plasma CRP level, which resulted in only a modest attenuation of hazard ratios (Table 4 and Table 5). The CVD risk estimates remained essentially unchanged when we limited analyses to the 3256 white participants with 2 plasma CRP measurements, and incorporated both measurements into the regression models (data available from authors).
Haplotype-based analyses gave results similar to the single SNP analyses in both white and black participants (Table 6). Statistically significant (event-specific) haplotype-based associations were observed for 3 (MI risk in black participants; stroke risk and CVD mortality in white participants) of the 6 ethnicity-event combinations and a trend was observed for a fourth event (stroke risk in black participants). The haplotype-based association with CVD mortality in white participants remained statistically significant (P<.05) after applying a Bonferroni correction to account for multiple testing with respect to the 3 events (white participants).
Our results suggest that a genetic basis may underlie, in part, the relationship between CRP concentration and CVD risk in older adults. The association we observed between CRP tag SNP genotype and plasma CRP concentration confirms previous findings in younger adults.12- 14,27 Together these results provide strong evidence that genetic variation within the CRP gene influences plasma CRP levels over the adult lifespan.
We additionally observed associations between CRP SNP genotypes and risk of future CVD events. Specifically, the 790 T allele was associated with increased risk of MI in black participants, the 1919 T allele was associated with increased risk of stroke and CVD mortality in white participants, and the minor alleles of 2667 and 3872 SNPs were associated with decreased risk of CVD mortality in white participants. The direction of the CVD risk estimates tended to be consistent with associations with plasma CRP concentration measured late in life. The specific CVD outcomes associated with CRP genotype differed by ethnicity. Possible explanations include differences in allele frequencies between white and black participants, reduced statistical power in the smaller black cohort, or chance association. In contrast with the observed association with clinical CVD, there was no evidence of association between the same CRP gene variants and CIMT, a measure of extent of subclinical atherosclerotic disease, consistent with lack of association of plasma CRP levels with this same measurement.28,29
The direction of the observed CRP genotype/haplotype–phenotype associations are also supported by previous functional data on polymorphisms within the CRP promoter region.12,13 The common CRP promoter haplotype in black participants tagged by the 790 T allele was associated with increased basal CRP transcriptional activity in human hepatocytes in vitro.12 Moreover, the 790 T allele and the 1919 T allele (located within the first intron) are both closely linked with high-activity CRP haplotype groups that are defined by another common promoter polymorphism at site 1440.12 In contrast, the minor alleles of the 2667 and 3872 polymorphisms are linked with CRP haplotypes associated with decreased promoter activity.
C-reactive protein may be both a marker of CVD and an active participant in the disease process.30 It is plausible that CRP genotype influences CRP synthesis, which in turn could mediate the onset of clinical CVD events. The association between CRP genotype and clinical CVD, together with the absence of association between CRP genotype and CIMT, suggests a greater involvement of CRP in the transition from subclinical to clinical disease than in atherosclerosis progression. This hypothesis is also consistent with the observation that CRP levels predict risk of stroke in CHS independently of CIMT.5 C-reactive protein is present in atherosclerotic plaques31 and has been associated both experimentally and clinically with plaque instability and thrombosis,32- 34 which might increase the propensity toward acute ischemic events, such as MI or stroke. C-reactive protein may be capable of maintaining a procoagulant phenotype through induction of tissue factor and plasminogen activator inhibitor 1 expression in vascular cells or blood monocytes.35 The complement-dependent ability of CRP to enhance ischemia-induced myocardial tissue damage36 and cerebral infarct size37 is another potential mechanism linking CRP and CVD events.
We found that the hazard ratios of CVD for CRP genotypes were only slightly attenuated on adjustment for plasma CRP concentration. If the mechanism of the association between CVD events and CRP genotype was solely through the gene's influence on CRP concentration, we might expect the attenuation to be more pronounced. However, within our study, individual CRP SNP genotypes account for less than 2% of the interindividual variation in plasma CRP concentration. It may not be that surprising that both CRP concentration and CRP genotype are independently associated with CVD events. Moreover, if chronic inflammation plays a causal role in the development of CVD, plasma CRP measured at a single point in time may not adequately reflect an individual's cumulative inflammatory burden. It is known that there is significant intraindividual and analytical variability in plasma CRP levels.19,38
It also remains remotely possible that the genetic variant responsible for the CRP genotype–CVD event association lies outside of the CRP gene locus and that the observed association is not mediated by plasma CRP levels. The low level of linkage disequilibrium between the structural CRP gene transcript and either flanking gene (serum amyloid P gene or the dual-specificity kinase 23 locus) observed in the HapMap data makes it unlikely that the genetic variant responsible for the CRP genotype–CVD event association lies outside the CRP gene locus. Thus, we consider it substantially more likely that the observed association is mediated directly through the CRP gene than through a linked secondary locus. Finally, the current findings imply that the interpretation of data from CRP genotype–CVD association studies with respect to mendelian randomization and causality of CRP levels may prove to be more complex than previously appreciated.39- 41
Although CRP genotype was strongly associated with plasma CRP concentration, there was little association between CRP genotype and risk of MI or stroke in the Physician's Heath Study, the Framingham Heart Study, or the Rotterdam Study.14,27,42,43 There are a number of possible explanations for the differences in CRP genotype–clinical CVD association between CHS and these other studies. Our analysis included a larger sample size and greater numbers of clinical events; therefore, we had greater power to detect association. In CHS, the association between CRP genotype and CVD risk appeared to be strongest for fatal events. In the Physician's Health Study and the Rotterdam Study, the results were not examined according to fatal vs nonfatal outcomes. It is possible that through an effect on acute phase response CRP genotype is more strongly associated with more severe events. This concept is supported by a recent study performed on 1827 white patients with acute coronary syndrome, in which CRP genotype was observed to be strongly associated with CRP levels measured within 10 days of an event.44 In addition, the participants in the Physician's Health Study were younger and 100% were men. Although it is generally expected that the influence of genetic factors is more important in younger individuals, it may be that older persons have had greater exposure to environmental triggers, thus allowing genetic susceptibilities to become more evident. Indeed, plasma markers of inflammation are under substantial genetic influence, even in older adults.45 Therefore, confirmation of our findings in other large population studies of older adults with large numbers of fatal events will be important for clarifying the role of the CRP gene and risk of CVD.
Corresponding Author: Leslie A. Lange, PhD, Department of Genetics, University of North Carolina, Campus Box 7264, 103 Mason Farm Rd, Chapel Hill, NC 27599-7264 (email@example.com).
Author Contributions: Dr L. Lange had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: L. Lange, Carlson, Durda, Nickerson, Psaty, Tracy, Reiner.
Acquisition of data: Durda, Cushman, Kuller, Psaty, Tracy, Reiner.
Analysis and interpretation of data: L. Lange, Carlson, Hindorff, E. Lange, Walston, Durda, Bis, Zeng, Lin, Kuller, Tracy, Reiner.
Drafting of the manuscript: L. Lange, Carlson, Lin, Nickerson, Reiner.
Critical revision of the manuscript for important intellectual content: L. Lange, Carlson, Hindorff, E. Lange, Walston, Durda, Cushman, Bis, Zeng, Lin, Kuller, Psaty, Tracy, Reiner.
Statistical analysis: L. Lange, Carlson, E. Lange, Zeng, Lin, Psaty, Reiner.
Obtained funding: Kuller, Nickerson, Psaty, Tracy, Reiner.
Administrative, technical, or material support: Hindorff, Durda, Kuller, Nickerson, Psaty, Tracy, Reiner.
Study supervision: Durda, Cushman, Reiner.
Financial Disclosures: None reported.
Funding/Support: This study was supported by contracts N01-HC-85079 through N01-HC-85086, N01-HC-35129, and N01 HC-15103, and grant R01 HL-071862-03 (Dr Reiner) from the National Heart, Lung, and Blood Institute.
Role of the Sponsor: The funding organization participated in the design and conduct of CHS and monitored the study's data collection efforts and management. The funding organization did not participate in analysis and interpretation of the data for the manuscript but reviewed it for accuracy of the CHS study description and acknowledgment of study investigators and funding sources.
Additional Information: A full list of participating CHS investigators and institutions can be found at http://www.chs-nhlbi.org.
Acknowledgment: We thank Jill Perrotte, BS, and Kate M. Durda, Department of Pathology, University of Vermont College of Medicine, Burlington, for assistance with genotyping, and Jean Yee, BS, Department of Epidemiology, University of Washington, Seattle, for data management. All received support from the National Heart, Lung, and Blood Institute grant R01 HL-071862-03 for their contribution.