Genome-wide association study of circulating interleukin 6 levels identifies novel loci

Effects of genetically predicted posttraumatic stress disorder on autoimmune phenotypes

Observational studies suggest that posttraumatic stress disorder (PTSD) increases risk for various autoimmune diseases. Insights into shared biology and causal relationships between these diseases may inform intervention approaches to PTSD and co-morbid autoimmune conditions. We investigated the shared genetic contributions and causal relationships between PTSD, 18 autoimmune diseases, and 3 immune/inflammatory biomarkers. Univariate MiXeR was used to contrast the genetic architectures of phenotypes. Genetic correlations were estimated using linkage disequilibrium score regression. Bi-directional, two-sample Mendelian randomization (MR) was performed using independent, genome-wide significant single nucleotide polymorphisms; inverse variance weighted and weighted median MR estimates were evaluated. Sensitivity analyses for uncorrelated (MR PRESSO) and correlated horizontal pleiotropy (CAUSE) were also performed. PTSD was considerably more polygenic (10,863 influential variants) than autoimmune diseases (median 255 influential variants). However, PTSD evidenced significant genetic correlation with nine autoimmune diseases and three inflammatory biomarkers. PTSD had putative causal effects on autoimmune thyroid disease (p = 0.00009) and C-reactive protein (CRP) (p = 4.3 × 10-7). Inferences were not substantially altered by sensitivity analyses. Additionally, the PTSD-autoimmune thyroid disease association remained significant in multivariable MR analysis adjusted for genetically predicted inflammatory biomarkers as potential mechanistic pathway variables. No autoimmune disease had a significant causal effect on PTSD (all p values > 0.05). Although causal effect models were supported for associations of PTSD with CRP, shared pleiotropy was adequate to explain a putative causal effect of CRP on PTSD (p = 0.18). In summary, our results suggest a significant genetic overlap between PTSD, autoimmune diseases, and biomarkers of inflammation. PTSD has a putative causal effect on autoimmune thyroid disease, consistent with existing epidemiologic evidence. A previously reported causal effect of CRP on PTSD is potentially confounded by shared genetics. Together, results highlight the nuanced links between PTSD, autoimmune disorders, and associated inflammatory signatures, and suggest the importance of targeting related pathways to protect against disease and disability.

Causal effect of interleukin (IL)-6 on blood pressure and hypertension: A mendelian randomization study

To examine whether circulating interleukin-6 (IL-6) levels (CirIL6) have a causal effect on blood pressure using Mendelian randomization (MR) methods. We used data from genome-wide association studies (GWAS) of European ancestry to obtain genetic instruments for circulating IL-6 levels and blood pressure measurements. We applied several robust MR methods to estimate the causal effects and to test for heterogeneity and pleiotropy. We found that circulating IL-6 had a significant positive causal effect on systolic blood pressure (SBP) and pulmonary arterial hypertension (PAH), but not on diastolic blood pressure (DBP) or hypertension. We found that as CirIL6 genetically increased, SBP increased using Inverse Variance Weighted (IVW) method (for ukb-b-20175, β = 0.082 with SE = 0.032, P = 0.011; for ukb-a-360, β = 0.075 with SE = 0.031, P = 0.014) and weighted median (WM) method (for ukb-b-20175, β = 0.061 with SE = 0.022, P = 0.006; for ukb-a-360, β = 0.065 with SE = 0.027, P = 0.014). Moreover, CirIL6 may be associated with an increased risk of PAH using WM method (odds ratio (OR) = 15.503, 95% CI, 1.025-234.525, P = 0.048), but not with IVW method. Our study provides novel evidence that circulating IL-6 has a causal role in the development of SBP and PAH, but not DBP or hypertension. These findings suggest that IL-6 may be a potential therapeutic target for preventing or treating cardiovascular diseases and metabolic disorders. However, more studies are needed to confirm the causal effects of IL-6 on blood pressure and to elucidate the underlying mechanisms and pathways.

A genetic exploration of the relationship between Posttraumatic Stress Disorder and cardiovascular diseases

Background and Aims

Experiencing a traumatic event may lead to Posttraumatic Stress Disorder (PTSD), including symptoms such as flashbacks and hyperarousal. Individuals suffering from PTSD are at increased risk of cardiovascular disease (CVD), but it is unclear why. This study assesses shared genetic liability and potential causal pathways between PTSD and CVD.

Methods

