Summary Data-based Mendelian Randomization Research Articles

Association of lipid-lowering drug targets with risk of cutaneous melanoma: a mendelian randomization study

BackgroundMelanoma proliferation is partly attributed to dysregulated lipid metabolism. The effectiveness of lipid-lowering drugs in combating cutaneous melanoma (CM) is a subject of ongoing debate in both in vitro and clinical studies.MethodThis study aims to evaluate the causal relationship between various lipid-lowering drug targets, namely 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR, targeted by statins), Proprotein convertase subtilisin/kexin type 9 (PCSK9, targeted by alirocumab and evolocumab), and Niemann-Pick C1-like 1 (NPC1L1, targeted by ezetimibe), and the outcomes of cutaneous melanoma. To mimic the effects of lipid-lowering drugs, we utilized two genetic tools: analysis of polymorphisms affecting the expression levels of drug target genes, and genetic variations linked to low-density lipoprotein cholesterol levels and drug target genes. These variations were sourced from genome-wide association studies (GWAS). We applied Summary-data-based Mendelian Randomization (SMR) and Inverse Variance Weighted Mendelian Randomization (IVW-MR) to gauge the effectiveness of these drugs.ResultsOur findings, with SMR results showing an odds ratio (OR) of 1.44 (95% CI: 1.08–1.92; P = 0.011) and IVW-MR results indicating an OR of 1.56 (95% CI: 1.10–2.23; P = 0.013), demonstrate a positive correlation between PCSK9 expression and increased risk of CM. However, no such correlations were observed in other analyses.ConclusionThe study concludes that PCSK9 plays a significant role in the development of CM, and its inhibition is linked to a reduced risk of the disease.

Identifying potential drug targets for varicose veins through integration of GWAS and eQTL summary data

BackgroundVaricose veins (VV) are a common chronic venous disease that is influenced by multiple factors. It affects the quality of life of patients and imposes a huge economic burden on the healthcare system. This study aimed to use integrated analysis methods, including Mendelian randomization analysis, to identify potential pathogenic genes and drug targets for VV treatment.MethodsThis study conducted Summary-data-based Mendelian Randomization (SMR) analysis and colocalization analysis on data collected from genome-wide association studies and cis-expression quantitative trait loci databases. Only genes with PP.H4 > 0.7 in colocalization were chosen from the significant SMR results. After the above analysis, we screened 12 genes and performed Mendelian Randomization (MR) analysis on them. After sensitivity analysis, we identified four genes with potential causal relationships with VV. Finally, we used transcriptome-wide association studies and The Drug-Gene Interaction Database data to identify and screen the remaining genes and identified four drug targets for the treatment of VV.ResultsWe identified four genes significantly associated with VV, namely, KRTAP5-AS1 [Odds ratio (OR) = 1.08, 95% Confidence interval (CI): 1.05–1.11, p = 1.42e-10] and PLEKHA5 (OR = 1.13, 95% CI: 1.06–1.20, p = 6.90e-5), CBWD1 (OR = 1.05, 95% CI: 1.01–1.11, p = 1.42e-2) and CRIM1 (OR = 0.87, 95% CI: 0.81–0.95, p = 3.67e-3). Increased expression of three genes, namely, KRTAP5-AS1, PLEKHA5, and CBWD1, was associated with increased risk of the disease, and increased expression of CRIM1 was associated with decreased risk of the disease. These four genes could be targeted for VV therapy.ConclusionWe identified four potential causal proteins for varicose veins with MR. A comprehensive analysis indicated that KRTAP5-AS1, PLEKHA5, CBWD1, and CRIM1 might be potential drug targets for varicose veins.

Genetic variations in anti-diabetic drug targets and COPD risk: evidence from mendelian randomization

