A generalized linear mixed model association tool for biobank-scale data

Clinical significance of glomerular IgM deposit in IgA nephropathy: a 5-year follow-up study

The significance of glomerular IgM deposit intensity in IgA Nephropathy (IgAN) remained ambiguous and requires further research. Patients with biopsy-proven IgAN in our hospital from January 2018 to May 2023 were recruited into this retrospective single-center study. Patients who presented with positive IgM deposit were included in IgM + cohort while patients with negative IgM deposit were included in IgM- cohort. Of the IgM+, patients whose IF intensity of IgM deposits exceeded 1+ formed IgM-H cohort while patients whose IF intensity of IgM deposits was equal to 1+ consisted IgM-L cohort. Pairwise comparisons were performed among these cohorts to determine clinical disparities, following the propensity score matching process. Among 982 IgAN patients, 539 patients presented with positive IgM deposit. The Kaplan-Meier analysis showed that the IgM deposit did not contribute adversely to the outcomes (eGFR decreased from the baseline ≥ 50% continuously or reached end-stage renal disease). However, the Cox regression analysis showed that increased intensity of IgM deposit was an independent risk factor (p = 0.03) in IgM+. The IgM-H exhibited more pronounced segmental glomerulosclerosis (p = 0.02) than the IgM-L, which may also be associated more directly with higher urine protein levels (p = 0.02). Moreover, our generalized linear mixed model demonstrated a remarkably higher urine albumin/creatinine ratio (p < 0.01) and serum creatinine (p = 0.04) levels as well as lower serum albumin (p < 0.01) level in IgM-H persistently during the 5-year follow-up. This study concluded that increased intensity of glomerular IgM deposits may contribute adversely to clinicopathologic presentation and outcome in those IgM + patients.

Causal relationships between neuropsychiatric disorders and nonalcoholic fatty liver disease: A bidirectional Mendelian randomization study

BACKGROUND

Increasing evidences suggest that nonalcoholic fatty liver disease (NAFLD) is associated with neuropsychiatric disorders. Nevertheless, whether there were causal associations between them remained vague. A causal association between neuropsychiatric disorders and NAFLD was investigated in this study.

METHODS

We assessed the published genome-wide association study summary statistics for NAFLD, seven mental disorder-related diseases and six central nervous system dysfunction-related diseases. The causal relationships were first assessed using two-sample and multivariable Mendelian randomization (MR). Then, sensitivity analyses were performed, followed by a reverse MR analysis to determine whether reverse causality is possible. Finally, we performed replication analyses and combined the findings from the above studies.

RESULTS

Our meta-analysis results showed NAFLD significantly increased the risk of anxiety disorders (OR = 1.016, 95% CI = 1.010-1.021, P value < 0.0001). In addition, major depressive disorder was the potential risk factor for NAFLD (OR = 1.233, 95% CI = 1.063-1.430, P value = 0.006). Multivariable MR analysis showed that the causal effect of major depressive disorder on NAFLD remained significant after considering body mass index, but the association disappeared after adjusting for the effect of waist circumference. Furthermore, other neuropsychiatric disorders and NAFLD were not found to be causally related.

CONCLUSIONS

These results implied causal relationships of NAFLD with anxiety disorders and Major Depressive Disorder. This study highlighted the need to recognize and understand the connection between neuropsychiatric disorders and NAFLD to prevent the development of related diseases.

Minor Genetic Overlap Among Rheumatoid Arthritis, Myocardial Infarction, and Myocardial Infarction Risk Determinants

OBJECTIVE

The aim of this study was to investigate whether a shared genetic susceptibility exists between individuals with rheumatoid arthritis (RA) and individuals with myocardial infarction (MI)-including major MI risk factors-and to quantify the degree of any such overlap.

METHODS

Genome-wide association study (GWAS) data for individuals with RA were constructed from a sample of 26,637 Swedish patients with RA and controls without RA. For patients with MI, GWAS data were obtained from a previously published meta-analysis. Genome-wide genetic correlation was estimated via linkage disequilibrium score regression. LAVA was employed to estimate local genetic correlations in ~2,500 nonoverlapping loci, including the major histocompatibility complex. The controls without RA were used for reference panel data. We also assessed stratified estimates of both genome-wide and local genetic correlation based on subsamples of individuals with seropositive RA and those with seronegative RA. Furthermore, genome-wide genetic correlation was estimated between RA and selected cardiovascular risk factors to elucidate pleiotropic relationships.

RESULTS

Following quality control, our GWAS of patients with RA consisted of 25,826 individuas. Genome-wide genetic correlation between patients with RA and MI was estimated to 0.13 (95% confidence interval -0.03 to 0.29). Six regions exhibited significant local genetic correlation, though none harbored any known risk single-nucleotide polymorphisms for either of the two traits. Estimates were similar in both individuals with seropositive RA and those with seronegative RA. No statistically significant genetic correlations were observed between RA risk factors and any of the MI risk factors.

CONCLUSION

Our findings indicate that genetic overlap between patients with RA and MI is minor. Furthermore, genetic overlap between RA and MI risk factors seem unlikely to provide a major contribution to the increased risk of MI observed in patients with RA.

Causal relationships between gut microbiota, plasma metabolites, and HIV infection: insights from Mendelian randomization and mediation analysis

OBJECTIVE

Gut dysbiosis and metabolic abnormalities have been implicated in HIV infection. However, the exact causal relationships among the gut microbiota, metabolites, and HIV infection remain poorly understood. Our study involving Mendelian randomization (MR) and mediation analysis aims to unveil these causalities.

METHODS

Genetic instrumental variables for the gut microbiota were retrieved from MiBioGen consortium (n = 18,340). Metabolism-related genetic variants were sourced from the CLSA cohort (n = 8299). GWAS summary statistics for symptomatic HIV infection were derived from the FinnGen study (n = 309,154), and the UK Biobank (n = 208,808). We performed the bidirectional two-sample MR to assess causalities with the inverse-variance weighted (IVW) method as the primary analysis. Moreover, we executed a mediation analysis using two-step MR methods.

RESULTS

Compared to the causal effects of HIV infection on gut microbiota (or metabolites), those of gut microbiota (or plasma metabolites) on the risk of HIV infection were more substantial. Phylum Proteobacteria (OR: 2.114, 95% CI 1.042-4.288, P = 0.038), and genus Ruminococcaceae UCG013 (OR: 2.127, 95% CI 1.080-4.191, P = 0.029) exhibited an adverse causal effect on HIV infection, whereas genus Clostridium sensu stricto 1(OR: 0.491, 95% CI 0.252-0.956, P = 0.036) and family Erysipelotrichaceae (OR: 0.399, 95% CI 0.193-0.827, P = 0.013) acted as significant protective factors for HIV. The salicyluric glucuronide level (OR = 2.233, 95% CI 1.120-4.453, P = 0.023) exhibited a considerably adverse causal effect on HIV infection. Conversely, the salicylate-to-citrate ratio (OR: 0.417, 95% CI 0.253-0.688, P = 0.001) was identified as a protective factor for HIV. We identified only one bidirectional causality between 1-palmitoyl-GPI and HIV infection. Mechanistically, genus Haemophilus mediated the causal effects of three phospholipids on HIV infection risk: 1-palmitoyl-GPI (mediation proportion = 33.7%, P = 0.018), 1-palmitoyl-2-arachidonoyl-GPI (mediation proportion = 18.3%, P = 0.019), and 1-linoleoyl-2-linolenoyl-GPC (mediation proportion = 20.3%, P = 0.0216). Additionally, 5-Dodecenoylcarnitine (C12:1) mediated the causal effect of genus Sellimonas on the risk of HIV infection (mediation proportion = 13.7%, P = 0.0348).

CONCLUSION

Our study revealed that gut microbiota and metabolites causally influence HIV infection risk more substantially than the reverse. We identified the bidirectional causality between 1-palmitoyl-GPI (16:0) and HIV infection, and elucidated four mediation relationships. These findings provide genetic insights into prediction, prevention, and personalized medicine of HIV infection.

A comprehensive analysis of the health effects associated with smoking in the largest population using UK Biobank genotypic and phenotypic data

Background

Smoking is a widespread behavior, while the relationship between smoking and various diseases remains a topic of debate.

Objective

We conducted analysis to further examine the identified associations and assess potential causal relationships.

Methods

We utilized seven single nucleotide polymorphisms (SNPs) known to be linked to smoking extracting genotype data from the UK Biobank, a large-scale biomedical repository encompassing comprehensive health-related and genetic information of European descent. Phenome-wide association study (PheWAS) analysis was conducted to map the association of genetically predicted smoking status with 1,549 phenotypes. The associations identified in the PheWAS were then meticulously examined through two-sample Mendelian randomization (MR) analysis, utilizing data from the UK Biobank (n = 487,365) and the Sequencing Consortium of Alcohol and Nicotine Use (GSCAN) (n = 337,334). This approach allowed us to comprehensively characterize the links between smoking and disease patterns.

Results

The PheWAS analysis produced 34 phenotypes that demonstrated significant associations with smoking (P = 0.05/1460). Importantly, sickle cell anemia and type 2 diabetes exhibited the most significant SNPs (both 85.71% significant SNPs). Furthermore, the MR analyses provided compelling evidence supporting causal associations between smoking and the risk of following diseases: obstructive chronic bronchitis (IVW: Beta = 0.48, 95% confidence interval (CI) 0.36-0.61, P = 1.62×10-13), cancer of the bronchus (IVW: Beta = 0.92, 95% CI 0.68-1.17, P = 2.02×10-13), peripheral vascular disease (IVW: Beta = 1.09, 95% CI 0.71-1.46, P = 1.63×10-8), emphysema (IVW: Beta = 1.63, 95% CI 0.90-2.36, P = 1.29×10-5), pneumococcal pneumonia (IVW: Beta = 0.30, 95% CI 0.11-0.49, P = 1.60×10-3), chronic airway obstruction (IVW: Beta = 0.83, 95% CI 0.30-1.36, P = 2.00×10-3) and type 2 diabetes (IVW: Beta = 0.53, 95% CI 0.16-0.90, P = 5.08×10-3).

Conclusion