We leveraged summary-level data of genome-wide association studies (PTSD: N= 1,222,882; atrial fibrillation (AF): N=482,409; coronary artery disease (CAD): N=1,165,690; hypertension: N=458,554; heart failure (HF): N=977,323). First, we estimated genetic correlations and utilized genomic structural equation modeling to identify a common genetic factor for PTSD and CVD. Next, we assessed biological, behavioural, and psychosocial factors as potential mediators. Finally, we employed multivariable Mendelian randomization to examine causal pathways between PTSD and CVD, incorporating the same potential mediators.

Results

Significant genetic correlations were found between PTSD and CAD, HT, and HF (rg =0.21-0.32, p≤ 3.08 · 10-16), but not between PTSD and AF. Insomnia, smoking, alcohol dependence, waist-to-hip ratio, and inflammation (IL6, C-reactive protein) partly mediated these associations. Mendelian randomization indicated that PTSD causally increases CAD (IVW OR=1.53, 95% CIs=1.19-1.96, p=0.001), HF (OR=1.44, CIs=1.08-1.92, p=0.012), and to a lesser degree hypertension (OR=1.25, CIs=1.05-1.49, p=0.012). While insomnia, smoking, alcohol, and inflammation were important mediators, independent causal effects also remained.

Conclusions

In addition to shared genetic liability between PTSD and CVD, we present strong evidence for causal effects of PTSD on CVD. Crucially, we implicate specific lifestyle and biological mediators (insomnia, substance use, inflammation) which has important implications for interventions to prevent CVD in PTSD patients.

Genetic associations and potential mediators between psychiatric disorders and irritable bowel syndrome: a Mendelian randomization study with mediation analysis

Objective

Potential causal associations between psychiatric disorders and irritable bowel syndrome have been demonstrated in observational studies; however, these studies are susceptible to underlying confounding and reverse causation biases. We aimed to assess the causal effects of psychiatric disorders on irritable bowel syndrome (IBS) and the potential mediators from a genetic perspective by conducting a Mendelian randomization (MR) study with mediation analysis.

Method

Genetic instruments associated with psychiatric disorders, potential mediators, and IBS were obtained from large-scale genome-wide association studies (GWAS). Three MR methods - the inverse-variance weighted (IVW) method, MR-Egger method, and weighted median method, were used to investigate causal association estimates. Heterogeneity among different genetic instrumental variables (IVs) was assessed using Q tests. Additionally, the MR-PRESSO and MR-Pleiotropy methods were used to verify horizontal pleiotropy and detect outliers that might bias the results, which were removed from further analysis. Consequently, we used MR mediation analysis to investigate potential mediators in the causal associations between psychiatric disorders and IBS.

Results

MR provided evidence of the causal effects of genetically predicted broad depression, major depressive disorder (MDD), anxiety disorder, post-traumatic stress disorder (PTSD), and schizophrenia on IBS. The results of MR mediation analysis demonstrated that the reduction in acetate levels mediated 12.6% of the effects of broad depression on IBS; insomnia mediated 16.00%, 16.20%, and 27.14% of the effects of broad depression, MDD, and PTSD on IBS, respectively; and the increase in blood β-hydroxybutyrate levels mediated 50.76% of the effects of schizophrenia on IBS.

Conclusion

Our study confirmed the brain-gut axis involvement and potential modulators in the pathophysiology of psychiatric disorder-induced IBS from a genetic perspective, and suggests potential therapeutic targets for the disrupted brain-gut axis.

Distinct genomic signatures and modifiable risk factors underly the comorbidity between major depressive disorder and cardiovascular disease

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Using genomic data, this study elucidates biological mechanisms, key risk factors, and causal pathways underlying their comorbidity. We show that CVDs share a large proportion of their genetic risk factors with MDD. Multivariate genome-wide association analysis of the shared genetic liability between MDD and atherosclerotic CVD (ASCVD) revealed seven novel loci and distinct patterns of tissue and brain cell-type enrichments, suggesting a role for the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic, and psychosocial/lifestyle risk factors. Finally, we found support for causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and demonstrated that the causal effects were partly explained by metabolic and psychosocial/lifestyle factors. The distinct signature of MDD-ASCVD comorbidity aligns with the idea of an immunometabolic sub-type of MDD more strongly associated with CVD than overall MDD. In summary, we identify plausible biological mechanisms underlying MDD-CVD comorbidity, as well as key modifiable risk factors for prevention of CVD in individuals with MDD.

SumVg: Total Heritability Explained by All Variants in Genome-Wide Association Studies Based on Summary Statistics with Standard Error Estimates

