Guidelines for performing Mendelian randomization investigations

Causal association between uric acid levels and the risk of aortic aneurysm and aortic dissection: A two-sample Mendelian randomization study

BACKGROUND AND AIMS

Increasing numbers of studies have sought to uncover the relationship between serum uric acid (UA) levels and the risk of aortic aneurysm (AA) or aortic dissection (AD), but the causality of the associations has not been established yet. To assess this potential relationship, we conducted a two-sample Mendelian randomization (MR) analysis.

METHODS AND RESULTS

We performed two-sample Mendelian randomization (MR) analysis using independent genetic variants for UA levels from a published genome-wide association study (GWAS). Summary statistics for instrument-outcome associations from FinnGen database for AA and AD. Various sensitivity analyses were performed using different MR methods including random effects inverse variance weighting, fix effects inverse variance weighting, MR-Egger, weighted median/mode, and the contamination mixture method. Genetically predicted UA levels was associated with a higher AA risk (odds ratio (OR), 1.223; 95 % confidence interval (CI), 1.058-1.388; p = 0.017) in a simple size of 209,366 individuals. No association was found between uric acid levels and the risk of AD (OR,0.812; 95 % CI, 0.423-1.200; p = 0.293).

CONCLUSION

Our study suggests a significant and robust causal association between UA levels and risk of AA but did not support such a relationship between UA levels and AD risk, which might be interpreted with caution and further confirmed. These findings support a clinically relevant causal effect of serum urate levels on the risk of AA.

Visceral and ectopic fat are more predictively associated with primary liver cancer than overall obesity from genetic sights: A Mendelian randomization study

Several observational studies have reported an association between obesity and primary liver cancer (PLC), while the causality behind this association and the comparison of the risk effects of different obesity indicators on PLC remain unclear. In this study, we performed two-sample Mendelian randomization (MR) analyses to assess the associations of genetically determined liver fat, visceral adipose tissue (VAT), and body mass index (BMI) with the risk of PLC. The summary statistics of exposures were obtained from two genome-wide association studies (GWASs) based on the UK Biobank (UKB) imaging cohort and the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. GWAS summary statistics for PLC were obtained from FinnGen consortium R7 release data, including 304 PLC cases and 218 488 controls. Inverse-variance weighted (IVW) was used as the primary analysis, and a series of sensitivity analyses were performed to further verify the robustness of these findings. IVW analysis highlighted a significant association of genetically determined liver fat (OR per SD increase: 7.14; 95% CI: 5.10-9.99; P = 2.35E-30) and VAT (OR per SD increase: 5.70; 95% CI: 1.32-24.72; P = .020) with PLC but not of BMI with PLC. The findings were further confirmed by a series of MR methods. No evidence of horizontal pleiotropy between these associations existed. Our study suggested that genetically determined liver fat and VAT rather than BMI were associated with an increased risk of PLC, which suggested that visceral fat distribution is more predictive of the clinical risk of PLC than common in vitro measures.

Causal effects of osteoporosis on structural changes in specific brain regions: a Mendelian randomization study

Observational studies have reported that osteoporosis is associated with cortical changes in the brain. However, the inherent limitations of observational studies pose challenges in eliminating confounding factors and establishing causal relationships. And previous observational studies have not reported changes in specific brain regions. By employing Mendelian randomization, we have been able to infer a causal relationship between osteoporosis and a reduction in the surficial area (SA) of the brain cortical. This effect is partially mediated by vascular calcification. We found that osteoporosis significantly decreased the SA of global brain cortical (β = -1587.62 mm2, 95%CI: -2645.94 mm2 to -529.32 mm2, P = 0.003) as well as the paracentral gyrus without global weighted (β = - 19.42 mm2, 95%CI: -28.90 mm2 to -9.95 mm2, P = 5.85 × 10-5). Furthermore, we estimated that 42.25% and 47.21% of the aforementioned effects are mediated through vascular calcification, respectively. Osteoporosis leads to a reduction in the SA of the brain cortical, suggesting the presence of the bone-brain axis. Vascular calcification plays a role in mediating this process to a certain extent. These findings establish a theoretical foundation for further investigations into the intricate interplay between bone, blood vessels, and the brain.

Assessing the genetic contribution of cumulative behavioral factors associated with longitudinal type 2 diabetes risk highlights adiposity and the brain-metabolic axis

While genetic factors, behavior, and environmental exposures form a complex web of interrelated associations in type 2 diabetes (T2D), their interaction is poorly understood. Here, using data from ~500K participants of the UK Biobank, we identify the genetic determinants of a "polyexposure risk score" (PXS) a new risk factor that consists of an accumulation of 25 associated individual-level behaviors and environmental risk factors that predict longitudinal T2D incidence. PXS-T2D had a non-zero heritability (h2 = 0.18) extensive shared genetic architecture with established clinical and biological determinants of T2D, most prominently with body mass index (genetic correlation [rg] = 0.57) and Homeostatic Model Assessment for Insulin Resistance (rg = 0.51). Genetic loci associated with PXS-T2D were enriched for expression in the brain. Biobank scale data with genetic information illuminates how complex and cumulative exposures and behaviors as a whole impact T2D risk but whose biology have been elusive in genome-wide studies of T2D.

Causal relationship between diet and knee osteoarthritis: A Mendelian randomization analysis

BACKGROUND

Knee osteoarthritis (KOA) is a common disabling joint disease that affects millions of people worldwide. Diet may play a role in the etiology and progression of KOA, but evidence for a causal relationship is limited. We aimed to investigate the causal impact of dietary intake on KOA risk using Mendelian randomization (MR).

METHODS

We used summary-level data from genome-wide association studies (GWAS) including dietary intake (n = 335, 394-462, 342), and KOA (n = 403, 124). We selected 6-77 genetic variants as instrumental variables for 18 dietary factors, including processed meat, poultry, beef, oily fish, non-oily fish, pork, lamb, frequency of alcohol intake, alcoholic beverages, tea, coffee, dried fruit, cereals, cheese, bread, cooked vegetables, salad/raw vegetables, and fresh fruit. We performed univariate and multivariate MR analyses to estimate the causal effect of each dietary factor on KOA risk. We also performed some sensitivity analyses to assess the validity of the MR hypothesis.

RESULTS

We found that higher coffee intake was associated with increased KOA risk, whereas higher intake of dried fruits, grains, cheese, and oily fish was associated with reduced KOA risk. After multivariate adjustment, we found that coffee and oily fish intake may affect KOA through obesity, body mass index (BMI), diabetes, hypertension, and prolonged standing. Sensitivity analyses did not reveal any evidence of pleiotropy.

CONCLUSIONS

Our study provides new causal evidence that dietary intake may influence KOA risk. Specifically, we suggest that increased intake of dried fruits, grains, cheese, and oily fish and decreased coffee intake may be beneficial in preventing and mitigating KOA. further studies are needed to elucidate the underlying mechanisms and to confirm our findings in different populations.

Effect of the gut microbiome in glaucoma risk from the causal perspective

OBJECTIVE

Evidence from observational studies has reported possible associations between the gut microbiome (GM) and glaucoma. However, the causal effect of GM on glaucoma risk remains to be determined.

METHODS AND ANALYSIS

We conducted two-sample bidirectional Mendelian randomisation (MR) analyses to explore the causal association between GM and glaucoma. Genome-wide association study summary statistics of 196 GM taxa (n=18 340) and glaucoma (18 902 cases and 358 375 controls) were obtained from MiBioGen and FinnGen Consortium. Inverse variance weighted, MR-Egger, weighted median, weighted mode, Mendelian Randomisation Pleiotropy Residual Sum and Outlier, MR-Egger intercept and Cochran's Q statistical analyses were used to supplement MR results and sensitivity analysis. An independent cohort from the Medical Research Council (MRC) Integrative Epidemiology Unit at the University of Bristol (MRC-IEU) Consortium (1715 cases and 359 479 controls) was used to validate causal effects.

RESULTS

Results of the MR analysis suggested that the family Oxalobacteraceae (OR 0.900, 95% CI 0.843 to 0.961, p=0.002) and the genus Eggerthella (OR 0.881, 95% CI 0.811 to 0.957, p=0.003) had a negative effect on glaucoma, whereas the genus Bilophila (OR 1.202, 95% CI 1.074 to 1.346, p=0.001), LachnospiraceaeUCG010 (OR 1.256, 95% CI 1.109 to 1.423, p=0.0003) and Ruminiclostridium 9 (OR 1.258, 95% CI 1.083 to 1.461, p=0.003) had a positive effect on glaucoma. Among these, the positive causal effect of LachnospiraceaeUCG010 (OR 1.002, 95% CI 1.000 to 1.004, p=0.033) on glaucoma was replicated in an independent cohort.

CONCLUSION

This MR analysis from large population studies demonstrated the causal effect of GM on glaucoma risk and supported the role of GM in influencing glaucoma susceptibility.

Interaction between systemic iron parameters and left ventricular structure and function in the preserved ejection fraction population: a two-sample bidirectional Mendelian randomization study

BACKGROUND

Left ventricular (LV) remodeling and diastolic function in people with heart failure (HF) are correlated with iron status; however, the causality is uncertain. This Mendelian randomization (MR) study investigated the bidirectional causal relationship between systemic iron parameters and LV structure and function in a preserved ejection fraction population.

METHODS

Transferrin saturation (TSAT), total iron binding capacity (TIBC), and serum iron and ferritin levels were extracted as instrumental variables for iron parameters from meta-analyses of public genome-wide association studies. Individuals without myocardial infarction history, HF, or LV ejection fraction (LVEF) < 50% (n = 16,923) in the UK Biobank Cardiovascular Magnetic Resonance Imaging Study constituted the outcome dataset. The dataset included LV end-diastolic volume, LV end-systolic volume, LV mass (LVM), and LVM-to-end-diastolic volume ratio (LVMVR). We used a two-sample bidirectional MR study with inverse variance weighting (IVW) as the primary analysis method and estimation methods using different algorithms to improve the robustness of the results.

RESULTS

In the IVW analysis, one standard deviation (SD) increased in TSAT significantly correlated with decreased LVMVR (β = -0.1365; 95% confidence interval [CI]: -0.2092 to -0.0638; P = 0.0002) after Bonferroni adjustment. Conversely, no significant relationships were observed between other iron and LV parameters. After Bonferroni correction, reverse MR analysis showed that one SD increase in LVEF significantly correlated with decreased TSAT (β = -0.0699; 95% CI: -0.1087 to -0.0311; P = 0.0004). No heterogeneity or pleiotropic effects evidence was observed in the analysis.

CONCLUSIONS

We demonstrated a causal relationship between TSAT and LV remodeling and function in a preserved ejection fraction population.

