Genome-wide association study identifies seven novel loci associating with circulating cytokines and cell adhesion molecules in Finns

Genetic Associations with Placental Proteins in Maternal Serum Identify Biomarkers for Hypertension in Pregnancy

Background

Preeclampsia is a complex syndrome that accounts for considerable maternal and perinatal morbidity and mortality. Despite its prevalence, no effective disease-modifying therapies are available. Maternal serum placenta-derived proteins have been in longstanding use as markers of risk for aneuploidy and placental dysfunction, but whether they have a causal contribution to preeclampsia is unknown.

Objective

We aimed to investigate the genetic regulation of serum placental proteins in early pregnancy and their potential causal links with preeclampsia and gestational hypertension.

Study design

This study used a nested case-control design with nulliparous women enrolled in the nuMoM2b study from eight clinical sites across the United States between 2010 and 2013. The first- and second-trimester serum samples were collected, and nine proteins were measured, including vascular endothelial growth factor (VEGF), placental growth factor, endoglin, soluble fms-like tyrosine kinase-1 (sFlt-1), a disintegrin and metalloproteinase domain-containing protein 12 (ADAM-12), pregnancy-associated plasma protein A, free beta-human chorionic gonadotropin, inhibin A, and alpha-fetoprotein. This study used genome-wide association studies to discern genetic influences on these protein levels, treating proteins as outcomes. Furthermore, Mendelian randomization was used to evaluate the causal effects of these proteins on preeclampsia and gestational hypertension, and their further causal relationship with long-term hypertension, treating proteins as exposures.

Results

A total of 2,352 participants were analyzed. We discovered significant associations between the pregnancy zone protein locus and concentrations of ADAM-12 (rs6487735, P= 3.03×10 -22 ), as well as between the vascular endothelial growth factor A locus and concentrations of both VEGF (rs6921438, P= 7.94×10 -30 ) and sFlt-1 (rs4349809, P= 2.89×10 -12 ). Our Mendelian randomization analyses suggested a potential causal association between first-trimester ADAM-12 levels and gestational hypertension (odds ratio=0.78, P= 8.6×10 -4 ). We also found evidence for a potential causal effect of preeclampsia (odds ratio=1.75, P =8.3×10 -3 ) and gestational hypertension (odds ratio=1.84, P =4.7×10 -3 ) during the index pregnancy on the onset of hypertension 2-7 years later. The additional mediation analysis indicated that the impact of ADAM-12 on postpartum hypertension could be explained in part by its indirect effect through gestational hypertension (mediated effect=-0.15, P= 0.03).

Conclusions

Our study discovered significant genetic associations with placental proteins ADAM-12, VEGF, and sFlt-1, offering insights into their regulation during pregnancy. Mendelian randomization analyses demonstrated evidence of potential causal relationships between the serum levels of placental proteins, particularly ADAM-12, and gestational hypertension, potentially informing future prevention and treatment investigations.

Common variation in a long non-coding RNA gene modulates variation of circulating TGF-β2 levels in metastatic colorectal cancer patients (Alliance)

BACKGROUND

Herein, we report results from a genome-wide study conducted to identify protein quantitative trait loci (pQTL) for circulating angiogenic and inflammatory protein markers in patients with metastatic colorectal cancer (mCRC). The study was conducted using genotype, protein marker, and baseline clinical and demographic data from CALGB/SWOG 80405 (Alliance), a randomized phase III study designed to assess outcomes of adding VEGF or EGFR inhibitors to systemic chemotherapy in mCRC patients. Germline DNA derived from blood was genotyped on whole-genome array platforms. The abundance of protein markers was quantified using a multiplex enzyme-linked immunosorbent assay from plasma derived from peripheral venous blood collected at baseline. A robust rank-based method was used to assess the statistical significance of each variant and protein pair against a strict genome-wide level. A given pQTL was tested for validation in two external datasets of prostate (CALGB 90401) and pancreatic cancer (CALGB 80303) patients. Bioinformatics analyses were conducted to further establish biological bases for these findings.

RESULTS

The final analysis was carried out based on data from 540,021 common typed genetic variants and 23 protein markers from 869 genetically estimated European patients with mCRC. Correcting for multiple testing, the analysis discovered a novel cis-pQTL in LINC02869, a long non-coding RNA gene, for circulating TGF-β2 levels (rs11118119; AAF = 0.11; P-value < 1.4e-14). This finding was validated in a cohort of 538 prostate cancer patients from CALGB 90401 (AAF = 0.10, P-value < 3.3e-25). The analysis also validated a cis-pQTL we had previously reported for VEGF-A in advanced pancreatic cancer, and additionally identified trans-pQTLs for VEGF-R3, and cis-pQTLs for CD73.

CONCLUSIONS

This study has provided evidence of a novel cis germline genetic variant that regulates circulating TGF-β2 levels in plasma of patients with advanced mCRC and prostate cancer. Moreover, the validation of previously identified pQTLs for VEGF-A, CD73, and VEGF-R3, potentiates the validity of these associations.

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.

Genomic analysis of venous thrombosis in children with acute lymphoblastic leukemia from diverse ancestries

