Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults

Exploring the role of CD13 and inflammatory factors in radiation enteritis: insights from high-throughput proteomics and Mendelian randomization analysis

BACKGROUND

Radiation enteritis (RE) is an unavoidable complication during radiotherapy for pelvic malignancies, characterized by chronic inflammation, fibrosis, and vascular injury in the intestinal tissue. Currently, there is a lack of research that delves into the relationship between inflammatory factors and key proteins in RE.

METHODS

This study employed high-throughput proteomics to analyze intestinal tissues from RE rats and healthy controls, identifying differentially expressed key proteins. The degree of intestinal damage was validated through HE staining. Furthermore, five Mendelian randomization methods were used to analyze the causal relationship between 70 serum circulating inflammatory factors and CD13 levels. Sensitivity analyses, including heterogeneity tests, leave-one-out tests, and horizontal pleiotropy tests, were performed to ensure the robustness and reliability of the results.

RESULTS

CD13 was identified as a key differentially expressed protein, with its expression significantly upregulated in RE rats and positively correlated with disease severity. Bidirectional Mendelian randomization analysis revealed causal relationships between CD13 and four inflammatory factors: increased levels of CCL28 and EN-RAGE may promote the rise in CD13, while increased levels of TAM-binding protein may be associated with decreased CD13 levels. Additionally, higher CD13 levels were found to be associated with increased levels of interleukin-12. Sensitivity analyses indicated good consistency and reliability in terms of heterogeneity and pleiotropy for these exposure variables.

CONCLUSION

This study reveals the potential mechanistic role of CD13 in RE. Moreover, the identified CD13-associated inflammatory factors offer potential targets for the development of new prevention and treatment strategies, with significant clinical implications.

Influence of Altered Immune-Inflammatory Axis on the Risk of Osteomyelitis and Its Network Interaction Effect in European Population

Background: Osteomyelitis (OM) is a severe bone infection with rising incidence rates, particularly among elderly and diabetic patients. The immune-inflammatory axis is implicated in OM pathogenesis, but its complex interplay remains poorly defined. This study aims to explore the causal relationships between immunophenotypes, plasma inflammatory proteins, and OM using Mendelian randomization (MR) and bioinformatics. Methods: Utilizing publicly available genetic data, we undertook a series of quality control measures to identify instrumental variables (IVs) associated with exposure. Subsequently, we conducted MR using inverse variance weighting to explore the causal relationships between immunophenotypes, plasma inflammatory proteins, and OM. Bioinformatics tools were applied to explore the functional enrichment and protein-protein interaction networks of the implicated genes. Results: The MR analysis identified 13 immune cell phenotypes and 2 plasma inflammatory proteins associated with risk of OM. Notably, higher levels of HLA DR on plasmacytoid dendritic cells, memory B cell absolute counts, and CD8dim T cell percentages were associated with increased OM risk. Additionally, elevated levels of CD6 and IL-12 subunit B were correlated with OM risk. Bioinformatics analysis revealed the enrichment of related genes in immune-related pathways and highlighted the complex interaction networks of the implicated proteins. Conclusions: This study provides novel insights into the immune-inflammatory axis in OM and identifies potential biomarkers for risk assessment. The findings warrant further validation in diverse populations and pave the way for developing targeted preventive and therapeutic strategies for OM.

DNA methylation-predicted plasma protein levels and breast cancer risk

BACKGROUND

Blood DNA methylation (DNAm) profiles have been used to show that changes in circulating leukocyte composition occur during breast cancer development, suggesting that peripheral immune system alterations are markers of breast cancer risk. Blood DNAm profiles have recently been used to predict plasma protein concentrations ("Protein EpiScores"), but their associations with breast cancer risk have not been examined in detail.

METHODS

Whole blood DNAm profiles were obtained for a case-cohort sample of participants in the Sister Study and used to calculate 109 Protein EpiScores. Of the 4,479 women included, 2,151 (48%) were diagnosed with breast cancer within 15 years of their baseline blood draw (median time to diagnosis: 8.6 years; 1,673 invasive cancer and 478 ductal carcinomas in situ). Protein EpiScores associations with breast cancer incidence were estimated using weighted Cox regression models, overall and stratified by time and participant characteristics.

RESULTS

Protein EpiScores for RARRES2, IGFBP4, and CCL21 were positively associated with invasive breast cancer risk (hazard ratios from 1.17 to 1.24), while those for F7, SELL, CXCL9, CD48, and IL19 were inversely associated (hazard ratios from 0.82 to 0.86) (all FDR < 0.10). Eight immune response-related Protein EpiScores (CXCL9, CD48, FCGR3B, CXCL11, CCL21, CRTAM, VCAM1, GZMA) were associated with invasive cancers diagnosed within five years of enrollment. Protein EpiScore associations were consistently stronger for estrogen receptor-negative tumors.

CONCLUSIONS

Several Protein EpiScores, including many related to immune response, were associated with breast cancer risk, highlighting novel changes to the peripheral immune system that occur during breast cancer development.

Bivariate genome-wide association study of circulating fibrinogen and C-reactive protein levels

BACKGROUND

Fibrinogen and C-reactive protein (CRP) play an important role in inflammatory pathways and share multiple genetic loci reported in previously published genome-wide association studies (GWAS), highlighting their common genetic background. Leveraging the shared biology may identify further loci pleiotropically associated with both fibrinogen and CRP.

OBJECTIVES

To identify novel genetic variants that are pleiotropic and associated with both fibrinogen and CRP, by integrating both phenotypes in a bivariate GWAS by using a multitrait GWAS.

METHODS

We performed a bivariate GWAS to identify further pleiotropic genetic loci, using summary statistics of previously published GWAS on fibrinogen (n = 120 246) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium, consisting of European ancestry samples and CRP (n = 363 228) from UK Biobank, including 5 different population groups. The main analysis was performed using metaUSAT and N-GWAMA. We conducted replication for novel CRP associations to test the robustness of the findings using an independent GWAS for CRP (n = 148 164). We also performed colocalization analysis to compare the associations in identified loci for the 2 traits and Genotype-Tissue Expression data.

RESULTS

