Conducting a Reproducible Mendelian Randomization Analysis Using the R Analytic Statistical Environment

Phenotypic and genetic characteristics of retinal vascular parameters and their association with diseases

Fundus images allow for non-invasive assessment of the retinal vasculature whose features provide important information on health. Using a fully automated image processing pipeline, we extract 17 different morphological vascular phenotypes, including median vessels diameter, diameter variability, main temporal angles, vascular density, central retinal equivalents, the number of bifurcations, and tortuosity, from over 130,000 fundus images of close to 72,000 UK Biobank subjects. We perform genome-wide association studies of these phenotypes. From this, we estimate their heritabilities, ranging between 5 and 25%, and genetic cross-phenotype correlations, which mostly mirror the corresponding phenotypic correlations, but tend to be slightly larger. Projecting our genetic association signals onto genes and pathways reveals remarkably low overlap suggesting largely decoupled mechanisms modulating the different phenotypes. We find that diameter variability, especially for the veins, associates with diseases including heart attack, pulmonary embolism, and age of death. Mendelian Randomization analysis suggests a causal influence of blood pressure and body mass index on retinal vessel morphology, among other results. We validate key findings in two independent smaller cohorts. Our analyses provide evidence that large-scale analysis of image-derived vascular phenotypes has sufficient power for obtaining functional and causal insights into the processes modulating the retinal vasculature.

Associations of intestinal diseases with anal diseases: a Mendelian randomization study

Although observational clinical studies have established an association between Intestinal Diseases (IDS) and Anal Diseases (ADS), the causal relationship is still not fully understood due to the limitations of observational studies. Genome-wide association study (GWAS) statistical data for IDS and ADS were obtained from publicly available databases. To assess the causal effects of IDS on ADS, we conducted Mendelian randomization analysis. The inverse variance weighted method indicated that Inflammatory bowel disease (IBD) had a significant causal relationship with three kinds of ADS: Anorectal abscess (ARB), Haemorrhoidal disease (HEM), and Fissure and fistula of anal and rectal regions (FISSANAL). Crohn's disease (CD) and Ulcerative colitis (UC) also showed significant causal effects with three ADS: ARB, HEM, and FISSANAL. Furthermore, a potential link between CD and BNA(Benign neoplasm of anus and anal canal), Irritable bowel syndrome (IBS) and HEM, Colorectal cancer (CRC) and BNA, and Celiac disease and MNA (Malignant neoplasm of anus and anal canal) was observed. This comprehensive MR analysis highlight the significant and increased risk of common Anal Diseases (ARB, FISSANAL, and HEM) in patients with IBD, CD, and UC. Additionally, potential positive causal associations emerged between IBS and HEM, CRC and BNA, as well as between celiac disease and MNA.

Causal relationships of infection with Helicobacter pylori and herpesvirus on periodontitis: A Mendelian randomization study

Background

To explore the causal association between Helicobacter pylori (H. pylori) infection, herpesvirus infection and periodontitis (PD) from a genetic perspective using Mendelian randomization (MR).

Methods

The PD data were derived from genome-wide association study (GWAS) from the Dental Endpoints (GLIDE) consortium, and the FinnGen Biobank provided data on H. pylori and herpesvirus infections. In addition, we examined GWAS data for subtypes of H. pylori and herpesvirus infection. Inverse variance weighting (IVW) was selected as a major analysis technique, and weighted median (WM), weighted model, simple model, and MR-Egger regression were added as supplementary methods. To verify the findings, the effects of pleiotropy and heterogeneity were assessed.

Results

Genetically predicted H. pylori infection (OR = 0.914, 95%CI = 0.693-1.205, P = 0.523), anti-H. pylori VacA (OR = 0.973, 95%CI = 0.895-1.057, P = 0.515), anti-H. pylori CagA (OR = 1.072, 95%CI = 0.986-1.164; P = 0.102), Epstein-Barr virus (EBV) infection (OR = 1.026, 95%CI = 0.940-1.120, P = 0.567), Herpes simplex virus (HSV) infection (OR = 0.962, 95%CI = 0.883-1.048, P = 0.372), cytomegalovirus (CMV) infection (OR = 1.025, 95%CI = 0.967-1.088, P = 0.415), EBV nuclear antigen-1 (EBNA1) (OR = 1.061, 95%CI = 0.930-1.209, P = 0.378), EBV virus capsid antigen (VCA) (OR = 1.043, 95CI% = 0.890-1.222, P = 0.603), HSV-1 (OR = 1.251, 95%CI = 0.782-2.001, P = 0.351), HSV-2 (OR = 1.020, 95%CI = 0.950-1.096, P = 0.585), CMV IgG (OR = 0.990, 95CI% = 0.882-1.111, P = 0.861) were not associated with PD, indicated that H. pylori and herpesvirus infection had no causal relationship to PD. Reverse studies also found no cause effect of PD on H. pylori or herpesvirus infection. The results of the sensitivity analysis suggested the robustness of the MR results.

Conclusion

This study offered preliminary proof that H. pylori and herpesvirus infections were not causally linked to PD, and vice versa. However, more robust instrumental variables (IVs) and larger samples of GWAS data were necessary for further MR analysis.

Unraveling NEK4 as a Potential Drug Target in Schizophrenia and Bipolar I Disorder: A Proteomic and Genomic Approach

BACKGROUND AND HYPOTHESIS

Investigating the shared brain protein and genetic components of schizophrenia (SCZ) and bipolar I disorder (BD-I) presents a unique opportunity to understand the underlying pathophysiological processes and pinpoint potential drug targets.

STUDY DESIGN

To identify overlapping susceptibility brain proteins in SCZ and BD-I, we carried out proteome-wide association studies (PWAS) and Mendelian Randomization (MR) by integrating human brain protein quantitative trait loci with large-scale genome-wide association studies for both disorders. We utilized transcriptome-wide association studies (TWAS) to determine the consistency of mRNA-protein dysregulation in both disorders. We applied pleiotropy-informed conditional false discovery rate (pleioFDR) analysis to identify common risk genetic loci for SCZ and BD-I. Additionally, we performed a cell-type-specific analysis in the human brain to detect risk genes notably enriched in distinct brain cell types. The impact of risk gene overexpression on dendritic arborization and axon length in neurons was also examined.

STUDY RESULTS

Our PWAS identified 42 proteins associated with SCZ and 14 with BD-I, among which NEK4, HARS2, SUGP1, and DUS2 were common to both conditions. TWAS and MR analysis verified the significant risk gene NEK4 for both SCZ and BD-I. PleioFDR analysis further supported genetic risk loci associated with NEK4 for both conditions. The cell-type specificity analysis revealed that NEK4 is expressed on the surface of glutamatergic neurons, and its overexpression enhances dendritic arborization and axon length in cultured primary neurons.

CONCLUSIONS

These findings underscore a shared genetic origin for SCZ and BD-I, offering novel insights for potential therapeutic target identification.

Assessment of causal association between the socio-economic status and osteoporosis and fractures: a bidirectional Mendelian randomization study in European population

