SEARCH · Search PubMed

Results for “Druggable gene”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

Decoding Primary Open-Angle Glaucoma: A Multi-Omics Approach to Identify Druggable Effector Genes.

PURPOSE: Genomewide association studies (GWAS) have identified numerous primary open angle glaucoma (POAG) risk loci, yet most reside in non-coding regions with unclear function. Mapping these loci to effector genes can elucidate disease mechanisms, identify functionally conserved variants, improve cross-ancestry risk prediction by reducing population-specific noise, and uncover shared therapeutic targets. METHODS: Here, we integrate European POAG GWAS with six types of multi-omics molecular Quantitative Trait Locis (xQTLs) using multi-trait colocalization to identify candidate effector variants and evaluate their cross-population relevance using genetic risk score (GRS) analysis, and their therapeutic potential through drug target prioritization. RESULTS: We identified 25 POAG effector variants colocalized with at least one xQTLs. In non-European populations, effector variants showed stronger effect size correlations with Europeans than non-colocalized variants (Pearson r2 = African 0.85 vs. 0.71; East Asian 0.81 vs. 0.69; and Latin American 0.91 vs. 0.75). Effector variants also had smaller allele frequency variations across populations (average interquartile range [IQR] = 0.15 vs. 0.20). The genetic risk score based on effector variants performed comparably to the genome-wide significant single-nucleotide polymorphism (SNP)-based GRS in non-European populations. Drug prioritization identified zinc, copper, sunitinib, probucol, and astemizole as potential common therapeutic agents for POAG and its subtypes. CONCLUSIONS: Our findings offer deeper insight into the molecular mechanisms underlying glaucoma and effector variants for developing more robust GRS models and broadly effective therapeutic strategies for POAG.

Humans↗

PubMed2025

Reproducible autosomal gene expression changes with loss of typical X and Y complement across tumor types.

Although there are known sex differences in cancer incidence, severity, and treatment, the sex chromosomes are typically excluded from genomic analyses because of the unique technical challenges associated with assessing their copy number, sequence variation, and expression. Here we assess sex chromosome complement in three widely-used human genomics datasets from normal (non-cancerous) tissues, primary tumors, and cancer cell lines and study the effects on genome-wide gene expression. Expected sex chromosome complements based on reported patient sex were observed in non-cancerous tissues, but about half of tumors and cancer cell lines showed loss of typical sex chromosome gene expression across tissue types with three categories: loss of chromosome Y (LOY), loss of chromosome X (LOX) and reactivation of the inactive X chromosome (XaXa). Genes consistently differentially expressed in tumors with loss of chromosome X, loss of chromosome Y, or loss of X chromosome inactivation are associated with the hallmarks of cancer and include both sex-linked and autosomal genes from nearly all chromosomes, druggable genes, and genes with molecular functions relevant to cancer signaling, such as kinase activity. Strikingly, tumors that are X0, including tumors from female patients that have lost an X chromosome and tumors from male patients that have lost a Y chromosome, cluster together by gene expression profile. Patients with tumors that have LOX or LOY had poorer survival outcomes compared to those with tumors that had maintained their sex chromosome complement. Further, LOX and LOY eliminates nearly all of the differential gene expression between tumors from different patient sexes, affecting sex chromosomal and autosomal gene expression. Going forward, considering patient sex as well as the entire genome, including assessment of the sex chromosome complement, will provide additional insights into personalized tumor etiology, progression, treatment, and patient outcome.

Journal Article↗

PubMed2026

Systematic Identification of Therapeutic Targets and Repurposed Drugs for Stroke: From Genome Causal Analysis to Multilevel Validation.