Genome-wide association studies (GWAS) are commonly employed to study the genetic basis of complex traits/diseases, and a key question is how much heritability could be explained by all single nucleotide polymorphisms (SNPs) in GWAS. One widely used approach that relies on summary statistics only is linkage disequilibrium score regression (LDSC); however, this approach requires certain assumptions about the effects of SNPs (e.g., all SNPs contribute to heritability and each SNP contributes equal variance). More flexible modeling methods may be useful. We previously developed an approach recovering the "true" effect sizes from a set of observed z-statistics with an empirical Bayes approach, using only summary statistics. However, methods for standard error (SE) estimation are not available yet, limiting the interpretation of our results and the applicability of the approach. In this study, we developed several resampling-based approaches to estimate the SE of SNP-based heritability, including two jackknife and three parametric bootstrap methods. The resampling procedures are performed at the SNP level as it is most common to estimate heritability from GWAS summary statistics alone. Simulations showed that the delete-d-jackknife and parametric bootstrap approaches provide good estimates of the SE. In particular, the parametric bootstrap approaches yield the lowest root-mean-squared-error (RMSE) of the true SE. We also explored various methods for constructing confidence intervals (CIs). In addition, we applied our method to estimate the SNP-based heritability of 12 immune-related traits (levels of cytokines and growth factors) to shed light on their genetic architecture. We also implemented the methods to compute the sum of heritability explained and the corresponding SE in an R package SumVg. In conclusion, SumVg may provide a useful alternative tool for calculating SNP heritability and estimating SE/CI, which does not rely on distributional assumptions of SNP effects.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Glycoprotein acetyls and depression: Testing for directionality and potential causality using longitudinal data and Mendelian randomization analyses

BACKGROUND

Inflammation is associated with depression, but causality remains unclear. We investigated potential causality and direction of effect between inflammation and depression.

METHODS

Using data from the ALSPAC birth cohort (n = 4021; 42.18 % male), we used multivariable regression to investigate bidirectional longitudinal associations of GlycA and depression and depression symptoms, assessed at ages 18y and 24y. We used two-sample Mendelian randomization (MR) to investigate potential causality and directionality. Genetic variants for GlycA were obtained from UK Biobank (UKB) (N = 115,078); for depression from the Psychiatric Genomics Consortium and UKB (N = 500,199); and for depressive symptoms (N = 161,460) from the Social Science Genetic Association Consortium. In addition to the Inverse Variance Weighted method, we used sensitivity analyses to strengthen causal inference. We conducted multivariable MR adjusting for body mass index (BMI) due to known genetic correlation between inflammation, depression and BMI.

RESULTS

In the cohort analysis, after adjusting for potential confounders we found no evidence of associations between GlycA and depression symptoms score or vice versa. We observed an association between GlycA and depression (OR = 1∙18, 95 % CI: 1∙03-1∙36). MR suggested no causal effect of GlycA on depression, but there was a causal effect of depression on GlycA (mean difference in GlycA = 0∙09; 95 % CI: 0∙03-0∙16), which was maintained in some, but not all, sensitivity analyses.

LIMITATIONS

The GWAS sample overlap could incur bias.

CONCLUSION

We found no consistent evidence for an effect of GlycA on depression. There was evidence that depression increases GlycA in the MR analysis, but this may be confounded/mediated by BMI.

The proteome and phosphoproteome of circulating extracellular vesicle-enriched preparations are associated with characteristic clinical features in type 1 diabetes

Introduction

There are no validated clinical or laboratory biomarkers to identify and differentiate endotypes of type 1 diabetes (T1D) or the risk of progression to chronic complications. Extracellular vesicles (EVs) have been studied as biomarkers in several different disease states but have not been well studied in T1D.

Methods

As the initial step towards circulating biomarker identification in T1D, this pilot study aimed to provide an initial characterization of the proteomic and phosphoproteomic landscape of circulating EV-enriched preparations in participants with established T1D (N=10) and healthy normal volunteers (Controls) (N=7) (NCT03379792) carefully matched by age, race/ethnicity, sex, and BMI. EV-enriched preparations were obtained using EVtrap® technology. Proteins were identified and quantified by LC-MS analysis. Differential abundance and coexpression network (WGCNA), and pathway enrichment analyses were implemented.

Results

The detected proteins and phosphoproteins were enriched (75%) in exosomal proteins cataloged in the ExoCarta database. A total of 181 proteins and 8 phosphoproteins were differentially abundant in participants with T1D compared to controls, including some well-known EVproteins (i.e., CD63, RAB14, BSG, LAMP2, and EZR). Enrichment analyses of differentially abundant proteins and phosphoproteins of EV-enriched preparations identified associations with neutrophil, platelet, and immune response functions, as well as prion protein aggregation. Downregulated proteins were involved in MHC class II signaling and the regulation of monocyte differentiation. Potential key roles in T1D for C1q, plasminogen, IL6ST, CD40, HLA-DQB1, HLA-DRB1, CD74, NUCB1, and SAP, are highlighted. Remarkably, WGCNA uncovered two protein modules significantly associated with pancreas size, which may be implicated in the pathogenesis of T1D. Similarly, these modules showed significant enrichment for membrane compartments, processes associated with inflammation and the immune response, and regulation of viral processes, among others.