BackgroundPrevious research has emphasized the potential benefits of anti-diabetic medications in inhibiting the exacerbation of Chronic Obstructive Pulmonary Disease (COPD), yet the role of anti-diabetic drugs on COPD risk remains uncertain.MethodsThis study employed a Mendelian randomization (MR) approach to evaluate the causal association of genetic variations related to six classes of anti-diabetic drug targets with COPD. The primary outcome for COPD was obtained from the Global Biobank Meta-analysis Initiative (GBMI) consortium, encompassing a meta-analysis of 12 cohorts with 81,568 cases and 1,310,798 controls. Summary-level data for HbA1c was derived from the UK Biobank, involving 344,182 individuals. Positive control analysis was conducted for Type 2 Diabetes Mellitus (T2DM) to validate the choice of instrumental variables. The study applied Summary-data-based MR (SMR) and two-sample MR for effect estimation and further adopted colocalization analysis to verify evidence of genetic variations.ResultsSMR analysis revealed that elevated KCNJ11 gene expression levels in blood correlated with reduced COPD risk (OR = 0.87, 95% CI = 0.79–0.95; p = 0.002), whereas an increase in DPP4 expression corresponded with an increased COPD incidence (OR = 1.18, 95% CI = 1.03–1.35; p = 0.022). Additionally, the primary method within MR analysis demonstrated a positive correlation between PPARG-mediated HbA1c and both FEV1 (OR = 1.07, 95% CI = 1.02–1.13; P = 0.013) and FEV1/FVC (OR = 1.08, 95% CI = 1.01–1.14; P = 0.007), and a negative association between SLC5A2-mediated HbA1c and FEV1/FVC (OR = 0.86, 95% CI = 0.74–1.00; P = 0.045). No colocalization evidence with outcome phenotypes was detected (all PP.H4 < 0.7).ConclusionThis study provides suggestive evidence for anti-diabetic medications' role in improving COPD and lung function. Further updated MR analyses are warranted in the future, following the acquisition of more extensive and comprehensive data, to validate our results.

Genetically determined blood pressure, antihypertensive drug classes, and frailty: A Mendelian randomization study.

Observational studies have suggested that the use of antihypertensive drugs was associated with the risk of frailty; however, these findings may be biased by confounding and reverse causality. This study aimed to explore the effect of genetically predicted lifelong lowering blood pressure (BP) through different antihypertensive medications on frailty. One-sample Mendelian randomization (MR) and summary data-based MR (SMR) were applied. We utilized two kinds of genetic instruments to proxy the antihypertensive medications, including genetic variants within or nearby drugs target genes associated with systolic/diastolic BP, and expression level of the corresponding gene. Among 298,618 UK Biobank participants, one-sample MR analysis observed that genetically proxied BB use (relative risk ratios, 0.76; 95% CI, 0.65-0.90; p = 0.001) and CCB use (0.83; 0.72-0.95; p = 0.007), equivalent to a 10-mm Hg reduction in systolic BP, was significantly associated with lower risk of pre-frailty. In addition, although not statistically significant, the effect directions of systolic BP through ACEi variants (0.72; 0.39-1.33; p = 0.296) or thiazides variants (0.74; 0.53-1.03; p = 0.072) on pre-frailty were also protective. Similar results were obtained in analyses for diastolic BP. SMR of expression in artery showed that decreased expression level of KCNH2, a target gene of BBs, was associated with lower frailty index (beta -0.02, p = 2.87 × 10-4). This MR analysis found evidence that the use of BBs and CCBs was potentially associated with reduced frailty risk in the general population, and identified KCNH2 as a promising target for further clinical trials to prevent manifestations of frailty.

Genetic regulation of human brain proteome reveals proteins implicated in psychiatric disorders.

Psychiatric disorders are highly heritable yet polygenic, potentially involving hundreds of risk genes. Genome-wide association studies have identified hundreds of genomic susceptibility loci with susceptibility to psychiatric disorders; however, the contribution of these loci to the underlying psychopathology and etiology remains elusive. Here we generated deep human brain proteomics data by quantifying 11,608 proteins across 268 subjects using 11-plex tandem mass tag coupled with two-dimensional liquid chromatography-tandem mass spectrometry. Our analysis revealed 788 cis-acting protein quantitative trait loci associated with the expression of 883 proteins at a genome-wide false discovery rate <5%. In contrast to expression at the transcript level and complex diseases that are found to be mainly influenced by noncoding variants, we found protein expression level tends to be regulated by non-synonymous variants. We also provided evidence of 76 shared regulatory signals between gene expression and protein abundance. Mediation analysis revealed that for most (88%) of the colocalized genes, the expression levels of their corresponding proteins are regulated by cis-pQTLs via gene transcription. Using summary data-based Mendelian randomization analysis, we identified 4 proteins and 19 genes that are causally associated with schizophrenia. We further integrated multiple omics data with network analysis to prioritize candidate genes for schizophrenia risk loci. Collectively, our findings underscore the potential of proteome-wide linkage analysis in gaining mechanistic insights into the pathogenesis of psychiatric disorders.