Venous thrombosis is a common adverse effect of modern therapy for acute lymphoblastic leukemia (ALL). Prior studies to identify risks of thrombosis in pediatric ALL have been limited by genetic screens of pre-identified genetic variants or genome- wide association studies (GWAS) in ancestrally uniform populations. To address this, we performed a retrospective cohort evaluation of thrombosis risk in 1,005 children treated for newly diagnosed ALL. Genetic risk factors were comprehensively evaluated from genome-wide single nucleotide polymorphism (SNP) arrays and were evaluated using Cox regression adjusting for identified clinical risk factors and genetic ancestry. The cumulative incidence of thrombosis was 7.8%. In multivariate analysis, older age, T-lineage ALL, and non-O blood group were associated with increased thrombosis while non-low-risk treatment and higher presenting white blood cell count trended toward increased thrombosis. No SNP reached genome-wide significance. The SNP most strongly associated with thrombosis was rs2874964 near RFXAP (G risk allele; P=4x10-7; hazard ratio [HR] =2.8). In patients of non-European ancestry, rs55689276 near the α globin cluster (P=1.28x10-6; HR=27) was most strongly associated with thrombosis. Among GWAS catalogue SNP reported to be associated with thrombosis, rs2519093 (T risk allele, P=4.8x10-4; HR=2.1), an intronic variant in ABO, was most strongly associated with risk in this cohort. Classic thrombophilia risks were not associated with thrombosis. Our study confirms known clinical risk features associated with thrombosis risk in children with ALL. In this ancestrally diverse cohort, genetic risks linked to thrombosis risk aggregated in erythrocyte-related SNP, suggesting the critical role of this tissue in thrombosis risk.

Common variation in a long non-coding RNA gene modulates variation of circulating TGF-β2 levels in metastatic colorectal cancer patients (Alliance)

Background

Herein, we report results from a genome-wide study conducted to identify protein quantitative trait loci (pQTL) for circulating angiogenic and inflammatory protein markers in patients with metastatic colorectal cancer (mCRC).The study was conducted using genotype, protein marker, and baseline clinical and demographic data from CALGB/SWOG 80405 (Alliance), a randomized phase III study designed to assess outcomes of adding VEGF or EGFR inhibitors to systemic chemotherapy in mCRC patients. Germline DNA derived from blood was genotyped on whole-genome array platforms. The abundance of protein markers was quantified using a multiplex enzyme-linked immunosorbent assay from plasma derived from peripheral venous blood collected at baseline. A robust rank-based method was used to assess the statistical significance of each variant and protein pair against a strict genome-wide level. A given pQTL was tested for validation in two external datasets of prostate (CALGB 90401) and pancreatic cancer (CALGB 80303) patients. Bioinformatics analyses were conducted to further establish biological bases for these findings.

Results

The final analysis was carried out based on data from 540,021 common typed genetic variants and 23 protein markers from 869 genetically estimated European patients with mCRC. Correcting for multiple testing, the analysis discovered a novel cis-pQTL in LINC02869, a long non-coding RNA gene, for circulating TGF-β2 levels (rs11118119; AAF = 0.11; P-value < 1.4e-14). This finding was validated in a cohort of 538 prostate cancer patients from CALGB 90401 (AAF = 0.10, P-value < 3.3e-25). The analysis also validated a cis-pQTL we had previously reported for VEGF-A in advanced pancreatic cancer, and additionally identified trans-pQTLs for VEGF-R3, and cis-pQTLs for CD73.

Conclusions

This study has provided evidence of a novel cis germline genetic variant that regulates circulating TGF-β2 levels in plasma of patients with advanced mCRC and prostate cancer. Moreover, the validation of previously identified pQTLs for VEGF-A, CD73, and VEGF-R3, potentiates the validity of these associations.

Circulating Inflammatory Cytokines and Female Reproductive Diseases: A Mendelian Randomization Analysis

CONTEXT

Extensive studies have provided considerable evidence suggesting the role of inflammation in the development of female reproductive diseases. However, causality has not been established.

OBJECTIVE

To explore whether genetically determined circulating levels of cytokines are causally associated with female reproductive diseases and discover potential novel drug targets for these diseases.

METHODS

Instrumental variables (IVs) for 47 circulating cytokines were obtained from a genome-wide association study (GWAS) meta-analysis of 31 112 European individuals. Protein quantitative trait loci and expression quantitative trait loci close to genes served as our IVs. Summary data of 9 female reproductive diseases were mainly derived from GWAS meta-analysis of the UK biobank and FinnGen. We elevated the association using the Wald ratio or inverse variance-weighted Mendelian randomization (MR) with subsequent assessments for MR assumptions in several sensitivity and colocalization analyses. We consider a false discovery rate <0.05 as statistical significance in MR analyses. Replication studies were conducted for further validation, and phenome-wide association studies were designed to explore potential side effects.

RESULTS

Our results indicated that high levels of macrophage colony-stimulating factor (MCSF), growth-regulated oncogene-alpha (GROα), and soluble intercellular adhesion molecule-1 were associated with increased risks of endometriosis, female infertility, and pre-eclampsia, respectively. High platelet-derived growth factor-BB (PDGF-BB) levels that reduced the risk of ovarian aging were also supported. Replication analysis supported the relationship between GROα and female infertility, and between MCSF and endometriosis.

CONCLUSION

We identified 4 correlated pairs that implied potential protein drug targets. Notably, we preferred highlighting the value of PDGF-BB as a drug target for ovarian aging, and MCSF as a drug target for endometriosis.

The influence of gut microbiota on circulating inflammatory cytokines and host: A Mendelian randomization study with meta-analysis

AIMS

The gut microbiota has been found to be altered in different inflammatory disorders, but its involvement in the regulation of inflammatory cytokines remains unclear. Therefore, this study aimed to investigate the impacts of gut microbiota on circulating inflammatory cytokines and their potential roles in host diseases.

MAIN METHODS

Two-sample Mendelian randomization (MR) analyses were conducted using summary-level data from genome-wide association studies (GWAS) to identify significant causal associations between 196 gut microbiota and 41 inflammatory cytokines. Meta-analysis was applied to test the robustness of the results. Enrichment analyses of identified cytokines were further utilized to infer the effects of gut microbiota on the host.

KEY FINDINGS