This study affirms causal relationships between smoking and obstructive chronic bronchitis, cancer of the bronchus, peripheral vascular disease, emphysema, pneumococcal pneumonia, chronic airway obstruction, type 2 diabetes, in the European population. These findings highlight the broad health impacts of smoking and support smoking cessation efforts.

Genome-wide association study based on clustering by obesity-related variables uncovers a genetic architecture of obesity in the Japanese and the UK populations

Whether all obesity-related variants contribute to the onset of obesity or one or a few variants cause obesity in genetically heterogeneous populations remains obscure. Here, we investigated the genetic architecture of obesity by clustering the Japanese and British populations with obesity using obesity-related factors. In Step-1, we conducted a genome-wide association study (GWAS) with body mass index (BMI) as the outcome for eligible participants. In Step-2, we assigned participants with obesity (BMI ≥25 kg/m2) to five clusters based on obesity-related factors. Subsequently, participants from each cluster and those with a BMI <25 kg/m2 were combined. A GWAS was conducted for each cluster. Several previously identified obesity-related genes were verified in Step-1. Of the genes detected in Step-1, unique obesity-related genes were detected separately for each cluster in Step-2. Our novel findings suggest that a smaller sample size with increased homogeneity may provide insights into the genetic architecture of obesity.

Identifying novel proteins underlying bipolar disorder via integrating pQTLs of the plasma, CSF, and brain with GWAS summary data

Bipolar disorder (BD) presents a significant challenge due to its chronic and relapsing nature, with its underlying pathogenesis remaining elusive. This study employs Mendelian randomization (MR), a widely recognized genetic approach, to unveil intricate causal associations between proteins and BD, leveraging protein quantitative trait loci (pQTL) as key exposures. We integrate pQTL data from brain, cerebrospinal fluid (CSF), and plasma with genome-wide association study (GWAS) findings of BD within a comprehensive systems analysis framework. Our analyses, including two-sample MR, Steiger filtering, and Bayesian colocalization, reveal noteworthy associations. Elevated levels of AGRP, FRZB, and IL36A in CSF exhibit significant associations with increased BD_ALL risk, while heightened levels of CTSF and LRP8 in CSF, and FLRT3 in plasma, correlate with decreased BD_ALL risk. Specifically for Bipolar I disorder (BD_I), increased CSF AGRP levels are significantly linked to heightened BD_I risk, whereas elevated CSF levels of CTSF and LRP8, and plasma FLRT3, are associated with reduced BD_I risk. Notably, genes linked to BD-related proteins demonstrate substantial enrichment in functional pathways such as "antigen processing and presentation," "metabolic regulation," and "regulation of myeloid cell differentiation." In conclusion, our findings provide beneficial evidence to support the potential causal relationship between IL36A, AGRP, FRZB, LRP8 in cerebrospinal fluid, and FLRT3 in plasma, and BD and BD_I, providing insights for future mechanistic studies and therapeutic development.

Biological aging mediates the association between periodontitis and cardiovascular disease: results from a national population study and Mendelian randomization analysis

BACKGROUND

The relationship between periodontitis and cardiovascular disease (CVD) has been extensively studied, but the role of biological aging in this relationship remains poorly understood. This study is dedicated to investigating the effect of periodontitis on the incidence of CVD and to elucidating the potential mediating role of biological aging. Furthermore, this study will seek to elucidate the causal association between periodontitis, CVD, and biological aging.

METHODS

We included 3269 participants from the National Health and Nutrition Examination Survey (2009-2014) with diagnostic information on periodontitis and composite CVD events. Biological aging was evaluated by utilizing both the Klemera-Doubal method's calculated biological age (KDMAge) and phenotypic age (PhenoAge). Logistic regression, restricted cubic spline (RCS) analysis, and subgroup analysis were used for data analysis. Mediation analysis was employed to explore the mediating role of biological aging. Subsequently, Mendelian randomization (MR) analyses were performed using genome-wide association study databases to explore potential causal relationships between periodontitis, CVD, and biological aging.

RESULTS

Periodontitis was associated with a higher risk of CVD. Participants with periodontitis were found to have increased levels of biological aging, and elevated levels of biological aging were associated with increased CVD risk. Mediation analyses showed a partial mediating effect of biological aging (PhenoAge: 44.6%; KDMAge: 22.9%) between periodontitis and CVD risk. MR analysis showed that periodontitis played a causal role in increasing the risk of small vessel stroke, while myocardial infarction was found to increase the risk of periodontitis. In addition, reverse MR analysis showed that phenotypic aging can increase the risk of periodontitis, and there is a two-way causal relationship between CVD and biological aging.

CONCLUSIONS

Periodontitis is associated with an increased CVD risk, partially mediated by biological aging, with a complex causal interrelationship. Targeted interventions for periodontal health may slow the biological aging processes and reduce CVD risk.

Causal relationship between genetically predicted antibody-Mediated Immune Responses and female infertility

OBJECTIVE

This study aims to explore the causal relationship between Antibody-Mediated Immune Response (AIR) and female infertility using Mendelian randomization analysis, providing new targets and insights for treating female infertility.

METHODS

Exposure data from a genome-wide association study (GWAS) involving 9724 participants provided information on 46 types of AIR. Outcome data were sourced from a large-scale meta-analysis of female infertility, which included 899 cases and 246,641 controls. Five different Mendelian randomization analysis methods were used to assess the genetic causal relationship between AIR and female infertility, with the inverse variance-weighted method as the primary approach. The robustness of the study results was evaluated through heterogeneity analysis, pleiotropy tests, and leave-one-out validation. Replication and meta-analysis were conducted using FinnGen data to further identify AIR.

RESULTS

Through rigorous genetic variant screening, IVW, sensitivity analysis, replication, and meta-analysis, the results indicate that only seropositivity of anti-polyomavirus 2 IgG shows significant evidence of association with female infertility. The reverse MR analysis did not reveal a genetic causal association between female infertility and this antibody-influencing relationship. Sensitivity analysis showed no evidence of pleiotropy or heterogeneity. Replication and meta-analysis using FinnGen data further supported that only anti-polyomavirus 2 IgG seropositivity has a significant association with female infertility.

CONCLUSION

Using Mendelian randomization analysis in two independent samples, this study presents evidence of a genetic association between anti-polyomavirus 2 IgG seropositivity and female infertility. These findings carry significant implications for the prevention and treatment of female infertility in clinical practice.

Causal relationship between gut microbiota and diabetic complications: a two-sample Mendelian randomization study

BACKGROUND

Imbalances in gut microbiota (GM) have been proposed as a potential contributing factor to diabetic complications; however, the causal relationship remains incompletely understood.

METHODS

Summary statistics were obtained from genome-wide association studies (GWAS) of 196 gut microbial taxa, including 9 phyla, 16 classes, 20 orders, 32 families, and 119 genera. These data were then analyzed using mediation Mendelian randomization (MR) analyses to explore the potential mediating effect of diabetes complications risk factors on the relationship between gut microbiota and specific diabetic complications such as diabetic kidney disease (DKD), ketoacidosis, and diabetic retinopathy (DR).

RESULTS

In our Mendelian analysis, we observed negative associations between Bifidobacterial order and Actinomycete phylum with DKD in type 1 diabetes (T1D) as well as early DKD in T1D. Conversely, these taxa showed positive associations with ketoacidosis in type 2 diabetes (T2D). In reverse Mendelian analysis, we found that DR in both T1D and T2D as well as ketoacidosis in T2D affected the abundance of Eubacterium fissicaten genus and LachnospiraceaeUCG010 family within the gut microbiota.

CONCLUSIONS

Our findings provide compelling evidence for causal relationships between specific GM taxa and various diabetes complications. These insights contribute valuable knowledge for developing treatments targeting diabetes-related complications.

Dissecting the mediating role of inflammatory factors in the interaction between metabolites and sepsis: insights from bidirectional Mendelian randomization

Sepsis, a life-threatening condition, involves complex interactions among metabolic alterations, inflammatory mediators, and host responses. This study utilized a bidirectional Mendelian randomization approach to investigate the causal relationships between 1400 metabolites and sepsis, and the mediating role of inflammatory factors. We identified 36 metabolites significantly associated with sepsis (p < 0.05), with AXIN1, FGF-19, FGF-23, IL-4, and OSM showing an inverse association, suggesting a protective role, while IL-2 exhibited a positive correlation, indicating a potential risk factor. Among these metabolites, Piperine and 9-Hydroxystearate demonstrated particularly interesting protective effects against sepsis. Piperine's protective effect was mediated through its interaction with AXIN1, contributing to a 16.296% reduction in sepsis risk. This suggests a potential pathway where Piperine influences sepsis outcomes by modulating AXIN1 levels. 9-Hydroxystearate also exhibited a protective role against sepsis, mediated through its positive association with FGF-19 and negative association with IL-2, contributing 9.436% and 12.565%, respectively, to its protective effect. Experimental validation confirmed significantly elevated IL-2 levels and reduced FGF-19, AXIN1, piperine, and 9-hydroxyoctadecanoic acid levels in sepsis patients compared to healthy controls. Piperine levels positively correlated with AXIN1, while 9-hydroxyoctadecanoic acid levels negatively correlated with IL-2 and positively correlated with FGF-19, supporting the Mendelian randomization findings. Our findings provide insights into the molecular mechanisms of sepsis, highlighting the unique roles and contributions of specific metabolites and their interactions with inflammatory mediators. This study enhances our understanding of sepsis pathophysiology and opens avenues for targeted therapeutic interventions and biomarker development for sepsis management. However, further research is essential to validate these pathways across diverse populations and fully explore the roles of these metabolites in sepsis.

Causal relationships between rheumatoid arthritis and neurodegenerative diseases: a two-sample univariable and multivariable Mendelian randomization study

Background

Observational research has highlighted a potential relationship between rheumatoid arthritis (RA) and neurodegenerative diseases (NDs). However, the confirmation of a causal connection is impeded by the inherent limitations of such studies, including vulnerability to confounding factors and the possibility of reverse causality. This study employs a two-sample Mendelian randomization (MR) approach to assess the causal impact of RA on three NDs, including Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS).

Methods

