Original Research ARTICLE
Exposing the Causal Effect of C-Reactive Protein on the Risk of Type 2 Diabetes Mellitus: A Mendelian Randomization Study
- 1College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- 2Department of Radiology, The Second Affiliated Hospital of Harbin Medical University, Harbin, China
- 3Department of Radiology, Longhua Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 4Department of Rehabilitation and Pharmacy, Heilongjiang Province Land Reclamation Headquarters General Hospital, Harbin, China
As a biomarker of inflammation, C-reactive protein (CRP) has attracted much attention due to its role in the incidence of type 2 diabetes mellitus (T2DM). Prospective studies have observed a positive correlation between the level of serum CRP and the incidence of T2DM. Recently, studies have reported that drugs for curing T2DM can also decrease the level of serum CRP. However, it is not yet clear whether high CRP levels cause T2DM. To evaluate this, we conducted a Mendelian randomization (MR) analysis using genetic variations as instrumental variables (IVs). Significantly associated single nucleotide polymorphisms (SNPs) of CRP were obtained from a genome-wide study and a replication study. Therein, 17,967 participants were utilized for the genome-wide association study (GWAS), and another 14,747 participants were utilized for the replication of identifying SNPs associated with CRP levels. The associations between SNPs and T2DM were from the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) consortium. After removing SNPs in linkage disequilibrium (LD) and T2DM-related SNPs, the four remaining CRP-related SNPs were deemed as IVs. To evaluate the pooled influence of these IVs on the risk of developing T2DM through CRP, the penalized robust inverse-variance weighted (IVW) method was carried out. The combined result (OR 1.114048; 95% CI 1.058656 to 1.172338; P = 0.024) showed that high levels of CRP significantly increase the risk of T2DM. In the subsequent analysis of the relationship between CRP and type 1 diabetes mellitus (T1DM), the pooled result (OR 1.017145; 95% CI 0.9066489 to 1.14225; P = 0.909) supported that CRP levels cannot determine the risk of developing T1DM.
According to the 8th edition of the International Diabetes Federation (IDF) Diabetes Atlas, the number of Diabetes Mellitus (DM) patients has continued to increase across the globe. DM includes a group of chronic metabolic diseases involving hyperglycemia (Olokoba et al., 2012; Pan et al., 2013; Shi and Hu, 2014), which, over prolonged periods, leads to injury in various tissues of the body (Atar and Hanssen, 2017; Gentile et al., 2017). The corresponding abundant fatal complications occur for patients suffering from DM for over 10 years. The IDF Diabetes Atlas reported that 425 million DM patients are adults, two-thirds (327 million) of whom are of working age. Over 90% of these patients are type 2 DM (T2DM) patients. Thus, it is urgent for us to explore a way to reduce the number of individuals with T2DM (Zou et al., 2018).
Accumulated evidence shows that intermediate phenotypes, such as body mass index (BMI) (Zhang et al., 2017; Cheng et al., 2018a), systolic blood pressure (Svensson et al., 2017), circulating uric acid (Xu et al., 2016), etc., are responsible for the onset of this type of DM (Cheng et al., 2018c). To prevent T2DM, it would be valuable to influence these intermediate phenotypes. Although there have been recent successes, a large number of potential intermediate phenotypes causing T2DM need to be identified.
Since chronic low-grade inflammation is associated with T2DM (Kohn et al., 2005), one of its markers, C-reactive protein (CRP), is frequently investigated as an intermediate phenotype that increases the risk of T2DM. Doi et al. (2005) conducted an observational study for a mean of 9.0 years on 1,759 Japanese patients without diabetes. They observed that elevated CRP concentration is a significant predictor of diabetes, independent of obesity and insulin resistance. Analogous conclusions were also observed in American and European individuals (Freeman et al., 2002; Nakanishi et al., 2003). On the other hand, after adjusting for sex, adiposity, and insulin resistance, no association between CRP and the risk of T2DM was shown in other studies (Festa et al., 2002; Krakoff et al., 2003; Thorand et al., 2007).
A recent prospective analysis and meta-analysis study (Lee et al., 2009) aimed to elucidate the role of CRP in T2DM. The high association (OR 1.49; 95% CI 1.03–2.15, p = 0.03) between serum CRP and the risk of diabetes was observed in a Norfolk cohort including 293 diabetes cases and 708 controls. A consistent result (RR 1.72; 95% CI 1.54–1.92) was shown in the further meta-analysis of 16 published studies involving 3,920 diabetes cases and 24,914 controls. However, no significant association (OR 1.00; 95% CI 0.66–1.51, p = 1.0) was obtained after adjustment for Waist-Hip Ratio (WHR), serum γ-glutamyltransferase and adiponectin. In addition, the heterogeneity of these publications (I2 = 52.8%) limited the conclusion in the meta-analysis.
The main challenge for current prospective studies is that the potential confounding factors have a great effect on the observations (Cheng et al., 2017; Zeng et al., 2018). This lack of clarity has hindered researchers in determining the association between CRP and T2DM. In addition, without evidence from randomized controlled trials (RCTs), it is hard to be sure of the role of CRP in the risk of T2DM. To this end, the Mendelian randomization (MR) method (Lawlor et al., 2008), revealing the role of a risk factor in disease etiology, is utilized here. MR is an instrumental variable (IV)-based method to infer the causality between intermediate phenotypes and disease. Genetic variants that are associated with intermediate phenotypes are introduced as IVs in MR to estimate the effect of phenotypic exposures on disease outcome (Figure 1). Due to random distribution of gene variants during gametogenesis, the confounding effect can be extremely low.
Materials and Methods
Here, we introduced MR analysis for analyzing the causal effect of CRP on the risk of T2DM, where single nucleotide polymorphisms (SNPs) are deemed as IVs. The Mendelian randomization method is a method for testing whether a changeable exposure has a causal effect on the development of a disease. In this study, we applied a two-sample MR method. Through the explanation of the larger proportion of the exposure variance by general instrumental variables (IVs), this method can lead to more accurate causal effect estimation (Bowden et al., 2015; Hartwig et al., 2017). The basic principle of MR analysis for this purpose is shown in Figure 1. According to MR analysis, SNPs influence T2DM outcome through intermediate phenotypes of CRP. Therefore, the SNPs should be significantly associated with CRP but not with T2DM. At the same time, in the population of genes associated with the phenotype, there may be a linkage disequilibrium (LD) effect. Due to the non-random linkage of alleles among different loci in a given population, loci are thought to be in LD when the frequency of association of their different alleles is higher or lower than what would be expected if the loci were independent and associated randomly (Slatkin, 2008; Aissani, 2014). Thus, a notable issue is that the assessment would be biased for those SNPs in LD, since their effects would be expanded when combining all of those SNPs (Noyce et al., 2017).
Figure 2 shows the basic idea of dissecting the causal effect of CRP on the risk of T2DM. First, SNP-CRP, and SNP-T2DM association information were extracted from two types of samples, including major and minor alleles for each SNP, minor allele frequency (MAF), standard error (SE), and beta coefficients for each effect allele, both in the T2DM and CRP GWAS databases. Then, suitable SNPs were selected as IVs based on the following principles: (1) SNPs should be significantly associated with CRP but not with T2DM; (2) SNPs with LD should be removed. Finally, a pooled analysis for those IVs based on MR analysis was conducted to evaluate the influence of CRP on T2DM. To assess the effect of a single SNP on the result, a sensitivity analysis was carried out using the leave-one-out method. Pleiotropy in a MR study occurs when the effect of a SNP of IVs on T2DM (the outcome) could be independent of CRP (the intermediate phenotype). Here MR-Egger regression (Bowden et al., 2015) was used to explore the possibility of pleiotropy.
Figure 2. Overview of our method demonstrating the basic ideas for dissecting the influence of CRP on T2DM.
Samples for Associations Between Genetic Variants and CRP
Significantly associated SNPs of CRP were sourced from (Elliott et al., 2009). They carried out a genome-wide study and a replication study to identify loci associated with plasma CRP concentrations. Their data collection took place from 1989 to 2008, and genotyping occurred from 2003 to 2008. In total, 17,967 participants from five studies (Wyszynski et al., 2005; Bouatia-Naji et al., 2008; Chambers et al., 2008; Yuan et al., 2008; Sabatti et al., 2009) were collected for the genome-wide association study (GWAS), and another 14,747 participants who were not included in the GWAS were collected for the replication study of SNPs associated with CRP. The participants were from West London in the United Kingdom, the Northern Finnish Birth Cohort, the northernmost provinces of Finland, the Lausanne Cohort, and France. All of these participants were 30–75 years of age. Genetic loci associated with CRP were identified from GWAS. The single most closely associated SNP of CRP for each genetic locus was further replicated.
Samples for Associations Between Genetic Variants and T2DM
Associations between SNPs and T2DM were obtained from the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) consortium (Zeggini et al., 2008; Lango Allen et al., 2010; Morris et al., 2012), which is a grouping of researchers with shared interest in performing large-scale studies to characterize the genetic basis of T2DM in individuals of European descent. The initial instance of DIAGRAM (DIAGRAM v1) combined data from the Diabetes Genetics Initiative, the Finland-United States Investigation of NIDDM Genetics, and the Wellcome Trust Case Control Consortium (Zeggini et al., 2008; Lango Allen et al., 2010; Morris et al., 2012). An incremental meta-analysis (DIAGRAM v2) added another five GWAS studies of European-descent samples. The recent meta-analysis (DIAGRAM v3) collected 12,171 cases and 56,862 controls.
We extracted the significantly associated SNPs (P-value < 5 × 10-8) of CRP from (Elliott et al., 2009). Then, we extracted the association between these SNPs and T2DM to remove T2DM-associated SNPs (P-value < 5 × 10-2). The remaining SNPs were further assessed using the web tool SNAP1 (Johnson et al., 2008) to remove SNPs in LD (International HapMap Consortium et al., 2007). The remaining SNPs were selected as IVs for further analysis.
Figure 2 shows the process of pooled analysis based on the selected SNPs (IVs), which includes MR analysis, sensitivity analysis, and pleiotropy analysis.
We first calculated the influence of individual SNPs on T2DM through CRP as individual-level data and then calculated the pooled influence of these individual-level data using the penalized robust inverse-variance weighted (IVW) method (Bowden et al., 2015; Burgess et al., 2016). Assuming X, Y, and Z are CRP, T2DM, and selected SNPs (IVs), respectively, the Wald ratio (βXY) of CRP to T2DM for a specified SNP is defined as the following.
where βZY represents the per-allele log(OR) of T2DM based on Morris et al.’s (2012) study. βZX is the per-allele log(OR) of CRP based on Locke et al.’s (2015) study. In addition, the SE of the BMI-T2DM association of each Wald ratio is defined as the following.
where SEZY and SEZX represent the SE of the variant-T2DM and variant-CRP associations from corresponding studies, respectively. Based on Equations 1 and 2, the individual-level data of the influence of each SNP (IV) on T2DM through CRP was obtained. These data were then analyzed using the penalized robust IVW method (Burgess et al., 2016). The IVW estimate is the same as the two-stage least squares estimate using individual-level data. The penalization of weights from candidate instruments produces heterogeneity in causal estimates, which generates the penalized IVW method, robust IVW method, and penalized robust IVW methods; these methods have been demonstrated to improve the robustness of the findings.
Sensitivity Analysis Based on Leave-One-Out Validation
We conducted a leave-one-out validation to test the sensitivity of the selected SNPs (IVs). First, a SNP among the IVs was removed from the IVs to carry out a penalized robust IVW estimate. Then, the fluctuation of the results before and after removing the SNP was observed as the sensitivity. This process was repeated for each of the IVs.
Pleiotropy Analysis Using MR-Egger Regression
MR-Egger regression (Bowden et al., 2015), an adaption of Egger regression, was utilized here to detect the directional pleiotropy of IVs. This regression permits variants to show their pleiotropy and assumes that each variant is valid. Only when every variant used for IVs is independent of the alternative path affecting the disease (i.e., horizontal pleiotropy) can yield a consistent estimate via MR-Egger regression (Bowden et al., 2015). The pleiotropy bias of MR analysis is regarded as analogous to small-study bias. Additionally, the average pleiotropic effect across the genetic variants was captured by the estimated value of the intercept in Egger regression. Thus, an intercept that differs from zero is indicative of overall directional pleiotropy.
IVs for Pooled Analysis
Five SNPs of five genetic loci were significantly associated with CRP (P-value < 5 × 10-8) based on the previous genome-wide study and the replication study (Elliott et al., 2009). One of them (rs4420638) was associated with T2DM and so was removed. Since the remaining four SNPs were assessed without LD, these SNPs were eventually selected for the pooled analysis in Table 1. Each line of the table documents 13 items involving the SNP, EA and its frequencies, the beta coefficients of the SNP on the risk of BMI and T2DM, and SEs.
Pooled Analysis Result
Figure 3 shows the effect of each SNP and the combined effect of all four SNPs on T2DM through CRP using the penalized robust IVW method. The pooled result (OR 1.114048; 95% CI 1.058656 to 1.172338; P = 0.024) showed that high CRP significantly increases the risk of T2DM.
Figure 4 shows the sensitivity analysis result of SNPs based on leave-one-out validation. The ORs undergo substantial change after removing rs4537545 or rs7553007. In comparison, the ORs undergo a smaller change after removing rs6700896 or rs1183910. These results demonstrated that rs4537545 or rs7553007 drive the penalized robust IVW estimate.
The further estimate of the horizontal pleiotropic effect of these SNPs (IVs) was conducted via MR-Egger analysis. As a result, we got an intercept of 0.017 (95% CI 0.001 to 0.033, P = 0.039). This means a potential horizontal pleiotropic effect could bias our estimates.
In this study, we conducted a MR analysis to explore the causal effect of CRP on the risk of T2DM. To reduce the confounding effects, genetic variants were selected as IVs for MR analysis. Five SNPs associated significantly with CRP were extracted from Elliott et al.’s (2009) study. One of them (rs4420638) was associated with T2DM and was removed (Zeggini et al., 2008; Lango Allen et al., 2010; Morris et al., 2012). As a result, the remaining four SNPs without LD were deemed as IVs. The combined effect of the four SNPs on T2DM through CRP (OR 1.114048; 95% CI 1.058656 to 1.172338; P = 0.024) was determined using the penalized robust IVW method. This showed that high CRP significantly increases the risk of T2DM.
Brunner et al. (2008) utilized MR analysis for an analogous purpose. They reached the conclusion that associations between CRP and diabetes are likely to be noncausal. In that study, three SNPs were utilized as IVs for MR analysis. The associations between SNPs and CRP were measured on 5,274 men and women, and the associations between SNPs and diabetes were observed among 1,923 patients and 2,932 controls. In comparison, our study included more samples (17,967 and 69,033 participants for measuring SNPs-CRP and SNPs-diabetes associations, respectively). In addition, the diabetes cases in Brunner et al.’s (2008) study include type 1 Diabetes Mellitus (T1DM), T2DM, and others. In contrast, our study focused on T2DM patients. Overall, our study is more specific, and the corresponding conclusion should be more reliable.
To further validate that the different conclusion of Brunner et al. (2008) may be caused by their use of multiple types of diabetes samples, we conducted a MR analysis to evaluate the causal effect of CRP on the risk of T1DM. Associations between SNPs and T1DM were sourced from a meta-analysis of genome-wide association studies with 5,913 T1DM cases and 8,828 reference samples of European ancestry (Barrett et al., 2009), which identified three SNPs (rs6700896, rs7553007, rs1183910) of the five CRP-associated SNPs identified in our study. The associations between these three SNPs and T1DM were then extracted in Table 2 as IVs. The pooled result (OR 1.017145; 95% CI 0.9066489 to 1.14225; P = 0.909) did not support a causal role of CRP in the onset of T1DM using the penalized robust IVW method. Therefore, the T1DM samples in Brunner et al.’s (2008) study may be responsible for confusing the causality of CRP in T2DM.
The advantage of MR analysis is that no confounding factor should be considered when using genetic variants as IVs, since the genetic variants are free and unaffected by confounding factors. In comparison, confounding factors can heavily affect the results of observational studies. Therefore, most of these observational studies should be adjusted for potential confounding factors. However, abundant phenotypes of T2DM make it hard to do so. This issue occurs in Lee et al.’s (2009) prospective study about whether adiponectin is a confounder. In their study, observed that CRP was associated with T2DM after adjusting for age, sex, body mass index (BMI), etc. However, the association was completely attenuated after further adjustment for serum adiponectin. The recent evidence suggests that CRP inhibits adiponectin gene expression in adipocytes (Yuan et al., 2007). If so, adjustment for adiponectin would be over-adjustment.
Current studies on the treatment of T2DM have also reflected the potential causal effect of CRP on the risk of T2DM. During the treatment of T2DM patients with Metformin and Silymarin, researchers observed a significant decrease of CRP levels (Koujan et al., 2015; Shi et al., 2016). They even used these decreased serum levels of CRP as a sensitive predictor in T2DM patients being treated with drugs (Koujan et al., 2015; Shi et al., 2016). To explore effective curative options for T2DM, some researchers mainly refer to the effect of treatment on the serum CRP levels of T2DM patients. Kim (2014) observed that training modes can help to decrease the serum CRP level of T2DM patients. Rutten et al. conducted intensive care and routine care on 235 and 189 T2DM patients, respectively. A significant decrease of CRP levels based on intensified multifactorial treatment was observed in T2DM patients after six years (den Ouden et al., 2015).
The incidence of T2DM can be associated with multiple phenotypic factors. The potential linkages among these factors can help us to comprehensively understand the mechanisms of T2DM. The role of CRP on the risk of T2DM necessitates a measurable effect for us to mine these linkages. To the best of our knowledge, CRP is an annular, pentameric protein in blood plasma, the levels of which rise in response to inflammation. Thus, long-term inflammation can also increase the risk of T2DM through CRP. In the previous MR analysis, body mass index (BMI), and waist circumference (WC) are two immediate phenotypes of T2DM (Corbin et al., 2016). Stronger associations between these two phenotypes and CRP were observed in a recent study (Priyanka et al., 2013). This means that either CRP affects T2DM through BMI and WC, or BMI and WC affect T2DM through CRP.
Although our samples were restricted to those of European ancestry and IVs with LD were removed to reduce bias, our analysis also has limitations worth consideration. First, IV-CRP and IV-T2DM associations were sourced from two different samples. In theory, a single sample source may be more reliable than multiple sample sources. Fortunately, the increase of the number of samples reduces this type of bias. Next, horizontal pleiotropy of IVs may influence the conclusion of MR method. Here we evaluated the relevance between genes that IVs were located in (Table 1) with T2DM using an enrichment tool BLAT2DOLite (Cheng et al., 2016). No associations were detected, which means that these IVs cannot impact T2DM directly. Nevertheless, they could take effect in other immediate phenotypes of T2DM, such as BMI, WC, etc. The lack of associations between IVs and other phenotypes limits our further validation. Finally, without a specified indication, these results are ill-suited for clinical application.
In summary, we validated that the increase of CRP levels enhances the risk of T2DM and has no effect on T1DM using MR analysis. Here, genetic variations are deemed as IVs to reduce the potential confounding effect. The further validation relies on RCTs. The result has a very positive guiding role in finding new therapeutic strategies and therapeutic targets. This is related to the early prevention of T2DM, which has the potential to indirectly, from the intervention phenotype, regulate the development of disease. In the future, we will construct a dataset including the relationships between genetic variations and the phenotypes of DM. The data set may be beneficial in designing computational tools for phenotype interaction (Peng et al., 2017, 2018a) and function association prediction (Cheng et al., 2018b; Hu et al., 2018; Peng et al., 2018b).
LC, HJ, JZ, and SW conceived and designed the experiments. LC, HZ, and SY analyzed the data. LC wrote the manuscript. All authors read and approved the final manuscript.
This work was supported by the Major State Research Development Program of China (No. 2016YFC1202302), the National Natural Science Foundation of China (Grant Nos. 61871160 and 61502125), the Heilongjiang Postdoctoral Fund (Grant Nos. LBH-TZ20 and LBH-Z15179), the China Postdoctoral Science Foundation (Grant Nos. 2018T110315 and 2016M590291), and Scientific Research Project of Health and Family Planning Commission of Heilong Jiang Province (Grant No. 2018415 to HJ).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Guiyou Liu for the improvement of this manuscript.
Barrett, J. C., Clayton, D. G., Concannon, P., Akolkar, B., Cooper, J. D., Erlich, H. A., et al. (2009). Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat. Genet. 41, 703–707. doi: 10.1038/ng.381
Bouatia-Naji, N., Rocheleau, G., Van, Lommel L, Lemaire, K., Schuit, F., Cavalcanti-Proença, C., et al. (2008). A polymorphism within the G6PC2 gene is associated with fasting plasma glucose levels. Science 320, 1085–1088. doi: 10.1126/science.1156849
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080
Brunner, E. J., Kivimäki, M., Witte, D. R., Lawlor, D. A., Davey, Smith G, Cooper, J. A., et al. (2008). Inflammation, insulin resistance, and diabetes–Mendelian randomization using CRP haplotypes points upstream. PLoS Med. 5:e155. doi: 10.1371/journal.pmed.0050155
Burgess, S., Bowden, J., Dudbridge, F., and Thompson, S. G. (2016). Robust instrumental variable methods using multiple candidate instruments with application to Mendelian randomization. arXiv 2016:1606.03729.
Chambers, J. C., Elliott, P., Zabaneh, D., Zhang, W., Li, Y., Froguel, P., et al. (2008). Common genetic variation near MC4R is associated with waist circumference and insulin resistance. Nat. Genet. 40, 716–718. doi: 10.1038/ng.156
Cheng, L., Hu, Y., Sun, J., Zhou, M., and Jiang, Q. (2018a). DincRNA: a comprehensive web-based bioinformatics toolkit for exploring disease associations and ncRNA function. Bioinformatics 34, 1953–1956. doi: 10.1093/bioinformatics/bty002
Cheng, L., Jiang, Y., Ju, H., Sun, J., Peng, J., Zhou, M., et al. (2018b). InfAcrOnt: calculating cross-ontology term similarities using information flow by a random walk. BMC Genomics 19:919. doi: 10.1186/s12864-017-4338-6
Cheng, L., Wang, P., Tian, R., Wang, S., Guo, Q., Luo, M., et al. (2018c). LncRNA2Target v2.0: a comprehensive database for target genes of lncRNAs in human and mouse. Nucleic Acids Res. doi: 10.1093/nar/gky1051 [Epub ahead of print].,
Cheng, L., Yang, H., Zhao, H., Pei, X., Shi, H., Sun, J., et al. (2017). MetSigDis: a manually curated resource for the metabolic signatures of diseases. Brief. Bioinform. doi: 10.1093/bib/bbx103 [Epub ahead of print].
Cheng, L., Zhang, S., and Hu, Y. (2016). BLAT2DOLite: an online system for identifying significant relationships between genetic sequences and diseases. PLoS One 11:e0157274. doi: 10.1371/journal.pone.0157274
Corbin, L. J., Richmond, R. C., Wade, K. H., Burgess, S., Bowden, J., Smith, G. D., et al. (2016). Body mass index as a modifiable risk factor for type 2 diabetes: refining and understanding causal estimates using Mendelian randomisation. Diabetes 65:db160418. doi: 10.2337/db16-0418
den, Ouden H, Berends, J., Stellato, R. K., Beulens, J. W., and Rutten, G. E. (2015). Effect of six years intensified multifactorial treatment on levels of hs-CRP and adiponectin in patients with screen detected type 2 diabetes: the ADDITION-Netherlands randomized trial. Diabetes Metab. Res. Rev. 31, 758–766. doi: 10.1002/dmrr.2669
Doi, Y., Kiyohara, Y., Kubo, M., Ninomiya, T., Wakugawa, Y., Yonemoto, K., et al. (2005). Elevated C-reactive protein is a predictor of the development of diabetes in a general Japanese population: the Hisayama Study. Diabetes Care 28:2497. doi: 10.2337/diacare.28.10.2497
Elliott, P. P., Chambers, D. J. C., Zhang, D. W., Clarke, R., Hopewell, J. C., Peden, J. F., et al. (2009). Genetic loci influencing c-reactive protein levels and risk of coronary heart disease. JAMA J. Am. Med. Assoc. 302:37. doi: 10.1001/jama.2009.954
Festa, A., D’Agostino, R Jr, Tracy, R. P., Haffner, S. M., and Insulin, Resistance Atherosclerosis Study (2002). Elevated levels of acute-phase proteins and plasminogen activator inhibitor-1 predict the development of type 2 diabetes: the insulin resistance atherosclerosis study. Diabetes 51, 1131–1137. doi: 10.2337/diabetes.51.4.1131
Freeman, D. J., Norrie, J., Caslake, M. J., Gaw, A., Ford, I., Lowe, G. D., et al. (2002). C-reactive protein is an independent predictor of risk for the development of diabetes in the West of Scotland Coronary Prevention Study. Diabetes 51, 1596–1600. doi: 10.2337/diabetes.51.5.1596
Gentile, S., Piscitelli, P., Viazzi, F., Russo, G., Ceriello, A., Giorda, C., et al. (2017). Antihyperglycemic treatment in patients with type 2 diabetes in Italy: the impact of age and kidney function. Oncotarget 8, 62039–62048. doi: 10.18632/oncotarget.18816
Hartwig, F. P., Davies, N. M., Hemani, G., and Davey Smith, G. (2017). Two-sample Mendelian randomization: avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int. J. Epidemiol. 45, 1717–1726. doi: 10.1093/ije/dyx028
International HapMap Consortium, Frazer, K. A., Ballinger, D. G., Cox, D. R., Hinds, D. A., Stuve, L. L., et al. (2007). A second generation human haplotype map of over 3.1 million SNPs. Nature 449, 851–861. doi: 10.1038/nature06258
Johnson, A. D., Handsaker, R. E., Pulit, S. L., Nizzari, M. M., O’Donnell, C. J., and de, Bakker PI (2008). SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 24:2938. doi: 10.1093/bioinformatics/btn564
Kim, K. B. (2014). Effect of different training mode on Interleukin-6 (IL-6) and C-reactive protein (CRP) in type 2 diabetes mellitus (T2DM) patients. J. Exerc. Nutr. Biochem. 18, 371–378. doi: 10.5717/jenb.2014.18.4.371
Kohn, L. D., Wallace, B., Schwartz, F., and McCall, K. (2005). Is type 2 diabetes an autoimmune-inflammatory disorder of the innate immune system? Endocrinology 146, 4189–4191. doi: 10.1210/en.2005-0920
Koujan, S. E., Gargari, B. P., Mobasseri, M., Valizadeh, H., and Asghari-Jafarabadi, M. (2015). Effects of Silybum marianum (L.) Gaertn. (silymarin) extract supplementation on antioxidant status and hs-CRP in patients with type 2 diabetes mellitus: a randomized, triple-blind, placebo-controlled clinical trial. Phytomedicine 22, 290–296. doi: 10.1016/j.phymed.2014.12.010
Krakoff, J., Funahashi, T., Stehouwer, C. D., Schalkwijk, C. G., Tanaka, S., Matsuzawa, Y., et al. (2003). Inflammatory markers, adiponectin, and risk of type 2 diabetes in the Pima Indian. Diabetes Care 26, 1745–1751. doi: 10.2337/diacare.26.6.1745
Lango Allen H, Estrada, K., Lettre, G., Berndt, S. I., Weedon, M. N., Rivadeneira, F., et al. (2010). Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 467, 832–838. doi: 10.1038/nature09410
Lawlor, D. A., Harbord, R. M., Sterne, J. A., Timpson, N., and Davey Smith, G. (2008). Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat. Med. 27, 1133–1163. doi: 10.1002/sim.3034
Lee, C. C., Adler, A. I., Sandhu, M. S., Sharp, S. J., Forouhi, N. G., Erqou, S., et al. (2009). Association of C-reactive protein with type 2 diabetes: prospective analysis and meta-analysis. Diabetologia 52, 1040–1047. doi: 10.1007/s00125-009-1338-3
Locke, A. E., Kahali, B., Berndt, S. I., Justice, A. E., Pers, T. H., Day, F. R., et al. (2015). Genetic studies of body mass index yield new insights for obesity biology. Nature 518, 197–206. doi: 10.1038/nature14177
Morris, A. P., Voight, B. F., Teslovich, T. M., Ferreira, T., Segrè, A. V., Steinthorsdottir, V., et al. (2012). Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat. Genet. 44, 981–990. doi: 10.1038/ng.2383
Nakanishi, S., Yamane, K., Kamei, N., Okubo, M., and Kohno, N. (2003). Elevated C-reactive protein is a risk factor for the development of type 2 diabetes in Japanese Americans. Diabetes Care 26, 2754–2757. doi: 10.2337/diacare.26.10.2754
Noyce, A. J., Kia, D. A., and Hemani, G. (2017). Estimating the causal influence of body mass index on risk of Parkinson disease: a mendelian randomisation study. PLoS Med. 14:e1002314. doi: 10.1371/journal.pmed.1002314
Peng, J., Zhang, X., Hui, W., Lu, J., Li, Q., Liu, S., et al. (2018b). Improving the measurement of semantic similarity by combining gene ontology and co-functional network: a random walk based approach. BMC Syst. Biol. 12:18. doi: 10.1186/s12918-018-0539-0
Peng, J., Xue, H., Shao, Y., Wang, Y., Chen, J., and Shang, X. (2017). A novel method to measure the semantic similarity of HPO terms. Int. J. Data Min. Bioinform. 17, 173–188. doi: 10.1504/IJDMB.2017.084268
Priyanka, N., Kumari, M., Kalra, N., Arjun, P., Naik, S. B., and Pradeep, A. R. (2013). Crevicular fluid and serum concentrations of progranulin and high sensitivity CRP in chronic periodontitis and type 2 diabetes. Dis. Mark. 35, 389–394. doi: 10.1155/2013/803240
Sabatti, C., Service, S. K., Hartikainen, A. L., Pouta, A., Ripatti, S., Brodsky, J., et al. (2009). Genome-wide association analysis of metabolic traits in a birth cohort from a founder population. Nat. Genet. 41, 35–46. doi: 10.1038/ng.271
Shi, L., Tan, G., and Zhang, K. (2016). Relationship of the serum CRP level with the efficacy of metformin in the treatment of type 2 diabetes mellitus: a meta-analysis. J. Clin. Lab. Anal. 30, 13–22. doi: 10.1002/jcla.21803
Svensson, M. K., Afghahi, H., Franzen, S., Björk, S., Gudbjörnsdottir, S., Svensson, A. M., et al. (2017). Decreased systolic blood pressure is associated with increased risk of all-cause mortality in patients with type 2 diabetes and renal impairment: a nationwide longitudinal observational study of 27,732 patients based on the Swedish National Diabetes Register. Diab. Vasc. Dis. Res. 14, 226–235. doi: 10.1177/1479164116683637
Thorand, B., Baumert, J., Kolb, H., Meisinger, C., Chambless, L., Koenig, W., et al. (2007). Sex differences in the prediction of type 2 diabetes by inflammatory markers: results from the MONICA/KORA Augsburg case-cohort study, 1984-2002. Diabetes Care 30, 854–860. doi: 10.2337/dc06-1693
Wyszynski, D. F., Waterworth, D. M., Barter, P. J., Cohen, J., Kesäniemi, Y. A., Mahley, R. W., et al. (2005). Relation between atherogenic dyslipidemia and the adult treatment program-III definition of metabolic syndrome (genetic epidemiology of metabolic syndrome project). Am. J. Cardiol. 95, 194–198. doi: 10.1016/j.amjcard.2004.08.091
Xu, Y. L., Xu, K. F., Bai, J. L., Liu, Y., Yu, R. B., Liu, C. L., et al. (2016). Elevation of serum uric acid and incidence of type 2 diabetes: a systematic review and meta-analysis. Chronic Dis. Transl. Med. 2, 81–91. doi: 10.1016/j.cdtm.2016.09.003
Yavorska, O. O., and Burgess, S. (2017). MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int. J. Epidemiol. 46, 1734–1739. doi: 10.1093/ije/dyx034
Yuan, G., Chen, X., Ma, Q., Qiao, J., Li, R., Li, X., et al. (2007). C-reactive protein inhibits adiponectin gene expression and secretion in 3T3-L1 adipocytes. J. Endocrinol. 194, 275–281. doi: 10.1677/JOE-07-0133
Yuan, X., Waterworth, D., Perry, J. R., Lim, N., Song, K., Chambers, J. C., et al. (2008). Population-based genome-wide association studies reveal six loci influencing plasma levels of liver enzymes. Am. J. Hum. Genet. 83, 520–528. doi: 10.1016/j.ajhg.2008.09.012
Zeggini, E., Scott, L. J., Saxena, R., Voight, B. F., Marchini, J. L., Hu, T., et al. (2008). Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes. Nat. Genet. 40, 638–645. doi: 10.1038/ng.120
Zeng, X., Liu, L., Lü, L., and Zou, Q. (2018). Prediction of potential disease-associated microRNAs using structural perturbation method. Bioinformatics 34, 2425–2432. doi: 10.1093/bioinformatics/bty112
Keywords: C-reactive protein, type 2 diabetes mellitus, causal effect, genome-wide association studies, mendelian randomization, type 1 diabetes mellitus
Citation: Cheng L, Zhuang H, Yang S, Jiang H, Wang S and Zhang J (2018) Exposing the Causal Effect of C-Reactive Protein on the Risk of Type 2 Diabetes Mellitus: A Mendelian Randomization Study. Front. Genet. 9:657. doi: 10.3389/fgene.2018.00657
Received: 14 September 2018; Accepted: 03 December 2018;
Published: 20 December 2018.
Edited by:Yan Huang, Harvard Medical School, United States
Reviewed by:Xinjun Zhang, Thermo Fisher Scientific, United States
Jorge L. Del-Aguila, Washington University in St. Louis, United States
Juan Wang, Inner Mongolia University, China
Copyright © 2018 Cheng, Zhuang, Yang, Jiang, Wang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.