Skip to main content
  • Research article
  • Open access
  • Published:

Association of high sensitive C-reactive protein with coronary heart disease: a Mendelian randomization study

Abstract

Objectives

Whether high sensitivity C-reactive protein (hs-CRP) has a causal effect on coronary heart disease (CHD) is unclear. This study investigated the causal effect of hs-CRP on CHD risk using Mendelian Randomization (MR) analysis.

Methods

A total of 3802 subjects were recruited in the follow-up study. Linear regression model was used to evaluate the relationship between CRP polymorphisms and hs-CRP. Survival receiver operator characteristic curve method was used to explore the cut-off of hs-CRP on CHD incidence. Cox regression model was applied to detect the association of hs-CRP with CHD by calculating the hazard ratio (HR) and 95% confidence interval (CI). Rs1205 and rs876537 in CRP were selected as instrumental variables in MR analysis.

Results

During a median follow-up time of 5.01 years, 98 CHD incidence was identified (47.03/104 person-years). Hs-CRP was significantly increased among rs1205 and rs876537 genotypes with r values of 0.064 and 0.066, respectively. Hs-CRP 1.08 mg/L was identified as the cut-off value with a maximum value of sensitivity and specificity on prediction of CHD. Participants with ≥1.08 mg/L of hs-CRP has a higher risk of CHD incidence than that of participants with < 1.08 mg/L, the adjusted HR (95% CI) was 1.69 (1.11–2.60) with a P value of 0.016. No significant casual association was observed between hs-CRP and CHD with a P value of 0.777.

Conclusions

The association between hs-CRP and CHD is unlikely to be causal, hs-CRP might be a predictor for incidence of CHD in general population.

Peer Review reports

Background

Coronary heart disease (CHD) is one of the leading global cause of death and disability. The current situation of rapid aging in China makes CHD become the 4th leading cause of mortality [1]. The increased adiposity, type 2 diabetes mellitus (T2DM), hypertension, sedentary lifestyle, as well as genetic factors, are the well-known determinants of CHD [2].

Inflammation promotes endothelial cell damage and atherogenesis, is an important risk factor for CHD in clinical [3]. As a famous parameter of inflammation, C-reactive protein (CRP) is also considered a valuable predictor of CHD risk [4]. To meet the super-precision of clinical practice, the measurement of high sensitivity C-reactive protein (hs-CRP) is proposed for its better capacity in identifying the individuals at high risk for CHD [5]. Hs-CRP has commonly been used for cardiovascular diseases (CVDs) risk stratification. The elevated level of serum hs-CRP was significantly increased the risk of myocardial infarction [6]. Randomized clinical trials have observed that the use of statins in individuals with elevated hs-CRP was associated with a reduction in hs-CRP and a decreased risk of vascular events [7]. Even numerous of studies support that the hs-CRP is an established marker for future risk of CVDs, however, inconsistent findings have been reported on CHD [8, 9]. Therefore, it has considerable interest in establishing whether hs-CRP has a causal role in CHD.

Mendelian randomization (MR) studies utilize genetic variants (such as he-CRP) as instrumental variables (IVs) to investigate possible causal relationship between exposure and outcomes through an intermediate trait [10]. The dominant advantage of MR studies is that they are not as vulnerable to confounding and reverse causality [11, 12]. If the intermediate trait is causally linked to disease, then genetic variants influencing the trait should also influence disease risk. Indeed, MR has been regarded as nature’s analogue of randomized controlled trials, thus, it has also been used in cardiovascular research to explore novel potential etiologic mechanisms and enhance our understanding of current therapies [13, 14].

The aim of the current study was to identify common single nucleotide polymorphisms (SNPs) in CRP that influence hs-CRP levels, and use the concept of MR to improve understanding of the possible causal relationship of hs-CRP levels with CHD.

Methods

Study population

A total of 4222 participants were drawn from the community hypertension survey in Yixing city, China, 2009 [15]. After we excluded the individuals with a history of CHD (n = 50), elder participants (n = 94) and missing measurement of hs-CRP (n = 373), 3802 subjects were finally enrolled in the follow-up study.

To access the incidence of disease status in the follow-up study, both concentrated and household surveys were conducted. Data collected from the local hospitals, Centers for Disease Control and Prevention, community health service centers and social security center was further inspected to reduce the information bias. During a median follow-up time of 5.01 years, 98 new-onset CHDs were recorded.