The correlation between socio-economic status (SES) and bone-related diseases garners increasing attention, prompting a bidirectional Mendelian randomization (MR) analysis in this study. Genetic data on SES indicators (average total household income before tax, years of schooling completed, and Townsend Deprivation Index at recruitment), femoral neck bone mineral density (FN-BMD), heel bone mineral density (eBMD), osteoporosis, and five different sites of fractures (spine, femur, lower leg-ankle, foot, and wrist-hand fractures) were derived from genome-wide association summary statistics of European ancestry. The inverse variance weighted method was employed to obtain the causal estimates, complemented by alternative MR techniques, including MR-Egger, weighted median, and MR-pleiotropy residual sum and outlier (MR-PRESSO). Furthermore, sensitivity analyses and multivariable MR were performed to enhance the robustness of our findings. Higher educational attainment exhibited associations with increased eBMD (β: .06, 95% confidence interval [CI]: 0.01-0.10, P = 7.24 × 10-3), and reduced risks of osteoporosis (OR: 0.78, 95% CI: 0.65-0.94, P = 8.49 × 10-3), spine fracture (OR: 0.76, 95% CI: 0.66-0.88, P = 2.94 × 10-4), femur fracture (OR: 0.78, 95% CI: 0.67-0.91, P = 1.33 × 10-3), lower leg-ankle fracture (OR: 0.79, 95% CI: 0.70-0.88, P = 2.05 × 10-5), foot fracture (OR: 0.78, 95% CI: 0.66-0.93, P = 5.92 × 10-3), and wrist-hand fracture (OR: 0.83, 95% CI: 0.73-0.95, P = 7.15 × 10-3). Material deprivation appeared to increase the risk of spine fracture (OR: 2.63, 95% CI: 1.43-4.85, P = 1.91 × 10-3). A higher FN-BMD level positively affected increased household income (β: .03, 95% CI: 0.01-0.04, P = 6.78 × 10-3). All these estimates were adjusted for body mass index, type 2 diabetes, smoking initiation, and frequency of alcohol intake. The MR analyses show that higher educational levels is associated with higher eBMD, reduced risk of osteoporosis and fractures, while material deprivation is positively related to spine fracture. Enhanced FN-BMD correlates with increased household income. These findings provide valuable insights for health guideline formulation and policy development.

Gastroesophageal reflux disease and risk for arrhythmias: a Mendelian randomization analysis

Background

Clinical observations and epidemiological studies suggest a potential linkage between gastroesophageal reflux disease (GERD) and arrhythmias, yet the underlying mechanism remains elusive. This study investigates the causal relationship between GERD and four types of arrhythmia through a genetic lens, employing Mendelian randomization analysis to elucidate the directionality of these associations.

Methods

Selected single nucleotide polymorphisms (SNPs) from genome-wide association study (GWAS) data were utilized as instrumental variables. The inverse variance weighting (IVW) method, MR-Egger regression analysis, and the weighted median method were employed in two-sample Mendelian randomization analysis. Horizontal pleiotropy was detected and corrected using the MR-PRESSO test and MR-Egger regression. The stability and reliability of the Mendelian randomization results were assessed using the leave-one-out method, Cochran's Q test, and funnel plots. The causal relationship between GERD and four types of arrhythmias was evaluated using the odds ratio (OR).

Results

IVW results indicated that GERD could increase the risk of arrhythmias. A one standard deviation increases in the logarithmically transformed GERD score resulted in a 34% increase in the risk of arrhythmia (OR = 1.34; 95% CI 1.19-1.51; p = 1.66E-06). No significant correlation was found between GERD and other arrhythmias.

Conclusion

A causal relationship exists between GERD and arrhythmias, suggesting that GERD increases the risk of developing these arrhythmias.

Inflammatory bowel disease and risk for hemorrhoids: a Mendelian randomization analysis

Observational studies have reported an association between inflammatory bowel disease (IBD), which includes Crohn's disease (CD) and ulcerative colitis (UC), and hemorrhoids (HEM). However, the presence of a causal relationship within this observed association remains to be confirmed. Consequently, we utilized the Mendelian randomization (MR) method to assess the causal effects of IBD on hemorrhoids. We validated the association between IBD and hemorrhoids in humans based on genome-wide association studies (GWAS) data. To investigate the causal relationship between IBD and hemorrhoids, we performed a two-sample Mendelian randomization study using training and validation sets. The genetic variation data for IBD, CD, UC, and hemorrhoids were derived from published genome-wide association studies (GWAS) of individuals of European. Two-sample Mendelian randomization and Multivariable Mendelian randomization (MVMR) were employed to determine the causal relationship between IBD (CD or UC) and hemorrhoids. Genetically predicted overall IBD was positively associated with hemorrhoids risk, with ORs of 1.02 (95% CIs 1.01-1.03, P = 4.39 × 10-4) and 1.02 (95% CIs 1.01-1.03, P = 4.99 × 10-5) in the training and validation sets, respectively. Furthermore, we found that CD was positively associated with hemorrhoids risk, with ORs of 1.02 (95% CIs 1.01-1.03, P = 4.12 × 10-6) and 1.02 (95% CIs 1.01-1.02, P = 3.78 × 10-5) for CD in the training and validation sets, respectively. In addition, we found that UC in the training set was positively associated with hemorrhoids risk (ORs 1.02, 95% CIs 1.01-1.03, P = 4.65 × 10-3), while no significant causal relationship between UC and hemorrhoids was shown in the validation set (P > 0.05). However, after MVMR adjustment, UC in the training set was not associated with an increased risk of hemorrhoids. Our study showed that there is a causal relationship between CD and hemorrhoids, which may suggest that clinicians need to prevent the occurrence of hemorrhoids in CD patients.

Identification and validation of potential genes for the diagnosis of sepsis by bioinformatics and 2-sample Mendelian randomization study

This integrated study combines bioinformatics, machine learning, and Mendelian randomization (MR) to discover and validate molecular biomarkers for sepsis diagnosis. Methods include differential expression analysis, weighted gene co-expression network analysis (WGCNA) for identifying sepsis-related modules and hub genes, and functional enrichment analyses to explore the roles of hub genes. Machine learning algorithms identify 3 diagnostic genes - CD177, LDHA, and MCEMP1 - consistently highly expressed in sepsis patients. The nomogram model effectively predicts sepsis risk, supported by receiver operator characteristic (ROC) curves. Correlations between diagnostic genes and immune cell infiltration are observed. MR analysis reveals a positive causal relationship between MCEMP1 and sepsis risk. In conclusion, this study presents potential sepsis diagnostic biomarkers, highlighting the genetic association of MCEMP1 with sepsis for insights into early diagnosis.

CAUSAL ROLES OF SERUM URIC ACID LEVELS AND GOUT IN SEPSIS: A MENDELIAN RANDOMIZATION STUDY

ABSTRACT