An integrative framework to prioritize genes in more than 500 loci associated with body mass index.

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

​Comprehensive mendelian randomization analysis of plasma proteomics to identify new therapeutic targets for the treatment of coronary heart disease and myocardial infarction.

Ischemic heart disease is one of the leading causes of mortality worldwide, and thus calls for development of more effective therapeutic strategies. This study aimed to identify potential therapeutic targets for coronary heart disease (CHD) and myocardial infarction (MI) by investigating the causal relationship between plasma proteins and these conditions. A two-sample Mendelian randomization (MR) study was performed to evaluate more than 1600 plasma proteins for their causal associations with CHD and MI. The MR findings were further confirmed through Bayesian colocalization, Summary-data-based Mendelian Randomization (SMR), and Transcriptome-Wide Association Studies (TWAS) analyses. Further analyses, including enrichment analysis, single-cell analysis, MR analysis of cardiovascular risk factors, phenome-wide Mendelian Randomization (Phe-MR), and protein-protein interaction (PPI) network construction were conducted to verify the roles of selected causal proteins. Thirteen proteins were causally associated with CHD, seven of which were also causal for MI. Among them, FES and PCSK9 were causal proteins for both diseases as determined by several analytical methods. PCSK9 was a risk factor of CHD (OR = 1.25, 95% CI: 1.13-1.38, P = 7.47E-06) and MI (OR = 1.36, 95% CI: 1.21-1.54, P = 2.30E-07), whereas FES was protective against CHD (OR = 0.68, 95% CI: 0.59-0.79, P = 6.40E-07) and MI (OR = 0.65, 95% CI: 0.54-0.77, P = 5.38E-07). Further validation through enrichment and single-cell analysis confirmed the causal effects of these proteins. Moreover, MR analysis of cardiovascular risk factors, Phe-MR, and PPI network provided insights into the potential drug development based on the proteins. This study investigated the causal pathways associated with CHD and MI, highlighting the protective and risk roles of FES and PCSK9, respectively. FES. Specifically, the results showed that these proteins are promising therapeutic targets for future drug development.

Integration of multi-omics summary data reveals the role of N6-methyladenosine in neuropsychiatric disorders.

N6-methyladenosine (m6A) methylation regulates gene expression/protein by influencing numerous aspects of mRNA metabolism and contributes to neuropsychiatric diseases. Here, we integrated multi-omics data and genome-wide association study summary data of schizophrenia (SCZ), bipolar disorder (BP), attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD) to reveal the role of m6A in neuropsychiatric disorders by using transcriptome-wide association study (TWAS) tool and Summary-data-based Mendelian randomization (SMR). Our investigation identified 86 m6A sites associated with seven neuropsychiatric diseases and then revealed 7881 associations between m6A sites and gene expressions. Based on these results, we discovered 916 significant m6A-gene associations involving 82 disease-related m6A sites and 606 genes. Further integrating the 58 disease-related genes from TWAS and SMR analysis, we obtained 61, 8, 7, 3, and 2 associations linking m6A-disease, m6A-gene, and gene-disease for SCZ, BP, AD, MDD, and PD separately. Functional analysis showed the m6A mapped genes were enriched in "response to stimulus" pathway. In addition, we also analyzed the effect of gene expression on m6A and the post-transcription effect of m6A on protein. Our study provided new insights into the genetic component of m6A in neuropsychiatric disorders and unveiled potential pathogenic mechanisms where m6A exerts influences on disease through gene expression/protein regulation.

XWAS analysis in neuropsychiatric disorders by integrating multi-molecular phenotype quantitative trait loci and GWAS summary data