We aggregated data from genome-wide association studies (GWASs) targeting RA or NDs within populations of European descent. Single nucleotide polymorphisms (SNPs) with robust associations to RA were identified as instrumental variables (IVs). To estimate the association between RA and AD, PD, and ALS, we utilized the inverse variance weighted (IVW) method in our univariable MR (UVMR) analysis. Validation of the IVW results ensued through supplementary analyses using MR-Egger and weighted median methods. The multivariable MR (MVMR) analysis was conducted, adjusting for body mass index (BMI), alcohol drinking, and type 2 diabetes mellitus (T2DM).

Results

The UVMR analysis, based on the IVW method, revealed a significantly positive causal association between RA and late-onset (LO) AD (OR [95% CI] = 1.084 [1.020-1.153]; p = 9.980 × 10-3), while suggesting a possible inverse relationship with PD (OR [95% CI] = 0.727 [0.563-0.938]; p = 0.014). Our study did not detect any causal connections between RA and early-onset (EO) AD, atypical or mixed (AM) AD, and ALS (all p > 0.05). The MVMR analysis results indicated that after adjusting for alcohol drinking, RA remains a risk factor for LOAD (OR [95% CI] = 1.094 [1.024-1.169]; p = 0.008). However, MVMR analysis revealed no causal connections between RA and PD after adjustments for BMI, alcohol drinking, or T2DM (all p > 0.05). Sensitivity analyses showed no evidence of heterogeneity and horizontal pleiotropy.

Conclusions

This research provides genetic evidence indicating that RA potentially causes an increased risk of developing LOAD and PD. Such a revelation underscores the importance for individuals suffering from RA to be vigilant about the potential emergence of LOAD and PD. Ongoing monitoring and prompt detection are essential for successfully managing and intervening in this possible risk.

A two sample mendelian randomization analysis investigates causal effects between gut microbiome and immune related Vasculitis

Observational data suggest a link between gut microbiota and immune-related vasculitis, but causality remains unclear. A bidirectional mendelian randomization study was conducted using public genome-wide data. The inverse-variance-weighted (IVW) method identified associations and addressed heterogeneity.Families Clostridiaceae 1 and Actinomycetaceae correlated positively with granulomatosis with polyangiitis risk, while classes Lentisphaeria and Melainabacteria, and families Lachnospiraceae and Streptococcaceae showed negative associations. Behçet's disease was positively associated with the risk of family Streptococcaceae abundance. And other several gut microbiota constituents were identified as potential risk factors for immune-related vasculitis. Furthermore, combining positive association results from the IVW analysis revealed numerous shared gut microbiota constituents associated with immune-related vasculitis. MR analysis demonstrated a causal association between the gut microbiota and immune-related vasculitis, offering valuable insights for subsequent mechanistic and clinical investigations into microbiota-mediated immune-related vasculitis.

Exploring genetic associations in systemic lupus erythematosus through Mendelian randomization: implications for novel biomarkers and therapeutic targets

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by a significant health burden. There is an essential need for novel biomarkers and therapeutic targets to improve diagnosis and management. Mendelian randomization (MR) was applied to explore causal links between SLE and various biomarkers like immune cells, metabolites, and inflammatory cytokines using multiple databases. Initially, biomarkers significantly associated with SLE were identified. Bidirectional MR helped clarify these relationships, and a two-step mediation MR examined their effects on SLE risk. Intersection analysis was used to identify biomarkers with consistent effects across datasets. Four biomarkers were identified as having significant associations with SLE risk: 1-palmitoyl-2-arachidonoyl-GPI levels [odds ratio (OR), 1.379; 95% confidence interval (CI), 1.180 to 1.613; FDR, 0.046], IL-17A levels (OR, 2.197; 95% CI, 1.412 to 3.418; FDR, 0.044), N-acetyl-aspartyl-glutamate (NAAG) levels (OR, 0.882; 95% CI, 0.831 to 0.936; FDR, 0.030), and ribitol levels (OR, 0.743; 95% CI, 0.644 to 0.857; FDR, 0.012). Bidirectional MR showed an inverse effect of NAAG on IL-17A levels (OR, 0.978; 95% CI, 0.962 to 0.994; p = 0.006). Mediation analysis indicated that NAAG influenced SLE risk both directly (beta =  - 0.108) and indirectly through IL-17A (beta =  - 0.018), highlighting the potential mediating role of IL-17A. After expanding the significance criteria to p < 0.05, intersection analysis across multiple datasets revealed 29 biomarkers with consistent beta directions, including 19 potential risk factors (beta > 0) and 10 protective factors (beta < 0) for SLE. This research has revealed significant genetic associations with SLE and demonstrated that IL-17A mediates the relationship between NAAG levels and SLE risk, highlighting potential new targets for personalized therapeutic interventions. Key Points • This study employs MR to identify significant genetic associations between various biomarkers and SLE, providing novel insights into potential biomarkers and therapeutic targets. • Four key biomarkers were identified as significantly associated with SLE risk: 1-palmitoyl-2-arachidonoyl-GPI, IL-17A, N-acetyl-aspartyl-glutamate (NAAG), and ribitol. • The findings suggest that NAAG levels have a protective effect against SLE, partly mediated through IL-17A, indicating a complex interplay between these biomarkers in the pathogenesis of SLE. • Intersectional analysis across multiple datasets revealed 29 biomarkers with consistent effects on SLE risk, highlighting new directions for future research and potential personalized therapeutic strategies.

MOCHA's advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human cohorts

Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is being increasingly used to study gene regulation. However, major analytical gaps limit its utility in studying gene regulatory programs in complex diseases. In response, MOCHA (Model-based single cell Open CHromatin Analysis) presents major advances over existing analysis tools, including: 1) improving identification of sample-specific open chromatin, 2) statistical modeling of technical drop-out with zero-inflated methods, 3) mitigation of false positives in single cell analysis, 4) identification of alternative transcription-starting-site regulation, and 5) modules for inferring temporal gene regulatory networks from longitudinal data. These advances, in addition to open chromatin analyses, provide a robust framework after quality control and cell labeling to study gene regulatory programs in human disease. We benchmark MOCHA with four state-of-the-art tools to demonstrate its advances. We also construct cross-sectional and longitudinal gene regulatory networks, identifying potential mechanisms of COVID-19 response. MOCHA provides researchers with a robust analytical tool for functional genomic inference from scATAC-seq data.

Mendelian randomization study reveals causal effects of specific gut microbiota on the risk of interstitial cystitis/bladder pain syndrome (IC/BPS)

Evidence from previous studies have demonstrated that gut microbiota are closely associated with occurrence of interstitial cystitis/bladder pain syndrome (IC/BPS), yet the causal link between the two is not well known. In this study, we performed a two-sample Mendelian randomization (MR) analysis to determine the possible causal association between gut microbiota with IC/BPS. Gut microbiota summary level data were derived from the genome-wide association study (GWAS) conducted by MiBioGen and the IC/BPS GWAS summary level data were obtained from the GWAS Catalog. Next, we performed an MR study to investigate the causal link between gut microbiota and IC/BPS. The primary method for causal analysis was the inverse variance weighted (IVW), and the MR results were validated through multiple sensitivity analyses. A positive association was found between IC/BPS and eight gut microbial taxa, including genus Bacteroides, genus Haemophilus, genus Veillonella, genus Coprococcus1, genus Butyricimonas, family Bacteroidaceae, family Christensenellaceae, and order Lactobacillales. Sensitivity analysis revealed lack of significant pleiotropy or heterogeneity in the obtained results. This MR analysis reveals that a causal association exists between some gut microbiota with IC/BPS. This finding may is expected to guide future research and development of IC/BPS preventions and treatments based on the bladder-gut axis. However, given the clinical complexity and diagnostic challenges of IC/BPS, along with the limitations of using large-scale GWAS summary data for analysis, our MR results require further validation through additional research.

Multisystem health comorbidity networks of metabolic dysfunction-associated steatotic liver disease

BACKGROUND

The global burden of metabolic dysfunction-associated steatotic liver disease (MASLD) is growing, but its subsequent health consequences have not been thoroughly examined.

METHODS

A phenome-wide association study was conducted to map the associations of MASLD with 948 unique clinical outcomes among 361,021 Europeans in the UK Biobank. Disease trajectory and comorbidity analyses were applied to visualize the sequential patterns of multiple comorbidities related to the occurrence of MASLD. The associations jointly verified by observational and polygenic phenome-wide analyses were further replicated by two-sample Mendelian randomization analysis using data from the FinnGen study and international consortia.

FINDINGS

The observational and polygenic phenome-wide association study revealed the associations of MASLD with 96 intrahepatic and extrahepatic diseases, including circulatory, metabolic, genitourinary, neurological, gastrointestinal, and hematologic diseases. Sequential patterns of MASLD-related extrahepatic comorbidities were primarily found in circulatory, metabolic, and inflammatory diseases. Mendelian randomization analyses supported the causal associations between MASLD and the risk of several intrahepatic disorders, metabolic diseases, cardio-cerebrovascular disease, and ascites but found no associations with neurological diseases.

CONCLUSIONS

This study elucidated multisystem comorbidities and health consequences of MASLD, contributing to the development of combination interventions targeting distinct pathways for health promotion among patients with MASLD.

FUNDING

X.L. was funded by the Natural Science Fund for Distinguished Young Scholars of Zhejiang Province (LR22H260001) and the National Nature Science Foundation of China (82204019) and Y.D. was funded by the Key Project of Traditional Chinese Medicine Science and Technology Plan of Zhejiang Province (GZY-ZJ-KJ-24077) and the National Natural Science Foundation of China (82001673 and 82272860).

Association of gut microbiota with lactose intolerance and coeliac disease: a two-sample Mendelian randomization study

Background and objectives

Lactose intolerance and coeliac disease are common clinical nutrient malabsorption disorders, with an unclear pathogenesis and limited therapeutic options. It is widely believed that the gut microbiota plays an important role in many digestive disorders, but its role in lactose intolerance and coeliac disease is not yet clear. This study aimed to investigate the correlation between gut microbiota and lactose intolerance and coeliac disease.

Materials and methods

This study utilized the genome-wide association study database to investigate the association between gut microbiota and lactose intolerance and coeliac disease using Mendelian randomization (MR). The robustness of our findings was confirmed through subsequent analyses including Cochrane's Q statistic, MR-Egger Intercept Regression, MR-PRESSO Global Test and Leave-one-out methods.