Objective: Several epidemiological studies have identified a potential link between serum uric acid (UA), gout, and sepsis. The primary objective of this study is to delve deeper into this connection, investigating the causal effect of UA and gout on sepsis by applying Mendelian randomization (MR). Methods: The causal relationship was analyzed using data from Genome-Wide Association Study (GWAS). Inverse variance weighting (IVW) was used as the main analysis method. Three complementary methods were used for our MR analysis, which included the MR-Egger regression method, the weighted median method, the simple median method. Horizontal pleiotropy was identified by MR-Egger intercept test. Cochran's Q statistics were employed to assess the existence of instrument heterogeneity. The leave-one-out method was used as a sensitivity analysis. Results: The IVW results indicated that there was a positive causal relationship between UA and sepsis (critical care) (odds ratio [OR] = 0.24, 95% confidence interval [CI]: 0.04 to 0.43, P = 0.018, F = 4,291.20). There was no significant association between UA and sepsis (28-day death in critical care) (OR = 0.10, 95% CI = -0.29 to 0.50, P = 0.604). There was no significant association between gout and sepsis (critical care) (OR = 0.85, 95% CI = -4.87 to 6.57, P = 0.771), and sepsis (28-day death in critical care) (OR = -6.30, 95% CI = -17.41 to 4.81, P = 0.267). Horizontal pleiotropy was absent in this study. The results were robust under all sensitivity analyses. Conclusion: The study revealed that elevated UA levels were causally linked with sepsis (critical care). No causal relationship had been found between UA and sepsis (28-day death in critical care), as well as between gout and sepsis.

Causal relationship between obesity and anorectal abscess: a Mendelian randomization study

Background

Observational studies have indicated that obesity is a risk factor for anorectal abscess (ARB). However, it remains unclear whether a causal genetic relationship exists between obesity and ARB.

Methods

Univariate and multivariate Mendelian randomization (MR) were conducted using data from a large, published genome-wide association study (GWAS) of European ancestry to infer a causal relationship between obesity and ARB. Inverse variance weighted (IVW) analysis served as the primary analysis method, with results reported as odds ratios (OR).

Results

MR analysis revealed that body mass index (BMI) positively affects ARB (OR 1.974, 95% confidence interval (CI) 1.548-2.519, p = 4.34 × 10-8). The weighted median method (OR = 1.879, 95% CI 1.248-2.829, p = 0.002) and Bayesian model averaging (BMA) (OR = 1.88, 95% CI 1.477-2.392, p = 2.85 × 10-7) also demonstrated consistent results. Subsequently, the impact of several obesity-related characteristics on ARB was assessed. Body fat percentage (BF), whole body fat mass (FM), waist circumference (WC), and hip circumference (HC) were found to be causally associated with an increased risk of ARB. However, these associations vanished after adjusting for BMI effects.

Conclusion

The study confirms a positive causal effect of obesity on ARB, highlighting that reasonable weight control is an important strategy to reduce the incidence of ARB.

Employing Bidirectional Two-Sample Mendelian Randomization Analysis to Verify the Potential of Polyunsaturated Fatty Acid Levels in the Prevention of Pancreatic Cancer

Polyunsaturated fatty acids (PUFAs), specifically Omega-3 (FAω3) and docosahexaenoic acid (DHA), have been studied for their potential role in modulating pancreatic cancer (PC) risk. Although observational studies suggest a beneficial effect in reducing this risk, their findings are often limited by confounding variables and issues of reverse causation. This study used a two-way two-sample Mendelian randomization (MR) method to test the hypothesized genetic causal relationship between PUFAs and PC risk. Data from an extensive genome-wide association study (GWAS) were analyzed, focusing on FAω3 and FAω6 levels, their ratios, and DHA as variables and PC incidence as outcomes. This relationship was comprehensively evaluated using related MR methods, such as inverse variance weighting (IVW), MR Egger, and weighted median (WM). This study finds a significant negative correlation between FAω3 and DHA levels and PC risk, while FAω6 levels show no significant correlation. Interestingly, the ratio of FAω6 to FAω3 was positively associated with increased risk of PC. Neither the MR Egger nor the MR-PRESSO tests detected significant pleiotropy, nor did the Cochrane's Q test show significant heterogeneity. Leave-one-out analyzes further confirmed the robustness of these results. Using MR analysis of two samples, this study provides genetic causal evidence that FAω3 and DHA levels reduce the risk of PC, whereas the ratio of FAω6 to FAω3 increases the risk of PC. These insights highlight the potential utility of supplementing FAω3 and DHA or altering PUFAs in developing PC prevention strategies.

IL27 and IL1RN are causally associated with acute pancreatitis: a Mendelian randomization study

BACKGROUND

The interleukin (IL) plays a role in the development of acute pancreatitis (AP). However, the specific IL in AP has not been fully revealed. Therefore, the association between prospective IL and AP was studied via Mendelian randomization (MR).

METHODS

The HUGO Gene nomenclature committee (HGNC) database provided 47 interleukin related genes (ILRGs). ILRGs and differentially expressed genes (DEGs) from GSE194331 were overlapped to create differently expressed ILRGs (DE-ILRGs). The integrative epidemiology unit (IEU) open genome-wide association study (GWAS) database provided exposure and outcome datasets. Univariate MR (UVMR) analysis using MR-Egger, IVW, simple mode, and weighted mode was done. UVMR results were verified using sensitivity analysis. Drug prediction, MVMR analysis, and PPI network development were also performed.

RESULTS

Six DE-ILRGs were obtained. IL27 and IL1RN were substantially causally linked with AP by UVMR analysis (OR = 0.926, P < 0.001 and OR = 1.031, P = 0.023). Our sensitivity analysis showed the dependability of our results. Direct effect of IL27 was suggested by MVMR analysis. In the cytokine receptor binding pathway, IL27 and IL1RN interacted with IL36G and IL1R2. TAE-684, ARQ-680, and 12 other IL1RN and 14 IL27 medications were predicted.

CONCLUSIONS

IL1RN was identified as a risk factor for acute pancreatitis (AP), but IL27 was found to be a protective factor for AP.

A phenome-wide association and Mendelian randomisation study of alcohol use variants in a diverse cohort comprising over 3 million individuals

BACKGROUND

Alcohol consumption is associated with numerous negative social and health outcomes. These associations may be direct consequences of drinking, or they may reflect common genetic factors that influence both alcohol consumption and other outcomes.

METHODS

We performed exploratory phenome-wide association studies (PheWAS) of three of the best studied protective single nucleotide polymorphisms (SNPs) in genes encoding ethanol metabolising enzymes (ADH1B: rs1229984-T, rs2066702-A; ADH1C: rs698-T) using up to 1109 health outcomes across 28 phenotypic categories (e.g., substance-use, mental health, sleep, immune, cardiovascular, metabolic) from a diverse 23andMe cohort, including European (N ≤ 2,619,939), Latin American (N ≤ 446,646) and African American (N ≤ 146,776) populations to uncover new and perhaps unexpected associations. These SNPs have been consistently implicated by both candidate gene studies and genome-wide association studies of alcohol-related behaviours but have not been investigated in detail for other relevant phenotypes in a hypothesis-free approach in such a large cohort of multiple ancestries. To provide insight into potential causal effects of alcohol consumption on the outcomes significant in the PheWAS, we performed univariable two-sample and one-sample Mendelian randomisation (MR) analyses.