The MR analyses and meta-analyses identified the following significant causal associations: phylum Euryarchaeota on interleukin-2 (IL-2) (βIVW = 0.085, P = 1.5 × 10-2) and interleukin-8 (IL-8) (βIVW = 0.065, P = 4.1 × 10-2), phylum Tenericutes and class Mollicutes on macrophage inflammatory protein 1a (MIP1a) (βIVW = -0.142, P = 7.0 × 10-3), class Bacilli on hepatocyte growth factor (HGF) (βIVW = -0.106, P = 2.5 × 10-2), order Enterobacteriales on monocyte chemoattractant protein-1 (MCP1) (βIVW = 0.182, P = 1.8 × 10-2), and genus Lachnospiraceae NC2004 group on TNF-related apoptosis-inducing ligand (TRAIL) (βIVW = -0.207, P = 6.0 × 10-4). Enrichment analyses suggested that phylum Euryarchaeota and order Enterobacteriales might be risk factors for certain autoimmune diseases and neoplasms, while the phylum Tenericutes may have a protective effect.

SIGNIFICANCE

This study represents the first evidence confirming the causal effect of specific gut microbial taxa on circulating inflammatory cytokines and sheds light on their potential roles in the development and progression of various host diseases.

Circulating inflammatory cytokines in relation to the risk of renal cell carcinoma: A gender-specific two-sample Mendelian randomization study

BACKGROUND

Currently there is no specific molecular biomarker for the diagnosis and treatment of renal cell carcinoma (RCC). Here we performed a gender-specific two-sample Mendelian randomization analysis to systematically assess the effects of circulating cytokines on RCC.

METHODS

We have employed cis-quantitative trait loci as instrumental variables for the protein levels and expression of circulating cytokines. We estimated the causal effects of circulating cytokines on RCC risk in males and females with several Mendelian randomization methods.

RESULTS

We observed a significant causal effect of Eotaxin on the increased risk of RCC in males (Odds ratio [OR] = 2.546, 95% confidence interval [CI] = 1.617-4.010, p value = 5.496 × 10-5), but not in females (OR = 1.352, 95% CI = 0.766-2.388, p value = 0.298). Besides, we also identified several cytokines as potentially associated with RCC in males including RANTES, MCP3, PDGFbb, TRAIL, and several other cytokines as potentially associated with RCC in females including sICAM and SCGFb.

CONCLUSION

Our study highlighted that a higher level of circulating Eotaxin is causally associated with an increased risk of RCC in males but not in females. Further studies are needed to elucidate the exact mechanism and its potential application in the prognosis and treatment of RCC.

Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease characterized by insulin deficiency and loss of pancreatic islet β-cells. The objective of this study is to identify de novo mutations in 13 trios from singleton families that contribute to the genetic basis of T1DM through the application of whole-exome sequencing (WES). Of the 13 families sampled for this project, 12 had de novo variants, with Family 7 having the highest number (nine) of variants linked to T1DM/autoimmune pathways, whilst Family 4 did not have any variants past the filtering steps. There were 10 variants of 7 genes reportedly associated with T1DM (MST1; TDG; TYRO3; IFIHI; GLIS3; VEGFA; TYK2). There were 20 variants of 13 genes that were linked to endocrine, metabolic, or autoimmune diseases. Our findings demonstrate that trio-based WES is a powerful approach for identifying new candidate genes for the pathogenesis of T1D. Genotyping and functional annotation of the discovered de novo variants in a large cohort is recommended to ascertain their association with disease pathogenesis.

The interplay between inflammatory cytokines and cardiometabolic disease: bi-directional mendelian randomisation study

Objective

To leverage large scale genetic association data to investigate the interplay between circulating cytokines and cardiometabolic traits, and thus identifying potential therapeutic targets.

Design

Bi-directional Mendelian randomisation study.

Setting

Genome-wide association studies from three Finnish cohorts (Northern Finland Birth Cohort 1966, Young Finns Study, or FINRISK study), and genetic association summary statistics pooled from observational studies for expression quantitative trait loci and cardiometabolic traits.

Participants

Data for 47 circulating cytokines in 13 365 individuals from genome-wide association studies, summary statistic data for up to 21 735 individuals on circulating cytokines, summary statistic gene expression data across 49 tissues in 838 individuals, and summary statistic data for up to 1 320 016 individuals on cardiometabolic traits.

Interventions

Relations between circulating cytokines and cardiovascular, anthropometric, lipid, or glycaemic traits (coronary artery disease, stroke, type 2 diabetes mellitus, body mass index, waist circumference, waist to hip ratio, systolic blood pressure, glycated haemoglobin, high density lipoprotein cholesterol, low density lipoprotein cholesterol, total cholesterol, triglycerides, C reactive protein, glucose, fasting insulin, and lifetime smoking).

Main outcome methods

Genetic instrumental variables that are biologically plausible for the circulating cytokines were generated. The effects of cardiometabolic risk factors on concentrations of circulating cytokines, circulating cytokines on other circulating cytokines, and circulating cytokines on cardiometabolic outcomes were investigated.

Results

Genetic evidence (mendelian randomisation P<0.0011) suggests that higher body mass index, waist circumference, smoking, higher concentrations of lipids, and systolic blood pressure increase circulating concentrations of several inflammatory cytokines and C reactive protein. Evidence for causal relations (mendelian randomisation P<0.0011) were noted between circulating cytokines, including a key role of vascular endothelial growth factor on influencing the concentrations of 10 other cytokines. Both mendelian randomisation (P<0.05) and colocalisation (posterior probability >0.5) suggested that coronary artery disease risk is increased by higher concentrations of circulating tumour necrosis factor related apoptosis-inducing ligand (TRAIL), interleukin-1 receptor antagonist (IL1RA), and macrophage colony-stimulating factor (MCSF).

Conclusion

This study offers insight into inflammatory mediators of cardiometabolic risk factors, cytokine signalling cascades, and effects of circulating cytokines on different cardiometabolic outcomes.

The role of VCAM-1 in diabetic retinopathy: A systematic review and meta-analysis

