Genetic causal assessment between major depression and hypertension: A two-sample bidirectional Mendelian randomization study

Background:To understand a causal role of modifiable lifestyle factors in angiotensin-converting enzyme 2 (ACE2) expression (a putative severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2] receptor) across 44 human tissues/organs, and in coronavirus disease 2019 (COVID-19) susceptibility and severity, we conducted a phenome-wide two-sample Mendelian randomization (MR) study.Methods:More than 500 genetic variants were used as instrumental variables to predict smoking and alcohol consumption. Inverse-variance weighted approach was adopted as the primary method to estimate a causal association, while MR-Egger regression, weighted median, and MR pleiotropy residual sum and outlier (MR-PRESSO) were performed to identify potential horizontal pleiotropy.Results:We found that genetically predicted smoking intensity significantly increased ACE2 expression in thyroid (β=1.468, p=1.8×10−8), and increased ACE2 expression in adipose, brain, colon, and liver with nominal significance. Additionally, genetically predicted smoking initiation significantly increased the risk of COVID-19 onset (odds ratio=1.14, p=8.7×10−5). No statistically significant result was observed for alcohol consumption.Conclusions:Our work demonstrates an important role of smoking, measured by both status and intensity, in the susceptibility to COVID-19.Funding:XJ is supported by research grants from the Swedish Research Council (VR-2018–02247) and Swedish Research Council for Health, Working Life and Welfare (FORTE-2020–00884).

Decision letter: Mendelian randomization analysis provides causality of smoking on the expression of ACE2, a putative SARS-CoV-2 receptor

ObjectivePrior research has revealed a heightened prevalence of neoplasms in individuals diagnosed with rheumatoid arthritis (RA). The primary objective of this study is to delve into the causal association between RA and two distinct types of neoplasms: benign neoplasm of bone and articular cartilage (BNBAC) and malignant neoplasm of bone and articular cartilage (MNBAC).MethodsWe employed summary data from genome-wide association analyses (GWAS) to investigate the causal relationship between RA and two neoplasms, BNBAC and MNBAC, using a two-sample bidirectional Mendelian randomization (MR) study design. The IEU OpenGWAS database provided the GWAS summary data for RA, while the Finnish consortium supplied the GWAS summary data for BNBAC and MNBAC. Our analysis involved the utilization of eight distinct MR methods, namely random-effects inverse variance weighted (IVW), MR Egger, weighted median, simple mode, weighted mode, maximum likelihood, penalized weighted median, and fixed effects IVW. Subsequently, we conducted assessments to evaluate heterogeneity, horizontal pleiotropy, outliers, the impact of a single-nucleotide polymorphism (SNP), and adherence to the assumption of normal distribution in the MR analysis.ResultsThe results from the MR analysis revealed that there was no significant genetic association between RA and BNBAC (P = 0.427, odds ratio [OR] 95% confidence interval [CI] = 0.971 [0.904–1.044]). However, a positive genetic association was observed between RA and MNBAC (P = 0.001, OR 95% CI = 1.413 [1.144–1.745]). Conducting a reverse MR analysis, we found no evidence to support a genetic causality between BNBAC (P = 0.088, OR 95% CI = 1.041 [0.994–1.091]) or MNBAC (P = 0.168, OR 95% CI = 1.013 [0.995–1.031]) and RA. Our MR analysis demonstrated the absence of heterogeneity, horizontal pleiotropy, and outliers and confirmed that the effect was not driven by a single SNP. Additionally, the data exhibited a normal distribution.ConclusionThe findings of this study demonstrate that RA constitutes a significant risk factor for MNBAC. In the context of clinical application, it is advisable to conduct MNBAC screening in RA patients and remain vigilant regarding its potential manifestation. Importantly, the outcomes of this investigation introduce a fresh vantage point into the understanding of the tumorigenesis associated with RA.

