Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study

Background C-reactive protein (CRP) is associated with immune, cardiometabolic, and psychiatric traits and diseases. Yet it is inconclusive whether these associations are causal. Methods and Findings We performed Mendelian randomization (MR) analyses using two genetic risk scores (GRSs) as instrumental variables (IVs). The first GRS consisted of four single nucleotide polymorphisms (SNPs) in the CRP gene (GRSCRP), and the second consisted of 18 SNPs that were significantly associated with CRP levels in the largest genome-wide association study (GWAS) to date (GRSGWAS). To optimize power, we used summary statistics from GWAS consortia and tested the association of these two GRSs with 32 complex somatic and psychiatric outcomes, with up to 123,865 participants per outcome from populations of European ancestry. We performed heterogeneity tests to disentangle the pleiotropic effect of IVs. A Bonferroni-corrected significance level of less than 0.0016 was considered statistically significant. An observed p-value equal to or less than 0.05 was considered nominally significant evidence for a potential causal association, yet to be confirmed. The strengths (F-statistics) of the IVs were 31.92–3,761.29 and 82.32–9,403.21 for GRSCRP and GRSGWAS, respectively. CRP GRSGWAS showed a statistically significant protective relationship of a 10% genetically elevated CRP level with the risk of schizophrenia (odds ratio [OR] 0.86 [95% CI 0.79–0.94]; p < 0.001). We validated this finding with individual-level genotype data from the schizophrenia GWAS (OR 0.96 [95% CI 0.94–0.98]; p < 1.72 × 10−6). Further, we found that a standardized CRP polygenic risk score (CRPPRS) at p-value thresholds of 1 × 10−4, 0.001, 0.01, 0.05, and 0.1 using individual-level data also showed a protective effect (OR < 1.00) against schizophrenia; the first CRPPRS (built of SNPs with p < 1 × 10−4) showed a statistically significant (p < 2.45 × 10−4) protective effect with an OR of 0.97 (95% CI 0.95–0.99). The CRP GRSGWAS showed that a 10% increase in genetically determined CRP level was significantly associated with coronary artery disease (OR 0.88 [95% CI 0.84–0.94]; p < 2.4 × 10−5) and was nominally associated with the risk of inflammatory bowel disease (OR 0.85 [95% CI 0.74–0.98]; p < 0.03), Crohn disease (OR 0.81 [95% CI 0.70–0.94]; p < 0.005), psoriatic arthritis (OR 1.36 [95% CI 1.00–1.84]; p < 0.049), knee osteoarthritis (OR 1.17 [95% CI 1.01–1.36]; p < 0.04), and bipolar disorder (OR 1.21 [95% CI 1.05–1.40]; p < 0.007) and with an increase of 0.72 (95% CI 0.11–1.34; p < 0.02) mm Hg in systolic blood pressure, 0.45 (95% CI 0.06–0.84; p < 0.02) mm Hg in diastolic blood pressure, 0.01 ml/min/1.73 m2 (95% CI 0.003–0.02; p < 0.005) in estimated glomerular filtration rate from serum creatinine, 0.01 g/dl (95% CI 0.0004–0.02; p < 0.04) in serum albumin level, and 0.03 g/dl (95% CI 0.008–0.05; p < 0.009) in serum protein level. However, after adjustment for heterogeneity, neither GRS showed a significant effect of CRP level (at p < 0.0016) on any of these outcomes, including coronary artery disease, nor on the other 20 complex outcomes studied. Our study has two potential limitations: the limited variance explained by our genetic instruments modeling CRP levels in blood and the unobserved bias introduced by the use of summary statistics in our MR analyses. Conclusions Genetically elevated CRP levels showed a significant potentially protective causal relationship with risk of schizophrenia. We observed nominal evidence at an observed p < 0.05 using either GRSCRP or GRSGWAS—with persistence after correction for heterogeneity—for a causal relationship of elevated CRP levels with psoriatic osteoarthritis, rheumatoid arthritis, knee osteoarthritis, systolic blood pressure, diastolic blood pressure, serum albumin, and bipolar disorder. These associations remain yet to be confirmed. We cannot verify any causal effect of CRP level on any of the other common somatic and neuropsychiatric outcomes investigated in the present study. This implies that interventions that lower CRP level are unlikely to result in decreased risk for the majority of common complex outcomes.

statistically significant. An observed p-value equal to or less than 0.05 was considered nominally significant evidence for a potential causal association, yet to be confirmed.
Recent large-scale MR studies, focusing mainly on cardiovascular disease and metabolic traits, failed to show a causal association between CRP level and these outcomes (S1 Table). This has led to the notion that elevated CRP levels do not causally contribute to these traits and disorders. However, these studies used a single CRP-associated single nucleoid polymorphism (SNP) or a very limited set of CRP-associated SNPs (S1 Table). Common SNPs serving as proxies for CRP level represent only a small effect on CRP level per se and thus require a large enough sample size to detect causal effects on the outcome. Moreover, most studies have generally included a limited range of common complex diseases, often not more than two or three outcomes, or they have been performed in a single or small population, yielding inadequate study power (S1 Table). In other words, the evidence for a causal relationship between CRP and a broad range of common traits or diseases remains inconclusive. This is mostly due to the lack of well-powered MR studies that use optimally informative genetic IVs for CRP. Here, we sought to comprehensively examine the hypothesis that genetically determined CRP level directly contributes to common somatic and psychiatric outcomes. To optimize IV power, we applied a MR approach using summary statistics from large-scale genome-wide association study (GWAS) consortia of 32 somatic and psychiatric phenotypes for the four CRP variants representing 98% of the common variation in the CRP gene and for the largest known set of independent SNPs known to be associated with CRP. We further aimed to confirm the identified association between CRP and schizophrenia using a CRP polygenic risk score (CRP PRS ) from individual-level genotype data from the largest consortium of schizophrenia to date. We performed an in silico pathway analysis (see Discussion) to provide insights into the possible mechanism underlying the observed association of CRP level with schizophrenia.