Results

By employing the inverse variance weighted method, we identified that family Veillonellaceae, genus Oxalobacter and Senegalimassilia were protective against lactose intolerance, whereas genus Anaerotruncus, Eubacterium rectale group and Ruminococcus2 were found to be risk factors for lactose intolerance. Regarding coeliac disease, class Bacilli and Gammaproteobacteria, family FamilyXIII and Veillonellaceae, genus Eisenbergiella, Lachnoclostridium, RuminococcaceaeUCG014 and Ruminococcus2 were identified as protective factors, while class Betaproteobacteria, genus Eubacterium xylanophilum group and Blautia were risk factors. Furthermore, reverse the MR analysis did not reveal any evidence of a causal relationship between lactose intolerance or coeliac disease and the bacteria identified in our study.

Conclusion

This study provides novel insights into exploring the role of gut microbiota in lactose intolerance and coeliac disease; however, further experiments investigations are required to elucidate the specific underlying mechanisms.

Impact of childhood maltreatment on aging: a comprehensive Mendelian randomization analysis of multiple age-related biomarkers

BACKGROUND

Childhood maltreatment (CM) is linked to long-term adverse health outcomes, including accelerated biological aging and cognitive decline. This study investigates the relationship between CM and various aging biomarkers: telomere length, facial aging, intrinsic epigenetic age acceleration (IEAA), GrimAge, HannumAge, PhenoAge, frailty index, and cognitive performance.

METHODS

We conducted a Mendelian randomization (MR) study using published GWAS summary statistics. Aging biomarkers included telomere length (qPCR), facial aging (subjective evaluation), and epigenetic age markers (HannumAge, IEAA, GrimAge, PhenoAge). The frailty index was calculated from clinical assessments, and cognitive performance was evaluated with standardized tests. Analyses included Inverse-Variance Weighted (IVW), MR Egger, and Weighted Median (WM) methods, adjusted for multiple comparisons.

RESULTS

CM was significantly associated with shorter telomere length (IVW: β = - 0.1, 95% CI - 0.18 to - 0.02, pFDR = 0.032) and increased HannumAge (IVW: β = 1.33, 95% CI 0.36 to 2.3, pFDR = 0.028), GrimAge (IVW: β = 1.19, 95% CI 0.19 to 2.2, pFDR = 0.040), and PhenoAge (IVW: β = 1.4, 95% CI 0.12 to 2.68, pFDR = 0.053). A significant association was also found with the frailty index (IVW: β = 0.31, 95% CI 0.13 to 0.49, pFDR = 0.006). No significant associations were found with facial aging, IEAA, or cognitive performance.

CONCLUSIONS

CM is linked to accelerated biological aging, shown by shorter telomere length and increased epigenetic aging markers. CM was also associated with increased frailty, highlighting the need for early interventions to mitigate long-term effects. Further research should explore mechanisms and prevention strategies.

Pourcain B, Robinson EB, Bourgeron T, Baron-Cohen S, Børglum AD, Martin HC, Warrier V. An axis of genetic heterogeneity in autism is indexed by age at diagnosis and is associated with varying developmental and mental health profiles

There is growing recognition that earliest signs of autism need not clearly manifest in the first three years of life. To what extent is this variation in developmental trajectories associated with age at autism diagnosis? Does the genetic profile of autism vary with age at autism diagnosis? Using longitudinal data from four birth cohorts, we demonstrate that two different trajectories of socio-emotional behaviours are associated with age at diagnosis. We further demonstrate that the age at autism diagnosis is partly heritable (h2 SNP = 0.12, s.e.m = 0.01), and is associated with two moderately correlated (rg = 0.38, s.e.m = 0.07) autism polygenic factors. One of these factors is associated with earlier diagnosis of autism, lower social and communication abilities in early childhood. The second factor is associated with later autism diagnosis, increased socio-emotional difficulties in adolescence, and has moderate to high positive genetic correlations with Attention-Deficit/Hyperactivity Disorder, mental health conditions, and trauma. Overall, our research identifies an axis of heterogeneity in autism, indexed by age at diagnosis, which partly explains heterogeneity in autism and the profiles of co-occurring neurodevelopmental and mental health profiles. Our findings have important implications for how we conceptualise autism and provide one model to explain some of the diversity within autism.

Estimating the direct effects of the genetic liabilities to bipolar disorder, schizophrenia, and behavioral traits on suicide attempt using a multivariable Mendelian randomization approach

Bipolar disorder (BD) and schizophrenia (SZ) are associated with higher odds of suicide attempt (SA). In this study, we aimed to explore the effect of BD and SZ genetic liabilities on SA, also considering the contribution of behavioral traits, socioeconomic factors, and substance use disorders. Leveraging large-scale genome-wide association data from the Psychiatric Genomics Consortium (PGC) and the UK Biobank (UKB), we conducted a two-sample Mendelian randomization (MR) analysis to evaluate the putative causal effect of BD (41,917 cases, 371,549 controls) and SZ (53,386 cases, 77,258 controls) on SA (26,590 cases, 492,022 controls). Then, we assessed the putative causal effect of BD and SZ on behavioral traits, socioeconomic factors, and substance use disorders. Considering the associations identified, we evaluated the direct causal effect of behavioral traits, socioeconomic factors, and substance use disorders on SA using a multivariable MR approach. The genetic liabilities to BD and SZ were associated with higher odds of SA (BD odds ratio (OR) = 1.24, p = 3.88 × 10-12; SZ OR = 1.09, p = 2.44 × 10-20). However, while the effect of mental distress (OR = 1.17, p = 1.02 × 10-4) and risk-taking (OR = 1.52, p = 0.028) on SA was independent of SZ genetic liability, the BD-SA relationship appeared to account for the effect of these risk factors. Similarly, the association with loneliness on SA was null after accounting for the effect of SZ genetic liability. These findings highlight the complex interplay between genetic risk of psychiatric disorders and behavioral traits in the context of SA, suggesting the need for a comprehensive mental health assessment for high-risk individuals.

Association between cathepsins and skin cancers: A bidirectional two-sample Mendelian randomization study

BACKGROUND

Several cathepsins have been identified as being involved in the development of cancer. Nevertheless, the connection between cathepsins and skin cancers remained highly elusive.

METHODS

A bidirectional Mendelian randomization (MR) analysis was performed to investigate the causal association between cathepsins and skin malignancies. The genome-wide association studies (GWAS) data for cathepsins, malignant melanoma (MM), and basal cell carcinoma (BCC) were obtained from European research. The primary method employed was inverse variance weighted. In addition, MR-Egger, weighted median, weighted mode, and simple mode were also executed. Sensitivity analysis was performed using Cochran's Q test, MR-Egger, and MR-PRESSO.

RESULTS

From univariable MR (UVMR), cathepsin H, and S were determined to have a causal relationship with BCC. Additionally, cathepsin H was identified as associated with MM. Multivariable MR (MVMR) showed that after correcting for risk factors of skin carcinoma, cathepsin H was detected to be protective against BCC, whereas cathepsin S has been observed as a risk factor for BCC. No substantial pleiotropy and heterogeneity were identified in the sensitivity analysis.

CONCLUSION

This study was the first to establish a direct link between cathepsins and skin malignancies. Cathepsin H and S have the potential to serve as new biomarkers for BCC, offering valuable assistance in the prompt identification, treatment, and prevention of the disease. Nevertheless, additional clinical trials are required to validate our findings.

Association of glucagon-like peptide-1 receptor agonists with risk of cancers-evidence from a drug target Mendelian randomization and clinical trials

BACKGROUND

Glucagon-like peptide-1 receptor (GLP1R) agonists have been approved by Food and Drug Administration for management of obesity. However, the causal relationship of GLP1R agonists (GLP1RA) with cancers still unclear.

METHODS

The available cis-eQTLs for drugs target genes (GLP1R) were used as proxies for exposure to GLP1RA. Mendelian randomizations (MR) were performed to reveal the association of genetically-proxied GLP1RA with 14 common types cancer from large-scale consortia. Type 2 diabetes was used as positive control, and the GWASs data including 80 154 cases and 853 816 controls. Replicating the findings in the FinnGen study and then pooled with meta-analysis. Finally, all the related randomized controlled trails (RCTs) on GLP1RA were systematically searched from PubMed, Embase, and the Cochrane Library to comprehensively synthesize the evidence to validate any possible association with cancers.

RESULT

A total of 22 significant cis-eQTL single-nucleotide polymorphisms were included as genetic instrument. The association of genetically-proxied GLP1RA with significantly decreased type 2 diabetes risk [OR (95%)=0.82 (0.79-0.86), P <0.001], which ensuring the effectiveness of identified genetic instruments. The authors found favorable evidence to support the association of GLP1RA with reduced breast cancer and basal cell carcinoma risk [0.92 (0.88-0.96), P <0.001, 0.92 (0.85-0.99), P =0.029, respectively], and with increased colorectal cancer risk [1.12 (1.07-1.18), P <0.001]. In addition, there was no suggestive evidence to support the association of GLP1RA with ovarian cancer [0.99 (0.90-1.09), P =0.827], lung cancer [1.01 (0.93-1.10), P =0760], and thyroid cancer [0.83 (0.63-1.10), P =0.187]. Our findings were consistent with the meta-analysis. Finally, 80 RCTs were included in the systematic review, with a low incidence of different kinds of cancer.

CONCLUSIONS

Our study suggests that GLP1RA may decrease the risk of breast cancer and basal cell carcinoma, but increase the risk of colorectal cancer. However, according to the systematic review of RCTs, the incidence of cancer in patients treated with GLP1RA is low. Larger sample sizes of RCTs with long-term follow-up are necessary to establish the incidence of cancers and evaluate the risk-benefit ratios.

Exploring the bidirectional causal associations between pain and circulating inflammatory proteins: A Mendelian randomization study