No causal association between allergic rhinitis and migraine: a Mendelian randomization study

PURPOSE

Allergic rhinitis (AR) and migraine are among the most common public health problems worldwide. Observational studies on the correlation between AR and migraine have reported inconsistent results. This study aimed to investigate the causal relationship of AR with migraine and its subtypes, including migraine with aura (MA) and migraine without aura (MO).

METHODS

Bidirectional two-sample Mendelian randomization (MR) analysis was performed with publicly available summary-level statistics of large genome-wide association studies to estimate the possible causal effects. The inverse variance-weighted method was selected for primary analysis and was supplemented with the weighted median, weighted mode, and MR-Egger methods. The causal analysis using summary effect estimates (CAUSE) were further performed to verify the causality. Several sensitivity tests, including the leave-one-out, Cochran's Q, MR-Egger intercept, and MR-PRESSO tests, were performed to assess the robustness of the results.

RESULTS

AR did not exhibit a significant causal correlation with the elevated risk of any migraine (odd ratio (OR), 0.816; 95% confidence interval (CI), 0.511-1.302; P = 0.394), MA (OR, 0.690; 95% CI 0.298-1.593; P = 0.384), or MO (OR, 1.022; 95% CI 0.490-2.131; P = 0.954). Consistently, reverse MR analysis did not reveal causal effects of any migraine or its subtypes on AR. Almost all sensitivity analyses supported the robustness of the results.

CONCLUSIONS

This MR study did not reveal a clear causal association between AR and migraine risk. More research is warranted to reveal the complex association between AR and migraine.

Parental genetically predicted liability for coronary heart disease and risk of adverse pregnancy outcomes: a cohort study

BACKGROUND

Adverse pregnancy outcomes (APO) may unmask or exacerbate a woman's underlying risk for coronary heart disease (CHD). We estimated associations of maternal and paternal genetically predicted liability for CHD with lifelong risk of APOs. We hypothesized that associations would be found for women, but not their male partners (negative controls).

METHODS

We studied up to 83,969‬ women (and up to 55,568‬ male partners) from the Norwegian Mother, Father and Child Cohort Study or the Trøndelag Health Study with genotyping data and lifetime history of any APO in their pregnancies (1967-2019) in the Medical Birth Registry of Norway (miscarriage, stillbirth, hypertensive disorders of pregnancy, gestational diabetes, small for gestational age, large for gestational age, and spontaneous preterm birth). Maternal and paternal genetic risk scores (GRS) for CHD were generated using 148 gene variants (p-value < 5 × 10-8, not in linkage disequilibrium). Associations between GRS for CHD and each APO were determined using logistic regression, adjusting for genomic principal components, in each cohort separately, and combined using fixed effects meta-analysis.

RESULTS

One standard deviation higher GRS for CHD in women was related to increased risk of any hypertensive disorders of pregnancy (odds ratio [OR] 1.08, 95% confidence interval [CI] 1.05-1.10), pre-eclampsia (OR 1.08, 95% CI 1.05-1.11), and small for gestational age (OR 1.04, 95% CI 1.01-1.06). Imprecise associations with lower odds of large for gestational age (OR 0.98, 95% CI 0.96-1.00) and higher odds of stillbirth (OR 1.04, 95% CI 0.98-1.11) were suggested. These findings remained consistent after adjusting for number of total pregnancies and the male partners' GRS and restricting analyses to stable couples. Associations for other APOs were close to the null. There was weak evidence of an association of paternal genetically predicted liability for CHD with spontaneous preterm birth in female partners (OR 1.02, 95% CI 0.99-1.05), but not with other APOs.

CONCLUSIONS

Hypertensive disorders of pregnancy, small for gestational age, and stillbirth may unmask women with a genetically predicted propensity for CHD. The association of paternal genetically predicted CHD risk with spontaneous preterm birth in female partners needs further exploration.

Blood DNA methylation signature of diet quality and association with cardiometabolic traits

AIMS

Diet quality might influence cardiometabolic health through epigenetic changes, but this has been little investigated in adults. Our aims were to identify cytosine-phosphate-guanine (CpG) dinucleotides associated with diet quality by conducting an epigenome-wide association study (EWAS) based on blood DNA methylation (DNAm) and to assess how diet-related CpGs associate with inherited susceptibility to cardiometabolic traits: body mass index (BMI), systolic blood pressure (SBP), triglycerides, type 2 diabetes (T2D), and coronary heart disease (CHD).

METHODS AND RESULTS

Meta-EWAS including 5274 participants in four cohorts from Spain, the USA, and the UK. We derived three dietary scores (exposures) to measure adherence to a Mediterranean diet, to a healthy plant-based diet, and to the Dietary Approaches to Stop Hypertension. Blood DNAm (outcome) was assessed with the Infinium arrays Human Methylation 450K BeadChip and MethylationEPIC BeadChip. For each diet score, we performed linear EWAS adjusted for age, sex, blood cells, smoking and technical variables, and BMI in a second set of models. We also conducted Mendelian randomization analyses to assess the potential causal relationship between diet-related CpGs and cardiometabolic traits. We found 18 differentially methylated CpGs associated with dietary scores (P < 1.08 × 10-7; Bonferroni correction), of which 12 were previously associated with cardiometabolic traits. Enrichment analysis revealed overrepresentation of diet-associated genes in pathways involved in inflammation and cardiovascular disease. Mendelian randomization analyses suggested that genetically determined methylation levels corresponding to lower diet quality at cg02079413 (SNORA54), cg02107842 (MAST4), and cg23761815 (SLC29A3) were causally associated with higher BMI and at cg05399785 (WDR8) with greater SBP, and methylation levels associated with higher diet quality at cg00711496 (PRMT1) with lower BMI, T2D risk, and CHD risk and at cg0557921 (AHRR) with lower CHD risk.

CONCLUSION

Diet quality in adults was related to differential methylation in blood at 18 CpGs, some of which related to cardiometabolic health.

SARS-CoV-2 neurovascular invasion supported by Mendelian randomization

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is known to affect vessels and nerves and can be easily visualized in the retina. However, the effect of SARS-CoV-2 on retinal morphology remains controversial. In the present research, we applied Mendelian randomization (MR) analysis to estimate the association between SARS-CoV-2 and changes in the thickness of the inner retina.

METHODS

Two-sample MR analysis was conducted using summary-level data from 3 open genome-wide association study databases concerning COVID-19 infection (2,942,817 participants) and COVID-19 hospitalization (2,401,372 participants); moreover, the dataset of inner retina thickness, including the macular retinal nerve fiber layer (mRNFL) and macular ganglion cell-inner plexiform layer (mGCIPL), included 31,434 optical coherence tomography (OCT) images derived from healthy UK Biobank participants. All the participants were of European ancestry. The inverse variance weighted (IVW) meta-analysis was used as our primary method. Various complementary MR approaches were established to provide robust causal estimates under different assumptions.

RESULTS

According to our MR analysis, genetically predicted COVID-19 infection was associated with an increased risk of mRNFL and mGCIPL thickness (OR = 1.74, 95% CI 1.20-2.52, P = 3.58 × 10-3; OR = 2.43, 95% CI 1.49-3.96, P = 3.6 × 10-4). The other MR methods produced consistent results. However, genetically predicted COVID-19 hospitalization did not affect the thickness of the inner retina (OR = 1.11, 95% CI 0.90-1.37, P = 0.32; OR = 1.28, 95% CI 0.88-1.85, P = 0.19).

CONCLUSION

This work provides the first genetically predictive causal evidence between COVID-19 infection and inner retinal thickness in a European population. These findings will contribute to further understanding of the pathogenesis of COVID-19 and stimulate improvements in treatment modalities.

A Mendelian analysis of the relationships between immune cells and breast cancer

Background

Emerging evidence showed immune cells were associated with the development of breast cancer. Nonetheless, the causal link between them remains uncertain. Consequently, the objective of this study was to investigate the causal connection between immune traits and the likelihood of developing breast cancer.

Methods

A two-sample Mendelian randomization (MR) analysis was conducted to establish the causal relationship between immune cells and breast cancer in this study. Utilizing publicly accessible genetic data, we investigated causal connections between 731 immune cells and the occurrence of breast cancer. The primary approach for exploring this relationship was the application of the inverse-variance-weighted (IVW) method. Furthermore, sensitivity analyses, encompassing the leave-one-out analysis, Cochran Q test, and Egger intercept test were performed to validate the reliability of the Mendelian randomization results. Finally, we used Bayesian Weighted Mendelian Randomization (BWMR) approach to test the results of MR study.

Results

According to the Bonferroni correction, no immune trait was identified with a decreased or increased risk of overall breast cancer risk. As for the ER+ breast cancer, 6 immune trait was identified after the Bonferroni method. the IVW method results showed that CD45RA- CD4+ %CD4+ (p-value:1.37×10-6), CD8dim %T cell (p-value:4.62×10-43), BAFF-R on IgD+ CD38- unsw mem (p-value:6.93×10-5), CD27 on PB/PC (p-value:2.72×10-18) lowered the risk of breast cancer. However, CD19 on IgD- CD38br (p-value:1.64×10-6), CD25 on IgD+ CD38dim (p-value: - ∞) were associated with a higher risk of developing breast cancer. As for the CX3CR1 on CD14+ CD16- monocyte (p-value: 1.15×10-166), the IVW method clearly demonstrated a protective effect against ER- breast cancer. For the above positive results, BAFF-R on IgD+ CD38- unsw mem was the sole association linked to reduced breast cancer risk using the BWMR method. The intercept terms' p-values in MR-Egger regression all exceeded 0.05, indicating the absence of potential horizontal pleiotropy.

Conclusion

Through genetic approaches, our study has illustrated the distinct correlation between immune cells and breast cancer, potentially paving the way for earlier diagnosis and more efficient treatment alternatives.

Genetic determinants of complement activation in the general population

Complement is a fundamental innate immune response component. Its alterations are associated with severe systemic diseases. To illuminate the complement's genetic underpinnings, we conduct genome-wide association studies of the functional activity of the classical (CP), lectin (LP), and alternative (AP) complement pathways in the Cooperative Health Research in South Tyrol study (n = 4,990). We identify seven loci, encompassing 13 independent, pathway-specific variants located in or near complement genes (CFHR4, C7, C2, MBL2) and non-complement genes (PDE3A, TNXB, ABO), explaining up to 74% of complement pathways' genetic heritability and implicating long-range haplotypes associated with LP at MBL2. Two-sample Mendelian randomization analyses, supported by transcriptome- and proteome-wide colocalization, confirm known causal pathways, establish within-complement feedback loops, and implicate causality of ABO on LP and of CFHR2 and C7 on AP. LP causally influences collectin-11 and KAAG1 levels and the risk of mouth ulcers. These results build a comprehensive resource to investigate the role of complement in human health.

