Combining the strengths of inverse-variance weighting and Egger regression in Mendelian randomization using a mixture of regressions model

Mendelian-randomization study reveals causal relationships between nitrogen dioxide and gut microbiota

Exposure to nitrogen dioxide might potentially change the makeup and operation of gut microbes. Nitrogen dioxide data was procured from the IEU Open GWAS (N = 456 380). Subsequently, a two-sample Mendelian randomization study was executed, utilizing summary statistics of gut microbiota sourced from the most expansive available genome-wide association study meta-analysis, conducted by the MiBioGen consortium (N = 13 266). The causal relationship between nitrogen dioxide and gut microbiota was determined using inverse variance weighted, maximum likelihood, MR-Egger, Weighted Median, Weighted Model, Mendelian randomization pleiotropy residual sum and outlier, and constrained maximum likelihood and model averaging and Bayesian information criterion. The level of heterogeneity of instrumental variables was quantified by utilizing Cochran's Q statistic. The colocalization analysis was used to examine whether nitrogen dioxide and the identified gut microbiota shared casual variants. Inverse variance weighted estimate suggested that nitrogen dioxide was causally associated with Akkermansia (β = -1.088, 95% CI: -1.909 to -0.267, P = 0.009). In addition, nitrogen dioxide presented a potential association with Bacteroides (β = -0.938, 95% CI: -1.592 to -0.284, P = 0.005), Barnesiella (β = -0.797, 95% CI: -1.538 to -0.055, P = 0.035), Coprococcus 3 (β = 1.108, 95% CI: 0.048-2.167, P = 0.040), Eubacterium hallii group (E. hallii) (β = 0.776, 95% CI: 0.090-1.463, P = 0.027), Holdemania (β = -1.354, 95% CI: -2.336 to -0.372, P = 0.007), Howardella (β = 1.698, 95% CI: 0.257-3.139, P = 0.021), Olsenella (β = 1.599, 95% CI: 0.151-3.048, P = 0.030) and Sellimonas (β = -1.647, 95% CI: -3.209 to -0.086, P = 0.039). No significant heterogeneity of instrumental variables or horizontal pleiotropy was found. The associations of nitrogen dioxide with Akkermansia (PH4 = 0.836) and E. hallii (PH4 = 0.816) were supported by colocalization analysis. This two-sample Mendelian randomization study found that increased exposure to nitrogen dioxide had the potential to impact the human gut microbiota.

Metabolic disorders and risk of cardiovascular diseases: a two-sample mendelian randomization study

BACKGROUND

Metabolic disorders are increasing worldwide and are characterized by various risk factors such as abdominal obesity, insulin resistance, impaired glucose metabolism, and dyslipidemia. Observational studies suggested a bidirectional association between cardiovascular diseases and metabolic disorders and its components. However, the causal associations between them remained unclear. This study aims to investigate the causal relationship between metabolic disorders and cardiovascular disease through Mendelian randomization (MR) analysis.

METHODS

A two-sample MR analysis based on publicly available genome-wide association studies were used to infer the causality. The single-nucleotide polymorphisms with potential pleiotropy were excluded by MR-PRESSO. The effect estimates were constructed using the random-effects inverse-variance-weighted method as the primary estimate. Furthermore, MR-Egger and weighted median were also performed to detect heterogeneity and pleiotropy.

RESULTS

Genetically predicted metabolic disorders increased the risk for coronary heart disease (OR = 1.77, 95% CI: 1.55-2.03, p < 0.001), myocardial infarction (OR = 1.75, 95% CI: 1.52-2.03, p < 0.001), heart failure (OR = 1.26, 95% CI: 1.14-1.39, p < 0.001), hypertension (OR = 1.01, 95% CI: 1.00-1.02, p = 0.002), and stroke (OR = 1.19, 95% CI: 1.08-1.32, p < 0.001). The concordance of the results of various complementary sensitivity MR methods reinforces the causal relationship further.

CONCLUSION

This study provides evidence of a causal relationship between metabolic disorders and increased risk of coronary heart disease, myocardial infarction, heart failure, hypertension, and stroke. Special attention should be paid to improving metabolic disorders to reduce the development of cardiovascular diseases.

Association between gut microbiota and sensorineural hearing loss: a Mendelian randomization study

Background

Several recent studies speculated that the gut microbiota is associated with sensorineural hearing loss (SNHL) and proposed the concept of the gut-inner ear axis. However, the causal effect of gut microbiota on SNHL is still unknown. In this study, we performed a two-sample Mendelian randomization (MR) analysis to estimate the causal effect of gut microbiota on SNHL.

Methods

Gut microbiota data were obtained from the largest available genome-wide association study (n = 18,340) conducted by the MiBioGen consortium. The summary statistics of SNHL were obtained from the FinnGen consortium R8 release data (28,310 cases and 302,750 controls). The causal effects were estimated with inverse-variance weighted, MR-Egger, and weighted median. Reverse Mendelian randomization analysis was performed on the bacteria that were found to be associated with SNHL in forward Mendelian randomization analysis. We then performed sensitivity analyses, including Cochran's Q-test, MR-Egger intercept test, MR-PRESSO, cML-MA-BIC, and leave-one-out analysis, to detect heterogeneity and pleiotropy.

Results

The inverse-variance weighted results suggested that Lachnospiraceae (UCG001) had a significant protective effect against SNHL (odds ratio = 0.85, 95% confidence interval: 0.78-0.93, P = 6.99 × 10-4). In addition, Intestinimonas (odds ratio = 0.89, 95% confidence interval: 0.82-0.97, P = 8.53 × 10-3) presented a suggestively protective effect on SNHL. Rikenellaceae (RC9gutgroup) (odds ratio = 1.08, 95% confidence interval: 1.02-1.15, P = 0.01) and Eubacterium (hallii group) (odds ratio = 1.12, 95% confidence interval: 1.00-1.24, P = 0.048) suggestively increase the risk of SNHL. The results of the reverse MR analysis showed that there is no significant causal effect of SNHL on the gut microbiota. No significant heterogeneity of instrumental variables or pleiotropy was detected.

Conclusion

The evidence that the four genera mentioned above are associated with SNHL supports the hypothesis of a gut-inner ear axis. Our study provides microbial markers for the prevention and treatment of SNHL, and further studies are needed to explore the mechanisms of the gut microbiome-inner ear axis in health and diseases.

SGLT2 inhibition, circulating metabolites, and atrial fibrillation: a Mendelian randomization study

BACKGROUND

Sodium-glucose cotransporter 2 (SGLT2) inhibitors have shown promise in reducing the risk of atrial fibrillation (AF). However, the results are controversial and the underlying metabolic mechanism remains unclear. Emerging evidence implied that SGLT2 inhibitors have extra beneficial metabolic effects on circulating metabolites beyond glucose control, which might play a role in reducing the risk of AF. Hence, our study aimed to investigate the effect of circulating metabolites mediating SGLT2 inhibition in AF by Mendelian randomization (MR).