To determine whether a bidirectional causal relationship exists between major depressive disorder (MDD) and heart failure (HF). Our two-sample bidirectional Mendelian randomization (MR) study consisted of two parts. In the first part, we conducted a forward MR analysis where MDD was considered as the exposure and HF as the outcome. In the second part, a reverse MR analysis was performed, treating HF as the exposure and MDD as the outcome. Summary data on MDD and HF were obtained from the IEU Open GWAS database. Based on the results of the MR-Egger regression intercept test, there was no evidence of horizontal pleiotropy in this study. Furthermore, the IVW results consistently suggested estimates of causal effect values. The findings revealed that individuals with MDD had a 16.9% increased risk of HF compared to those without MDD (OR = 1.169, 95%CI: 1.044-1.308, P = 0.007). However, there was no evidence to support that HF would increase the risk of MDD (OR = 1.012, 95%CI: 0.932-1.099, P = 0.773). Heterogeneity in SNPs of MDD and HF was observed through the heterogeneity test and funnel plot. Additionally, the leave-one-out method did not identify any instances where a single SNP was biased toward or dependent on causation. Our study provides evidence supporting a one-way causal relationship between MDD and HF. Specifically, MDD increases the risk of developing HF. However, our findings did not provide any evidence suggesting that HF increases the risk of developing MDD.

ObjectivesTo assess the association between thyroid dysfunction and diabetic retinopathy (DR), a two-sample bidirectional Mendelian randomization (MR) study utilizing the Genome-wide Association Study (GWAS) database was conducted to investigate the causal relationship between these two variables.MethodsIn this study, GWAS of 48,328,151 single nucleotide polymorphisms(SNP) in the European population from the IEU open GWAS database were utilized as genetic tools for investigating thyroid dysfunction. The total sample size for the study on hyperthyroidism was 460,499 (case group: 3557; control group: 456,942). The total sample size for hypothyroidism was 410,141 (case group: 30,155; control group: 37,986). In addition, the data on DR were extracted from the FinnGen Biobank, comprising a total sample size of 319,046 individuals (10,413 cases and 308,633 controls). For the forward MR analysis, hyperthyroidism and hypothyroidism were considered as exposures with DR as the outcome. Reverse MR analysis was conducted using DR as exposure and hyperthyroidism and hypothyroidism as outcomes. Methods: The main analytical approach employed inverse variance weighting(IVW), supplemented by MR-Egger, Weighted mode method, weighted median, and Simple mode. Cochran's Q test, MR-PRESSO, MR-Egger and leave-one-out analysis were used to evaluate the sensitivity and pleiotropy.ResultsTwo-sample bidirectional MR analysis revealed a significant association between the presence of hyperthyroidism and hypothyroidism and an increased risk of DR in the forward MR analysis (IVW: OR = 1.29, 95% [CI] = 1.12–1.49, P < 0.001; OR = 1.17, 95% CI = 1.10–1.25, P < 0.001). In the reverse MR analysis, DR was found to be associated with an elevated risk of developing hyperthyroidism and hypothyroidism (IVW: OR = 1.56, 95% CI 1.38–1.76, P < 0.001; OR = 1.41, 95% CI 1.25–1.59, P < 0.001). Furthermore, most supplementary MR methods also demonstrated statistically significant differences and exhibited effect sizes consistent with those obtained from IVW. The sensitivity analysis confirmed the relative reliability of our causal findings.ConclusionsOur findings provide genetic evidence supporting a bidirectional causal relationship between thyroid function and DR.

The clinical incidence of sjogren's syndrome combined with gastroesophageal reflux disease is high. Existing observational studies have shown inconsistent results in the association between gastroesophageal reflux disease (GERD) and Sjogren's syndrome (SS).We observed that the symptoms of SS patients also improved after receiving GERD-related treatment. Therefore, we aimed to investigate the relationship between GERD and SS through a bidirectional two-sample Mendelian randomization (MR) study. Independent SNPs associated with GERD and SS were selected from a genome-wide association study (GWAS) as instrumental variables to conduct a bidirectional two-sample Mendelian analysis of GERD and SS. Genetic data were obtained from two databases for the following two outcomes: Gastroesophageal reflux (IEU Open GWAS) [sample size = 602,604 (patients = 129,080; nonpatients = 473,524)] and SS (FinnGen) [sample size = 392,423 (patients = 2,495; nonpatients = 389,928)]. Statistical methods for the MR analysis included the inverse-variance weighting method, weighted median, simple mode and weighted mode, as well as heterogeneity and sensitivity analyses using the Cochran Q statistic, MR‒Egger regression, outlier detection methods (MR-PRESSO). In addition, Steiger Test was conducted to test the direction of causality. MR analysis showed a positive correlation between GERD and SS risk [odds ratio (OR) = 1.3279 (95% confidence interval 1.0312–1.7099, P = 0.0280)]. However, in contrast, no significant causal effect of SS on GERD was observed [OR = 1.0024 (95% CI 0.9651–1.0412; P = 0.8995)]. This bidirectional two-sample Mendelian randomization study confirmed a causal relationship between SS and GERD, and suggested that GERD is a risk factor for SS, while SS does not affect GERD.