BACKGROUND: Stroke is a severe cerebrovascular disease characterized by narrow time windows and complications. This study aimed to identify novel drug targets and repurposed drugs for stroke. METHODS: This study used expression quantitative trait loci data from druggable genes in brain and blood as instrumental variables. Mendelian randomization, colocalization, and phenome-wide Mendelian randomization were applied to evaluate causal relationships and potential side effects, with stroke and ischemic stroke as primary outcomes. Preclinical validation used oxygen-glucose deprivation/reperfusion and middle cerebral artery occlusion/reperfusion models. Pharmacological and behavioral assessments evaluated the therapeutic potential of candidate targets and drugs. Additionally, proteomic sequencing was performed following GGCX (γ-glutamyl carboxylase) overexpression to explore its biological functions. RESULTS: Elevated GGCX expression in brain and blood was potentially causally associated with reduced risk of stroke and ischemic stroke, supported by colocalization evidence, although potential cardiovascular risks could not be excluded. Drug repositioning identified ifenprodil as a candidate agent that reduced infarction volume, improved motor and cognitive functions, and reversed GGCX downregulation in mice. Ifenprodil treatment and GGCX overexpression alleviated oxygen-glucose deprivation/reperfusion-induced injury and upregulated GGCX expression. Mechanistically, GGCX conferred neuroprotection by regulating protein homeostasis, suppressing inflammation, promoting metabolic recovery, and modulating nuclear transcriptional regulation. CONCLUSIONS: This study established a potential causal link between GGCX and stroke risk, particularly ischemic stroke. GGCX represents a promising therapeutic target for ischemic stroke. Targeted GGCX expression upregulation and drug repurposing, particularly ifenprodil, may offer novel therapeutic avenues. Further validation is warranted to assess clinical efficacy and safety.

Animals↗

PubMed2026

Multi-omics integration and colocalization analyses prioritize candidate molecular loci associated with hypothermia.

BACKGROUND: Hypothermia is a life-threatening condition lacking specific pharmacological treatments. This study aimed to prioritize genetically supported molecular loci associated with hypothermia and to explore their pharmacological tractability using multi-omics data. METHODS: Initially, 2532 druggable genes were curated from the Drug-Gene Interaction Database and established literature. These were cross-referenced with cis-eQTL and cis-pQTL datasets, encompassing 870,655 and 114,281 SNPs for blood, respectively, alongside 2379 shared SNPs across adipose, skeletal muscle, and heart tissues. Matched instrumental variables were integrated with hypothermia GWAS summary statistics for two-sample Mendelian randomization (MR) and Bayesian colocalization. Transcriptomic differential expression analysis (DEA) was subsequently conducted as an exploratory analysis of cold-exposure-associated expression changes. Database-derived compound annotations were systematically re-evaluated according to target specificity, established pharmacological mechanism, and concordance with the direction of the MR estimates. RESULTS: Among 671 gene-level MR tests, 36 genes reached nominal significance, whereas only ABCC8 remained significant after FDR correction. Colocalization was evaluable for 8 of these 36 genes, and 4 loci (COL18A1, SLC1A7, ADIPOQ, and MERTK) met the prespecified PP.H4>0.90 threshold. The remaining 28 loci were not evaluable because sufficient overlapping regional variants were unavailable after harmonization. Transcriptomic analysis identified altered expression of SLC1A3 and SLCO4A1 under cold exposure, although these findings did not directly validate the colocalization-supported loci. Re-evaluation of database-derived compound annotations did not identify any direct, selective, and directionally concordant drug-repurposing candidate for hypothermia. CONCLUSIONS: COL18A1, SLC1A7, ADIPOQ, and MERTK showed colocalization support among the 8 evaluable nominal MR-associated loci. Because colocalization coverage was limited, these genes should be regarded as preliminary candidate loci rather than established therapeutic targets. The pharmacological annotations were indirect, non-selective, unsupported, or directionally inconsistent and should be interpreted solely as hypothesis-generating information.

Bayesian colocalization↗

PubMed2025

RNAi screening of uncharacterized genes identifies promising druggable targets in Schistosoma japonicum.

Schistosomiasis affects more than 250 million people worldwide and is one of the neglected tropical diseases. Currently, the treatment of schistosomiasis relies on a single drug-praziquantel-which has led to increasing pressure from drug resistance. Therefore, there is an urgent need to find new treatments. The development of genome sequencing has provided valuable information for understanding the biology of schistosomes. In the genome of Schistosoma japonicum, approximately 11% of the protein-coding sequences are uncharacterized genes (UGs) annotated as "hypothetical protein" or "protein of unknown function." These poorly understood genes have been unjustifiably neglected, although some may be essential for the survival of the parasites and serve as potential drug targets. In this study, we systematically mined the highly expressed UGs in both genders of this parasite throughout key developmental stages in their mammalian host, using our previously published S. japonicum genome and RNA-seq data. By employing in vitro RNA interference (RNAi), we screened 126 UGs that lack homologs in Homo sapiens and identified 8 that are essential for the parasite vitality. We further investigated two UGs, Sjc_0002003 and Sjc_0009272, which resulted in the most severe phenotypes. Fluorescence in situ hybridization demonstrated that both genes were expressed throughout the body without sex bias. Silencing either Sjc_0002003 or Sjc_0009272 reduced the cell proliferation in the body. Furthermore, in vivo RNAi indicated both genes are required for the growth and survival of the parasites in the mammalian host. For Sjc_0002003, we further characterize the underlying molecular cause of the observed phenotype. Through RNA-seq analysis and functional studies, we revealed that silencing Sjc_0002003 reduces the expression of a series of intestinal genes, including Sjc_0007312 (hypothetical protein), Sjc_0008276 (vha-17), Sjc_0002942 (PLA2G15), and Sjc_0003646 (SJCHGC09134 protein), leading to gut dilation. Our work highlights the importance of UGs in schistosomes as promising targets for drug development in the treatment of the schistosomiasis.