BackgroundIntegrating quantitative trait loci (QTL) data related to molecular phenotypes with genome-wide association study (GWAS) data is an important post-GWAS strategic approach employed to identify disease-associated molecular features. Various types of molecular phenotypes have been investigated in neuropsychiatric disorders. However, these findings pertaining to distinct molecular features are often independent of each other, posing challenges for having an overview of the mapped genes.MethodsIn this study, we comprehensively summarized published analyses focusing on four types of risk-related molecular features (gene expression, splicing transcriptome, protein abundance, and DNA methylation) across five common neuropsychiatric disorders. Subsequently, we conducted supplementary analyses with the latest GWAS dataset and corresponding deficient molecular phenotypes using Functional Summary-based Imputation (FUSION) and summary data-based Mendelian randomization (SMR). Based on the curated and supplemented results, novel reliable genes and their functions were explored.ResultsOur findings revealed that eQTL exhibited superior ability in prioritizing risk genes compared to the other QTL, followed by sQTL. Approximately half of the genes associated with splicing transcriptome, protein abundance, and DNA methylation were successfully replicated by eQTL-associated genes across all five disorders. Furthermore, we identified 436 novel reliable genes, which enriched in pathways related with neurotransmitter transportation such as synaptic, dendrite, vesicles, axon along with correlations with other neuropsychiatric disorders. Finally, we identified ten multiple molecular involved regulation patterns (MMRP), which may provide valuable insights into understanding the contribution of molecular regulation network targeting these disease-associated genes.ConclusionsThe analyses prioritized novel and reliable gene sets related with five molecular features based on published and supplementary results for five common neuropsychiatric disorders, which were missed in the original GWAS analysis. Besides, the involved MMRP behind these genes could be given priority for further investigation to elucidate the pathogenic molecular mechanisms underlying neuropsychiatric disorders in future studies.

Prioritization of therapeutic targets for cancers using integrative multi-omics analysis

BackgroundThe integration of transcriptomic, proteomic, druggable genetic and metabolomic association studies facilitated a comprehensive investigation of molecular features and shared pathways for cancers’ development and progression.MethodsComprehensive approaches consisting of transcriptome-wide association studies (TWAS), proteome-wide association studies (PWAS), summary-data-based Mendelian randomization (SMR) and MR were performed to identify genes significantly associated with cancers. The results identified in above analyzes were subsequently involved in phenotype scanning and enrichment analyzes to explore the possible health effects and shared pathways. Additionally, we also conducted MR analysis to investigate metabolic pathways related to cancers.ResultsTotally 24 genes (18 transcriptomic, 1 proteomic and 5 druggable genetic) showed significant associations with cancers risk. All genes identified in multiple methods were mainly enriched in nuclear factor erythroid 2-related factor 2 (NRF2) pathway. Additionally, biosynthesis of ubiquinol and urate were found to play an important role in gastrointestinal tumors.ConclusionsA set of putatively causal genes and pathways relevant to cancers were identified in this study, shedding light on the shared biological processes for tumorigenesis and providing compelling genetic evidence to prioritize anti-cancer drugs development.

Association of Cadherin-Related Family Member 1 with Traumatic Brain Injury

The cadherin family plays a pivotal role in orchestrating synapse formation in the central nervous system. Cadherin-related family member 1 (CDHR1) is a photoreceptor-specific calmodulin belonging to the expansive cadherin superfamily. However, its role in traumatic brain injury (TBI) remains largely unknown. CDHR1 expression across various brain tissue sites was analyzed using the GSE104687 dataset. Employing a summary-data-based Mendelian Randomization (SMR) approach, integrated analyses were performed by amalgamating genome-wide association study abstracts from TBI with public data on expressed quantitative trait loci and DNA methylation QTL from both blood and diverse brain tissues. CDHR1 expression and localization in different brain tissues were meticulously delineated using western blotting, immunohistochemistry, and enzyme-linked immunosorbent assay. CDHR1 expression was consistently elevated in the TBI group compared to that in the sham group across multiple tissues. The inflammatory response emerged as a crucial biological mechanism, and pro-inflammatory and anti-inflammatory factors were not expressed in either group. Integrated SMR analyses encompassing both blood and brain tissues substantiated the heightened CDHR1 expression profiles, with methylation modifications emerging as potential contributing factors for increased TBI risk. This was corroborated by western blotting and immunohistochemistry, confirming augmented CDHR1 expression following TBI. This multi-omics-based genetic association study highlights the elevated TBI risk associated with CDHR1 expression coupled with putative methylation modifications. These findings provide compelling evidence for future targeted investigations and offer promising avenues for developing interventional therapies for TBI.Graphical Cadherin-related family member 1 as a potential risk factor for traumatic brain injury.

Dissecting the shared genetic landscape of anxiety, depression, and schizophrenia