To investigate the genetic causality for the insomnia and common orthopedic diseases, such as rheumatoid arthritis (RA), ankylosing spondylitis (AS), osteoporosis (OP), and gout (GT). The genome-wide association study (GWAS) summary data on insomnia were obtained from a published study, while the GWAS summary data on RA, AS, OP, and GT were sourced from the FinnGen consortium. We utilized the TwoSampleMR package of the R software (version 4.1.2) to conduct a two-sample Mendelian randomization (MR) analysis. Our primary method of analysis was the random-effects inverse variance weighted (IVW) approach. Subsequently, we conducted a series of sensitivity analyses for the MR analysis. The MR analysis revealed a positive genetic causal relationship between insomnia and RA (P=0.016, odds ratio [OR] 95% confidence interval [CI]=1.112 [1.020-1.212]). However, no significant genetic causal relationship was observed between insomnia and AS (P=0.194, OR 95% CI=1.121 [0.944-1.331]), OP (P=0.788, OR 95% CI=1.016 [0.904-1.142]), and GT (P=0.757, OR 95% CI=1.018 [0.912-1.136]). The MR analysis did not exhibit heterogeneity, horizontal pleiotropy, outlier effects, or dependence on a single SNP, and demonstrated normal distribution, which guaranteed the robustness of the results. The results of this study suggest that insomnia may be a significant risk factor for RA, and controlling insomnia may represent a promising strategy for preventing RA. While insomnia was not observed to be associated with AS, OP, and GT at the genetic level, other levels of association cannot be excluded.

BackgroundSeveral observational studies have suggested a potential relationship between gut microbiome and psoriatic arthritis (PsA). However, the causality of this relationship still remains unclear. We aim to explore if the specific gut microbiome is causally associated with PsA at the genetic level and offer valuable insights into the etiology of PsA.MethodsIn this study, we employed a bidirectional two-sample Mendelian randomization (MR) analysis to investigate the causal effects of the gut microbiome on PsA. Publicly accessible genome-wide association study summary data of gut microbiome were obtained from the MiBioGen consortium (n = 14,306), while the summary statistics of psoriatic arthropathies were sourced from the FinnGen consortium R8 release data (2,776 cases and 221,323 controls). The primary analytical method employed was inverse variance weighted (IVW), complemented by supplementary methods including MR-Egger, weighted median, weighted mode, maximum likelihood, MR-PRESSO, and cML-MA. Reverse MR analysis was performed on the bacteria that were found to be causally associated with PsA in forward MR analysis. Cochran’s IVW Q statistic was utilized to assess the heterogeneity of instrumental variables among the selected single nucleotide polymorphisms.ResultsIVW estimates revealed that Ruminococcaceae_UCG-002 (odds ratio (OR) = 0.792, 95% confidence interval (CI), 0.643–0.977, p = 0.029) exhibited a protective effect on PsA. Conversely, Blautia (OR = 1.362, 95% CI, 1.008–1.842, p = 0.044), Eubacterium_fissicatena_group (OR = 1.28, 95% CI, 1.075–1.524, p = 0.006), and Methanobrevibacter (OR = 1.31, 95% CI, 1.059–1.621, p = 0.013) showed a positive correlation with the risk of PsA. No significant heterogeneity, horizontal pleiotropy, or outliers were observed, and the results of the MR analysis remained unaffected by any single nucleotide polymorphisms. According to the results of reverse MR analysis, no significant causal effect of PsA was found on gut microbiome.ConclusionThis study establishes for the first time a causal relationship between the gut microbiome and PsA, providing potential valuable strategies for the prevention and treatment of PsA. Further randomized controlled trials are urgently warranted to support the targeted protective mechanisms of probiotics on PsA.

