Multivariable Mendelian Randomization and Mediation

Causal association between the peripheral immunity and the risk and disease severity of multiple sclerosis

Background

Growing evidence links immunological responses to Multiple sclerosis (MS), but specific immune factors are still unclear.

Methods

Mendelian randomization (MR) was performed to investigate the association between peripheral hematological traits, MS risk, and its severity. Then, further subgroup analysis of immune counts and circulating cytokines and growth factors were performed.

Results

MR revealed higher white blood cell count (OR [95%CI] = 1.26 [1.10,1.44], P = 1.12E-03, P adjust = 3.35E-03) and lymphocyte count (OR [95%CI] = 1.31 [1.15,1.50], P = 5.37E-05, P adjust = 3.22E-04) increased the risk of MS. In further analysis, higher T cell absolute count (OR [95%CI] = 2.04 [1.36,3.08], P = 6.37E-04, P adjust = 2.19E-02) and CD4+ T cell absolute count (OR [95%CI] = 2.11 [1.37,3.24], P = 6.37E-04, P adjust = 2.19E-02), could increase MS risk. While increasing CD25++CD4+ T cell absolute count (OR [95%CI] = 0.75 [0.66,0.86], P = 2.12E-05, P adjust = 1.72E-03), CD25++CD4+ T cell in T cell (OR [95%CI] = 0.79[0.70,0.89], P = 8.54E-05, P adjust = 5.29E-03), CD25++CD4+ T cell in CD4+ T cell (OR [95%CI] = 0.80[0.72,0.89], P = 1.85E-05, P adjust = 1.72E-03), and CD25++CD8+ T cell in T cell (OR [95%CI] = 0.68[0.57,0.81], P = 2.22E-05, P adjust = 1.72E-03), were proved to be causally defensive for MS. For the disease severity, the suggestive association between some traits related to CD4+ T cell, Tregs and MS severity were demonstrated. Moreover, elevated levels of IL-2Ra had a detrimental effect on the risk of MS (OR [95%CI] = 1.22 [1.12,1.32], P = 3.20E-06, P adjust = 1.34E-04).

Conclusions

This study demonstrated a genetically predicted causal relationship between elevated peripheral immune cell counts and MS. Subgroup analysis revealed a specific contribution of peripheral immune cells, holding potential for further investigations into the underlying mechanisms of MS and its severity.

Smoking and BMI mediate the causal effect of education on lower back pain: observational and Mendelian randomization analyses

Objective

Low back pain (LBP) has been associated with education in previous observational studies, but the causality remains unclear. This study aims to assess the impact of education on LBP and to explore mediation by multiple lifestyle factors.

Design

Univariable Mendelian randomization (MR) was performed to examine the overall effect of education on LBP. Subsequently, multivariable MR was conducted to assess both the direct effect of education on LBP and the influence of potential mediators. Indirect effects were estimated using either the coefficient product method or the difference method, and the proportion of mediation was calculated by dividing the indirect effect by the total effect. The observational study utilized data from the NHANES database collected between 1999 and 2004, and included 15,580 participants aged 20 years and above.

Results

Increasing education by 4.2 years leads to a 48% reduction in the risk of LBP (OR=0.52; 95% CI: 0.46 to 0.59). Compared to individuals with less than a high school education, those with education beyond high school have a 28% lower risk of LBP (OR=0.72; 95% CI: 0.63 to 0.83). In the MR study, smoking accounts for 12.8% (95% CI: 1.04% to 20.8%) of the total effect, while BMI accounts for 5.9% (95% CI: 2.99% to 8.55%). The combined mediation effect of smoking and BMI is 27.6% (95% CI: 23.99% to 32.7%). In the NHANES study, only smoking exhibits a mediating effect, accounting for 34.3% (95% CI: 21.07% to 41.65%) of the effect, while BMI does not demonstrate a mediating role.

Conclusions

Higher levels of education provide a protective effect against the risk of LBP. Additionally, implementing interventions to reduce smoking and promote weight loss among individuals with lower levels of education can also decrease this risk.

Associations between genetically predicted iron status and cardiovascular disease risk: A Mendelian randomization study

Background

Mendelian randomization (MR) studies suggest a causal effect of iron (Fe) status on cardiovascular disease (CVD) risk, but it is unknown if these associations are confounded by pleiotropic effects of the instrumental variables (IV) on CVD risk factors. We aimed to investigate the effect of Fe status on CVD risk controlling for CVD risk factors.

Methods

Fe biomarker IVs (total Fe binding capacity (TIBC, n=208,422), transferrin saturation (TSAT, n=198,516), serum Fe (SI, n=236,612), ferritin (n=257,953)) were selected from a European GWAS meta-analysis. We performed two-sample univariate (UV) MR of each Fe trait on CVD outcomes (all-cause ischemic stroke (IS), cardioembolic IS (CES), large artery IS (LAS), small vessel IS (SVS), and coronary heart disease (CHD)) from MEGASTROKE (n=440,328) and CARDIoGRAMplusC4D (n=183,305). We then implemented multivariate (MV) MR conditioning on six CVD risk factors from independent European samples to evaluate their potential confounding and/or mediating effects on the observed Fe-CVD associations.

Results

With UVMR analyses, we found higher genetically predicted Fe status to be associated with a greater risk of CES (TSAT: OR 1.17 [95%CI 1.03, 1.33], SI: OR 1.21 [ 95%CI 1.02, 1.44]; TIBC: OR 0.81 [95%CI 0.69, 0.94]). The detrimental effects of Fe status on CES risk remained unaffected when adjusting for CVD risk factors (all P<0.05). Additionally, we found diastolic blood pressure (DBP) to mediate between 7.1-8.8% of the total effect of Fe status on CES incidence. While UVMR initially suggested a protective effect of Fe status on LAS and CHD, MVMR analyses factoring CVD risk factors revealed a complete annulment of this perceived protective effect (all P>0.05).

Discussion

Higher Fe status was associated with a greater risk of CES independent of CVD risk factors, and this effect was partly mediated by DBP. These findings support a role of Fe status as a modifiable risk factor for CES.

Causal Relationship between Immune Cells and Gynecological Cancers through Bidirectional and Multivariable Mendelian Randomization Analyses

Background: Evidence suggests potential associations between gynecological malignancies and various immune cell chemicals and systems. However, establishing a causal relationship remains uncertain. Methods: This work employed Wald ratio for one single-nucleotide polymorphism (SNP) or the inverse-variance weighted method (IVW) for multiple SNPs to conduct bidirectional two-sample Mendelian randomization (MR) analysis by utilizing genome-wide association study (GWAS) data. We employed supplementary methods, including MR-Egger and weighted median methods, to detect and correct for the influence of horizontal pleiotropy. In addition, we also use colocalization analysis for further validation. Results: In IVW analysis, increases in relative count of circulating CD11c+ HLA-DR++ conventional dendritic cells (cDC) were associated with an elevated risk of breast cancer (OR [95% CI], 1.1295 [1.0632-1.2000], P = 8.044 × 10-5), while elevated levels of HLA-DR on plasmacytoid dendritic cells (DC) and HLA-DR on DC were protective against breast cancer. In addition, actual count of CD39+ resting Treg AC was also shown to be causally associated with the development of ovarian cancer, whereas a high relative count of CD28+ CD45RA- CD8+ T cells reduced the risk of cervical cancer. Sensitivity analysis revealed almost no evidence of bias in the current study. Multivariable MR (MVMR) analyses further confirmed a direct impact of the CD11c+ HLA-DR++ cDC immune phenotype on breast cancer. Colocalization analysis showed the lead SNP, rs780094, suggesting HLA-DR GWAS shared a common genetic mechanism with breast cancer. Conclusions: The MR study identified significant causal relationships between multiple immunophenotypes and breast cancer, aiming to provide clinicians with some reference for cancer prediction and explore further potential associations between immune phenotypes and gynecologic tumors.

Mendelian randomization analysis reveals an independent causal relationship between four gut microbes and acne vulgaris

Background

Numerous studies have suggested a correlation between gut microbiota and acne vulgaris; however, no specific causal link has been explored.

Materials and methods

To investigate the possible causal relationship between gut microbiota and acne vulgaris, this study employed a large-scale genome-wide association study (GWAS) summary statistic. Initially, a two-sample Mendelian randomization (MR) analysis was utilized to identify the specific gut microflora responsible for acne vulgaris. We used the Inverse Variance Weighted (IVW) method as the main MR analysis method. Additionally, we assessed heterogeneity and horizontal pleiotropy, while also examining the potential influence of individual single-nucleotide polymorphisms (SNPs) on the analysis results. In order to eliminate gut microbiota with reverse causal associations, we conducted reverse MR analysis. Multivariate Mendelian randomization analysis (MVMR) was then employed to verify the independence of the causal associations. Finally, we performed SNP annotation on the instrumental variables of independent gut microbiota and acne vulgaris to determine the genes where these genetic variations are located. We also explored the biological functions of these genes through enrichment analysis.

Result

The IVW method of forward MR identified nine gut microbes with a causal relationship with acne vulgaris (p < 0.05). The findings from the sensitivity analysis demonstrate the absence of heterogeneity or horizontal pleiotropy, and leave-one-out analysis indicates that the results are not driven by a single SNP. Additionally, the Reverse MR analysis excluded two reverse-correlated pathogenic gut microbes. And then, MVMR was used to analyze seven gut microbes, and it was found that Cyanobacterium and Family XIII were risk factors for acne vulgaris, while Ruminococcus1 and Ruminiclostridium5 were protective factors for acne vulgaris. After conducting biological annotation, we identified six genes (PLA2G4A, FADS2, TIMP17, ADAMTS9, ZC3H3, and CPSF4L) that may be associated with the pathogenic gut microbiota of acne vulgaris patients. The enrichment analysis results indicate that PLA2G4A/FADS2 is associated with fatty acid metabolism pathways.

Conclusion

Our study found independent causal relationships between four gut microbes and acne vulgaris, and revealed a genetic association between acne vulgaris patients and gut microbiota. Consider preventing and treating acne vulgaris by interfering with the relative content of these four gut microbes.

Exploring Causal Relationships between Leukocyte Telomere Length, Sex Hormone-Binding Globulin Levels, and Osteoporosis Using Univariable and Multivariable Mendelian Randomization

OBJECTIVE

Recent evidence supports that leukocyte telomere length (LTL) may be positively associated with healthy living and inversely correlated with the risk of age-related diseases, including osteoporosis. Furthermore, it is important to note that sex hormone-binding globulin (SHBG) levels play a crucial role in the regulation of osteoporosis by influencing the availability of sex hormones. Hence, this study holds significant importance as it aims to unravel the roles of LTL and SHBG levels and determine which one acts as a predominant intermediary factor in influencing osteoporosis. Using Mendelian randomization (MR), we can gain valuable insights into the intricate relationships between aging, sex hormones, and bone health.

METHODS