We identified 87 pleiotropic loci that overlapped between metaUSAT and N-GWAMA, including 23 previously known for either fibrinogen or CRP, 58 novel loci for fibrinogen, and 6 novel loci for both fibrinogen and CRP. Overall, there were 30 pleiotropic and novel loci for both traits, and 7 of these showed evidence of colocalization, located in or near ZZZ3, NR1I2, RP11-72L22.1, MICU1, ARL14EP, SOCS2, and PGM5. Among these 30 loci, 13 replicated for CRP in an independent CRP GWAS.

CONCLUSION

Bivariate GWAS identified additional associated loci for fibrinogen and CRP. This analysis suggests fibrinogen and CRP share a common genetic architecture with many pleiotropic loci.

Evaluating the association between immunological proteins and common intestinal diseases using a bidirectional two-sample Mendelian randomization study

Dysregulation of intestinal homeostasis, characterized by imbalanced immunological proteins, contributes to the pathogenesis of common intestinal diseases, e.g., irritable bowel syndrome (IBS), inflammatory bowel disease (IBD), and colorectal cancer (CRC). However, the potential causal relationships between specific immunological proteins and these diseases remain to be fully elucidated. In this study, we employed the bidirectional two-sample Mendelian randomization analysis to infer potential causal relationships between representative immunological proteins and these intestinal diseases. Genome-wide association study (GWAS) summary statistics of IBS, IBD, and CRC were obtained from public databases and utilized in MR analysis. Multiple sensitivity analyses were performed to evaluate the robustness, with p-values adjusted using the Benjamini-Hochberg method for multiple comparisons. Our findings revealed a significant association between IL-1β (OR = 0.783, 95 % CI: 0.676 to 0.908, adjusted P = 0.048) and a decreased risk of IBS. Furthermore, genetic predisposition to IBS was related to the reduced levels of IL-25 (β =  - 0.233, 95 % CI: -0.372 to -0.094, adjusted P = 0.047). Additionally, genetic predisposition to IBD was correlated with elevated levels of IL-6 (β = 0.046, 95 % CI: 0.022-0.069, adjusted P = 0.010). The levels of TNF-α (OR = 1.252, 95 % CI: 1.102 to 1.423, adjusted P = 0.047) were associated with an increased risk of CRC. Our study suggests associations between specific immunological proteins and intestinal diseases, which would provide valuable insights for developing targeted immunomodulation therapies for these conditions. Further investigation into underlying mechanisms remains a research priority in the future.

The Genetic Variants Influencing Hypertension Prevalence Based on the Risk of Insulin Resistance as Assessed Using the Metabolic Score for Insulin Resistance (METS-IR)

Insulin resistance is a major indicator of cardiovascular diseases, including hypertension. The Metabolic Score for Insulin Resistance (METS-IR) offers a simplified and cost-effective way to evaluate insulin resistance. This study aimed to identify genetic variants associated with the prevalence of hypertension stratified by METS-IR score levels. Data from the Korean Genome and Epidemiology Study (KoGES) were analyzed. The METS-IR was calculated using the following formula: ln [(2 × fasting blood glucose (FBG) + triglycerides (TG)) × body mass index (BMI)]/ ln [high-density lipoprotein cholesterol (HDL-C)]. The participants were divided into tertiles 1 (T1) and 3 (T3) based on their METS-IR scores. Genome-wide association studies (GWAS) were performed for hypertensive cases and non-hypertensive controls within these tertile groups using logistic regression adjusted for age, sex, and lifestyle factors. Among the METS-IR tertile groups, 3517 of the 19,774 participants (17.8%) at T1 had hypertension, whereas 8653 of the 20,374 participants (42.5%) at T3 had hypertension. A total of 113 single-nucleotide polymorphisms (SNPs) reached the GWAS significance threshold (p < 5 × 10-8) in at least one tertile group, mapping to six distinct genetic loci. Notably, four loci, rs11899121 (chr2p24), rs7556898 (chr2q24.3), rs17249754 (ATP2B1), and rs1980854 (chr20p12.2), were significantly associated with hypertension in the high-METS-score group (T3). rs10857147 (FGF5) was significant in both the T1 and T3 groups, whereas rs671 (ALDH2) was significant only in the T1 group. The GWASs identified six genetic loci significantly associated with hypertension, with distinct patterns across METS-IR tertiles, highlighting the role of metabolic context in genetic susceptibility. These findings underscore critical genetic factors influencing hypertension prevalence and provide insights into the metabolic-genetic interplay underlying this condition.

Exploring causal associations between interleukins, their receptors, and acute respiratory distress syndrome: a bidirectional-Mendelian Randomization study

Background

Acute Respiratory Distress Syndrome (ARDS) is a severe lung condition posing significant health risks. Observational studies have indicated a potential linkage between ARDS and various interleukins, such as Interleukin-8, -2, -6, and -1β. However, the causal relationships between specific interleukins, their receptors, and ARDS remain unclear.

Methods

This study employed bidirectional Mendelian Randomization (MR) to investigate the causal associations between 197 different interleukins, interleukin receptors, and ARDS. GWAS summary data for interleukins and their receptors were sourced from publicly available studies, while ARDS data were derived from the FinnGen database. Wald Ratio and Inverse Variance Weighted (IVW) methods were primarily utilized in the MR analyses. Reverse MR was also conducted to explore reverse causations.

Results

The study identified significant causal relationships between several interleukins, their receptors and ARDS. Specifically, Interleukin-17 receptor D (IL-17RD) (OR = 0.67, 95%CI = 0.51-0.89, p = 0.006), Interleukin-22 receptor subunit alpha-1 (IL22RA-1) (OR = 1.45, 95%CI = 1.08-1.96, p = 0.014), Interleukin-20 receptor subunit alpha (IL20RA) (OR = 0.53, 95%CI = 0.32-0.89, p = 0.016), Interleukin-22 (IL-22) (OR = 3.60, 95%CI = 1.20-10.83, p = 0.022), Interleukin-23 receptor (IL-23R) (OR = 2.14, 95%CI = 1.10-4.17, p = 0.025), Single Ig IL-1-related receptor (SIGIRR) (OR = 1.22, 95%CI = 1.00-1.48, p = 0.047) showed notable associations with ARDS. No reverse causal relationships were found on results above, and neither heterogeneity nor pleiotropy was detected in the analysis.

Conclusion