Multisite chronic pain (MCP) and site-specific chronic pain (SSCP) may be influenced by circulating inflammatory proteins, but the causal relationship remains unknown. To overcome this limitation, two-sample bidirectional Mendelian randomization (MR) analysis was used to analyse data for 91 circulating inflammatory proteins, MCP and SSCP encompassing headache, back pain, shoulder pain, hip pain, knee pain, stomach abdominal pain and facial pain. The primary MR method used was inverse variance weighting, sensitivity analyses included weighted median, MR pleiotropy residual sum and outlier and the Egger intercept method. Heterogeneity was also detected using Cochrane's Q test and leave-one-out analyses. Finally, a causal relationship between 29 circulating inflammatory proteins and chronic pain was identified. Among these proteins, 14 exhibited a protective effect, including MCP (T-cell surface glycoprotein cluster of differentiation 5), headache (4E-binding protein 1 [4EBP1], cluster of differentiation 40, cluster of differentiation 6 and C-X-C motif chemokine [CXCL] 11), back pain (leukaemia inhibitory factor), shoulder pain (fibroblast growth factor [FGF]-5 and interleukin [IL]-18R1), stomach abdominal pain (tumour necrosis factor [TNF]-α), hip pain (CXCL1, IL-20 and signalling lymphocytic activation molecule 1) and knee pain (IL-7 and TNF-β). Additionally, 15 proteins were identified as risk factors for MCP and SSCP: MCP (colony-stimulating factor 1, human glial cell line-derived neurotrophic factor and IL-17C), headache (fms-related tyrosine kinase 3 ligand, IL-20 receptor subunit α [IL-20RA], neurotrophin-3 and tumour necrosis factor receptor superfamily member 9), facial pain (CXCL1), back pain (TNF), shoulder pain (IL-17C and matrix metalloproteinase-10), stomach abdominal pain (IL-20RA), hip pain (C-C motif chemokine 11/eotaxin-1 and tumour necrosis factor ligand superfamily member 12) and knee pain (4EBP1). Importantly, in the opposite direction, MCP and SSCP did not exhibit a significant causal impact on circulating inflammatory proteins. Our study identified potential causal influences of various circulating inflammatory proteins on MCP and SSCP and provided promising treatments for the clinical management of MCP and SSCP.

Genetic and causal relationship between chronic gastrointestinal diseases and erectile dysfunction: a Mendelian randomization study

Background

Studies based on observations have indicated potential associations between chronic gastrointestinal diseases and an increased risk of erectile dysfunction (ED). However, the causality of these connections remains ambiguous.

Methods

Summary data for chronic gastrointestinal diseases were extracted from public data. Summary data on ED were extracted from three distinct sources. The genetic correlations between chronic gastrointestinal diseases and ED were explored using linkage disequilibrium score regression (LDSC). The causal associations between chronic gastrointestinal diseases and ED were evaluated using Mendelian randomization (MR) analysis, followed by a meta-analysis to determine the ultimate causal effect.

Results

The LDSC results suggested significant genetic correlations between Crohn's disease (CD) and ED. Inflammatory bowel disease (IBD), ulcerative colitis (UC), and liver cirrhosis (LC) were found to have potential genetic correlations with ED. The combined multiple MR results indicate that IBD and CD have significant causal relationships with ED, while colorectal cancer (CRC) may have a potential causal effect on ED.

Conclusion

This research provided evidence supporting a causal association between IBD, CD, CRC, and ED. The impact of chronic gastrointestinal diseases on ED warrants greater attention in clinical practice.

The role of hypertension in the relationship between leisure screen time, physical activity and migraine: a 2-sample Mendelian randomization study

BACKGROUND

The relationship between lifestyle and migraine is complex, as it remains uncertain which specific lifestyle factors play the most prominent role in the development of migraine, or which modifiable metabolic traits serve as mediators in establishing causality.

METHODS

Independent genetic variants strongly associated with 20 lifestyle factors were selected as instrumental variables from corresponding genome-wide association studies (GWASs). Summary-level data for migraine were obtained from the FinnGen consortium (18,477 cases and 287,837 controls) as a discovery set and the GWAS meta-analysis data (26,052 cases and 487,214 controls) as a replication set. Estimates derived from the two datasets were combined using fixed-effects meta-analysis. Two-step univariable MR (UVMR) and multivariable Mendelian randomization (MVMR) analyses were conducted to evaluate 19 potential mediators of association and determine the proportions of these mediators.

RESULTS

The combined effect of inverse variance weighted revealed that a one standard deviation (SD) increase in genetically predicted Leisure screen time (LST) was associated with a 27.7% increase (95% CI: 1.14-1.44) in migraine risk, while Moderate or/and vigorous physical activity (MVPA) was associated with a 26.9% decrease (95% CI: 0.61-0.87) in migraine risk. The results of the mediation analysis indicated that out of the 19 modifiable metabolic risk factors examined, hypertension explains 24.81% of the relationship between LST and the risk of experiencing migraine. Furthermore, hypertension and diastolic blood pressure (DBP) partially weaken the association between MVPA and migraines, mediating 4.86% and 4.66% respectively.

CONCLUSION

Our research findings indicated that both LST and MVPA in lifestyle have independent causal effects on migraine. Additionally, we have identified that hypertension and DBP play a mediating role in the causal pathway between these two factors and migraine.

Causality between six psychiatric disorders and digestive tract cancers risk: a two-sample Mendelian randomization study

Associations between psychiatric disorders and digestive tract cancers have been proposed. However, the causal link between these factors remains unclear. This study pioneers Mendelian randomization (MR) analysis to explore the genetic link between psychiatric disorders and digestive tract cancers risk. We analysed data on six psychiatric disorders [schizophrenia, bipolar disorder, major depressive disorder (MDD), attention deficit hyperactivity disorder, autism spectrum disorder, and panic disorder (PD)] and digestive tract cancers [esophagus cancer (EC), gastric cancer (GC), and colorectal cancer (CRC)] from genome-wide association studies databases. Using instrumental variables identified from significant single nucleotide polymorphism associations, we employed the inverse variance weighted (IVW) method alongside the weighted median (WM) method and MR-Egger regression. The results revealed no causal link between psychiatric disorders and the risk of EC or GC. Psychiatric disorders were not identified as risk factors for CRC. Notably, PD demonstrated a lower CRC risk (OR = 0.79, 95% CI 0.66-0.93, P = 0.01). This MR analysis underscores the lack of a causal association between psychiatric disorders and digestive tract cancers risk while suggesting a potential protective effect of PD against CRC.

Genome-wide meta-analysis identifies ancestry-specific loci for Alzheimer's disease

INTRODUCTION

Alzheimer's disease (AD) is a devastating neurological disease with complex genetic etiology. Yet most known loci have only identified from the late-onset type AD in populations of European ancestry.

METHODS

We performed a two-stage genome-wide association study (GWAS) of AD totaling 6878 Chinese and 63,926 European individuals.

RESULTS

In addition to the apolipoprotein E (APOE) locus, our GWAS of two independent Chinese samples uncovered three novel AD susceptibility loci (KIAA2013, SLC52A3, and TCN2) and a novel ancestry-specific variant within EGFR (rs1815157). More replicated variants were observed in the Chinese (31%) than in the European samples (15%). In combining genome-wide associations and functional annotations, EGFR and TCN2 were prioritized as two of the most biologically significant genes. Phenome-wide Mendelian randomization suggests that high mean corpuscular hemoglobin concentration might protect against AD.

DISCUSSION

The current study reveals novel AD susceptibility loci, emphasizes the importance of diverse populations in AD genetic research, and advances our understanding of disease etiology.

HIGHLIGHTS

Loci KIAA2013, SLC52A3, and TCN2 were associated with Alzheimer's disease (AD) in Chinese populations. rs1815157 within the EGFR locus was associated with AD in Chinese populations. The genetic architecture of AD varied between Chinese and European populations. EGFR and TCN2 were prioritized as two of the most biologically significant genes. High mean corpuscular hemoglobin concentrations might have protective effects against AD.

Relationships of omega-3 and omega-6 polyunsaturated fatty acids with esophageal diseases: a two-sample Mendelian randomization analysis

Introduction

Omega-3 polyunsaturated fatty acids (PUFAs) have been widely studied and used as nutritional supplements because of their anti-inflammatory effects. Previous studies have shown an association between polyunsaturated fatty acids such as omega-3 and omega-6 PUFAs with the development of malignant tumors. However, the relationships of omega-3 and omega-6 PUFAs with esophageal diseases have not been characterized.

Methods