Schistosoma japonicum↗

PubMed2026

Multi-omics uncovers the pleiotropic genetic mechanisms linking MASLD and cardiometabolic syndromes.

BACKGROUND: Metabolic dysfunction-associated steatotic liver disease (MASLD) and cardiovascular-kidney-metabolic (CKM) syndrome are interrelated conditions with shared pathophysiological features; however, the genetic architecture underlying their relationship has not been fully elucidated. Deciphering this shared genetic basis holds promise for advancing mechanistic insights and therapeutic discovery. METHODS: We performed an integrated genome-wide cross-trait analysis using GWAS summary statistics for MASLD and 38 CKM traits. Our analysis estimated genetic correlations, inferred causal relationships, and identified pleiotropic variants. Candidate causal genes and druggable targets were subsequently prioritized through integrating multi-omics data. RESULTS: MASLD exhibited significant genetic correlations with 16 CKM traits, especially metabolic and cardiovascular conditions. Bidirectional causal relationships were observed between MASLD and T2D, adiposity, and lipid traits. We discovered 116 pleiotropic loci, including 65 shared causal variants such as rs429358 near APOE, which exerted influence across multiple traits. Gene-based analyses prioritized 152 unique candidate pleiotropic genes, enriched in lipid and cholesterol metabolism, and highly expressed in the liver, adipose, and immune-related cell types, such as macrophages and endothelial cells. Multi-omics integration validated 131 genes using eQTL and pQTL data from multiple tissues and cohorts. Notably, FTO and APOE emerged as central pleiotropic hubs, and druggability evaluation highlighted APOE, LPL, PPARG, and GPBAR1 as established therapeutic targets for metabolic diseases. CONCLUSION: This study provides a comprehensive map of the shared genetic architecture between MASLD and CKM syndrome, reveals novel causal genes and repurposable drug targets, and offers insights into precision medicine approaches for cardiometabolic and liver diseases.

Humans↗

PubMed2026

Integrative proteogenomic and observational analysis identifies potential biomarkers for latent autoimmune diabetes in adults.

BACKGROUND: Latent autoimmune diabetes in adults (LADA) shares core genetic and immunological features with type 1 diabetes (T1D) but is frequently misdiagnosed as type 2 diabetes (T2D). With few biomarkers for its timely diagnosis and management, this study integrated proteome-wide Mendelian randomisation (MR) and observational clinical analysis to identify potential LADA biomarkers. METHODS: We performed proteome-wide MR using cis-protein quantitative trait loci (cis-pQTLs) for 1,389 plasma proteins from the deCODE study (n = 35,559) and genome-wide association study (GWAS) data for LADA (2,634 cases and 5,947 controls, European ancestry). Robustness was enhanced via multiple sensitivity analyses. Pathway enrichment analysis, druggability evaluation, phenome-wide MR, and interaction analyses were performed to investigate the clinical relevance and biological context of candidate proteins. Candidate proteins were further evaluated using enzyme-linked immunosorbent assays in a matched Chinese clinical study (n = 241) to assess their discriminative ability for LADA. RESULTS: Proteome-wide MR and colocalisation analyses indicated associations between genetically predicted plasma levels of C-X-C motif chemokine ligand 10 (CXCL10; OR [95% CI] per 1-SD increase in protein levels: 5.49 [1.74,17.32]), serum amyloid A1 (SAA1; 1.28 [1.14,1.45]), and SAA2 (1.22 [1.11,1.34]) with LADA risk. Replication, multi-tissue eQTL, and multivariable MR supported CXCL10's association. Druggability evaluation suggested CXCL10 as a drug target under investigation, and phenome-wide MR of 1,006 diseases and traits indicated no major safety concerns for CXCL10 as a potential biomarker. In the observational clinical study, CXCL10 differentiated LADA from healthy controls (area under the receiver operating characteristic curve [ROC-AUC]: 0.889; precision-recall area under the curve [PR-AUC]: 0.919) and T2D (ROC-AUC: 0.838; PR-AUC: 0.921), with both models showing adequate calibration. CONCLUSIONS: This study suggests that CXCL10 is a putative biomarker associated with LADA, demonstrating discriminative ability to distinguish LADA from T2D in an observational clinical cohort. These findings contribute to understanding the autoimmune molecular aetiology of LADA and support its diagnostic potential in resolving the clinical ambiguity between LADA and T2D.