This study elucidates the causal connections between specific interleukins, their receptors and ARDS, contributing significantly to the understanding of ARDS pathogenesis. These findings offer a foundation for further research and potential therapeutic interventions targeting these interleukins and receptors in ARDS management. The absence of reverse causation and pleiotropy, heterogeneity reinforces the robustness of these associations.

Mendelian randomization analysis using multiple biomarkers of an underlying common exposure

Mendelian randomization (MR) analysis is increasingly popular for testing the causal effect of exposures on disease outcomes using data from genome-wide association studies. In some settings, the underlying exposure, such as systematic inflammation, may not be directly observable, but measurements can be available on multiple biomarkers or other types of traits that are co-regulated by the exposure. We propose a method for MR analysis on latent exposures (MRLE), which tests the significance for, and the direction of, the effect of a latent exposure by leveraging information from multiple related traits. The method is developed by constructing a set of estimating functions based on the second-order moments of GWAS summary association statistics for the observable traits, under a structural equation model where genetic variants are assumed to have indirect effects through the latent exposure and potentially direct effects on the traits. Simulation studies show that MRLE has well-controlled type I error rates and enhanced power compared to single-trait MR tests under various types of pleiotropy. Applications of MRLE using genetic association statistics across five inflammatory biomarkers (CRP, IL-6, IL-8, TNF-α, and MCP-1) provide evidence for potential causal effects of inflammation on increasing the risk of coronary artery disease, colorectal cancer, and rheumatoid arthritis, while standard MR analysis for individual biomarkers fails to detect consistent evidence for such effects.

Phenotype wide association study links bronchopulmonary dysplasia with eosinophilia in children

Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth. Despite this, genetic drivers of BPD are poorly understood. The objective of this study is to better understand the impact of single nucleotide polymorphisms (SNPs) previously associated with BPD by examining associations with other phenotypes. We drew pediatric subjects from the biorepository at the Center for Applied Genomics to identify associations between these SNPs and 2,146 imputed phenotypes. Methylation data, external cohorts, and in silico validation methods were used to corroborate significant associations. We identified 60 SNPs that were previously associated with BPD. We found a significant association between rs3771150 and rs3771171 and mean eosinophil percentage in a European cohort of 6,999 patients and replicated this in external cohorts. Both SNPs were also associated with asthma, COPD and FEV1/FVC ratio. These SNPs displayed associations with methylation probes and were functionally linked to ST2 (IL1RL1) levels in blood and lung tissue. Our findings support a genetic justification for the epidemiological link between BPD and asthma. Given the well-established link between ST2 and type 2 inflammation in asthma, these findings provide a rationale for future studies exploring the role of type 2 inflammation in the pathogenesis of BPD.

Exploring the relationship between the interleukin family and lung adenocarcinoma through Mendelian randomization and RNA sequencing analysis

BACKGROUND

Lung adenocarcinoma (LUAD) is still one of the most prevalent malignancies. Interleukin factors are closely associated with the initiation and progression of cancer. However, the relationship between interleukin factors and LUAD has not been fully elucidated. This study aimed to use Mendelian randomization (MR) and RNA sequencing (RNA-seq) analyses to identify the interleukin factors associated with the onset and progression of LUAD.

METHODS

Exposure-related instrumental variables were selected from interleukin factor summary datasets. The LUAD summary dataset from FINGENE served as the outcome. MR and sensitivity analyses were conducted to screen for interleukin factors associated with LUAD occurrence. Transcriptome analyses revealed the role of interleukin factors in lung tissues. The results were validated through Western blotting and further confirmed with driver gene-negative patients from multiple centers. Potential mechanisms influencing LUAD occurrence and development were explored using bulk RNA-seq and single-cell RNA-seq data.

RESULTS

MR analysis indicated that elevated plasma levels of IL6RB, IL27RA, IL22RA1, and IL16 are causally associated with increased LUAD risk, while IL18R1 and IL11RA exhibit the opposite effect. Transcriptome analyses revealed that IL11RA, IL18R1, and IL16 were downregulated in tumor tissues compared with normal lung tissue, but only higher expression of IL11RA correlated with improved prognosis in patients with LUAD from different centers and persisted even in driver-gene negative patients. The IL11RA protein level was lower in various LUAD cell lines than in human bronchial epithelial cells. The genes co-expressed with IL11RA were enriched in the Ras signaling pathway and glycosylation processes. Fibroblasts were the primary IL11RA-expressing cell population, with IL11RA+fibroblasts exhibiting a more immature state. The genes differentially expressed between IL11RA+and IL11RA- fibroblasts were involved in the PI3K-Akt/TNF signaling pathway.

CONCLUSION

According to the MR and transcriptome analyses, the downregulation of IL11RA was closely related to the occurrence and development of LUAD.

DNAm scores for serum GDF15 and NT-proBNP levels associate with a range of traits affecting the body and brain

BACKGROUND

Plasma growth differentiation factor 15 (GDF15) and N-terminal proB-type natriuretic peptide (NT-proBNP) are cardiovascular biomarkers that associate with a range of diseases. Epigenetic scores (EpiScores) for GDF15 and NT-proBNP may provide new routes for risk stratification.

RESULTS

In the Generation Scotland cohort (N ≥ 16,963), GDF15 levels were associated with incident dementia, ischaemic stroke and type 2 diabetes, whereas NT-proBNP levels were associated with incident ischaemic heart disease, ischaemic stroke and type 2 diabetes (all PFDR < 0.05). Bayesian epigenome-wide association studies (EWAS) identified 12 and 4 DNA methylation (DNAm) CpG sites associated (Posterior Inclusion Probability [PIP] > 95%) with levels of GDF15 and NT-proBNP, respectively. EpiScores for GDF15 and NT-proBNP were trained in a subset of the population. The GDF15 EpiScore replicated protein associations with incident dementia, type 2 diabetes and ischaemic stroke in the Generation Scotland test set (hazard ratios (HR) range 1.36-1.41, PFDR < 0.05). The EpiScore for NT-proBNP replicated the protein association with type 2 diabetes, but failed to replicate an association with ischaemic stroke. EpiScores explained comparable variance in protein levels across both the Generation Scotland test set and the external LBC1936 test cohort (R2 range of 5.7-12.2%). In LBC1936, both EpiScores were associated with indicators of poorer brain health. Neither EpiScore was associated with incident dementia in the LBC1936 population.