All participants were interviewed and underwent physical examinations and laboratory tests. The demographic characteristics of the participants were obtained by trained research staff though a standard questionnaire. The weight, height and 3 blood pressure measurements were obtained from each participant by trained and certified observers according to a standard protocol.

Written informed consent was obtained from all participants. The research protocol was approved by the ethics committee of Nanjing Medical University and written informed consent was obtained from all subjects during epidemiological interviews.

SNP selection

The CRP gene was located on chromosome 1q23.2 (Gene ID: 1401; NC_000001.11) and spans 2.3 kbps and contains 4 exons. We searched the SNPs covered CRP gene from the upstream 5 kb to the downstream 2 kb and selected tagging SNPs (tagSNPs) from the database of the Chinese Han population in Beijing, China of the International Hap MAP Project (HapMap Data Rel 24/phase II Nov08, on NCBI B36 assembly, dbSNPb126). All tagSNPs were selected with a minor allele frequency (MAF) ≥ 0.05 and linkage disequilibrium (LD) r2 ≥ 0.8. We also applied a functional candidate strategy to select potential functional SNPs on the bioinformatics effect prediction website (SNPINFO, https://snpinfo.niehs.nih.gov/). Finally, rs1205(C > T), rs1073715(C > T), rs876537(C > T) and rs2808630(C > T) were selected. The specific biological information was summarized in the Additional file 1: Table S1.

Blood sampling

Blood samples were taken from participants after an over-night fasting (> 10 h), adding ethylenediamine tetraacetic acid (EDTA)-containing receptacles. The blood samples were centrifuged at 3500 rpm to separate the serum and stored at − 20 °C. Serum hs-CRP, total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C) and glucose (GLU) were measured by an automatic biochemical analysis system of Siemens (Advia 1200, German).

Following recommendations from the American Centers for Disease Control and the American Heart Association [16], the participants with hs-CRP ≥1.0 mg/L and hs-CRP < 1.0 mg/L correspond to high risk group (HRG) and normal group (NG), respectively.

SNP genotyping

DNA was extracted using a standard phenole-chloroform method. Genotyping was performed using the TaqMan allelic discrimination assay in 384-well plates on the platform of 7900HT Real-time polymerase chain reaction (PCR) System (Applied Biosystems, Foster City, CA). The primers and probes were designed using Primer Express Oligo Design software v2.0 (ABI PRISM). Genotyping results were determined using SDS 2.3 Allelic Discrimination Software (Applied Biosystems). Meanwhile, each plate was included blank samples as negative controls for the genotyping quality confirmation. The successful call rates of SNPs rs10737175, rs1205, rs2808630 and rs2966449 were 98.44, 98.31, 98.39 and 98.41%, respectively.

Statistical analysis

Unpaired Student’s t-test was used to test the differences in of all the quantitative variables among groups presented as means ± standard deviation (SD). Chi-square (χ2) test was performed to compare the proportion of hypertension, T2DM, smoking and drinking between HRG and NG group. Linear regression model was used to analyze the correlation between the variations of CRP gene and hs-CRP. Besides, we used survival receiver operator characteristic curve (SROC) to explore the cut-off values of hs-CRP with CHD incidence. Cox proportional hazard regression model was applied to evaluate the association of CRP gene with CHD by calculating hazard ratio (HR) and 95% confidence interval (CI).

The MR analysis for exploring the causal effect of hs-CRP on the risk of CHD applying inverse-variance weighted (IVW) methods, where the SNPs are deemed as IVs. The basic principle of MR analysis for the current study was shown in Fig. 1. Effects are given as odds ratios (OR) and 95%CI. A two-tailed P value of 0.05 was defined as the cut-off for statistical significance. The SROC and MR analysis were undertaken using the R Packages of “survivalROC” and “Mendelian Randomization”.

Fig. 1
figure 1

The basic principle of Mendelian Randomization analysis for the current study. BMI, body mass index; CRP, C-reactive protein; CHD, coronary heart disease; IVs, instrumental variables; SNP, single nucleotide polymorphisms

Results

Demographic characteristics at baseline