Discussion

This study demonstrates the potential of proteomic and phosphoproteomic signatures of EV-enriched preparations to provide insight into the pathobiology of T1D. The WGCNA analysis could be a powerful tool to discriminate signatures associated with different pathobiological components of the disease.

Assessing the causal relationship between genetically determined inflammatory biomarkers and low back pain risk: a bidirectional two-sample Mendelian randomization study

Background

Observational studies have suggested an association between inflammatory markers and low back pain (LBP), but the causal relationship between these factors remains uncertain.

Methods

We conducted a bidirectional two-sample Mendelian randomization analysis (MR) study to investigate whether there is a causal relationship between inflammatory markers and low back pain. We obtained genetic data for CRP, along with its upstream inflammatory markers IL-6, IL-8, and IL-10, as well as low back pain from publicly available genome-wide association studies (GWAS). We applied several MR methods, including inverse variance weighting, weighted median, MR-Egger, Wald Ratio, and MR-PRESSO, to test for causal relationships. Sensitivity analyses were also conducted to assess the robustness of the results.

Results

Our analyses utilizing the Inverse Variance Weighted (IVW) method, the MR-Egger method, and the weighted median method indicated that IL-6 may be associated with an increased risk of LBP (Effect Size: -0.009, 95% Confidence Interval: -0.013-0.006, p = 9.16e-08); however, in the reverse direction, there was no significant causal effect of LBP on inflammatory markers.

Conclusion

Our study used a Mendelian randomization approach and found that elevated IL-6 levels may reduce the risk of LBP.

Causal Relationship Between Inflammation and Preeclampsia: Genetic Evidence from a Mendelian Randomization Study

Preeclampsia (PE) is a hypertensive disorder of pregnancy. PE patients were reported to have higher serum levels of C-reactive protein (CRP), interleukin-6 (IL-6) and tumor necrosis factor α (TNF-α) than those in healthy controls. However, whether the expressions of these inflammation biomarkers have a causal relationship with PE is unspecified. We applied the Mendelian randomization method to infer the causal relationship between inflammation biomarkers (e.g., CRP, IL-6, interleukin 1 receptor antagonist [IL-1ra] and TNF-α) and PE. Single nucleotide polymorphisms (SNPs) strongly related to inflammation biomarkers were used as instrumental variables. CRP, IL-1ra and IL-6 levels showed no significant effect on PE progression, while the genetic predicted higher level of TNF-α significantly increased the risk of PE (OR per 1-SD increase in TNF-α: 4.33; 95% CI [1.99, 9.39]; p = .00021). The findings suggest that pro-inflammatory activity of TNF-α could be a determinant for PE progression. More antenatal care should be given to those pregnant women with higher level of inflammation biomarkers, especially TNF-α.

Association between inflammation and cognition: Triangulation of evidence using a population-based cohort and Mendelian randomization analyses

BACKGROUND

Inflammation is associated with cognitive functioning and dementia in older adults, but whether inflammation is related to cognitive functioning in youth and whether these associations are causal remains unclear.

METHODS

In a population-based cohort (Avon Longitudinal Study of Parents and Children; ALSPAC), we investigated cross-sectional associations of inflammatory markers (C-reactive protein [CRP], Interleukin-6 [IL-6] and Glycoprotein acetyls [GlycA]) with measures of cold (working memory, response inhibition) and hot (emotion recognition) cognition at age 24 (N = 3,305 in multiple imputation models). Furthermore, we conducted one-sample and two-sample bidirectional Mendelian randomization (MR) analyses to examine potential causal effects of genetically-proxied inflammatory markers (CRP, GlycA, IL-6, IL-6 receptor, soluble IL-6 receptor) on cognitive measures (above) and on general cognitive ability.

RESULTS