Autoimmune diseases and risk of non-Hodgkin lymphoma: A Mendelian randomisation study

Objective

To examine whether genetically predicted susceptibility to ten autoimmune diseases (Behçet's disease, coeliac disease, dermatitis herpetiformis, lupus, psoriasis, rheumatoid arthritis, sarcoidosis, Sjögren's syndrome, systemic sclerosis, and type 1 diabetes) is associated with risk of non-Hodgkin lymphoma (NHL).

Design

Two sample Mendelian randomization (MR) study.

Setting

Genome wide association studies (GWASs) of ten autoimmune diseases, NHL, and four NHL subtypes (i.e., follicular lymphoma, mature T/natural killer-cell lymphomas, non-follicular lymphoma, and other and unspecified types of NHL).

Analysis

We used data from the largest publicly available GWASs of European ancestry for each autoimmune disease, NHL, and NHL subtypes. For each autoimmune disease, we extracted single nucleotide polymorphisms (SNPs) strongly associated (P < 5×10-8) with that disease and that were independent of one another (R2 < 1×10-3) as genetic instruments. SNPs within the human leukocyte antigen region were not considered due to potential pleiotropy. Our primary MR analysis was the inverse-variance weighted analysis. Additionally, we conducted MR-Egger, weighted mode, and weighted median regression to address potential bias due to pleiotropy, and robust adjusted profile scores to address weak instrument bias. We carried out sensitivity analysis limited to the non-immune pathway for nominally significant findings. To account for multiple testing, we set the thresholds for statistical significance at P < 5×10-3.

Participants

The number of cases and controls identified in the relevant GWASs were 437 and 3,325 for Behçet's disease, 4,918 and 5,684 for coeliac disease, 435 and 341,188 for dermatitis herpetiformis, 4,576 and 8,039 for lupus, 11,988 and 275,335 for psoriasis, 22,350 and 74,823 for rheumatoid arthritis, 3,597 and 337,121 for sarcoidosis, 2,735 and 332,115 for Sjögren's syndrome, 9,095 and 17,584 for systemic sclerosis, 18,942 and 501,638 for type 1 diabetes, 2,400 and 410,350 for NHL; and 296 to 2,340 cases and 271,463 controls for NHL subtypes.

Exposures

Genetic variants predicting ten autoimmune diseases: Behçet's disease, coeliac disease, dermatitis herpetiformis, lupus, psoriasis, rheumatoid arthritis, sarcoidosis, Sjögren's syndrome, systemic sclerosis, and type 1 diabetes.

Main outcome measures

Estimated associations between genetically predicted susceptibility to ten autoimmune diseases and the risk of NHL.

Results

The variance of each autoimmune disease explained by the SNPs ranged from 0.3% to 3.1%. Negative associations between type 1 diabetes and sarcoidosis and the risk of NHL were observed (odds ratio [OR] 0.95, 95% confidence interval [CI]: 0.92 to 0.98, P = 5×10-3, and OR 0.92, 95% CI: 0.85 to 0.99, P = 2.8×10-2, respectively). These findings were supported by the sensitivity analyses accounting for potential pleiotropy and weak instrument bias. No significant associations were found between the other eight autoimmune diseases and NHL risk. Of the NHL subtypes, type 1 diabetes was most strongly associated with follicular lymphoma (OR 0.91, 95% CI: 0.86 to 0.96, P = 1×10-3), while sarcoidosis was most strongly associated with other and unspecified NHL (OR 0.86, 95% CI: 0.75 to 0.97, P = 1.8×10-2).

Conclusions

These findings suggest that genetically predicted susceptibility to type 1 diabetes, and to some extent sarcoidosis, might reduce the risk of NHL. However, future studies with different datasets, approaches, and populations are warranted to further examine the potential associations between these autoimmune diseases and the risk of NHL.

Causal relationship between multiple sclerosis and cortical structure: a Mendelian randomization study

BACKGROUND

Observational studies have suggested an association between multiple sclerosis (MS) and cortical structure, but the results have been inconsistent.

OBJECTIVE

We used two-sample Mendelian randomization (MR) to assess the causal relationship between MS and cortical structure.

METHODS

MS data as the exposure trait, including 14,498 cases and 24,091 controls, were obtained from the International Multiple Sclerosis Genetics Consortium. Genome-wide association study (GWAS) data for cortical surface area (SAw/nw) and thickness (THw/nw) in 51,665 individuals of European ancestry were obtained from the ENIGMA Consortium. The inverse-variance weighted (IVW) method was used as the primary analysis for MR. Sensitivity analyses were conducted to evaluate heterogeneity and pleiotropy. Enrichment analysis was performed on MR analyses filtered by sensitivity analysis.

RESULTS

After IVW and sensitivity analysis filtering, only six surviving MR results provided suggestive evidence supporting a causal relationship between MS and cortical structure, including lingual SAw (p = .0342, beta (se) = 5.7127 (2.6969)), parahippocampal SAw (p = .0224, beta (se) = 1.5577 (0.6822)), rostral middle frontal SAw (p = .0154, beta (se) = - 9.0301 (3.7281)), cuneus THw (p = .0418, beta (se) =  - 0.0020 (0.0010)), lateral orbitofrontal THw (p = .0281, beta (se) = 0.0025 (0.0010)), and lateral orbitofrontal THnw (p = .0417, beta (se) = 0.0029 (0.0014)). Enrichment analysis suggested that leukocyte cell-related pathways, JAK-STAT signaling pathway, NF-kappa B signaling pathway, cytokine-cytokine receptor interaction, and prolactin signaling pathway may be involved in the effect of MS on cortical morphology.

CONCLUSION

Our results provide evidence supporting a causal relationship between MS and cortical structure. Enrichment analysis suggests that the pathways mediating brain morphology abnormalities in MS patients are mainly related to immune and inflammation-driven pathways.

Large-Scale Mendelian Randomization Study Reveals Circulating Blood-based Proteomic Biomarkers for Psychopathology and Cognitive Task Performance

Background

Research on peripheral (e.g., blood-based) biomarkers for psychiatric illness has typically been low-throughput in terms of both the number of subjects and the range of assays performed. Moreover, traditional case-control studies examining blood-based biomarkers are subject to potential confounds of treatment and other exposures common to patients with psychiatric illnesses. Our research addresses these challenges by leveraging large-scale, high-throughput proteomics data and Mendelian Randomization (MR) to examine the causal impact of circulating proteins on psychiatric phenotypes and cognitive task performance.

Methods

We utilized plasma proteomics data from the UK Biobank (3,072 proteins assayed in 34,557 European-ancestry individuals) and deCODE Genetics (4,719 proteins measured across 35,559 Icelandic individuals). Significant proteomic quantitative trait loci (both cis-pQTLs and trans-pQTLs) served as MR instruments, with the most recent GWAS for schizophrenia, bipolar disorder, major depressive disorder, and cognitive task performance (all excluding overlapping UK Biobank participants) as phenotypic outcomes.

Results

MR revealed 109 Bonferroni-corrected causal associations (44 novel) involving 88 proteins across the four phenotypes. Several immune-related proteins, including interleukins and complement factors, stood out as pleiotropic across multiple outcome phenotypes. Drug target enrichment analysis identified several novel potential pharmacologic repurposing opportunities, including anti-inflammatory agents for schizophrenia and bipolar disorder and duloxetine for cognitive performance.

Conclusions

Identification of causal effects for these circulating proteins suggests potential biomarkers for these conditions and offers insights for developing innovative therapeutic strategies. The findings also indicate substantial evidence for the pleiotropic effects of many proteins across different phenotypes, shedding light on the shared etiology among psychiatric conditions and cognitive ability.

The Impact of Different Intensities of Physical Activity on Serum Urate and Gout: A Mendelian Randomization Study

Physical activity is a potential protective factor against gout, but the role of exercise intensity in this context remains unclear. To overcome the limitations of observational studies in causal inference, this study employed a two-sample Mendelian randomization approach to explore the impact of different genetically proxied/predicted intensities of physical activity on serum urate concentration and the incidence of gout. Our data related to physical activity, serum urate, and gout were obtained from the UK Biobank, the Global Urate Genetics Consortium (GUGC), and the FinnGen dataset, respectively. Walking was included as representative of typical low-intensity physical activity in the analysis, and the other two types were moderate and vigorous physical activities. The estimation methods we used included the inverse-variance-weighted (IVW) method, MR-Egger regression, weighted-median method, simple-mode method, and weighted-mode method. Sensitivity analyses involved Rucker's framework, Cochran's Q test, funnel plots, MR-PRESSO outlier correction, and leave-one-out analysis. We found suggestive evidence from the inverse-variance-weighted method that moderate physical activity was a potential factor in reducing the incidence of gout (OR = 0.628, p = 0.034), and this association became more substantial in our subsequent sensitivity analysis (OR = 0.555, p = 0.006). However, we observed no distinctive effects of physical activity on serum urate concentration. In conclusion, our study supports some findings from observational studies and emphasizes the preventive role of moderate physical activity against gout. Given the limitations of the existing datasets, we call for future reexamination and expansion of our findings using new GWAS data.

The causal relationship between immune cells and ankylosing spondylitis: a bidirectional Mendelian randomization study

BACKGROUND

Ankylosing spondylitis (AS) is one of several disorders known as seronegative spinal arthritis (SpA), the origin of which is unknown. Existing epidemiological data show that inflammatory and immunological factors are important in the development of AS. Previous research on the connection between immunological inflammation and AS, however, has shown inconclusive results.

METHODS

To evaluate the causal association between immunological characteristics and AS, a bidirectional, two-sample Mendelian randomization (MR) approach was performed in this study. We investigated the causal connection between 731 immunological feature characteristic cells and AS risk using large, publically available genome-wide association studies.

RESULTS

After FDR correction, two immunophenotypes were found to be significantly associated with AS risk: CD14 - CD16 + monocyte (OR, 0.669; 95% CI, 0.544 ~ 0.823; P = 1.46 × 10-4; PFDR = 0.043), CD33dim HLA DR + CD11b + (OR, 0.589; 95% CI = 0.446 ~ 0.780; P = 2.12 × 10-4; PFDR = 0.043). AS had statistically significant effects on six immune traits: CD8 on HLA DR + CD8 + T cell (OR, 1.029; 95% CI, 1.015 ~ 1.043; P = 4.46 × 10-5; PFDR = 0.014), IgD on IgD + CD24 + B cell (OR, 0.973; 95% CI, 0.960 ~ 0.987; P = 1.2 × 10-4; PFDR = 0.021), IgD on IgD + CD38 - unswitched memory B cell (OR, 0.962; 95% CI, 0.945 ~ 0.980; P = 3.02 × 10-5; PFDR = 0.014), CD8 + natural killer T %lymphocyte (OR, 0.973; 95% CI, 0.959 ~ 0.987; P = 1.92 × 10-4; PFDR = 0.021), CD8 + natural killer T %T cell (OR, 0.973; 95% CI, 0.959 ~ 0.987; P = 1.65 × 10-4; PFDR = 0.021).