OBJECTIVE

Vascular cell adhesion molecule-1 (VCAM-1) plays a regulatory role in inflammatory diseases. However, the exact role of VCAM-1 in diabetic retinopathy (DR) remains unclear, and there is a lack of meta-analyses.

METHODS

The role of VCAM-1 in DR was screened by database searching. A random effects model was used, and the estimated mean difference was evaluated.

RESULTS

Twenty articles were included. The level of VCAM-1 increased significantly in the DR group compared with the control group (SMD: 0.67, 95 % CI: 0.34-1.01, P < 0.0001). VCAM-1 levels correlated with sample size and DR type, method and severity based on subgroup analysis.

CONCLUSION

A high level of VCAM-1 is present in DR patients and is related to the severity of DR. Therefore, VCAM-1 is a potential detection biomarker for DR.

Ten interleukins and risk of prostate cancer

Background

Interleukins (ILs) have been reported to be related to prostate cancer. The aims of this study were to estimate the levels for several key interleukins in prostate cancer and the causal effects between them.

Methods

We conducted a bi-directional two-sample Mendelian randomization (MR) study to assess the causal associations between ILs and prostate cancer. Genetic instruments and summary-level data for 10 ILs were obtained from three genome-wide association meta-analyses. Prostate cancer related data were obtained from the PRACTICAL (79,148 cases and 61,106 controls), UK Biobank (7,691 cases and 169,762 controls) and FinnGen consortium (10,414 cases and 124,994 controls), respectively.

Results

The odds ratio of prostate cancer was 0.92 (95% confidence interval (CI), 0.89, 0.96; P=1.58×10^-05) and 1.12 (95% CI, 1.07, 1.17; P=6.61×10^-07) for one standard deviation increase in genetically predicted IL-1ra and IL-6 levels, respectively. Genetically predicted levels of IL-1ß, IL-2a, IL-6ra, IL-8, IL-16, IL-17, IL-18, and IL-27 were not associated with the risk of prostate cancer. Reverse MR analysis did not find the associations between genetic liability to prostate cancer and higher levels of IL-1ra (β, -0.005; 95% CI, -0.010, 0.001; P=0.111) and IL-6 (β, 0.002; 95% CI, -0.011, 0.014; P=0.755).

Conclusion

This MR study suggests that long-term IL-6 may increase the risk of prostate cancer and IL-1ra may reduce it.

Germline genetic variation and predicting immune checkpoint inhibitor induced toxicity

Immune checkpoint inhibitor (ICI) therapy has revolutionised the treatment of various cancer types. ICIs reinstate T-cell function to elicit an anti-cancer immune response. The resulting immune response can however have off-target effects which manifest as autoimmune type serious immune-related adverse events (irAE) in ~10-55% of patients treated. It is currently challenging to predict both who will experience irAEs and to what severity. Identification of patients at high risk of serious irAE would revolutionise patient care. While the pathogenesis driving irAE development is still unclear, host genetic factors are proposed to be key determinants of these events. This review presents current evidence supporting the role of the host genome in determining risk of irAE. We summarise the spectrum and timing of irAEs following treatment with ICIs and describe currently reported germline genetic variation associated with expression of immuno-modulatory factors within the cancer immunity cycle, development of autoimmune disease and irAE occurrence. We propose that germline genetic determinants of host immune function and autoimmune diseases could also explain risk of irAE development. We also endorse genome-wide association studies of patients being treated with ICIs to identify genetic variants that can be used in polygenic risk scores to predict risk of irAE.

Genetic influences of the effect of circulating inflammatory cytokines on osteoarthritis in humans

OBJECTIVE

The causal relationship between inflammatory cytokines and Osteoarthritis (OA) has not been well investigated. This study investigated the causal role of inflammatory cytokines in the risk of OA and total joint arthroplasty using the Mendelian randomization (MR) method.

METHOD

Single nucleotide polymorphisms (SNPs) robustly associated with inflammatory cytokines were used as instrumental variables. The inverse-variance weighted (IVW) method with false discovery rate (FDR) adjusted P-value (q-value) for multiple comparisons were used as the main MR method to estimate causal effects based on the summary-level data for OA (knee and hip OA, respectively) and total joint arthroplasty (TJA). Sensitivity analyses validated the robustness of the results and ensured the absence of heterogeneity and horizontal pleiotropy.

RESULTS

After FDR adjustment, macrophage colony-stimulating factor (MCSF) and vascular endothelial growth factor (VEGF) were identified as causally associated with knee OA (MCSF, odds ratio [OR]: 1.16, 95% confidence interval [CI]: 1.09-1.23, q = 5.05 × 10^-5; VEGF, OR: 1.09, 95% CI: 1.04-1.15, q = 0.011). We also observed that genetically predicted MCSF and VEGF were positively associated with the risk of TJA, and MCP3 was negatively associated with for the risk of TJA, although the effects seem fairly modest. Sensitivity analysis further excluded the influence of heterogeneity and horizontal pleiotropy.

CONCLUSIONS

Inflammatory cytokines, namely MCSF and VEGF, were causally associated with knee OA, which could enhance our understanding of inflammation in OA pathology.

Connecting aging biology and inflammation in the omics era

Aging is characterized by the accumulation of damage to macromolecules and cell architecture that triggers a proinflammatory state in blood and solid tissues, termed inflammaging. Inflammaging has been implicated in the pathogenesis of many age-associated chronic diseases as well as loss of physical and cognitive function. The search for mechanisms that underlie inflammaging focused initially on the hallmarks of aging, but it is rapidly expanding in multiple directions. Here, we discuss the threads connecting cellular senescence and mitochondrial dysfunction to impaired mitophagy and DNA damage, which may act as a hub for inflammaging. We explore the emerging multi-omics efforts that aspire to define the complexity of inflammaging - and identify molecular signatures and novel targets for interventions aimed at counteracting excessive inflammation and its deleterious consequences while preserving the physiological immune response. Finally, we review the emerging evidence that inflammation is involved in brain aging and neurodegenerative diseases. Our goal is to broaden the research agenda for inflammaging with an eye on new therapeutic opportunities.