Univariable and multivariable and MR analyses were employed in this study. First, we used genetic variants associated with both LTL, as determined from a study involving 472,174 European participants by Codd et al., and SHBG levels, as identified in a study conducted by Ruth et al. with 370,125 participants, as instrumental variables (IVs). Then we aimed to establish a causal relationship between LTL and SHBG levels and their potential impact on osteoporosis using univariable MR. Finally, we conducted multivariable MR to provide insights into the independent and combined effects of LTL, SHBG levels on osteoporosis risk. We used various MR methods, with the primary analysis employing the inverse-variance weighted (IVW) model.

RESULTS

Univariable MR analysis reveals a potential causal effect of longer LTL on reduced risk of osteoporosis [odds ratio (OR): 0.85; 95% confidence interval (CI): 0.73-0.99; p = 0.03]. Conversely, higher genetically determined SHBG levels affect the risk of osteoporosis positively. (OR: 1.38; 95% CI: 1.09-1.75; p < 0.01). We observed a negative causal effect for LTL on the occurrence of SHBG (OR: 0.96; 95% CI 0.94-0.98, p < 0.01). After adjustment of using multivariable MR, the causal effect of LTL on osteoporosis (OR: 0.92; 95% CI: 0.84-1.03; p = 0.14), and the effect of SHBG on osteoporosis (OR: 1.43; 95% CI: 1.16-1.75; p < 0.01) were observed.

CONCLUSION

Longer LTL may confer a protective effect against osteoporosis. Additionally, the levels of SHBG appear to play a crucial role in mediating the relationship between LTL and osteoporosis. By understanding the interplay between these factors, we can gain valuable insights into the mechanisms underlying bone health and aging and potentially identify new avenues for prevention and intervention strategies.

Causal relationship between inflammatory factors and cerebral small vessel disease: Univariate, multivariate, and summary-data-based mendelian randomization analysis

OBJECTIVE

To explore the impact of inflammatory factors on the incidence of cerebral small vessel disease (CSVD), we performed a mendelian randomization (MR) study to analyze the causal relationship between multiple inflammatory factors and CSVD imaging markers and utilized summary-data-based mendelian randomization (SMR) analysis to infer whether the impact of instrumental variables (IVs) on disease is mediated by gene expression or DNA methylation.

METHODS

Using public databases such as UKB and IEU, and original genome-wide association studies, we obtained IVs related to exposure (inflammatory factors) and outcome (CSVD imaging markers). We performed the inverse variance weighted, weighted median, and MR-Egger methods to assess causal effects between exposure and outcome in univariate MR analysis. To evaluate their heterogeneity, a series of sensitivity analyses were conducted, including the Cochrane Q test, MR-Egger intercept test, MR-Presso, and leave-one-out analysis. We also applied mediation and multivariate MR analysis to explore the interactions between positive exposures on the same outcome. Additionally, we conducted the SMR, which utilizes instruments within or near relevant genes in blood or brain tissues, to elucidate the causal associations with CSVD markers.

RESULTS

ABO Univariate MR of multiple cohorts revealed that the risk of small vessel stroke (SVS) increases with elevated levels of TNF-related apoptosis-inducing ligand (TRAIL, OR, 1.23, 95% CI, 1.08-1.39) and interleukin-1 receptor-like 2, (IL-1RL2, OR, 1.29, 95% CI, 1.04-1.61). IL-18 was a potential risk factor for extensive basal ganglia perivascular space burden (BGPVS, OR, 1.02, 95% CI, 1.00-1.05). Moreover, the risk of extensive white matter perivascular space burden (WMPVS) decreased with rising levels of E-selectin (OR, .98, 95% CI, .97-1.00), IL-1RL2 (OR, .97, 95% CI, .95-1.00), IL-3 receptor subunit alpha (IL-3Ra, OR, .98, 95% CI, .97-1.00), and IL-5 receptor subunit alpha (IL-5Ra, OR, .98, 95% CI, .97-1.00). Mediation and multivariate MR analysis indicated that E-selectin and IL-3Ra might interact during the pathogenesis of WMPVS. SMR estimates showed that TRAIL-related IVs rs5030044 and rs2304456 increased the risk of SVS by increasing the expression of gene Kininogen-1 (KNG1) in the cerebral cortex, particularly in the frontal cortex (βsmr = .10, Psmr = .003, FDR = .04). Instruments (rs507666 and rs2519093) related to E-selectin and IL-3Ra could increase the risk of WMPVS by enhancing DNA methylation of the gene ABO in blood tissue (βsmr = .01-.02, Psmr = .001, FDR = .01-.03).

CONCLUSION

According to MR and SMR analysis, higher levels of TRAIL increased the risk of SVS by upregulating gene expression of KNG1 in brain cortex tissues. In addition, protective effects of E-selectin and IL-3a levels on WMPVS were regulated by increased DNA methylation of gene ABO in blood tissue.

Genetic prediction of modifiable lifestyle factors for erectile dysfunction

Background

The causal relationship between certain lifestyle factors and erectile dysfunction (ED) is still uncertain.

Aim

The study sought to investigate the causal effect of 9 life factors on ED through 2-sample single-variable Mendelian randomization (SVMR) and multivariable Mendelian randomization (MVMR).

Methods

Genetic instruments to proxy 9 risk factors were identified by genome-wide association studies. The genome-wide association studies estimated the connection of these genetic variants with ED risk (n = 223 805). We conducted SVMR, inverse variance-weighting, Cochran's Q, weighted median, MR-Egger, MR-PRESSO (Mendelian Randomization Pleiotropy RESidual Sum and Outlier), and MVMR analyses to explore the total and direct relationship between life factors and ED.

Outcomes

The primary outcome was defined as self or physician-reported ED, or using oral ED medication, or a history of surgery related to ED.

Results

In SVMR analyses, suggestive associations with increased the risk of ED were noted for ever smoked (odds ratio [OR], 5.894; 95% confidence interval [CI], 0.469 to 3.079; P = .008), alcohol consumption (OR, 1.495; 95% CI, 0.044 to 0.760; P = .028) and body mass index (BMI) (OR, 1.177; 95% CI, 0.057 to 0.268; P = .003). Earlier age at first intercourse was significantly related to reduced ED risk (OR, 0.659; 95% CI, -0.592 to -0.244; P = 2.5 × 10-6). No strong evidence was found for the effect of coffee intake, time spent driving, physical activity, and leisure sedentary behaviors on the incidence of ED (All P > .05). The result of MVMR analysis for BMI (OR, 1.13; 95% CI, 1.01 to 1.25; P = .045) and earlier age at first intercourse (OR, 0.77; 95% CI, 0.56 to 0.99; P = .018) provided suggestive evidence for the direct impact on ED, while no causal factor was detected for alcoholic drinks per week and ever smoked.

Clinical implications

This study provides evidence for the impact of certain modifiable lifestyle factors on the development of ED.

Strengths and limitations

We performed both SVMR and MVMR to strengthen the causal relationship between exposures and outcomes. However, the population in this study was limited to European ancestry.

Conclusion

Ever smoked, alcoholic drinks per week, BMI, and age first had sexual intercourse were causally related to ED, while the potential connection between coffee intake, physical activity, recreational sedentary habits, and increased risk of ED needs to be further confirmed.

Causal relationship between prostate cancer and 12 types of cancers: multivariable and bidirectional Mendelian randomization analyses

BACKGROUND

Previous observational studies have shown an association between certain cancers and the subsequent risk of prostate cancer (PCa). However, the causal relationship between these cancers and PCa is still unclear. This study aimed to investigate the causal relationship between 12 common cancers and the risk of PCa.

METHODS

We employed genome-wide association studies (GWAS) to perform forward and reverse Mendelian randomization (MR) within two-sample frameworks. Furthermore, we conducted multivariable MR analyses to investigate the relationships between different types of cancer. In addition, multiple sensitivity analysis methods were employed to assess the robustness of our findings.

RESULTS

Our univariable MR analysis showed that genetically predicted hematological cancer was associated with a reduced risk of PCa (OR: 0.911, 95% CI 0.89-0.922, P = 0.03). Furthermore, MR analysis demonstrates that genetically predicted occurrence of thyroid gland and endocrine gland cancer also raised the risk of PCa (all P < 0.05). Multivariable analysis showed that thyroid gland cancer exhibited a higher incidence of PCa (OR: 1.12, 95% CI: 1.08-1.16, P = 0.008). In the reverse MR analysis, we found no significant inverse causal associations between PCa and 12 types of cancers.

CONCLUSION

In summary, this study provided insights into the causal relationships between various types of cancer and PCa. Hematological cancer was suggested to associate with a lower risk of PCa, while thyroid gland cancer and endocrine gland cancer might increase the risk. These findings contribute to the understanding of genetic factors related to PCa and its potential associations with other cancers.

Correlations between genetically predicted lipid-lowering drug targets and inflammatory bowel disease

BACKGROUND

Millions of individuals globally suffer from Inflammatory bowel diseases (IBDs). There is a dearth of large population-based investigations on lipid metabolism and IBDs, and it is unclear whether lipid-lowering drugs target IBDs causally. Consequently, the aim of this study was to investigate the effects of lipid-lowering medication targets on the occurrence and progression of IBDs.

METHODS