CONCLUSIONS

EpiScores for serum levels of GDF15 and Nt-proBNP associate with body and brain health traits. These EpiScores are provided as potential tools for disease risk stratification.

The role of serum interleukins in Cancer: A Multi-center Mendelian Randomization study

The occurrence of cancer is often accompanied by immune evasion and tumor-promoting inflammation, with interleukins (IL) playing a pivotal role in the immune-inflammatory mechanism. However, the precise contribution of serum interleukins in cancer remains elusive. We obtained GWAS summary data for 35 interleukins from eight independent large-scale serum proteome studies of European ancestry populations and for 23 common cancers from the FinnGen Consortium. We then conducted a multicenter Mendelian Randomization (MR) study to explore the relationship between systemic inflammatory status and cancers. 24 causal associations between interleukins and cancers were supported by multicenter data, 18 of which were reported for the first time. Our results indicated that IL-1α (Hodgkin lymphoma), IL-5 (bladder cancer), IL-7 (prostate cancer), IL-11 (bone malignant tumor), IL-16 (lung cancer), IL-17A (pancreatic cancer), IL-20 (bladder cancer), IL-22 (lymphocytic leukemia), IL-34 (breast cancer), IL-36β (prostate cancer), and IL-36γ (liver cancer) were risk factors for related cancers. Conversely, IL-9 (malignant neoplasms of the corpus uteri), IL-17C (liver cancer), and IL-31 (colorectal cancer, bladder cancer, pancreatic cancer, and cutaneous melanoma) exhibited protective effects against related cancers. Notably, the dual effects of serum interleukins were also observed. IL-18 acted as a risk factor for prostate cancer, however, was a protective factor against laryngeal cancer. Similarly, IL-19 promoted the development of lung cancer and myeloid leukemia, while conferring protection against Breast, cervical, and thyroid cancers. Our study confirmed the genetic association between multiple serum interleukins and cancers. Immune and anti-inflammatory strategies targeting these associations provide opportunities for prevention and treatment.

Causality between various cytokines and asthma: a bidirectional two-sample Mendelian randomization analysis

Background

Many studies have shown that cytokines play an important role in the pathogenesis of asthma, but their biological effects on asthma remain unclear. The Mendelian randomization (MR) method was used to evaluate the causal relationship between various cytokines [such as interleukins (ILs), interferons (IFNs), tumor necrosis factors (TNFs), colony-stimulating factors (CSFs), transforming growth factor (TGF), etc.,] and asthma.

Methods

In this study, inverse variance weighting was used to evaluate the causal relationship between asthma and cytokines. In addition, the reliability of the results is ensured by multiple methods such as MR-Egger, weighted median, MR-Raps, MR-Presso, and RadialMR, as well as sensitivity analysis.

Results

The results showed that none of the 11 cytokines was associated with the risk of asthma. In contrast, asthma can increase levels of IL-5 [odds ratio (OR) = 1.112, 95% confidence interval (CI): 1.009-1.224, P = 0.032] and IL-9 (OR = 1.111, 95% CI: 1.013-1.219, P = 0.025).

Conclusion

Genetically predicted asthma was positively associated with elevated levels of IL-5 and IL-9, indicating the downstream effects of IL-5 and IL-9 on asthma. Medical treatments can thus be designed to target IL-5 and IL-9 to prevent asthma exacerbations.

The interplay between cytokines and stroke: a bi-directional Mendelian randomization study

Stroke, the second leading cause of death and disability, causes massive cell death in the brain followed by secondary inflammatory injury initiated by disease associated molecular patterns released from dead cells. Nonetheless, the evidence regarding the causal relationship between inflammatory cytokines and stroke subtypes is obscure. To leverage large scale genetic association data to investigate the interplay between circulating cytokines and stroke, we adopted a two-sample bi-directional Mendelian randomization (MR) analysis. Firstly, we performed a forward MR analysis to examine the associations of genetically determined 31 cytokines with 6 stroke subtypes. Secondly, we conducted a reverse MR analysis to check the associations of 6 stroke subtypes with 31 cytokines. In the forward MR analysis, genetic evidence suggests that 21 cytokines were significantly associated with certain stroke subtype risk with |β| ranging from 1.90 × 10-4 to 0.74. In the reverse MR analysis, our results found that five stroke subtypes (intracerebral hemorrhage (ICH), large artery atherosclerosis ischemic stroke (LAAS), lacunar stroke (LS), cardioembolic ischemic stroke (CEI), small-vessel ischemic stroke (SV)) caused significantly changes in 16 cytokines with |β| ranging from 1.08 × 10-4 to 0.69. In particular, those five stroke subtypes were statistically significantly associated with C-reactive protein (CRP). In addition, ICH, LAAS, LS and SV were significantly correlated with vascular endothelial growth factor (VEGF), while LAAS, LS, CEI and SV were significantly related to fibroblast growth factor (FGF). Moreover, integrated bi-directional MR analysis, these factors (IL-3Rα, IL-6R, IL-6Rα, IL-1Ra, insulin-like growth factor-1(IGF-1), IL-12Rβ2) can be used as predictors of some specific stroke subtypes. As well as, IL-16 and C-C motif chemokine receptor 7 (CCR7) can be used as prognostic factors of stroke. Our findings prognostic identify potential pharmacological opportunities, including perturbation of circulating cytokines for both predicting stroke risk and post stroke treatment effects. As we conducted a comprehensive search and analysis of stroke subtype and cytokines in the existing publicly available GWAS database, the results have good population-generalizability.

Gut Subdoligranulum variabile ameliorates rheumatoid arthritis by promoting TSG-6 synthesis from joint cells

Background

A burgeoning body of evidence has substantiated the association between alterations in the composition of the gut microbiota and rheumatoid arthritis (RA). Nevertheless, our understanding of the intricate mechanisms underpinning this association is limited.

Methods

To investigate whether the gut microbiota influences the pathogenesis of RA through metabolism or immunity, we performed rigorous synthesis analyses using aggregated statistics from published genome-wide association studies (GWAS) using two-sample Mendelian randomization (MR) and mediated MR techniques, including two-step MR and multivariate MR analyses. Subsequently, we conducted in vitro cellular validation of the analyzed Microbial-Cytokine-RA pathway. We determined the optimal culture conditions through co-culture experiments involving concentration and time. Cell Counting Kit-8 (CCK-8) assays were employed to assess cellular viability, and enzyme-linked immunosorbent assays (ELISA) were performed to assess tumor necrosis factor-inducible gene 6 protein (TSG-6) and tumor necrosis factor-α (TNF-α) levels.