Humans↗

PubMed2026

Genome-wide CRISPR screens map synthetic lethal interactions across recurrent cancer driver alterations.

Synthetic lethality (SL) provides a treatment paradigm for targeting cancer with alterations in driver genes that are not conventionally druggable, including tumor suppressor genes. We execute a series of genome-wide CRISPR screens using functionally validated isogenic cell lines and conduct a large-scale SL analysis using data from the cancer dependency map (DepMap). We chart SL interactions across 15 driver alterations: FBXW7, CCNE1, CDK12, ARID1A, KMT2D, DNMT3A, TET2, KEAP1, STK11, IDH1, SF3B1, SRSF2, U2AF1, chromosome 18q loss, and chromosome 13q loss. We show validation of several SL interactions, including ARID1A and the hexosamine biosynthetic pathway aminotransferase GFPT1, STK11 with CAMK protein kinase MARK2, FBXW7 and the CDK1 regulatory kinase PKMYT1, and CCNE1 amplification and the anaphase-promoting complex or cyclosome (APC/C). In summary, this study offers a rich resource of genetic interactions across cancer drivers enabling the discovery of biological insights and drug targets for future therapeutic development.

CP: cancer↗

PubMed2025

Genomics-informed drug-repurposing strategy identifies two therapeutic targets for preventing liver disease associated with metabolic dysfunction.

Identification of drug-repurposing targets with genetic and biological support is an economically and temporally efficient strategy for improving the treatment of diseases. We employed a cross-disciplinary approach to identify potential therapeutics for the prevention of metabolic-dysfunction-associated steatotic liver disease (MASLD) in at-risk individuals by using humans as a model organism. We identified 212 putative candidate genes associated with MASLD by using data from a large multi-ancestry genetic association study, of which 158 (74.5%) were previously unreported. From this set, we identified 57 genes that encode for druggable protein targets and for which the effects of increasing genetically predicted gene expression on MASLD risk align with the function of that drug on the protein target. We then used We then evaluated these potential targets for evidence of efficacy by using Mendelian randomization, pathway analysis, and protein structural modeling. Through these approaches, we present compelling evidence to suggest that the activation of FADS1 by icosapent ethyl, as well as S1PR2 by fingolimod, could be a promising therapeutic strategy for MASLD prevention.

Humans↗

PubMed2026

CRISPR/Cas9 loss-of-function screen in a neuronal model of AP-4 deficiency identifies ATG9A trafficking modulators.

Biallelic loss-of-function variants in adaptor protein complex 4 (AP-4) disrupt trafficking of transmembrane proteins at the trans-Golgi network, including autophagy-related protein 9A (ATG9A), leading to childhood-onset hereditary spastic paraplegia (AP-4-HSP). AP-4-HSP is characterized by features of both a neurodevelopmental and a degenerative neurological disease. To investigate the molecular mechanisms underlying AP-4-HSP and identify potential therapeutic targets, we conducted an arrayed CRISPR/Cas9 loss-of-function screen of 8,478 genes, targeting the "druggable genome," in a human neuronal model of AP-4 deficiency. Through this phenotypic screen and subsequent experiments, key modulators of ATG9A trafficking were identified, and complementary pathway analyses provided insights into the regulatory landscape of ATG9A transport. Knockdown of ANPEP and NPM1 enhanced ATG9A availability outside the trans-Golgi network, suggesting that they regulate ATG9A localization. These findings deepen our understanding of ATG9A trafficking in the context of AP-4 deficiency and offer a framework for the development of targeted interventions for AP-4-HSP.

Humans↗

PubMed2025

Possible linking and treatment between Parkinson's disease and inflammatory bowel disease: a study of Mendelian randomization based on gut-brain axis.