Mendelian randomization (MR) is a statistical method for identifying instrumental variables (IVs) from genome-wide association study (GWAS) data, and is associated with little confounding by environmental or other disease-related factors. We used genome-wide association study (GWAS) data from previously published studies on circulating concentrations of omega-3, omega-6, docosahexaenoic acid (DHA) and linoleic acid (LA), as well as esophageal cancer and other esophageal diseases, which were downloaded from the IEU OpenGwas database (https://gwas.mrcieu.ac.uk/) and the GWAS Catalog database (https://www.ebi.ac.uk/). The inverse variance-weighted approach was used as the principal analysis, and the MR-Egger and weighted median methods were used alongside. A series of sensitivity analyses were used to ensure the robustness of the causality estimates.

Results

We found that the circulating omega-3 PUFAs concentration was positively associated with esophageal cancer (p = 8 × 10-4), and circulating DHA concentration (the main component of omega-3 in food), was also positively associated with esophageal cancer (p = 2 × 10-2), but no significant association was found between circulating omega-6 PUFAs and esophageal cancer (p = 0.17), and circulating LA concentration (the main component of omega-6 in food), was also no significant associated with esophageal cancer (p = 0.32). We found no significant relationships of circulating omega-3 and omega-6 PUFAs concentration with four other esophageal diseases.

Conclusion

This study indicates that higher levels of circulating omega-3 PUFAs and DHA concentrations may be a risk factor for the development of esophageal cancer. Conversely, an increased omega-6/omega-3 ratio may serve as a protective factor against esophageal cancer. These findings have significant implications for the clinical application of omega-3 PUFAs and the prevention and treatment of esophageal cancer.

Multi-omics characterization of type 2 diabetes associated genetic variation

Discerning the mechanisms driving type 2 diabetes (T2D) pathophysiology from genome-wide association studies (GWAS) remains a challenge. To this end, we integrated omics information from 16 multi-tissue and multi-ancestry expression, protein, and metabolite quantitative trait loci (QTL) studies and 46 multi-ancestry GWAS for T2D-related traits with the largest, most ancestrally diverse T2D GWAS to date. Of the 1,289 T2D GWAS index variants, 716 (56%) demonstrated strong evidence of colocalization with a molecular or T2D-related trait, implicating 657 cis-effector genes, 1,691 distal-effector genes, 731 metabolites, and 43 T2D-related traits. We identified 773 of these cis- and distal-effector genes using either expression QTL data from understudied ancestry groups or inclusion of T2D index variants enriched in underrepresented populations, emphasizing the value of increasing population diversity in functional mapping. Linking these variants, genes, metabolites, and traits into a network, we elucidated mechanisms through which T2D-associated variation may impact disease risk. Finally, we showed that drugs targeting effector proteins were enriched in those approved to treat T2D, highlighting the potential of these results to prioritize drug targets for T2D. These results represent a leap in the molecular characterization of T2D-associated genetic variation and will aid in translating genetic findings into novel therapeutic strategies.

Causal effects of gut microbiota on risk of interstitial cystitis: a two-sample Mendelian randomization study

Background

The correlation between gut microbiota and interstitial cystitis has garnered significant attention in previous studies. Nevertheless, the causal relationship between them remains to be clarified.

Methods

Genetic variation serves as a tool in Mendelian randomization analyses, facilitating the inference of causal relationships between exposure factors and disease outcomes. In this study, summary statistics derived from a comprehensive genome-wide association study conducted by the MiBioGen consortium were utilized as exposure factors, while interstitial cystitis data sourced from the GWAS Catalog served as the disease outcome. Then, a two-sample Mendelian randomization analysis was performed by applying inverse variance-weighted, MR-Egger, Weighted Median, Simple Mode, and Weighted Mode. In addition, heterogeneity and horizontal pleiotropy were excluded by sensitivity analysis.

Results

IVW results confirmed that genus Haemophilus (OR = 2.20, 95% CI: 1.16-4.15, p = 0.015), genus Butyricimonas (OR = 2.26, 95% CI: 1.15-4.45, p = 0.018), genus Bacteroides (OR = 4.27, 95% CI: 1.36-13.4, p = 0.013) and Coprococcus1 (OR = 3.39, 95% CI: 1.28-8.99, p = 0.014) had a risk effect on interstitial cystitis. Sensitivity analysis did not find outlier SNPs.

Conclusion

Our analysis has identified a causal relationship between specific genera and interstitial cystitis. However, further validation through randomized controlled trials is essential to substantiate these findings.

Polygenic risk score of metabolic dysfunction-associated steatotic liver disease amplifies the health impact on severe liver disease and metabolism-related outcomes

BACKGROUND

Although the inherited risk factors associated with fatty liver disease are well understood, little is known about the genetic background of metabolic dysfunction-associated steatotic liver disease (MASLD) and its related health impacts. Compared to non-alcoholic fatty liver disease (NAFLD), MASLD presents significantly distinct diagnostic criteria, and epidemiological and clinical features, but the related genetic variants are yet to be investigated. Therefore, we conducted this study to assess the genetic background of MASLD and interactions between MASLD-related genetic variants and metabolism-related outcomes.

METHODS

Participants from the UK Biobank were grouped into discovery and replication cohorts for an MASLD genome-wide association study (GWAS), and base and target cohorts for polygenic risk score (PRS) analysis. Autosomal genetic variants associated with NAFLD were compared with the MASLD GWAS results. Kaplan-Meier and Cox regression analyses were used to assess associations between MASLD and metabolism-related outcomes.

RESULTS

Sixteen single-nucleotide polymorphisms (SNPs) were identified at genome-wide significance levels for MASLD and duplicated in the replication cohort. Differences were found after comparing these SNPs with the results of NAFLD-related genetic variants. MASLD cases with high PRS had a multivariate-adjusted hazard ratio of 3.15 (95% confidence interval, 2.54-3.90) for severe liver disease (SLD), and 2.81 (2.60-3.03) for type 2 diabetes mellitus. The high PRS amplified the impact of MASLD on SLD and extrahepatic outcomes.

CONCLUSIONS

High PRS of MASLD GWAS amplified the impact of MASLD on SLD and metabolism-related outcomes, thereby refining the process of identification of individuals at high risk of MASLD. Supplementation of this process with relevant genetic backgrounds may lead to more effective MASLD prevention and management.

Associations between modifiable risk factors and hepatocellular carcinoma: a trans-ancestry Mendelian randomization study

BACKGROUND

Potentially modifiable risk factors for hepatocellular carcinoma (HCC) have been investigated in observational epidemiology studies in East Asian and European populations, whereas the causal associations of most of these risk factors remain unclear.

METHODS

We collected genome-wide association summary statistics of 22 modifiable risk factors in East Asians and 33 risk factors in Europeans. Genetic summary statistics of HCC were sourced from the Biobank Japan study (1,866 cases and 195,745 controls) for East Asians, and the deCODE genetics study (406 cases and 49,302 controls) and the UK Biobank (168 cases and 372 016 controls) for Europeans. Two-sample Mendelian randomization (MR) analyses were performed independently for East Asian and European populations.

RESULTS

In East Asians, genetically predicted alcohol frequency, ever drinkers, aspartate aminotransferase (AST), hypothyroidism, chronic hepatitis B, and chronic hepatitis C, metabolic dysfunction-associated steatotic liver disease (MASLD), and autoimmune hepatitis were significantly associated with an increased HCC risk (P < 0.05/22). Among European population, alanine transaminase, AST, MASLD, percent liver fat, and liver iron content were significantly associated with a higher risk of HCC (P < 0.05/33). The replication dataset and meta-analysis further confirmed these results.

CONCLUSIONS

Although East Asian and European populations have different factors for HCC, their common modifiable risk factors AST and MASLD for HCC, offer valuable insights for targeted intervention strategies to mitigate society burden of HCC.

Association between blood metabolites and basal cell carcinoma risk: a two-sample Mendelian randomization study

Background

Circulating metabolites, which play a crucial role in our health, have been reported to be disordered in basal cell carcinoma (BCC). Despite these findings, evidence is still lacking to determine whether these metabolites directly promote or prevent BCC's progression. Therefore, our study aims to examine the potential effects of circulating metabolites on BCC progression.

Material and methods

We conducted a two-sample Mendelian randomization (MR) analysis using data from two separate genome-wide association studies (GWAS). The primary study included data for 123 blood metabolites from a GWAS with 25,000 Finnish individuals, while the secondary study had data for 249 blood metabolites from a GWAS with 114,000 UK Biobank participants.GWAS data for BCC were obtained from the UK Biobank for the primary analysis and the FinnGen consortium for the secondary analysis. Sensitivity analyses were performed to assess heterogeneity and pleiotropy.

Results

In the primary analysis, significant causal relationships were found between six metabolic traits and BCC with the inverse variance weighted (IVW) method after multiple testing [P < 4 × 10-4 (0.05/123)]. Four metabolic traits were discovered to be significantly linked with BCC in the secondary analysis, with a significance level of P < 2 × 10-4 (0.05/249). We found that all the significant traits are linked to Polyunsaturated Fatty Acids (PUFAs) and their degree of unsaturation.

Conclusion

Our research has revealed a direct link between the susceptibility of BCC and Polyunsaturated Fatty Acids and their degree of unsaturation. This discovery implies screening and prevention of BCC.

Relationships between nine neuropsychiatric disorders and cervical cancer: insights from genetics, causality and shared gene expression patterns

BACKGROUND

Neuropsychiatric disorders and cervical cancer exert substantial influences on women's health. Furthermore, neuropsychiatric disorders frequently manifest as common symptoms in cancer patients, potentially increasing the risk of malignant neoplasms. This study aimed to identify neuropsychiatric disorders that are genetically and causally related to cervical cancer and to investigate the molecular mechanisms underlying these associations.

METHODS

GWAS data related to nine neuropsychiatric disorders, namely, schizophrenia, bipolar disorder, autism spectrum disorder, Parkinson's disease, anxiety, Alzheimer's disease, mood disorders, depression, and alcohol dependence, were obtained to calculate heritability (h2) and genetic correlation (rg) with cervical cancer using linkage disequilibrium score regression (LDSC). Mendelian randomization (MR) analysis of the two cohorts was employed to assess the causal effects. Shared gene expression pattern analysis was subsequently conducted to investigate the molecular mechanism underlying these significant associations.

RESULTS

Anxiety, mood disorders, depression, and alcohol dependence were genetically correlated with cervical cancer (all adjusted P < 0.05). Only depression was causally related to cervical cancer in both the discovery (ORIVW: 1.41, PIVW = 0.02) and replication cohorts (ORIVW: 1.80, PIVW = 0.03) in the MR analysis. Gene expression pattern analysis revealed that 270 genes related to depression and cervical cancer, including tumour necrosis factor (TNF), were significantly upregulated in cervical cancer patients, while vascular endothelial growth factor A (VEGFA), transcription factor AP-1 (JUN), and insulin-like growth factor I (IGF-I) were associated with prognosis in cervical cancer patients (all P < 0.05). These overlapping genes implicated the involvement of multiple biological mechanisms, such as neuron death, the PI3K-Akt signalling pathway, and human papillomavirus infection.

CONCLUSIONS

Genetic, causal and molecular evidence indicates that depression increases the risk of cervical cancer. The TNF, VEGFA, JUN, and IGF-1 genes and the neuron death, PI3K-Akt, and human papillomavirus infection signalling pathways may possibly explain this association.

Socioeconomic disparity in the natural history of cutaneous melanoma: evidence from two large prospective cohorts

BACKGROUND

Previous studies on the associations between socioeconomic status (SES) and cutaneous malignant melanoma (CMM) failed to distinguish the effects of different SES factors under an individual-data-based prospective study design.

METHODS

Based on UK Biobank (UKB) and China Kadoorie Biobank (CKB), we estimated the effects of four SES factors on transitions from baseline to CMM in situ, subsequently to invasive CMM and further CMM mortality by applying multistate models. We further explored to which extent the associations between SES and CMM incidence could be explained by potential mediators including sun exposure, lifestyle and ageing in UKB.

RESULTS

In multistate analyses, good household income was independently associated with an increased risk of CMM in situ (HR=1.38, 95% CI: 1.21 to 1.58) and invasive CMM (HR=1.34, 95% CI: 1.22 to 1.48) in UKB. These findings were partly validated in CKB. Especially in UKB, we observed an increased risk of CMM in situ and invasive CMM among participants with good type of house; only good education was independently associated with lower risk of evolving to invasive CMM among patients with CMM in situ (HR=0.69, 95% CI: 0.52 to 0.92); only good household income was independently associated with lower risk of CMM mortality among patients with CMM (HR=0.65, 95% CI: 0.45 to 0.95). In mediation analysis, the proportions attributable to the mediating effect were <6% for all selected variables, including self-reported sun exposure-related factors.

CONCLUSION

SES factors have different effects on the incidence and progression of CMM. The association between SES and incident CMM is neither causal nor well explained by selected mediators.

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

Association between glioma and neurodegenerative diseases risk: a two-sample bi-directional Mendelian randomization analysis

Background

Earlier observational studies have demonstrated a correlation between glioma and the risk of neurodegenerative diseases (NDs), but the causality and direction of their associations remain unclear. The objective of this study was to ascertain the causal link between glioma and NDs using Mendelian randomization (MR) methodology.

Methods

Genome-wide association study (GWAS) data were used in a two-sample bi-directional MR analysis. From the largest meta-analysis GWAS, encompassing 18,169 controls and 12,488 cases, summary statistics data on gliomas was extracted. Summarized statistics for NDs, including Alzheimer's disease (AD), multiple sclerosis (MS), amyotrophic lateral sclerosis (ALS) and Parkinson's disease (PD) were obtained from the GWAS of European ancestry. Inverse variance weighted (IVW) method was elected as the core MR approach with weighted median (WM) method and MR-Egger method as complementary methods. In addition, sensitivity analyses were performed. A Bonferroni correction was used to correct the results.

Results

Genetically predicted glioma had been related to decreased risk of AD. Specifically, for all glioma (IVW: OR = 0.93, 95% CI = 0.90-0.96, p = 4.88 × 10-6) and glioblastoma (GBM) (IVW: OR = 0.93, 95% CI = 0.91-0.95, p = 5.11 × 10-9). We also found that genetically predicted all glioma has a suggestive causative association with MS (IVW: OR = 0.90, 95% CI = 0.81-1.00, p = 0.045). There was no evidence of causal association between glioma and ALS or PD. According to the results of reverse MR analysis, no discernible causal connection of NDs was found on glioma. Sensitivity analyses validated the robustness of the above associations.

Conclusion

We report evidence in support of potential causal associations of different glioma subtypes with AD and MS. More studies are required to uncover the underlying mechanisms of these findings.

ENPP2/Autotaxin: The potential drug target for alcoholic liver disease identified through Mendelian randomization analysis

BACKGROUND AND AIMS

At present, there is still a lack of radical drug targets for intervention in alcoholic liver disease (ALD), and drug discovery through randomized controlled trials is a lengthy, risky, and expensive undertaking, so we aimed to identify effective drug targets based on human genetics.

METHODS

We used Mendelian randomization (MR) and Bayesian colocalization analysis to investigate 2639 genes encoding druggable proteins and examined the causal effects on ALD (PMID 34737426: 456348 European with 451 cases and 455 897 controls). In addition, we conducted the mediation analysis to explore the potential mechanism using the genome-wide association study (GWAS) data of blood biomarkers as mediators.

RESULTS

We finally identified the drug target: ENPP2/Autotaxin and genetically proxied ENPP2/Autotaxin was causally associated with the risk of ALD (OR = 2.28, 95% CI: 1.64 to 3.16, p = 7.49E-7). In addition, we found that the effect of ENPP2/Autotaxin on ALD may be partly mediated by effector memory CD8+ T cell (the proportion of mediation effect: 8.49%).

CONCLUSIONS

Our integrative analysis suggested that genetically determined levels of circulating ENPP2/Autotaxin have a causal effect on ALD risk and are a promising drug target.

Summary data-based Mendelian randomization and single-cell RNA sequencing analyses identify immune associations with low-level LGALS9 in sepsis

Sepsis is one of the major challenges in intensive care units, characterized by the complexity of the host immune status. To gain a deeper understanding of the pathogenesis of sepsis, it is crucial to study the phenotypic changes in immune cells and their underlying molecular mechanisms. We conducted Summary data-based Mendelian randomization analysis by integrating genome-wide association studies data for sepsis with expression quantitative trait locus data, revealing a significant decrease in the expression levels of 17 biomarkers in sepsis patients. Furthermore, based on single-cell RNA sequencing data, we elucidated potential molecular mechanisms at single-cell resolution and identified that LGALS9 inhibition in sepsis patients leads to the activation and differentiation of monocyte and T-cell subtypes. These findings are expected to assist researchers in gaining a more in-depth understanding of the immune dysregulation in sepsis.

Casual relationships between circulating metabolites and rheumatoid arthritis: A mendelian randomization study

Background

Blood metabolites serve as pivotal indicators in identifying and predicting the course of rheumatoid arthritis (RA). However, empirical substantiation of a direct causal link between these serum biomarkers and the development of RA is still lacking comprehensive support.

Method

In pursuit of a thorough exploration of the causal links between circulating blood metabolites and RA, we deployed a two-sample Mendelian randomization (MR) approach during our initial investigative phase. This method was utilized to examine the potential connections between 249 distinct circulating metabolites and the prevalence of RA. In the validation phase, we conducted replication analyses with a new metabolic dataset consisting of 123 metabolites. Furthermore, we employed the Mendelian randomization based on Bayesian model averaging (MR-BMA) technique to pinpoint key metabolic characteristics that have significant causal implications.

Results

In our primary analysis, we found that acetate, acetoacetate and pyruvate exhibited a consistent protective causal association with rheumatoid arthritis, while lactate demonstrated a positive correlation with rheumatoid arthritis risk. It is also noteworthy that a substantial subset of traits related to both saturated and unsaturated fatty acids showed causal influences. Subsequent secondary analyses substantiated these observations, revealing that traits associated with the average number of methylene groups in a fatty acid chain exhibited protective effects. Ultimately, our MR-BMA analyses unveiled that the ratio of polyunsaturated fatty acids (PUFAs) to total fatty acids assumes a paramount role in increasing the susceptibility to rheumatoid arthritis.

Conclusions

By employing systemic MR analyses, our study has successfully generated an all-encompassing atlas elucidating the intricate connections between circulating metabolites and the susceptibility to rheumatoid arthritis. Our results indicate the high unsaturation degree is a dominant risk factors correlated with rheumatoid arthritis.

Genetic susceptibility to breast cancer increases the risk of neutropenia and agranulocytosis: insights from Mendelian randomization

PURPOSE

The causal relationship between breast cancer and its estrogen receptor (ER) subtypes and neutropenia and agranulocytosis is unclear.

METHODS

In two-sample Mendelian randomization (MR), we used inverse variance weighting (IVW), Bayesian weighted MR (BWMR), MR-Egger, weighted median, simple mode, and weighted mode methods to analyze causality for ER-positive breast cancer, ER-negative breast cancer, overall breast cancer, and drug-induced neutropenia and agranulocytosis. To validate the results, we performed the analysis again using GWAS data on neutropenia from different databases. In multivariable MR (MVMR), we assessed the independent effects of ER-positive and ER-negative breast cancer on causality.

RESULTS

Two-sample MR analysis showed a causal relationship between ER-positive breast cancer (IVW odds ratio (OR) = 1.319, P = 7.580 × 10-10), ER-negative breast cancer (OR = 1.285, P = 1.263 × 10-4), overall breast cancer (OR = 1.418, P = 2.123 × 10-13), and drug-induced neutropenia and a causal relationship between ER-positive breast cancer (OR = 1.349, P = 1.402 × 10-7), ER-negative breast cancer (OR = 1.235, P = 7.615 × 10-3), overall breast cancer (OR = 1.429, P = 9.111 × 10-10), and neutropenia. Similarly, ER-positive breast cancer (OR = 1.213, P = 5.350 × 10-8), ER-negative breast cancer (OR = 1.179, P = 1.300 × 10-3), and overall breast cancer (OR = 1.275, P = 8.642 × 10-11) also had a causal relationship with agranulocytosis. MVMR analysis showed that ER-positive breast cancer remained causally associated with drug-induced neutropenia (OR = 1.233, P = 4.188 × 10-4), neutropenia (OR = 1.283, P = 6.363 × 10-4), and agranulocytosis (OR = 1.142, P = 4.549 × 10-3). Heterogeneity analysis and pleiotropy test showed that our results were reliable.

CONCLUSION

Our study provides genetic evidence for a causal association between breast cancer and its estrogen receptor subtypes and neutropenia. In clinical practice, in addition to focusing on therapeutic factors, additional attention should be given to breast cancer patients to avoid severe neutropenia.

The relationship between obstructive sleep apnea and osteoarthritis: evidence from an observational and Mendelian randomization study

Objectives

Obstructive sleep apnea (OSA) and osteoarthritis (OA) are common comorbidities that significantly impact individuals' quality of life. However, the relationship between OSA and OA remains unclear. This study aims to explore the connection between OSA and OA and evaluate causality using Mendelian randomization (MR).

Methods

A total of 12,454 participants from the National Health and Nutrition Examination Survey (2009-2012) were included. OSA participants were identified based on self-reported interviews. The association between OA and OSA was assessed through multivariable logistic regression analysis. A two-sample MR was employed to investigate the relationship between OSA and OA, specifically hip OA and knee OA, utilizing the inverse variance-weighted (IVW) approach.

Results

Based on the observational study, individuals with OSA exhibited a higher risk of OA (OR = 1.67, 95% CI = 1.40-1.98). IVW demonstrated that the risk of OA (OR = 1.13, 95% CI: 1.05-1.21, p = 0.001), hip OA (OR = 1.11, 95% CI: 1.04-1.18, p = 0.002), and knee OA (OR = 1.08, 95% CI: 1.02-1.14, p = 0.005) was significantly associated with OSA. Reverse MR analyses indicated no effect of OA on OSA. Additionally, body mass index (BMI) was found to mediate 36.9% (95% CI, 4.64-73.2%, p = 0.026) of the OSA effects on OA risk.

Conclusion

The cross-sectional observational analysis unveiled noteworthy associations between OSA and OA. Meanwhile, findings from the MR study provide support for a causal role.

Prioritization of drug targets for thyroid cancer: a multi-omics Mendelian randomization study

OBJECTIVES

Recurrence or tumor metastasis and drug resistance remain the major challenge in the treatment of thyroid cancer. It is needed to identify novel drug targets for thyroid cancer.

METHODS

Summary data-based Mendelian randomization (SMR) and colocalization analysis were performed to evaluate the associations between gene methylation, expression, protein levels with thyroid cancer. We additionally performed protein-protein interaction (PPI) network, gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) analyses to further explore the potential roles of identified genes in thyroid cancer.