METHODS

A two-sample and two-step MR study was conducted to evaluate the association of SGLT2 inhibition with AF and the mediation effects of circulating metabolites linking SGLT2 inhibition with AF. Genetic instruments for SGLT2 inhibition were identified as genetic variants, which were both associated with the expression of SLC5A2 gene and glycated hemoglobin level (HbA1c). Positive control analysis on type 2 diabetes mellitus (T2DM) was conducted to validate the selection of genetic instruments.

RESULTS

Genetically predicted SGLT2 inhibition (per 1 SD decrement in HbA1c) was associated with reduced risk of T2DM (odds ratio [OR] = 0.63 [95% CI 0.45, 0.88], P = 0.006) and AF (0.51 [0.27, 0.97], P = 0.039). Among 168 circulating metabolites, two metabolites were both associated with SGLT2 inhibition and AF. The effect of SGLT2 inhibition on AF through the total concentration of lipoprotein particles (0.88 [0.81, 0.96], P = 0.004) and the concentration of HDL particles (0.89 [0.82, 0.97], P = 0.005), with a mediated proportion of 8.03% (95% CI [1.20%, 14.34%], P = 0.010) and 7.59% ([1.09%, 13.34%], P = 0.011) of the total effect, respectively.

CONCLUSIONS

This study supported the association of SGLT2 inhibition with a reduced risk of AF. The total concentration of lipoprotein particles and particularly the concentration of HDL particles might mediate this association. Further mechanistic and clinical studies research are needed to understand the mediation effects of circulating metabolites especially blood lipids in the association between SGLT2 inhibition and AF.

Type 1 diabetes, glycemic traits, and risk of dental caries: a Mendelian randomization study

Background: Regarding past epidemiological studies, there has been disagreement over whether type 1 diabetes (T1DM) is one of the risk factors for dental caries. The purpose of this study was to determine the causative links between genetic susceptibility to T1DM, glycemic traits, and the risk of dental caries using Mendelian randomization (MR) approaches. Methods: Summary-level data were collected on genome-wide association studies (GWAS) of T1DM, fasting glucose (FG), glycated hemoglobin (HbA1c), fasting insulin (FI), and dental caries. MR was performed using the inverse-variance weighting (IVW) method, and sensitivity analyses were conducted using the MR-Egger method, weighted median, weighted mode, replication cohort, and multivariable MR conditioning on potential mediators. Results: The risk of dental caries increased as a result of genetic susceptibility to T1DM [odds ratio (OR) = 1.044; 95% confidence interval (CI) = 1.015-1.074; p = 0.003], with consistent findings in the replication cohort. The relationship between T1DM and dental caries was stable when adjusted for BMI, smoking, alcohol intake, and type 2 diabetes (T2DM) in multivariable MR. However, no significant correlations between the risk of dental caries and FG, HbA1c, or FI were found. Conclusion: These results indicate that T1DM has causal involvement in the genesis of dental caries. Therefore, periodic reinforcement of oral hygiene instructions must be added to the management and early multidisciplinary intervention of T1DM patients, especially among adolescents and teenagers, who are more susceptible to T1DM.

Integrated Analyses of Single-Cell Transcriptome and Mendelian Randomization Reveal the Protective Role of Resistin in Sepsis Survival in Intensive Care Unit

The high morbidity and mortality rates associated with sepsis highlight the challenges of finding specific remedies for this condition in the intensive care unit (ICU). This study aimed to explore the differentially expressed genes (DEGs) specific to cell types in sepsis and investigate the role of resistin in the survival of sepsis patients through Mendelian randomization (MR) analyses. We used single-cell and bulk transcriptome data to identify cell type-specific DEGs between sepsis and healthy controls. MR analyses were then conducted to investigate the causal relationships between resistin (one of the identified DEGs) levels and the survival of sepsis patients. Additionally, we utilized meQTL (methylation quantitative trait loci) to identify cytosine-phosphate-guanine (CpG) sites that may directly affect sepsis. We identified 560 cell type-specific DEGs between sepsis and healthy controls. Notably, we observed the upregulation of resistin levels in macrophages during sepsis. In bulk transcriptome, RETN is also upregulated in sepsis samples compared with healthy controls. MR analyses revealed a negative association existed between the expression of resistin, at both gene and protein levels, and the mortality or severity of sepsis patients in ICU. Moreover, there were no associations observed between resistin levels and death or organ failure due to other causes. We also identified three methylation CpG sites, located in RETN or its promoter region-cg06633066, cg22322184, and cg02346997-that directly affected both resistin protein levels and sepsis death in the ICU. Our findings suggest that resistin may provide feasible protection for sepsis patients, particularly those with severe cases, without serious side effects. Therefore, resistin could be a potential drug candidate for sepsis treatment. Additionally, we identified two CpG sites, cg06633066 and cg22322184, that were associated with RETN protein levels and sepsis death, providing novel insights into the underlying mechanisms of sepsis.

Investigating the Impact of Sun/UV Protection and Ease of Skin Tanning on the Risk of Pseudoexfoliation Glaucoma: A Mendelian Randomization Study

Purpose

To investigate the potential causal associations between the use of sun/ultraviolet (UV) protection and ease of skin tanning and the risk of pseudoexfoliation glaucoma (PXG) in European populations.

Methods

Single nucleotide polymorphisms (SNPs) associated with the use of sun/UV protection and ease of skin tanning were selected from the UK Biobank genome-wide association study database consisting of 498,751 European participants. SNPs of PXG were obtained from the FinnGen study including 3424 PXG cases and 326,434 controls. Two-sample Mendelian randomization (MR) analyses were performed to assess the association between the use of sun/UV protection and ease of skin tanning and risk of PXG.

Results

Inverse variance weighted regression of genetic susceptibility predicted that both use of sun/UV protection and ease of skin tanning were potentially positively associated with the decreased risk of PXG in the European ancestry (use of sun/UV protection: odds ratio [OR] = 0.47; 95% confidence interval [CI], 0.24-0.92; P = 0.028; ease of skin tanning: OR = 0.81; 95% CI, 0.67-0.97; P = 0.025).

Conclusions

We found genetic evidence supporting a potential causal association between UV protection and a decreased risk of PXG in European population. Further research will help elucidate the underlying mechanisms and promote UV protection for eyes, especially in people with a high risk of PXG.

Causality between depression and ankylosing spondylitis in a European population: Results from a Mendelian randomization analysis