BackgroundNumerous studies highlight the genetic underpinnings of mental disorders comorbidity, particularly in anxiety, depression, and schizophrenia. However, their shared genetic loci are not well understood. Our study employs Mendelian randomization (MR) and colocalization analyses, alongside multi-omics data, to uncover potential genetic targets for these conditions, thereby informing therapeutic and drug development strategies.MethodsWe utilized the Consortium for Linkage Disequilibrium Score Regression (LDSC) and Mendelian Randomization (MR) analysis to investigate genetic correlations among anxiety, depression, and schizophrenia. Utilizing GTEx V8 eQTL and deCODE Genetics pQTL data, we performed a three-step summary-data-based Mendelian randomization (SMR) and protein–protein interaction analysis. This helped assess causal and comorbid loci for these disorders and determine if identified loci share coincidental variations with psychiatric diseases. Additionally, phenome-wide association studies, drug prediction, and molecular docking validated potential drug targets.ResultsWe found genetic correlations between anxiety, depression, and schizophrenia, and under a meta-analysis of MR from multiple databases, the causal relationships among these disorders are supported. Based on this, three-step SMR and colocalization analyses identified ITIH3 and CCS as being related to the risk of developing depression, while CTSS and DNPH1 are related to the onset of schizophrenia. BTN3A1, PSMB4, and TIMP4 were identified as comorbidity loci for both disorders. Molecules that could not be determined through colocalization analysis were also presented. Drug prediction and molecular docking showed that some drugs and proteins have good binding affinity and available structural data.ConclusionsOur study indicates genetic correlations and shared risk loci between anxiety, depression, and schizophrenia. These findings offer insights into the underlying mechanisms of their comorbidities and aid in drug development.

Mendelian randomization analysis identifies druggable genes and drugs repurposing for chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is a prevalent condition that significantly impacts public health. Unfortunately, there are few effective treatment options available. Mendelian randomization (MR) has been utilized to repurpose existing drugs and identify new therapeutic targets. The objective of this study is to identify novel therapeutic targets for COPD. Cis-expression quantitative trait loci (cis-eQTL) were extracted for 4,317 identified druggable genes from genomics and proteomics data of whole blood (eQTLGen) and lung tissue (GTEx Consortium). Genome-wide association studies (GWAS) data for doctor-diagnosed COPD, spirometry-defined COPD (Forced Expiratory Volume in one second [FEV1]/Forced Vital Capacity [FVC] <0.7), and FEV1 were obtained from the cohort of FinnGen, UK Biobank and SpiroMeta consortium. We employed Summary-data-based Mendelian Randomization (SMR), HEIDI test, and colocalization analysis to assess the causal effects of druggable gene expression on COPD and lung function. The reliability of these druggable genes was confirmed by eQTL two-sample MR and protein quantitative trait loci (pQTL) SMR, respectively. The potential effects of druggable genes were assessed through the phenome-wide association study (PheWAS). Information on drug repurposing for COPD was collected from multiple databases. A total of 31 potential druggable genes associated with doctor-diagnosed COPD, spirometry-defined COPD, and FEV1 were identified through SMR, HEIDI test, and colocalization analysis. Among them, 22 genes (e.g., MMP15, PSMA4, ERBB3, and LMCD1) were further confirmed by eQTL two-sample MR and protein SMR analyses. Gene-level PheWAS revealed that ERBB3 expression might reduce inflammation, while GP9 and MRC2 were associated with other traits. The drugs Montelukast (targeting the MMP15 gene) and MARIZOMIB (targeting the PSMA4 gene) may reduce the risk of spirometry-defined COPD. Additionally, an existing small molecule inhibitor of the APH1A gene has the potential to increase FEV1. Our findings identified 22 potential drug targets for COPD and lung function. Prioritizing clinical trials that target these identified druggable genes with existing drugs or novel medications will be beneficial for the development of COPD treatments.

Identification of prospective aging drug targets via Mendelian randomization analysis.