The demographic characteristics of participants were summarized in Table 1. The age, BMI, SBP, DBP, TC, TG, LDL-C, GLU and proportion of hypertension in HRS were significantly higher than NG, whereas HDL-C, the proportion of T2DM, smoking and drinking were significantly lower (P < 0.05). Overall, these characteristics were adjusted as confounding factors in the multivariate regression analysis.

Table 1 Baseline characteristics for different hs-CRP groups

Association of CRP polymorphisms with CHD incidence

During a median follow-up time of 5.01 years, 98 CHD incidence was identified with an incidence density of 47.03/104 person-years. No significant associations were observed between each of the four SNPs with CHD in the whole population. After adjusted for age, gender, BMI, HDL-C, LDL-C, TC, TG, smoking, drinking, hypertension and T2DM, the results were still not significant (Table 2). Stratification analysis by age (55 years), gender, smoking and drinking status were further conducted, however, no significant associations were found between CRP and CHD (Additional file 1: Table S2).

Table 2 Association analysis of CRP polymorphisms with CHD incidence

Association of hs-CRP with CHD

The incidence density of CHD for HRG and NG group were 71.08/104 and 33.79/104 person-years, respectively. The elevated hs-CRP was significantly increased CHD risk, the HR (95% CI) was 2.05 (1.36–3.09) with a P value of 0.001. But, the result was still not significant after adjusted for cofounders, the adjusted HR (95% CI) was 1.54 (0.99–2.36) with a P value of 0.051.

The SROC analysis showed that the highest sensitivity and specificity of cut-off point for hs-CRP to predict CHD incidence was 1.08 mg/L, the sensitivity and specificity were 60.2 and 59.1%, respectively (Additional file 1: Figure S1). Participants with ≥1.08 mg/L of hs-CRP has a higher risk of CHD incidence than that of participants with < 1.08 mg/L (74.92/104 vs 32.97/104 person-years), the adjusted HR (95% CI) was 1.69 (1.11–2.60) with a P value of 0.016. Detailed results were displayed in Table 3.

Table 3 Association analysis of hs-CRP with CHD in the follow-up study

Correlation of CRP variation with hs-CRP

CRP rs1205 and rs876537 C > T variations were significantly correlated with hs-CRP, after adjusted for age, gender, BMI, TC, TG, HDL-C, LDL-C, smoking, drinking, T2DM and hypertension, the standard regression coefficients were 0.064 and 0.066, respectively, P < 0.001. No significant correlations were observed between rs10737175, rs2808630 and hs-CRP, with P values of 0.713 and 0.271 (Additional file 1: Table S3). Thus, we selected rs1205 and rs876537 were selected as IVs in the MR analysis.

MR analysis between hs-CRP and CHD

The causal effect was evaluated by using MR analysis assessing the association of genetically predicted hs-CRP with CHD risk (Table 4). The crude model showed that hs-CRP was not statistically associated with CHD, the cruder OR (95% CI) was 1.04 (0.67–1.62) with a P value of 0.855. The result remained negative after adjustment for age, gender, BMI, TC, TG, HDL-C, LDL-C, smoking, drinking, T2DM and hypertension, the adjust OR (95% CI) was 1.07 (0.69–2.18) with a P value of 0.777.

Table 4 The Mendelian randomization analysis of hs-CRP with CHD

Discussion

The present study firstly applied MR design to explore the casual effect of hs-CRP and CHD risk in Chinese Han population. The main results showed that CRP rs1205 and rs876537 were correlated with serum hs-CRP level, but no causal effect of elevated hs-CRP on CHD risk was observed. Although the association between hs-CRP and CHD is unlikely to be causal, hs-CRP might be a predictor for incidence of CHD in general population.

Epidemiological studies have shown that plasma levels of hs-CRP are a strong independent predictor of risk of future CVDs [17, 18]. As a well-known inflammatory biomarker, hs-CRP is significantly elevated in patients dying suddenly with severe coronary artery diseases, both with and without acute coronary thrombosis [19]. There is considerable interest in establishing whether hs-CRP has a causal role in CHD. MR study, an IV-based method to infer the causality between intermediate phenotypes and disease, has been widely conducted in cardiovascular disease research [20]. Genetic variants that are associated with intermediate phenotypes are introduced as IVs in MR to estimate the effect of phenotypic exposures on disease outcome.