The aim of this study was to explore the application of Mendelian randomization (MR) Egger and inverse variance weighted (IVW) in a causal effect on depression and ankylosing spondylitis (AS). Instrumental variables (IVs) were determined using genome-wide association studies. The 2-sample MR analysis was conducted by MR Egger to test the causal effect between depression and AS. The pleiotropy of potential instrumental variables was evaluated. The results of MR Egger and IVW were further compared. A total of 3 single nucleotide polymorphisms as the construct IVs were included. IVW results showed a significant causal effect between depression and AS (P < .001). Depression could promote the risk of AS (odds ratio = 1.060, 95% confidence interval: 1.026-1.094). However, the MR Egger showed no causal effect (P = .311). Heterogeneity statistics suggested that no heterogeneity was existed (P > .05). It was also suggested that there was no horizontal pleiotropy in IVs (MR Egger intercept: -0.0004, P = .471). Reverse MR analysis suggested that there was no causal effect between AS and depression (P > .05). Gene expression quantitative trait locus (QTLs) suggested that rs2517601 and RNF39 were positively correlated (beta = 1.066, P < .001). Depression may be one of the causes of AS by MR analysis in a European population. We can estimate the causal effect based on IVW when horizontal pleiotropy is very tiny.

Causal relationship between leisure sedentary behaviors and low back pain risk: a Mendelian randomization study

PURPOSE

This study aims to assess the causal associations of leisure sedentary behaviors with low back pain (LBP).

METHODS

A Mendelian randomization (MR) study was carried out utilizing genetic instruments to determine whether leisure sedentary behaviors (including leisure television watching, leisure computer use, and driving) are causally associated with LBP. All instrumental variables were selected from publicly available genetic summary data. The inverse variance weighted (IVW) was used as the main method to conduct univariable MR analyses. Further sensitivity analyses were utilized to test the stability of the results. Moreover, multivariable MR was performed to evaluate the independent causal relationship between leisure sedentary behaviors and LBP when body mass index (BMI), waist circumference, smoking initiation, and vigorous physical activity were taken into account.

RESULTS

The MR analyses showed evidence that television watching increased the risk for LBP (OR: 1.97, 95% CI 1.45, 2.66; P = 1.19 × 10-5). Genetically determined computer use is causally associated with a decreased risk of LBP (OR: 0.53, 95% CI 0.41, 0.68; P = 4.79 × 10-7). However, no evidence was found of a causal relationship between driving and LBP (OR: 2.27, 95% CI 0.75, 6.81; P = 0.145). After adjusting for BMI, waist circumference, smoking initiation, and vigorous physical activity, only television maintained its causal effect on LBP.

CONCLUSIONS

This study indicated that genetically predicted television watching was a risk factor for LBP independent of BMI, waist circumference, smoking initiation, and vigorous physical activity. This finding may be helpful for the diagnosis and management of LBP.

Identifying the causal relationship between sedentary behavior and heart failure: Insights from a Mendelian randomization study and mediation analysis

BACKGROUND

Observational studies have revealed that a lack of physical exercise may be linked to a higher risk of heart failure (HF). Here, the causal relationship between sedentary behavior (SB) and HF was investigated using Mendelian randomization (MR).

HYPOTHESIS

SB was considered as an important risk factor of HF.

METHODS

Single nucleotide polymorphisms with a genome-wide statistical significance threshold of <5 × 10-8 among the SB-proxied phenotypes (TV screen time, computer use, and driving) from genome-wide association study (GWAS) datasets were identified as instrumental variables (IVs). The MR study was performed using the inverse-variance weighting (IVW) model as a primary standard to evaluate causal relationships. Simultaneously, MR-Egger regression, weighted median, and maximum likelihood models were used as supplements. Sensitivity analysis, consisting of a heterogeneity and horizontal pleiotropy test, was performed using Cochran's Q, MR-Egger intercept, and MR-PRESSO tests to ensure the reliability of conclusions.

RESULTS

The IVW model results showed that increased TV screen time correlated with a higher genetic susceptibility for HF in both HF-associated GWAS datasets, which was also supported by weighted median and maximum likelihood model results. The odds ratios with 95% confidence intervals were 1.418 (1.182-1.700) and 1.486 (1.136-1.943), respectively. Although the results of Cochran's Q test indicated certain heterogeneity among the IVs. The MR-Egger intercept and MR-PRESSO tests suggested no horizontal pleiotropy and verified the reliability of the conclusion.

CONCLUSIONS

This MR study identified that increased TV screen time may predispose individuals to the development of HF.

COVID-19 hospitalization increases the risk of developing glioblastoma: a bidirectional Mendelian-randomization study

Background

As a result of the COVID-19 pandemic, patients with glioblastoma (GBM) are considered a highly vulnerable population. Despite this, the extent of the causative relationship between GBM and COVID-19 infection is uncertain.

Methods

Genetic instruments for SARS-CoV-2 infection (38,984 cases and 1,644,784 control individuals), COVID-19 hospitalization (8,316 cases and 1,549,095 control individuals), and COVID-19 severity (4,792 cases and 1,054,664 control individuals) were obtained from a genome-wide association study (GWAS) from European populations. A total of 6,183 GBM cases and 18,169 controls from GWAS were enrolled in our study. Their associations were evaluated by applying Mendelian randomization (MR) including IVW meta-analysis, MR-Egger regression, and weighted-median analysis. To make the conclusions more robust and reliable, sensitivity analyses were performed.

Results

Our results showed that genetically predicted COVID-19 hospitalization increases the risk of GBM (OR = 1.202, 95% CI = 1.035-1.395, p = 0.016). In addition, no increased risk of SARS-CoV-2 infection, COVID-19 hospitalization and severity were observed in patients with any type of genetically predicted GBM.

Conclusion

Our MR study indicated for the first time that genetically predicted COVID-19 hospitalization was demonstrated as a risk factor for the development of GBM.

Genetically predicted circulating levels of cytokines and the risk of depression: a bidirectional Mendelian-randomization study

Objective: Inflammatory cytokines disturbance is the main result of immune dysregulation, which is widely described in major depressive disorder (MDD). However, the potential causal relationship between these two factors has not been discovered. Therefore, the purpose of this study was to investigate the causal relationship between inflammatory cytokines and MDD risk by using the two-sample Mendelian randomization (MR) analysis. Method: Two genetic instruments obtained from publicly available gene profile data were utilized for the analysis. We obtained the genetic variation data of 41 inflammatory cytokines from genome-wide association studies (GWAS) meta-analysis of 8293 individuals of Finnish descent. The MDD data, including 135,458 MDD cases and 344,901 controls, were obtained from the Psychiatric Genomics Consortium Database. For the Mendelian randomization (MR) estimation, several methods were employed, namely, MR-Egger regression, inverse-variance weighted (IVW), weighted median, and MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) methods. Result: A causal relationship was identified between the genetically proxied levels of Interleukin (IL) -18, IL-1β, and Regulated upon activation normal T cell expressed and secreted (RANTES) and the risk of MDD (OR = 0.968, 95%CI = 0.938, 0.998, p = 0.036; OR = 0.875, 95%CI = 0.787, 0.971, p = 0.012; OR = 0.947, 95%CI = 0.902, 0.995, p = 0.03; respectively). However, our Mendelian randomization (MR) estimates provided no causality of MDD on inflammatory cytokines. Conclusion: Our study elucidates the connection between inflammatory cytokines and MDD by using MR analysis, thereby enhancing our comprehension of the potential mechanisms. By identifying these associations, our findings hold substantial implications for the development of more effective treatments aimed at improving patient outcomes. However, further investigation is required to fully comprehend the exact biological mechanisms involved.