Proteome-wide Mendelian randomization identifies causal links between blood proteins and severe COVID-19

In November 2021, the COVID-19 pandemic death toll surpassed five million individuals. We applied Mendelian randomization including >3,000 blood proteins as exposures to identify potential biomarkers that may indicate risk for hospitalization or need for respiratory support or death due to COVID-19, respectively. After multiple testing correction, using genetic instruments and under the assumptions of Mendelian Randomization, our results were consistent with higher blood levels of five proteins GCNT4, CD207, RAB14, C1GALT1C1, and ABO being causally associated with an increased risk of hospitalization or respiratory support/death due to COVID-19 (ORs = 1.12-1.35). Higher levels of FAAH2 were solely associated with an increased risk of hospitalization (OR = 1.19). On the contrary, higher levels of SELL, SELE, and PECAM-1 decrease risk of hospitalization or need for respiratory support/death (ORs = 0.80-0.91). Higher levels of LCTL, SFTPD, KEL, and ATP2A3 were solely associated with a decreased risk of hospitalization (ORs = 0.86-0.93), whilst higher levels of ICAM-1 were solely associated with a decreased risk of respiratory support/death of COVID-19 (OR = 0.84). Our findings implicate blood group markers and binding proteins in both hospitalization and need for respiratory support/death. They, additionally, suggest that higher levels of endocannabinoid enzymes may increase the risk of hospitalization. Our research replicates findings of blood markers previously associated with COVID-19 and prioritises additional blood markers for risk prediction of severe forms of COVID-19. Furthermore, we pinpoint druggable targets potentially implicated in disease pathology.

A plasma metabolite score of three eicosanoids predicts incident type 2 diabetes: a prospective study in three independent cohorts

INTRODUCTION

Peptide markers of inflammation have been associated with the development of type 2 diabetes. The role of upstream, lipid-derived mediators of inflammation such as eicosanoids, remains less clear. The aim of this study was to examine whether eicosanoids are associated with incident type 2 diabetes.

RESEARCH DESIGN & METHODS

In the FINRISK (Finnish Cardiovascular Risk Study) 2002 study, a population-based sample of Finnish men and women aged 25-74 years, we used directed, non-targeted liquid chromatography-mass spectrometry to identify 545 eicosanoids and related oxylipins in the participants' plasma samples (n=8292). We used multivariable-adjusted Cox regression to examine associations between eicosanoids and incident type 2 diabetes. The significant independent findings were replicated in the Framingham Heart Study (FHS, n=2886) and DIetary, Lifestyle and Genetic determinants of Obesity and Metabolic syndrome (DILGOM) 2007 (n=3905). Together, these three cohorts had 1070 cases of incident type 2 diabetes.

RESULTS

In the FINRISK 2002 cohort, 76 eicosanoids were associated individually with incident type 2 diabetes. We identified three eicosanoids independently associated with incident type 2 diabetes using stepwise Cox regression with forward selection and a Bonferroni-corrected inclusion threshold. A three-eicosanoid risk score produced an HR of 1.56 (95% CI 1.41 to 1.72) per 1 SD increment for risk of incident diabetes. The HR for comparing the top quartile with the lowest was 2.80 (95% CI 2.53 to 3.07). In the replication analyses, the three-eicosanoid risk score was significant in FHS (HR 1.24 (95% CI 1.10 to 1.39, p<0.001)) and directionally consistent in DILGOM (HR 1.12 (95% CI 0.99 to 1.27, p=0.07)). Meta-analysis of the three cohorts yielded a pooled HR of 1.31 (95% CI 1.05 to 1.56).

CONCLUSIONS

Plasma eicosanoid profiles predict incident type 2 diabetes and the clearest signals replicate in three independent cohorts. Our findings give new information on the biology underlying type 2 diabetes and suggest opportunities for early identification of people at risk.

Circulating inflammatory cytokines and risk of five cancers: a Mendelian randomization analysis

BACKGROUND

Epidemiological and experimental evidence has linked chronic inflammation to cancer aetiology. It is unclear whether associations for specific inflammatory biomarkers are causal or due to bias. In order to examine whether altered genetically predicted concentration of circulating cytokines are associated with cancer development, we performed a two-sample Mendelian randomisation (MR) analysis.

METHODS

Up to 31,112 individuals of European descent were included in genome-wide association study (GWAS) meta-analyses of 47 circulating cytokines. Single nucleotide polymorphisms (SNPs) robustly associated with the cytokines, located in or close to their coding gene (cis), were used as instrumental variables. Inverse-variance weighted MR was used as the primary analysis, and the MR assumptions were evaluated in sensitivity and colocalization analyses and a false discovery rate (FDR) correction for multiple comparisons was applied. Corresponding germline GWAS summary data for five cancer outcomes (breast, endometrial, lung, ovarian, and prostate), and their subtypes were selected from the largest cancer-specific GWASs available (cases ranging from 12,906 for endometrial to 133,384 for breast cancer).

RESULTS

There was evidence of inverse associations of macrophage migration inhibitory factor with breast cancer (OR per SD = 0.88, 95% CI 0.83 to 0.94), interleukin-1 receptor antagonist with endometrial cancer (0.86, 0.80 to 0.93), interleukin-18 with lung cancer (0.87, 0.81 to 0.93), and beta-chemokine-RANTES with ovarian cancer (0.70, 0.57 to 0.85) and positive associations of monokine induced by gamma interferon with endometrial cancer (3.73, 1.86 to 7.47) and cutaneous T-cell attracting chemokine with lung cancer (1.51, 1.22 to 1.87). These associations were similar in sensitivity analyses and supported in colocalization analyses.

