|Year : 2018 | Volume
| Issue : 21 | Page : 2630-2633
Plasma C-Reactive Protein and Abdominal Aortic Aneurysm: A Mendelian Randomization Analysis
Xue-Ying Qin1, Juan Juan1, Xiao Xiang1, Ying-Qi Wei2, Shang-Wei Zuo3, Tao Huang1, Da-Fang Chen1, Roger Marshall4, Jiang Xiong3, Wei Guo3, Yong-Hua Hu1
1 Department of Epidemiology and Biostatistics, School of Public Health, Peking University Health Science Center, Beijing 100191, China
2 Beijing Center for Disease Prevention and Control, Beijing 100013, China
3 Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing 100853, China
4 Section of Epidemiology and Biostatistics, School of Population Health, University of Auckland, Private Bag 92019, Auckland, New Zealand
|Date of Submission||25-May-2018|
|Date of Web Publication||25-Oct-2018|
Prof. Yong-Hua Hu
Department of Epidemiology and Biostatistics, School of Public Health, Peking University Health Science Center, 38 Xueyuan Road, Haidian District, Beijing 100191
Source of Support: None, Conflict of Interest: None
|How to cite this article:|
Qin XY, Juan J, Xiang X, Wei YQ, Zuo SW, Huang T, Chen DF, Marshall R, Xiong J, Guo W, Hu YH. Plasma C-Reactive Protein and Abdominal Aortic Aneurysm: A Mendelian Randomization Analysis. Chin Med J 2018;131:2630-3
|How to cite this URL:|
Qin XY, Juan J, Xiang X, Wei YQ, Zuo SW, Huang T, Chen DF, Marshall R, Xiong J, Guo W, Hu YH. Plasma C-Reactive Protein and Abdominal Aortic Aneurysm: A Mendelian Randomization Analysis. Chin Med J [serial online] 2018 [cited 2018 Nov 15];131:2630-3. Available from: http://www.cmj.org/text.asp?2018/131/21/2630/244124
To the Editor: Abdominal aortic aneurysm (AAA), a disease defined as a 50% increase in the aortic diameter or the diameter larger than 30 mm compared to the normal abdominal aortic, is a severe and chronic disorder occurring mostly among men older than 65 years of age. Understanding the risk factors and the biomarkers for AAA formation is imperative for creating targeted interventions to prevent, screen, or slow the disease progression.
Inflammatory processes are important in the formation of AAA. C-reactive protein (CRP) is an acute phase reactant produced by the liver that has long been used as a marker of systemic inflammation. Compelling epidemiological studies have suggested CRP is associated with the development of AAA. However, it is unclear whether this association is causal or because of confounding and reverse causality existing in observational studies., Moreover, limited randomized trials have provided data that can precisely demonstrate the underlying causal relationships between CRP and AAA.
Mendelian randomization (MR) analysis, a natural randomized trial, has been widely used to study the causal relevance of emerging risk factors with diseases. Genetic variants which are randomly allocated at birth, are used as instrumental variables in MR. In this study, using established genetic variants that determine CRP level, we performed an MR to investigate the causal role of CRP in AAA.
The study included 465 Chinese participants aged 40 years and above from an age- and gender-matched case–control study and details of the design were described previously. AAA patients (n = 155) were diagnosed at the Vascular and Endovascular Surgery Department of the Chinese PLA General Hospital. The patients were diagnosed by abdominal Doppler ultrasound or computed tomography. Patients with inflammatory AAA, mental disorders, or pregnant were excluded from the study. The control group came from two sources, control Group 1 (n = 155) were selected from the same hospital and control Group 2 (n = 155) were selected from communities in Fangshan District in Beijing. The control participants were age- and gender-matched with AAA patients but were negative for AAA by abdominal Doppler ultrasound [Figure 1].
|Figure 1: The flow chart of the study design. *rs1130864 was excluded from the data analysis for its high linkage disequilibrium with rs1417938; †rs1417938 was not used as an instrumental variable individually because its extremely low F-statistic in the first stage of MR analysis, but it was included into the genetic score. MR: Mendelian randomization.|
Click here to view
All participants were interviewed in person using a standardized questionnaire to collect data related to sociodemographics, history of chronic diseases and lifestyle factors, including smoking status, alcohol consumption, and physical activity. Participants underwent a standard physical examination including height, weight, systolic, and diastolic blood pressure and abdominal ultrasound or computed tomography. Peripheral venous fasting blood samples were obtained to detect biochemical parameters. Biochemical parameters were analyzed at the clinical laboratory of the Chinese PLA Hospital.
Four common single-nucleotide polymorphisms (SNPs) (rs1205, rs1130864, rs1800947, and rs1417938) on this gene were selected from the International HapMap Project collection of Han Chinese data and because they were reported to be significantly related to CRP levels in a Chinese population. Genomic DNA was extracted from the peripheral blood using the salt fractionation method. Genotypes were analyzed with Sequenom's MassARRAY® genotyping platform. All SNPs passed the quality control with genotyping call rate of >95%. We obtained the linkage disequilibrium estimates using data from SNP Annotation and Proxy Search, and rs1130864 was excluded from the analysis for its high linkage disequilibrium with rs1417938 (r2 = 0.962) [Table S1 [Additional file 1]].
Both separate SNPs and the constructed genetic score were used as the genetic instruments in our study. Instrumental variables were selected according to its relationship with CRP level and the F-statistics in the first stage regression of MR analysis; therefore, rs1800947, rs1205, and genetic score were involved for their F-statistics larger than 10 or approximated to 10. rs1417938 was not used as an instrumental variable individually because its F-statistic was extremely low, but it was included into the genetic score to maximize power and avoid the weak instrument bias in MR analysis [Figure 1] and [Table S2 [Additional file 2]].
Before calculating the genetic score, we first defined the exposure allele (or the risk allele) as the allele related to elevate circulating CRP levels according to our data and previous report, that is, T allele in rs1417938, G allele in rs1800947, and C allele in rs1205. We then calculated an unweighted CRP genetic score by adding the number of CRP-increasing alleles that the person had inherited at each variant. We treated genetic score as instrumental variable in two ways – one way was to do the median split to obtain a binary variable (below or above the median genetic score), and the other way was to use the original values of the genetic score.
All statistical analysis was undertaken using STATA software package (version 13.0, StataCorp, College Station, Texas, USA). Power analysis was undertaken using PASS software (version 14.0, NCSS, LLC. Kaysville, Utah, USA, ncss.com/software/pass). Hardy–Weinberg Equilibrium was tested in control group using the Chi-square test. Concentrations of CRP were natural log transformed for its skewed distribution (LnCRP). Variables are expressed as numbers and percentages (categorical variables), mean ± standard deviation (continuous variables with normal distribution), median and interquartile range for continuous variables with abnormal distribution. First, we compared the basic characteristics of the participants as divided by genetic variants and genetic score, respectively, to assess whether the potential confounding factors differed across the instrumental variables. Differences between the groups were analyzed using t-test and Kruskal–Wallis rank test for continuous variables and Chi-square tests for categorical variables. Second, we measured the difference of CRP level for each instrumental variable using linear regression. Third, we used two-stage least squares regression for the instrumental analysis. The first stage was to do the linear regression of Ln-transferred CRP levels on the genetic instruments, obtaining the predicted values for CRP levels. The second stage was to do the logistic regression of AAA on the predicted CRP levels. The standard errors were corrected by bootstrapping both first and second stage together 1000 times. Fourth, we compared the associations derived from the instrumental variable analysis to the observational analysis by conducting multivariable analysis. In this step, the original genetic score was used as instrumental variable. Potential confounding factors were adjusted gradually, including age and gender, lifestyle factors including smoking, drinking, and physical activity, history of diseases, and body mass index (BMI). P = 0.05 was considered as the statistical significance threshold.
[Table 1] shows the basic characteristics of the study participants. There were significant differences for smoking, drinking, history of hypertension and hyperlipidemia, and the concentration of plasma CRP between cases and controls, and no differences were found for age, gender, physical activity, history of diabetes, BMI, and genotypes of rs1417938, rs1800947, and rs1205.
All the basic characteristics were comparable when we grouped participants based on the genotypes of rs1800947 and rs1205, and the genetic score, respectively, thus indicating that these genetic variants are reliable instrumental variables. Only in rs1417938, we found that age and BMI were a significant difference in AT and TT groups, but did not find the significant difference in CRP [Table S3 [Additional file 3]],[Table S4 [Additional file 4]],[Table S5 [Additional file 5]],[Table S6 [Additional file 6]]. All the SNPs were in Hardy–Weinberg equilibrium (P > 0.01), and the allele frequencies in the control group from the community were shown in [Table S7 [Additional file 7]].
[Figure 2] shows the genetic association with LnCRP and estimated causal relationship between CRP and AAA. We found that G allele of rs1800947 (0.782 unit per G allele with 95% confidence interval[CI], 0.204–1.360, P = 0.008), C allele of rs1205 (0.358 unit per C allele with 95% CI, 0.171–0.545, P < 0.001), higher genetic score (0.423 unit per allele with 95% CI: 0.247–0.599, P < 0.001) were significantly associated with elevated CRP concentrations. MR analyses showed that genetically elevated CRP level were not significantly associated with higher risk of AAA (odds ratio [OR] = 1.163 for rs1800947, 95% CI, 0.174–7.748, P = 0.876; OR = 0.767 for rs1205, 95% CI, 0.249–2.367, P = 0.645; OR = 1.130 for binary genetic score, 95% CI, 0.524–2.439, P = 0.755; OR = 0.196 for continuous genetic score, 95% CI, 0.605–2.364, P = 0.606).
|Figure 2: Effect of genetically elevated CRP levels on abdominal aortic aneurysm using MR approach in a Chinese population. *Mean difference in LnCRP was assessed using regression model with genetic variants and genetic score as repressors, respectively; †OR was for the association and was estimated from the instrumental variable analysis by using genetic variants and genetic score as instrumental variables, respectively. No other covariates were adjusted in the MR model. CRP: C-reactive protein; MR: Mendelian randomization; OR: Odds ratio.|
Click here to view
[Table 2] compares the conventional observational multivariable results and the MR estimates for the association between CRP and AAA in three models. In model 1, age and gender were adjusted, and we observed that observational analysis obtained significant directionally association of CRP with risk of AAA (OR = 1.357, 95% CI: 1.163–1.583, P < 0.001), but instrumental variable analysis found nonsignificant association (OR = 1.203, 95% CI, 0.574–2.519, P = 0.625). When further adjusted for smoking, drinking, and physical activity in model 2, we observed significant associations in observational analysis for CRP and AAA, but not in instrumental variable analysis (OR in observational analysis, 1.276, 95% CI, 1.085–1.502, P = 0.003; OR in instrumental variable analysis, 1.291, 95% CI, 0.606–2.754, P = 0.508). In model 3, BMI and history of chronic diseases were additionally adjusted, and similar results were observed for the association between CRP and AAA as in model 2.
|Table 2: The effect size between LnCRP and AAA estimated from the conventional multivariate model and instrumental variable analysis|
Click here to view
Our findings showed that the genetically elevated CRP concentrations were not associated with AAA, indicating that this association was not causal.
By using genetic variants as the instrumental variables, this approach could provide evidence of causality between exposure and outcome. To be effective and strong instrumental variables, the genetic variants must satisfy these assumptions: (1) the genetic variants should be related to the exposure we studied; (2) the instrument should not be associated with confounding factors of the exposure and outcome association; and (3) the genetic variants must affect the outcome only through the exposure pathway, and no other pathways. Both separate CRP genetic variants (except rs1417938) and genetic score were found to be related significantly to the plasma CRP levels, and the F-statistics were reasonably large for these variants and genetic score, which satisfied the first assumption of the instrumental variables. Our results showed that the basic characteristics were equally distributed across each instrumental variable so that the second assumption was satisfied according to the current data although they could not represent all the confounders. For the third assumption, the functions of the genetic variants we chose located in CRP gene and their relationships with CRP level have been widely studied. In addition, the genetic score, we constructed could combine information from multiple variants that associated with the CRP, and could maximize power and avoid the weak instrument bias.
Our finding did not find plasma CRP was a determinant of the AAA. This was inconsistent with the results of most observational studies, including one study based on the same population used in the present study., Atherosclerosis and inflammation, with increased levels of related markers, such as CRP, are considered the common pathophysiologic process for vascular-related disease. However, some studies in vivo did not find the significant association between plasma CRP and cardiovascular events, which agreed with our finding. In the pathway of AAA, we suppose that plasma CRP may be an early predictor of AAA, but not a causal factor of the AAA development; and therefore, the observed association may derive from reverse causality. Irrespective of the causal effect of plasma CRP on AAA, however, there are considerable evidence persisting that inflammatory activation and degradation of the elastic media are the important pathophysiologic process in the development of AAA. More evidence is needed to identify specific genetic, environmental, and biochemical factors for AAA development.
There were several limitations in our study. One limitation is the low power of the study mainly led by small sample size and the lack of enough common genetic variants in CRP genes to construct a complete genetic score. The concentration of CRP in the circulation is influenced by many polymorphisms in CRP gene and other genes, more power would be obtained if more genetic variants could be involved. Therefore, larger sample size and more genotyping data are expected to confirm the CRP and AAA through MR studies. Moreover, our study was conducted in Chinese population; therefore, the results cannot be generalized to other ethnic populations.
Supplementary information is linked to the online version of the paper on the Chinese Medical Journal website.
We would like to thank Professor Tao Wu from School of Public Health in Peking University for her advices and thank all participants in the study.
Financial support and sponsorship
The study was supported by the National Natural Science Foundation of China (81230066, 81502874) and The Eleventh Five-year Plan in Health Care Foundation of PLA (09BJZ04).
Conflicts of interest
There are no conflicts of interest.
| References|| |
Bird AN, Davis AM. Screening for abdominal aortic aneurysm. JAMA 2015;313:1156-7. doi: 10.1001/jama.2015.0996.
Folsom AR, Yao L, Alonso A, Lutsey PL, Missov E, Lederle FA, et al.
Circulating biomarkers and abdominal aortic aneurysm incidence: The atherosclerosis risk in communities (ARIC) study. Circulation 2015;132:578-85. doi: 10.1161/circulationaha.115.016537.
Stather PW, Sidloff DA, Dattani N, Gokani VJ, Choke E, Sayers RD, et al.
Meta-analysis and meta-regression analysis of biomarkers for abdominal aortic aneurysm. Br J Surg 2014;101:1358-72. doi: 10.1002/bjs.9593.
Mojtahedzadeh M, Chelkeba L, Ranjvar-Shahrivar M, Najafi A, Moini M, Najmeddin F, et al
. Randomized trial of the effect of magnesium sulfate continuous infusion on IL-6 and CRP serum levels following abdominal aortic aneurysm surgery. Iran J Pharm Res 2016;15:951-6.
Burgess S, Timpson NJ, Ebrahim S, Davey Smith G. Mendelian randomization: Where are we now and where are we going? Int J Epidemiol 2015;44:379-88. doi: 10.1093/ije/dyv108.
Wei Y, Xiong J, Zuo S, Chen F, Chen D, Wu T, et al.
Association of polymorphisms on chromosome 9p21.3 region with increased susceptibility of abdominal aortic aneurysm in a Chinese Han population. J Vasc Surg 2014;59:879-85. doi: 10.1016/j.jvs. 2013.10.095.
Kong H, Qian YS, Tang XF, Zhang J, Gao PJ, Zhang Y, et al.
C-reactive protein (CRP) gene polymorphisms, CRP levels and risk of incident essential hypertension: Findings from an observational cohort of Han Chinese. Hypertens Res 2012;35:1019-23. doi: 10.1038/hr.2012.89.
Bradić M, Costa J, Chelo IM. Genotyping with Sequenom. Methods Mol Biol 2011;772:193-210. doi: 10.1007/978-1-61779-228-1_11.
Lee CC, You NC, Song Y, Hsu YH, Manson J, Nathan L, et al.
Relation of genetic variation in the gene coding for C-reactive protein with its plasma protein concentrations: Findings from the women's health initiative observational cohort. Clin Chem 2009;55:351-60. doi: 10.1373/clinchem.2008.117176.
Crawford DC, Sanders CL, Qin X, Smith JD, Shephard C, Wong M, et al.
Genetic variation is associated with C-reactive protein levels in the third national health and nutrition examination survey. Circulation 2006;114:2458-65. doi: 10.1161/circulationaha.106.615740.
Hopewell JC, Clarke R. Emerging risk factors for stroke: What have we learned from Mendelian randomization studies? Stroke 2016;47:1673-8. doi: 10.1161/strokeaha.115.010646.
Shangwei Z, Yingqi W, Jiang X, Zhongyin W, Juan J, Dafang C, et al.
Serum high-sensitive C-reactive protein level and CRP genetic polymorphisms are associated with abdominal aortic aneurysm. Ann Vasc Surg 2017;45:186-92. doi: 10.1016/j.avsg.2017.05.024.
Pai JK, Pischon T, Ma J, Manson JE, Hankinson SE, Joshipura K, et al.
Inflammatory markers and the risk of coronary heart disease in men and women. N
Engl J Med 2004;351:2599-610. doi: 10.1056/NEJMoa040967.
Koike T, Kitajima S, Yu Y, Nishijima K, Zhang J, Ozaki Y, et al.
Human C-reactive protein does not promote atherosclerosis in transgenic rabbits. Circulation 2009;120:2088-94. doi: 10.1161/circulationaha.109.872796.
Dehghan A, Dupuis J, Barbalic M, Bis JC, Eiriksdottir G, Lu C, et al.
Meta-analysis of genome-wide association studies in >80 000 subjects identifies multiple loci for C-reactive protein levels. Circulation 2011;123:731-8. doi: 10.1161/circulationaha.110.948570.
[Figure 1], [Figure 2]
[Table 1], [Table 2]