Mendelian randomization indicates that atopic dermatitis contributes to the occurrence of diabetes

BACKGROUND

An association has been indicated between atopic dermatitis (AD), a prevalent chronic inflammatory skin disease, and diabetes mellitus. However, the exact causal relationship between AD and both type 1 diabetes (T1D) and type 2 diabetes (T2D) remains controversial. This study aimed to explore the causal association between AD and diabetes by Mendelian Randomization (MR) approaches.

METHODS

Public genetic summary data for AD was obtained from EAGLE study. Single nucleotide polymorphisms of diabetes were retrieved from four genome-wide association studies that had been performed in European populations. Inverse variance weighted (IVW) in MR analysis was used as the primary means of causality estimation. Several complementary analyses and sensitivity analyses were performed to calculate MR estimates and to enhance the causal inference, respectively. The R package 'TwoSampleMR' was used for analysis.

RESULTS

Genetically predicted AD led to a higher risk of T1D (OR, 1.19; 95% CI, 1.05, 1.34; P = 0.006) and T2D (OR, 1.07; 95% CI, 1.02, 1.11; P = 0.003) based on random-effect IVW method. The complementary analyses provided similar positive results. Cochran's Q test and I2 statistics indicated moderate heterogeneity between AD and both T1D and T2D. No significant horizontal pleiotropy was detected by MR-Egger Intercept p except summary data from FinnGen consortium.

CONCLUSION

Genetically predicted AD is a risk factor for both T1D and T2D. These findings imply potential shared pathological mechanisms between AD and diabetes, thus suggesting the significance of early clinical diagnosis and prevention of AD in reducing the incidence of diabetes.

Association between periodontitis and breast cancer: two-sample Mendelian randomization study

OBJECTIVES

The purpose of this study was to investigate whether there is a causal relationship between periodontitis and breast cancer by Mendelian randomization analysis.

MATERIALS AND METHODS

We performed a two-sample bidirectional Mendelian randomization (MR) analysis using publicly released genome-wide association studies (GWAS) statistics. The inverse-variance weighted (IVW) method was used as the primary analysis. We applied complementary methods, including weighted median, weighted mode, simple mode, MR-Egger regression, and MR-pleiotropy residual sum and outlier (MR-PRESSO) to detect and correct for the effect of horizontal pleiotropy.

RESULTS

IVW MR analysis showed no effect of periodontitis on breast cancer (IVW OR=0.99, P =0.14). Similarly, no significant causal relationship between breast cancer and periodontitis was found in reverse MR analysis (IVW OR=0.95, P =0.83). The results of MR-Egger regression, weighted median, and weighted mode methods were consistent with those of the IVW method. Based on sensitivity analyses, horizontal pleiotropy is unlikely to distort causal estimates.

CONCLUSIONS

Although observational studies have reported an association between periodontitis and breast cancer, the results of our MR analysis do not support a causal relationship between periodontitis and breast cancer.

CLINICAL RELEVANCE

Mendelian randomization study can more clearly analyze the causal relationship between periodontitis and breast cancer, in order to provide a certain reference for clinicians and deepen the understanding of the relationship between periodontitis and breast cancer, to explore more possible associations between periodontitis and systemic diseases.

Combining Mendelian randomization and network deconvolution for inference of causal networks with GWAS summary data