RESULTS

SDCCAG8 and VCAM1 genes were associated with risk of thyroid cancer with tier 1 evidence, while TCN2 gene was with tier 3 evidence. SDCCAG8 gene was associated with risk of papillary thyroid cancer with tier 1 evidence. At the level of circulating proteins, genetically predicted higher levels of SDCCAG8 (OR = 0.46, 95% CI 0.34-0.64) and VCAM1 (OR = 0.21, 95% CI 0.10-0.45) were inversely associated with thyroid cancer risk; higher level of TCN2 was associated with an increased risk of thyroid cancer (OR = 1.30, 95% CI 1.15-1.47); and the higher level of SDCCAG8 (OR = 0.40, 95% CI 0.28-0.58) was associated with a decreased risk of papillary thyroid cancer. The bioinformatics analysis showed that SDCCAG8, VCAM1 and TCN2 might play roles in immune-related pathways.

CONCLUSION

SDCCAG8, VCAM1 and TCN2 genes were associated with thyroid cancer risk with evidence at multi-omics levels. There were potential roles of SDCCAG8, VCAM1 and TCN2 in immune-related pathways. Our findings might improve the understanding of the pathogenesis of thyroid cancer and discovery of novel potential drug targets for this disease.

Genetic influences on brain and cognitive health and their interactions with cardiovascular conditions and depression