CONCLUSION

Our findings extend genetic research into the intimate link between immune cells and AS, which can help guide future clinical and basic research.

The peripheral immune cell counts and mouth ulcers: A two-sample Mendelian randomization study

Objective

This study explored the causal association of peripheral immune cell counts with mouth ulcers (MUs) by two-sample Mendelian Randomization.

Design

The counts of 12 circulating immune cell types (leukocytes, lymphocytes, monocytes, eosinophils, neutrophils, basophils, CD4+ cells, CD8+ cells, unswitched memory B cells, NK cells, B cells and a derived ratio (CD4+/CD8+)) were determined as the exposure. MUs were the outcome. The analysis was conducted mostly using the inverse-variance weighted (IVW) approach. MR Egger, weighted median, weighted mode and simple mode were used to detect the horizontal pleiotropy.

Results

The IVW results for leukocytes and lymphocyte counts were OR = 0.93, 95 % CI = 0.88-0.98, p = 0.0115 and OR = 0.91, 95 % CI: 0.84-0.98, p = 0.0150, respectively. The Wald ratio result for CD4+ cell and CD8+ cell counts were OR = 0.70, 95 % CI: 0.65-0.75, p = 1.05 × 10-20 and OR = 1.25, 95 % CI: 1.19-1.31, p = 9.99 × 10-21, respectively.

Conclusions

This study supports a causal effect of peripheral immune cell counts on MUs. Higher leukocyte, lymphocyte and CD4+ cell counts can protect against MUs, but higher CD8+ cell counts enhance the risk of MUs. This finding confirms host immune factors play a crucial role in the aetiology of MUs.

Elucidating the relationship between breast cancer and brain cortical structure: a Mendelian randomization study

Cancer-associated cognitive impairment is a significant challenge for individuals who have survived breast cancer, affecting their quality of life. In this study, we conducted an inaugural comprehensive Mendelian randomization analysis discerning the causal relationship between breast cancer, including its two subtypes, and the cerebral cortical structure. Our analysis indicated that estrogen receptor-negative breast cancer significantly decreased surface area (β = -593.01 mm2, 95% CI: -1134.9 to -51.1 mm2, P = 0.032). At the regional level, estrogen receptor-negative breast cancer showed a significant association with surface area and thickness in 17 cortical regions. These regions included the insula, posterior cingulate, superior frontal, precuneus, fusiform, lateral occipital, and rostral middle frontal. Specifically, estrogen receptor-negative breast cancer had a significant impact on decreasing the surface area of the insula without considering global weight (β = -14.09 mm2, 95% CI: -22.91 to -5.27 mm2, P = 0.0017). The results from meta-analysis and LD Score Regression provide support for our findings. This investigation unveils the correlations between breast cancer, its various subcategories, and the cerebral cortical structure. Notably, breast cancer of the estrogen receptor-negative variety may elicit more widespread cerebral atrophy.

Genetically predicted chronic rhinosinusitis and the risk of stroke: a two-sample Mendelian randomization study

Objective

The causal association between chronic rhinosinusitis (CRS) and stroke remains uncertain due to the susceptibility of observational studies to confounding and the possibility of reverse causality. This study aims to examine the potential causal relationship between CRS and the risk of stroke, encompassing various subtypes.

Methods

In this research, we utilized genome-wide association study (GWAS) data for CRS from FinnGen. We identified significant single-nucleotide polymorphisms (SNPs) associated with CRS and used them as instrumental variables (IVs). GWAS data for any ischemic stroke (AIS), ischemic stroke (IS), large-artery atherosclerotic stroke (LAS), small-vessel strokes (SVS), cardioembolic strokes (CES), intracerebral hemorrhage (ICH), lobar ICH, and non-lobar ICH came from multi-ancestry GWAS datasets. We conducted two-sample Mendelian randomization (MR) analyses using inverse variance weighting (IVW), weighted median, and MR-Egger regression methods to investigate potential causal relationships between CRS and stroke. Both heterogeneity and pleiotropy were evaluated by sensitivity analyses.

Result

The IVW analysis revealed no significant associations between CRS and AIS (OR = 0.99, 95% CI [0.93-1.05], p = 0.73), IS (OR = 0.97, 95% CI [0.81-1.17], p = 0.09), SVS (OR = 0.96, 95% CI [0.82-1.12], p = 0.58), LAS (OR = 0.91, 95% CI [0.77-1.08], p = 0.09), CES (OR = 0.97, 95% CI [0.81-1.17], p = 0.79), ICH (OR = 1.28, 95% CI [0.74-2.22], p = 0.28), lobar ICH (OR = 1.22, 95% CI [0.60-2.50], p = 0.28), and non-lobar ICH (OR = 1.25, 95% CI [0.65-2.40], p = 0.79). Sensitivity analysis found no evidence of horizontal pleiotropy.

Conclusion

According to genetic evidence, this Mendelian randomization (MR) study does not indicate a causal relationship between CRS and stroke in European populations. However, further studies are necessary to comprehensively evaluate the potential association between CRS and stroke.

Potential disease biomarkers for diabetic retinopathy identified through Mendelian randomization analysis

Background

Diabetic retinopathy (DR), a leading cause of vision loss, has limited options for effective prevention and treatment. This study aims to utilize genomics and proteomics data to identify potential drug targets for DR.

Methods

We utilized plasma protein quantitative trait loci data from the Atherosclerosis Risk in Communities Study and the Icelandic Decoding Genetics Study for discovery and replication, respectively. Genetic associations with DR, including its subtypes, were derived from the FinnGen study. Mendelian Randomization (MR) analysis estimated associations between protein levels and DR risk, complemented by colocalization analysis to examine shared causal variants.

Results

Our MR analysis identified significant associations of specific plasma proteins with DR and proliferative DR (PDR). Elevated genetically predicted levels of WARS (OR = 1.16; 95% CI = 0.095-0.208, FDR = 1.31×10-4) and SIRPG (OR = 1.15; 95% CI = 0.071-0.201, FDR = 1.46×10-2) were associated with higher DR risk, while increased levels of ALDOC (OR = 1.56; 95% CI = 0.246-0.637, FDR = 5.48×10-3) and SIRPG (OR = 1.15; 95% CI = 0.068-0.208, FDR = 4.73×10-2) were associated with higher PDR risk. These findings were corroborated by strong colocalization evidence.

Conclusions

Our study highlights WARS, SIRPG, and ALDOC as significant proteins associated with DR and PDR, providing a basis for further exploration in drug development. Additional studies are needed to validate these proteins as disease biomarkers across diverse populations.

Visceral adipose tissue and osteoarthritis, a two-sample Mendelian randomized study

Background

The relationship between visceral adipose tissue and osteoarthritis is not clear. The purpose of our study was to explore the relationship between visceral adipose tissue and osteoarthritis.

Methods

We used a two-sample Mendelian randomization method to select single-nucleotide polymorphisms (SNPs) significantly associated with visceral adipose tissue as instrumental variables to explore the relationship between visceral adipose tissue and all osteoarthritis, hand osteoarthritis, hip osteoarthritis, knee osteoarthritis, and spine osteoarthritis. The reliability of the results was tested using sensitivity analysis.

Results

Our findings indicated that visceral adipose tissue was associated with all osteoarthritis, hip osteoarthritis, knee osteoarthritis, and spine osteoarthritis (all osteoarthritis: OR = 1.399, 95% CI: 1.335-1.467, p = 7.95e-44; hip osteoarthritis: OR = 1.399, 95% CI: 1.284-1.524, p = 1.41e-14; knee osteoarthritis: OR = 1.794, 95% CI: 1.662-1.937, p = 1.33e-50; and spine osteoarthritis: OR = 1.445, 95% CI: 1.314-1.589, p = 2.89e-14). Sensitivity analysis demonstrated the reliability of these results.

Conclusion

Our study suggests that genetically predicted visceral adipose tissue is associated with osteoarthritis. Reducing the accumulation of visceral adipose tissue could potentially have an impact on the incidence of osteoarthritis.

Exploring the intestinal ecosystem: from gut microbiota to associations with subtypes of inflammatory bowel disease

Objective

Significant differences have been discovered between subtypes of Crohn's disease (CD) and ulcerative colitis (UC). The role of gut microbiota in promoting the onset of UC and CD is established, but conclusions regarding subtype-specific analyses remain limited.

Methods

This study aims to explore the influence of gut microbiota on subtypes of UC and CD, offering novel insights into the pathogenesis and treatment of UC and CD.Two-sample Mendelian randomization (MR) analysis was employed to examine the causal relationship between subtypes of UC and CD and gut microbiota composition. Gut microbiota data were sourced from the International Consortium MiBioGen, while UC and CD data were obtained from FINNGEN. Eligible single nucleotide polymorphisms (SNPs) were selected as instrumental variables. Multiple analytical approaches such as inverse variance-weighted (IVW), MR-Egger regression, weighted median, weighted mode, and MR-RAPS were utilized. Sensitivity analyses including MR-Egger intercept test, Cochran's Q test, and leave-one-out analysis were conducted for quality control. Subsequently, we employed multivariable IVW, MR-Egger, weighted median, and LASSO regression methods to identify independently significant genera or families and conducted sensitivity analyses.

Results

We have determined that Hungatella, Acidaminococcaceae, and 15 other microbial taxa act as protective factors for various CD and UC subtypes, while Terrisporobacter, Anaerostipes, and 23 other microbial taxa are associated with increased risk for different CD and UC subtypes. Furthermore, through multivariable MR analysis, we have identified significant genera or families with independent effects.

Conclusion

Our study confirms a causal relationship between dysbiosis of gut microbiota and the occurrence of CD and UC subtypes. Furthermore, it validates etiological distinctions among different subtypes of CD and UC. A novel approach to adjunctive therapy involving distinct UC or CD subtypes may involve the use of probiotics and represents a potential avenue for future treatments.

Assessment of bidirectional relationships between multiple sclerosis and epilepsy: A two-sample Mendelian randomization study

BACKGROUND AND OBJECTIVE

Epidemiological studies indicate that multiple sclerosis (MS) is associated with epilepsy. However, the causality and directionality of this association remain under-elucidated. This study aimed to reveal the causality between MS and epilepsy.

METHODS

