Introduction
Depression and anxiety are common contributors to population ill health. Clinical depression and anxiety disorders are responsible for an estimated 8.8% of the total years lived with disability globally (Global Burden of Disease Collaborative Network, 2019). Prevention of depression and anxiety requires sound knowledge about their aetiology. Elevated serum levels of C-reactive protein (CRP) and the biologically linked interleukin-6 (IL-6), as seen in low-grade inflammation, are among the many proposed causes, and observational studies have reported higher serum CRP levels among people with depression (Gimeno et al., Reference Gimeno, Kivimäki, Brunner, Elovainio, De Vogli, Steptoe and Ferrie2009; Howren, Lamkin, & Suls, Reference Howren, Lamkin and Suls2009; Miller & Raison, Reference Miller and Raison2016; Osimo, Baxter, Lewis, Jones, & Khandaker, Reference Osimo, Baxter, Lewis, Jones and Khandaker2019) and anxiety (Michopoulos, Powers, Gillespie, Ressler, & Jovanovic, Reference Michopoulos, Powers, Gillespie, Ressler and Jovanovic2017; Pitsavos et al., Reference Pitsavos, Panagiotakos, Papageorgiou, Tsetsekou, Soldatos and Stefanadis2006; Vogelzangs, Beekman, de Jonge, & Penninx, Reference Vogelzangs, Beekman, de Jonge and Penninx2013). However, it is uncertain whether these associations reflect causal effects of serum CRP, causality in the reverse direction (i.e. depression and anxiety causing increased serum CRP levels), or confounding due to common causes of elevated serum CRP and anxiety and depression.
Mendelian Randomisation (MR) may resolve the causal nature of these associations, by using genetically determined variation in serum CRP and serum IL-6 levels as instrumental variables to inform about their causal effects, while avoiding confounding and reverse causation bias (Haycock et al., Reference Haycock, Burgess, Wade, Bowden, Relton and Smith2016). However, MR studies have yielded inconsistent evidence, suggesting that higher genetically predicted serum CRP levels may either increase (Khandaker et al., Reference Khandaker, Zuber, Rees, Carvalho, Mason, Foley and Burgess2020) or have no effect (Kim Wium-Andersen, Dynnes Ørsted, & Grønne Nordestgaard, Reference Kim Wium-Andersen, Dynnes Ørsted and Grønne Nordestgaard2014) on the risk of depression. In apparent contrast to those results, a recent MR study in the UK Biobank population by Ye et al., yielded evidence that serum CRP may reduce, whereas serum IL-6 may increase anxiety and depression symptoms (Ye et al., Reference Ye, Kappelmann, Moser, Davey Smith, Burgess, Jones and Khandaker2021). However, UK Biobank has a participation rate of only 5.5% (Fry et al., Reference Fry, Littlejohns, Sudlow, Doherty, Adamska, Sprosen and Allen2017), and is therefore prone to collider bias (Munafò, Tilling, Taylor, Evans, & Davey Smith, Reference Munafò, Tilling, Taylor, Evans and Davey Smith2017). We aimed to investigate the effects of serum CRP and serum IL-6 on anxiety and depression symptoms and life satisfaction by using MR analyses in a study population less prone to collider bias, the population-based Trøndelag Health Study (HUNT) in Norway with 54% (HUNT3, 2006–08) to 69% (HUNT2, 1995–97) participation rates (Krokstad et al., Reference Krokstad, Langhammer, Hveem, Holmen, Midthjell, Stene and Holmen2013).
Method
Study design and study population
Using data from the HUNT Study, we estimated the cross-sectional associations of serum CRP with symptoms of depression and anxiety and life satisfaction and performed one-sample MR analyses to estimate the causal effects of serum CRP levels on these mental health outcomes. We used a two-sample MR approach to examine the corresponding causal effects of serum IL-6. For details about the HUNT study and its variables, see supplementary data, section ‘1. The HUNT study and variables’.
For our main analyses, we used similar genetic instruments for serum CRP as applied by Hartwig et al., in 2017 (Hartwig, Borges, Horta, Bowden, & Davey Smith, Reference Hartwig, Borges, Horta, Bowden and Davey Smith2017). This included two instruments: ‘CRP-conservative’ contains single nucleotide polymorphisms (SNPs) in the CRP gene that showed strong associations with serum CRP levels in a meta-analysis of 194 418 individuals (Wensley et al., Reference Wensley, Gao, Burgess, Kaptoge, Di Angelantonio, Shah and Danesh2011). ‘CRP-liberal’ includes all other SNPs that were genome-wide significantly associated with serum CRP levels in a meta-analysis of 80 000 participants of European descent (Dehghan et al., Reference Dehghan, Dupuis, Barbalic, Bis, Eiriksdottir, Lu and Chasman2011). We included the ‘CRP-conservative’ SNPs in our ‘CRP-liberal’ instrument. We put most emphasis on the results using the conservative instrument, as the restriction to the CRP gene makes horizontal pleiotropy less likely. The advantage of the liberal instrument is that the higher number of variants across the genome may increase statistical power, though at the cost of higher possibility of horizontal pleiotropy. As instruments for serum IL-6 levels we selected 3 SNPs within 55 kb of the IL6R gene that have been used in a recent MR study (Ye et al., Reference Ye, Kappelmann, Moser, Davey Smith, Burgess, Jones and Khandaker2021), originating from a 2012 genetic association study (Swerdlow et al., Reference Swerdlow, Holmes, Kuchenbaecker, Engmann, Shah, Sofat and Casas2012) where HUNT was not a part of the meta-analysis.
As part of the sensitivity analyses, we conducted one-sample MR with an alternative CRP-instrument (denoted as ‘CRP-sensitivity’). This instrument has recently been used in several MR studies (Georgakis et al., Reference Georgakis, Malik, Gill, Franceschini, Sudlow, Dichgans and Slagboom2020; Kappelmann et al., Reference Kappelmann, Arloth, Georgakis, Czamara, Rost, Ligthart and Binder2021; Ye et al., Reference Ye, Kappelmann, Moser, Davey Smith, Burgess, Jones and Khandaker2021) and contains 24 SNPs in the CRP region, based on a GWAS by Ligthart et al. (Ligthart et al. Reference Ligthart, Vaez, Võsa, Stathopoulou, de Vries, Prins and Alizadeh2018). For more details about SNPs and genetic instruments, see supplementary data, section ‘2. Genetic instruments’.
Among 78 962 participants in HUNT2 or HUNT3, we excluded 9751 participants without genotype information and 442 participants with no information on either serum CRP, Hospital Anxiety and Depression Scale (HADS) or life satisfaction, leaving 68 769 participants of European decent for the MR-analysis. However, since we calculated SNP-exposure and SNP-outcome associations separately the population size differs for each phenotype according to how many missing values there were for that particular phenotype, see the ‘n exposure’ and ‘n outcome’ column in online Supplementary Figs S3–S6. In the cross-sectional analysis, we required that all variables and covariates for each individual were collected at the same time point. Therefore we additionally excluded 29 935 participants that did not have data on serum CRP, HADS, life satisfaction, BMI, diabetes, alcohol use, cardiovascular disease and smoking status from the same survey (HUNT2 or HUNT3), leaving 38 834 participants for analysis.
Statistical methods
HADS and life satisfaction scores were natural log transformed due to skewed distributions, and serum CRP was log transformed with base 2 so that all outcomes are reported per doubling in CRP. We added 0.1 to all serum CRP values and 1 to all HADS values before transformation to avoid log(0) values. For participants with data from both HUNT2 and HUNT3 we used the earliest measurement, as these measurements were likely less influenced by confounding factors occurring more commonly in older age.
To estimate the cross-sectional associations of the measured serum CRP levels with HADS and life satisfaction, we used multivariable linear regression to estimate the percentage difference in HADS and life satisfaction scores per doubling in serum CRP, and used logistic regression to estimate the odds ratio (OR) of elevated HADS-A or HADS-D per doubling in serum CRP. We adjusted for age, sex, BMI, alcohol use, cardiovascular disease, smoking status and diabetes.
For the one-sample MR analysis we estimated the SNP-CRP and SNP-outcome associations using individual-level HUNT information on genotypes, serum CRP and HADS/life satisfaction. To account for relatedness in the HUNT population, we used a linear mixed effects model, with BOLT-LMM software (Loh et al., Reference Loh, Tucker, Bulik-Sullivan, Vilhjálmsson, Finucane, Salem and Price2015), to control for sex, birth year, genotyping batch and the first 4 genetic principal components. The output from BOLT-LMM represents the linear relationship between dosage of the individual SNP and the phenotype. Since the outcome-scores were log transformed, we converted the beta estimate to percentage change by the following formula:
For the binary phenotypes we used the following formula from the BOLM-LMM User manual to approximate the log (OR):
The standard errors for the log (OR) were approximated by the following formula:
For the one-sample MR we used the SNP-exposure and SNP-outcome associations from the BOLT-LMM software and calculated inverse variance weighted (IVW) MR estimates using the ‘TwoSampleMR’ package (Hemani et al., Reference Hemani, Zheng, Elsworth, Wade, Haberland, Baird and Haycock2018) version 0.5.6, combined with the ‘MendelianRandomisation’ package (Olena Yavorska, Reference Olena Yavorska2021) version 0.5.1, which allowed us to account for linkage disequilibrium between the SNPs.
For the two-sample MR analyses we used SNP-IL-6 associations as reported in Ye et al. (Ye et al. Reference Ye, Kappelmann, Moser, Davey Smith, Burgess, Jones and Khandaker2021), and SNP-outcome associations from the BOLT-LMM analysis using individual genotype and outcome information from HUNT. In the same manner as for the one-sample MR we combined the ‘TwoSampleMR’ package (Hemani et al., Reference Hemani, Zheng, Elsworth, Wade, Haberland, Baird and Haycock2018) and the ‘MendelianRandomisation’ package (Olena Yavorska, Reference Olena Yavorska2021) to obtain IVW MR estimates accounting for linkage disequilibrium.
Except for the mixed effect model analysis, all statistical analyses were performed in R (R Core Team, 2021) version 4.0.5. Mixed effect model analysis was performed with the BOLT-LMM software (Loh et al., Reference Loh, Tucker, Bulik-Sullivan, Vilhjálmsson, Finucane, Salem and Price2015) (version 2.3.4).
Assessment of MR assumptions and other sensitivity analyses
MR builds on three main assumptions: (1) the relevance assumption that the genetic instrument is associated with the exposure; (2) the exclusion restriction assumption that there should be no other mechanism in which the SNP affects the outcome other than through the exposure (a violation of which would cause horizontal pleiotropy), and (3) the independence assumption that the SNP should be independent from any confounder of the exposure-outcome association. To evaluate the plausibility of these assumptions we used several different methods including calculation of F-statistics and R 2, Phenoscanner, regression analyses between SNPs and possible confounders, Cochran Q-statistic, MR Egger, median weighted MR, weighted mode MR and MR with an alternative genetic instrument. For details on how this was conducted, see supplementary data, section ‘3. Method for assessment of MR assumptions and other sensitivity analyses’.
Results
Descriptive data
The characteristics of the study population are described in Table 1.
Table 1. Characteristics of study populations
IQR, interquartile range; s.d., standard deviation; n, number.
Cross-sectional analysis
Adjusted for age, sex, BMI, alcohol use, cardiovascular disease, smoking status and diabetes, each doubling in serum CRP was associated with −0.14% difference in HADS-T (95% CI −0.64 to 0.37, p value: 0.599), 0.27% difference in HADS-D (95% CI −0.20 to 0.75; p value: 0.260) and −0.77% difference in HADS-A (95% CI −1.24 to −0.29; p value: 0.002) (Fig. 1). For binary outcomes, each doubling in serum CRP was associated with an OR of 1.03 (95% CI 0.996 to 1.06; p value: 0.085) for HADS-D ≥ 8 and an OR of 0.996 (95% CI 0.97 to 1.02, p value: 0.701) for HADS-A ≥ 8 (Fig. 2). There was a difference of −0.10% in life satisfaction score per doubling in CRP (95% CI −0.41 to 0.21; p value: 0.523) (Fig. 1).
Figure 1. Results from cross-sectional analysis with continuous outcomes. Adjusted for age, sex, BMI, alcohol use, cardiovascular disease, smoking status and diabetes. CRP, C-reactive protein; HADS, hospital anxiety and depression scale; HADS-A, HADS-anxiety score; HADS-D, HADS-depression score; HADS-T, total HADS score; 95% CI, 95% confidence interval; n, number of participants.
Figure 2. Results from cross-sectional analysis with binary outcomes. Adjusted for age, sex, BMI, alcohol use, cardiovascular disease, smoking status and diabetes. CRP, C-reactive protein; HADS, hospital anxiety and depression scale; HADS-A, HADS-anxiety score; HADS-D, HADS-depression score; HADS-T, total HADS score; 95% CI, 95% confidence interval; n, number of participants.
One-sample MR analysis of CRP serum levels
In one-sample MR analysis using the conservative genetic instrument for CRP, each doubling in genetically predicted serum CRP was associated with 2.63% (95% CI −0.02 to 5.35, p value: 0.051) higher HADS-T score, 2.43% (95% CI −0.11 to 5.03, p value: 0.061) higher HADS-D score, 1.94% (95% CI −0.58 to 4.52, p value: 0.132) higher HADS-A score, and 2.00% (95% CI 0.45–3.59, p value: 0.011) higher life satisfaction score. The associations were in the same direction, but weaker for the liberal instrument (Fig. 3). A doubling in genetically predicted serum CRP gave an OR of 0.94 (95% CI 0.82–1.08, p value: 0.392) for HADS-D ≥ 8, and 1.03 (95% CI 0.92–1.15, p value: 0.591) for HADS-A ≥ 8 (Fig. 4).
Figure 3. Results from one-sample MR with continuous outcomes. CRP, C-reactive protein; CRP-lib, CRP-liberal genetic instrument; CRP-con, CRP-conservative genetic instrument; SNP, single nucleotide polymorphism; HADS, hospital anxiety and depression scale; HADS-A, HADS-anxiety score; HADS-D, HADS-depression score; HADS-T, total HADS score; 95% CI, 95% confidence interval; p-Qstat, p value for Q-statistic; n exposure, number of participants in SNP-exposure analysis; n outcome, number of participants in SNP-outcome analysis.
Figure 4. Results from one-sample MR with binary outcomes. CRP, C-reactive protein; CRP-lib, CRP-liberal genetic instrument; CRP-con, CRP-conservative genetic instrument; SNP, single nucleotide polymorphism; HADS, hospital anxiety and depression scale; HADS-A, HADS-anxiety score; HADS-D, HADS-depression score; 95% CI, 95% confidence interval; p-Qstat, p value for Q-statistic; n exposure, number of participants in SNP-exposure analysis; n outcome, number of participants in SNP-outcome analysis.
Two-sample MR analysis of IL-6 serum levels
Genetically predicted higher serum IL-6 were associated with slightly lower HADS-T, HADS-D, HADS-A and life satisfaction scores, but all estimates were hampered by large imprecision (Fig. 5). For clinical depression and anxiety, signs of an inverse association were only seen for depression (Fig. 6); however, no association with IL-6 reached conventional criteria for statistical significance.
Figure 5. Results from two-sample MR with continuous outcomes. CRP, C-reactive protein; IL-6, interleukin-6; SNP, single nucleotide polymorphism; HADS, hospital anxiety and depression scale; HADS-A, HADS-anxiety score; HADS-D, HADS-depression score; HADS-T, total HADS score; 95% CI, 95% confidence interval; p-Qstat, p value for Q-statistic; n exposure, number of participants in SNP-exposure analysis; n outcome, number of participants in SNP-outcome analysis.
Figure 6. Results from two-sample MR with binary outcomes. CRP, C-reactive protein; IL-6, interleukin-6; SNP, single nucleotide polymorphism; HADS, hospital anxiety and depression scale; HADS-A, HADS anxiety score; HADS-D, HADS depression score; 95% CI, 95% confidence interval; p-Qstat, p value for Q-statistic; n exposure, number of participants in SNP-exposure analysis; n outcome, number of participants in SNP-outcome analysis.
Assessment of assumptions
Relevance assumption
F-statistics for the SNP-exposure associations were >214 for all SNPs included in the conservative CRP instrument, >10 for 14 of the 18 SNPs additionally included in the liberal CRP instrument, and >22 for all SNPs used in the IL-6 instrument (online Supplementary Tables S5–S7). This indicates that 4 of the SNPs in the liberal instrument were too weak, which in a one-sample setting would bias the MR estimate towards the cross-sectional estimate (Davies, Holmes, & Davey Smith, Reference Davies, Holmes and Davey Smith2018).
To investigate this further, we ran linear regression with the allele count for the conservative and the liberal CRP-instrument. This yielded an adjusted R 2 of 0.019 for the conservative and 0.038 for the liberal instrument. The F-statistics were 989.9 and 2025.2 respectively (online Supplementary Table S8). This indicates that even if some SNPs may be too weak, the general strength of the liberal instrument is satisfactory and the risk for weak instrument bias is small.
Exclusion restriction assumption and independence assumption
Using the R-package Phenoscanner (Kamat et al., Reference Kamat, Blackshaw, Young, Surendran, Burgess, Danesh and Staley2019; Staley et al., Reference Staley, Blackshaw, Kamat, Ellis, Surendran, Sun and Butterworth2016) to search through the NHGRI-EBI GWAS-catalog (Buniello et al., Reference Buniello, MacArthur, Cerezo, Harris, Hayhurst, Malangone and Parkinson2019), we found no genome-wide significant associations between the four SNPs in the conservative CRP instrument and other phenotypes than CRP (online Supplementary Table S9). However, many of the 18 additional SNPs in the liberal CRP instrument were associated with traits that could be potential confounders or pleiotropic factors, in particular CVD and CVD risk factors (online Supplementary Table S10). Among SNPs used in the instrument for IL-6, we found a few associations with potential confounders such as CVD, eczema and rheumatoid arthritis (online Supplementary Table S11). These associations may reflect that IL-6 may increase the risk of these diseases, which would not imply a violation of the independence assumption. Analyses within the HUNT study population showed no strong associations between the allele scores for the genetic instruments for CRP and the potential pleiotropic factors such as BMI, smoking status, or a history of myocardial infarction, heart failure or rheumatoid arthritis (online Supplementary Tables S12 and S13). For the IL-6 instrument there were no statistically significant associations between allele score and likely confounders or pleiotropic factors (online Supplementary Table S14).
The MR-Egger, weighted median MR and weighted mode MR analyses were quite consistent with the IVW analysis in that they did not point to large effects of serum CRP on the outcomes (online Supplementary Figs S1 and S2). The MR Egger test (pleiotropy p value) did not indicate any large directional pleiotropy. The Q-statistics for the conservative and liberal CRP instruments (online Supplementary Figs S3 and S4) indicated no strong heterogeneity, except for the liberal CRP instrument with HADS-D as outcome, with a Q-statistic of 39.48 (p value 0.009) suggesting a violation of the exclusion restriction assumption. Since there only were three SNPs in the IL-6 instrument, we did not run MR-Egger, weighted median or weighted mode MR for IL-6. The Q-statistics for the IL-6 SNPs were relatively small, all below 5 with p values ≥ 0.120, suggesting little heterogeneity (Figs 5 and 6).
Sensitivity analysis
The results from IVW MR using the alternative genetic instrument from Ligthart et al. (Ligthart et al. Reference Ligthart, Vaez, Võsa, Stathopoulou, de Vries, Prins and Alizadeh2018) combined with the CRP-conservative instrument were consistent with the main results (online Supplementary Figs S3 and S4). The Q-statistics for this instrument showed no strong signs of heterogeneity and MR-Egger did not indicate any problems with pleiotropy with all pleiotropy p values above 0.05 (online Supplementary Figs S5 and S6). However, with F-statistics below 10 for 5 of the 26 SNPs and an overall F-statistic for the entire instrument of 723, this instrument was weaker than the CRP liberal instrument in HUNT (online Supplementary Tables S15 and S16).
Discussion
Key results
We observed weak cross-sectional associations where each doubling in serum CRP was associated with approximately 0.3% higher HADS-D score and −0.1% difference in life satisfaction score (where lower scores indicate higher life satisfaction), but −0.8% difference in HADS-A score. Only the association with HADS-A reached conventional criteria for statistical significance. MR analyses provided some evidence that serum CRP may modestly increase anxiety and depression symptoms and reduce life satisfaction, with ~2% higher scores per doubling in serum CRP, but only the association with life satisfaction reached conventional criteria for statistical significance. For IL-6, causal estimates were in the opposite direction, but imprecise and far from conventional criteria for statistical significance.
Limitations
Our sensitivity analyses indicate that the core MR assumptions are held for the conservative CRP instrument, whereas associations to possible confounders and signs of heterogeneity suggest that the estimates for the liberal CRP instrument should be interpreted with more caution. Applying a different set of SNPs from a recent GWAS (Ligthart et al., Reference Ligthart, Vaez, Võsa, Stathopoulou, de Vries, Prins and Alizadeh2018), the CRP-sensitivity instrument revealed results consistent with the liberal CRP instrument. This instrument showed no signs of heterogeneity or pleiotropy, but had weaker F-statistics and R 2 than the CRP liberal instrument. Weak instrument bias could potentially bias the estimates for the CRP liberal and CRP sensitivity instruments towards the cross-sectional estimate (Davies et al., Reference Davies, Holmes and Davey Smith2018). Future studies could use a two-sample design to avoid this, but also to gain more statistical power to precisely estimate weak causal estimates like those we observed. Another limitation of this study is that the measurement of life satisfaction only included a single question, making it difficult to evaluate what aspects of life satisfaction were captured.
Interpretations
CRP has been suggested to influence mental health in previous studies. Even though CRP does not pass the blood brain barrier freely, CRP could affect the central nervous system via passage through leaky regions, active uptake, activation of endothelial cells and immunity cells or binding to receptors on peripheral nerves (Felger et al., Reference Felger, Haroon, Patel, Goldsmith, Wommack, Woolwine and Miller2020). Our findings, however, do not support a major causal role of serum CRP as a determinant of anxiety and depression symptoms and life satisfaction. Nonetheless, our MR analyses provide weak evidence that serum CRP may modestly increase anxiety and depression symptoms and lower life satisfaction, which would be in agreement with conventional observational studies showing serum CRP to be elevated in a quarter of patience with depression (Osimo et al., Reference Osimo, Baxter, Lewis, Jones and Khandaker2019) and serum CRP to be elevated in subgroups of patients with anxiety disorders (Vogelzangs et al., Reference Vogelzangs, Beekman, de Jonge and Penninx2013). Also, such a causal effect would be in agreement with the hypothesis that inflammation may increase sickness behaviour and fatigue which again would affect mental health (Miller & Raison, Reference Miller and Raison2016).
Among MR studies, Khandaker et al. (Khandaker et al. Reference Khandaker, Zuber, Rees, Carvalho, Mason, Foley and Burgess2020) found increased risk of lifetime major depression at higher genetically predicted levels of serum CRP and IL-6 in UK Biobank. In contrast, Kappelman et al. (Kappelmann et al. Reference Kappelmann, Arloth, Georgakis, Czamara, Rost, Ligthart and Binder2021) and Wium Andersen et al. (Kim Wium-Andersen et al. Reference Kim Wium-Andersen, Dynnes Ørsted and Grønne Nordestgaard2014) found no effect of serum CRP on depression, whereas Perry et al. (Perry et al. Reference Perry, Upthegrove, Kappelmann, Jones, Burgess and Khandaker2021) found no effect of serum CRP, but evidence of an effect of higher serum IL-6 on depression. In opposite direction of our findings, evidence from the recent UK Biobank MR study by Ye et al. (Ye et al. Reference Ye, Kappelmann, Moser, Davey Smith, Burgess, Jones and Khandaker2021) suggested that serum CRP may lower, while serum IL-6 may increase anxiety and depression symptoms. The reasons for the discrepancy between our results and those of Ye et al., are not known, but the low participation rate in UK Biobank could give rise to collider bias (Munafò et al., Reference Munafò, Tilling, Taylor, Evans and Davey Smith2017). Ye et al., comprehensively evaluated the risk for potential collider bias with inverse probability weighting techniques and found no evidence of bias. However, the ability of such analyses to detect and correct for selection bias depends on how accurately the probability of selection can be modelled (Gkatzionis & Burgess, Reference Gkatzionis and Burgess2019). We also note that Ye et al., assessed symptoms using the Patient Health Questionnaire (PHQ-9) and the Generalised Anxiety Disorder (GAD-7) Questionnaire, whereas we used HADS-D and HADS-A. PHQ-9 assesses 9 symptoms from the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) criteria for depression, while GAD-7 asseses 7 symptoms related to general anxiety disorder based on the DSM-IV criteria (Spitzer, Kroenke, Williams, & Löwe, Reference Spitzer, Kroenke, Williams and Löwe2006). One important difference between these questionaires and HADS, is that HADS does not meassure somatic symptoms of depression and anxiety such as dizziness, headaches, insomnia and fatigue. HADS also excludes symptoms related to serious mental disorders (Bjelland, Dahl, Haug, & Neckelmann, Reference Bjelland, Dahl, Haug and Neckelmann2002). However, studies show that HADS-D and HADS-A are comparable with PHQ-9 and GAD-7 in their psychometric properties (Esser et al., Reference Esser, Hartung, Friedrich, Johansen, Wittchen, Faller and Mehnert2018; Hansson, Chotai, Nordstöm, & Bodlund, Reference Hansson, Chotai, Nordstöm and Bodlund2009). Hence, comparing our results from HADS with PHQ-9 and GAD-7 is sensible.
Generalisability
Our study is based on a Norwegian population from Trøndelag County with low ethnic diversity. However, it seems unlikely that serum CRP and serum IL-6 should have major different effects on anxiety and depression symptoms in this group compared with other populations.
Conclusion
In conclusion our findings do not support a major causal role of serum CRP on anxiety and depression symptoms and life satisfaction, but provide weak evidence that serum CRP may modestly decrease life satisfaction and increase anxiety and depression symptoms. For IL-6, causal estimates were in the opposite direction, but impresise and therefore inconclusive. Our results do not support the recent suggestion that serum CRP may lower anxiety and depression symptoms.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291723001290.
Acknowledgements
The Trøndelag Health Study (HUNT) is a collaboration between HUNT Research Centre (Faculty of Medicine and Health Sciences, Norwegian University of Science and Technology NTNU), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health.
Financial support
Department of Public Health and Nursing, Faculty of Medicine and Health Sciences, Norwegian University of Science and Technology (NTNU).
Conflict of interest
None.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The study was ethically approved by the Regional Committee for Medical and Health Research Ethics in Mid-Norway (REK# 2018/2487). All study participants have given informed consent.
Open access