Among the more than 400,000 participants in the UK Biobank cohort and the more than 170,000 participants in the Global Lipids Genetics Consortium, a total of nine genes linked to lipid-lowering drug targets were obtained (ABCG5/ABCG8, APOB, APOC3, LDLR, LPL, HMGCR, NPC1L1, PCSK9, and PPARA). IBD data were acquired from de Lange et al. (patients/sample size of IBDs: 25042/59957; ulcerative colitis (UC): 12366/45,975; Crohn's disease (CD): 12194/40,266) and the FinnGen cohort (patients/total sample size of IBDs: 4420/176,899; CD: 1520/171,906; UC: 3325/173,711). All four datasets were cross-combined for validation via Mendelian randomization analysis, and potential mediating factors were explored via mediation analysis.

RESULTS

Genetically proxied APOC3 inhibition was related to increased IBD risk (odds ratio (95% confidence interval): 0.87 (0.80-0.95); P < 0.01) and UC risk (0.83 (0.73-0.94); P < 0.01). IBD and CD risk were reduced by genetic mimicry of LDLR and LPL enhancements, respectively (odds ratioLDLR: 1.18 (1.03-1.36); P = 0.018; odds ratioCD: 1.26 (1.11-1.43); P = 2.60E-04). Genetically proxied HMGCR inhibition was associated with increased CD risk (0.68 (0.50-0.94); P = 0.018). These findings were confirmed through Mendelian analysis of the cross-combination of four separate datasets. APOC3-mediated triglyceride levels may contribute to IBDs partly through mediated triglycerides, Clostridium sensu stricto 1, Clostridiaceae 1, or the Lachnospiraceae FCS020 group. LDLR enhancement may contribute to IBDs partly through increasing Lactobacillaceae.

CONCLUSION

Vigilance is required to prevent adverse effects on IBDs (UC) for patients receiving volanesorsen (an antisense oligonucleotide targeting ApoC3 mRNA) and adverse effects on CD for statin users. LPL and LDLR show promise as candidate drug targets for CD and IBD, respectively, with mechanisms that are potentially independent of their lipid-lowering effects.

Association of thrombopoietin-related drugs with thromboembolic events: Mendelian randomization and a real-world study

Background

Studies have shown conflicting results when using thrombopoietin-related drugs (TPORD) for thromboembolic events (TEEs). Our study aimed to explore the correlation between TPORDs and TEEs.

Method

Drug-targeted Mendelian randomization (MR) and multivariate MR (MVMR) analysis were used to explore the causal relationship between TPORDs and TEEs such as venous thromboembolism (VTE), deep vein thrombosis (DVT), pulmonary embolism (PE), myocardial infarction (MI) and ischemic stroke (STR). At the same time, a real-world study was conducted by extracting adverse events (AEs) from the FDA Adverse Event Reporting System database included in AERSMine to further validate our findings.

Outcome

In drug-target MR, TPORDs were associated with VTE (OR = 1.193, 95% confidence interval (CI): 1.001-1.423, p = 0.049], DVT (OR = 1.321, 95% CI: 1.027-1.700, p = 0.030), MI (OR = 1.216, 95% CI: 1.010-1.464, p = 0.039), STR (OR = 1.224, 95% CI: 1.021-1.468, p = 0.029). VTE/DVT/STR remained stable in MVMR (VTE: OR = 1.3, 95% CI: 1.187-1.422, p < 0.001; DVT: OR = 1.465,95% CI:1.285-1.671, p < 0.001; STR: OR = 1.119, 95% CI: 1.018-1.229, p = 0.019) and real-world studies [lower bound of proportional reporting ratio (ROR) greater than 1]. The significance of myocardial infarction disappeared in MVMR (OR = 0.996, 95% CI: 0.894-1.109, p = 0.942) and in real-world studies (lower ROR lower than 1). There was no evidence of a causal relationship between TPORD and PE (OR = 1.244, 95% CI: 0.969-1.597, p = 0.087), but it generated a signal from a real-world study (lower bound of ROR greater than 1).

Conclusion

This study suggests that TPORDs may be associated with an increased risk of TEEs, particularly AEs leading to VTE/DVT/STR. In addition, the relationship between TPORDs and PE/MI is debatable and requires more research.

Gut microbiota, circulating metabolites, and gallstone disease: a Mendelian randomization study

Background

The link between Gut microbiota (GM) and Gallstone disease (GSD) is well established, but it is not clear whether there is a causal relationship between the two associations.

Methods

We conducted bidirectional Mendelian randomization (MR) analyses, leveraging aggregated data from the Genome-Wide Association Study (GWAS) of GM and Circulating Metabolites. Our primary objective was to investigate the causal interplay between intestinal flora and GSD. Additionally, we performed mediational analyses, two-step MR, and multivariate MR to uncover the potential mediating effect of circulating metabolites in this relationship.

Result

Our study has revealed a causal relationship between GSD and six distinct bacterial groups. Genetically predicted Class Bacilli (Odds Ratio (OR): 0.901, 95% Confidence Interval (95% CI): 0.825-0.985; p = 0.021), Order Lactobacillales (OR: 0.895, 95% CI: 0.816-0.981; p = 0.017), and Genus Coprococcus 2 (OR: 0.884, 95% CI: 0.804-0.973; p = 0.011) were inversely associated with the risk of GSD. Conversely, the Genus Clostridiumsensustricto1 (OR: 1.158, 95% CI: 1.029-1.303; p = 0.015), Genus Coprococcus3 (OR: 1.166, 95% CI: 1.024-1.327; p = 0.020), and Genus Peptococcus (OR: 1.070, 95% CI: 1.017-1.125; p = 0.009) were positively associated with the risk of GSD. Moreover, our findings suggest that the positive influence of the Genus Peptococcus on GSD may be mediated through Omega-3 polyunsaturated fatty acids (PUFA).

Conclusion

This study reinforces the connection between the gut microbiome and the risk of GSD while also unveiling the mediating role of Omega-3 PUFA in the causal relationship between these factors.

The Influence of Autoimmune Thyroid Diseases on Viral Pneumonia Development, Including COVID-19: A Two-Sample Mendelian Randomization Study

The association between thyroid function and viral pneumonia has undergone extensive examination, yet the presence of a causal link remains uncertain. The objective of this paper was to employ Two-Sample Mendelian Randomization (MR) analysis to investigate the connections between three thyroid diseases and thyroid hormone indicators with viral pneumonia and COVID-19. We obtained summary statistics datasets from seven genome-wide association studies (GWASs). The primary method used for estimating relationships was inverse-variance weighting (IVW). In addition, we employed weighted median, weighted mode, MR-Egger, and MR-PRESSO as supplementary analytical tools. Sensitivity analyses encompassed Cochran's Q test, MR-Egger intercept test, and MR-PRESSO. Our study revealed significant causal relationships between having a genetic predisposition to autoimmune thyroid disease (AITD) and an increased susceptibility to viral pneumonia (odds ratio [OR]: 1.096; 95% confidence interval [CI]: 1.022-1.176). Moreover, it demonstrated a heightened susceptibility and severity of COVID-19 (OR for COVID-19 susceptibility, COVID-19 hospitalization, and COVID-19 critical illness, with 95% CIs of 1.016, 1.001-1.032; 1.058, 1.003-1.116; 1.045, 1.010-1.081). However, no statistically significant associations were found between TSH, FT4, subclinical hypo- or hyperthyroidism, and the risk of viral pneumonia incidence, or the susceptibility and severity of COVID-19 (all p > 0.05). This study establishes a cause-and-effect relationship between AITD and the development of viral pneumonia, as well as the susceptibility and severity of COVID-19.

Causal association of leisure sedentary behavior and cervical spondylosis, sciatica, intervertebral disk disorders, and low back pain: a Mendelian randomization study

Background

Some studies suggest sedentary behavior is a risk factor for musculoskeletal disorders. This study aimed to investigate the potential causal association between leisure sedentary behavior (LSB) (including television (TV) viewing, computer use, and driving) and the incidence of sciatica, intervertebral disk degeneration (IVDD), low back pain (LBP), and cervical spondylosis (CS).

Methods

We obtained the data of LSB, CS, IVDD, LBP, sciatica and proposed mediators from the gene-wide association studies (GWAS). The causal effects were examined by Inverse Variance Weighted (IVW) test, MR-Egger, weighted median, weighted mode and simple mode. And sensitivity analysis was performed using MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) and MR-Egger intercept test. Multivariable MR (MVMR) was conducted to investigate the independent factor of other LSB; while two-step MR analysis was used to explore the potential mediators including Body mass index (BMI), smoking initiation, type 2 diabetes mellitus (T2DM), major depressive disorder (MDD), schizophrenia, bipolar disorder between the causal association of LSB and these diseases based on previous studies.

Results

Genetically associated TV viewing was positively associated with the risk of CS (OR = 1.61, 95%CI = 1.25 to 2.07, p = 0.002), IVDD (OR = 2.10, 95%CI = 1.77 to 2.48, p = 3.79 × 10-18), LBP (OR = 1.84, 95%CI = 1.53 to 2.21, p = 1.04 × 10-10) and sciatica (OR = 1.82, 95% CI = 1.45 to 2.27, p = 1.42 × 10-7). While computer use was associated with a reduced risk of IVDD (OR = 0.66, 95%CI = 0.55 to 0.79, p = 8.06 × 10-6), LBP (OR = 0.49, 95%CI = 0.40 to 0.59, p = 2.68 × 10-13) and sciatica (OR = 0.58, 95%CI = 0.46 to 0.75, p = 1.98 × 10-5). Sensitivity analysis validated the robustness of MR outcomes. MVMR analysis showed that the causal effect of TV viewing on IVDD (OR = 1.59, 95%CI = 1.13 to 2.25, p = 0.008), LBP (OR = 2.15, 95%CI = 1.50 to 3.08, p = 3.38 × 10-5), and sciatica (OR = 1.61, 95%CI = 1.03 to 2.52, p = 0.037) was independent of other LSB. Furthermore, two-step MR analysis indicated that BMI, smoking initiation, T2DM may mediate the causal effect of TV viewing on these diseases.

Conclusion

This study provides empirical evidence supporting a positive causal association between TV viewing and sciatica, IVDD and LBP, which were potentially mediated by BMI, smoking initiation and T2DM.

MRSL: a causal network pruning algorithm based on GWAS summary data

Causal discovery is a powerful tool to disclose underlying structures by analyzing purely observational data. Genetic variants can provide useful complementary information for structure learning. Recently, Mendelian randomization (MR) studies have provided abundant marginal causal relationships of traits. Here, we propose a causal network pruning algorithm MRSL (MR-based structure learning algorithm) based on these marginal causal relationships. MRSL combines the graph theory with multivariable MR to learn the conditional causal structure using only genome-wide association analyses (GWAS) summary statistics. Specifically, MRSL utilizes topological sorting to improve the precision of structure learning. It proposes MR-separation instead of d-separation and three candidates of sufficient separating set for MR-separation. The results of simulations revealed that MRSL had up to 2-fold higher F1 score and 100 times faster computing time than other eight competitive methods. Furthermore, we applied MRSL to 26 biomarkers and 44 International Classification of Diseases 10 (ICD10)-defined diseases using GWAS summary data from UK Biobank. The results cover most of the expected causal links that have biological interpretations and several new links supported by clinical case reports or previous observational literatures.

Causal relationship between COVID-19 and chronic pain: A mendelian randomization study

OBJECTIVE

COVID-19 is a highly transmissible disease that can result in long-term symptoms, including chronic pain. However, the mechanisms behind the persistence of long-COVID pain are not yet fully elucidated, highlighting the need for further research to establish causality. Mendelian randomization (MR), a statistical technique for determining a causal relationship between exposure and outcome, has been employed in this study to investigate the association between COVID-19 and chronic pain.

MATERIAL AND METHODS

The IVW, MR Egger, and weighted median methods were employed. Heterogeneity was evaluated using Cochran's Q statistic. MR Egger intercept and MR-PRESSO tests were performed to detect pleiotropy. The Bonferroni method was employed for the correction of multiple testing. R software was used for all statistical analyses.

RESULT

Based on the IVW method, hospitalized COVID-19 patients exhibit a higher risk of experiencing lower leg joint pain compared to the normal population. Meanwhile, the associations between COVID-19 hospitalization and back pain, headache, and pain all over the body were suggestive. Additionally, COVID-19 patients requiring hospitalization were found to have a suggestive higher risk of experiencing neck or shoulder pain and pain all over the body compared to those who did not require hospitalization. Patients with severe respiratory-confirmed COVID-19 showed a suggestive increased risk of experiencing pain all over the body compared to the normal population.

CONCLUSION

Our study highlights the link between COVID-19 severity and pain in different body regions, with implications for targeted interventions to reduce COVID-19 induced chronic pain burden.

Causal role of lipid metabolism in pulmonary alveolar proteinosis: an observational and mendelian randomisation study

BACKGROUND