BACKGROUND: Mounting evidence suggests that Parkinson's disease (PD) and inflammatory bowel disease (IBD) are closely associated and becoming global health burdens. However, the causal relationships and common pathogeneses between them are uncertain. Furthermore, they are uncurable. Thus, we aimed to identify the causal relationships and novel therapeutic targets shared between them based on their common pathophysiological mechanisms in gut-brain-axis (GBA). METHODS: A meta-analysis on bidirectional Mendelian randomization (MR) utilizing various datasets was performed to estimate their causal relationship. Then, pleiotropic analysis under the composite null hypothesis (PLACO) with functional mapping combined with annotation of genetic associations (FUMA) analysis were conducted to identify pleiotropic genes. Next, blood, brain and intestine expression quantitative trait locus (eQTL) were taken to perform drug-target MR finding common causal genes in two diseases. Colocalization analysis ensured the eQTLs of corresponding gene colocalized with disease. Enrichment analysis and protein‒protein interaction (PPI) network were done to explore common pathogenesis pathways. Genes passed all analysis were regarded as drug targets. RESULTS: Our MR meta-analysis revealed the bidirectional causal relationship between diseases, with combined ORs for PD on IBD, CD, UC (1.050 [95% CI 1.014-1.086], 1.044 [95% CI 0.995-1.095], 1.063 [95% CI 1.016-1.120]); for IBD, CD, UC on PD (1.003 [95% CI 0.973-1.034], 1.035 [95% CI 1.004-1.067], 1.008 [95% CI 0.977-1.040]). Overall, 277, 216 and 201 genes were identified as pleiotropic genes between PD and IBD, CD, UC. Total of 733 genes were classified as tier 3 (found in only one tissue) druggable targets, 57 as tier 2 (found in two tissues, 51 protein-coding genes) and 9 as tier 3 (found in three tissues). Among 60 protein-coding druggable targets over tier 2, 18 overlapped with pleiotropic genes and enriched in mitochondria, antigen presentation, processing and immune cell regulation pathways. Three druggable genes (LRRK2, RAB29 and HLA-DQA2) passed colocalization analysis. LRRK2 and RAB29 were reported to be pleiotropic genes, and RAB29 and HLA-DQA2 were reported for the first time as potential drug targets. CONCLUSIONS: This study established a reliable causal relationship, possible shared drug targets and common pathogenesis pathways of two diseases, which had important implications for intervention and treatment of two diseases simultaneously.

Humans↗

PubMed2026

Genomic insights into stroke recovery: cross-phenotype associations.

Stroke is a major cause of long-term disability with variable recovery. While clinical factors such as initial severity play a role, genetic factors are increasingly recognized as important contributors to stroke recovery. Genotype studies are generally focused on a single post-stroke behavioural domain, but some genes might relate to broad mechanisms of plasticity. This study therefore aimed to identify cross-phenotypic genetic variants associated across two or more stroke recovery domains. DNA from Stroke, Stress, Rehabilitation, and Genetics study participants was genotyped, resulting in 9 814 610 variants. In order to examine cross-phenotypic results, we first conducted genome-wide association studies on the six recovery domains: motor (grip force), cognition (Telephone Montreal Cognitive Assessment), depression (Patient Health Questionnaire-8), stress (Primary Care Post-Traumatic Stress Disorder Screen), functional status (Stroke Impact Scale-Activities of Daily Living), and disability (modified Rankin Scale 0-2 versus 3-6), some of which were tested longitudinally, yielding nine phenotypes. Models were adjusted for age, sex, initial severity (NIH Stroke Scale score), and ancestry. Cross-phenotype associations were identified by evaluating single nucleotide polymorphisms (SNPs) associated (P < 5e-5) with multiple phenotypes. To determine how these genetic variants may relate to biological mechanisms of recovery, we conducted gene enrichment analyses. Participants (n = 565, 59% male) had mild-moderate initial stroke severity (median acute NIH Stroke Scale score = 4). After accounting for the correlation structure among the nine phenotypes, we observed 319 cross-phenotypic SNPs, 3.45 times the expected number. Five of the cross-phenotypic SNPs were linked to genes relevant to neural development, function and plasticity, e.g. ERICH1 (rs11778883-C), FOX3 (rs55726768-G), LIFR-AS1 (rs76401391-T), RPS6KA2 (rs113518460-C) and TUBGCP2 (rs147150392-C), as were enrichments in RAB5-EEA1, CTNNA1-CTNNB1, CIN85-SH3GL2 and ELMO1-DOCK2 complexes. Multiple gene enrichments were found, e.g. Stroke Impact Scale-Activities of Daily Living and Patient Health Questionnaire 8 at 3 months were enriched for CREB phosphorylation, which is important for long-term potentiation. We identified cross-phenotypic SNPs associated with multiple behavioural domains of stroke recovery. Some of these genes encode, or regulate, druggable proteins. These genetic factors are not well captured by clinical or neuroimaging assessments and so provide a unique window into stroke recovery. These findings, if validated, suggest that some genes may be broadly important to stroke recovery.