Mendelian randomization (MR) has been increasingly applied for causal inference with observational data by using genetic variants as instrumental variables (IVs). However, the current practice of MR has been largely restricted to investigating the total causal effect between two traits, while it would be useful to infer the direct causal effect between any two of many traits (by accounting for indirect or mediating effects through other traits). For this purpose we propose a two-step approach: we first apply an extended MR method to infer (i.e. both estimate and test) a causal network of total effects among multiple traits, then we modify a graph deconvolution algorithm to infer the corresponding network of direct effects. Simulation studies showed much better performance of our proposed method than existing ones. We applied the method to 17 large-scale GWAS summary datasets (with median N = 256879 and median #IVs = 48) to infer the causal networks of both total and direct effects among 11 common cardiometabolic risk factors, 4 cardiometabolic diseases (coronary artery disease, stroke, type 2 diabetes, atrial fibrillation), Alzheimer's disease and asthma, identifying some interesting causal pathways. We also provide an R Shiny app (https://zhaotongl.shinyapps.io/cMLgraph/) for users to explore any subset of the 17 traits of interest.

Effects of epigenetic age acceleration on kidney function: a Mendelian randomization study

BACKGROUND

Previous studies have reported cross-sectional associations between measures of epigenetic age acceleration (EAA) and kidney function phenotypes. However, the temporal and potentially causal relationships between these variables remain unclear. We conducted a bidirectional two-sample Mendelian randomization study of EAA and kidney function. Genetic instruments for EAA and estimate glomerular filtration rate (eGFR) were identified from previous genome-wide association study (GWAS) meta-analyses of European-ancestry participants. Causal effects of EAA on kidney function and kidney function on EAA were assessed through summary-based Mendelian randomization utilizing data from the CKDGen GWAS meta-analysis of log-transformed estimated glomerular filtration rate (log-eGFR; n = 5,67,460) and GWAS meta-analyses of EAA (n = 34,710). An allele score-based Mendelian randomization leveraging individual-level data from UK Biobank participants (n = 4,33,462) further examined the effects of EAA on kidney function.

RESULTS

Using summary-based Mendelian randomization, we found that each 5 year increase in intrinsic EAA (IEAA) and GrimAge acceleration (GrimAA) was associated with - 0.01 and - 0.02 unit decreases in log-eGFR, respectively (P = 0.02 and P = 0.09, respectively), findings which were strongly supported by allele-based Mendelian randomization study (both P < 0.001). Summary-based Mendelian randomization identified 24% increased odds of CKD with each 5-unit increase in IEAA (P = 0.05), with consistent findings observed in allele score-based analysis (P = 0.07). Reverse-direction Mendelian randomization identified potentially causal effects of decreased kidney function on HannumAge acceleration (HannumAA), GrimAA, and PhenoAge acceleration (PhenoAA), conferring 3.14, 1.99, and 2.88 year decreases in HanumAA, GrimAA, and PhenoAA, respectively (P = 0.003, 0.05, and 0.002, respectively) with each 1-unit increase in log-eGFR.

CONCLUSION

This study supports bidirectional causal relationships between EAA and kidney function, pointing to potential prevention and therapeutic strategies.

Robust multivariable Mendelian randomization based on constrained maximum likelihood

Mendelian randomization (MR) is a powerful tool for causal inference with observational genome-wide association study (GWAS) summary data. Compared to the more commonly used univariable MR (UVMR), multivariable MR (MVMR) not only is more robust to the notorious problem of genetic (horizontal) pleiotropy but also estimates the direct effect of each exposure on the outcome after accounting for possible mediating effects of other exposures. Despite promising applications, there is a lack of studies on MVMR's theoretical properties and robustness in applications. In this work, we propose an efficient and robust MVMR method based on constrained maximum likelihood (cML), called MVMR-cML, with strong theoretical support. Extensive simulations demonstrate that MVMR-cML performs better than other existing MVMR methods while possessing the above two advantages over its univariable counterpart. An application to several large-scale GWAS summary datasets to infer causal relationships between eight cardiometabolic risk factors and coronary artery disease (CAD) highlights the usefulness and some advantages of the proposed method. For example, after accounting for possible pleiotropic and mediating effects, triglyceride (TG), low-density lipoprotein cholesterol (LDL), and systolic blood pressure (SBP) had direct effects on CAD; in contrast, the effects of high-density lipoprotein cholesterol (HDL), diastolic blood pressure (DBP), and body height diminished after accounting for other risk factors.

Causal association between obstructive sleep apnea and gastroesophageal reflux disease: A bidirectional two-sample Mendelian randomization study

Objectives: Correlations between obstructive sleep apnea (OSA) and gastroesophageal reflux disease (GERD) have been detected in previous observational studies. However, this association remains uncertain due to the potential presence of selection and confounding biases. Therefore, this bidirectional two-sample Mendelian randomization (MR) study was conducted to evaluate the causal relationship between OSA and GERD.

Methods: In this study, instrumental variables (IVs) for OSA were selected from publicly available genetic summary data (27,207 cases and 280,720 controls). Summary statistics for GERD were obtained from a genome-wide association study of 602,604 individuals. The inverse variance weighted (IVW) method was used as the main MR method. The MR-Egger intercept test, MR pleiotropy residual sum and outlier, and leave-one-out analysis were used to detect pleiotropy. Heterogeneity was detected by Cochran’s Q test.

Results: The IVW results revealed that OSA [odds ratio (OR): 1.19, 95% confidence interval (CI): 1.11–1.28, p = 8.88E-07] was causally associated with the incidence of GERD. Moreover, there was evidence of GERD leading to OSA in the IVW analysis (OR: 1.44, 95%CI: 1.33–1.57, p = 7.74E-19). No directional pleiotropy was detected by the MR-Egger intercept test (all p &gt; 0.05).

Conclusion: This study found that OSA is linked to a higher incidence of GERD, and vice versa. This finding might be helpful for the screening and prevention of these two diseases.

Association of smoking, alcohol, and coffee consumption with the risk of ovarian cancer and prognosis: a mendelian randomization study

OBJECTIVE

Currently, the association between smoking, alcohol, and coffee intake and the risk of ovarian cancer (OC) remains conflicting. In this study, we used a two-sample mendelian randomization (MR) method to evaluate the association of smoking, drinking and coffee consumption with the risk of OC and prognosis.

METHODS

Five risk factors related to lifestyles (cigarettes per day, smoking initiation, smoking cessation, alcohol consumption and coffee consumption) were chosen from the Genome-Wide Association Study, and 28, 105, 10, 36 and 36 single-nucleotide polymorphisms (SNPs) were obtained as instrumental variables (IVs). Outcome variables were achieved from the Ovarian Cancer Association Consortium. Inverse-variance-weighted method was mainly used to compute odds ratios (OR) and 95% confidence intervals (Cl).

RESULTS

The two-sample MR analysis supported the causal association of genetically predicted smoking initiation (OR: 1.15 per SD, 95%CI: 1.02-1.29, P = 0.027) and coffee consumption (OR: 1.40 per 50% increase, 95%CI: 1.02-1.93, P = 0.040) with the risk of OC, but not cigarettes per day, smoking cessation, and alcohol consumption. Subgroup analysis based on histological subtypes revealed a positive genetical predictive association between coffee consumption and endometrioid OC (OR: 3.01, 95%CI: 1.50-6.04, P = 0.002). Several smoking initiation-related SNPs (rs7585579, rs7929518, rs2378662, rs10001365, rs11078713, rs7929518, and rs62098013), and coffee consumption-related SNPs (rs4410790, and rs1057868) were all associated with overall survival and cancer-specific survival in OC.

CONCLUSION

Our findings provide the evidence for a favorable causal association of genetically predicted smoking initiation and coffee consumption with OC risk, and coffee consumption is linked to a greater risk of endometrioid OC.

Psychiatric disorders and Type 2 diabetes mellitus: A bidirectional Mendelian randomization

BACKGROUND

Extensive observational evidence put forward the association between psychiatric disorders and type 2 diabetes mellitus (T2DM). However, causal relationships between these two diseases required further research. Thus, we evaluated the bidirection casual effect between five psychiatric disorders and T2DM using two-sample mendelian randomization (MR).

METHODS

By selecting single nucleotide polymorphisms associated with T2DM and five psychiatric disorders (attention-deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), schizophrenia, anxiety disorder and panic disorder), a bidirectional two-sample MR was applied to evaluate causality between these diseases. The inverse-variance weighted (IVW) method was used as the primary analysing approach for estimating possible causal effects. MR-Egger and weighted median were also conducted to verify the results. The funnel plot, Cochran's Q test and MR-Egger intercept test were used for sensitivity analyses. In addition, potential mediators were investigated by risk factor analyses.

RESULTS

Genetic susceptibilities of ADHD and MDD would increase the risk of T2DM (ADHD: OR = 1.14, 95%CI 1.08-1.20; p = 5.7 × 10 - 6 ; MDD: OR = 1.22, 95%CI 1.09-1.36; p = 0.0004 ). In addition, genetic predisposition to T2DM was also associated with ADHD (OR = 1.09, 95%CI 1.04-1.14; p = 0.0004). Several risk factors of T2DM were implicated in the above causal associations, including smoking, high body mass index, waist-to-hip ratio and elevated serum triglycerides.

CONCLUSION

Our studies indicated a causal effect of ADHD and MDD on increasing the risk of T2DM, which was potentially mediated by smoking and obesity-related phenotypes. Meanwhile, we found a causal effect of T2DM on ADHD. Thus, prevention strategies for T2DM should also include mental health and vice versa.