In the ALSPAC cohort, there was limited evidence of an association between standardised inflammatory markers and standardised cognitive measures at age 24 after adjusting for potential confounders (N = 3,305; beta range, -0.02 [95 % confidence interval (CI) -0.06 to 0.02, p = 0.27] to 0.02 [95 % CI -0.02 to 0.05, p = 0.33]). Similarly, we found limited evidence of potential effects of 1-unit increase in genetically-proxied inflammatory markers on standardised working memory, emotion recognition or response inhibition in one-sample MR using ALSPAC data (beta range, -0.73 [95 % CI -2.47 to 1.01, p = 0.41] to 0.21 [95 % CI -1.42 to 1.84, p = 0.80]; or on standardised general cognitive ability in two-sample MR using the latest Genome-Wide Association Study (GWAS) datasets (inverse-variance weighted beta range, -0.02 [95 % CI -0.05 to 0.01, p = 0.12] to 0.03 [95 % CI -0.01 to 0.07, p = 0.19]).

CONCLUSIONS

Our MR findings do not provide strong evidence of a potential causal effect of inflammatory markers (CRP, IL-6, IL-6 receptor, GlycA) on the cognitive functions examined here. Given the large confidence intervals in the one-sample MR, larger GWAS of specific cognitive measures are needed to enable well-powered MR analyses to investigate whether inflammation causally influences specific cognitive domains.

Dissecting the causal association between inflammation and post-traumatic stress disorder: A bidirectional Mendelian randomization study

BACKGROUND

Accumulating evidence showed a bidirectional association between post-traumatic stress disorder and inflammation. However, whether the association is causal remains unclear. We aimed to evaluate the causal relationships between inflammatory cytokines and post-traumatic stress disorder using two-sample bi-directional Mendelian randomization analysis.

METHODS

Single nucleotide polymorphism from genome-wide association studies of inflammatory cytokines, C-reactive protein, and post-traumatic stress disorder (23,212 patients and 151,447 controls) was selected as instrumental variables. The causal associations were estimated by inverse variance weighting with sensitivity analyses using weighted median, MR-Egger, and MR-PRESSO methods.

RESULTS

We observed suggestive associations of genetically predicted interleukin-17 (IL-17) and RANTES with post-traumatic stress disorder. One standard deviation (SD) increase in genetically predicted level of IL-17 lowered the risk of post-traumatic stress disorder with an odds ratio (OR) of 0.902 (95 % CI = 0.828, 0.984, P = 0.02). One SD higher genetically predicted RANTES (CCL5) concentration increased post-traumatic stress disorder risk (OR = 1.067, 95 % CI = 1.005, 1.133, P = 0.032). However, we found no evidence of causal associations of post-traumatic stress disorder with the selected inflammatory cytokines and biomarkers. We observed no evidence supporting the presence of pleiotropy. The results of sensitivity analyses demonstrated the same directions and similar effect sizes as the primary findings.

LIMITATIONS

Potential pleiotropy, possible weak instruments, and low statistical power limited our findings.

CONCLUSION

Inflammation was suggestively causally associated with the risk of post-traumatic stress disorder, and inflammatory cytokines had no downstream effect on post-traumatic stress disorder. Further studies are needed to explain the mechanisms of systemic inflammation and neuroinflammation in post-traumatic stress disorder.

Cellular and molecular features of COVID-19 associated ARDS: therapeutic relevance

The severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) infection can be asymptomatic or cause a disease (COVID-19) characterized by different levels of severity. The main cause of severe COVID-19 and death is represented by acute (or acute on chronic) respiratory failure and acute respiratory distress syndrome (ARDS), often requiring hospital admission and ventilator support.The molecular pathogenesis of COVID-19-related ARDS (by now termed c-ARDS) is still poorly understood. In this review we will discuss the genetic susceptibility to COVID-19, the pathogenesis and the local and systemic biomarkers correlated with c-ARDS and the therapeutic options that target the cell signalling pathways of c-ARDS.

Gene expression profiles in mesenchymal stromal cells from bone marrow, adipose tissue and lung tissue of COPD patients and controls

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is characterized by irreversible lung tissue damage. Novel regenerative strategies are urgently awaited. Cultured mesenchymal stem/stromal cells (MSCs) have shown promising results in experimental models of COPD, but differences between sources may impact on their potential use in therapeutic strategies in patients.

AIM

To assess the transcriptome of lung-derived MSCs (LMSCs), bone marrow-derived MSCs (BM-MSC) and adipose-derived MSCs (AD-MSCs) from COPD patients and non-COPD controls.

METHODS

We studied differences in gene expression profiles between the MSC-subtypes, as well as between COPD and control using RNA sequencing (RNA-seq).

RESULTS