Results

Our univariable MR results confirmed 15 microbial traits, 7 metabolites and 2 cytokines that may be causally associated with RA (P FDR < 0.05). Mediation analysis revealed that microbial traits influence the risk of RA through metabolite or cytokine (proportion mediated: 7.75% - 58.22%). In vitro experiments demonstrated that TSG-6 was highly expressed in the Subdoligranulum variabile treatment group and was correlated with decreased RA severity (reduced TNF-α expression). Silencing the TSG-6 gene significantly increased TNF-α expression, regardless of treatment with S. variabile. Additionally, S. variabile-secreted exosomes exhibited the same effect.

Conclusion

The results of this study suggest that S. variabile has the potential to promote TSG-6 secretion, thereby reducing RA inflammation.

Single nucleotide variants in the CCL2, OAS1 and DPP9 genes and their association with the severity of COVID-19 in an Ecuadorian population

COVID-19 has a broad clinical spectrum, ranging from asymptomatic-mild form to severe phenotype. The severity of COVID-19 is a complex trait influenced by various genetic and environmental factors. Ethnic differences have been observed in relation to COVID-19 severity during the pandemic. It is currently unknown whether genetic variations may contribute to the increased risk of severity observed in Latin-American individuals The aim of this study is to investigate the potential correlation between gene variants at CCL2, OAS1, and DPP9 genes and the severity of COVID-19 in a population from Quito, Ecuador. This observational case-control study was conducted at the Carrera de Biologia from the Universidad Central del Ecuador and the Hospital Quito Sur of the Instituto Ecuatoriano de Seguridad Social (Quito-SUR-IESS), Quito, Ecuador. Genotyping for gene variants at rs1024611 (A>G), rs10774671 (A>G), and rs10406145 (G>C) of CCL2, OAS1, and DPP9 genes was performed on 100 COVID-19 patients (43 with severe form and 57 asymptomatic-mild) using RFLP-PCR. The genotype distribution of all SNVs throughout the entire sample of 100 individuals showed Hardy Weinberg equilibrium (P=0.53, 0.35, and 0.4 for CCL2, OAS1, and DPP9, respectively). The HWE test did not find any statistically significant difference in genotype distribution between the study and control groups for any of the three SNVs. The multivariable logistic regression analysis showed that individuals with the GG of the CCL2 rs1024611 gene variant had an increased association with the severe COVID-19 phenotype in a recessive model (P = 0.0003, OR = 6.43, 95% CI 2.19-18.89) and for the OAS1 rs10774671 gene variant, the log-additive model showed a significant association with the severe phenotype of COVID-19 (P=0.0084, OR=3.85, 95% CI 1.33-11.12). Analysis of haplotype frequencies revealed that the coexistence of GAG at CCL2, OAS1, and DPP9 variants, respectively, in the same individual increased the presence of the severe COVID-19 phenotype (OR=2.273, 95% CI: 1.271-4.068, P=0.005305). The findings of the current study suggests that the ethnic background affects the allele and genotype frequencies of genes associated with the severity of COVID-19. The experience with COVID-19 has provided an opportunity to identify an ethnicity-based approach to recognize genetically high-risk individuals in different populations for emerging diseases.

Proteogenomic Data Integration Reveals CXCL10 as a Potentially Downstream Causal Mediator for IL-6 Signaling on Atherosclerosis

BACKGROUND

Genetic and experimental studies support a causal involvement of IL-6 (interleukin-6) signaling in atheroprogression. Although trials targeting IL-6 signaling are underway, any benefits must be balanced against an impaired host immune response. Dissecting the mechanisms that mediate the effects of IL-6 signaling on atherosclerosis could offer insights about novel drug targets with more specific effects.

METHODS

Leveraging data from 522 681 individuals, we constructed a genetic instrument of 26 variants in the gene encoding the IL-6R (IL-6 receptor) that proxied for pharmacological IL-6R inhibition. Using Mendelian randomization, we assessed its effects on 3281 plasma proteins quantified with an aptamer-based assay in the INTERVAL cohort (n=3301). Using mediation Mendelian randomization, we explored proteomic mediators of the effects of genetically proxied IL-6 signaling on coronary artery disease, large artery atherosclerotic stroke, and peripheral artery disease. For significant mediators, we tested associations of their circulating levels with incident cardiovascular events in a population-based study (n=1704) and explored the histological, transcriptomic, and cellular phenotypes correlated with their expression levels in samples from human atherosclerotic lesions.

RESULTS

We found significant effects of genetically proxied IL-6 signaling on 70 circulating proteins involved in cytokine production/regulation and immune cell recruitment/differentiation, which correlated with the proteomic effects of pharmacological IL-6R inhibition in a clinical trial. Among the 70 significant proteins, genetically proxied circulating levels of CXCL10 (C-X-C motif chemokine ligand 10) were associated with risk of coronary artery disease, large artery atherosclerotic stroke, and peripheral artery disease, with up to 67% of the effects of genetically downregulated IL-6 signaling on these end points mediated by decreases in CXCL10. Higher midlife circulating CXCL10 levels were associated with a larger number of cardiovascular events over 20 years, whereas higher CXCL10 expression in human atherosclerotic lesions correlated with a larger lipid core and a transcriptomic profile reflecting immune cell infiltration, adaptive immune system activation, and cytokine signaling.

CONCLUSIONS

Integrating multiomics data, we found a proteomic signature of IL-6 signaling activation and mediators of its effects on cardiovascular disease. Our analyses suggest the interferon-γ-inducible chemokine CXCL10 to be a potentially causal mediator for atherosclerosis in 3 vascular compartments and, as such, could serve as a promising drug target for atheroprotection.

Sex-biased genetic regulation of inflammatory proteins in the Dutch population

BACKGROUND