FINDINGS

The minor allele rs1229984-T, which is protective against alcohol behaviours, showed the highest number of PheWAS associations across the three cohorts (N = 232, European; N = 29, Latin American; N = 7, African American). rs1229984-T influenced multiple domains of health. We replicated associations with alcohol-related behaviours, mental and sleep conditions, and cardio-metabolic health. We also found associations with understudied traits related to neurological (migraines, epilepsy), immune (allergies), musculoskeletal (fibromyalgia), and reproductive health (preeclampsia). MR analyses identified evidence of causal effects of alcohol consumption on liability for 35 of these outcomes in the European cohort.

INTERPRETATION

Our work demonstrates that polymorphisms in genes encoding alcohol metabolising enzymes affect multiple domains of health beyond alcohol-related behaviours. Understanding the underlying mechanisms of these effects could have implications for treatments and preventative medicine.

FUNDING

MVJ, NCK, SBB, SSR and AAP were supported by T32IR5226 and 28IR-0070. SSR was also supported by NIDA DP1DA054394. NCK and RBC were also supported by R25MH081482. ASH was supported by funds from NIAAA K01AA030083. JLMO was supported by VA 1IK2CX002095. JLMO and JJMM were also supported by NIDA R21DA050160. JJMM was also supported by the Kavli Postdoctoral Award for Academic Diversity. EGA was supported by K01MH121659 from the NIMH/NIH, the Caroline Wiess Law Fund for Research in Molecular Medicine and the ARCO Foundation Young Teacher-Investigator Fund at Baylor College of Medicine. MSA was supported by the Instituto de Salud Carlos III and co-funded by the European Union Found: Fondo Social Europeo Plus (FSE+) (P19/01224, PI22/00464 and CP22/00128).

A multitrait genetic study of hemostatic factors and hemorrhagic transformation after stroke treatment

BACKGROUND

Thrombolytic recombinant tissue plasminogen activator (r-tPA) treatment is the only pharmacologic intervention available in the ischemic stroke acute phase. This treatment is associated with an increased risk of intracerebral hemorrhages, known as hemorrhagic transformations (HTs), which worsen the patient's prognosis.

OBJECTIVES

To investigate the association between genetically determined natural hemostatic factors' levels and increased risk of HT after r-tPA treatment.

METHODS

Using data from genome-wide association studies on the risk of HT after r-tPA treatment and data on 7 hemostatic factors (factor [F]VII, FVIII, von Willebrand factor [VWF], FXI, fibrinogen, plasminogen activator inhibitor-1, and tissue plasminogen activator), we performed local and global genetic correlation estimation multitrait analyses and colocalization and 2-sample Mendelian randomization analyses between hemostatic factors and HT.

RESULTS

Local correlations identified a genomic region on chromosome 16 with shared covariance: fibrinogen-HT, P = 2.45 × 10-11. Multitrait analysis between fibrinogen-HT revealed 3 loci that simultaneously regulate circulating levels of fibrinogen and risk of HT: rs56026866 (PLXND1), P = 8.80 × 10-10; rs1421067 (CHD9), P = 1.81 × 10-14; and rs34780449, near ROBO1 gene, P = 1.64 × 10-8. Multitrait analysis between VWF-HT showed a novel common association regulating VWF and risk of HT after r-tPA at rs10942300 (ZNF366), P = 1.81 × 10-14. Mendelian randomization analysis did not find significant causal associations, although a nominal association was observed for FXI-HT (inverse-variance weighted estimate [SE], 0.07 [-0.29 to 0.00]; odds ratio, 0.87; 95% CI, 0.75-1.00; raw P = .05).

CONCLUSION

We identified 4 shared loci between hemostatic factors and HT after r-tPA treatment, suggesting common regulatory mechanisms between fibrinogen and VWF levels and HT. Further research to determine a possible mediating effect of fibrinogen on HT risk is needed.

Causal role of metabolites in non-small cell lung cancer: Mendelian randomization (MR) study

On a global scale, lung cancer(LC) is the most commonly occurring form of cancer. Nonetheless, the process of screening and detecting it in its early stages presents significant challenges. Earlier research endeavors have recognized metabolites as potentially reliable biomarkers for LC. However, the majority of these studies have been limited in scope, featuring inconsistencies in terms of the relationships and levels of association observed.Moreover, there has been a lack of consistency in the types of biological samples utilized in previous studies. Therefore, the main objective of our research was to explore the correlation between metabolites and Non-small cell lung cancer (NSCLC).Thorough two-sample Mendelian randomization (TSMR) analysis, we investigated potential cause-and-effect relationships between 1400 metabolites and the risk of NSCLC.The analysis of TSMR revealed a significant causal impact of 61 metabolites on NSCLC.To ensure the reliability and validity of our findings, we perform FDR correction for P-values by Benjaminiand Hochberg(BH) method, Our results indicate that Oleate/vaccenate (18:1) levels and Caffeine to paraxanthine ratio may be causally associated with an increased risk of NSCLC [Oleate/vaccenate(18:1)levels: OR = 1.171,95%CI: 1.085-1.265, FDR = 0.036; Caffeine to paraxanthine ratio: OR = 1.386, 95%CI:1.191-1.612,FDR = 0.032].

Body mass index and risk of inflammatory breast disease: a Mendelian randomization study

Introduction

Introduction: in previous studies, obesity was identified as a risk factor for inflammatory breast disease, but its causality is uncertain. In the present study, we performed a two-sample Mendelian randomization (TSMR) analysis to investigate the causal relationship between obesity and inflammatory breast disease. Methods: we use body mass index (BMI) as a measure of obesity. Data for single nucleotide polymorphisms (SNPs) associated with BMI were obtained from UK Biobank. Data for single nucleotide polymorphisms (SNPs) associated with mastitis were obtained from FinnGen Biobank. We used several MR analysis methods, such as inverse-variance weighting (IVW), MR-Egger, weighted median, simple mode and weighted mode to make our results more convincing. We also performed MR-PRESSO test, MR-Egger test, heterogeneity test, pleiotropy test and leave-one-out analysis to make our analysis results more robust and credible. We used odds ratio (OR) to evaluate the causal relationship between BMI and mastitis. Results: based on the IVW random effects model, we found that a one-standard deviation (SD) increase in BMI increased the risk of mastitis by 62.1 % (OR = 1.621, 95 % CI: 1.262-2.083, p = 1.59E-4), which is almost consistent with the results of several other methods. Conclusions: in European individuals, an increase in the number of BMI increases the risk of inflammatory breast disease. People with high BMI need to control their weight to reduce the incidence of inflammatory breast disease.

Genetic causal association between rheumatoid arthritis and periodontitis: a bidirectional two-sample Mendelian randomization analysis

BACKGROUND

A link between rheumatoid arthritis (RA) and periodontitis (PD) has been established. However, their causal relationship remains unclear and the effect of different serotypes of RA on the PD development has not been investigated. This study aims to elucidate the causal association between PD and different serotypes of RA using Mendelian randomization (MR).