Study Design and Rationale
The present MR study consists of two key components. First, we used established gene variants associated with CRP level and combined them to build two genetic risk scores (GRSs) for CRP. The first GRS consisted of only four SNPs in the CRP gene (GRS CRP ) selected from the largest recent MR study of CRP [29], and the second consisted of 18 SNPs that were associated with CRP level at a genome-wide significance level in the largest GWAS for CRP to date (GRS GWAS ) [30]. Second, we obtained summary association statistics from GWAS consortia for a panel of 32 common somatic and psychiatric outcomes (Table 1). BPP and BZA selected the studies, and contacted each consortium with a standardized request for study data, including the name of the study or consortium, the number of cases and controls, the number of available CRP SNPs for GRS CRP and GRS GWAS , and the estimated effect for each SNP (or its proxy) on outcome, i.e., the per allele regression coefficient with standard error or the odds ratio (OR) and corresponding 95% confidence interval. Data were available for 32 different outcomes in five broad disease/trait classes (autoimmune/inflammatory, cardiovascular, metabolic, neurodegenerative, and psychiatric), including at least 1,566, and up to 184,305, participants per outcome from populations of European ancestry (Table 1). These outcomes were selected based on the following two inclusion criteria: (i) the outcome having been associated with CRP level in epidemiological studies and (ii) availability of large meta-GWAS analyses for the outcome (Table 1).

Genetic Instruments
Weak IVs yielding insufficient statistical power may have hampered estimation of causal effects of CRP on the outcomes in previous analyses (S1 Table). Our MR approach, by using GWAS data and combining multiple independent SNPs into a GRS (i.e., IV), has the potential to greatly increase power. The selected SNPs have been described elsewhere [30,55,56] and are further detailed in S2-S4 Tables. These IVs were used to test the combined effect of the associations of CRP-level-influencing alleles with the outcomes. Our approach was implemented in such a way that the effects of both independent SNPs in the CRP gene (GRS CRP ) [55,56] (S1 Methods) and independent SNPs known to be genome-wide significantly associated with CRP levels (GRS GWAS ) [30], as well as pleiotropic effects of SNPs, could be discriminated [57]. Pleiotropy exists if CRP SNPs influence exposures (risk factors) other than CRP level and therefore violate one of the key MR assumptions.

Statistical Analysis
All analyses were done using the GRS function implemented in the grs.summary module of the R package Genetics ToolboX (version 2.15.1 for Windows). The grs.summary module approximates the regression of an outcome onto an additive GRS, using only single SNP association summary statistics extracted from GWAS results. The method is described in more detail elsewhere [58]. In brief, we performed MR analyses using GRS IVs in two steps. First, we used four individual CRP gene SNPs (i.e., IVs) associated with CRP level [56,59] (S2 and S3 Tables) to create a weighted GRS, named GRS CRP, corresponding to the joint effect of the four SNPs within the CRP gene [55]. We extracted ω (the estimated coefficient, or weight) for individual SNPs from association results reported by the CRP Coronary Heart Disease Genetics Collaboration (CCGC) [29,55]; ω represents a one-unit (in mg/l) increase of the natural log of CRP level (lnCRP) per dose of the coded allele. The four tagging SNPs represent 98% of the common variation in the CRP gene, assuming a minor allele frequency of 0.05 and an r 2 threshold of 0.8, and aggregately explain~2% of the total variation (i.e., phenotypic variance) in serum CRP level in populations of European descent [55,59]. Second, we constructed a multilocus GRS, named GRS GWAS , that combined 18 SNPs associated with serum CRP level at a genomewide significance level (p < 5×10 −8 ; S2 and S3 Tables), derived from a large meta-GWAS analysis of CRP conducted by the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium [30]. This multilocus GRS explains approximately~5% of the total variation in serum CRP level [30]. We integrated ω for each CRP SNP from the reference data of CCGC [55] or meta-analysis of GWASs [30] for CRP level with the summary association statistics extracted from the GWAS consortium data for each outcome (S1 Data; S2 Methods). This MR approach using meta-GWAS summary statistics data is equivalent to an inverse-variance-weighted meta-analysis and has previously been validated in comparison to individual-level data [57,60]. To estimate the causal effect of CRP level on an outcome, we obtained the β values (estimated effects from regression analysis) for the effects of CRP SNPs on the outcome, with standard errors, se β , from the corresponding GWAS results. Where no summary statistics for a CRP SNP in the GRS IVs were available in the look-up dataset, we chose the proxy SNP that had the highest linkage disequilibrium with the initial SNP (r 2 > 0.9 in HapMap release 22; S3 Table). If several proxy SNPs had the exact same r 2 value, we chose the proxy nearest to the original SNP in the instrument. Separate regressions of outcomes on GRSs were performed to calculate α IV estimators (i.e., causal IV estimators) for each outcome. Correspondingly, the value of a GRS is the sum of the ω values multiplied by the allele dosage (i.e., 0, 1, or 2) for each CRP SNP in the CCGC or in the CHARGE Consortium data [30,55]. For uncorrelated SNPs, when maximizing the likelihood function, the α IV value and its standard error, se α , can be approximated with the formula α ffi (Sω × β × se β −2 ) /(Sω 2 × se β −2 ), with se α ffi p 1/ Sω 2 × se β −2 . lnCRP was used as the outcome in reference studies [30,55], so in obtaining the ω values (i.e., effect sizes) for each of the CRP SNPs, a unit increase in lnCRP equals a 10 symmetric percentage (s%) increase in CRP level, which corresponds to a unit change in the level of a continuous outcome or logit of risk estimate (i.e., beta coefficient) for a dichotomous outcome [61]. The α IV value (i.e., causal estimate) for each CRP SNP is, therefore, presented for each outcome as corresponding to a 10-s% increase in actual CRP level. During the course of this study, an updated, larger GWAS dataset for coronary artery disease (CAD) became publicly available (CARDIoGRAMplusC4D Consortium, release 2015 [41]); we therefore redid the analysis for CAD using the release 2015 data.
To assess which SNPs might have violated the key MR assumption regarding pleiotropy, we performed goodness-of-fit tests to correct both GRSs for the heterogeneity of their corresponding SNPs' effects on each outcome. Heterogeneity, which indicates the potential presence of pleiotropy, was measured using the Q statistic and was considered statistically significant at a conservative uncorrected p-value of <0.05. Although heterogeneity could be an indicator of pleiotropy, there are other factors that could introduce heterogeneity in the analyses. Even though the adjustments for heterogeneity that we have made could be overconservative, we have used this method in order to minimize false positives. After stepwise removal of SNPs with potential pleiotropic effects, we repeated the analyses until significant heterogeneity was no longer observed.
To further ensure the strength of these two GRSs as IVs, we generated an F-statistic for each outcome. We used variance in lnCRP explained by each set of CRP SNPs (2% and 5%, respectively, for GRS CRP and GRS GWAS ) to calculate the F-statistic using the formula where R 2 represents the proportion of variability in CRP level that is explained by the GRS, n represents sample size, and K represents the number of IVs included in model (i.e., for this study K = 1) [62]. As a rule of thumb, an F-value above ten indicates that a causal estimate is unlikely to be biased due to weak instruments [57].