CONCLUSIONS

Our study adds to current knowledge on the role of specific inflammatory biomarker pathways in cancer aetiology. Further validation is needed to assess the potential of these cytokines as pharmacological or lifestyle targets for cancer prevention.

Genetically Determined Inflammatory Biomarkers and the Risk of Heart Failure: A Mendelian Randomization Study

Background: Positive associations between inflammatory biomarkers and the risk of heart failure (HF) have been reported in conventional observational studies. However, the causal effects of inflammatory biomarkers on HF have not been fully elucidated. We conducted a Mendelian randomization (MR) study to examine the possible etiological roles of inflammatory biomarkers in HF. Methods: Summary statistical data for the associations between single nucleotide polymorphisms (SNPs) and C-reactive protein (CRP), fibrinogen, and components of the interleukin-1 (IL-1)-interleukin-6 (IL-6) inflammatory signaling pathway, namely, interleukin-1β (IL-1β), IL-1 receptor antagonist (IL-1ra), IL-6, and soluble IL-6 receptor (sIL-6r), were obtained from genome-wide association studies (GWASs) for individuals of European descent. The GWAS dataset of 977,323 participants of European ancestry, which included 47,309 HF cases and 930,014 controls, was collected to identify genetic variants underlying HF. A two-sample Mendelian randomization framework was implemented to examine the causality of the association between these inflammatory biomarkers and HF. Results: Our MR analyses found that genetically determined CRP and fibrinogen were not causally associated with HF risk (odds ratio [OR] = 0.93, 95% confidence interval [CI] = 0.84-1.02, p = 0.15; OR = 0.94, 95% CI = 0.55-1.58, p = 0.80, respectively). These findings remained consistent using different Mendelian randomization methods and in sensitivity analyses. For the IL-1-IL-6 pathway, causal estimates for IL-6 (OR = 0.86, 95% CI 0.81-0.91, p < 0.001), but not for IL-1β, IL-1ra, or sIL-6r, were significant. However, the association between genetically determined IL-6 and HF risk became non-significant after excluding SNPs with potential pleiotropy (OR = 0.89, 95% CI = 0.77-1.03, p = 0.12). Conclusion: Our study did not identify convincing evidence to support that CRP and fibrinogen, together with their upstream IL-1-IL-6 signaling pathway, were causally associated with HF risk.

Human immune diversity: from evolution to modernity

The extreme diversity of the human immune system, forged and maintained throughout evolutionary history, provides a potent defense against opportunistic pathogens. At the same time, this immune variation is the substrate upon which a plethora of immune-associated diseases develop. Genetic analysis suggests that thousands of individually weak loci together drive up to half of the observed immune variation. Intense selection maintains this genetic diversity, even selecting for the introgressed Neanderthal or Denisovan alleles that have reintroduced variation lost during the out-of-Africa migration. Variations in age, sex, diet, environmental exposure, and microbiome each potentially explain the residual variation, with proof-of-concept studies demonstrating both plausible mechanisms and correlative associations. The confounding interaction of many of these variables currently makes it difficult to assign definitive contributions. Here, we review the current state of play in the field, identify the key unknowns in the causality of immune variation, and identify the multidisciplinary pathways toward an improved understanding.

Genome-wide association study revealed novel candidate gene loci associated with soluble E-selectin levels in a Taiwanese population

BACKGROUND AND AIMS

Increase soluble E-selectin (sE-selectin) levels are associated with various inflammation and cardiometabolic disorders.

METHODS

This study aimed to investigate the genetic determinants of circulating sE-selectin levels by genome-wide association study (GWAS) in 4,525 Taiwan Biobank (TWB) participants and genotype-phenotype association analysis for sE-selectin level-determining alleles in over 80,000 TWB participants.

RESULTS

By GWAS, ABO, SELE, and FUT6 gene variants were identified as the determinants of sE-selectin levels, which reach genome-wide significance (maximum p = 3.25 × 10-271, 4.81 × 10-14, and 9.64 × 10-12, respectively). After further adjustment for the lead ABO rs2519093 genotypes, three novel gene loci, EVI5, FER and DMAC1, were associated with sE-selectin levels at p < 5 × 10-7. Three other previously reported gene loci, CELSR2, ST3GAL6-AS1, and HNF1A-AS1, also showed supportive evidence for the association with sE-selectin levels (maximum p < 0.0073). A multivariate analysis revealed age, body mass index, current smoking, hemoglobin A1C, hematocrit, leukocyte and platelet counts, serum alanine aminotransferase, triglycerides, and uric acid levels were independently associated with sE-selectin levels, in which the above ten gene loci contribute to 27.68% of the variance. For genotype-phenotype association analysis, a pleiotropic effect was demonstrated with genome-wide significant association between ABO gene variants and total and low-density-lipoprotein cholesterol levels, leukocyte counts and hematocrit.

CONCLUSIONS

Our data provide novel insight into the regulation of sE-selectin levels. These results may open new avenues in understanding the critical role of E-selectin on the pathogenesis of inflammatory and cardiometabolic disorders.

Genetic variant in 3' untranslated region of the mouse pycard gene regulates inflammasome activity

Quantitative trait locus mapping for interleukin-1β release after inflammasome priming and activation was performed on bone-marrow-derived macrophages (BMDM) from an AKRxDBA/2 mouse strain intercross. The strongest associated locus mapped very close to the Pycard gene on chromosome 7, which codes for the inflammasome adaptor protein apoptosis-associated speck-like protein containing a CARD (ASC). The DBA/2 and AKR Pycard genes only differ at a single-nucleotide polymorphism (SNP) in their 3' untranslated region (UTR). DBA/2 vs. AKR BMDM had increased levels of Pycard mRNA expression and ASC protein, and increased inflammasome speck formation, which was associated with increased Pycard mRNA stability without an increased transcription rate. CRISPR/Cas9 gene editing was performed on DBA/2 embryonic stem cells to change the Pycard 3'UTR SNP from the DBA/2 to the AKR allele. This single base change significantly reduced Pycard expression and inflammasome activity after cells were differentiated into macrophages due to reduced Pycard mRNA stability.