Aging represents a multifaceted process culminating in the deterioration of biological functions. Despite the introduction of numerous anti-aging strategies, their therapeutic outcomes have often been less than optimal. Consequently, discovering new targets to mitigate aging effects is of critical importance. We applied Mendelian randomization (MR) to identify potential pharmacological targets against aging, drawing upon summary statistics from both the Decode and FinnGen cohorts, with further validation in an additional cohort. To address potential reverse causality, bidirectional MR analysis with Steiger filtering was utilized. Additionally, Bayesian co-localization and phenotype scanning were implemented to investigate previous associations between genetic variants and traits. Summary-data-based Mendelian randomization (SMR) analysis was conducted to assess the impact of genetic variants on aging via their effects on protein expression. Additionally, mediation analysis was orchestrated to uncover potential intermediaries in these associations. Finally, we probed the systemic implications of drug-target protein expression across diverse indications by MR-PheWas analysis. Utilizing a Bonferroni-corrected threshold, our MR examination identified 10 protein-aging associations. Within this cohort of proteins, MST1, LCT, GMPR2, PSMB4, ECM1, EFEMP1, and ISLR2 appear to exacerbate aging risks, while MAX, B3GNT8, and USP8 may exert protective influences. None of these proteins displayed reverse causality except EFEMP1. Bayesian co-localization inferred shared variants between aging and proteins such as B3GNT8 (rs11670143), ECM1 (rs61819393), and others listed. Mediator analysis pinpointed 1,5-anhydroglucitol as a partial intermediary in the influence LCT exhibits on telomere length. Circulating proteins play a pivotal role in influencing the aging process, making them promising candidates for therapeutic intervention. The implications of these proteins in aging warrant further investigation in future clinical research.

Potential Benefits of Statin Therapy in Reducing Osteoarthritis Risk: A Mendelian Randomization Study.

The purpose of this study was to determine the causal effect of statins on osteoarthritis (OA) risk using Mendelian randomization (MR). Single nucleotide polymorphism-based genome-wide association analyses of statins were collected from the UK Biobank and FinnGen dataset, and OA data were collected from the UK Biobank and Arthritis Research UK Osteoarthritis Genetics (arcOGEN) study. Two-sample MR analyses were performed using the inverse-variance weighted (IVW) technique. MR-Egger, weighted median, and weighted mode served as supplementary analyses. MR-Egger regression, Cochran's Q test, and Mendelian Randomization Pleiotropy Residual Sum and Outlier analysis were performed as sensitivity analyses. Hydroxymethylglutaryl-coenzyme A reductase (HMGCR) expression and OA risk were evaluated using summary data-based MR (SMR). MR analyses consistently supported a causal connection between statin use and OA risk. A causal effect was observed for atorvastatin (IVW: β = -2.989, P = 0.003) and rosuvastatin (IVW: β = -14.141, P = 0.006) treatment on hip OA. Meta-analysis showed the association between atorvastatin and knee OA was statistically significant (odds ratio 0.15; P = 0.004). Simvastatin use exhibited a protective effect against knee (IVW: β = -1.056, P = 0.004) and hip OA (IVW: β = -1.405, P = 0.001). Statin medication showed a protective effect on hip OA (IVW: β = -0.054, P = 0.013). HMGCR correlated significantly with a reduced risk of knee OA (β = -0.193, PSMR = 0.017), rather than hip OA (β = 0.067, PSMR = 0.502), which suggested that statins' protective effect on OA may not be related to its lipid-lowering effect. This MR study provides compelling evidence that statin treatment may be a protective factor for OA. Further research is required to clarify its underlying mechanism.

Causal effect of gut microbiota on pancreatic cancer: A Mendelian randomization and colocalization study.

The causal relationship between gut microbiota (GM) and pancreatic cancer (PC) remains unclear. This study aimed to investigate the potential genes underlying this mechanism. GM Genome-wide association study (GWAS) summary data were from the MiBioGen consortium. PC GWAS data were from the National Human Genome Research Institute-European Bioinformatics Institute (NHGRI-EBI) GWAS Catalogue. To detect the causal relationship between GM and PC, we implemented three complementary Mendelian randomization (MR) methods: Inverse Variance Weighting (IVW), MR-Egger and Weighted Median, followed by sensitivity analyses. Furthermore, we integrated GM GWAS data with blood cis-expression quantitative trait loci (eQTLs) and blood cis-DNA methylation QTL (mQTLs) using Summary data-based Mendelian Randomization (SMR) methods. This integration aimed to prioritize potential GM-affecting genes through SMR analysis of two molecular traits. PC cis-eQTLs and cis-mQTLs were summarized from The Cancer Genome Atlas (TCGA) data. Through colocalization analysis of GM cis-QTLs and PC cis-QTLs data, we identified common genes that influence both GM and PC. Our study found a causal association between GM and PC, including four protective and five risk-associated GM [Inverse Variance Weighted (IVW), p < 0.05]. No significant heterogeneity of instrumental variables (IVs) or horizontal pleiotropy was found. The gene SVBP was identified as a GM-affecting gene using SMR analysis of two molecular traits (FDR<0.05, P_HEIDI>0.05). Additionally, two genes, MCM6 and RPS26, were implicated in the interaction between GM and PC based on colocalization analysis (PPH4>0.5). In summary, this study provides evidence for future research aimed at developing suitable therapeutic interventions and disease prevention.