A two-sample Mendelian randomization (MR) analysis was performed by using summarized statistics derived from large genome-wide association studies of MS and epilepsy. We used the inverse variance weighted method as the primary approach, and then four other MR methods to bidirectionally evaluate the causality of the association between MS and epilepsy. Additional sensitivity analyses were performed to measure the robustness of the findings.

RESULTS

Genetically predicted MS was positively correlated with developing all epilepsy [odds ratio (OR) = 1.027 (1.003-1.051), P  =  0.028] and generalized epilepsy [OR = 1.050 (1.008-1.094), P = 0.019]. In the reverse MR analysis, all epilepsy [OR = 1.310 (1.112-1.543), P = 0.001], generalized epilepsy [OR = 1.173 (1.010-1.363), P = 0.037], and focal epilepsy [OR = 1.264 (1.069-1.494), P  =  0.006] elevated the risk of developing MS. The result remained robust and congruous across all sensitivity analyses conducted.

CONCLUSIONS

MS is potentially associated with a higher risk of developing epilepsy. Furthermore, epilepsy may be a causal determinant of MS risk. These findings may further the understanding of the interaction of the two conditions.

Potential Biomarkers in Cerebrospinal Fluid and Plasma for Dementia

Background

The identification of biomarkers for different dementias in plasma and cerebrospinal fluid (CSF) has made substantial progress. However, they are observational studies, and there remains a lack of research on dementias with low incidence rates.

Objective

We performed a comprehensive Mendelian randomization to identify potential biomarkers for different dementia type.

Methods

The summary-level datasets encompassed 734 plasma and 154 cerebrospinal fluid proteins sourced from recently published genome-wide association studies (GWAS). Summary statistics for different dementias, including any dementia (refering to any type of dementia symptoms, 218,792 samples), Alzheimer's disease (AD, 63,926 samples), vascular dementia (212,389 samples), frontotemporal dementia (3,024 samples), dementia with Lewy bodies (DLB, 6,618 samples), and dementia in Parkinson's disease (216,895 samples), were collected from large GWAS. The primary method is inverse variance weighting, with additional sensitivity analyses conducted to ensure the robustness of the findings.

Results

The molecules released into CSF, namely APOE2 for any dementia, APOE2 and Siglec-3 for AD, APOE2 for vascular dementia, and APOE2 for DLB, might be potential biomarkers. CD33 for AD and SNCA for DLB in plasma could be promising biomarkers.

Conclusions

This is the first study to integrate plasma and CSF proteins to identify potential biomarkers for different dementias.

Primary sclerosing cholangitis causally affects kidney function decline: A Mendelian randomization study

BACKGROUND AND AIM

The causal linkage between primary sclerosing cholangitis (PSC) and kidney function is unexplored despite their potential for long-term detrimental effects on kidney function.

METHODS

Two-sample summary-level Mendelian randomization (MR) study was conducted to identify the association between PSC and kidney function. The genetic variants were extracted from the PSC-specific multi-trait analyzed genome-wide association study (GWAS) of European ancestry. Summary-level data for kidney function traits, including estimated glomerular filtration rate (eGFR), annual eGFR decline, and chronic kidney disease (CKD), were obtained from the CKDGen consortium. Multiplicative random-effects inverse-variance weighted (MR-IVW), and a series of pleiotropy-robust analyses were performed to investigate the causal effects and ascertain their robustness.

RESULTS

Significant causal associations between genetically predicted PSC and kidney function traits were identified. Genetically predicted PSC was associated with decreased log-transformed eGFR (MR-IVW; beta = -0.41%; standard error [SE] = 0.02%; P < 0.001), increased rate of annual eGFR decline (MR-IVW; beta = 2.43%; SE = 0.18%; P < 0.001), and higher risk of CKD (MR-IVW; odds ratio = 1.07; 95% confidence interval = 1.06-1.08; P < 0.001). The main findings were supported by pleiotropy-robust analysis, including MR-Egger with bootstrapped error and weighted median.

CONCLUSIONS

Our study demonstrates that genetically predicted PSC is causally associated with kidney function impairment. Further studies are warranted to identify the underlying mechanisms.

Exploring the impact of computer game playing on cognitive function, Alzheimer's disease risk, and brain-derived neurotrophic factor levels: Basic evidence from Mendelian randomization

Introduction

The potential positive impact of computer game playing on cognitive function and its potential role in reducing the risk of Alzheimer's disease (AD) has been suggested. However, current observational studies have certain limitations. We utilized Mendelian randomization (MR) alongside extensive genome-wide association study (GWAS) data to examine the relationship between computer game playing, cognitive function, risk of AD, and levels of brain-derived neurotrophic factor (BDNF).

Methods

We collected datasets on computer game playing, cognition function, risk of AD, and BDNF level from the IEU Open GWAS project. Causal effects were assessed using various MR methods, including inverse variance weighted (IVW), weighted median, MR-Egger, simple mode, and weighted mode. To ensure the accuracy of the results, sensitivity analyses were conducted.

Results

Our analysis revealed a significant association between computer game playing and cognitive function (β  =  0.801, 95% CI: 0.351, 1.328, P  =  0.001). There was no statistically significant association between computer game playing and either BDNF level or risk of AD (β  =  -0.112, 95%CI: -1.315, 1.091, P  =  0.855; OR = 1.000, 95% CI: 1.004, 0.997, P  =  0.891, respectively). We further confirmed the reliability of our evidence through the MR-Egger intercept test, MR-PRESSO global test, Cochran's Q test, and funnel plots.

Conclusion

The results of our study indicate that engaging in computer game playing may confer a safeguarding influence on cognitive function. This underscores the potential advantages associated with computer gaming. Nevertheless, given the constraints inherent in our research, further investigation is warranted to substantiate our findings and delve into the underlying mechanisms.

Exploring the causal relationship between inflammatory cytokines and inflammatory arthritis: A Mendelian randomization study

OBJECTIVES

Previous studies have reported an association between inflammatory cytokines and inflammatory arthritis, including Ankylosing spondylitis (AS), rheumatoid arthritis (RA), and psoriatic arthritis (PsA). This study aims to explore the causal relationship between inflammatory cytokines and AS, RA, and PsA using Mendelian randomization (MR).

METHODS

We conducted a bidirectional two-sample MR analysis using genetic summary data from a publicly available genome-wide association study (GWAS) that included 41 genetic variations of inflammatory cytokines, as well as genetic variant data for AS, RA, and PsA from the FinnGen consortium. The main analysis method used was Inverse variance weighted (IVW) to investigate the causal relationship between exposure and outcome. Additionally, other methods such as MR Egger, weighted median (WM), simple mode, and weighted mode were employed to strengthen the final results. Sensitivity analysis was also performed to ensure the reliability of the findings.

RESULTS

The results showed that macrophage colony-stimulating factor (MCSF) was associated with an increased risk of AS (OR = 1.163, 95 % CI = 1.016-1.33, p = 0.028). Conversely, high levels of TRAIL and beta nerve growth factor (β-NGF) were associated with a decreased risk of AS (OR = 0.892, 95 % CI = 0.81-0.982, p = 0.002; OR = 0.829, 95 % CI = 0.696-0.988, p = 0.036). Four inflammatory cytokines were found to be associated with an increased risk of PsA: vascular endothelial growth factor (VEGF) (OR = 1.161, 95 % CI = 1.057-1.275, p = 0.002); Interleukin 12p70 (IL12p70) (OR = 1.189, 95 % CI = 1.049-1.346, p = 0.007); IL10 (OR = 1.216, 95 % CI = 1.024-1.444, p = 0.026); IL13 (OR = 1.159, 95 % CI = 1.05-1.28, p = 0.004). Interleukin 1 receptor antagonist (IL-1rα) was associated with an increased risk of seropositive RA (OR = 1.181, 95 % CI = 1.044-1.336, p = 0.008). Similarly, genetic susceptibility to inflammatory arthritis was found to be causally associated with multiple inflammatory cytokines. Lastly, the sensitivity analysis supported the robustness of these findings.

CONCLUSIONS

This study provides additional insights into the relationship between inflammatory cytokines and inflammatory arthritis, and may offer new clues for the etiology, diagnosis, and treatment of inflammatory arthritis.

Analysing the causal relationship between potentially protective and risk factors and cutaneous melanoma: A Mendelian randomization study

BACKGROUND

Previous observational studies reported altered melanoma risks in relation to many potential factors, such as coffee intake, smoking habits and photodamage-related conditions. Considering the susceptibility of epidemiological studies to residual confounders, there remains uncertainty about the actual causal roles of these reported factors in melanoma aetiology.

OBJECTIVES

This study aims to investigate the causal association between cutaneous melanoma (CM) and previously reported factors: coffee intake, alcohol consumption, lifetime smoking, socioeconomic status (SES), ease of skin tanning, childhood sunburn and facial ageing, providing insight into its underlying aetiology and preventative strategies.

METHODS

We utilized a two-sample MR analysis on data from the largest meta-analysis summary statistics of confirmed cutaneous melanoma including 30,134 patients. Genetic instrumental variables were constructed by identifying single nucleotide polymorphisms (SNPs) that associate with corresponding factors. Inverse variance weighted (IVW) was the primary MR method. For sensitivity and heterogeneity, MR Egger, weighted median, simple mode, weighted mode and MR Egger intercept tests were examined.

RESULTS

Cutaneous melanoma risks were found to be elevated in association with a predisposition towards ease of skin tanning (IVW: OR = 2.842, 95% CI 2.468-3.274, p < 0.001) and with childhood sunburn history (IVW: OR = 6.317, 95% CI 4.479-8.909, p < 0.001). Repeated MR after removing potential confounders and outliers demonstrated resolved horizontal pleiotropy and statistically significant results that closely mirrored the initial findings. Other potential factors, such as coffee intake, alcohol consumption, smoking and socioeconomic status (SES), indicated insignificant effects on melanoma risk in the analysis, and therefore, our Mendelian randomization study does not support their roles in modifying melanoma risks.

CONCLUSIONS

Our extensive MR analysis provides strong evidence of the causative role of ease of skin tanning and childhood sunburn history in elevating melanoma risk. Curtailing ultraviolet radiation (UVR) exposure may be the single best preventative strategy to reduce melanoma risk.

Genetically Proxied Sarcopenia-Related Muscle Traits and Depression: Evidence from the FinnGen Cohort

BACKGROUND

Sarcopenia and depression are common and often coexist in the elderly. This study aims to determine the impact of sarcopenia-related muscle traits on depression.

METHODS

A two-sample Mendelian randomization (MR) study was performed on the summary-level data from the FinnGen cohort to estimate the causal association of appendicular lean mass (ALM), walking pace, or low hand grip strength with depression. Additionally, multivariable MR (MVMR) was performed to assess the dependence of each muscle trait in the causality and adjust the effect of body mass index (BMI). Supplementary backward MR analyses were performed to estimate the effect of depression on sarcopenia-related muscle traits.