Association between ABO and Duffy blood types and circulating chemokines and cytokines

Blood group antigens are inherited traits that may play a role in immune and inflammatory processes. We investigated associations between blood groups and circulating inflammation-related molecules in 3537 non-Hispanic white participants selected from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Whole-genome scans were used to infer blood types for 12 common antigen systems based on well-characterized single-nucleotide polymorphisms. Serum levels of 96 biomarkers were measured on multiplex fluorescent bead-based panels. We estimated marker associations with blood type using weighted linear or logistic regression models adjusted for age, sex, smoking status, and principal components of population substructure. Bonferroni correction was used to control for multiple comparisons, with two-sided p values < 0.05 considered statistically significant. Among the 1152 associations tested, 10 were statistically significant. Duffy blood type was associated with levels of CXCL6/GCP2, CXCL5/ENA78, CCL11/EOTAXIN, CXCL1/GRO, CCL2/MCP1, CCL13/MCP4, and CCL17/TARC, whereas ABO blood type was associated with levels of sVEGFR2, sVEGFR3, and sGP130. Post hoc pairwise t-tests showed that individuals with type Fy(a+b−) had the lowest mean levels of all Duffy-associated markers, while individuals with type A blood had the lowest mean levels of all ABO-associated markers. Additional work is warranted to explore potential clinical implications of these differences.

Associations of Arachidonic Acid Synthesis with Cardiovascular Risk Factors and Relation to Ischemic Heart Disease and Stroke: A Univariable and Multivariable Mendelian Randomization Study

Arachidonic acid (AA), a major long-chain omega-6 polyunsaturated fatty acid, is associated with ischemic heart disease (IHD) and stroke. We assessed bi-directional associations of AA synthesis reflected by plasma phospholipid AA with CVD risk factors, and identified mediators of associations of AA with IHD and stroke using Mendelian randomization (MR). We used two-sample MR to assess bi-directional associations of AA synthesis with lipids, blood pressure, adiposity, and markers of inflammation and coagulation. We used multivariable MR to assess mediators of associations of AA with IHD and stroke. Genetically predicted AA (% of total fatty acids increase) was positively associated with apolipoprotein B (ApoB, 0.022 standard deviations (SD), 95% confidence interval (CI) 0.010, 0.034), high-density (0.030 SD, 95% CI 0.012, 0.049) and low-density lipoprotein cholesterol (LDL-C, 0.016 SD, 95% CI 0.004, 0.027) and lower triglycerides (-0.031 SD, 95% CI -0.049, -0.012) but not with other traits. Genetically predicted these traits gave no association with AA. The association of AA with IHD was attenuated adjusting for ApoB or LDL-C. Genetically predicted AA was associated with lipids but not other traits. Given ApoB is thought to be the key lipid in IHD, the association of AA with IHD is likely mediated by ApoB.

Circulating Interleukins and Risk of Multiple Sclerosis: A Mendelian Randomization Study

Background

Previous research have implicated critical roles of systemic inflammation in the development of Multiple Sclerosis (MS). But the causal relationship between interleukins (ILs) and MS has not been fully elucidated.

Objective

In this study, we applied Mendelian randomization (MR) approaches to address the causal associations between genetically determined circulating levels of ILs and the risk of MS.

Methods

Genetic instruments for circulating IL-1 receptor antagonist (IL-1Ra), IL-2 receptor α subunit (IL-2Rα), IL-6, IL-16, IL-17, and IL-18 were obtained from recently published genome-wide association studies (GWAS). Summary-level data for MS were obtained from the International Multiple Sclerosis Genetics Consortium. MR analyses were performed using the R software (version 3.6.1, The R Foundation) and the TwoSampleMR package.

Results

Genetic predisposition to higher circulating levels of IL-2Rα were significantly associated with MS risk. The odds ratio (OR) was 1.22 (95% confidence interval [CI], 1.12-1.32; p < 0.001) per one standard deviation increase in circulating IL-2Rα levels. There was a suggestive association of circulating IL-1Ra with MS risk (OR, 0.94; 95% CI, 0.88-0.99; p = 0.027). The other ILs were not associated with the outcome.

Conclusion

Our results indicated that circulating IL-2Rα was causally associated with risk of MS.

A phenome-wide association study of ABO blood groups

BACKGROUND

ABO blood group is associated with differences in lifespan, cardiovascular disease, and some cancers, for reasons which are incompletely understood. To gain sex-specific additional insight about potential mechanisms driving these common conditions for future interventions, we characterized associations of ABO blood group antigen across the phenotype sex-specifically.

METHODS

We performed a phenome-wide association study (PheWAS) assessing the association of tag single nucleotide polymorphisms (SNPs) for ABO blood group antigens (O, B, A1, and A2) with 3873 phenotypes.

RESULTS

The tag SNP for the O antigen was inversely associated with diseases of the circulatory system (particularly deep vein thrombosis (DVT)), total cholesterol, low-density lipoprotein cholesterol (LDL-C), and ovarian cancer, and positively associated with erythrocyte traits, leukocyte counts, diastolic blood pressure (DBP), and healthy body composition; the tag SNP for the A1 antigen tended to have associations in reverse to O. Stronger associations were more apparent for men than women for DVT, DBP, leukocyte traits, and some body composition traits, whereas larger effect sizes were found for women than men for some erythrocyte and lipid traits.

CONCLUSION