Previous large sample size studies showed that no association of variants in the CRP locus and CHD were found, arguing against a causal role for CRP in atherosclerosis [21]. Another large-scale MR research selected CRP rs1205 as one of IVs and indicated that CRP was unlikely to be even a modest causal factor in CHD [13]. As a non-specific inflammatory biomarker of CRP, hs-CRP takes the advantages of a higher sensitivity, accuracy and reproducibility. Our novel results firstly indicate that genetically raised hs-CRP are unrelated to the risk of CHD, which are consistent with previous studies [20, 21]. These findings imply that the cardiologists should treat CRP or hs-CRP as genuine confounders and pay more attention to other inflammatory markers in high-risk individuals.

Although American Centers for Disease Control and the American Heart Association has recommended hs-CRP 1.0 mg/L as a cut-off value to identify individuals who is under a high risk of inflammatory reaction [16], no study has proved 1.0 mg/L is applicable to Chinese population. The SROC method was used to explore the hs-CRP cut-off value when evaluating the association of CHD. Our result showed hs-CRP 1.08 mg/L is the cut-off value with a relative high sensitivity and specificity to predict CHD incidence, inferring the criterion of 1.0 mg/L is still suitable in China.

The current study investigated the cause effect of hs-CRP and CHD risk in a prospective cohort study. The strengths of this study should be noted. First, for the expensive cost of hs-CRP measurement, limited MR studies have carried out to evaluate the association of hs-CRP with CHD in a prospective cohort study. All the participants with a measurement of hs-CRP was enrolled in the follow-up study, the relative large sample size in our study substantially enhanced the statistical power, providing an objective and credible result on he-CRP and CHD. Second, it was firstly to explore the cut-off value of hs-CRP on predicting CHD incidence in Chinese Han population.

Several limitations were existed in the current study. First of all, the common SNPs that significantly correlated with hs-CRP are not derived from genome-wide association study data. In addition, the causal effect of hs-CRP on sub-types of CHD were not investigated. Besides, the pleiotropy effect of SNPs in MR analysis was not evaluated. Last but not least, during acute events, the CRP gene expression can be up-regulated, leading to cascades that would increase the hs-CRP levels. However, our data might not accurately represent the fundamental of gene expression related to hs-CRP.

Conclusion

In conclusion, CRP rs1205 and rs876537 variation were correlated with hs-CRP level.

Although the association between hs-CRP and CHD is unlikely to be causal, hs-CRP might be a predictor for incidence of CHD in general population.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CHD:

Coronary heart disease

CI:

Confidence interval

CVDs:

Cardiovascular diseases

DL-C:

High-density lipoprotein cholesterol

GLU:

Glucose

HR:

Hazard ratio

HRG:

High risk group

Hs-CRP:

High-sensitivity C-reactive protein

IV:

Instrumental variable

LDL-C:

Low-density lipoprotein cholesterol

MAF:

Minor allele frequency

MR:

Mendelian Randomization

NG:

Normal group

OR:

Odds Ratios

PCR:

Polymerase chain reaction

SD:

Standard deviation

SNP:

Single nucleotide polymorphism

SROC:

Survival receiver operator characteristic curve

T2DM:

Type 2 diabetes mellitus

TC:

Total cholesterol

TG:

Triglycerides