Pulmonary alveolar proteinosis (PAP) is a rare interstitial lung disease characterised by the accumulation of lipoprotein material in the alveoli. Although dyslipidaemia is a prominet feature, the causal effect of lipid traits on PAP remains unclear. This study aimed to explore the role of lipid traits in PAP and evaluate the potential of lipid-lowering drug targets in PAP.

METHODS

Clinical outcomes, lipid profiles and lung function tests were analysed in a clinical cohort of diagnosed PAP patients and propensity score-matched healthy controls. Genome-wide association study data on PAP, lipid metabolism, blood cells and variants of genes encoding potential lipid-lowering drug targets were obtained for Mendelian randomisation (MR) and mediation analyses.

FINDINGS

Observational results showed that higher levels of total cholesterol (TC), triglycerides and low-density lipoprotein (LDL) were associated with increased risks of PAP. Higher levels of TC and LDL were also associated with worse PAP severity. In MR analysis, elevated LDL was associated with an increased risk of PAP (OR: 4.32, 95% CI: 1.63 to 11.61, p=0.018). Elevated monocytes were associated with a lower risk of PAP (OR 0.34, 95% CI: 0.18 to 0.66, p=0.002) and mediated the risk impact of LDL on PAP. Genetic mimicry of PCSK9 inhibition was associated with a reduced risk of PAP (OR 0.03, p=0.007).

INTERPRETATION

Our results support the crucial role of lipid and metabolism-related traits in PAP risk, emphasising the monocyte-mediated, causal effect of elevated LDL in PAP genetics. PCSK9 mediates the development of PAP by raising LDL. These finding provide evidence for lipid-related mechanisms and promising lipid-lowering drug target for PAP.

Causal associations between liver enzymes and diabetic microvascular complications: A univariable and multivariable Mendelian randomization

BACKGROUND

Observational studies show that liver enzymes are diabetes risk factors. However, previous observational investigations on the relationship between liver enzymes and diabetic microvascular complications produced contradictory results. The purpose of this research is to examine the independent causal effects of liver enzymes on diabetic microvascular complications.

METHODS

Univariable Mendelian randomization (UVMR) and multivariable Mendelian randomization (MVMR) were utilized to disentangle the causal effects. The genome-wide association study (GWAS) summary-level statistics were collected from the UK biobank and the FinnGen consortium. Single nucleotide polymorphisms (SNPs) were selected as genetic instruments with genome-wide significance (p < 5 ×10-8). Five UVMR approaches, including inverse variance weighted (IVW), Bayesian weighted Mendelian randomization, MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO), weighted median, and MR-Egger, and three MVMR approaches, including the extended versions of IVW, MR-Egger, and the Q-minimization methods, were performed to evaluate the causal effects. The robustness of the MR results was further confirmed using several sensitivity analyses.

RESULTS

UVMR revealed that a genetically predisposed per standard deviation increase in serum alanine aminotransferase (ALT) level increased the risk of diabetic retinopathy (DR) in type 2 diabetes mellitus (T2DM) (IVW OR = 1.489, 95% CI = 1.206-1.772, p = 0.006). Likewise, serum aspartate aminotransferase (AST) levels showed similar results (IVW OR = 1.376, 95% CI = 1.115-1.638, p = 0.017). Furthermore, these effects were consistent after controlling for glycemia and blood pressure using MVMR analysis. Additionally, sensitivity analyses further strengthened the causality. However, no significant associations were found between alkaline phosphatase (ALP), gamma-glutamyl transferase (GGT), and diabetic microvascular complications.

CONCLUSIONS

Robust evidence was demonstrated for an independent causal effect of serum ALT or AST concentration on the risk of DR in T2DM. Further investigations are necessary to elucidate the potential biological mechanisms and confirm their clinical significance for early prevention and intervention.

Gastroesophageal reflux disease and the risk of respiratory diseases: a Mendelian randomization study

BACKGROUND

Observational studies have suggested a suspected association between gastroesophageal reflux disease (GERD) and respiratory diseases, but the causality remains equivocal. The goal of this study was to evaluate the causal role of GERD in respiratory diseases by employing Mendelian randomization (MR) studies.

METHODS

We conducted Mendelian randomization analysis based on summary data of genome-wide association studies (GWASs) and three MR statistical techniques (inverse variance weighted, weighted median and MR-Egger) were employed to assess the probable causal relationship between GERD and the risk of respiratory diseases. Sensitivity analysis was also carried out to ensure more trustworthy results, which involves examining the heterogeneity, pleiotropy and leave-one-SNP-out method. We also identified 33 relevant genes and explored their distribution in 26 normal tissues.

RESULTS

In the analysis, for every unit increase in developing GERD, the odds ratio for developing COPD, bronchitis, pneumonia, lung cancer and pulmonary embolism rose by 72% (ORIVW = 1.72, 95% CI 1.50; 1.99), 19% (ORIVW = 1.19, 95% CI 1.11; 1.28), 16% (ORIVW = 1.16, 95% CI 1.07; 1.26), 0. 3% (ORIVW = 1.003, 95% CI 1.0012; 1.0043) and 33% (ORIVW = 1.33, 95% CI 1.12; 1.58), respectively, in comparison with non-GERD cases. In addition, neither heterogeneity nor pleiotropy was found in the study. This study also found that gene expression was higher in the central nervous system and brain tissue than in other normal tissues.

CONCLUSIONS

This study provided evidence that people who developed GERD had a higher risk of developing COPD, bronchitis, pneumonia, lung cancer and pulmonary embolism. Our research suggests physicians to give effective treatments for GERD on respiratory diseases. By exploring the gene expression, our study may also help to reveal the role played by the central nervous system and brain tissue in developing respiratory diseases caused by GERD.

The impact of age-specific childhood body-mass index on adult cardiometabolic traits: a Mendelian randomization study

Objective

To evaluate the causal relationship between childhood body-mass index (BMI) at different ages and adult cardiometabolic traits.

Methods

We retrieved genetic instrument variables (IVs) for exposures (standardized BMI at newborn, infant, toddler and late childhood), cardiometabolic traits and potential confounders or mediators (adult BMI, SHBG, testosterone and age at menarche) from the corresponding genome-wide association analysis. We performed univariate and multivariable Mendelian randomization (MR) to dissect associations between age-specific childhood BMI and adult cardiometabolic outcomes. Odds ratio was used to present the direction of the causal association.

Results

In univariate MR, higher newborn BMI was causally associated with reduced risk for type 2 diabetes in women. Late childhood BMI was associated with increased risk for female diabetes and coronary artery disease (CAD), myocardial infarction (MI), and chronic kidney disease (CKD) in general population. Among these associations, only association between late childhood BMI with MI remained significant after adjusting for adult male BMI and sex hormones, (OR = 1.120, 95% CI 1.023-1.226, p = 0.014). Besides, in multivariable MR, we found evidence for causal association between newborn BMI with reduced risk for CAD (OR = 0.862, 95% CI 0.751-0.989, p = 0.034) and MI (OR = 0.864, 95% CI 0.752-0.991, p = 0.037) in men. No obvious impact of infant or toddler BMI was identified on the above-mentioned diseases. For continuous cardiometabolic traits, in all age epochs except infant, higher BMI was associated with increased level of fasting glucose in women.

Conclusion

BMI at birth and late childhood exerts different impact on adult cardiometabolic diseases, while BMI at infant and toddler ages is not causally associated with these outcomes. The effect of childhood BMI may be influenced by sex disparities.

A Mendelian randomization study of the effect of serum 25-hydroxyvitamin D levels on autoimmune thyroid disease

Objective

The influence of vitamin D on autoimmune thyroid disease (AITD) remains a subject of ongoing debate. This study employs Mendelian randomization (MR) to investigate the causal correlations of serum 25-hydroxyvitamin D (25[OH]D) levels with autoimmune thyroiditis (AIT), autoimmune hyperthyroidism (AIH), and Graves disease (GD).

Methods

Data on single nucleotide polymorphisms related to serum 25(OH)D levels, AIT, AIH, and GD were sourced from UK Biobank and FinnGen. Inverse variance weighted, MR-Egger, and weighted median were employed to test the exposure-outcome causal relationship. Assessments of horizontal pleiotropy, heterogeneity, and stability were performed using the MR-Egger intercept, Cochran's Q test, and leave-one-out sensitivity analysis, respectively.

Results

The results of MR analysis showed increased serum 25(OH)D levels was associated with a reduced risk of AIT (OR 0.499, 95% CI 0.289 to 0.860, p = 0.012) but not causal associated with AIH (OR 0.935, 95% CI 0.695 to 1.256, p = 0.654) and GD (OR 0.813, 95% CI 0.635 to 1.040, p = 0.100). Intercept analysis showed no horizontal pleiotropy (p > 0.05), and Cochran's Q test showed no heterogeneity (p > 0.05). Sensitivity analysis suggested that these results were robust.

Conclusion

An increased serum 25(OH)D level is associated with AIT risk reduction but unrelated to AIH and GD. This finding suggests that vitamin D supplementation can be valuable for preventing and treating AIT.

Exploring the causal role of gut microbiota in giant cell arteritis: a Mendelian randomization analysis with mediator insights

Background

Giant Cell Arteritis (GCA) is a complex autoimmune condition. With growing interest in the role of gut microbiota in autoimmune diseases, this research aimed to explore the potential causal relationship between gut microbiota and GCA, and the mediating effects of specific intermediaries.

Methods

Using a bidirectional two-sample Mendelian randomization (MR) design, we investigated associations between 191 microbial taxa and GCA. A two-step MR technique discerned the significant mediators on this relationship, followed by Multivariable MR analyses to quantify the direct influence of gut microbiota on GCA and mediation effect proportion, adjusting for these mediators.

Results

Nine taxa displayed significant associations with GCA. Among them, families like Bacteroidales and Clostridiaceae1 had Odds Ratios (OR) of 1.48 (p=0.043) and 0.52 (p=5.51e-3), respectively. Genera like Clostridium sensu stricto1 and Desulfovibrio showed ORs of 0.48 (p=5.39e-4) and 1.48 (p=0.037), respectively. Mediation analyses identified 25 hydroxyvitamin D level (mediation effect of 19.95%), CD14+ CD16- monocyte counts (mediation effect of 27.40%), and CD4+ T cell counts (mediation effect of 28.51%) as significant intermediaries.

Conclusion

Our findings provide invaluable insights into the complex interplay between specific gut microbiota taxa and GCA. By highlighting the central role of gut microbiota in influencing GCA risk and long-term recurrence, and their interactions with vital immune mediators, this research paves the way for potential therapeutic interventions in GCA management.

Inverse association between type 2 diabetes and hepatocellular carcinoma in East Asian populations

Aims

To investigate the potential association between type 2 diabetes (T2D) and hepatocellular carcinoma (HCC) in East Asian populations using Mendelian randomization (MR) analyses.

Methods

Bidirectional Mendelian randomization (MR) studies were conducted using summary statistics from genome-wide association studies (GWAS) related to T2D and HCC. The potential effects of confounders such as chronic hepatitis B, chronic hepatitis C, body mass index, and alcohol intake frequency were corrected using a multivariate MR study. Various MR methods, including the inverse variance weighted (IVW) approach, were used to estimate the associations between T2D and HCC. Sensitivity analysis and assessment of heterogeneity were performed to ensure the robustness of the results.