We show that besides heterogeneity between donors, MSCs from different sources have strongly divergent gene signatures. The growth factors FGF10 and HGF were predominantly expressed in LMSCs. MSCs from all sources displayed altered expression profiles in COPD, with most pronounced significantly up- and downregulated genes in MSCs from adipose tissue. Pathway analysis revealed that the most differentially expressed genes in COPD-derived AD-MSCs are involved in extracellular matrix (ECM) binding and expression. In LMSCs, the gene that differed most strongly between COPD and control was CSGALNACT1, an ECM modulating gene.

CONCLUSION

Autologous MSCs from COPD patients display abnormalities with respect to their transcriptome, which were surprisingly most profound in MSCs from extrapulmonary sources. LMSCs may be optimally equipped for lung tissue repair because of the expression of specific growth factor genes.

Mendelian randomization accounting for complex correlated horizontal pleiotropy while elucidating shared genetic etiology

Mendelian randomization (MR) harnesses genetic variants as instrumental variables (IVs) to study the causal effect of exposure on outcome using summary statistics from genome-wide association studies. Classic MR assumptions are violated when IVs are associated with unmeasured confounders, i.e., when correlated horizontal pleiotropy (CHP) arises. Such confounders could be a shared gene or inter-connected pathways underlying exposure and outcome. We propose MR-CUE (MR with Correlated horizontal pleiotropy Unraveling shared Etiology and confounding), for estimating causal effect while identifying IVs with CHP and accounting for estimation uncertainty. For those IVs, we map their cis-associated genes and enriched pathways to inform shared genetic etiology underlying exposure and outcome. We apply MR-CUE to study the effects of interleukin 6 on multiple traits/diseases and identify several S100 genes involved in shared genetic etiology. We assess the effects of multiple exposures on type 2 diabetes across European and East Asian populations.

IL6 gene polymorphism association with calcific aortic valve stenosis and influence on serum levels of interleukin-6

Aortic valve stenosis is the most frequent valve disease in developed countries and its prevalence will increase with population aging. There is still no pharmaceutical treatment nor biomarker to determine the susceptibility to develop aortic stenosis. Therefore, we analyzed the association of polymorphisms in risk loci with calcific aortic stenosis. Patients with aortic valve disease were genotyped for PALMD rs6702619, LPA rs10455872, and IL6 rs1800795 polymorphisms and circulating levels of interleukin-6 (IL-6) were measured. Calcium content of leaflets obtained in valve replacement surgeries was determined by micro-computed tomography. In the genotyping of 578 individuals, we found significant association between PALMD and IL6 polymorphisms and aortic stenosis in patients with tricuspid aortic valve, independently of other potentially confounding variables such as age and dyslipidemia. There was no association of these polymorphisms with valve calcium content, but this value correlated with the mean aortic pressure gradient (r = 0.44; P < 0.001). The CC genotype of IL6 polymorphism was associated with higher levels of serum IL-6 compared to other genotypes (23.5 vs. 10.5 pg/ml, respectively; P = 0.029). Therefore, patients carrying the CC genotype of IL6 rs1800795 polymorphism present higher levels of circulating IL-6 and this could contribute to the severity of the aortic valve stenosis. Our results agree with the identification of IL6 as a locus risk for stenosis and also with the intervention of this cytokine in aortic valve calcification. A more exhaustive follow-up of those patients carrying risk genotypes is therefore recommended.

A Mendelian randomization study investigating the causal role of inflammation on Parkinson’s disease

There is increasing evidence for inflammation as a determinant in the pathogenesis of Parkinson’s disease, but its role in parkinsonian neurodegeneration remains elusive. It is not clear whether inflammatory cascades are causes or consequences of dopamine neuron death. In the present study, we aim to perform an in-depth statistical investigation of the causal relationship between inflammation and Parkinson’s disease using a two-sample Mendelian randomization design.

Genetic instruments were selected using summary-level data from the largest genome-wide association studies to date (sample size ranging from 13 955 to 204 402 individuals) conducted on a European population for the following inflammation biomarkers: C-reactive protein, interleukin-6, interleukin 1 receptor antagonist and tumour necrosis factor α. Genetic association data on Parkinson’s disease (56 306 cases and 1 417 791 controls) and age at onset of Parkinson’s disease (28 568 cases) were obtained from the International Parkinson’s Disease Genomics Consortium. On primary analysis, causal associations were estimated on sets of strong (P-value < 5 × 10⁻⁸; F-statistic > 10) and independent (linkage disequilibrium r² < 0.001) genetic instruments using the inverse-variance weighted method. In sensitivity analysis, we estimated causal effects using robust Mendelian randomization methods and after removing pleiotropic genetic variants. Reverse causation was also explored. We repeated the analysis on different data sources for inflammatory biomarkers to check the consistency of the findings.