Multiple Testing
The present study included 32 independent sample sets. For each sample set, we did one statistical test, for which a global nominal significance level of 0.05 was considered as satisfactory to derive conclusions. The need for correction for multiple testing is debatable. Nevertheless, to ensure the validity of our conclusions, we took a conservative approach and applied a Bonferroni-corrected significance threshold calculated as 0.05 divided by 32 (i.e., 0.0016). We considered a statistical test with an observed p-value more than 0.05 as a definitely nonsignificant result, i.e., no association; an observed p-value equal to or less than 0.05 as nominally significant evidence for a potential, but yet to be confirmed, causal association; and an observed pvalue equal to or less than 0.0016 as statistically significant evidence for a causal association.

CRP Polygenic Risk Score and Schizophrenia Using Individual-Level Data
In an ancillary follow-up study, inspired by comments by the editors and the reviewers, we aimed further to determine whether GRS GWAS was causally associated with schizophrenia using individual-level data retrieved from the Psychiatric Genomics Consortium (PGC) schizophrenia dataset (S3 Methods) [54]. This dataset consisted of 36 independent cohorts with a combined 25,629 cases and 30,976 controls for which we had ethics approval (S4 Methods). Three family-based samples of European ancestry (1,235 parent-affected offspring trios) were excluded from our analysis. To evaluate whether the observed protective causal association between GRS GWAS and schizophrenia was persistent, we investigated whether the CRP PRS was also protectively associated with schizophrenia. Briefly, CRP PRS values were calculated for each individual by summing the total effect of the SNP dosages by their effect size. In addition to the 18 genome-wide significant CRP SNPs, we grouped subthreshold CRP-associated SNPs at the following p-value thresholds: 1 × 10 −4 , 0.001, 0.01, 0.05, and 0.1. Standardized CRP PRS values were tested for association with schizophrenia case status in each cohort with adjustment for ten principal components (PCs). A fixed effects inverse-variance-weighted meta-analysis was performed across all 36 cohorts to obtain the overall effect size estimate as explained in S4 Methods and elsewhere [63]. The variance in schizophrenia case status explained by CRP PRS was estimated using the deviation in Nagelkerke's pseudo-R 2 between a null model (which included ten PCs) and the full model (which included GRS in addition to the ten PCs), calculated in R using the Functions for Medical Statistics Book with Some Demographic Data (fmsb) R package (S3 Methods). Similar to previous studies, the statistical significance of CRP PRS values was estimated based on their logistic regression coefficient [64], and reported CRP PRS ORs correspond to a 1-SD increase in CRP PRS [65].