Significant differences in immune responses, prevalence or susceptibility of diseases and treatment responses have been described between males and females. Despite this, sex-differentiation analysis of the genetic architecture of inflammatory proteins is largely unexplored. We performed sex-stratified meta-analysis after protein quantitative trait loci (pQTL) mapping using inflammatory biomarkers profiled using targeted proteomics (Olink inflammatory panel) of two population-based cohorts of Europeans.

RESULTS

Even though, around 67% of the pQTLs demonstrated shared effect between sexes, colocalization analysis identified two loci in the males (LINC01135 and ITGAV) and three loci (CNOT10, SRD5A2, and LILRB5) in the females with evidence of sex-dependent modulation by pQTL variants. Furthermore, we identified pathways with relevant functions in the sex-biased pQTL variants. We also showed through cross-validation that the sex-specific pQTLs are linked with sex-specific phenotypic traits.

CONCLUSION

Our study demonstrates the relevance of genetic sex-stratified analysis in the context of genetic dissection of protein abundances among individuals and reveals that, sex-specific pQTLs might mediate sex-linked phenotypes. Identification of sex-specific pQTLs associated with sex-biased diseases can help realize the promise of individualized treatment.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Structural brain morphometry differences and similarities between young patients with Crohn's disease in remission and healthy young and old controls

Introduction

Crohn's disease (CD), one of the main phenotypes of inflammatory bowel disease (IBD), can affect any part of the gastrointestinal tract. It can impact the function of gastrointestinal secretions, as well as increasing the intestinal permeability leading to an aberrant immunological response and subsequent intestinal inflammation. Studies have reported anatomical and functional brain changes in Crohn's Disease patients (CDs), possibly due to increased inflammatory markers and microglial cells that play key roles in communicating between the brain, gut, and systemic immune system. To date, no studies have demonstrated similarities between morphological brain changes seen in IBD and brain morphometry observed in older healthy controls..

Methods

For the present study, twelve young CDs in remission (M = 26.08 years, SD = 4.9 years, 7 male) were recruited from an IBD Clinic. Data from 12 young age-matched healthy controls (HCs) (24.5 years, SD = 3.6 years, 8 male) and 12 older HCs (59 years, SD = 8 years, 8 male), previously collected for a different study under a similar MR protocol, were analyzed as controls. T1 weighted images and structural image processing techniques were used to extract surface-based brain measures, to test our hypothesis that young CDs have different brain surface morphometry than their age-matched young HCs and furthermore, appear more similar to older HCs. The phonemic verbal fluency (VF) task (the Controlled Oral Word Association Test, COWAT) (Benton, 1976) was administered to test verbal cognitive ability and executive control.

Results/Discussion

On the whole, CDs had more brain regions with differences in brain morphometry measures when compared to the young HCs as compared to the old HCs, suggesting that CD has an effect on the brain that makes it appear more similar to old HCs. Additionally, our study demonstrates this atypical brain morphometry is associated with function on a cognitive task. These results suggest that even younger CDs may be showing some evidence of structural brain changes that demonstrate increased resemblance to older HC brains rather than their similarly aged healthy counterparts.

Multi-omics analysis of innate and adaptive responses to BCG vaccination reveals epigenetic cell states that predict trained immunity

Immune responses are tightly regulated yet highly variable between individuals. To investigate human population variation of trained immunity, we immunized healthy individuals with Bacillus Calmette-Guérin (BCG). This live-attenuated vaccine induces not only an adaptive immune response against tuberculosis but also triggers innate immune activation and memory that are indicative of trained immunity. We established personal immune profiles and chromatin accessibility maps over a 90-day time course of BCG vaccination in 323 individuals. Our analysis uncovered genetic and epigenetic predictors of baseline immunity and immune response. BCG vaccination enhanced the innate immune response specifically in individuals with a dormant immune state at baseline, rather than providing a general boost of innate immunity. This study advances our understanding of BCG's heterologous immune-stimulatory effects and trained immunity in humans. Furthermore, it highlights the value of epigenetic cell states for connecting immune function with genotype and the environment.

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.

Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

Background

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

Methods

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,162 cancer cases and up to 824,556 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 x 10-8) cis-acting SNPs (i.e. in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") < 0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) > 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes.

Results

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR 1.19, 95% CI 1.10-1.29, q-value=0.033, PPH4=84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR 1.42, 95% CI 1.20-1.69, q-value=0.055, PPH4=73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR 0.66, 95% CI 0.53-0.81, q-value=0.067, PPH4=81.8%), macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR 1.14, 95% CI 1.05-1.23, q-value=0.072, PPH4=76.1%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR 0.92, 95% CI 0.88-0.97, q-value=0.15), PPH4=85.6%). For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥ 0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

Conclusion

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 5 circulating inflammatory markers in risk of 5 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

Fibroblast Growth Factor-23 and Risk of Cardiovascular Diseases: A Mendelian Randomization Study

BACKGROUND

Fibroblast growth factor-23 (FGF-23) is associated with a range of cardiovascular and noncardiovascular diseases in conventional epidemiological studies, but substantial residual confounding may exist. Mendelian randomization approaches can help control for such confounding.

METHODS

SCALLOP Consortium data of 19,195 participants were used to generate an FGF-23 genetic score. Data from 337,448 UK Biobank participants were used to estimate associations between higher genetically predicted FGF-23 concentration and the odds of any atherosclerotic cardiovascular disease (n=26,266 events), nonatherosclerotic cardiovascular disease (n=12,652), and noncardiovascular diseases previously linked to FGF-23. Measurements of carotid intima-media thickness and left ventricular mass were available in a subset. Associations with cardiovascular outcomes were also tested in three large case-control consortia: CARDIOGRAMplusC4D (coronary artery disease, n=181,249 cases), MEGASTROKE (stroke, n=34,217), and HERMES (heart failure, n=47,309).

RESULTS

We identified 34 independent variants for circulating FGF-23, which formed a validated genetic score. There were no associations between genetically predicted FGF-23 and any of the cardiovascular or noncardiovascular outcomes. In UK Biobank, the odds ratio (OR) for any atherosclerotic cardiovascular disease per 1-SD higher genetically predicted logFGF-23 was 1.03 (95% confidence interval [95% CI], 0.98 to 1.08), and for any nonatherosclerotic cardiovascular disease, it was 1.01 (95% CI, 0.94 to 1.09). The ORs in the case-control consortia were 1.00 (95% CI, 0.97 to 1.03) for coronary artery disease, 1.01 (95% CI, 0.95 to 1.07) for stroke, and 1.00 (95% CI, 0.95 to 1.05) for heart failure. In those with imaging, logFGF-23 was not associated with carotid or cardiac abnormalities.