Results

In the forward MR study, the IVW approach of MR analysis suggested an inverse association between T2D and HCC, with a risk odds ratio of 0.8628 (95% confidence interval [CI], 0.7888-0.9438). Furthermore, even after adjusting for BMI, chronic hepatitis B, and alcohol intake frequency, this study still supports the inverse association between T2D and HCC. Additional MR methods provided further support for this relationship. Sensitivity analysis and assessment of heterogeneity confirmed the robustness of the results. The reverse MR analysis did not show a clear impact of genetic liability to HCC on reduced risk of T2D(OR=0.9788; 95% CI, 0.9061-1.0574).

Conclusion

This study provides evidence of an inverse association between T2D and HCC in East Asian populations using MR analysis. Further studies are warranted to validate these findings.

Leukocyte telomere length and sarcopenia-related traits: A bidirectional Mendelian randomization study

Accumulating evidence indicated that leukocyte telomere length (LTL) was related to sarcopenia. However, it is still not clear whether the association of changes in LTL with sarcopenia is likely to be causal, or could be explained by reverse causality. Thus, we carried on bidirectional Mendelian randomization (MR) and multivariable MR analyses to identify the causal relationship between LTL and sarcopenia-related traits. Summary-level data and independent variants used as instruments came from large genome-wide association studies of LTL (472,174 participants), appendicular lean mass (450,243 participants), low grip strength (256,523 participants), and walking pace (450,967 participants). We identified suggestive association of longer LTL with larger appendicular lean mass [odds ratio (OR) = 1.053; 95% confidence interval (CI), 1.009-1.099; P = 0.018], and causal association of longer LTL with a lower risk of low grip strength (OR = 0.915; 95% CI, 0.860-0.974; P = 0.005). In the reverse MR analysis, we also observed a positive causal association between walking pace and LTL (OR = 1.252; 95% CI, 1.121-1.397; P < 0.001). Similar results can be repeated in sensitivity analyses. While in the multivariable MR analysis, the estimate of the impact of walking pace on LTL underwent a transformation after adjusting for T2DM (OR = 1.141; 95%CI: 0.989-1.317; P = 0.070). The current MR analysis supported a causal relationship between shorter telomere length and both low muscle mass and strength. Additionally, walking pace may affect LTL through T2DM.

Causal Association Between Allergic Diseases and Dementia: Evidence from Multivariate Mendelian Randomization Study

Background

The link between allergic diseases and dementia remains controversial, and the genetic causality of this link is unclear.

Objective

This study investigated the causal relationship between allergic diseases and dementia using univariate and multivariate Mendelian randomization (MR) methods.

Methods

We selected genome-wide association studies including 66,645 patients with allergic diseases and 12,281 patients with dementia, with statistical datasets derived from the FinnGen Consortium of European origin. After a rigorous screening process for single nucleotide polymorphisms to eliminate confounding effects, MR estimation was performed mainly using the inverse variance weighting method and the MR-Egger method. Sensitivity analyses were performed using Cochran's Q test, MR-PRESSO test, MR Pleiotropy residuals and leave-one-out analysis.

Results

Univariate and multivariate MR together demonstrated a causal relationship between atopic dermatitis and reduced vascular dementia (VaD) risk (OR = 0.89, 95% CI: 0.81-0.99, p = 0.031; OR = 0.85, 95% CI: 0.76-0.95, p = 0.003). MVMR confirmed asthma was associated with a reduction in the risk of Alzheimer's disease (AD) (OR = 0.82, 95% CI: 0.71-0.94, p = 0.005) and may be associated with a reduction in the risk of VaD (OR = 0.80, 95% CI: 0.65-0.99, p = 0.042); allergic rhinitis may be causally associated with an increased risk of AD (OR = 1.16, 95% CI: 1.00-1.35, p = 0.046) and VaD (OR = 1.29, 95% CI: 1.03-1.62, p = 0.027). In sensitivity analyses, these findings were reliable.

Conclusions

MR methods have only demonstrated that allergic rhinitis dementia is associated with an increased risk of developing dementia. Previously observed associations between other allergic diseases and dementia may be influenced by comorbidities and confounding factors rather than causality.

Ankylosing spondylitis status and risk of secondary systemic amyloidosis: A two-sample mendelian randomization study

OBJECTIVES

There is still controversy regarding the causal relationship between ankylosing spondylitis (AS) and secondary systemic amyloidosis (SSA). This study utilized aggregated data from genome-wide association studies (GWAS) on population cohorts to investigate whether a causal relationship exists between AS and SSA.

METHODS

The genetic causal relationship between AS status and SSA was analyzed utilizing a two-sample Mendelian randomization (TSMR). The analyses were conducted using the weighted mode method (WM2), inverse variance weighted method (IVW), simple mode (SM), weighted median method (WM1), and Mendelian randomization Egger regression (MR-Egger). Among these methods, the primary results were based on the IVW approach. The association was evaluated using the odds ratio (OR) along with a 95% confidence interval (95% CI).

RESULTS

The IVW analysis revealed a positive causal relationship between AS status and SSA (OR = 1.411, 95 % CI = 1.069, 1.862, P = 0.015). Meanwhile, the WM1 (OR = 1.394, 95 % CI = 1.115, 1.742, P = 0.004) and WM2 (OR = 1.393, 95 % CI = 1.112, 1.743, P = 0.045) methods also identified a positive causal relationship between AS status and SSA. The MR-Egger method did not identify a causal relationship between AS and SSA (OR = 1.175, 95 % CI = 0.888, 1.555, P = 0.342). The SM results demonstrated that the observed genotypes did not exhibit statistically significant differences between AS and SSA (OR = 1.184, 95 % CI = 0.416, 3.366, P = 0.767). The results of the MR-Egger regression suggested that the results were unaffected by bias caused by genetic pleiotropy (Intercept = 0.283, SE = 0.134, P = 0.126). Cochran's Q test did not reveal any significant heterogeneity (Q = 1.759, P = 0.624). The "leave-one-out" analysis further confirmed that the absence of any single SNP did not impact the robustness of our results.

CONCLUSION

This study revealed a positive causal relationship between AS status and the occurrence of SSA, providing new insights into the genetic analysis of SSA.

Causal Relationships between Air Pollutant Exposure and Bone Mineral Density and the Risk of Bone Fractures: Evidence from a Two-Stage Mendelian Randomization Analysis

A number of studies from the literature have suggested that exposure to air pollutants is associated with a declined bone mineral density (BMD), and increased risks of osteoporosis (OP) and bone fractures. This study was performed to systemically assess the genetically causal associations of air pollutants with site-/age-specific BMD and risk of bone fractures with the implementation of two-sample Mendelian randomization (TSMR) and multivariate Mendelian randomization (MVMR). The TSMR analysis was implemented to infer the causal associations between air pollutants and BMD and the risk of bone fractures, additional MVMR analysis was used to further estimate the direct causal effects between air pollutants and BMD, the occurrence of OP, and bone fractures. The results showed that NOx exposure contributed to lower femoral neck BMD (FN-BMD) (β = -0.71, 95%CI: -1.22, -0.20, p = 0.006) and total body BMD (TB-BMD) (β = -0.55, 95%CI: -0.90, -0.21, p = 0.002). Additionally, exposure to PM10 was found to be associated with a decreased TB-BMD (B β = -0.42, 95%CI: -0.66, -0.18, p = 0.001), further age-specific subgroup analysis demonstrated the causal effect of PM10 exposure on the decreased TB-BMD in a subgroup aged 45 to 60 years (β = -0.70, 95%CI: -1.12, -0.29, p = 0.001). Moreover, the findings of the MVMR analysis implied that there was a direct causal effect between PM10 exposure and the decreased TB-BMD (45 < age < 60), after adjusting for PM2.5 and PM2.5 -10 exposure. Our study provides additional evidence to support the causal associations of higher concentrations of air pollutant exposure with decreased BMD, especially in those populations aged between 45 to 60 years, suggesting that early intervention measures and public policy should be considered to improve public health awareness and promote bone health.

Association between type 1 diabetes mellitus and ankylosing spondylitis: a two-sample Mendelian randomization study

Objective

The development of ankylosing spondylitis (AS) is closely related to autoimmune system dysfunction. Type 1 diabetes mellitus (T1DM) is an autoimmune disease that is a risk factor for many diseases. This study aimed to investigate the causal relationship between T1DM mellitus and AS genetically.

Methods

A genome-wide association study (GWAS) of causal relationships between exposure (T1DM) and outcome (AS) was performed using summary data from the GWAS database. We conducted a two-sample Mendelian randomization (MR) study of these two diseases. Inverse variance weighting (IVW) was used as the primary analysis method, with MR Egger, weighted median, and weighted mode used as supplementary methods. Sensitivity analyses were performed using Cochran's Q test, MR-Egger intercept, MR-Pleiotropy RESidual Sum and outlier methods, leave-one-out analysis, and funnel plots.

Results