METHODS

A bidirectional two-sample MR analysis was performed using available large-scale genome-wide association studies statistics. The inverse-variance weighted (IVW) or multiplicative random-effects IVW was used to determine causality, depending on the heterogeneity of instrumental variables. Additional sensitivity analyses were also performed.

RESULTS

The forward MR analysis identified that seropositive RA (odds ratio (OR), 1.26; 95% confidence interval (CI), 1.07-1.44; p = 0.0018), but not seronegative RA (OR, 1.01; 95% CI, 0.95-1.06; p = 0.9098), was associated with an increased risk of PD. The reverse MR analysis did not show any significant causal effect of PD on RA, independent of the serotypes. The sensitivity tests (p > 0.05) confirmed the robustness and accuracy of these findings.

CONCLUSION

This study revealed that there was a genetic causal effect of seropositive RA on PD, suggesting that this subtype of RA patients may require specific clinical attention to prevent the development of PD.

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.

Causal associations between prodromal infection and neuromyelitis optica spectrum disorder: A Mendelian randomization study

BACKGROUND AND PURPOSE

Prodromal infections are associated with neuromyelitis optica spectrum disorder (NMOSD), but it remains unclear which type of infection has a causal association with NMOSD. We aimed to explore the causal associations between four herpesvirus infections (chickenpox, cold sores, mononucleosis and shingles) and NMOSD, as well as between other types of infections and NMOSD.

METHODS

For data on infections, we used the genome-wide association study (GWAS) summary statistics from the 23andMe cohort. For outcomes, we used the GWAS data of participants of European ancestry, including 215 NMOSD patients (132 anti-aquaporin-4 antibody [AQP4-ab]-positive patients and 83 AQP4-ab-negative patients) and 1244 normal controls. Single-nucleotide polymorphism (SNP) identification and two-sample Mendelian randomization (MR) analyses were then performed.

RESULTS

In the 23andMe cohort, we identified one SNP for chickenpox (rs9266089 in HLA-B gene), one SNP for cold scores (rs885950 in the POU5F1 gene), one SNP for mononucleosis (rs2596465 in the HCP5 gene), and three SNPs for shingles (rs2523591 in the HLA-B gene; rs7047299 in the IFNA21 gene; rs9260809 in the MICD gene). The association between cold sores and AQP4-ab-positive NMOSD reached statistical significance (odds ratio [OR] 745.318; 95% confidence interval [CI] 22.176, 25,049.53 [p < 0.001, Q < 0.001]). The association between shingles and AQP4-ab-positive NMOSD was also statistically significant (OR 21.073; 95% CI 4.271, 103.974 [p < 0.001, Q < 0.001]). No significant association was observed between other infections and AQP4-ab-positive or AQP4-ab-negative NMOSD.

CONCLUSION

These findings suggest there are positive associations between cold sores and shingles and AQP4-ab-positive NMOSD, indicating there may be causal links between herpes simplex virus and varicella-zoster virus infection and AQP4-ab-positive NMOSD.

Association of BMI and waist circumference with diabetic microvascular complications: A prospective cohort study from the UK Biobank and Mendelian randomization analysis

AIMS

To investigate the precise association between BMI and waist circumference (WC) and diabetic complications, including retinopathy (DR), nephropathy (DN) and peripheral neuropathy (DPN).

METHODS

A multivariable-adjusted Cox proportional hazard model was used to evaluate the observed association from 30,541 UK Biobank participants with diabetes. A two-sample Mendelian randomization (MR) framework was applied to summary-level GWASs of BMI and WC comprising a total of 461,460 and 462,166 participants from UK Biobank to explore the potential causal association.

RESULTS

Higher BMI and WC were associated with increased risks of DR, DN, and DPN (HR (95% CI), per-SD increase: BMI: DR 1.09 (1.04-1.13), DN 1.37 (1.33-1.41), DPN 1.27 (1.20-1.34); WC: DR 1.11 (1.07-1.16), DN 1.41 (1.36-1.46), DPN 1.38 (1.30-1.45)) in the UK Biobank cohort. Univariate MR indicated that increased BMI and WC were causal risk factors for these complications (OR (95% CI), per-SD increase: BMI: DR 1.33 (1.22-1.45), DN 1.74 (1.47-2.07), DPN 2.20 (1.67-2.90); WC: DR 1.43 (1.27-1.61), DN 2.03 (1.62-2.55), DPN 2.80 (1.99-3.92)), and the effect sizes remained significant after adjustment for glycated hemoglobin.

CONCLUSIONS

Prospective observational and MR analyses provided evidence that high BMI and WC may represent potential causal risk factors for diabetic microvascular complications. Weight control might modify the risks of these complications independently of glycemic control and should be considered as a therapeutic recommendation.

A Mendelian randomization-based approach to explore the relationship between leukocyte counts and breast cancer risk in European ethnic groups

Exploring the potential association between peripheral blood leukocyte counts and breast cancer risk by Mendelian randomization (MR) analysis methods. Genetic data related to peripheral blood sorting counts of leukocytes were collected from a genome-wide association study by Blood Cell Consortium (BCX). Single nucleotide polymorphic loci predicting peripheral blood sorting counts of these leukocytes were selected as instrumental variables according to the correlation assumption, independence assumption and exclusivity assumption of MR. The data on breast cancer and its subtypes were obtained from Breast Cancer Association Consortium (BCAC) and FinnGen Consortium. In this study, the Inverse-Variance Weighted (IVW), Weighted Median, MR-Egger, Maximum Likelihood (ML), MR-PRESSO and Constrained Maximum Likelihood and Model Averaging (cML-MA) methods of random effects models were used for MR analysis. Cochran's Q analysis, and MR-Egger intercept analysis were applied for sensitivity analysis. IVW and cML-MA were considered the primary analytical tools, and the results of the other 4 MRs were used as complementary and validation. The results suggest that there is no significant causal relationship between leukocyte count and breast cancer risk (IVW OR = 0.98 [95% CI: 0.93-1.03], p-value = 0.35; CML-MA OR = 1.01 [95% CI: 0.98-1.05], p-value = 0.51). In addition, we analyzed whether there was a potential correlation between the five main types of categorized leukocyte counts and different breast cancer subtypes. We did not find significant evidence to support a significant correlation between leukocyte counts and breast cancer subtypes.

Correlations between schizophrenia and lichen planus: a two-sample bidirectional Mendelian randomization study

Background

Several existing studies have shown a correlation between schizophrenia and lichen planus (LP). However, the causality of this relationship remains uncertain. Thus, this study aimed to examine the causal association between schizophrenia and LP.

Methods

A two-sample Mendelian randomization (MR) study was carried out to investigate whether schizophrenia is causally related to LP and vice versa, and genetic variants in this study were taken from previous genome-wide association studies. We used the inverse variance weighted (IVW) method as the main analysis. Furthermore, several sensitivity analyses were performed to assess heterogeneity, horizontal pleiotropy, and stability.

Results