Results
Using GRS CRP , we first tested whether a CRP-gene-determined increase in lnCRP was associated with each outcome. In Table 2, the causal effects of lnCRP estimated for each outcome are summarized. We found no heterogeneity in the IV analyses (p heterogeneity 0.11 for all outcomes), and GRS CRP was a strong instrument ( GRS GWAS showed a statistically significant protective effect of lnCRP on the risk of schizophrenia (per 10-s% increase in CRP level, OR 0.86 [95% CI 0.79-0.94]; p < 0.0010) (Figs 1 and S1; Table 3). In a follow-up analysis using the individual-level PGC data, we found that a GRS incorporating the same 18 CRP SNPs used to construct the GRS GWAS was again significantly associated with a lower risk of schizophrenia (OR 0.96 [95% CI 0.94-0.98]; p < 1.72 × 10 −6 ). This signal persisted when we included all SNPs meeting a less stringent p-value threshold of 1 × 10 −4 (OR 0.97 [95% CI 0.95-0.99]; p < 2.45 × 10 −4 ). At less stringent p-value thresholds, less variance was explained by the logistic model, and the protective effect of CRP risk scores became less significant, but across all p-value thresholds, the direction of the effect was consistently protective (Figs 2 and 3). To ensure that the association between risk alleles for CRP and schizophrenia was not driven by a small number of genome-wide significant SNPs, we performed a leave-one-out sensitivity analysis of the 18 genome-wide SNPs. In the 18 sets of 17 SNPs, the variance explained (Nagelkerke's pseudo-R 2 ) ranged from 0.012% to 0.034%, with pvalues ranging from 9.3 × 10 −5 to 1.6 × 10 −2 , suggesting that the protective effect observed between risk alleles for CRP and schizophrenia was not driven by a small number of SNPs with large effects.
GRS GWAS also showed moderate but nominally significant effects of lnCRP on the risk of IBD (OR 0. 85   Effect size (95% CI) per 1-mg/l increase in lnCRP. For risk of disease, effect size is given as an OR, otherwise given in the specific units in which the outcome was measured. Derived from the IV causal estimator α. 2 Effect size unit is millimeters of mercury per 1-mg/l increase in lnCRP. 3 Effect size unit is standard deviations per 1-mg/l increase in lnCRP (the body mass index results were inverse normal transformed to a distribution with μ = 0 and σ = 1). 4 Effect size unit is milliliters/minute/1.73 m 2 per 1-mg/l increase in lnCRP. 5 Effect size unit is grams/deciliter per 1-mg/l increase in lnCRP. eGFR cr , estimated glomerular filtration rate from serum creatinine; F-value, F-statistic value for the genetic instrument; M, number of markers used in the genetic instrument; N, number of samples in the disease/trait meta-analysis; p-het, p-value of heterogeneity of effect test.  Effect size (95% CI) per 1-mg/l increase in lnCRP. For risk of disease, effect size is given as an OR, otherwise given in the specific units in which the outcome was measured. Derived from the IV causal estimator α. 2 Effect size unit is millimeters of mercury per 1-mg/l increase in lnCRP. 3 Effect size unit is standard deviations per 1-mg/l increase in lnCRP (the body mass index results were inverse normal transformed to a distribution with μ = 0 and σ = 1). 4 Effect size unit is milliliters/minute/1.73 m 2 per 1-mg/l increase in lnCRP. 5 Effect size unit is grams/deciliter per 1-mg/l increase in lnCRP.  (Table 3; S1 Fig).
Using GRS GWAS , there was no significant evidence of heterogeneity of the effect size for knee osteoarthritis, bipolar disorder, schizophrenia, or SBP, while the heterogeneity test was statistically significant for psoriatic arthritis, IBD, Crohn disease, CAD, DBP, eGFR cr , serum albumin, and serum protein. These heterogeneities in the effects of GRS GWAS may be attributable to pleiotropic effects of the SNPs used to build the GRS GWAS . We subsequently performed a stepwise removal of SNPs from GRS GWAS until no significant heterogeneity remained (Table 4). This adjustment in the GRS GWAS resulted in the removal of three SNPs from the GRS GWAS for IBD (in GCKR, IRF1, and PTPN2), five SNPs from the GRS GWAS for Crohn disease (in GCKR, IL6R, IRF1, PABPC4, and PTPN2), one SNP from the GRS GWAS for psoriatic arthritis (in IRF1), three SNPs for CAD (in APOC1, HNF1A, and IL6R), one SNP from the GRS GWAS for DBP (in PABPC4), two SNPs from the GRS GWAS for eGFR cr (in LEPR and GCKR), six SNPs from the GRS GWAS for serum albumin level (in APOC1, BCL7B, GCKR, PPP1R3B, PTPN2, and IRF1), and one SNP from the GRS GWAS for serum protein level (in GCKR). After removal of these variants from the GRS GWAS , we found no statistically significant (at p < 0.0016) association between genetically increased lnCRP level and any of these outcomes (Table 4). However, the effect estimate of CRP on DBP, serum albumin, and psoriatic arthritis showed nominal association at p < 0.05. For example, for DBP, 17 SNPs remained in the GRS GWAS and yielded a slightly lower causal estimate (compared to the values before adjustment) of a 0.39 (95% CI −0.01 to 0.78) mm Hg increase in DBP per 10-s% increase in lnCRP level, with a nominal significance of p < 0.05.
Likewise, we hypothesized that the fact that GRS GWAS showed a nonsignificant effect of CRP on celiac disease, ulcerative colitis, rheumatoid arthritis, type 1 diabetes, and type 2 diabetes can be to some extent explained by the significant heterogeneity observed for these outcomes (Table 3). The stepwise adjustment in the GRS GWAS resulted in the removal of two SNPs from the GRS GWAS for celiac disease (in PABPC4 and PTPN2), one SNP from the GRS GWAS for ulcerative colitis (in GCKR), five SNPS from the GRS GWAS for rheumatoid arthritis (in HNF4A, IL6R, SALL1, NLRP3, and PTPN2), one SNP from the GRS GWAS for type 1 diabetes (in PTPN2), and one SNP from the GRS GWAS for type 2 diabetes (in APOC1). After adjusting for heterogeneity, the association of GRS GWAS with these outcomes remained statistically nonsignificant (Table 4).

Discussion
In this large-scale cross-consortium MR study of 32 complex outcomes, we found evidence for a potential protective causal relationship between elevated CRP level and schizophrenia in both genetic IVs (i.e., GRS CRP and GRS GWAS ) and confirmed this protective relationship in followup analyses using individual-level genotype data from the schizophrenia GWAS. We also found a statistically significant association of CRP level with CAD, and nominally significant evidence for a predisposing causal association of CRP level with IBD, Crohn disease, psoriatic arthritis, knee osteoarthritis, SBP, DBP, eGFR cr , serum albumin level, serum protein level, and bipolar disorder, using GRS GWAS as an IV. However, after adjustment for heterogeneity, neither GRS showed a significant effect (at p < 0.0016) of CRP level on any of these outcomes, including CAD, nor on the 20 other common somatic and psychiatric outcomes we investigated, including celiac disease, ulcerative colitis, psoriasis (all types), rheumatoid arthritis, systemic lupus erythematous, systemic sclerosis, type 1 and 2 diabetes, stroke (all types), body mass index, chronic kidney disease, amyotrophic lateral sclerosis, Alzheimer disease, Parkinson disease, autism, and major depressive disorder.

CRP Protection against Schizophrenia
Strikingly, as opposed to the current literature and previous inconclusive small-scale studies [66][67][68], our findings suggest that genetically elevated levels of CRP are not predisposing but in fact protective for schizophrenia. The significant causal protective role of CRP for schizophrenia was consistent in both IVs using summary statistics, i.e., GRS CRP and GRS GWAS . When incorporating 18 genome-wide CRP-associated SNPs using individual-level data, we confirmed a modest, but significant, protective effect of CRP level for schizophrenia. This signal persisted when we included all SNPs meeting a less stringent p-value threshold of 1×10 −4 . Notably, the leave-one-out sensitivity analysis revealed that the genetic overlap between CRP level and schizophrenia we observed at genome-wide and 1×10 −4 significance thresholds was not driven by a few major SNPs. In contrast, others have previously shown that CRP levels are significantly elevated in patients with schizophrenia [69,70], with a recent meta-analysis concluding Effect size unit is millimeters of mercury per 1-mg/l increase in lnCRP. 3 Effect size unit is milliliters/minute/1.73 m 2 per 1-mg/l increase in lnCRP. 4 Effect size unit is grams/deciliter per 1-mg/l increase in lnCRP. M, number of markers used in the genetic instrument; p-het, p-value of heterogeneity of effect test. that the association between elevated CRP and schizophrenia is indeed robust [71]. Given that clinical studies report elevated CRP levels in schizophrenia, one would expect to find that alleles for elevated CRP would confer an increased risk for schizophrenia. The fact that we found the completely opposite effect-in a cohort of over 25,000 cases and 30,000 controlsshould give one pause when deriving clinical meaning from these results. Our observation that a genetically determined marginal increase in the level of CRP is likely to be protective for schizophrenia may fuel the debate about whether the observed CRP elevation in schizophrenia is a by-product of the pathogenesis of schizophrenia or directly contributing to clinical manifestations of the disorder [6]. Our finding may also point out potential biases in previous studies regarding the causes of elevated CRP levels in patients with schizophrenia, such as reverse causality and/or pleiotropic effects within chosen instruments. The exact mechanism for how elevated CRP levels are linked to schizophrenia requires a well-defined experimental analysis. In addition to CRP variants, other recent studies have identified several inflammatory genetic variants associated with schizophrenia and bipolar disorder, which include variants in the major histocompatibility complex (MHC) region on Chromosome 6p21 [72]-harboring many cytokine genes [54,[73][74][75][76]-and in the IL10 promoter [77], TNF promoter [78], IL1B [79], and C4 [80].

Biological Annotation
Following comments made by the reviewers, we explored the possible underlying pathways that may explain the potential protective causal association between CRP and schizophrenia. We performed a follow-up in silico functional pathway analysis using a previously reported approach [81] as summarized in S5 Methods and S4-S13 Tables. In brief, our results show that pathways associated with the interferon response are significantly enriched amongst genes harbored by CRP loci and their associated expression quantitative trait loci (eQTLs) and that there are differentially expressed genes between schizophrenia cases and controls. Previous studies showed that the induction of T cell IFN cytokine release stimulates microglia and astrocytes to facilitate glutamate clearance in neuronal cells without evoking inflammatory mediators [82,83]. One could speculate that CRP-interferon pathways may induce neuroprotection by contributing to glutamate clearance, leading to the protection of neurons against the oxidative stress associated with an excess of glutamate [84,85], and thereby offering a protective effect against schizophrenia.

CRP GRS GWAS Association with Bipolar Disorder
As for bipolar disorder, we found a nominal effect of a 1.21-fold increase in risk for bipolar disorder with a 10-s% increase in CRP level. Though this nominal predisposing effect needs to be confirmed, our finding corroborates epidemiological observations suggesting that elevated CRP is associated with the disease and supports a potential causal influence of general inflammation in bipolar disorder [86]. We note that, though it may be biologically sensible, this result failed to pass multiple testing correction. Confirmation by replication in independent cohorts, functional follow-up analyses, or the use of a stronger CRP GRS GWAS in upcoming studies is required to draw a definitive conclusion.

CRP GRS GWAS Association with Blood Pressure and Hypertension
We found nominally significant evidence for an up to~0.70-mm Hg increase in blood pressure with a 10-s% increase in CRP level and no evidence of heterogeneity for SBP. Additionally, there was nominally borderline significance for a causal association between CRP and DBP after adjustment for heterogeneity. These nominally significant findings, on the one hand, are in line with numerous epidemiological studies that have highlighted an association between elevated CRP and an increased risk of hypertension. For instance, one study found an association between CRP loci and hypertension in Asian individuals [87]. An additional line of support for a possible causal association of CRP and blood pressure comes from an experimental study in which an increase in CRP gene expression in mice, and subsequently CRP protein levels, led to a rise in SBP particularly [88]. Moreover, an ex vivo study by Zhou et al. showed that combining IL6 treatment and mechanical strain leads to a consistent increase in CRP expression at the protein and mRNA levels in smooth muscle cells [89]. Both inflammatory factors and local mechanical strains are abundant in blood vessels and are well-known risk factors for high blood pressure. Our finding did not reach a statistically significant level after correction for multiple testing; thus, it may echo previous MR studies that have failed to find a causal relationship between CRP level and blood pressure or hypertension in Europeans [90,91]. However, our systematic literature review showed that previous studies had some limitations (S1 Table). For instance, no study used a refined GWAS set of 18 CRP-associated SNPs; instead, they tested single or a limited set of CRP SNPs. Using such instruments might have led to biased estimates as their corresponding effects on CRP levels have been found to be small [30,57]. A combination of weak instruments and small sample sizes might have led to type II error [28,57] and hence to a conclusion of no causal association between CRP and blood pressure traits in previous studies. When all of the evidence is taken together, a direct link between CRP and blood pressure remains to be elucidated, though our nominal associations between GRS CRP and GRS GWAS and blood pressure do add to a line of findings from experimental studies suggesting a potential causal relationship between CRP and blood pressure.

CRP GRS GWAS Association with Osteoarthritis
Our nominally significant finding that CRP might be a potential causal factor for knee osteoarthritis (using GRS GWAS ) should be interpreted with caution. In line with our findings, we have previously shown that levels of CRP were higher in women with early radiological knee osteoarthritis (i.e., Kellgren-Lawrence grade 2+) and in women whose disease progressed [92]. Additionally, another study showed that genetically elevated CRP levels contribute to osteoarthritis severity [93]. However, other studies have found contrasting results [71,72,94]. One systematic review provided evidence that the relationship between CRP and osteoarthritis does exist but is dependent on body mass index [95]. It remains to be further investigated whether weight gain over the lifetime mediates the potential causal association between genetically elevated CRP and knee osteoarthritis.

CRP GRS GWAS Shows No Association with Other Remaining Outcomes
The present study was able to calculate nominal causal estimates for IBD, Crohn disease, psoriatic arthritis, CAD, eGFR cr , serum albumin level, and serum protein level using CRP GRS GWAS , but the estimates were altered by removal of SNPs from GRS GWAS based on heterogeneity tests, resulting in nominal or nonsignificant associations. These outcomes appeared therefore to have heterogeneity in the causal estimates, suggesting that these observed estimates were biased, likely due to pleiotropic effects of CRP loci. These results corroborate negative findings of previous studies (S1 Table), suggesting that a causal role of CRP in these traits and diseases is unlikely.

Methodological Concerns and Advantages
Pleiotropic biases in Mendelian randomization analyses using CRP GRS GWAS . A detailed evaluation of pleiotropic SNPs in our study showed that the method applied to identify heterogeneity sources was able to indicate and exclude several already known pleiotropic loci from the GRS GWAS IV. For instance, the use of a SNP in IL6R (rs4129267), amongst others, resulted in heterogeneity of effects on CAD risk. The same variant contributed to heterogeneity of effects for Crohn disease in our study, and it has been shown that this SNP is associated with levels of biomarkers other than CRP [56]. Further, a MR study found that IL6R SNPs, specifically the nonsynonymous SNP rs8192284, are associated with CAD risk and CRP levels [96]. Our selected IL6R SNPs, namely rs4537545 and rs4129267, are in extremely high linkage disequilibrium with rs8192284 (r 2 0.96 for both SNPs in HapMap data, CEU population). Carriers of the risk allele of rs8192284 have higher CRP, IL6, and fibrinogen levels [96]. Fibrinogen is also a well-known risk factor for CAD. Therefore, it is unclear so far which biomarker(s) mediates the effect of IL6R SNPs on CAD. Besides the IL6 locus, APOC1 and PABPC4 have been indicated as pleiotropic in three out of 32 our investigated outcomes, and PTPN2 and GCKR in six. With this information taken together, we were able to disentangle at least part of the pleiotropy regarding the causal estimates of CRP for outcomes. Again, we found no significant association of CRP GRS GWAS with IBD, Crohn disease, psoriatic arthritis, CAD, eGFR cr , serum albumin level, and serum protein level after adjustment for heterogeneity.
Using summary statistics of large-scale consortia. It is of utmost interest whether the observed effect of CRP as a risk predictor for human disease is causal, and thus whether reduction of CRP levels will lower the risk of disease. Here, we investigated the causality of CRP in 32 phenotypes by leveraging very large sample sizes collected by GWAS consortia, an approach that was much better powered than most previous MR studies. We found that genetically elevated CRP levels approximated by powerful instruments did not appear to contribute directly to most of the studied somatic and psychiatric outcomes. Our findings are consistent with previous MR studies reporting null associations of genetically elevated CRP levels with inflammation-related outcomes including CAD [56,59,97], type 2 diabetes [98], high body mass index [99], Alzheimer disease, and depression [100]. All previous MR studies were substantially limited to a single or a few outcomes, used only SNPs in the CRP gene, or had sample sizes much smaller than that of the present study (S1 Table). In addition to these studies, the current GWAS data do not corroborate epidemiological observations suggesting that elevated CRP levels are associated with amyotrophic lateral sclerosis [101], Alzheimer disease [102], Parkinson disease [103], and major depressive disorder [104]. Furthermore, patients with immunityrelated disorders frequently have a very high CRP level (as high as 100 mg/l) due to their disease status. Our findings may therefore more favorably indicate reverse causality. Taken together, these results show that CRP is highly unlikely to contribute causally to most of the major common somatic and neuropsychiatric outcomes that were investigated in the present study, with the possible exception of schizophrenia.
Strength of instrumental variables. The results presented in Table 2 show that our GRS CRP is not a weak instrument, as indicated by its high F-values owing to the large sample sizes of available outcomes from GWASs for the phenotypes under study. The strength of our instrument increased considerably in all disease classes when we used variants of multiple loci associated with CRP in GWASs. However, the variants comprising the CRP GRS GWAS explain on average only a moderate~5% of the total variance in baseline CRP levels [30]. Moreover, the possibility of effect modification by nongenetic CRP-related factors on the outcomes remains to be investigated. We may be able to create even stronger instruments based on ongoing efforts to identify additional variation influencing CRP levels. Even if larger sample sizes and stronger instruments can be realized, the overwhelming lack of causal effects observed for most outcomes in our study implies that therapies targeted at lowering CRP will not directly result in decreased risk of the investigated outcomes, or in better symptom management [105,106].
Using summary statistics instead of individual-level data. Here we used summary association statistics obtained from previously conducted meta-GWASs in order to maximize our study power. One may argue this may induce bias compared to when one uses individual-level data. Nevertheless, previous studies showed high agreement in results from MR methods using GWAS summary data and individual-level data [60,107]; Furthermore, our analyses of individual-level data for schizophrenia led to the same conclusion as our analyses using summary statistics data, confirming the robustness of our methodological approach.
Other potential sources of bias. An important rationale for MR is that the gene variants do not change over time and are inherited randomly. Thus, the genetic variants are considered free from confounding and reverse causation [108]. However, one cannot completely control for the possibility of confounding of genotype-intermediate phenotype-disease associations. For instance, there could be a confounding effect by ethnic/racial group (i.e., population stratification), but this is unlikely to be a major problem in most situations [108]. In the present study, we included summary statistics data from highly credible results of meta-GWASs. All the original meta-GWASs corrected for population stratification in cohort-level analyses and at meta-GWAS level.
Another caveat of MR is that developmental compensation might occur, through a genotype being expressed during fetal development that in turn buffers the effects of either environmental or genetic factors, a process called canalization [108,109]. Therefore, buffering mechanisms could hamper the associations between genetic variants and the outcome of interest. As opposed to this, a lifetime exposure to a risk factor may enhance its effects on the disease [109]. However, it is not clear to what extent genetically determined small changes in any given exposure would be sufficient to induce compensation [108].
All 32 of the meta-GWASs from which instrument summary estimates were taken were performed in individuals of European descent in Europe and the US and included thousands of samples for each outcome (S1 Table), which was also the case for our previous CRP meta-GWAS from which we chose the CRP-associated SNPs to calculate GRS GWAS . Therefore, the results of this MR study are applicable to individuals of European descent and are not necessarily generalizable to other ethnic groups.

Conclusion
We showed that elevated CRP levels driven by genetic factors are causally associated with protection against schizophrenia, suggesting that CRP may be one important puzzle piece that leads to an improved understanding of the pathogenesis of schizophrenia. We observed nominal evidence that genetically elevated CRP is causally associated with SBP, DBP, knee osteoarthritis, and bipolar disorder. Based on current GWAS data, we cannot verify any causal effect of CRP on the other 27 common somatic and neuropsychiatric outcomes investigated in the present study. Therefore, disease-associated rise in CRP levels may be a response to the disease process rather than a cause for these 27 outcomes. This implies that interventions to lower CRP levels are unlikely to result in decreased risk for the majority of common complex outcomes.   Table. List of genes at 108 genome-wide significant loci associated with schizophrenia. (XLSX) S11 Table. Brain and blood eQTL for credible sets of SNPs of the 108 schizophrenia loci. (XLSX) S12 Table. List of genes at 108 genome-wide significant loci associated with schizophrenia (yellow), brain eQTL (red), and blood eQTL (blue). (XLSX) S13 Table. Pathway enrichment results for the list of genes at 108 genome-wide significant loci associated with schizophrenia (yellow) and those associated with brain eQTL (red) and blood eQTL (blue). (XLSX)

Acknowledgments
We thank for Dan Arking for providing the genome-wide summary statistics data on CRP polymorphisms for the autism GWAS. We would also like to acknowledge the rheumatoid arthritis and systemic lupus erythematosus consortia, the Genetic Investigation of Anthropometric Traits (GIANT) consortium, and the PGC for making publicly available the data from which we extracted summary association statistics for CRP SNPs. The PAGE Consortium is grateful to Alex Tsoi for adding the consortium information and membership list.
The Type 1 Diabetes Genetics Consortium would like to thank all participants and their families, and also thank Smari Kristinsson, Larus Arni Hermannsson, and Asbjörn Krisbjörnsson of Raförninn for their extensive software design and contribution. We would like to thank the people with motor neuron disease and their families for their participation in this project. We acknowledge sample management undertaken by Biobanking Solutions funded by the Medical Research Council at the Centre for Integrated Genomic Medical Research, University of Manchester.

Consortia Coauthors
Please refer to S1 Consortia for the full names and affiliations of coauthors from the following five consortia: the Schizophrenia Working Group of the Psychiatric Genomics Consortium, the Autism Spectrum Disorder Working Group of the Psychiatric Genomics Consortium, the GERAD1 Consortium, the International Consortium for Blood Pressure Genome-Wide Association Studies, and the PAGE Consortium.

Editors' Summary
Background Inflammation is an important part of the human immune response, the network of cells and molecules that protects the body from attack by pathogens (infectious organisms) and from harmful substances and foreign particles (for example, splinters). When human cells are attacked by pathogens or injured by trauma or chemicals, molecules called inflammatory mediators induce fluid leakage from the blood vessels into the damaged tissue and attract "phagocytes" (a type of immune cell) to the site of infection or injury to "eat" the germs and dead or damaged cells. The end result is inflammation, which is characterized by swelling, redness, heat, and pain. The inflammatory response, although unpleasant, limits the damage caused by foreign invaders or chemicals by preventing further contact with body tissues. Sometimes, however, inflammation can be harmful. Persistent dysregulation of the inflammatory response is implicated in numerous somatic disorders (diseases that affect the body, such as cardiovascular disease) and neuropsychiatric disorders (mental disorders attributable to diseases of the nervous system, such as schizophrenia).

Why Was This Study Done?
Observational studies suggest that increased blood levels of C-reactive protein (CRP, an inflammatory protein) are associated with certain somatic and neuropsychiatric disorders. But observational studies cannot prove that changes in CRP levels actually cause any of these disorders. It could be that the individuals who develop a specific disease and who have a high CRP level also share another unknown characteristic that is actually responsible for disease development (confounding). Alternatively, it could be that the disease itself increases CRP levels (reverse causation). It is important to know whether CRP is causally involved in the development of specific diseases because it might then be possible to prevent or treat these diseases using drugs that control CRP levels. Here, the researchers undertake a Mendelian randomization study to determine whether CRP has a causal relationship with 32 common complex somatic and neuropsychiatric outcomes. Because gene variants are inherited randomly, they are not prone to confounding and are free from reverse causation. So, if CRP levels actually cause a specific somatic or neuropsychiatric disease, genetic variants that affect CRP levels should be associated with an altered risk for that disease.

What Did the Researchers Do and Find?
The researchers used data collected by several consortia involved in large genome-wide association studies (studies that ask whether specific genetic changes across the whole human genome, or blueprint, are associated with specific diseases) to look for associations between 32 somatic and neuropsychiatric outcomes and two genetic risk scores (GRSs) for CRP level. GRS CRP consisted of four single nucleotide polymorphisms (SNPs; a type of genetic variant) in the gene encoding CRP; GRS GWAS consisted of 18 SNPs that were associated with CRP level in a genome-wide association study. The researchers report that a genetically increased CRP level was significantly associated with a reduced risk of schizophrenia (a significant association is one unlikely to have arisen by chance). In addition, they found a nominally significant association (an association that needs to be confirmed) between genetically increased CRP levels and an increased risk of knee osteoarthritis, raised diastolic and systolic blood pressure, and bipolar disorder. Notably, there was no evidence for an effect of genetically increased CRP levels on any of the other 27 outcomes studied.
What Do These Findings Mean?
These findings suggest that genetically raised levels of CRP are causally associated with protection against schizophrenia, an unexpected finding given other recent studies that suggest that raised CRP levels and brain inflammation predispose individuals to schizophrenia. The findings also provide preliminary evidence that genetically raised levels of CRP may be causally associated with an increased risk of raised blood pressure, knee arthritis, and bipolar disorder. The lack of any association between genetically raised levels of CRP and the other outcomes studied suggests, however, that many previously identified disease-associated rises in CRP levels might be a response to the disease process rather than a cause of these diseases. Like all Mendelian randomization studies, the reliability of these findings depends on the validity of several assumptions made by the researchers and on the ability of the GRSs used in the study to explain variations in CRP level. Importantly, however, these findings suggest that interventions designed to lower CRP level are unlikely to decrease the risk of people developing the majority of common complex somatic and neuropsychiatric outcomes.

Additional Information
This list of resources contains links that can be accessed when viewing the PDF on a device or via the online version of the article at http://dx.doi.org/10.1371/journal.pmed.1001976.
• Wikipedia has pages on inflammation, C-reactive protein, and Mendelian randomization (note: Wikipedia is a free online encyclopedia that anyone can edit; available in several languages) • The MedlinePlus encyclopedia has a page on C-reactive protein (in English and Spanish) • The American Heart Association provides a short article on inflammation and heart disease • A UK National Health Service "Behind the Headlines" article explains a recent study that found an association between immune activity in the brain and schizophrenia