RESULTS

Univariable MR analyses demonstrated that there were causal associations of ALM (odds ratio [OR]: 0.94; 95% confidence interval [CI]: 0.88-0.99), walking pace (OR: 0.53; 95% CI: 0.32-0.88), and low hand grip strength (OR: 1.20; 95% CI: 1.05-1.38) with depression. MVMR analyses showed that ALM was the only trait that had a significant causal relationship with depression (OR: 0.91; 95% CI: 0.85-0.98) after accounting for the other two muscle traits. Moreover, the independent association of ALM with depression remained (OR: 0.92; 95% CI: 0.85-0.99) after being adjusted by BMI. The backward MR analyses showed no causal associations of depression with any sarcopenia-related muscle traits.

CONCLUSION

Low muscle mass independently increases the risk of depression. This study determined the muscle-related risk factors of depression, which may help establish the causality between sarcopenia and depression and provide evidence-based recommendations for improving mental health in the elderly.

Causal associations between dietary habits and CVD: a Mendelian randomisation study

Over the years, numerous observational studies have substantiated that various dietary choices have opposing effects on CVD. However, the causal effect has not yet been established. Thus, we conducted a Mendelian randomisation (MR) analysis to reveal the causal impact of dietary habits on CVD. Genetic variants strongly associated with 20 dietary habits were selected from publicly available genome-wide association studies conducted on the UK Biobank cohort (n 449 210). Summary-level data on CVD were obtained from different consortia (n 159 836-977 323). The inverse-variance weighted method (IVW) was the primary outcome, while MR-Egger, weighted median and MR Pleiotropy RESidual Sum and Outlier were used to assess heterogeneity and pleiotropy. We found compelling evidence of a protective causal effect of genetic predisposition towards cheese consumption on myocardial infarction (IVW OR = 0·67; 95 % CI = 0·544, 0·826; P = 1·784 × 10-4) and heart failure (IVW OR = 0·646; 95 % CI = 0·513, 0·814; P = 2·135 × 10-4). Poultry intake was found to be a detrimental factor for hypertension (IVW OR = 4·306; 95 % CI = 2·158, 8·589; P = 3·416 × 10-5), while dried fruit intake was protective against hypertension (IVW OR = 0·473; 95 % CI = 0·348, 0·642; P = 1·683 × 10-6). Importantly, no evidence of pleiotropy was detected. MR estimates provide robust evidence for a causal relationship between genetic predisposition to 20 dietary habits and CVD risk, suggesting that well-planned diets may help prevent and reduce the risk of CVD.

A Positive Causal Effect of Shrimp Allergy on Major Depressive Disorder Mediated by Allergy- and Immune-Related Pathways in the East Asian Population

BACKGROUND

Observational studies have implied a potential correlation between allergic diseases and major depressive disorder (MDD). However, the relationship is still inconclusive as it is likely to be interfered with by substantial confounding factors and potential reverse causality. The present study aimed to investigate causal correlation of the two diseases by a Mendelian randomization (MR) study and further elucidate the underlying molecular mechanisms.

METHODS

With the biggest summary datasets of a genome-wide association study (GWAS) in the East Asian population, we conducted a two-sample, bidirectional MR study to assess the causal correlation between shrimp allergy (SA) and MDD. Subsequently, we identified the pleiotropic genes' susceptibility to the two diseases at whole-genome and tissue-specific levels, respectively. Enriched GO sets and KEGG pathways were also discovered to elucidate the potential underlying mechanisms.

RESULTS

With the most suitable MR method, SA was identified as a causal risk factor for MDD based on three different groups of independent genetic instruments, respectively (p < 2.81 × 10-2). In contrast, we did not observe a significant causal effect of MDD on SA. The GWAS-pairwise program successfully identified seven pleiotropic genetic variants (PPA3 > 0.8), indicating that the two diseases indeed have a shared genetic basis. At a whole-genome level, the MAGMA program identified 44 pleiotropic genes, which were enriched in allergy-related pathways, such as antigen processing and presentation pathway (p = 1.46 × 10-2). In brain-specific tissue, the S-MultiXcan program found 17 pleiotropic genes that were significantly enriched in immune-related pathways and GO sets, including asthma-related pathway, T-cell activation-related, and major histocompatibility complex protein-related GO sets. Regarding whole-blood tissue, the program identified six pleiotropic genes that are significantly enriched in tolerance induction-related GO sets.

CONCLUSIONS

The present study for the first time indicated a significant causal effect of SA on the occurrence of MDD, but the reverse was not true. Enrichment analyses of pleiotropic genes at whole-genome and tissue-specific levels implied the involvement of allergy and immune-related pathways in the shared genetic mechanism of the two diseases. Elucidating the causal effect and the acting direction may be beneficial in reducing the incidence rate of MDD for the massive group of SA patients in the East Asian region.

Unlocking the Medicinal Mysteries: Preventing Lacunar Stroke with Drug Repurposing

Currently, only the general control of the risk factors is known to prevent lacunar cerebral infarction, but it is unknown which type of medication for controlling the risk factors has a causal relationship with reducing the risk of lacunar infarction. To unlock this medical mystery, drug-target Mendelian randomization analysis was applied to estimate the effect of common antihypertensive agents, hypolipidemic agents, and hypoglycemic agents on lacunar stroke. Lacunar stroke data for the transethnic analysis were derived from meta-analyses comprising 7338 cases and 254,798 controls. We have confirmed that genetic variants mimicking calcium channel blockers were found to most stably prevent lacunar stroke. The genetic variants at or near HMGCR, NPC1L1, and APOC3 were predicted to decrease lacunar stroke incidence in drug-target MR analysis. These variants mimic the effects of statins, ezetimibe, and antisense anti-apoC3 agents, respectively. Genetically proxied GLP1R agonism had a marginal effect on lacunar stroke, while a genetically proxied improvement in overall glycemic control was associated with reduced lacunar stroke risk. Here, we show that certain categories of drugs currently used in clinical practice can more effectively reduce the risk of stroke. Repurposing several drugs with well-established safety and low costs for lacunar stroke prevention should be given high priority when doctors are making decisions in clinical practice. This may contribute to healthier brain aging.

Mendelian randomization study supports the causal effects of air pollution on longevity via multiple age-related diseases

Growing evidence suggests that exposure to fine particulate matter (PM2.5) may reduce life expectancy; however, the causal pathways of PM2.5 exposure affecting life expectancy remain unknown. Here, we assess the causal effects of genetically predicted PM2.5 concentration on common chronic diseases and longevity using a Mendelian randomization (MR) statistical framework based on large-scale genome-wide association studies (GWAS) (>400,000 participants). After adjusting for other types of air pollution and smoking, we find significant causal relationships between PM2.5 concentration and angina pectoris, hypercholesterolaemia and hypothyroidism, but no causal relationship with longevity. Mediation analysis shows that although the association between PM2.5 concentration and longevity is not significant, PM2.5 exposure indirectly affects longevity via diastolic blood pressure (DBP), hypertension, angina pectoris, hypercholesterolaemia and Alzheimer's disease, with a mediated proportion of 31.5, 70.9, 2.5, 100, and 24.7%, respectively. Our findings indicate that public health policies to control air pollution may help improve life expectancy.

Mendelian randomization reveals association between retinal thickness and non-motor symptoms of Parkinson's disease

Retinal thickness is related to Parkinson's disease (PD), but its association with the severity of PD is still unclear. We conducted a Mendelian randomized (MR) study to explore the association between retinal thickness and PD. For the two-sample MR analysis, the summary statistics obtained from genome-wide association studies on the thickness of Retinal nerve fiber layer (RNFL) and ganglion cell inner plexiform layer (GCIPL) were employed as exposure, while the summary statistics associated with PD were used as the outcome. The primary approach utilized was inverse variance weighted. To correct for multiple testing, the false discovery rate (FDR) was employed. For sensitivity analysis, an array of robust MR methods was utilized. We found genetically predicted significant association between reduced RNFL thickness and a reduced risk of constipation in PD (odds ratio [OR] = 0.854, 95% confidence interval [CI] (0.782, 0.933), P < 0.001, FDR-corrected P = 0.018). Genetically predicted reduced RNFL thickness was associated with a reduced Unified Parkinson's Disease Rating Scale total score (β = -0.042, 95% CI (-0.079, 0.005), P = 0.025), and reduced GCIPL thickness was associated with a lower risk of constipation (OR = 0.901, 95% CI (0.821, 0.988), P = 0.027) but a higher risk of depression (OR = 1.103, 95% CI (1.016, 1.198), P = 0.020), insomnia (OR = 1.090, 95% CI (1.013, 1.172), P = 0.021), and rapid eye movement sleep behaviour disorder (RBD) (OR = 1.198, 95% CI (1.061, 1.352), P = 0.003). In conclusion, we identify an association between retinal thickness and non-motor symptoms (constipation, depression, insomnia and RBD) in PD, highlighting the potential of retinal thickness as a biomarker for PD nonmotor symptoms.

Proteomic insights into the associations between obesity, lifestyle factors, and coronary artery disease

BACKGROUND

We aimed to investigate the protein pathways linking obesity and lifestyle factors to coronary artery disease (CAD).

METHODS

Summary-level genome-wide association statistics of CAD were obtained from the CARDIoGRAMplusC4D consortium (60,801 cases and 123,504 controls) and the FinnGen study (R8, 39,036 cases and 303,463 controls). Proteome-wide Mendelian randomization (MR) analysis was conducted to identify CAD-associated blood proteins, supplemented by colocalization analysis to minimize potential bias caused by linkage disequilibrium. Two-sample MR analyses were performed to assess the associations of genetically predicted four obesity measures and 13 lifestyle factors with CAD risk and CAD-associated proteins' levels. A two-step network MR analysis was conducted to explore the mediating effects of proteins in the associations between these modifiable factors and CAD.

RESULTS

Genetically predicted levels of 41 circulating proteins were associated with CAD, and 17 of them were supported by medium to high colocalization evidence. PTK7 (protein tyrosine kinase-7), RGMB (repulsive guidance molecule BMP co-receptor B), TAGLN2 (transgelin-2), TIMP3 (tissue inhibitor of metalloproteinases 3), and VIM (vimentin) were identified as promising therapeutic targets. Several proteins were found to mediate the associations between some modifiable factors and CAD, with PCSK9, C1S, AGER (advanced glycosylation end product-specific receptor), and MST1 (mammalian Ste20-like kinase 1) exhibiting highest frequency among the mediating networks.

CONCLUSIONS

This study suggests pathways explaining the associations of obesity and lifestyle factors with CAD from alterations in blood protein levels. These insights may be used to prioritize therapeutic intervention for further study.