Statins as a risk factor for diabetic retinopathy: a Mendelian randomization and cross-sectional observational study

BackgroundDiabetic retinopathy (DR) is the foremost cause of vision loss among the global working-age population, and statins are among the most frequently prescribed drugs for lipid management in patients with DR. The exact relationship between statins and DR has not been determined. This study sought to validate the causal association between statins usage and diabetic retinopathy.MethodsThe summary-data-based Mendelian randomization (SMR) method and inverse-variance-weighted Mendelian randomization (IVW-MR) were used to identify the causal relationship between statins and DR via the use of expression quantitative trait loci (eQTL) data for 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) (31,684 blood samples), low density lipoprotein cholesterol-related GWAS data (sample size: 440,546), and DR-related GWAS data (14,584 cases and 176,010 controls). Additionally, a cross-sectional observational study based on the data from the National Health and Nutrition Examination Survey (NHANES) was conducted to supplement the association between DR and statins (sample size: 106,911). The odds ratios (ORs) with corresponding 95% confidence intervals (CIs) was employed to evaluate the results.ResultsBased on the results of the MR analysis, HMGCR inhibitors were causally connected with a noticeably greater incidence of DR (IVW: OR = 0.54, 95% CI [0.42, 0.69], p = 0.000002; SMR: OR = 0.66, 95% CI [0.52, 0.84], p = 0.00073). Subgroup analysis revealed that the results were not affected by the severity of DR. The sensitivity analysis revealed the stability and reliability of the MR analysis results. The results from the cross-sectional study based on NHANES also support the association between not taking statins and a decreased risk of DR (OR = 0.54, 95% CI [0.37, 0.79], p = 0.001).ConclusionsThis study revealed that a significant increase in DR risk was causally related to statins use, providing novel insights into the role of statins in DR. However, further investigations are needed to verify these findings.

Genetic overlap and causality between COVID-19 and multi-site chronic pain: the importance of immunity.

The existence of chronic pain increases susceptibility to virus and is now widely acknowledged as a prominent feature recognized as a major manifestation of long-term coronavirus disease 2019 (COVID-19) infection. Given the ongoing COVID-19 pandemic, it is imperative to explore the genetic associations between chronic pain and predisposition to COVID-19. We conducted genetic analysis at the single nucleotide polymorphism (SNP), gene, and molecular levels using summary statistics of genome-wide association study (GWAS) and analyzed the drug targets by summary data-based Mendelian randomization analysis (SMR) to alleviate the multi-site chronic pain in COVID-19. Additionally, we performed a latent causal variable (LCV) method to investigate the causal relationship between chronic pain and susceptibility to COVID-19. The cross-trait meta-analysis identified 19 significant SNPs shared between COVID-19 and chronic pain. Coloc analysis indicated that the posterior probability of association (PPH4) for three loci was above 70% in both critical COVID-19 and COVID-19, with the corresponding top three SNPs being rs13135092, rs7588831, and rs13135092. A total of 482 significant overlapped genes were detected from MAGMA and CPASSOC results. Additionally, the gene ANAPC4 was identified as a potential drug target for treating chronic pain (P=7.66E-05) in COVID-19 (P=8.23E-03). Tissue enrichment analysis highlighted that the amygdala (P=7.81E-04) and prefrontal cortex (P=8.19E-05) as pivotal in regulating chronic pain of critical COVID-19. KEGG pathway enrichment further revealed the enrichment of pleiotropic genes in both COVID-19 (P=3.20E-03,Padjust=4.77E-02,hsa05171) and neurotrophic pathways (P=9.03E-04,Padjust =2.55E-02,hsa04621). Finally, the latent causal variable (LCV) model was applied to find the genetic component of critical COVID-19 was causal for multi-site chronic pain (P=0.015), with a genetic causality proportion (GCP) of was 0.60. In this study, we identified several functional genes and underscored the pivotal role of the inflammatory system in the correlation between the paired traits. Notably, heat shock proteins emerged as potential objective biomarkers for chronic pain symptoms in individuals with COVID-19. Additionally, the ubiquitin system might play a role in mediating the impact of COVID-19 on chronic pain. These findings contribute to a more comprehensive understanding of the pleiotropy between COVID-19 and chronic pain, offering insights for therapeutic trials.