Breast cancer is a common cancer type that leads to cancer-related deaths among women. HER2-positive breast cancer, in particular, is associated with poor prognosis due to its high aggressiveness, increased risk of recurrence, and metastasis potential. Previous observational studies have explored potential associations between inflammatory cytokines and the risk of two breast cancer subtypes (HER2-positive and HER2-negative), but the results have been inconsistent. To further elucidate the causal relationship between inflammatory cytokines and the two breast cancer subtypes, we conducted a two-sample Mendelian randomization (MR) study. We employed a two-sample bidirectional MR analysis using publicly available genome-wide association study (GWAS) statistics. After obtaining instrumental variables, we conducted MR analyses using five different methods to ensure the reliability of our results. Additionally, we performed tests for heterogeneity and horizontal pleiotropy. Subsequently, we conducted a reverse MR study by reversing exposure and outcome variables. Evidence from our IVW analysis revealed that genetically predicted levels of IL-5 [odds ratio (OR): 1.18, 95% confidence interval (CI): 1.04-1.35, P = 0.012], IL-7 (OR: 1.11, 95% CI: 1.01-1.22, P = 0.037), and IL-16 (OR: 1.13, 95% CI: 1.02-1.25, P = 0.025) were associated with an increased risk of HER2-positive breast cancer. Conversely, IL-10 (OR: 1.14, 95% CI: 1.03-1.26, P = 0.012) was associated with an increased risk of HER2-negative breast cancer. These results showed no evidence of heterogeneity or horizontal pleiotropy (P > 0.05). Results from the reverse MR analysis indicated no potential causal association between breast cancer and inflammatory cytokines (P > 0.05). Our findings demonstrate that IL-5, IL-7, and IL-16 are risk factors for HER2-positive breast cancer, with varying degrees of increased probability of HER2-positive breast cancer associated with elevated levels of these inflammatory cytokines. Conversely, IL-10 is a risk factor for HER2-negative breast cancer. Reverse studies have confirmed that breast cancer is not a risk factor for elevated levels of inflammatory cytokines. This series of results clarifies the causal relationship between different types of inflammatory cytokines and different subtypes of breast cancer. Based on this research, potential directions for the mechanism research of different inflammatory cytokines and different subtypes of breast cancer have been provided, and potential genetic basis for identifying and treating different subtypes of breast cancer have been suggested.

ObjectiveThe immune response assumes a pivotal role in the underlying mechanisms of urticaria pathogenesis. The present study delves into an investigation of the genetic causal connections between urticaria and prevalent autoimmune afflictions, notably rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), ulcerative colitis (UC), and Crohn’s disease (CD).MethodsA bidirectional two-sample Mendelian randomization (MR) analysis was conducted to investigate the causal relationships involving four autoimmune diseases and urticaria. The genome-wide association study (GWAS) summary data of four autoimmune disease were sourced from the IEU OpenGWAS database. The GWAS summary data for urticaria were derived from the Finnish consortium dataset. The principal analytical approach employed in this study was the random-effects inverse variance weighted (IVW) method. Subsequently, a series of sensitivity analyses were performed, encompassing assessments of heterogeneity, horizontal pleiotropy, outliers, “Leave-one-out” analyses, and tests for adherence to the assumption of normal distribution.ResultsThe random-effects IVW analysis indicate a positive genetic causal association between RA and urticaria (P < 0.001, OR 95% CI = 1.091 [1.051-1.133]). Conversely, SLE, UC, and CD do not exhibit a significant genetic causal relationship with urticaria. The reverse MR analysis reveals a positive genetic causal linkage between urticaria and SLE (P = 0.026, OR 95% CI = 1.289 [1.031-1.612]). However, the analysis demonstrates no substantial genetic causal relationship between urticaria and RA, UC, or CD. Importantly, the genetic causal assessment absence of heterogeneity, horizontal pleiotropy, and outliers. Furthermore, it remains unaffected by any individual single nucleotide polymorphism (SNP), demonstrating adherence to a normal distribution.ConclusionThis investigation establishing RA as a predisposing factor for urticaria. Moreover, urticaria as a plausible risk determinant for SLE. Heightened vigilance is recommended among RA patients to monitor the manifestation of urticaria within clinical settings. Similarly, individuals afflicted by urticaria should duly acknowledge the prospective susceptibility to SLE.