Circulating Metabolites and Dental Traits: A Mendelian Randomization Study

It is of great importance to uncover causal biomarkers to gain insight into the pathogenesis of oral diseases and identify novel treatment targets for prevention and treatment thereof. This study aimed to systematically evaluate the causal effects of hundreds of metabolites on 10 dental traits using a 2-sample Mendelian randomization (MR) approach. Genetic variants from genome-wide association studies of 309 known metabolites were used as instrumental variables. We selected 10 dental traits, including clinical measures of dental diseases, from the Gene-Lifestyle Interactions in Dental Endpoints Consortium and self-reported oral health data from the UK Biobank. The causal relationships between metabolites and dental traits were inferred using the inverse variance-weighted approach and further controlled for horizontal pleiotropy using 5 additional MR methods. After correcting for multiple tests, 5 metabolites were identified as causal biomarkers. Genetically predicted increased levels of mannose were associated with lower risk of bleeding gums (odds ratio [OR] = 0.72; 95% confidence interval [CI], 0.61-0.85; P = 9.9 × 10-5). MR also indicated 4 metabolites on the causal pathway to dentures, with fructose (OR = 0.50; 95% CI, 0.36-0.70; P = 5.2 × 10-5) and 1-palmitoleoyl-glycerophosphocholine (OR = 0.67; 95% CI, 0.56-0.81; P = 4.8 × 10-5) as potential protective factors and glycine (OR = 1.22; 95% CI, 1.11-1.35; P = 5.6×10-5) and 1,5-anhydroglucitol (OR = 1.32; 95% CI, 1.14-1.52; P = 1.5 × 10-4) as risk factors. The causal associations were robust in various sensitivity analyses. We further observed some shared metabolites among different dental traits, implying similar biological mechanisms underlying the pathogenic processes. Finally, the pathway analysis revealed several significant metabolic pathways that may be involved in the development of dental disorders. Our study provides novel insights into the combination of metabolomics and genomics to reveal the pathogenesis of and therapeutic strategies for dental disorders. It highlighted 5 metabolites and several pathways as causal candidates, warranting further investigation.

Mendelian Randomisation Study on Association of Gut Microbiota and Periodontitis

OBJECTIVE

Several studies have demonstrated the possible association between gut microbiota and periodontitis. The mechanism by which gut microbiota contribute to periodontitis remains unknown.

METHODS

A 2-sample Mendelian randomisation (MR) study was conducted using publicly available Genome-Wide Association Studies (GWAS) data of European ancestry. The relationships between gut microbiota and tooth loss and periodontitis were assessed using summary-level data. Moreover, inverse variance weighted (IVW), MR-Egger, weighted median, and simple Mendelian were used. The results were further validated using sensitivity analyses.

RESULTS

A total of 211 gut microbiota were studied, including 9 phyla, 16 classes, 20 orders, 35 families, and 131 genera. The IVW method identified 16 bacterial genera related to the risk of periodontitis and tooth loss. Lactobacillaceae was associated with an increased risk of periodontitis (odds ratio [OR], 1.40, 95% confidence interval [CI], 1.03-1.91, P<.001) and tooth loss (OR, 1.12; 95% CIs, 1.02-1.24, P = .002), whereas Lachnospiraceae UCG008 was linked to a lower risk of tooth loss (P = .041). There was no heterogeneity and horizontal pleiotropy in the sensitivity analysis.

CONCLUSIONS

Several microorganisms were identified to be linked to the risk of periodontitis. Furthermore, the findings improved our understanding of gut microbiota and periodontitis pathology.

Causal effects of rheumatoid arthritis or ankylosing spondylitis on membranous nephropathy: a two-sample Mendelian randomization study

Background

Membranous nephropathy (MN) is the leading cause of adult-onset nephrotic syndrome, with primary MN of unclear cause accounting for 80% of cases. Retrospective clinical research reported that MN occurring in rheumatoid arthritis (RA) and ankylosing spondylitis (AS) patients was triggered by nephrotoxic drugs or of unknown cause. However, whether RA or AS itself increases the risk of developing MN is unknown.

Methods

We conducted mendelian randomization (MR) analysis to evaluate the causal effects of RA or AS on MN using genome-wide association study (GWAS) statistics. The inverse variance weighted (IVW) method was the primary analysis, and several supplementary analyses and sensitivity analyses were performed to test the causal estimates.

Results

We obtained 30 valid instrumental variables (IVs) of RA and 16 valid IVs of AS from large-scale open-access GWASs. The genetically predicted RA significantly increased the risk of MN [IVW odds ratios (OR) = 1.327, 95% confidence interval (CI) = (1.124, 1.565), P = 8.051 × 10-4]. Three supplementary MR analyses provided the consistent positive causal effect of RA on MN (all P < 0.05). No horizontal pleiotropy was detected by MR Egger intercept analysis (P = 0.411). However, the genetically predicted AS had no causal effect on MN by IVW and supplementary analysis (all P > 0.05).

Conclusions

Genetically predicted RA could increase the risk of MN, but genetically predicted AS was not associated with MN. Screening for kidney involvement in RA patients should be noted, and active treatment of RA will reduce the public health burden of MN.

Atopic dermatitis and risk of 14 site-specific cancers: A Mendelian randomization study

BACKGROUND

Atopic dermatitis (AD) accounts for a large proportion of the burden of skin disease, with a prevalence of around 10% among adults worldwide. In addition, systematic reviews and meta-analyses have found that AD is associated with cancer risk at several sites; if found to be causal this could highlight potential treatment targets to reduce cancer risk.

OBJECTIVES

To assess the potential causative link between AD and 14 site-specific cancers in a two-sample randomization study.

METHODS

From the largest genome-wide association study (GWAS) of AD (10,788 cases and 30,047 non-cases), genetic variants highly associated (p < 5 × 10-8 ) with AD in the European population were selected as instrumental variables (IVs). Data from large cancer consortia, as well as the UK Biobank study (n = 442,239) and the FinnGen study (n = 218,792), were employed to assess genetic associations with 14 site-specific cancers and overall cancer. A set of complementary approaches and sensitivity analyses were carried out to examine the robustness of our results. In addition, associations for the same cancer site from different data sources were combined using meta-analyses.

RESULTS

We discovered no strong causal evidence of AD on the risk of overall cancer, with effect estimates close to zero. After the Benjamini-Hochberg correction, the inverse-variance weighted method indicated no association between AD and overall cancer risk in both the UK Biobank (OR, 1.00; 95% CI, 0.94-1.06; FDR, 0.98) and FinnGen studies (OR, 0.96; 95% CI, 0.92-1.02; FDR, 0.68). No strong evidence of an association was found between genetically predicted AD and the risk of any site-specific cancers.

CONCLUSIONS

Our MR investigation does not support a causal effect of AD on cancer risk. This finding has important implications for the prevention and management of both AD and cancer, as it reduces the concern of potential adverse effects of AD on cancer outcomes.

Causal relationship between sleep traits and cognitive impairment: A Mendelian randomization study

OBJECTIVE

Observational studies had demonstrated a link between sleep disturbances and cognitive decline. Here, we aimed to investigate the causal association between genetically predicted sleep traits and cognitive impairment using Mendelian randomization (MR).

METHODS

Using strict criteria, we selected genetic variants from European ancestry Genome-wide association studies (GWAS) from the Sleep Disorders Knowledge Portal and UK Biobank as instrumental variables for several sleep traits, including insomnia, sleep duration, daytime sleepiness, daytime napping, and chronotype. Summary statistics related to cognitive impairment were derived from five different GWAS, including the Social Science Genetic Association Consortium. The role of self-reported sleep trait phenotypes in the etiology of cognitive impairment was explored using inverse-variance weighted (IVW) tests, MR-Egger tests, and weighted medians, and sensitivity analyses were conducted to ensure robustness.

RESULTS