In all the three data sources selected for interleukin-6, we found statistical evidence for an earlier age at onset of Parkinson’s disease associated with increased interleukin-6 concentration [years difference per 1 log-unit increase = −2.364, 95% confidence interval (CI) = −4.789–0.060; years difference per 1 log-unit increase = −2.011, 95% CI = −3.706 to −0.317; years difference per 1 log-unit increase = −1.569, 95% CI = −2.891 to −0.247]. We did not observe any statistical evidence for causal effects of C-reactive protein, interleukin 1 receptor antagonist and tumour necrosis factor α on both Parkinson’s disease and its age at onset. Results after excluding possible pleiotropic genetic variants were consistent with findings from primary analyses. When investigating reverse causation, we did not find evidence for a causal effect of Parkinson’s disease or age at onset on any biomarkers of inflammation.

We found evidence for a causal association between the onset of Parkinson’s disease and interleukin-6. The findings of this study suggest that the pro-inflammatory activity of the interleukin-6 cytokine could be a determinant of prodromal Parkinson’s disease.

Body Mass Index, Interleukin-6 Signaling and Multiple Sclerosis: A Mendelian Randomization Study

Objectives

We explored whether genetically predicted increased body mass index (BMI) modulates multiple sclerosis (MS) risk through interleukin-6 (IL-6) signaling.

Methods

We performed a two-sample Mendelian randomization (MR) study using multiple genome-wide association studies (GWAS) datasets for BMI, IL-6 signaling, IL-6 levels and c-reactive protein (CRP) levels as exposures and estimated their effects on risk of MS from GWAS data from the International Multiple Sclerosis Genetics Consortium (IMSGC) in 14,802 MS cases and 26,703 controls.

Results

In univariable MR analyses, genetically predicted increased BMI and IL-6 signaling were associated with higher risk of MS (BMI: odds ratio (OR) = 1.30, 95% confidence interval (CI) = 1.15-1.47, p = 3.76 × 10-5; IL-6 signaling: OR = 1.51, 95% CI = 1.11-2.04, p = 0.01). Furthermore, higher BMI was associated with increased IL-6 signaling (β = 0.37, 95% CI = 0.32,0.41, p = 1.58 × 10-65). In multivariable MR analyses, the effect of IL-6 signaling on MS risk remained after adjusting for BMI (OR = 1.36, 95% CI = 1.11-1.68, p = 0.003) and higher BMI remained associated with an increased risk for MS after adjustment for IL-6 signaling (OR = 1.16, 95% CI =1.00-1.34, p = 0.046). The proportion of the effect of BMI on MS mediated by IL-6 signaling corresponded to 43% (95% CI = 25%-54%). In contrast to IL-6 signaling, there was little evidence for an effect of serum IL-6 levels or CRP levels on risk of MS.

Conclusion

In this study, we identified IL-6 signaling as a major mediator of the association between BMI and risk of MS. Further explorations of pathways underlying the association between BMI and MS are required and will, together with our findings, improve the understanding of MS biology and potentially lead to improved opportunities for targeted prevention strategies.

Interleukin 6 in diabetes, chronic kidney disease and cardiovascular disease: mechanisms and therapeutic perspectives

INTRODUCTION

Diabetes, chronic kidney disease (CKD) and cardiovascular disease (CVD) are cardiometabolic diseases that remain amongst the leading causes of morbidity and premature mortality. Here, we review the current understanding of how anti-inflammatory intervention via inhibition of the pro-inflammatory but pleiotropic cytokine interleukin (IL) 6 may benefit patients with these or related diseases or complications.

AREAS COVERED

Based on a PubMed literature search, this review integrates and contextualizes evidence regarding the clinical utility of anti-IL-6 intervention in the treatment of cardiometabolic diseases, as well as of the associated condition non-alcoholic hepatosteatosis.

EXPERT OPINION

Evidence implicates the pro-inflammatory effects of IL-6 in the pathophysiology of diabetes, CKD and CVD. Thus, targeting the IL-6 pathway holds a therapeutic potential in these cardiometabolic disorders. However, because IL-6 has multiple homeostatic roles, antagonizing this cytokine may be associated with side effects such as increased risk of infection as seen with other anti-inflammatory drugs. Additional studies are required to establish the benefit-risk profile of anti-IL-6 intervention in the cardiometabolic diseases, whilst also considering alternative interventions such as lifestyle changes. IL-6 is also elevated in NASH, but the clinical usefulness of targeting IL-6 in this hepatic disorder remains largely unexplored.