Proteome-Wide Mendelian Randomization and Colocalization Analysis Identify Therapeutic Targets for Knee and Hip Osteoarthritis.

Osteoarthritis (OA) is a common degenerative disease. Although some biomarkers and drug targets of OA have been discovered and employed, limitations and challenges still exist in the targeted therapy of OA. Mendelian randomization (MR) analysis has been regarded as a reliable analytic method to identify effective therapeutic targets. Thus, we aimed to identify novel therapeutic targets for OA and investigate their potential side effects based on MR analysis. In this study, two-sample MR, colocalization analysis, summary-data-based Mendelian randomization (SMR) and Mendelian randomization phenome-wide association study (MR-PheWAS) were conducted. We firstly analyzed data from 4907 plasma proteins to identify potential therapeutic targets associated with OA. In addition, blood expression quantitative trait loci (eQTLs) data sources were used to perform additional validation. A protein-protein interaction (PPI) network was also constructed to delve into the interactions among identified proteins. Then, MR-PheWASs were utilized to assess the potential side effects of core therapeutic targets. After MR analysis and FDR correction, we identified twelve proteins as potential therapeutic targets for knee OA or hip OA. Colocalization analysis and additional validation supported our findings, and PPI networks revealed the interactions among identified proteins. Finally, we identified MAPK3 (OR = 0.855, 95% CI: 0.791-0.923, p = 6.88 × 10-5) and GZMK (OR = 1.278, 95% CI: 1.131-1.444, p = 8.58 × 10-5) as the core therapeutic targets for knee OA, and ITIH1 (OR = 0.847, 95% CI: 0.784-0.915, p = 2.44 × 10-5) for hip OA. A further MR phenome-wide association study revealed the potential side effects of treatments targeting MAPK3, GZMK, and ITIH1. This comprehensive study indicates twelve plasma proteins with potential roles in knee and hip OA as therapeutic targets. This advancement holds promise for the progression of OA drug development, and paves the way for more efficacious treatments of OA.

Integrated analysis reveals FLI1 regulates the tumor immune microenvironment via its cell-type-specific expression and transcriptional regulation of distinct target genes of immune cells in breast cancer

BackgroundImmunotherapy is a practical therapeutic approach in breast cancer (BRCA), and the role of FLI1 in immune regulation has gradually been unveiled. However, the specific role of FLI1 in BRCA was conflicted; thus, additional convincing evidence is needed.MethodsWe explored the upstream regulation of FLI1 expression via summary data-based Mendelian randomization (SMR) analysis and ncRNA network construction centering on FLI1 using BRCA genome-wide association study (GWAS) summary data with expression quantitative trait loci (eQTLs) and DNA methylation quantitative trait loci (mQTLs) from the blood and a series of in silico analyses, respectively. We illuminated the downstream function of FLI1 in immune regulation by integrating a series of analyses of single-cell RNA sequence data (scRNA-seq).ResultsWe verified a causal pathway from FLI1 methylation to FLI1 gene expression to BRCA onset and demonstrated that FLI1 was downregulated in BRCA. FLI1, a transcription factor, served as myeloid and T cells’ communication regulator by targeting immune-related ligands and receptor transcription in BRCA tissues. We constructed a ceRNA network centering on FLI1 that consisted of three LncRNAs (CKMT2-AS1, PSMA3-AS1, and DIO3OS) and a miRNA (hsa-miR-324-5p), and the expression of FLI1 was positively related to a series of immune-related markers, including immune cell infiltration, biomarkers of immune cells, and immune checkpoints.ConclusionLow-methylation-induced or ncRNA-mediated downregulation of FLI1 is associated with poor prognosis, and FLI1 might regulate the tumor immune microenvironment via a cell-type-specific target genes manner in BRCA.