Associations between genetically predicted levels of blood metabolites and pancreatic cancer risk

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive solid malignancies, which is featured by systematic metabolism. Thus, a better understanding of metabolic dysregulation in PDAC is important to better characterize its etiology. Here, we performed a large metabolome-wide association study (MWAS) to systematically explore associations between genetically predicted metabolite levels in blood and PDAC risk. Using data from 881 subjects of European descent in the TwinsUK Project, comprehensive genetic models were built to predict serum metabolite levels. These prediction models were applied to the genetic data of 8275 cases and 6723 controls included in the PanScan (I, II and III) and PanC4 consortia. After assessing the metabolite-PDAC risk associations by a slightly modified TWAS/FUSION framework, we identified five metabolites (including two dipeptides) showing significant associations with PDAC risk at false discovery rate (FDR) <0.05. Integrated with gut microbial information, two-sample Mendelian randomization (MR) analyses were further performed to investigate the relationship among serum metabolites, gut microbiome features and PDAC. The flavonoid-degrading bacteria Flavonifractor sp90199495 was found to be associated with metabolite X-21849 and it was also shown to be associated with PDAC risk. Collectively, our study identified novel candidate metabolites for PDAC risk, which could lead to new insights into the etiology of PDAC and improved treatment options.

Systemic inflammatory regulators and risk of acute-on-chronic liver failure: A bidirectional mendelian-randomization study

Inflammation plays a role in the pathogenesis of acute-on-chronic liver failure (ACLF), however, whether there is a causal relationship between inflammation and ACLF remains unclear. A two-sample Mendelian randomization (MR) approach was used to investigate the causal relationship between systemic inflammatory regulators and ACLF. The study analyzed 41 cytokines and growth factors from 8,293 individuals extracted from a genome-wide association study (GWAS) meta-analysis database involving 253 ACLF cases and 456,095 controls. Our results showed that lower stem cell factor (SCF) levels, lower basic fibroblast growth factor (bFGF) levels and higher Interleukin-13 (IL-13) levels were associated with an increased risk of ACLF (OR = 0.486, 95% CI = 0.264-0.892, p = 0.020; OR = 0.323, 95% CI = 0.107-0.972, p = 0.044; OR = 1.492, 95% CI = 1.111-2.004, p = 0.008, respectively). In addition, genetically predicted ACLF did not affect the expression of systemic inflammatory regulators. Our results indicate that cytokines play a crucial role in the pathogenesis of ACLF. Further studies are needed to determine whether these biomarkers can be used to prevent and treat ACLF.

Age at Menarche Mediating Visceral Adipose Tissue's Influence on Pre-eclampsia: A Mendelian Randomization Study

CONTEXT

The association between visceral adipose tissue (VAT) and pre-eclampsia (PE) shows inconsistent results and the underlying mediator remains unknown.

OBJECTIVE

We aimed to explore the causal effect of VAT on PE risks and the mediation role of age at menarche (AAM) in explaining this relationship.

METHODS

Summary data for PE were obtained from the FinnGen genome-wide association study (3556 cases and 114 735 controls). For exposure data, 70 genetic variants associated with the predicted VAT in 161 149 European women from UK Biobank were used as instrumental variables. Inverse variance weighted and multiple sensitivity analyses were applied. We also conducted multivariable Mendelian randomization (MR) analyses to test the association between VAT-associated single-nucleotide variations and PE. Next, mediation analyses were performed to study whether the association between VAT and PE was mediated via AAM.

RESULTS

In univariable MR analysis, higher volume of VAT was associated with the advancement of AAM and increased PE risk (beta = -0.33; 95% CI, -0.49 to -0.16 for AAM; odds ratio 1.65, 95% CI, 1.23 to 2.20 for PE). After adjusting for waist circumference, waist to hip ratio, and hip circumference, the multivariable MR results presented the consistent positive causality of VAT on PE. Two-step MR analysis proved an estimated 14.3% of the positive effect of VAT on PE was mediated by AAM.

CONCLUSION

Our findings provided evidence of the causal relationship between VAT and PE and proved VAT could accelerate AAM and then contribute to the risk of incident PE.

Constipation and cardiovascular disease: A two-sample Mendelian randomization analysis

Background

Although several observational studies have suggested positive associations between constipation and cardiovascular disease (CVD), a solid causal association has not been demonstrated. Therefore, a two-sample Mendelian randomization (MR) study was performed to investigate the causal associations between constipation and CVD.

Methods

Independent genetic variants strongly associated with constipation were obtained from the FinnGen consortium. Summary-level data for CVD, including coronary artery disease (CAD), myocardial infarction (MI), heart failure (HF), atrial fibrillation (AF), stroke, and its subtypes, were collected from a few extensive genome-wide association studies (GWASs). The inverse-variance weighted methods, weighted median, and MR-Egger were used for the MR estimates. The Cochran's Q test, MR-Egger intercept tests, MR-PRESSO, MR Steiger test, leave-one-out analyses, and funnel plot were used in the sensitivity analysis.

Results

Genetically determined constipation was suggestively associated with AF risk (odds ratio (OR), 1.07; 95% confidence interval (CI), 1.01, 1.14; p = 0.016). Constipation and other CVD do not appear to be causally related. It was demonstrated that the results were robust through sensitivity analyses.

Conclusion

This MR study demonstrated suggestive causal associations of constipation on AF, despite no associations achieving a significance value after multiple testing corrections. There was no evidence of an association between constipation and the risk of CAD, MI, HF, stroke, or stroke subtypes.

Sleep traits, fat accumulation, and glycemic traits in relation to gastroesophageal reflux disease: A Mendelian randomization study

Background

Sleep traits, fat accumulation, and glycemic traits are associated with gastroesophageal reflux disease (GERD) in observational studies. However, whether their associations are causal remains unknown. We performed a Mendelian randomization (MR) study to determine these causal relationships.

Methods

Independent genetic variants associated with insomnia, sleep duration, short sleep duration, body fat percentage, visceral adipose tissue (VAT) mass, type 2 diabetes, fasting glucose, and fasting insulin at the genome-wide significance level were selected as instrumental variables. Summary-level data for GERD were derived from a genome-wide association meta-analysis including 78,707 cases and 288,734 controls of European descent. Inverse variance weighted (IVW) was used for the main analysis, with weighted median and MR-Egger as complements to IVW. Sensitivity analyses were performed using Cochran's Q test, MR-Egger intercept test, and leave-one-out analysis to estimate the stability of the results.

Results

The MR study showed the causal relationships of genetically predicted insomnia (odds ratio [OR] = 1.306, 95% confidence interval [CI] 1.261 to 1.352; p = 2.24 × 10^-51), short sleep duration (OR = 1.304, 95% CI: 1.147 to 1.483, p = 4.83 × 10^-5), body fat percentage (OR = 1.793, 95% CI 1.496 to 2.149; p = 2.68 × 10^-10), and visceral adipose tissue (OR = 2.090, 95% CI 1.963 to 2.225; p = 4.42 × 10^-117) with the risk of GERD. There was little evidence for causal associations between genetically predicted glycemic traits and GERD. In multivariable analyses, genetically predicted VAT accumulation, insomnia, and decreased sleep duration were associated with an increased risk of GERD.

Conclusion

This study suggests the possible roles of insomnia, short sleep, body fat percentage, and visceral adiposity in the development of GERD.