Quantile-Dependent Expressivity of Serum Interleukin-6 Concentrations as a Possible Explanation of Gene-Disease Interactions, Gene-Environment Interactions, and Pharmacogenetic Effects

Interleukin 6 (IL-6) is a moderately heritable pleiotropic cytokine whose elevated concentrations in coronary artery disease, peripheral arterial disease, pulmonary arterial hypertension, Eales' disease, Sjògren's syndrome, osteoarthritis, adenocarcinoma, neuroblastoma, polymyalgia rheumatica, pulmonary tuberculosis, and enterovirus 71 infection, and following coronary artery bypass graft show larger genetic effects than in unaffected low IL-6 controls. We hypothesize that genetic effects may depend upon whether average IL-6 concentrations are high or low, i.e., quantile-dependent expressivity. Quantile-specific offspring-parent (β_OP) and full-sib regression slopes (β_FS) were estimated by applying quantile regression to the age- and sex-adjusted serum IL-6 concentrations in families surveyed in the Framingham Heart Study. Quantile-specific heritabilities were calculated as h² = 2β_OP / (1 + r_spouse) and h² = {(1 + 8r_spouseβ_FS)^0.5 -1} / (2r_spouse)). Heritability (h² ± SE) of IL-6 concentrations increased from 0.01 ± 0.01 at the 10th percentile (NS), 0.02 ± 0.01 at the 25th (P = 0.009), 0.03 ± 0.01 at the 50th (P = 0.007), 0.04 ± 0.02 at the 75th (P = 0.004), and 0.13 ± 0.05 at the 90th percentile (P = 0.03), or 0.0005 ± 0.0002 for each 1% increase in the offspring's phenotype distribution (P_{linear trend} = 0.02) when estimated from β_OP and from 0.02 ± 0.02 at the 10th (NS), 0.02 ± 0.02 at the 25th (NS), 0.06 ± 0.02 at the 50th (P = 0.01), 0.12 ± 0.04 at the 75th (P = 0.001), and 0.30 ± 0.03 at the 90th percentile (P < 10^-16), or 0.0015 ± 0.0007 for each 1% increase in the sibling phenotype distribution (P_{linear trend} = 0.02) when estimated from β_FS. Thus the heritability of serum IL-6 concentrations is quantile dependent, which may contribute in part to the larger genetic effect size reported for diseases and environmental conditions that elevate IL-6 concentrations vis-à-vis unaffected controls.

Dissecting the IL-6 pathway in cardiometabolic disease: a Mendelian randomization study on both IL6 and IL6R

BACKGROUND

Chronic inflammation is a risk factor for cardiovascular disease. IL-6 signaling perturbation through IL-6 or IL-6R blockade may have potential benefit on cardiovascular risk. It is unknown whether targeting either IL-6 or IL-6 receptor may result in similar effects on CVD and adverse events. We compared the anticipated effects of targeting IL-6 and IL-6 receptor on cardiometabolic risk and potential side effects.

METHODS

We constructed four instruments: two main instruments with genetic variants in the IL6 and IL6R loci weighted for their association with CRP, and two after firstly filtering variants for their association with IL-6 or IL-6R expression. Analyses were performed for coronary artery disease (CAD), ischemic stroke, atrial fibrillation (AF), heart failure, type 2 diabetes (T2D), rheumatoid arthritis (RA), infection endpoints, and quantitative hematological, metabolic, and anthropometric parameters.

RESULTS

A 1 mg/L lower CRP by the IL6 instrument was associated with lower CAD (OR 0.86, 95% CI 0.77;0.96), AF, and T2D risk. A 1mg/L lower CRP by the IL6R instrument was associated with lower CAD (OR 0.90, 95% CI 0.86;0.95), any stroke and ischemic stroke, AF, RA risk and higher pneumonia risk. The eQTL filtered results were in concordance with the main results, but with wider confidence intervals.

CONCLUSIONS

IL-6 signalling perturbation by either IL6 or IL6R genetic instruments is associated with a similar risk reduction for multiple cardiometabolic diseases, suggesting that both IL-6 and IL-6R are potential therapeutic targets to lower CVD. Moreover, IL-6 rather than IL-6R inhibition might have a more favorable pneumonia risk.

The Association of Human Leukocyte Antigen and COVID-19 in Southern China

Human leukocyte antigen (HLA) polymorphism is hypothesized to be associated with diverse immune responses toward infectious diseases. Herein, by comparing against multiple subpopulation groups as control, we confirmed that HLA-B*15:27 and HLA-DRB1*04:06 were associated with coronavirus disease 2019 susceptibility in China. Both alleles were predicted to have weak binding affinities toward viral proteins.