CONCLUSIONS

Genetically predicted FGF-23 levels are not associated with atherosclerotic and nonatherosclerotic cardiovascular diseases, suggesting no important causal link.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/CJASN/2023_01_10_CJN05080422.mp3.

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.

Integrated methylome and phenome study of the circulating proteome reveals markers pertinent to brain health

Characterising associations between the methylome, proteome and phenome may provide insight into biological pathways governing brain health. Here, we report an integrated DNA methylation and phenotypic study of the circulating proteome in relation to brain health. Methylome-wide association studies of 4058 plasma proteins are performed (N = 774), identifying 2928 CpG-protein associations after adjustment for multiple testing. These are independent of known genetic protein quantitative trait loci (pQTLs) and common lifestyle effects. Phenome-wide association studies of each protein are then performed in relation to 15 neurological traits (N = 1,065), identifying 405 associations between the levels of 191 proteins and cognitive scores, brain imaging measures or APOE e4 status. We uncover 35 previously unreported DNA methylation signatures for 17 protein markers of brain health. The epigenetic and proteomic markers we identify are pertinent to understanding and stratifying brain health.

Genome- and epigenome-wide studies of plasma protein biomarkers for Alzheimer's disease implicate TBCA and TREM2 in disease risk

Introduction

The levels of many blood proteins are associated with Alzheimer's disease (AD) or its pathological hallmarks. Elucidating the molecular factors that control circulating levels of these proteins may help to identify proteins associated with disease risk mechanisms.

Methods

Genome-wide and epigenome-wide studies (nindividuals ≤1064) were performed on plasma levels of 282 AD-associated proteins, identified by a structured literature review. Bayesian penalized regression estimated contributions of genetic and epigenetic variation toward inter-individual differences in plasma protein levels. Mendelian randomization (MR) and co-localization tested associations between proteins and disease-related phenotypes.

Results

Sixty-four independent genetic and 26 epigenetic loci were associated with 45 proteins. Novel findings included an association between plasma triggering receptor expressed on myeloid cells 2 (TREM2) levels and a polymorphism and cytosine-phosphate-guanine (CpG) site within the MS4A4A locus. Higher plasma tubulin-specific chaperone A (TBCA) and TREM2 levels were significantly associated with lower AD risk.

Discussion

Our data inform the regulation of biomarker levels and their relationships with AD.

DNA methylome-wide association study of genetic risk for depression implicates antigen processing and immune responses

BACKGROUND

Depression is a disabling and highly prevalent condition where genetic and epigenetic, such as DNA methylation (DNAm), differences contribute to disease risk. DNA methylation is influenced by genetic variation but the association between polygenic risk of depression and DNA methylation is unknown.

METHODS

We investigated the association between polygenic risk scores (PRS) for depression and DNAm by conducting a methylome-wide association study (MWAS) in Generation Scotland (N = 8898, mean age = 49.8 years) with replication in the Lothian Birth Cohorts of 1921 and 1936 and adults in the Avon Longitudinal Study of Parents and Children (ALSPAC) (Ncombined = 2049, mean age = 79.1, 69.6 and 47.2 years, respectively). We also conducted a replication MWAS in the ALSPAC children (N = 423, mean age = 17.1 years). Gene ontology analysis was conducted for the cytosine-guanine dinucleotide (CpG) probes significantly associated with depression PRS, followed by Mendelian randomisation (MR) analysis to infer the causal relationship between depression and DNAm.

RESULTS

Widespread associations (NCpG = 71, pBonferroni < 0.05, p < 6.3 × 10-8) were found between PRS constructed using genetic risk variants for depression and DNAm in CpG probes that localised to genes involved in immune responses and neural development. The effect sizes for the significant associations were highly correlated between the discovery and replication samples in adults (r = 0.79) and in adolescents (r = 0.82). Gene Ontology analysis showed that significant CpG probes are enriched in immunological processes in the human leukocyte antigen system. Additional MWAS was conducted for each lead genetic risk variant. Over 47.9% of the independent genetic risk variants included in the PRS showed associations with DNAm in CpG probes located in both the same (cis) and distal (trans) locations to the genetic loci (pBonferroni < 0.045). Subsequent MR analysis showed that there are a greater number of causal effects found from DNAm to depression than vice versa (DNAm to depression: pFDR ranged from 0.024 to 7.45 × 10-30; depression to DNAm: pFDR ranged from 0.028 to 0.003).

CONCLUSIONS

PRS for depression, especially those constructed from genome-wide significant genetic risk variants, showed methylome-wide differences associated with immune responses. Findings from MR analysis provided evidence for causal effect of DNAm to depression.

The genetic and epigenetic profile of serum S100β in the Lothian Birth Cohort 1936 and its relationship to Alzheimer's disease

Background: Circulating S100 calcium-binding protein (S100β) is a marker of brain inflammation that has been associated with a range of neurological conditions. To provide insight into the molecular regulation of S100β and its potential causal associations with Alzheimer's disease, we carried out genome- and epigenome-wide association studies (GWAS/EWAS) of serum S100β levels in older adults and performed Mendelian randomisation with Alzheimer's disease. Methods: GWAS (N=769, mean age 72.5 years, sd = 0.7) and EWAS (N=722, mean age 72.5 years, sd = 0.7) of S100β levels were performed in participants from the Lothian Birth Cohort 1936. Conditional and joint analysis (COJO) was used to identify independent loci. Expression quantitative trait locus (eQTL) analyses were performed for lead loci that had genome-wide significant associations with S100β. Bidirectional, two-sample Mendelian randomisation was used to test for causal associations between S100β and Alzheimer's disease. Colocalisation between S100β and Alzheimer's disease GWAS loci was also examined. Results: We identified 154 SNPs from chromosome 21 that associated (P<5x10 -8) with S100β protein levels. The lead variant was located in the S100β gene (rs8128872, P=5.0x10 -17). We found evidence that two independent causal variants existed for both transcription of S100β and S100β protein levels in our eQTL analyses . No CpG sites were associated with S100β levels at the epigenome-wide significant level (P<3.6x10 -8); the lead probe was cg06833709 (P=5.8x10 -6), which mapped to the LGI1 gene. There was no evidence of a causal association between S100β levels and Alzheimer's disease or vice versa and no evidence for colocalisation between S100β and Alzheimer's disease loci. Conclusions: These data provide insight into the molecular regulators of S100β levels. This context may aid in understanding the role of S100β in brain inflammation and neurological disease.