GWAS↗

PubMed2026

Druggable genome CRISPRi screen in hydrogels reveals regulators of cortactin-driven actin remodeling promoting glioblastoma invasion.

To identify therapeutic targets limiting glioblastoma invasion, we applied druggable genome CRISPRi screens and multiomic analysis to patient-derived glioblastoma cells in micro-dissectible biomimetic 3D hydrogels that permitted separation and analysis of core versus invasive fractions. Of 2,550 genes screened, 12 encoded druggable targets whose suppression limited invasion, of which AURKB (encoding aurora kinase B) and ACP1 (encoding low molecular weight protein tyrosine phosphatase, LMW-PTP) were validated in neurosphere assays and in vivo. Proximity labeling identified cortactin as a link between LMW-PTP and aurora B, and we observed that cortactin underwent serine phosphorylation by aurora B and tyrosine dephosphorylation by LMW-PTP. Targeting ACP1 or AURKB via CRISPRi or inhibitors in culture and in vivo shifted the cortactin phosphorylation balance in glioblastoma, reducing levels of cortactin and the actin-related protein 2/3 (Arp2/3) complex that mediates cortactin-induced actin stabilization, thereby reducing actin-cortactin-Arp2/3 colocalization and subsequent actin polymerization. AURKB or ACP1 targeting shifted actin from cytoplasm to the nucleus, reducing mesenchymal gene expression. Biophysical analysis implicated AURKB in glioblastoma cell adhesion and stiffness needed for initial migration and ACP1 in mechanical stress resistance required for later migration. These findings revealed a targetable axis balancing kinase and phosphatase activities to regulate actin polymerization during glioblastoma invasion.

Humans↗

PubMed2026

Searching the druggable genome using large language models.

SUMMARY: The druggable genome encompasses the genes that are known or predicted to interact with drugs. The Drug-Gene Interaction Database (DGIdb) provides an integrated resource for discovering and contextualizing these interactions, supporting a broad range of research and clinical applications. DGIdb is currently accessed through structured web interfaces and API calls, requiring users to translate natural-language questions into database-specific query patterns. To allow for the use of DGIdb through natural language, we developed the DGIdb Model Context Protocol (MCP) server, which allows large language models (LLMs) access to up-to-date information through the DGIdb API. We demonstrate that the MCP server improves an LLM's ability to answer questions requiring accurate, up-to-date biomedical knowledge drawn from structured external resources. AVAILABILITY AND IMPLEMENTATION: The DGIdb MCP server is detailed at https://github.com/dgidb/dgidb-mcp-server and includes instructions for accessing the server through the Claude desktop app.

Large Language Models↗

PubMed2026

Searching the Druggable Genome using Large Language Models.

SUMMARY: The druggable genome encompasses the genes that are known or predicted to interact with drugs. The Drug-Gene Interaction Database (DGIdb) provides an integrated resource for discovering and contextualizing these interactions, supporting a broad range of research and clinical applications. DGIdb is currently accessed through structured web interfaces and API calls, requiring users to translate natural-language questions into database-specific query patterns. To allow for the use of DGIdb through natural language, we developed the DGIdb Model Context Protocol (MCP) server, which allows large language models (LLMs) access to up-to-date information through the DGIdb API. We demonstrate that the MCP server greatly enhances an LLM's ability to answer questions requiring accurate, up-to-date biomedical knowledge drawn from structured external resources. AVAILABILITY AND IMPLEMENTATION: The DGIdb MCP server is detailed at https://github.com/griffithlab/dgidb-mcp-server and includes instructions for accessing the server through the Claude desktop app.

Journal Article↗

PubMed2026

A genome-wide cross-trait analysis characterizes the shared genetic architecture between rheumatoid arthritis and psychiatric disorders.