Blood group has a complex association with cardiovascular diseases and its major risk factors, including blood pressure and lipids, as well as with blood cell traits and body composition, with some differences by sex. Lower LDL-C may underlie some of the benefits of blood group O, but the complexity of associations with blood group antigen suggests overlooked drivers of common chronic diseases.

Genetics meets proteomics: perspectives for large population-based studies

Proteomic analysis of cells, tissues and body fluids has generated valuable insights into the complex processes influencing human biology. Proteins represent intermediate phenotypes for disease and provide insight into how genetic and non-genetic risk factors are mechanistically linked to clinical outcomes. Associations between protein levels and DNA sequence variants that colocalize with risk alleles for common diseases can expose disease-associated pathways, revealing novel drug targets and translational biomarkers. However, genome-wide, population-scale analyses of proteomic data are only now emerging. Here, we review current findings from studies of the plasma proteome and discuss their potential for advancing biomedical translation through the interpretation of genome-wide association analyses. We highlight the challenges faced by currently available technologies and provide perspectives relevant to their future application in large-scale biobank studies.

Could vitamin D reduce obesity-associated inflammation? Observational and Mendelian randomization study

BACKGROUND

Obesity is associated with inflammation but the role of vitamin D in this process is not clear.

OBJECTIVES

We aimed to assess the associations between serum 25-hydroxyvitamin D [25(OH)D], BMI, and 16 inflammatory biomarkers, and to assess the role of vitamin D as a potential mediator in the association between higher BMI and inflammation.

METHODS

Northern Finland Birth Cohort 1966 (NFBC1966) 31-y data on 3586 individuals were analyzed to examine the observational associations between BMI, 25(OH)D, and 16 inflammatory biomarkers. Multivariable regression analyses and 2-sample regression-based Mendelian randomization (MR) mediation analysis were performed to assess any role of vitamin D in mediating a causal effect of BMI on inflammatory biomarkers [soluble intercellular adhesion molecule 1 (sICAM-1), high sensitivity C-reactive protein (hs-CRP), and α1-acid glycoprotein (AGP)] for which observational associations were detected. For MR, genome-wide association study summary results ranging from 5163 to 806,834 individuals were used for biomarkers, 25(OH)D, and BMI. Findings were triangulated with a literature review of vitamin D supplementation trials.

RESULTS

In NFBC1966, mean BMI (kg/m2) was 24.8 (95% CI: 24.7, 25.0) and mean 25(OH)D was 50.3 nmol/L (95% CI: 49.8, 50.7 nmol/L). Inflammatory biomarkers correlated as 4 independent clusters: interleukins, adhesion molecules, acute-phase proteins, and chemokines. BMI was positively associated with 9 inflammatory biomarkers and inversely with 25(OH)D (false discovery rate < 0.05). 25(OH)D was inversely associated with sICAM-1, hs-CRP, and AGP, which were positively associated with BMI. The MR analyses showed causal association of BMI on these 3 inflammatory biomarkers. There was no observational or MR evidence that circulating 25(OH)D concentrations mediated the association between BMI and these 3 inflammatory markers. Review of randomized controlled trials (RCTs) supported our findings showing no impact of vitamin D supplementation on inflammatory biomarkers.

CONCLUSIONS

The findings from our observational study and causal MR analyses, together with data from RCTs, do not support a beneficial role of vitamin D supplementation on obesity-related inflammation.

Circulating interleukins in relation to coronary artery disease, atrial fibrillation and ischemic stroke and its subtypes: A two-sample Mendelian randomization study

Whether thyroid dysfunction plays a causal role in the development of cancer remains inconclusive. We conducted a two-sample Mendelian randomization study to investigate the associations between genetic predisposition to thyroid dysfunction and 22 site-specific cancers. Single-nucleotide polymorphisms associated with four traits of thyroid function were selected from a genome-wide association meta-analysis with up to 72,167 European-descent individuals. Summary-level data for breast cancer and 21 other cancers were extracted from the Breast Cancer Association Consortium (122,977 breast cancer cases and 105,974 controls) and UK Biobank (367,643 individuals). For breast cancer, a meta-analysis was performed using data from both sources. Genetically predicted thyroid dysfunction was associated with breast cancer, with similar patterns of associations in the Breast Cancer Association Consortium and UK Biobank. The combined odds ratios of breast cancer were 0.94 (0.91–0.98; p = 0.007) per genetically predicted one standard deviation increase in TSH levels, 0.96 (0.91–1.00; p = 0.053) for genetic predisposition to hypothyroidism, 1.04 (1.01–1.07; p = 0.005) for genetic predisposition to hyperthyroidism and 1.07 (1.02–1.12; p = 0.003) per genetically predicted one standard deviation increase in free thyroxine levels. Genetically predicted TSH levels and hypothyroidism were inversely with thyroid cancer; the odds ratios were 0.47 (0.30-0.73; p = 0.001) and 0.70 (0.51-0.98; p = 0.038), respectively. Our study provides evidence of a causal association between thyroid dysfunction and breast cancer (mainly ER-positive tumors) risk. The role of TSH and hypothyroidism for thyroid cancer and the associations between thyroid dysfunction and other cancers need further exploration.

Multivariate Genome-wide Association Analysis of a Cytokine Network Reveals Variants with Widespread Immune, Haematological, and Cardiometabolic Pleiotropy

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

Mendelian randomization indicates that TNF is not causally associated with Alzheimer's disease

Epidemiological research has suggested that inhibition of tumor necrosis factor (TNF)-α in patients with rheumatoid arthritis (RA) reduces the overall risk of Alzheimer's disease (AD). TNF-α antagonists have been suggested as a potential treatment for AD. We used a two-sample Mendelian randomization design to examine the causal relationship between blood TNF expression, serum TNF-α levels, and RA on AD risk. Our results do not support a causal relationship between TNF expression, serum TNF-α levels, and RA on AD risk.