Our results show that schizophrenia has a protective effect on LP (OR = 0.881, 95%CI = 0.795-0.975, p = 0.015). Conversely, we observed no significant relationship between LP and schizophrenia in reverse MR analysis (OR = 0.934, 95%CI = 0.851-1.026, p = 0.156).

Conclusion

Our two-sample Mendelian randomization study supports a significant causal relationship between LP and schizophrenia and finds that schizophrenia can reduce the incidence of LP. This is in contrast to previous findings and provides new insights into the relationship between LP and schizophrenia, but the exact mechanism needs further investigation.

Among men, androgens are associated with a decrease in Alzheimer's disease risk

INTRODUCTION

Increased levels of sex hormones have been hypothesized to decrease Alzheimer's disease (AD) risk. We assessed the association between sex steroid hormones with AD using a Mendelian randomization (MR) approach.

METHODS

An inverse-variance weighting (IVW) MR analysis was performed using effect estimates from external genome-wide association study (GWAS) summary statistics. We included independent variants (linkage disequilibrium R2  < 0.001) and a p-value threshold of 5 × 10-8 .

RESULTS

An increase in androgens was associated with a decreased AD risk among men: testosterone (odds ratio [OR]: 0.53; 95% confidence interval [CI]: 0.32-0.88; p-value: 0.01; false discovery rate [FDR] p-value: 0.03); dehydroepiandrosterone sulfate (DHEAS; OR: 0.56; 95% CI: 0.38-0.85; p-value: 0.01; FDR p-value: 0.03); and androsterone sulfate (OR: 0.69; 95% CI: 0.46-1.02; p-value: 0.06; FDR p-value: 0.10). There was no association between sex steroid hormones and AD among women, although analysis for estradiol had limited statistical power.

DISCUSSION

A higher concentration of androgens was associated with a decreased risk of AD among men of European ancestry, suggesting that androgens among men might be neuroprotective and could potentially prevent or delay an AD diagnosis.

HIGHLIGHTS

Sex hormones are hypothesized to play a role in developing Alzheimer's disease (AD). The effect of sex hormones on AD was assessed using Mendelian randomization (MR) analysis. Among women, genetically determined effects of sex hormones were limited or null. Among men, a higher concentration of androgens decreased AD risk. This study suggests a causal relationship between androgens and AD among men.

Participation bias in the UK Biobank distorts genetic associations and downstream analyses

While volunteer-based studies such as the UK Biobank have become the cornerstone of genetic epidemiology, the participating individuals are rarely representative of their target population. To evaluate the impact of selective participation, here we derived UK Biobank participation probabilities on the basis of 14 variables harmonized across the UK Biobank and a representative sample. We then conducted weighted genome-wide association analyses on 19 traits. Comparing the output from weighted genome-wide association analyses (neffective = 94,643 to 102,215) with that from standard genome-wide association analyses (n = 263,464 to 283,749), we found that increasing representativeness led to changes in SNP effect sizes and identified novel SNP associations for 12 traits. While heritability estimates were less impacted by weighting (maximum change in h2, 5%), we found substantial discrepancies for genetic correlations (maximum change in rg, 0.31) and Mendelian randomization estimates (maximum change in βSTD, 0.15) for socio-behavioural traits. We urge the field to increase representativeness in biobank samples, especially when studying genetic correlates of behaviour, lifestyles and social outcomes.

Integrating GWAS and proteome data to identify novel drug targets for MU

Mouth ulcers have been associated with numerous loci in genome wide association studies (GWAS). Nonetheless, it remains unclear what mechanisms are involved in the pathogenesis of mouth ulcers at these loci, as well as what the most effective ulcer drugs are. Thus, we aimed to screen hub genes responsible for mouth ulcer pathogenesis. We conducted an imputed/in-silico proteome-wide association study to discover candidate genes that impact the development of mouth ulcers and affect the expression and concentration of associated proteins in the bloodstream. The integrative analysis revealed that 35 genes play a significant role in the development of mouth ulcers, both in terms of their protein and transcriptional levels. Following this analysis, the researchers identified 6 key genes, namely BTN3A3, IL12B, BPI, FAM213A, PLXNB2, and IL22RA2, which were related to the onset of mouth ulcers. By combining with multidimensional data, six genes were found to correlate with mouth ulcer pathogenesis, which can be useful for further biological and therapeutic research.

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.

Bayesian Genetic Colocalization Test of Two Traits Using coloc

Genetic colocalization is an approach for determining whether a genetic variant at a particular locus is shared across multiple phenotypes. Genome-wide association studies (GWAS) have successfully mapped genetic variants associated with thousands of complex traits and diseases. However, a large proportion of GWAS signals fall in non-coding regions of the genome, making functional interpretation a challenge. Colocalization relies on a Bayesian framework that can integrate summary statistics, for example those derived from GWAS and expression quantitative trait loci (eQTL) mapping, to assess whether two or more independent association signals at a region of interest are consistent with a shared causal variant. The results from a colocalization analysis may be used to evaluate putative causal relationships between omics-based molecular measurements and a complex disease, and can generate hypotheses that may be followed up by tailored experiments. In this article, we present an easy and straightforward protocol for conducting a Bayesian test for colocalization of two traits using the 'coloc' package in R with summary-level results derived from GWAS and eQTL studies. We also provide general guidelines that can assist in the interpretation of findings generated from colocalization analyses. © 2022 Wiley Periodicals LLC. Basic Protocol: Performing a genetic colocalization analysis using the 'coloc' package in R and summary-level data Support Protocol: Installing the 'coloc' R package.

Identification of novel proteins for lacunar stroke by integrating genome-wide association data and human brain proteomes

BACKGROUND

Previous genome-wide association studies (GWAS) have identified numerous risk genes for lacunar stroke, but it is challenging to decipher how they confer risk for the disease. We employed an integrative analytical pipeline to efficiently transform genetic associations to identify novel proteins for lacunar stroke.

METHODS

We systematically integrated lacunar stroke genome-wide association study (GWAS) (N=7338) with human brain proteomes (N=376) to perform proteome-wide association studies (PWAS), Mendelian randomization (MR), and Bayesian colocalization. We also used an independent human brain proteomic dataset (N=152) to annotate the new genes.

RESULTS

We found that the protein abundance of seven genes (ICA1L, CAND2, ALDH2, MADD, MRVI1, CSPG4, and PTPN11) in the brain was associated with lacunar stroke. These seven genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and astrocytes. Three genes (ICA1L, CAND2, ALDH2) were causal in lacunar stroke (P < 0.05/proteins identified for PWAS; posterior probability of hypothesis 4 ≥ 75 % for Bayesian colocalization), and they were linked with lacunar stroke in confirmatory PWAS and independent MR. We also found that ICA1L is related to lacunar stroke at the brain transcriptome level.

CONCLUSIONS

Our present proteomic findings have identified ICA1L, CAND2, and ALDH2 as compelling genes that may give key hints for future functional research and possible therapeutic targets for lacunar stroke.

Causal Association Between Periodontitis and Type 2 Diabetes: A Bidirectional Two-Sample Mendelian Randomization Analysis