The causal relationship between adiponectin (ADPN) and estimated glomerular filtration rate (eGFR) is unclear. This study adopts a two-sample bidirectional Mendelian randomization (MR) study to explore the causal relationship between ADPN and eGFR. Using eight single nucleotide polymorphisms (SNP) of ADPN and 26 SNP of eGFR as instrumental variables, the study performs a two-sample bidirectional MR study using MR inverse-variance weighted (IVW), MR-Egger and weighted median approach to evaluate the causal relationship between ADPN and eGFR. Using the genetic risk score (GRS) of ADPN and eGFR as instrumental variables, the study performs a second MR analysis to assess the association between ADPN and eGFR. In ADPN to eGFR MR analysis, the IVW, weighted median and GRS analysis all showed that ADPN had a causal effect on eGFR after removing potential confounders of the ADPN-eGFR relation (IVW: β = .016, P = .002; weighted median: β = .012, P = .022; GRS: β = .016, P = 1.48E-05). As both ADPN and eGFR were natural log-transformed in the corresponding GWAS, eGFR increased by 0.15% for any 10% increase in ADPN. In eGFR to ADPN MR analysis, eGFR had no causal effect on ADPN after removing potential confounders of the eGFR-ADPN relation (All P values > 0.05). The heterogeneity test and sensitivity analysis indicated some heterogeneity, but no directional pleiotropy. Adiponectin has a causal effect on eGFR, while eGFR has no causal effect on ADPN. ADPN may be a clinical target for improving eGFR and treating chronic kidney disease caused by decreased eGFR.

Background Our objective is to investigate the potential causal relationship between telomere length (TL) and aneurysmal subarachnoid hemorrhage (aSAH) and intracranial aneurysms (IAs) by conducting a bidirectional two-sample Mendelian Randomization (MR) study. Methods We utilized publicly available summary data from genome-wide association studies (GWAS) for comprehensive analysis. Telomere length-associated data were sourced from the Epidemiology Unit (IEU) GWAS database (n = 472,174), while data pertaining to intracranial aneurysms were derived from a GWAS meta-analysis conducted by Bakker et al, encompassing aneurysmal subtypes including aSAH (n = 77,074), IAs (n = 79,429), and unruptured intracranial aneurysms (uIA) (n = 74,004), all sampled from European populations. The primary method for MR analysis employed was the Inverse Variance Weighted (IVW) method. Additionally, we conducted various sensitivity analyses to assess the heterogeneity and pleiotropy of study findings. Reverse MR analysis was employed to explore potential reverse causality. Results In the forward MR analysis, the IVW method indicated a negative association between TL and aSAH (OR = 0.636, 95% CI: 0.459–0.883, p = 0.006) as well as IAs (OR = 0.670, 95% CI: 0.499–0.900, p = 0.0079). There was no evidence of heterogeneity or horizontal pleiotropy in the forward MR analysis. Reverse MR analysis did not reveal any causal relationship between aSAH, IAs, uIA and TL. Conclusions In European populations, there exists a causal relationship between longer TL and reduced risks of aSAH and IAs Further research is warranted to elucidate the underlying mechanisms and the potential of TL as an intervention target for lowering the incidence of aSAH and IAs.

BackgroundNicotine dependence is a key factor influencing the diversity of gut microbiota, and targeting gut microbiota may become a new approach for the prevention and treatment of nicotine dependence. However, the causal relationship between the two is still unclear. This study aims to investigate the causal relationship between nicotine dependence and gut microbiota.MethodsA two-sample bidirectional Mendelian randomization (MR) study was conducted using the largest existing gut microbiota and nicotine dependence genome-wide association studies (GWAS). Causal relationships between genetically predicted nicotine dependence and gut microbiota abundance were examined using inverse variance weighted, MR-Egger, weighted median, simple mode, weighted mode, and MR-PRESSO approaches. Cochrane’s Q test, MR-Egger intercept test, and leave-one-out analysis were performed as sensitivity analyses to assess the robustness of the results. Multivariable Mendelian randomization analysis was also conducted to eliminate the interference of smoking-related phenotypes. Reverse Mendelian randomization analysis was then performed to determine the causal relationship between genetically predicted gut microbiota abundance and nicotine dependence.ResultsGenetically predicted nicotine dependence had a causal effect on Christensenellaceae (β: -0.52, 95% CI: -0.934–0.106, P = 0.014). The Eubacterium xylanophilum group (OR: 1.106, 95% CI: 1.004-1.218), Lachnoclostridium (OR: 1.118, 95% CI: 1.001-1.249) and Holdemania (OR: 1.08, 95% CI: 1.001-1.167) were risk factors for nicotine dependence. Peptostreptococcaceae (OR: 0.905, 95% CI: 0.837-0.977), Desulfovibrio (OR: 0.014, 95% CI: 0.819-0.977), Dorea (OR: 0.841, 95% CI. 0.731-0.968), Faecalibacterium (OR: 0.831, 95% CI: 0.735-0.939) and Sutterella (OR: 0.838, 95% CI: 0.739-0.951) were protective factor for nicotine dependence. The sensitivity analysis showed consistent results.ConclusionThe Mendelian randomization study confirmed the causal link between genetically predicted risk of nicotine dependence and genetically predicted abundance of gut microbiota. Gut microbiota may serve as a biomarker and offer insights for addressing nicotine dependence.