In the main IVW analysis, sleep duration (reaction time: β = -0.05, 95% CI -0.07 to -0.04, p = 1.93×10-12 ), daytime sleepiness (average cortical thickness: β = -0.12, 95% CI -0.22 to -0.02, p = 0.023), and daytime napping (fluid intelligence: β = -0.47, 95% CI -0.87 to -0.07, p = 0.021; hippocampal volume in Alzheimer's disease: β = -0.99, 95% CI -1.64 to -0.35, p = 0.002) were significantly negatively correlated with cognitive performance. However, any effects of insomnia and chronotype on cognitive impairment were not determined.

CONCLUSIONS

Our findings highlighted that focusing on sleep behaviors or distinct sleep patterns-particularly sleep duration, daytime sleepiness, and daytime napping, was a promising approach for preventing cognitive impairment. This study also shed light on risk factors for and potential early markers of cognitive impairment risk factors.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

An augmented Mendelian randomization approach provides causality of brain imaging features on complex traits in a single biobank-scale dataset

Mendelian randomization (MR) is an effective approach for revealing causal risk factors that underpin complex traits and diseases. While MR has been more widely applied under two-sample settings, it is more promising to be used in one single large cohort given the rise of biobank-scale datasets that simultaneously contain genotype data, brain imaging data, and matched complex traits from the same individual. However, most existing multivariable MR methods have been developed for two-sample setting or a small number of exposures. In this study, we introduce a one-sample multivariable MR method based on partial least squares and Lasso regression (MR-PL). MR-PL is capable of considering the correlation among exposures (e.g., brain imaging features) when the number of exposures is extremely upscaled, while also correcting for winner's curse bias. We performed extensive and systematic simulations, and demonstrated the robustness and reliability of our method. Comprehensive simulations confirmed that MR-PL can generate more precise causal estimates with lower false positive rates than alternative approaches. Finally, we applied MR-PL to the datasets from UK Biobank to reveal the causal effects of 36 white matter tracts on 180 complex traits, and showed putative white matter tracts that are implicated in smoking, blood vascular function-related traits, and eating behaviors.

Thyroid dysfunction and risk of cutaneous malignant melanoma: a bidirectional two-sample Mendelian randomization study

Background

Epidemiologic and observational data have found a risk association between thyroid dysfunction and cutaneous malignant melanoma (CMM), however, the cause and direction of these effects are yet unknown. By using a bidirectional two-sample Mendelian randomization (MR) methodology, we hoped to further investigate the causal link between thyroid dysfunction and CMM in this work.

Methods

A genome-wide association study (GWAS) of 9,851,867 single nucleotide polymorphisms (SNPs) in a European population was used to develop genetic tools for thyroid dysfunction. Hypothyroidism was linked to 22,687 cases and 440,246 controls. For hyperthyroidism, there were 3545 cases and 459,388 controls. A total of 3751 cases and 372016 controls were included in the genetic data for CMM from UK Biobank (http://www.nealelab.is/uk-biobank) (the Dataset: ieu - b - 4969). Among them, inverse variance weighting (IVW) is the main MR Analysis method for causality assessment. MR-Egger method, MR Pleiotropic residual and outlier test (MR-PRESSO), and simple and weighted median (VM) were used to supplement the IVW method. Sensitivity analyses, mainly Cochran's Q test, leave-one-out analysis, and MR Egger intercept test were performed to assess the robustness of the outcomes.

Results

The two-sample MR Analysis results revealed a negative correlation between genetically predicted hypothyroidism and the probability of CMM (OR=0.987, 95%CI =0.075-0.999, p=0.041). The supplemental MR Analysis did not reveal any statistically significant differences, although the direction of the effect sizes for the other approaches was consistent with the IVW effect sizes. The results of the causal analysis were relatively robust, according to a sensitivity analysis. The risk of CMM was unaffected by hyperthyroidism (p>0.05). No correlation between CMM and thyroid dysfunction was seen in the reverse MR analysis.

Conclusion

Although the magnitude of the causal association is weak and further investigation of the mechanism of this putative causal relationship is required, our findings imply that hypothyroidism may be a protective factor for CMM.

Causal association between dietary factors and esophageal diseases: A Mendelian randomization study

BACKGROUND

Using Mendelian randomization (MR) approach, our objective was to determine whether there was a causal association between dietary factors and gastroesophageal reflux disease (GERD), Barrett's esophagus (BE), or esophageal cancer (EC).

METHODS

Genome-wide association study (GWAS) data for eighteen types of dietary intake were obtained from the UK Biobank. GWAS data for GERD, BE, and EC were sourced from the FinnGen consortium. We performed univariable and multivariable MR analysis to assess the cause effect between dietary factors and esophageal diseases. MR results were expressed as odds ratios (OR) with 95% confidence intervals (CI).

RESULTS

Raw vegetable intake was associated with a lower risk of GERD (OR = 0.478; P = 0.011). On the contrary, cooked vegetable intake increased the risk of GERD (OR = 1.911; P = 0.024). Bread intake was associated with increased odds of BE (OR = 6.754; P = 0.007), while processed meat intake was associated with reduced risk of BE (OR = 0.210; P = 0.035). We also observed evidence that increased consumption of dried fruit (OR = 0.087; P = 0.022) and salt added to food (OR = 0.346; P = 0.045) could prevent EC. The results of multivariable MR showed that the protective effect of consumption of salt added to food on EC was no longer significant after adjusting for the consumption of dried fruit.

CONCLUSION

Vegetable consumption was associated with GERD, whereas consumption of bread and processed meat was associated with BE. Dried fruit intake was associated with a lower risk of EC, and the protective effect of consumption of salt added food on EC may also be mediated by consumption of dried fruit. Future research should be performed to investigate the mechanisms behind these cause-and-effect relationships to reduce the burden of disease caused by dietary habits.

Appraising the causal relationship between thyroid function and rheumatoid arthritis: a two-sample bidirectional Mendelian randomization study

Background

Hypothyroidism and hyperthyroidism are observationally associated with rheumatoid arthritis (RA), but causality is unclear. To evaluate the causal relationship between thyroid function and RA, we conducted a two-Sample bidirectional Mendelian Randomization (MR) study.

Methods

Single nucleotide polymorphisms associated with six phenotypes were selected from the FinnGen biobank database, The ThyroidOmics Consortium database, and the IEU Open GWAS database. For the forward MR analysis, we selected hypothyroidism (N=213,390), Graves' disease (GD) (N=199,034), other types of hyperthyroidism (N=190,799), free thyroxine (FT4, N=49,269), and thyroid-stimulating hormone (TSH, N=54,288) as the five related thyroid function phenotypes for exposure, with RA (N=58,284) as the outcome. Reverse MR analysis selected RA as the exposure and five phenotypes of thyroid function as the outcome. The Inverse variance weighting (IVW) method was used as the primary analysis method, supplemented by weighted median (WM) and MR-Egger methods. Cochran's Q test, MR-PRESSO, MR-Egger regression methods, and leave-one-out analysis were employed to assess sensitivity and pleiotropy.

Results

Forward MR evidence indicates that genetic susceptibility to hypothyroidism is associated with an increased risk of RA (ORIvw=1.758, P=7.61×10-5). Reverse MR evidence suggests that genetic susceptibility to RA is associated with an increased risk of hypothyroidism (ORIvw=1.274, P=3.88×10-20), GD (ORIvw=1.269, P=8.15×10-05), and other types of hyperthyroidism (ORIvw=1.141, P=1.80×10-03). There is no evidence to support a forward or reverse causal relationship between genetic susceptibility to RA and FT4, TSH.

Conclusion

Our results provide genetic evidence supporting bidirectional causal relationships between thyroid function and RA. These findings inform preventive strategies and interventions targeting RA and thyroid dysfunction.

Socioeconomic status and severe mental disorders: a bidirectional multivariable Mendelian randomisation study

BACKGROUND

Despite the evidence supporting the relationship between socioeconomic status (SES) and severe mental disorders (SMD), the directionality of the associations between income or education and mental disorders is still poorly understood.

OBJECTIVE

To investigate the potential bidirectional causal relationships between genetic liability to the two main components of SES (income and educational attainment (EA)) on three SMD: schizophrenia, bipolar disorder (BD) and depression.

METHODS

We performed a bidirectional, two-sample univariable Mendelian randomisation (UVMR) and multivariable Mendelian randomisation (MVMR) study using SES phenotypes (income, n=397 751 and EA, n=766 345) and SMD (schizophrenia, n=127 906; BD, n=51 710 and depression, n=500 119) genome-wide association studies summary-statistics to dissect the potential direct associations of income and EA with SMD.

FINDINGS

UVMR showed that genetic liability to higher income was associated with decreased risk of schizophrenia and depression, with a smaller reverse effect of schizophrenia and depression on income. Effects were comparable after adjusting for EA in the MVMR. UMVR showed bidirectional negative associations between genetic liability to EA and depression and positive associations between genetic liability to EA and BD, with no significant effects on schizophrenia. After accounting for income, MVMR showed a bidirectional positive direction between genetic liability to EA and BD and schizophrenia but not with depression.

CONCLUSIONS

Our results suggest a heterogeneous link pattern between SES and SMD. We found a negative bidirectional association between genetic liability to income and the risk of schizophrenia and depression. On the contrary, we found a positive bidirectional relationship of genetic liability to EA with schizophrenia and BD, which only becomes apparent after adjusting for income in the case of schizophrenia.

CLINICAL IMPLICATIONS

These findings shed light on the directional mechanisms between social determinants and mental disorders and suggest that income and EA should be studied separately in relation to mental illness.

The association between trimethylamine N-oxide levels and ischemic stroke occurrence: a meta-analysis and Mendelian randomization study

BACKGROUND

Trimethylamine-N-oxide (TMAO), an intestinal microbiota-derived choline metabolite, has been found to be associated with ischemic stroke (IS) in more and more studies. However, the causal role of TMAO on IS occurrence remains perplexing.

METHODS

We comprehensively screened the related clinical studies on PubMed, Web of Science, and Embase. Case-control and cohort studies that reported the TMAO levels of both IS patients and healthy controls were included, and the risk of bias was assessed according to the criteria by the Centre for Evidence-Based Medicine in Oxford, UK. A meta-analysis of the retrieved publications was performed with a random-effect model to analyze the connection between TMAO levels and IS events. Besides, a Mendelian randomization (MR) analysis was performed to study the causal effect of TMAO on IS, with pooled data of TMAO and IS obtained from genome-wide association studies (GWAS). The following methods were used: MR-Egger, weighted median, inverse-variance weighted, simple mode, and weighted mode. The study has been registered in INPLASY (Registration number: INPLASY2023100027).

RESULTS

Eight cohort or case-control studies covering 2444 cases and 1707 controls were identified. The pooled data indicated that the IS patients tended to have higher TMAO levels compared with the controls (mean difference: 1.97 μM; 95% confidence interval [CI]: 0.87, 3.07; P = 0.0005), while distinctive heterogeneity (I2 = 96%, P < 0.00001) was observed. Sub-group analysis revealed that the heterogeneity of the studies might be derived from the studies themselves. However, no causal effect of TMAO on IS was observed (P > 0.05) in the Mendelian randomization analysis of this study.

CONCLUSION

We confirmed that IS patients tend to have higher TMAO levels than healthy individuals, while our findings of MR analysis did not support the causal role of TMAO in IS occurrence. Therefore, more studies are required for a better understanding of the relationship between TMAO levels and IS onset.

The association between walking pace and hand grip strength with the risk of chronic obstructive pulmonary disease: a bidirectional Mendelian randomization study

BACKGROUND

Chronic Obstructive Pulmonary Disease (COPD) currently ranks as the third leading cause of mortality worldwide, imposing substantial burdens on societal and individual health. Amongst health research tools, walking pace (WP) and hand grip strength (HGS) are cornerstones, extensively associated with diverse health conditions. However, the intricate interplay between these factors and COPD risk remains ambiguous. This study aims to elucidate the causal association of WP, HGS, with COPD risk through a bidirectional Mendelian randomization (MR) approach.

METHODS

Bidirectional MR analysis was performed using Genome-wide association study (GWAS) data of European individuals for WP, HGS, and COPD. Inverse Variance Weighted (IVW) served as the primary MR analysis approach. To supplement the IVW findings, four additional MR methods [MR-Egger, weighted median, maximum likelihood, simple median] were used. To assess heterogeneity and pleiotropy, sensitivity analyses were performed. In addition, multivariate MR (MVMR) analysis was used to assess causality after adjustment for potential confounders.

RESULTS

IVW method results show a significant negative association between WP and COPD risk in both initial (genome-wide threshold, odds ratio (OR) = 0.21, 95% confidence interval (CI) 0.09-0.51, P = 5.06 × 10- 4) and secondary (locus-wide threshold, OR = 0.27, 95%CI: 0.18-0.41, P = 4.88 × 10- 10) MR analysis. The reverse MR analysis suggested that COPD also diminishes WP. Additionally, a causal risk reduction for COPD with right HGS (OR = 0.74, 95% CI: 0.58-0.94, P = 1.44 × 10- 2) was only found in secondary MR analysis. The outcomes of the four additional MR methods also suggested similar causal relationships, and sensitivity analyses endorsed their robustness. Lastly, the MVMR analysis demonstrated that the WP's effect on reducing COPD risk persisted independently of potential confounding variables.

CONCLUSION

A bidirectional causal relationship exists between typical WP and COPD risk. Conversely, a decrease in right HGS is unidirectionally associated with an increased risk of COPD. The study suggests that WP may serve as a predictive factor for COPD or as a simple evaluative indicator for prognosis.