Background: Previous observational studies have reported a bidirectional association between periodontitis and type 2 diabetes, but the causality of these relationships remains unestablished. We clarified the bidirectional causal association through two-sample Mendelian randomization (MR).

Methods: We obtained summary-level data for periodontitis and type 2 diabetes from several published large-scale genome-wide association studies (GWAS) of individuals of European ancestry. For the casual effect of periodontitis on type 2 diabetes, we used five independent single-nucleotide polymorphisms (SNPs) specific to periodontitis from three GWAS. The summary statistics for the associations of exposure-related SNPs with type 2 diabetes were drawn from the GWAS in the Diabetes Genetics Replication and Meta-analysis (DIAGRAM) consortium and the FinnGen consortium R5 release, respectively. For the reversed causal inference, 132 and 49 SNPs associated with type 2 diabetes from the DIAGRAM consortium and the FinnGen consortium R5 release were included, and the summary-level statistics were obtained from the Gene-Lifestyle Interactions in Dental Endpoints consortium. Multiple approaches of MR were carried out.

Results: Periodontitis was not causally related with the risk of type 2 diabetes (all p > 0.05). No causal effect of type 2 diabetes on periodontitis was found (all p > 0.05). Estimates were consistent across multiple MR analyses.

Conclusion: This study based on genetic data does not support a bidirectional causal association between periodontitis and type 2 diabetes.

Two-Sample Multivariable Mendelian Randomization Analysis Using R

Mendelian randomization is a framework that uses measured variation in genes for assessing and estimating the causal effect of an exposure on an outcome. Multivariable Mendelian randomization is an extension that can assess the causal effect of multiple exposures on an outcome, and can be advantageous when considering a set (>1) of potentially correlated candidate risk factors in evaluating the causal effect of each on a health outcome, accounting for measured pleiotropy. This can be seen, for example, in determining the causal effects of lipids and cholesterol on type 2 diabetes risk, where the correlated risk factors share genetic predictors. Similar to univariate Mendelian randomization, multivariable Mendelian randomization can be conducted using two-sample summary-level data where the gene-exposure and gene-outcome associations are derived from separate samples from the same underlying population. Here, we present a protocol for conducting a two-sample multivariable Mendelian randomization study using the 'MVMR' package in R and summary-level genetic data. We also provide a protocol for searching and obtaining instruments using available data sources in the 'MRInstruments' R package. Finally, we provide general guidelines and discuss the utility of performing a multivariable Mendelian randomization analysis for simultaneously assessing causality of multiple exposures. © 2021 Wiley Periodicals LLC. Basic Protocol: Performing a two-sample multivariable Mendelian randomization analysis using the 'MVMR' package in R and summarized genetic data Support Protocol 1: Installing the 'MVMR' R package Support Protocol 2: Obtaining instruments from the 'MRInstruments' R package.

Identification of novel drug targets for Alzheimer's disease by integrating genetics and proteomes from brain and blood

Genome-wide association studies (GWASs) have discovered numerous risk genes for Alzheimer's disease (AD), but how these genes confer AD risk is challenging to decipher. To efficiently transform genetic associations into drug targets for AD, we employed an integrative analytical pipeline using proteomes in the brain and blood by systematically applying proteome-wide association study (PWAS), Mendelian randomization (MR) and Bayesian colocalization. Collectively, we identified the brain protein abundance of 7 genes (ACE, ICA1L, TOM1L2, SNX32, EPHX2, CTSH, and RTFDC1) are causal in AD (P < 0.05/proteins identified for PWAS and MR; PPH4 >80% for Bayesian colocalization). The proteins encoded by these genes were mainly expressed on the surface of glutamatergic neurons and astrocytes. Of them, ACE with its protein abundance was also identified in significant association with AD on the blood-based studies and showed significance at the transcriptomic level. SNX32 was also found to be associated with AD at the blood transcriptomic level. Collectively, our current study results on genetic, proteomic, and transcriptomic approaches has identified compelling genes, which may provide important leads to design future functional studies and potential drug targets for AD.

Herpesvirus infections and Alzheimer's disease: a Mendelian randomization study

Background

Observational studies have suggested that herpesvirus infection increased the risk of Alzheimer’s disease (AD), but it is unclear whether the association is causal. The aim of the present study is to evaluate the causal relationship between four herpesvirus infections and AD.

Methods

We performed a two-sample Mendelian randomization analysis to investigate association of four active herpesvirus infections with AD using summary statistics from genome-wide association studies. The four herpesvirus infections (i.e., chickenpox, shingles, cold sores, mononucleosis) are caused by varicella-zoster virus, herpes simplex virus type 1, and Epstein-Barr virus (EBV), respectively. A large summary statistics data from International Genomics of Alzheimer’s Project was used in primary analysis, including 21,982 AD cases and 41,944 controls. Validation was further performed using family history of AD data from UK Biobank (27,696 cases of maternal AD, 14,338 cases of paternal AD and 272,244 controls).

Results

We found evidence of a significant association between mononucleosis (caused by EBV) and risk of AD after false discovery rates (FDR) correction (odds ratio [OR] = 1.634, 95% confidence interval [CI] = 1.092–2.446, P = 0.017, FDR-corrected P = 0.034). It has been verified in validation analysis that mononucleosis is also associated with family history of AD (OR [95% CI] = 1.392 [1.061, 1.826], P = 0.017). Genetically predicted shingles were associated with AD risk (OR [95% CI] = 0.867 [0.784, 0.958], P = 0.005, FDR-corrected P = 0.020), while genetically predicted chickenpox was suggestively associated with increased family history of AD (OR [95% CI] = 1.147 [1.007, 1.307], P = 0.039).

Conclusions

Our findings provided evidence supporting a positive relationship between mononucleosis and AD, indicating a causal link between EBV infection and AD. Further elucidations of this association and underlying mechanisms are likely to identify feasible interventions to promote AD prevention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-021-00905-5.

A Causal Web between Chronotype and Metabolic Health Traits

Observational and experimental evidence has linked chronotype to both psychological and cardiometabolic traits. Recent Mendelian randomization (MR) studies have investigated direct links between chronotype and several of these traits, often in isolation of outside potential mediating or moderating traits. We mined the EpiGraphDB MR database for calculated chronotype–trait associations (p-value < 5 × 10⁻⁸). We then re-analyzed those relevant to metabolic or mental health and investigated for statistical evidence of horizontal pleiotropy. Analyses passing multiple testing correction were then investigated for confounders, colliders, intermediates, and reverse intermediates using the EpiGraphDB database, creating multiple chronotype–trait interactions among each of the the traits studied. We revealed 10 significant chronotype–exposure associations (false discovery rate < 0.05) exposed to 111 potential previously known confounders, 52 intermediates, 18 reverse intermediates, and 31 colliders. Chronotype–lipid causal associations collided with treatment and diabetes effects; chronotype–bipolar associations were mediated by breast cancer; and chronotype–alcohol intake associations were impacted by confounders and intermediate variables including known zeitgebers and molecular traits. We have reported the influence of chronotype on several cardiometabolic and behavioural traits, and identified potential confounding variables not reported on in studies while discovering new associations to drugs and disease.