Blood-based epigenome-wide analyses of cognitive abilities

BACKGROUND

Blood-based markers of cognitive functioning might provide an accessible way to track neurodegeneration years prior to clinical manifestation of cognitive impairment and dementia.

RESULTS

Using blood-based epigenome-wide analyses of general cognitive function, we show that individual differences in DNA methylation (DNAm) explain 35.0% of the variance in general cognitive function (g). A DNAm predictor explains ~4% of the variance, independently of a polygenic score, in two external cohorts. It also associates with circulating levels of neurology- and inflammation-related proteins, global brain imaging metrics, and regional cortical volumes.

CONCLUSIONS

As sample sizes increase, the ability to assess cognitive function from DNAm data may be informative in settings where cognitive testing is unreliable or unavailable.

Epigenetic scores for the circulating proteome as tools for disease prediction

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

The genetic and epigenetic profile of serum S100β in the Lothian Birth Cohort 1936 and its relationship to Alzheimer's disease

Background: Circulating S100 calcium-binding protein (S100β) is a marker of brain inflammation that has been associated with a range of neurological conditions. To provide insight into the molecular regulation of S100β and its potential causal associations with Alzheimer's disease, we carried out genome- and epigenome-wide association studies (GWAS/EWAS) of serum S100β levels in older adults and performed Mendelian randomisation with Alzheimer's disease. Methods: GWAS (N=769, mean age 72.5 years, sd = 0.7) and EWAS (N=722, mean age 72.5 years, sd = 0.7) of S100β levels were performed in participants from the Lothian Birth Cohort 1936. Conditional and joint analysis (COJO) was used to identify independent loci. Expression quantitative trait locus (eQTL) analyses were performed for lead loci that had genome-wide significant associations with S100β. Bidirectional, two-sample Mendelian randomisation was used to test for causal associations between S100β and Alzheimer's disease. Colocalisation between S100β and Alzheimer's disease GWAS loci was also examined. Results: We identified 154 SNPs from chromosome 21 that associated (P<5x10 -8) with S100β protein levels. The lead variant was located in the S100β gene (rs8128872, P=5.0x10 -17). We found evidence that two independent causal variants existed for both transcription of S100β and S100β protein levels in our eQTL analyses . No CpG sites were associated with S100β levels at the epigenome-wide significant level (P<3.6x10 -8); the lead probe was cg06833709 (P=5.8x10 -6), which mapped to the LGI1 gene. There was no evidence of a causal association between S100β levels and Alzheimer's disease or vice versa and no evidence for colocalisation between S100β and Alzheimer's disease loci. Conclusions: These data provide insight into the molecular regulators of S100β levels. This context may aid in understanding the role of S100β in brain inflammation and neurological disease.

The 3p21.31 genetic locus promotes progression to type 1 diabetes through the CCR2/CCL2 pathway

Multiple cross-sectional and longitudinal studies have shown that serum levels of the chemokine ligand 2 (CCL-2) are associated with type 1 diabetes (T1D), although the direction of effect differs. We assessed CCL-2 serum levels in a longitudinal cohort to clarify this association, combined with genetic data to elucidate the regulatory role of CCL-2 in T1D pathogenesis. The Diabetes Autoimmunity Study in the Young (DAISY) followed 310 subjects with high risk of developing T1D. Of these, 42 became persistently seropositive for islet autoantibodies but did not develop T1D (non-progressors); 48 did develop T1D (progressors). CCL-2 serum levels among the three study groups were compared using linear mixed models adjusting for age, sex, HLA genotype, and family history of T1D. Summary statistics were obtained from the CCL-2 protein quantitative trait loci (pQTL) and CCR2 expression QTL (eQTL) studies. The T1D fine mapping association data were provided by the Type 1 Diabetes Genetics Consortium (T1DGC). Serum CCL-2 levels were significantly lower in both progressors (p = 0.004) and non-progressors (p = 0.005), compared to controls. Two SNPs (rs1799988 and rs746492) in the 3p21.31 genetic locus, which includes the CCL-2 receptor, CCR2, were associated with increased CCR2 expression (p = 8.2e-5 and 5.2e-5, respectively), decreased CCL-2 serum level (p = 2.41e-9 and 6.21e-9, respectively), and increased risk of T1D (p = 7.9e-5 and 7.9e-5, respectively). The 3p21.31 genetic region is associated with developing T1D through regulatory control of the CCR2/CCL2 immune pathway.

Connecting the epigenome, metabolome and proteome for a deeper understanding of disease

Epigenome-wide association studies (EWAS) identify genes that are dysregulated by the studied clinical endpoints, thereby indicating potential new diagnostic biomarkers, drug targets and therapy options. Combining EWAS with deep molecular phenotyping, such as approaches enabled by metabolomics and proteomics, allows further probing of the underlying disease-associated pathways. For instance, methylation of the TXNIP gene is associated robustly with prevalent type 2 diabetes and further with metabolites that are short-term markers of glycaemic control. These associations reflect TXNIP's function as a glucose uptake regulator by interaction with the major glucose transporter GLUT1 and suggest that TXNIP methylation can be used as a read-out for the organism's exposure to glucose stress. Another case is the association between DNA methylation of the AHRR and F2RL3 genes with smoking and a protein that is involved in the reprogramming of the bronchial epithelium. These examples show that associations between DNA methylation and intermediate molecular traits can open new windows into how the body copes with physiological challenges. This knowledge, if carefully interpreted, may indicate novel therapy options and, together with monitoring of the methylation state of specific methylation sites, may in the future allow the early diagnosis of impending disease. It is essential for medical practitioners to recognize the potential that this field holds in translating basic research findings to clinical practice. In this review, we present recent advances in the field of EWAS with metabolomics and proteomics and discuss both the potential and the challenges of translating epigenetic associations, with deep molecular phenotypes, to biomedical applications.