OBJECTIVES: Patients with RA have a 2- to 3-fold elevated risk of psychiatric disorders, suggesting an underlying genetic link between these phenotypes. However, the shared genetic architectures and pathological mechanisms driving RA-psychiatric disorder comorbidity remain to be fully elucidated. Herein, we performed cross-trait analysis to investigate the shared genetic architecture between RA and psychiatric disorders. METHODS: Leveraging European-ancestry genome-wide association studies (GWASs) datasets of RA (n = 1 026 690) and 10 major psychiatric disorders (n = 14 307-1 222 882), we performed cross-trait pleiotropic analysis to identify the shared pleiotropic loci and genes between RA and psychiatric disorders, followed by functional annotation and Mendelian randomization analysis to explore the pathological mechanisms underlying RA-psychiatric disorder comorbidity. RESULTS: Our analysis revealed significant positive genetic correlations between RA and seven psychiatric disorders, such as major depressive disorder. From these correlations, we identified 61 pleiotropic loci jointly influencing RA and psychiatric disorder risk, along with 208 pleiotropic genes predominantly involved in immune and inflammatory response biological processes. Druggable target exploration identified 21 drug-gene interactions involving pleiotropic genes, with two genes (RHOA and TRAF3) classified in the clinically actionable category, representing potential therapeutic targets for both RA and psychiatric disorders. Mendelian randomization further demonstrated a bidirectional causal relationship between RA and schizophrenia, while supporting the causal roles of attention-deficit/hyperactivity disorder, major depressive disorder and post-traumatic stress disorder in increasing RA risk. CONCLUSION: Our findings elucidate the shared genetic architecture between RA and psychiatric disorders, providing novel insights into the pathological mechanisms underlying their comorbidity and laying the groundwork for improved comorbidity management.

Arthritis, Rheumatoid↗

PubMed2026

Single-Cell Transcriptome-Wide Mendelian Randomization and Colocalization Uncover Potential Immunocytes-Related Therapeutic Targets for Obesity.

Weight-loss treatment is crucial for individuals with obesity to prevent various complications. The role of Immune cells in obesity has been recently recognized, whereas its translation into therapy requires identifying key target genes. We performed Mendelian randomization (MR) analysis to assess causal relationships between expression quantitative trait loci (eQTL) of 14 immune cells and obesity-related traits (obesity, body mass index and body fat percentage), and validated the results in colocalization analysis. For the putative causal genes identified by the MR and colocalization analyses, we conducted pathway enrichment, differential expressed gene (DEG) analysis and search of druggable evidence, and utilized a Tier system to prioritize drug targets for obesity. MR and colocalization evidence was observed for 1630 genes associated with one or more obesity-related traits, mainly expressed in CD4+ naive/central memory T cells and enriched in antigen processing and presentation pathways. Forty-one genes showed causal relationship with all three outcomes, among which 19 genes have not been reported for obesity previously. DEG analysis using single-cell RNA sequencing data of blood or adipose tissue indicated that the differential expression of UBE2Z in monocytes, ZCCHC7 in T cells, and FNBP4 in B cells between lean and obese individuals were consistent with the MR results. By searching drug-gene interaction databases, we found targeted drugs for PYGB and PRUNE1, and PYGB was the top gene ranked in the Tier system. This study provides evidence for the involvement of immune cells in obesity, and the potential cell-specific, immune-related targets for obesity treatment.

Obesity↗

PubMed2025

Network-based drug repurposing for psychiatric disorders using single-cell genomics.

Neuropsychiatric disorders lack effective treatments due to a limited understanding of the underlying cellular and molecular mechanisms. To address this, we integrated population-scale single-cell genomics data and analyzed 23 cell-type-level gene regulatory networks across schizophrenia, bipolar disorder, and autism. Our analysis revealed potential druggable transcription factors co-regulating known risk genes that converge into cell-type-specific co-regulated modules. We applied graph neural networks on those modules to prioritize novel risk genes and leveraged them in a network-based drug repurposing framework to identify 220 drug molecules with the potential for targeting specific cell types. We found evidence for 37 of these drugs in reversing disorder-associated transcriptional phenotypes. Additionally, we discovered 335 drug-cell quantitative trait loci (eQTLs), revealing genetic variation's influence on drug target expression at the cell-type level. Our results provide a single-cell network medicine resource that provides potential mechanistic insights for advancing treatment options for neuropsychiatric disorders.

Drug Repositioning↗