References

  1. He J, et al. Major causes of death among men and women in China. N Engl J Med. 2005;353(11):1124–34.

    Article  CAS  Google Scholar 

  2. Ripatti S, et al. A multilocus genetic risk score for coronary heart disease: case-control and prospective cohort analyses. Lancet. 2010;376(9750):1393–400.

    Article  Google Scholar 

  3. Ross R. Atherosclerosis--an inflammatory disease. N Engl J Med. 1999;340(2):115–26.

    Article  CAS  Google Scholar 

  4. Buckley DI, et al. C-reactive protein as a risk factor for coronary heart disease: a systematic review and meta-analyses for the U.S. preventive services task force. Ann Intern Med. 2009;151(7):483–95.

    Article  Google Scholar 

  5. Mirhafez SR, et al. Serum high-sensitivity C-reactive protein as a biomarker in patients with metabolic syndrome: evidence-based study with 7284 subjects. Eur J Clin Nutr. 2016;70(11):1298–304.

    Article  CAS  Google Scholar 

  6. Sugiyama T, et al. Relationship between serum high-sensitivity C-reactive protein and myocardial infarction in a general Japanese population. J Clin Lab Anal. 2016;30(6):999–1002.

    Article  CAS  Google Scholar 

  7. Ridker PM, et al. Rosuvastatin to prevent vascular events in men and women with elevated C-reactive protein. N Engl J Med. 2008;359(21):2195–207.

    Article  CAS  Google Scholar 

  8. Pepys MB. C-reactive protein is neither a marker nor a mediator of atherosclerosis. Nat Clin Pract Nephrol. 2008;4(5):234–5.

    Article  Google Scholar 

  9. Zacho J, et al. Genetically elevated C-reactive protein and ischemic vascular disease. N Engl J Med. 2008;359(18):1897–908.

    Article  CAS  Google Scholar 

  10. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  Google Scholar 

  11. Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.

    Article  Google Scholar 

  12. Hartwig FP, et al. Two-sample Mendelian randomization: avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int J Epidemiol. 2016;45(6):1717–26.

    Article  Google Scholar 

  13. Collaboration, C.R.P.C.H.D.G, et al. Association between C reactive protein and coronary heart disease: mendelian randomisation analysis based on individual participant data. BMJ. 2011;342:d548.

    Article  Google Scholar 

  14. Swerdlow DI, et al. HMG-coenzyme a reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials. Lancet. 2015;385(9965):351–61.

    Article  CAS  Google Scholar 

  15. Chen J, et al. Common variants in TGFBR2 and miR-518 genes are associated with hypertension in the Chinese population. Am J Hypertens. 2014;27(10):1268–76.

    Article  Google Scholar 

  16. Ridker PM. Inflammation in atherothrombosis: how to use high-sensitivity C-reactive protein (hsCRP) in clinical practice. Am Heart Hosp J. 2004;2(4 Suppl 1):4–9.

    PubMed  Google Scholar 

  17. Ridker PM, et al. C-reactive protein and other markers of inflammation in the prediction of cardiovascular disease in women. N Engl J Med. 2000;342(12):836–43.

    Article  CAS  Google Scholar 

  18. Tracy RP, et al. Relationship of C-reactive protein to risk of cardiovascular disease in the elderly. Results from the cardiovascular health study and the rural health promotion project. Arterioscler Thromb Vasc Biol. 1997;17(6):1121–7.

    Article  CAS  Google Scholar 

  19. Burke AP, et al. Elevated C-reactive protein values and atherosclerosis in sudden coronary death: association with different pathologies. Circulation. 2002;105(17):2019–23.

    Article  CAS  Google Scholar 

  20. Kawashiri MA, et al. Mendelian randomization: its impact on cardiovascular disease. J Cardiol. 2018;72(4):307–13.

    Article  Google Scholar 

  21. Elliott P, et al. Genetic loci associated with C-reactive protein levels and risk of coronary heart disease. JAMA. 2009;302(1):37–48.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

First and foremost, I would like to show my deepest gratitude to Dr. Song Yang and Prof. Chong Shen, who are respectable, responsible and resourceful scholars, providing me with valuable guidance in every stage of the writing of this thesis. Their keen and vigorous academic observation enlightens me not only in this thesis but also in my future study. I shall extend my thanks to Mr. Shen for all his kindness and help. My sincere appreciation also goes to my colleagues from Department of Cardiology, Affiliated Yixing People’s Hospital of Jiangsu University, who participated this study with great cooperation.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No.81573232), The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

SY and CS conceived and designed the study. Clinical data was collected and collated by QZ, YC, XZ, JS, PW. XC and YJ analyzed the data. QZ wrote the paper. SY and CS reviewed and edited the manuscript. All authors read and approved the manuscript.

Corresponding author

Correspondence to Song Yang.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all participants. The research protocol was approved by the ethics committee of Nanjing Medical University and written informed consent was obtained from all subjects during epidemiological interviews.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Table S1. Bioinformatics Analysis of 4 selected SNPs in CRP Gene. Table S2. Stratification analysis of CRP gene polymorphisms and CHD. Table S3. Partial correlation analysis of CRP variation with hs-CRP. Figure S1. Survival ROC analysis for hs-CRP and CHD incidence.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhuang, Q., Shen, C., Chen, Y. et al. Association of high sensitive C-reactive protein with coronary heart disease: a Mendelian randomization study. BMC Med Genet 20, 170 (2019). https://doi.org/10.1186/s12881-019-0910-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12881-019-0910-z

Keywords