Potential causal association between leisure sedentary behaviors, physical activity and musculoskeletal health: A Mendelian randomization study

BACKGROUND

Increasing evidence shows that leisure sedentary behaviors (LSB) and physical activity (PA) are associated with various musculoskeletal disorders. However, the causality between LSB/PA and musculoskeletal health remained unknown. In this study, we aimed to evaluate the causal relationships between LSB/PA and lower back pain (LBP), intervertebral disc disorder (IVDD), rheumatoid arthritis (RA), and bone mineral density (BMD) by using a two-sample Mendelian randomization method.

METHODS

The exposure data were obtained from large-scale genome-wide association studies (GWAS), including the PA dataset (self-reported PA, n  =  377,234; accelerometer-assessed PA, n  = 91,084) and LSB dataset (n  =  422,218). The outcome data were derived from the FinnGen LBP dataset (n  =  248,528), FinnGen IVDD dataset (n  =  256,896), BMD GWAS dataset (n  =  56,284), and RA GWAS dataset (n  =  58,284). The causal relationships were estimated with inverse variance weighted (IVW), MR-Egger, and weighted median methods. Sensitivity analyses were performed with Cochran's Q test, MR-Egger intercept test, and leave-one-out analysis to estimate the robustness of our findings.

RESULTS

Genetically predicted leisure television watching increased the risk of LBP (OR = 1.68, 95% CI 1.41 to 2.01; P = 8.23×10-9) and IVDD (OR  =  1.62, 95% CI 1.37 to 1.91; P  =  2.13 × 10-8). In addition, this study revealed a potential causal relationship between computer use and a reduced risk of IVDD (OR  = 0.60, 95% CI 0.42 to 0.86; P  =  0.005) and RA (OR  = 0.28, 95% CI 0.13 to 0.60; P  =  0.001).

CONCLUSIONS

Our results suggest that leisure television watching is a risk factor for LBP and IVDD, whereas leisure computer use may act as a protective factor against IVDD and RA. These findings emphasized the importance of distinguishing between different sedentary behaviors in musculoskeletal disease studies.

Investigating the causal mediating effect of type 2 diabetes on the relationship between traits and systolic blood pressure: A two-step Mendelian randomization study

Background

Type 2 diabetes mellitus (T2DM) and hypertension commonly coexist, and we presumed that T2DM might mediate the relationship between some shared risk factors and systolic blood pressure (SBP).

Methods

The causal association between T2DM and SBP was first confirmed using Mendelian randomization (MR) analyses, and a two-step MR design was then used to test the causal mediating effect of T2DM on the relationship between 107 traits and SBP using summary statistics from genome-wide association studies.

Results

T2DM was causally associated with SBP. The univariable MR of the two-step causal mediation analyses suggested that 44 and 45 of the 107 traits had causal associations with T2DM and SBP, respectively. Five of the 27 traits that were significantly associated with both T2DM and SBP could not be reversely altered by T2DM and were included in the second step of the causal mediation analyses. The results indicated that most of the investigated traits causally altered SBP independent of T2DM, but the partial causal mediating effect of T2DM on the association between fasting insulin and SBP was successfully identified with a mediation proportion of 33.6%.

Conclusions

Our study provides novel insights into the role of risk factors in the comorbidity of T2DM and high blood pressure, which is important for long-term disease prevention and management.

A Mendelian randomization-based exploration of red blood cell distribution width and mean corpuscular volume with risk of hemorrhagic strokes

Red blood cell distribution width (RCDW) and mean corpuscular volume (MCV) are associated with different risk factors for hemorrhagic stroke. However, whether RCDW and MCV are causally related to hemorrhagic stroke remains poorly understood. Therefore, we explored the causality between RCDW/MCV and nontraumatic hemorrhagic strokes using Mendelian randomization (MR) methods. We extracted exposure and outcome summary statistics from the UK Biobank and FinnGen. We evaluated the causality of RCDW/MCV on four outcomes (subarachnoid hemorrhage [SAH], intracerebral hemorrhage [ICH], nontraumatic intracranial hemorrhage [nITH], and a combination of SAH, cerebral aneurysm, and aneurysm operations) using univariable MR (UMR) and multivariable MR (MVMR). We further performed colocalization and mediation analyses. UMR and MVMR revealed that higher genetically predicted MCV is protective of ICH (UMR: odds ratio [OR] = 0.89 [0.8-0.99], p = 0.036; MVMR: OR = 0.87 [0.78-0.98], p = 0.021) and nITH (UMR: OR = 0.89 [0.82-0.97], p = 0.005; MVMR: OR = 0.88 [0.8-0.96], p = 0.004). There were no strong causal associations between RCDW/MCV and any other outcome. Colocalization analysis revealed a shared causal variant between MCV and ICH; it was not reported to be associated with ICH. Proportion mediated via diastolic blood pressure was 3.1% (0.1%,14.3%) in ICH and 3.4% (0.2%,15.8%) in nITH. The study constitutes the first MR analysis on whether genetically elevated RCDW and MCV affect the risk of hemorrhagic strokes. UMR, MVMR, and mediation analysis revealed that MCV is a protective factor for ICH and nITH, which may inform new insights into the treatments for hemorrhagic strokes.

Thyroid Stimulating Hormone Levels Are Associated With Genetically Predicted Nonalcoholic Fatty Liver Disease

CONTEXT

Findings from observational studies indicate an association of thyroid hormone levels with the risk of nonalcoholic fatty liver disease (NAFLD); however, conflicting results remain and reverse causality may be a possibility.

OBJECTIVE

This study aimed to evaluate the associations between NAFLD and both plasma thyroxine (T4) and thyroid stimulating hormone (TSH) at the phenotypic and genetic levels.

METHODS

We included 14 797 participants, aged 20 to 74 years who had undergone abdominal ultrasonography during the Third National Health and Nutrition Examination Survey (NHANES III). Multivariable logistic regression analyses were used to examine the observational associations of TSH and T4 with NAFLD. Mediation analyses were performed to study whether the relationship between NAFLD and TSH levels was mediated via potential confounders. A bidirectional, two-sample Mendelian randomization (MR) analysis was used to determine the potential causal relationship.

RESULTS

Multivariable logistic regression model suggested a "dose-response" relationship between TSH (Q4 vs Q1: OR = 1.29; 95% CI, 1.10-1.52; Ptrend = 0.001) and NAFLD. BMI and ALT partially mediated the association between TSH and NAFLD, while the proportion of the mediation effects of BMI and ALT were 39.1% and 22.3%, respectively. In MR analyses, the inverse-variance weighted method was selected as primary method and suggested a putative causal effect of NAFLD on serum TSH levels (OR = 1.022; 95% CI, 1.002-1.043). The result was further validated in the sensitivity analyses.

CONCLUSION

Circulating TSH levels were associated with the risk of NAFLD. MR analysis suggested a putative causal effect of NAFLD on TSH levels.