Previous observational studies reported that ankylosing spondylitis is closely related to hypertension. However, it is still controversial whether the association between ankylosing spondylitis and hypertension is causal. The effects of ankylosing spondylitis on diastolic and systolic blood pressure deserve further investigation. The objective of our study is to explore whether ankylosing spondylitis is causally associated with blood pressure. A bidirectional two-sample Mendelian randomization (MR) analysis was performed by employing five Mendelian randomization analysis methods. MR Egger regression, weighted median, inverse variance weighted, and weight mode methods were performed in the two-sample Mendelian randomization analysis. We performed Mendelian randomization to investigate the association between ankylosing spondylitis (finn-b-M13_ANKYLOSPON) and hypertension (ukb-b-14057), diastolic blood pressure (ebi-a-GCST90000062) and systolic blood pressure (ebi-a-GCST90000059). We also performed reverse Mendelian randomization between exposures and outcomes. Another new validation cohort (ukb-b-18194) was also performed. The heterogeneity, horizontal pleiotropy, and possible outliers were examined in the MR analysis results. The inverse variance weighted results showed that ankylosing spondylitis has no genetic causal relationship with hypertension (p = 0.441, OR = 1.001, 95% CI: 0.999-1.002). The inverse variance weighted results showed that ankylosing spondylitis has no genetic causal relationship with systolic blood pressure (p = 0.301, OR = 1.006, 95% CI: 0.995-1.018). The inverse variance weighted results showed that ankylosing spondylitis has no genetic causal relationship with diastolic blood pressure (p = 0.778, OR = 1.002, 95% CI: 0.988-1.016). The reverse Mendelian randomization between exposures and outcomes is negative. Another new validation cohort also confirmed the results. No heterogeneity was observed by the MR-pleiotropy residual sum and outlier test. The "leave-one-out" analysis indicated that the results of MR analysis were not affected by a single nucleotide polymorphism. This study represents the first two-sample Mendelian randomization analysis aimed at investigating the causal genetic relationship between ankylosing spondylitis and blood pressure. Our Mendelian randomization analysis results revealed a lack of causal association between ankylosing spondylitis and hypertension, diastolic blood pressure, as well as systolic blood pressure.

Psoriasis is a chronic, refractory inflammatory skin disease, with a high prevalence of psychiatric comorbidities, including depression, anxiety, and even suicidality, which may in turn initiate or exacerbate skin inflammation. However, the causal relationships between these comorbidities remain unclear. To investigate the cause-effect relationships between psoriasis and mental disorders including depression, anxiety, and suicidality, we conducted a bidirectional two-sample Mendelian randomization (MR) study utilizing summary statistics from the most comprehensive genome-wide association studies of psoriasis (n = 306 123), broad depression (n = 500 199), major depressive disorder (n = 173 005), anxiety (n = 17 310), and suicide attempts (n = 50 264). Using the random-effects inverse-variance weighted method as primary method, the forward MR analyses indicated that psoriasis was significantly associated with higher odds of broad depression (odds ratio [OR] 1.030, 95% confidence interval [CI] 1.010-1.051, P = 0.003) and suggestively associated with an increased risk of major depressive disorder (OR 1.054, 95% CI 1.002-1.109, P = 0.040), but not with the risk of anxiety (P = 0.160) or suicide attempts (P = 0.648). In reverse MR analyses, significant causal impact of broad depression (OR 1.363, 95% CI 1.103-1.684, P = 0.004) and major depressive disorder (OR 1.890, 95% CI 1.285-2.781, P = 0.001), but not anxiety (P = 0.787) and suicide attempts (P = 0.961) on psoriasis risk was observed. In addition, the results of primary analysis are consistent across sensitivity analyses, albeit the MR-Egger regression model produced wide CIs and negative results in several analyses. In conclusion, this MR study indicates a bidirectional causal relationship between psoriasis and depression that was previously unrecognized, which highlights the significance of screening for depression in psoriasis patients and initiating appropriate interventions. Further studies are required to elucidate the pathophysiology of the bidirectional causal relationship between these two conditions.