Approximately 40% of dementia cases could be prevented or delayed by modifiable risk factors related to lifestyle and environment. These risk factors, such as depression and vascular disease, do not affect all individuals in the same way, likely due to inter-individual differences in genetics. However, the precise nature of how genetic risk profiles interact with modifiable risk factors to affect brain health is poorly understood. Here we combine multiple data resources, including genotyping and postmortem gene expression, to map the genetic landscape of brain structure and identify 367 loci associated with cortical thickness and 13 loci associated with white matter hyperintensities (P < 5×10-8), with several loci also showing a significant association with cognitive function. We show that among 220 unique genetic loci associated with cortical thickness in our genome-wide association studies (GWAS), 95 also showed evidence of interaction with depression or cardiovascular conditions. Polygenic risk scores based on our GWAS of inferior frontal thickness also interacted with hypertension in predicting executive function in the Canadian Longitudinal Study on Aging. These findings advance our understanding of the genetic underpinning of brain structure and show that genetic risk for brain and cognitive health is in part moderated by treatable mid-life factors.

Genetically predicted 91 circulating inflammatory proteins in relation to risk of urological malignancies: a Mendelian randomization study

BACKGROUND

Urological malignancies, including kidney, bladder, and prostate cancer, are major health concerns worldwide. Inflammation has been implicated in the pathogenesis of these cancers, and circulating inflammatory proteins may play a role in their development. However, the causal relationship between specific plasma proteins and urological malignancies remains unclear.

METHODS

We performed a two-sample Mendelian randomization (MR) analysis using summary statistics from genome-wide association studies (GWAS). Instrumental variables representing genetic variants associated with circulating inflammatory proteins were used to infer causality on the risk of kidney, bladder, and prostate cancer. Four MR methods were utilized to provide robust effect estimates.

RESULTS

Our analysis identified several plasma proteins associated with a lower risk of kidney and bladder cancer, including Eukaryotic translation initiation factor 4E-binding protein 1, Caspase 8, Natural killer cell receptor 2B4, and Tumor necrosis factor ligand superfamily member 12. However, after adjusting for multiple testing, these associations did not remain statistically significant. For prostate cancer, CUB domain-containing protein 1 and Interleukin-10 receptor subunit beta were found to be protective, while Glial cell line-derived neurotrophic factor and SIR2-like protein 2 were identified as risk factors. After FDR adjustment, none of the inflammatory proteins were found to be significantly associated with a lower risk of prostate cancer.

CONCLUSION

Our findings suggest that certain plasma proteins may be involved in the development of urological malignancies. Mendelian randomization provides a useful framework for investigating causal relationships between inflammatory proteins and urological cancers, offering potential insights into their underlying biology and therapeutic targets.

Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration

OBJECTIVES

Single-cell and spatial transcriptomics analysis of human knee articular cartilage tissue to present a comprehensive transcriptome landscape and osteoarthritis (OA)-critical cell populations.

METHODS

Single-cell RNA sequencing and spatially resolved transcriptomic technology have been applied to characterise the cellular heterogeneity of human knee articular cartilage which were collected from 8 OA donors, and 3 non-OA control donors, and a total of 19 samples. The novel chondrocyte population and marker genes of interest were validated by immunohistochemistry staining, quantitative real-time PCR, etc. The OA-critical cell populations were validated through integrative analyses of publicly available bulk RNA sequencing data and large-scale genome-wide association studies.

RESULTS

We identified 33 cell population-specific marker genes that define 11 chondrocyte populations, including 9 known populations and 2 new populations, that is, pre-inflammatory chondrocyte population (preInfC) and inflammatory chondrocyte population (InfC). The novel findings that make this an important addition to the literature include: (1) the novel InfC activates the mediator MIF-CD74; (2) the prehypertrophic chondrocyte (preHTC) and hypertrophic chondrocyte (HTC) are potentially OA-critical cell populations; (3) most OA-associated differentially expressed genes reside in the articular surface and superficial zone; (4) the prefibrocartilage chondrocyte (preFC) population is a major contributor to the stratification of patients with OA, resulting in both an inflammatory-related subtype and a non-inflammatory-related subtype.

CONCLUSIONS

Our results highlight InfC, preHTC, preFC and HTC as potential cell populations to target for therapy. Also, we conclude that profiling of those cell populations in patients might be used to stratify patient populations for defining cohorts for clinical trials and precision medicine.

Mendelian Randomization Analysis of Circulating Cytokines and Risk of Autoimmune Neuroinflammatory Diseases

Background

Cytokines act a vital role in autoimmune neuroinflammatory diseases (ANDs) with undetermined causal relationships. Mendelian randomization (MR) analysis was performed to estimate the causal effects of circulating levels of cytokines on the risk of ANDs.

Methods

The causal relationship between 34 circulating cytokines and 4 kinds of ANDs, including multiple sclerosis (MS), neuromyelitis optica (NOM), chronic inflammatory demyelinating polyneuropathy (CIDP) and myasthenia gravis (MG) were explored using four methods of MR analysis. MR-PRESSO, MR-Egger regression methods and Cochran's Q statistic were utilized to identify the instrumental variables (IVs) with potential pleiotropy and heterogeneity. The Bonferroni correction was used for multiple group comparisons. P-value less than 3.68E-04 (0.05/ (34*4)) was considered statistically significant.

Results

Negative causal effects of circulating levels of interleukin (IL)-8 (OR = 0.648, 95% CI: 0.494-0.851, P = 0.002) on risk of MS, chemokine (C-C Motif) ligand (CCL)-5 (OR = 0.295, 95% CI: 0.103-0.841, P = 0.022) and stem cell growth factor-beta (SCGF-β) (OR = 0.745, 95% CI: 0.565-0.984, P = 0.038) on risk of CIDP, as well as positive causal effects of circulating levels of IL-2 receptor α (IL-2Rα) (OR = 1.216, 95% CI: 1.120-1.320, P = 3.20E-06) and chemokine C-X-C motif ligand (CXCL)-10 (OR = 1.404, 95% CI: 1.094-1.803, P = 0.008) on MS were observed. Nevertheless, only IL-2Rα still had a causal effect on MS after Bonferroni correction.

Conclusion

The results identify a genetically predicted causal effect of IL-2Rα, IL-8 and CXCL-10 on MS, CCL-5 and SCGF-β on CIDP.

Metformin treatment reduces the incidence of osteoporosis: a two-sample Mendelian randomized study

It remains unclear whether the association between metformin and osteoporosis (OP) risk is causal. This two-sample Mendelian randomization (MR) study suggests a causal relationship between metformin treatment and a decrease in OP and fracture incidence, as well as an increase in bone mineral density (BMD) in the lumbar spine, femoral neck, and heel. Nonetheless, no significant causal effect is observed on forearm BMD.

PURPOSE

We utilize a MR approach to investigate the association between metformin treatment and the risk of OP.

METHODS

Metformin treatment was selected as exposures. Outcomes included OP; BMD at the forearm (FA), femoral neck (FN), and lumbar spine (LS); estimated heel bone mineral density (eBMD); and fracture. Summary statistics for exposures and outcomes were obtained from corresponding genome-wide association studies. Inverse variance-weighted (IVW) analysis was mainly applied; the weighted median (WM), penalized weighted median (PWM), maximum likelihood (ML), and MR-Egger regression (MR-Egger) method were also used to obtain robust estimates. A series of sensitivity analyses including Cochran's Q test, MR-Egger regression, leave-one-out analysis, and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) were used to detect pleiotropy or heterogeneity.

RESULTS

In the main analysis, IVW estimates demonstrated that metformin treatment had a definite causal effect on the risk of OP (odds ratio (OR): 0.859, 95% CI: 0.774-0.953, P = 0.004), LS-BMD (OR: 1.063, 95% CI: 1.023-1.105, P = 0.002), FN-BMD (OR: 1.034, 95% CI: 1.000-1.069, P = 0.049), eBMD (OR: 1.035, 95% CI: 1.023-1.047, P ≤ 0.001), and fracture(OR: 0.958, 95% CI: 0.928-0.989, P = 0.008). However, it did not have an effect on FA-BMD(OR: 1.050, 95% CI: 0.969-1.138, P = 0.237).

CONCLUSIONS

This study indicated that metformin treatment is significantly associated with a reduction in the risk of OP, fracture and higher LS-BMD, FN-BMD, and eBMD. However, there was no significant association with FA-BMD.