An IFI6-based hydrogel promotes the healing of radiation-induced skin injury through regulation of the HSF1 activity

Radiation-induced skin injury (RISI) is a common complication of radiotherapy. Interferon-alpha inducible protein 6 (IFI6) significantly reduces the radiation sensitivity of HaCaT cells. Sodium alginate (SA) has substantial moisturizing properties. Graphene oxide (GO) is a suitable substrate with physical antibacterial properties. Therefore, we designed materials to modify IFI6 using the biogule of polydopamine (PDA) connected to GO/SA. The structure, size, morphology, and elemental compositions of IFI6-PDA@GO/SA were analyzed. Cytological studies suggested that IFI6-PDA@GO/SA is non-toxic to HaCaT cells, with antibacterial properties. It promotes migration and vascularization and inhibits apoptosis. These cells express IFI6 after irradiation. The mouse model suggested that IFI6-PDA@GO/SA promotes wound healing and reduces reactive oxygen species expression. IFI6-PDA@GO/SA accelerates RISI healing, possibly by initiating the SSBP1/HSF1 signaling pathway. In addition, IFI6-PDA@GO/SA improves the immune microenvironment. This study constitutes the first use of IFI6 as a RISI wound-healing material.

Causal relationship between physical activity, leisure sedentary behaviors and COVID-19 risk: a Mendelian randomization study

BACKGROUND

The 2019 coronavirus disease pandemic (COVID-19) poses an enormous threat to public health worldwide, and the ensuing management of social isolation has greatly decreased opportunities for physical activity (PA) and increased opportunities for leisure sedentary behaviors (LSB). Given that both PA and LSB have been established as major influencing factors for obesity, diabetes and cardiometabolic syndrome, whether PA/LSB in turn affects the susceptibility to COVID-19 by disrupting metabolic homeostasis remains to be explored. In this study, we aimed to systematically evaluate the causal relationship between PA/LSB and COVID-19 susceptibility, hospitalization and severity using a Mendelian randomization study.

METHODS

Data were obtained from a large-scale PA dataset (N = 377,000), LSB dataset (N = 422,218) and COVID-19 Host Genetics Initiative (N = 2,586,691). The causal effects were estimated with inverse variance weighted, MR-Egger, weighted median and MR-PRESSO. Sensitivity analyses were implemented with Cochran's Q test, MR-Egger intercept test, MR-PRESSO, leave-one-out analysis and the funnel plot. Risk factor analyses were further conducted to investigate the potential mediators.

RESULTS

Genetically predicted accelerometer-assessed PA decreased the risk for COVID-19 hospitalization (OR = 0.93, 95% CI 0.88-0.97; P = 0.002), while leisure television watching significantly increased the risk of COVID-19 hospitalization (OR = 1.55, 95% CI 1.29-1.88; P = 4.68 × 10-6) and disease severity (OR = 1.85, 95% CI 1.33-2.56; P = 0.0002) after Bonferroni correction. No causal effects of self-reported moderate to vigorous physical activity (MVPA), accelerometer fraction of accelerations > 425 milligravities, computer use or driving on COVID-19 progression were observed. Risk factor analyses indicated that the above causal associations might be mediated by several metabolic risk factors, including smoking, high body mass index, elevated serum triglyceride levels, insulin resistance and the occurrence of type 2 diabetes.

CONCLUSION

Our findings supported a causal effect of accelerometer-assessed PA on the reduced risk of COVID-19 hospitalization as well as television watching on the increased risk of COVID-19 hospitalization and severity, which was potentially mediated by smoking, obesity and type 2 diabetes-related phenotypes. Particular attention should be given to reducing leisure sedentary behaviors and encouraging proper exercise during isolation and quarantine for COVID-19.

A practical problem with Egger regression in Mendelian randomization

Mendelian randomization (MR) is an instrumental variable (IV) method using genetic variants such as single nucleotide polymorphisms (SNPs) as IVs to disentangle the causal relationship between an exposure and an outcome. Since any causal conclusion critically depends on the three valid IV assumptions, which will likely be violated in practice, MR methods robust to the IV assumptions are greatly needed. As such a method, Egger regression stands out as one of the most widely used due to its easy use and perceived robustness. Although Egger regression is claimed to be robust to directional pleiotropy under the instrument strength independent of direct effect (InSIDE) assumption, it is known to be dependent on the orientations/coding schemes of SNPs (i.e. which allele of an SNP is selected as the reference group). The current practice, as recommended as the default setting in some popular MR software packages, is to orientate the SNPs to be all positively associated with the exposure, which however, to our knowledge, has not been fully studied to assess its robustness and potential impact. We use both numerical examples (with both real data and simulated data) and analytical results to demonstrate the practical problem of Egger regression with respect to its heavy dependence on the SNP orientations. Under the assumption that InSIDE holds for some specific (and unknown) coding scheme of the SNPs, we analytically show that other coding schemes would in general lead to the violation of InSIDE. Other related MR and IV regression methods may suffer from the same problem. Cautions should be taken when applying Egger regression (and related MR and IV regression methods) in practice.

Investigating the Association Between Seven Sleep Traits and Nonalcoholic Fatty Liver Disease: Observational and Mendelian Randomization Study

Background and Aim: Aberrant sleep parameters are associated with the risk of nonalcoholic fatty liver disease (NAFLD). However, existing information is inconsistent among studies and involves reverse causation. Therefore, we aimed to investigate the observational associations and causations between sleep traits and NAFLD. Methods: We performed multivariable regression to assess observational associations of seven sleep traits (sleep duration, easiness of getting up in the morning, chronotype, nap during day, snoring, insomnia, and narcolepsy), and NAFLD in the UK Biobank (1,029 NAFLD). The Cox proportional hazards model was applied to derive hazard ratios and 95% confidence intervals (CIs). Furthermore, a bidirectional two-sample Mendelian randomization (MR) approach was used to explore the causal relationships between sleep traits and NAFLD. Results: In the multivariable regression model adjusted for potential confounders, getting up in the morning not at all easy (HR, 1.51; 95% CI, 1.27-1.78) and usually insomnia (HR, 1.46; 95% CI, 1.21-1.75) were associated with the risk of NAFLD. Furthermore, the easiness of getting up in the morning and insomnia showed a dose-response association with NAFLD (P_trend <0.05). MR analysis found consistent causal effects of NAFLD on easiness of getting up in the morning (OR, 0.995; 95% CI, 0.990-0.999; p = 0.033) and insomnia (OR, 1.006; 95% CI, 1.001-1.011; p = 0.024). These results were robust to weak instrument bias, pleiotropy, and heterogeneity. Conclusions: Findings showed consistent evidence of observational analyses and MR analyses that trouble getting up in the morning and insomnia were associated with an increased risk of NAFLD. Bidirectional MR demonstrated causal effects of NAFLD on sleep traits.