Investigation of glycaemic traits in psychiatric disorders using Mendelian randomisation revealed a causal relationship with anorexia nervosa

Data from observational studies have suggested an involvement of abnormal glycaemic regulation in the pathophysiology of psychiatric illness. This may be an attractive target for clinical intervention as glycaemia can be modulated by both lifestyle factors and pharmacological agents. However, observational studies are inherently confounded, and therefore, causal relationships cannot be reliably established. We employed genetic variants rigorously associated with three glycaemic traits (fasting glucose, fasting insulin, and glycated haemoglobin) as instrumental variables in a two-sample Mendelian randomisation analysis to investigate the causal effect of these measures on the risk for eight psychiatric disorders. A significant protective effect of a natural log transformed pmol/L increase in fasting insulin levels was observed for anorexia nervosa after the application of multiple testing correction (OR = 0.48 [95% CI: 0.33-0.71]-inverse-variance weighted estimate). There was no consistently strong evidence for a causal effect of glycaemic factors on the other seven psychiatric disorders considered. The relationship between fasting insulin and anorexia nervosa was supported by a suite of sensitivity analyses, with no statistical evidence of instrument heterogeneity or horizontal pleiotropy. Further investigation is required to explore the relationship between insulin levels and anorexia.

Shared genetic etiology between idiopathic pulmonary fibrosis and COVID-19 severity

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a complex lung disease, characterized by progressive lung scarring. Severe COVID-19 is associated with substantial pneumonitis and has a number of shared major risk factors with IPF. This study aimed to determine the genetic correlation between IPF and severe COVID-19 and assess a potential causal role of genetically increased risk of IPF on COVID-19 severity.

METHODS

The genetic correlation between IPF and COVID-19 severity was estimated with linkage disequilibrium (LD) score regression. We performed a Mendelian randomization (MR) study for IPF causality in COVID-19. Genetic variants associated with IPF susceptibility (P<5 × 10-8) in previous genome-wide association studies (GWAS) were used as instrumental variables (IVs). Effect estimates of those IVs on COVID-19 severity were gathered from the GWAS meta-analysis by the COVID-19 Host Genetics Initiative (4,336 cases & 623,902 controls).

FINDINGS

We detected a positive genetic correlation of IPF with COVID-19 severity (rg=0·31 [95% CI 0·04-0·57], P = 0·023). The MR estimates for severe COVID-19 did not reveal any genetic association (OR 1·05, [95% CI 0·92-1·20], P = 0·43). However, outlier analysis revealed that the IPF risk allele rs35705950 at MUC5B had a different effect compared with the other variants. When rs35705950 was excluded, MR results provided evidence that genetically increased risk of IPF has a causal effect on COVID-19 severity (OR 1·21, [95% CI 1·06-1·38], P = 4·24 × 10-3). Furthermore, the IPF risk-allele at MUC5B showed an apparent protective effect against COVID-19 hospitalization only in older adults (OR 0·86, [95% CI 0·73-1·00], P = 2·99 × 10-2) .

INTERPRETATION

The strongest genetic determinant of IPF, rs35705950 at MUC5B, seems to confer protection against COVID-19, whereas the combined effect of all other IPF risk loci seem to confer risk of COVID-19 severity. The observed effect of rs35705950 could either be due to protective effects of mucin over-production on the airways or a consequence of selection bias due to (1) a patient group that is heavily enriched for the rs35705950 T undertaking strict self-isolation and/or (2) due to survival bias of the rs35705950 non-IPF risk allele carriers. Due to the diverse impact of IPF causal variants on SARS-CoV-2 infection, with a possible selection bias as an explanation, further investigation is needed to address this apparent paradox between variance at MUC5B and other IPF genetic risk factors.

FUNDING

Novo Nordisk Foundation and Oak Foundation.

Genetically determined blood pressure, antihypertensive medications, and risk of Alzheimer's disease: a Mendelian randomization study

BACKGROUND

Observational studies suggest that the use of antihypertensive medications (AHMs) is associated with a reduced risk of Alzheimer's disease (AD); however, these findings may be biased by confounding and reverse causality. We aimed to explore the effects of blood pressure (BP) and lowering systolic BP (SBP) via the protein targets of different AHMs on AD through a two-sample Mendelian randomization (MR) approach.

METHODS

Genetic proxies from genome-wide association studies of BP traits and BP-lowering variants in genes encoding AHM targets were extracted. Estimates were calculated by inverse-variance weighted method as the main model. MR Egger regression and leave-one-out analysis were performed to identify potential violations.

RESULTS

There was limited evidence that genetically predicted SBP/diastolic BP level affected AD risk based on 400/398 single nucleotide polymorphisms (SNPs), respectively (all P > 0.05). Suitable genetic variants for β-blockers (1 SNP), angiotensin receptor blockers (1 SNP), calcium channel blockers (CCBs, 45 SNPs), and thiazide diuretics (5 SNPs) were identified. Genetic proxies for CCB [odds ratio (OR) = 0.959, 95% confidence interval (CI) = 0.941-0.977, P = 3.92 × 10-6] and overall use of AHMs (OR = 0.961, 95% CI = 0.944-0.978, P = 5.74 × 10-6, SNPs = 52) were associated with a lower risk of AD. No notable heterogeneity and directional pleiotropy were identified (all P > 0.05). Additional analyses partly support these results. No single SNP was driving the observed effects.

CONCLUSIONS

This MR analysis found evidence that genetically determined lowering BP was associated with a lower risk of AD and CCB was identified as a promising strategy for AD prevention.

Social Isolation, Social Interaction, and Alzheimer's Disease: A Mendelian Randomization Study

BACKGROUND

Social isolation and social interaction have been suggested to be associated with Alzheimer's disease. However, the causality cannot be unambiguously assessed as traditional epidemiological methods are easily subject to unmeasured confounders and potential bias.

OBJECTIVE

To examine bidirectional relationships between social isolation, social interaction, and Alzheimer's disease using Mendelian randomization method for assessing potential causal inference.

METHODS

This bidirectional two-sample Mendelian randomization study used independent genetic variants associated with social isolation and social interaction (n = 302,567-487,647), and Alzheimer's disease (n = 455,258). MR analyses were performed using the inverse-variance-weighted (IVW) as the main MR analytical method to estimate the causal effect. For sensitivity analyses, we applied weighted median, MR Egger to further assess the credibility of the causal effect.

RESULTS

Of the five types of social engagement examined in our study, only one showed evidence of an association with the risk of Alzheimer's disease. Attendance at a gym or sports club (IVW OR per SD change: 0.670; 95% CI: 0.463-0.970; p = 0.034) was inversely associated with the risk of Alzheimer's disease. We also found that AD may reduce the attendance at religious group (IVW OR per SD change: 1.017; 95% CI: 1.005-1.030; p = 0.004).

CONCLUSION

This study suggests that regular attendance at a gym or sports club is causally associated with reduced risk of Alzheimer's disease. Further studies are warranted to elucidate potential mechanisms.