A total of 11 single nucleotide polymorphisms (SNPs)were identified for instrumental variables(IVs) for MR analysis.IVW found that T1DM was causally associated with AS ((IVW: OR = 1.0006 (95% CI 1.0001, 1.0011), p = 0.0057; MR-Egger: OR = 1.0003 (95% CI 0.9995, 1.0012), p = 0.4147; weighted median: OR = 1.0006 (95% CI 1.0003, 1.0008), p = 0.0001; weighted mode: OR = 1.0007 (95% CI 1.0005, 1.0009), p = 0.0001). No horizontal pleiotropy was found for the MR-Egger intercept, and leave -one-out analysis found that the results remained stable after the removal of individual SNPs.

Conclusion

The results of the two-sample MR analysis supported a causal relationship between T1DM and AS risk.

Causal effects of potential risk factors on postpartum depression: a Mendelian randomization study

Background

Postpartum depression (PPD) is a type of depressive episode related to parents after childbirth, which causes a variety of symptoms not only for parents but also affects the development of children. The causal relationship between potential risk factors and PPD remains comprehensively elucidated.

Methods

Linkage disequilibrium score regression (LDSC) analysis was conducted to screen the heritability of each instrumental variant (IV) and to calculate the genetic correlations between effective causal factors and PPD. To search for the causal effect of multiple potential risk factors on the incidence of PPD, random effects of the inverse variance weighted (IVW) method were applied. Sensitivity analyses, including weighted median, MR-Egger regression, Cochrane's Q test, and MR Pleiotropy Residual Sum and Outlier (MR-PRESSO), were performed to detect potential Mendelian randomization (MR) assumption violations. Multivariable MR (MVMR) was conducted to control potential multicollinearity.

Results

A total of 40 potential risk factors were investigated in this study. LDSC regression analysis reported a significant genetic correlation of potential traits with PPD. MR analysis showed that higher body mass index (BMI) (Benjamini and Hochberg (BH) corrected p = 0.05), major depression (MD) (BH corrected p = 5.04E-19), and schizophrenia (SCZ) (BH corrected p = 1.64E-05) were associated with the increased risk of PPD, whereas increased age at first birth (BH corrected p = 2.11E-04), older age at first sexual intercourse (BH corrected p = 3.02E-15), increased average total household income before tax (BH corrected p = 4.57E-02), and increased years of schooling (BH corrected p = 1.47E-11) led to a decreased probability of PPD. MVMR analysis suggested that MD (p = 3.25E-08) and older age at first birth (p = 8.18E-04) were still associated with an increased risk of PPD.

Conclusion

In our MR study, we found multiple risk factors, including MD and younger age at first birth, to be deleterious causal risk factors for PPD.

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 for cardiovascular diseases: principles and applications

Large-scale genome-wide association studies conducted over the last decade have uncovered numerous genetic variants associated with cardiometabolic traits and risk factors. These discoveries have enabled the Mendelian randomization (MR) design, which uses genetic variation as a natural experiment to improve causal inferences from observational data. By analogy with the random assignment of treatment in randomized controlled trials, the random segregation of genetic alleles when DNA is transmitted from parents to offspring at gamete formation is expected to reduce confounding in genetic associations. Mendelian randomization analyses make a set of assumptions that must hold for valid results. Provided that the assumptions are well justified for the genetic variants that are employed as instrumental variables, MR studies can inform on whether a putative risk factor likely has a causal effect on the disease or not. Mendelian randomization has been increasingly applied over recent years to predict the efficacy and safety of existing and novel drugs targeting cardiovascular risk factors and to explore the repurposing potential of available drugs. This review article describes the principles of the MR design and some applications in cardiovascular epidemiology.

Genetic insights into the association of statin and newer nonstatin drug target genes with human longevity: a Mendelian randomization analysis

BACKGROUND

It remains controversial whether the long-term use of statins or newer nonstatin drugs has a positive effect on human longevity. Therefore, this study aimed to investigate the genetic associations between different lipid-lowering therapeutic gene targets and human longevity.

METHODS

Two-sample Mendelian randomization analyses were conducted. The exposures comprised genetic variants that proxy nine drug target genes mimicking lipid-lowering effects (LDLR, HMGCR, PCKS9, NPC1L1, APOB, CETP, LPL, APOC3, and ANGPTL3). Two large-scale genome-wide association study (GWAS) summary datasets of human lifespan, including up to 500,193 European individuals, were used as outcomes. The inverse-variance weighting method was applied as the main approach. Sensitivity tests were conducted to evaluate the robustness, heterogeneity, and pleiotropy of the results. Causal effects were further validated using expression quantitative trait locus (eQTL) data.

RESULTS

Genetically proxied LDLR variants, which mimic the effects of lowering low-density lipoprotein cholesterol (LDL-C), were associated with extended lifespan. This association was replicated in the validation set and was further confirmed in the eQTL summary data of blood and liver tissues. Mediation analysis revealed that the genetic mimicry of LDLR enhancement extended lifespan by reducing the risk of major coronary heart disease, accounting for 22.8% of the mediation effect. The genetically proxied CETP and APOC3 inhibitions also showed causal effects on increased life expectancy in both outcome datasets. The lipid-lowering variants of HMGCR, PCKS9, LPL, and APOB were associated with longer lifespans but did not causally increase extreme longevity. No statistical evidence was detected to support an association between NPC1L1 and lifespan.

CONCLUSION

This study suggests that LDLR is a promising genetic target for human longevity. Lipid-related gene targets, such as PCSK9, CETP, and APOC3, might potentially regulate human lifespan, thus offering promising prospects for developing newer nonstatin therapies.

Hypertension linked to Alzheimer's disease via stroke: Mendelian randomization

This study aimed to investigate the relationship between hypertension and Alzheimer's disease (AD) and demonstrate the key role of stroke in this relationship using mediating Mendelian randomization. AD, a neurodegenerative disease characterized by memory loss, cognitive impairment, and behavioral abnormalities, severely affects the quality of life of patients. Hypertension is an important risk factor for AD. However, the precise mechanism underlying this relationship is unclear. To investigate the relationship between hypertension and AD, we used a mediated Mendelian randomization method and screened for mediating variables between hypertension and AD by setting instrumental variables. The results of the mediated analysis showed that stroke, as a mediating variable, plays an important role in the causal relationship between hypertension and AD. Specifically, the mediated indirect effect value for stroke obtained using multivariate mediated MR analysis was 54.9%. This implies that approximately 55% of the risk of AD owing to hypertension can be attributed to stroke. The results suggest that the increased risk of AD owing to hypertension is mediated through stroke. The finding not only sheds light on the relationship between hypertension and AD but also indicates novel methods for the prevention and treatment of AD. By identifying the critical role of stroke in the link between hypertension and AD, this study provides insights into potential interventions that could mitigate the impact of hypertension on AD. This could help develop personalized treatments and help improve the quality of life of patients with AD who suffer from hypertension.

The mediating role of atrial fibrillation in causal associations between risk factors and stroke: a Mendelian randomization study

OBJECTIVES

Atrial fibrillation (AF) contributes to stroke development and progression. We aimed to quantify the mediating role of AF in the causal associations between a wide range of risk factors and stroke via a Mendelian randomization (MR) framework.

METHODS

We assessed the associations of 108 traits with stroke and its subtypes in a 2-sample univariable MR approach, then conducted a bidirectional MR analysis between these 108 traits and AF to evaluate the presence and direction of their causal associations. Finally, to further investigate the extent to which AF mediated the effects of eligible traits on stroke, we applied multivariable and 2-step MR techniques in a mediation analysis where outcomes were restricted to stroke types causally affected by AF (any stroke [AS], any ischemic stroke [AIS], and cardioembolic stroke [CES]).

RESULTS

Among 108 traits, 42 were putatively causal for at least 1 stroke type; of these 42 traits, 20 that had no bidirectional relationship with AF were retained. Finally, 33 associations of 15 eligible traits were examined in the mediation analysis. The mediation analyses for AS, AIS, and CES each included 11 eligible traits. After AF adjustment, the direct effects of all traits on CES were attenuated to null (all p>0.05), while the associations with AS and AIS persisted for most traits (AF-mediated proportion: from 6.6% [95% confidence interval, 2.7 to 0.6] to 52.0% [95% confidence interval, 39.8 to 64.3]).

CONCLUSIONS

The causal associations between all eligible traits and CES were largely mediated through AF, while most traits affected AS and AIS independently of AF.

Investigating linear and nonlinear associations of LDL cholesterol with incident chronic kidney disease, atherosclerotic cardiovascular disease and all-cause mortality: A prospective and Mendelian randomization study

BACKGROUND AND AIMS

Observational studies suggest potential nonlinear associations of low-density lipoprotein cholesterol (LDL-C) with cardio-renal diseases and mortality, but the causal nature of these associations is unclear. We aimed to determine the shape of causal relationships of LDL-C with incident chronic kidney disease (CKD), atherosclerotic cardiovascular disease (ASCVD) and all-cause mortality, and to evaluate the absolute risk of adverse outcomes contributed by LDL-C itself.

METHODS

Observational analysis and one-sample Mendelian randomization (MR) with linear and nonlinear assumptions were performed using the UK Biobank of >0.3 million participants with no reported prescription of lipid-lowering drugs. Two-sample MR on summary-level data from the Global Lipid Genetics Consortium (N = 296,680) and the CKDGen (N = 625,219) was employed to replicate the relationship for kidney traits. The 10-year probabilities of the outcomes was estimated by integrating the MR and Cox models.

RESULTS

Observationally, participants with low LDL-C were significantly associated with a decreased risk of ASCVD, but an increased risk of CKD and all-cause mortality. Univariable MR showed an inverse total effect of LDL-C on incident CKD (HR [95% CI]:0.84 [0.73-0.96]; p = 0.011), a positive effect on ASCVD (1.41 [1.29-1.53]; p<0.001), and no significant causal effect on all-cause mortality. Multivariable MR, controlling for high-density lipoprotein cholesterol (HDL-C) and triglycerides, identified a positive direct effect on ASCVD (1.32 [1.18-1.47]; p<0.001), but not on CKD and all-cause mortality. These results indicated that genetically predicted low LDL-C had an inverse indirect effect on CKD mediated by HDL-C and triglycerides, which was validated by a two-sample MR analysis using summary-level data from the Global Lipid Genetics Consortium (N = 296,680) and the CKDGen consortium (N = 625,219). Suggestive evidence of a nonlinear causal association between LDL-C and CKD was found. The 10-year probability curve showed that LDL-C concentrations below 3.5 mmol/L were associated with an increased risk of CKD.

CONCLUSIONS

In the general population, lower LDL-C was causally associated with lower risk of ASCVD, but appeared to have a trade-off for an increased risk of CKD, with not much effect on all-cause mortality. LDL-C concentration below 3.5 mmol/L may increase the risk of CKD.

The causal effects between selenium levels and the brain cortical structure: A two-sample Mendelian randomization study

Extensive research has demonstrated the critical role of selenium (Se) and selenoproteins in brain function and cognition. However, the impact of Se on brain cortical structure remains enigmatic. Therefore, this study used Mendelian randomization (MR) analysis to investigate the causal effect between Se levels and brain cortical structure.

METHODS

This study utilizes 11 genetic variants associated with Se level variations, extracted from a large-scale genome-wide association study (GWAS) encompassed circulating Se levels (n = 5477) and toenail Se levels (n = 4162) in the European population. Outcome data were sourced from the summary statistics of the ENIGMA Consortium, comprising GWAS data from 51,666 individuals. The variables include cortical surface area (SA), thickness (TH) at the global level, and 34 functional cortical regions evaluated by magnetic resonance imaging. The inverse-variance-weighted method was used as the primary estimate. Additionally, sensitivity analyses were conducted to detect potential violations of assumptions underlying MR.

RESULTS

At the global level, Se levels were not correlated with SA but showed a significant negative correlation with TH (β = -0.00485 mm, SE = 0.00192, p = .0115). Heterogeneity was observed across different brain regions, with positive correlations found between Se levels and the TH of the parahippocampal gyrus, superior frontal gyrus, and frontal pole, whereas negative correlations were found with the TH of the inferior parietal lobe and middle temporal lobe. Regarding SA, Se levels exhibit positive correlations with the pars triangularis, caudal anterior cingulate, inferior parietal lobe, and banks of the superior temporal sulcus. Conversely, negative correlations were observed with the medial orbitofrontal cortex, posterior cingulate gyrus, insula, and the middle, superior, and transverse gyrus of the temporal lobe. No pleiotropy was detected.

RESULTS

This MR study indicated that Se levels causally influence the brain cortical structure.

The causal association between smoking, alcohol consumption and risk of upper urinary calculi: insights from a Mendelian randomization study

Background: The causal link between smoking, alcohol consumption, and upper urinary calculi remains uncertain in observational studies due to confounding factors. To uncover potential causal associations, we utilized two-sample univariable and multivariable Mendelian randomization (MR) methods. Methods: Five risk factors related to lifestyles (cigarettes per day, lifetime smoking index, smoking initiation, drinks per week and alcohol intake frequency) were chosen from the Genome-Wide Association Study (GWAS). Upper urinary calculi were obtained from the FinnGen and United Kingdom Biobank consortium. Inverse-variance-weighted (IVW) was mainly used to compute odds ratios (OR) and 95% confidence intervals (Cl). While diligently scrutinizing potential sources of heterogeneity and horizontal pleiotropy via the rigorous utilization of Cochran's Q test, the MR-PRESSO method, and MR-Egger. Results: The summary OR for upper urinary calculi was 0.6 (IVW 95% CI: 0.49-0.74; p = 1.31 × 10-06) per standard deviation decrease in drinks per week. Interestingly, the genetically predicted alcohol intake frequency was associated with a significantly increased risk upper urinary calculi (OR = 1.27; 95% CI: 1.11-1.45; p = 0.0005). Our study found no association between smoking initiation, the number of cigarettes per day, and the lifetime smoking index and the risk of upper urinary calculi. By adjusting for body mass index and education, estimates of drinks per week remained consistent in multivariate MR analyses, while alcohol intake frequency became non-significant. Conclusion: MR analysis showed that drinks per week was negatively associated with upper urinary calculi, whereas the effect of tobacco on upper urinary calculi was not significant and the detrimental effect of alcohol intake frequency on upper urinary calculi became non-significant after adjusting for BMI and education.

Mitochondrial DNA copy number is associated with Crohn's disease: a comprehensive Mendelian randomization analysis

Mitochondrial DNA plays a critical role in the pathophysiological process of inflammation. However, the relationship between mitochondrial DNA copy number (mtDNA-CN) and inflammatory bowel diseases (IBD) remains poorly understood. We conducted a comprehensive Mendelian randomization (MR) using three instrumental variables (IVs) to explore the causal associations between mtDNA-CN and IBD, including Crohn's disease (CD), ulcerative colitis (UC). MR-Egger regression, weighted median, inverse-variance weighted (IVW), and weighted mode methods were used to evaluate the potential causal associations. The robustness of the IVW estimates was determined using the leave-one-out sensitivity test. A meta-analysis was conducted to pool the results from the three sets of IVs. Upon analysis, the findings of the current study revealed that genetically predicted mtDNA-CN was not associated with IBD (CD + UC) and UC. The results of MR analyses between mtDNA-CN and CD risk were inconsistent by using three sets of IVs. After a meta-analysis, we found that genetically predicted mtDNA-CN was associated with CD risk (odds ratio = 2.09; 95% confidence interval: 1.37-3.18). This finding was also confirmed by multivariable MR analyses and remained robust when tested with the leave-one-out sensitivity test. In conclusion, genetically predicted mtDNA-CN was found to be associated with CD risk. Therefore, mtDNA levels in the blood could potentially be used as a marker for CD risk assessment. Further studies are needed to elucidate the underlying mechanisms and validate the results of this study.

Evidence of shared genetic factors in the etiology of gastrointestinal disorders and endometriosis and clinical implications for disease management

In clinical practice, the co-existence of endometriosis and gastrointestinal symptoms is often observed. Using large-scale datasets, we report a genetic correlation between endometriosis and irritable bowel syndrome (IBS), peptic ulcer disease (PUD), gastro-esophageal reflux disease (GORD), and a combined GORD/PUD medicated (GPM) phenotype. Mendelian randomization analyses support a causal relationship between genetic predisposition to endometriosis and IBS and GPM. Identification of shared risk loci highlights biological pathways that may contribute to the pathogenesis of both diseases, including estrogen regulation and inflammation, and potential therapeutic drug targets (CCKBR; PDE4B). Higher use of IBS, GORD, and PUD medications in women with endometriosis and higher use of hormone therapies in women with IBS, GORD, and PUD, support the co-occurrence of these conditions and highlight the potential for drug repositioning and drug contraindications. Our results provide evidence of shared disease etiology and have important clinical implications for diagnostic and treatment decisions for both diseases.

Interleukin-17 and inflammatory bowel disease: a 2-sample Mendelian randomization study

Introduction

Observational studies have discovered a contradictory phenomenon between interleukin-17 (IL-17) and inflammatory bowel disease (IBD). The study aimed to confirm the causal association between each subtype of IL-17 and IBD.

Methods

We performed a 2-sample univariable and multivariable mendelian randomization (MR) to determine which subtype of IL-17 is causally related to IBD and its subtypes, and used a series of sensitivity analysis to examine the reliability of the main MR assumptions.

Results

We found that IL-17B, IL-17E and IL-17RB were significantly associated with an increased risk of UC (IL-17B: OR: 1.26, 95% CI, 1.09-1.46, P < 0.01; IL-17E: OR: 1.17, 95% CI, 1.05-1.30, P < 0.01; IL-17RB: OR: 1.30, 95% CI, 1.20-1.40, P < 0.0001) while IL-17C and IL-17RC showed causal effects on the increased risk of CD (IL-17C: OR: 1.23, 95% CI, 1.21-1.26, P < 0.0001; IL-17RC: OR: 2.01, 95% CI, 1.07-3.75, P=0.03). The results of multivariable MR (MVMR) showed that the causal effects of IL-17B and IL-17E on UC were unilaterally dependent on IL-17RB, while the effects of IL-17C and IL-17RC on CD were interdependent.

Discussion

Our study provided new genetic evidence for the causal relationships between each subtype of IL-17 and IBD, promoting future mechanistic research in IBD.

Interplay Between Chronic Kidney Disease, Hypertension, and Stroke: Insights From a Multivariable Mendelian Randomization Analysis

BACKGROUND AND OBJECTIVES

Chronic kidney disease (CKD) increases the risk of stroke, but the extent through which this association is mediated by hypertension is unknown. We leveraged large-scale genetic data to explore causal relationships between CKD, hypertension, and cerebrovascular disease phenotypes.

METHODS

We used data from genome-wide association studies of European ancestry to identify genetic proxies for kidney function (CKD diagnosis, estimated glomerular filtration rate [eGFR], and urinary albumin-to-creatinine ratio [UACR]), systolic blood pressure (SBP), and cerebrovascular disease (ischemic stroke and its subtypes and intracerebral hemorrhage). We then conducted univariable, multivariable, and mediation Mendelian randomization (MR) analyses to investigate the effect of kidney function on stroke risk and the proportion of this effect mediated through hypertension.

RESULTS

Univariable MR revealed associations between genetically determined lower eGFR and risk of all stroke (odds ratio [OR] per 1-log decrement in eGFR, 1.77; 95% CI 1.31-2.40; p < 0.001), ischemic stroke (OR 1.81; 95% CI 1.31-2.51; p < 0.001), and most strongly with large artery stroke (LAS) (OR 3.00; 95% CI 1.33-6.75; p = 0.008). These associations remained significant in the multivariable MR analysis, controlling for SBP (OR 1.98; 95% CI 1.39-2.82; p < 0.001 for all stroke; OR 2.16; 95% CI 1.48-3.17; p < 0.001 for ischemic stroke; OR 4.35; 95% CI 1.84-10.27; p = 0.001 for LAS), with only a small proportion of the total effects mediated by SBP (6.5% [0.7%-16.8%], 6.6% [0.8%-18.3%], and 7.2% [0.5%-24.8%], respectively). Total, direct and indirect effect estimates were similar across a number of sensitivity analyses (weighted median, MR-Egger regression).

DISCUSSION

Our results demonstrate an independent causal effect of impaired kidney function, as assessed by decreased eGFR, on stroke risk, particularly LAS, even when controlled for SBP. Targeted prevention of kidney disease could lower atherosclerotic stroke risk independent of hypertension.

Relationship Between the Ovarian Cyst and Depression: A Two-Sample Mendelian Randomization Study

Objective

Explore the causal relationship between the ovarian cyst and depression using a two-sample Mendelian randomization approach (MR).

Methods

Based on data pooled from genome-wide association studies, genetic variants of the ovarian cyst and depression were selected as instrumental variables, as well as the Mendelian randomization analysis was conducted using inverse variance weighted (IVW) as the main analysis method and MR-Egger regression analysis, MR-PRESSO and other sensitivity analysis methods as supplements.

Results

The IVW analysis showed a direct causal association between ovarian cysts and depression (OR=1.040; 95% CI: 1.003, 1.078; p=0.031). Meantime, there was a causal effect of genetically predicted depression on ovarian cysts (OR=1.327.; 95% CI: 1.197, 1.470; p<0.001). Sensitivity analyses such as MR-Egger regression analysis and MR-PRESSO indicated that the IVW results were robust and reliable.

Conclusion

This study suggested since ovarian cysts and female depression are mutually causal, the comorbidity of ovarian cysts and depression in women should be actively attended to and given appropriate prevention and treatment besides the diagnosis and treatment of depression or ovarian cysts.

Cannabis use and the risk of primary open-angle glaucoma: a Mendelian randomization study

Several observational studies have investigated the association between cannabis use and intraocular pressure, but its association with primary open-angle glaucoma (POAG) remains unclear. In this study, we leveraged human genetic data to assess through Mendelian randomization (MR) whether cannabis use affects POAG. We used five single-nucleotide polymorphisms (SNPs) associated with lifetime cannabis use (P-value < 5 × 10-8) from a genome-wide association study (GWAS) (N = 184,765) by the International Cannabis Consortium, 23andMe, and UK Biobank and eleven SNPs associated with cannabis use disorder (P-value < 5 × 10-7) from a GWAS meta-analysis of (17,068 cases and 357,219 controls of European descent) from Psychiatric Genomics Consortium Substance Use Disorders working group, Lundbeck Foundation Initiative for Integrative Psychiatric Research, and deCode. We associated the selected five SNPs from the GWAS of lifetime cannabis use and the eleven SNPs from the GWAS of cannabis use disorder, with the largest to date GWAS meta-analysis of POAG (16,677 cases and 199,580 controls). MR analysis suggested no evidence for a causal association of lifetime cannabis use and cannabis use disorder with POAG (odds ratio (OR) of outcome per doubling of the odds of exposure (95% confidence interval): 1.04 (0.88; 1.23) for lifetime cannabis use and 0.97 (0.92; 1.03) for cannabis use disorder). Sensitivity analyses to address pleiotropy and weak instrument bias yielded similar estimates to the primary analysis. In conclusion, our results do not support a causal association between cannabis use and POAG.

Causal association of rheumatoid arthritis with obstructive lung disease: Evidence from Mendelian randomization study

BACKGROUND

Observational studies have found an association between rheumatoid arthritis (RA) and risk of obstructive lung disease (ORDs). However, whether RA plays a role in ORDs development remains unclear.

OBJECTIVES

This study aimed to explore the causal association of RA with ORDs.

METHODS

Both univariable and multivariable Mendelian randomization (MR) analyses were employed. Summary statistics for RA were obtained from the genome-wide association study (GWAS) meta-analysis, and the GWAS data source of ORDs, including the chronic obstructive pulmonary disease (COPD) and asthma, was accessed from the FinnGen Biobank. Causal Analysis Using Summary Effect Estimates (CAUSE) method was used to improve statistical power. multivariable and two-step mediation MR was applied to calculate the independent and mediated effects.

RESULTS

The causal estimates by univariable and CAUSE results indicated genetic predisposition to RA had an effect on the increased risk of asthma/COPD (A/C) (ORCAUSE = 1.03; 95% CI: 1.02-1.04), COPD/asthma related infections (ACI) (ORCAUSE = 1.02; 95% CI: 1.01-1.03) and COPD/asthma related pneumonia or pneumonia derived septicemia (ACP) (ORCAUSE = 1.02; 95% CI: 1.01-1.03). Genetic predisposition to RA was significantly associated with early onset COPD (ORCAUSE = 1.02; 95% CI: 1.01-1.03) and asthma (ORCAUSE = 1.02; 95% CI: 1.01-1.03) risk and suggestively associated with non-allergic asthma (nAA) risk. After adjustment for confounders, independent causal effects remained for the associations of RA with risk of A/C, ACI, and ACP, as well as COPD, early-onset COPD, and asthma [total, nAA and allergic asthma (AA)] risk. Mediation analyses revealed no potential mediator.

CONCLUSION

This study indicates a causal effect of increased genetic predisposition to RA on an increased risk of ORDs, including COPD and asthma, especially early-onset COPD and nAA, and on asthma/COPD related infections, pneumonia or pneumonia derived septicemia.

Abnormal brain structure in atopic dermatitis: Evidence from Mendelian randomization study

BACKGROUND

Structural abnormalities in the brain of patients with atopic dermatitis (AD) have been reported; however, the cause has not been determined yet. Herein, we used Mendelian randomization (MR) to reveal the causal effect of AD on brain structure.

METHODS

This study utilized summary statistics from genome-wide association studies (GWASs) to investigate a collection of cerebral structural measures, encompassing cortical thickness (CT), cortical surface area (CA), and subcortical volumes in T1 images. A comprehensive GWAS meta-analysis identified a total of 20 independent single nucleotide polymorphisms linked to AD, surpassing the genome-wide significance threshold (p < 5 × 10⁻⁸). MR estimates were aggregated through the application of the inverse variance weighted method. Additional complementary analyses (i.e., MR-Egger and weighted median approaches) were conducted to further assess the robustness of the obtained results. Sensitivity analysis and multivariate MR (MVMR) while adjusting for brain structural changes risk factors (i.e., depression and anxiety) were performed to assess the reliability and stability of observed causality.

RESULTS

Genetically determined AD exhibited a causal link with reduced caudate volumes (IVW-MR: β = -0.186, p = 0.001, p-corrected = 0.009). Furthermore, we identified potential causal associations between AD and reduced CT in the cingulate region (posterior cingulate, IVW-MR: β = -0.065, p = 0.018, p-corrected = 0.551; isthmus cingulate, IVW-MR: β = -0.086, p = 0.003, p-corrected = 0.188), as well as abnormal cortical surface area (CA) in the supramarginal (IVW-MR: β = -0.047, p = 0.044, p-corrected = 0.714) and isthmus cingulate (IVW-MR: β = 0.053, p = 0.018, p-corrected = 0.714). Additional supplementary analyses yielded consistent outcomes. There was no evidence of horizontal pleiotropy. MVMR analysis showed that the causal effects of AD on abnormal brain structure remained significant while adjusting for depression and anxiety.

CONCLUSION

This MR study provided suggestive evidence that decreased caudate nucleus, posterior cingulate cortex, isthmus cingulate cortex and supramarginal gyrus are suggestively associated with higher AD risk. Future investigation into the brain regions is recommended, which helps to clarify the underlying mechanisms and point to new therapies against AD.

The association of metabolic syndrome with rotator cuff tendinopathy: a two-sample Mendelian randomization study

BACKGROUND

Observational research reported the underlying correlation of metabolic syndrome (MetS) and its components with rotator cuff tendinopathy (RCT), but their causality remained unclear. This study aimed to investigate whether genetically predicted MetS was related to the risk of RCT.

METHODS

Both univariable and multivariable Mendelian randomization (MR) analysis was applied using summary-level data from the most comprehensive genome-wide association studies to estimate the associations of MetS and its component with RCT, with the inverse variance weighted (IVW) as the primary method, and the method of Causal Analysis Using Summary Effect Estimates (CAUSE) as a supplement for false positives detection. The mediation analysis was furtherly used for the assessment of direct and indirect effects.

RESULTS

Univariable analysis revealed that genetically predicted MetS (OR: 1.0793; 95% CI 1.0311 to 1.1297), body mass index (BMI) (OR 1.2239; 95% CI 1.1357 to 1.3189), and waist circumference (WAC) (OR 1.3177; 95% CI 1.2015 to 1.4451) had a significant positive association with the risk of RCT. Triglycerides and systolic blood pressure were suggestively associated with RCT risk. These associations were also identified by CAUSE. There was independent causality of BMI (OR: 1.1806; 95% CI 1.0788 to 1.2920) and WAC (OR 1.3716; 95% CI 1.2076 to 1.5580) on RCT after adjustment for confounders. No mediator was found in the causal associations.

CONCLUSION

Our study revealed the genetic causality of MetS and its components, especially BMI and WAC, with RCT risk. Early prevention and diagnosis of excess central adiposity contributing to MetS are significant in the RCT risk management.

The associations of alcoholic liver disease and nonalcoholic fatty liver disease with bone mineral density and the mediation of serum 25-Hydroxyvitamin D: A bidirectional and two-step Mendelian randomization

BACKGROUND

Reduced bone mineral density (BMD) and osteoporosis are common in chronic liver diseases. However, the causal effect of alcoholic liver disease (ALD) and non-alcoholic fatty liver disease (NAFLD) on BMD remains uncertain.

OBJECTIVES

This study uses a two-sample Mendelian randomization (MR) design to evaluate the genetically predicted effect of ALD and NAFLD on BMDs using summary data from publically available genome-wide association studies (GWASs).

METHODS

The GWAS summary statistics of ALD (1416 cases and 213,592 controls) and NAFLD (894 cases and 217,898 controls) were obtained from the FinnGen consortium. BMDs of four sites (total body, n = 56,284; femoral neck, n = 32,735; lumbar spine, n = 28,498; forearm, n = 8143) were from the GEnetic Factors for OSteoporosis Consortium. Data for alcohol consumption (n = 112,117) and smoking (n = 33,299) and serum 25-Hydroxyvitamin D (25-OHD) level (n = 417,580) were from UK-biobank. We first performed univariate MR analysis with the Inverse Variance Weighted (IVW) method as the primary analysis to investigate the genetically predicted effect of ALD or NAFLD on BMD. Then, multivariate MR and mediation analysis were performed to identify whether the effect was mediated by alcohol consumption, smoking, or serum 25-OHD level.

RESULTS

The MR results suggested a robust genetically predicted effect of ALD on reduced BMD in the femoral neck (FN-BMD) (IVW beta = -0.0288; 95% CI: -0.0488, -0.00871; P = 0.00494) but not the other three sites. Serum 25-OHD level exhibited a significant mediating effect on the association between ALD and reduced FN-BMD albeit the proportion of mediation was mild (2.21%). No significant effects of NAFLD, alcohol consumption, or smoking on BMD in four sites, or reverse effect of BMD on ALD or NAFLD were detected.

CONCLUSION

Our findings confirm the genetically predicted effect of ALD on reduced FN-BMD, and highlight the importance of periodic BMD and serum 25-OHD monitoring and vitamin D supplementation as needed in patients with ALD. Future research is required to validate our results and investigate the probable underlying mechanisms.

Do Obesity-Related Traits Affect Prostate Cancer Risk through Serum Testosterone? A Mendelian Randomization Study

OBJECTIVE

This study aimed to investigate whether testosterone mediates or confounds the effect of obesity-related traits on prostate cancer (PCa) using Mendelian randomization (MR) analysis.

MATERIALS AND METHODS

Data of obesity-related traits (body mass index [BMI], waist-to-hip ratio [WHR], and waist-to-hip ratio adjusted for body mass index [WHRadjBMI]) were obtained from up to 806,834 people of European ancestry; data of testosterone (bioavailable testosterone [BT], total testosterone [TT], and sex hormone-binding globulin [SHBG]) were extracted from up to 194,453 participants in the UK Biobank; and the summary-level data of PCa (79,194 cases and 61,112 controls) were obtained from the PRACTICAL consortium.

RESULT

The results supported the causal relationship between higher BMI and a reduced risk of PCa (OR = 0.91, 95% confidence interval [CI]: 0.86-0.96). Furthermore, increased BT levels were associated with an elevated risk of PCa (OR = 1.15, 95% CI: 1.06-1.24). Importantly, our analysis revealed a unidirectional causal effect-higher BMI was linked to lower BT levels (beta = -0.27, 95% CI: -0.3--0.24), but not the other way around. This suggests that BT may mediate the effect of BMI on PCa rather than confound it. Our multivariable MR results further demonstrated that considering BT as a mediator led to the weakening of BMI's effect on PCa risk (OR = 0.97, 95% CI: 0.90-1.05), while the impact of BT on PCa remained unchanged when accounting for BMI. Moreover, we identified a significant indirect effect of BMI on PCa risk (OR = 0.96, 95% CI: 0.94-0.98).

CONCLUSION

Our study provided genetic evidence that serum BT can mediate the effect of BMI on the risk of PCa, indicating the possible mechanism by which obesity reduces PCa risk.

Evaluating the effect of body mass index and 25-hydroxy-vitamin D level on basal cell carcinoma using Mendelian randomization

Basal cell carcinoma (BCC) is the most common cancer with a rising incidence among white-skinned individuals. A number of epidemiological studies have suggested that obesity and serum 25-hydroxy-vitamin D (25(OH)D) levels may affect the arising of BCC. To address this, we selected 443 and 96 single nucleotide polymorphisms (SNPs) associated with body mass index (BMI) and serum level of 25(OH)D from large-scale genome-wide association studies (GWAS), respectively. The univariable and multivariable two-sample Mendelian randomization (MR) analyses were conducted with a series of sensitivity analyses to ensure the results were reliable and reproducible. The results of univariable two-sample MR analysis showed that higher BMI was related to lower risk for BCC (Odds ratio(OR) = 0.90; 95% confidence interval (CI),[0.81,0.99]; p = 0.02). In addition, this causal effect of BMI on BCC still remained (OR = 0.88; 95%CI,[- 0.22, - 0.03], p-value = 0.008) after adjusting for 25(OH)D level in the multivariable MR analysis. However, the results suggested that 25(OH)D level was not associated with BCC(OR = 1.02; 95%CI, [0.94,1.09], p-value = 0.67). In conclusion, similar to the conclusions of retrospective observational studies, the MR results indicate that high BMI is an independent protective factor for BCC. Meanwhile, vitamin D levels may not be causally associated with the risk of basal cell carcinoma and increasing vitamin D supplementation is unlikely to reduce the risk.