SEARCH · Search PubMed

Results for “Expression quantitative trait loci”

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 37 records · Page 2Linked to original sources

PubMed2026

Identification of putative causal associations between MicroRNAs and breast cancer via Mendelian randomization and bioinformatic analysis.

MicroRNAs (miRNAs) are implicated in breast cancer progression and prognosis. This study employed a Mendelian randomization (MR) framework to investigate causal relationships between plasma circulating miRNAs and breast cancer. miRNA expression quantitative trait loci were extracted from 2 independent cohorts. High-confidence miRNAs and their associated single-nucleotide polymorphisms were selected for 2-sample MR analyses using inverse-variance weighted and MR-Egger methods. Differential expression analysis and univariate Cox regression identified survival-associated genes in breast cancer, while enrichment analyses revealed pathways and biological processes linked to candidate targets. Pan-cancer analyses of miRNAs and targets were conducted via the ENCORI platform. Initial MR analyses in the discovery phase identified hsa-miR-100-5p, hsa-miR-125b-5p, and hsa-miR-339-5p as significantly associated with reduced breast cancer risk (P < .05), suggesting potential protective roles. A total of 1291 survival-associated differentially expressed genes were identified, with 39 overlapping targets implicated in miRNA-mediated breast cancer intervention. Enrichment analyses highlighted their involvement in cell cycle regulation and p53 signaling pathway. In the validation cohort, only hsa-miR-339-5p confirmed a protective effect on breast cancer risk, while hsa-miR-100-5p and hsa-miR-125b-5p did not reach significance. Pan-cancer profiling demonstrated aberrant miRNA expression across malignancies, prognostic relevance in multiple cancers, and significant negative correlations between miRNAs and target genes in breast tumors. Our findings provide novel insights into the causal roles of miRNAs in breast cancer pathogenesis and underscore their potential as noninvasive biomarkers and therapeutic targets. Future studies should prioritize functional validation and clinical translation of these miRNAs.

Humans↗

PubMed2026

Multistage Genetic, Transcriptomic, and Single-Cell Evidence Prioritizes MAP1LC3A among Ferroptosis-Related Genes in Glioblastoma.

Glioblastoma (GBM) remains a highly aggressive malignancy, and the contribution of ferroptosis-related genes to disease susceptibility remains incompletely understood. A genetically anchored, multistage framework was applied to prioritize ferroptosis-related genes associated with GBM. Among 483 genes curated from FerrDb V2, 315 had candidate cis-expression quantitative trait loci (cis-eQTLs) in eQTLGen, 250 retained at least three independent instruments after linkage disequilibrium clumping, and 226 yielded valid inverse-variance weighted (IVW) Mendelian randomization estimates using a GBM genome-wide association study comprising 6,183 cases and 18,169 controls. Thirty-four genes met the exploratory discovery criteria of P < 0.05 and a Benjamini-Hochberg false discovery rate (BH-FDR) < 0.20, with directionally concordant Bayesian weighted Mendelian randomization (BWMR) estimates. Replication-stage Mendelian randomization using GTEx V10 whole-blood cis-eQTLs supported four genes: ATG7, RPTOR, MAP1LC3A, and CHMP6. Evaluation across three independent tumor-control transcriptomic cohorts demonstrated that MAP1LC3A was consistently downregulated in tumor tissue and showed a significant random-effects pooled estimate (log₂ fold change, -1.273; 95% confidence interval, -1.625 to -0.920; false discovery rate = 0.016), whereas the other three genes lacked comparable cross-cohort statistical support. Single-cell virtual knockout analysis was subsequently performed in a patient-balanced subset of 2,400 malignant cells selected from 4,916 eligible cells across 20 adult IDH-wild-type GBM tumors. Across five independently seeded runs, 3, 15, 4, and 7 robust downstream genes were identified for ATG7, RPTOR, MAP1LC3A, and CHMP6, respectively. The resulting consensus sets comprised 17 unique genes, with RND3 shared across all four targets. Gene Ontology analysis indicated enrichment of cell-adhesion and cell-surface processes, whereas no KEGG or Reactome pathways remained significant after multiple-testing correction. Collectively, these findings prioritize MAP1LC3A for future experimental investigation while distinguishing genetic association, tumor-expression concordance, and computational perturbation from definitive evidence of causality or mechanism.

Humans↗

PubMed2025

Genetic effects on chromatin accessibility reveal the molecular mechanisms of complex traits in maize.

Cis-regulatory elements (CREs) are critical for modulating gene expression and phenotypic diversity in maize. While genome-wide association study (GWAS) hits and expression quantitative trait loci (eQTLs) are often enriched in CREs, their molecular mechanisms remain poorly understood. Characterizing CREs within accessible chromatin regions (ACRs) offers a powerful approach to link noncoding variants to chromatin structure alterations and phenotypic variation. Here, we generated ATAC-seq profiles from seedling leaves of 214 maize inbred lines, identifying 82 174 consensus ACRs. Notably, 39.55% of these ACRs exhibited significant population-wide chromatin accessibility variation. By mapping chromatin accessibility quantitative trait loci (caQTLs), we discovered 27 004 loci, including 1398 predicted to disrupt transcription factor (TF)-binding sites. Integration with multi-omics data revealed 7405 caACR-target gene pairs and linked 56 caACRs to GWAS signals for 51 agronomic traits, with significant enrichment in flowering-related pathways. Functional candidates such as ZmZIM30 - putatively regulated by caACRs - emerged as key regulators of flowering time. At the fad7 locus associated with linolenic acid content, allelic variants overlapping a caQTL showed differential chromatin accessibility. Our study provides a high-resolution cis-elements of maize leaves, deciphers the genetic basis of chromatin accessibility variation, and bridges noncoding caQTLs to molecular mechanisms underlying GWAS hits.

Zea mays↗

PubMed2026

Mapping the regulatory architecture of circadian clock adaptation: A genome-wide eQTL analysis in Drosophila melanogaster.

The circadian clock enables organisms to align internal daily rhythms with environmental cues, with major consequences for survival and fitness. Although the molecular framework of this system in Drosophila melanogaster is well characterized through transcription translation feedback loops involving ten core clock genes, the genetic basis of natural variation in their expression remains poorly understood. Here, we used natural expression variation to identify expression quantitative trait loci (eQTLs) through genome-wide association mapping. Using the Drosophila Genetic Reference Panel, we measured relative expression of all core clock genes at a single time point two hours after light onset. We identified 109 significant SNPs and 28 indels associated with expression variation across the clock network. Expression levels varied widely, with Pdp1ε showing the greatest variation (an 86-fold difference between extreme lines) and cyc the least (11.3-fold). Only three significant SNPs were located within clock genes themselves, all in Clk, whereas most associations represented trans-eQTLs in genes with diverse molecular functions. Candidate regulators included transcription factors such as Abd-B, tai, and E5; RNA binding proteins including Pum, Bru-3, and Mbl; and several long noncoding and antisense RNAs. Variants were also detected in gbb and the BMP pathway transcription factor Mad. Consistent with this, Mad knockdown reduced vri expression. Together, these results reveal a complex regulatory architecture underlying natural variation in circadian gene expression.

Journal Article↗

PubMed2026

Multi-omic characterization of the Hispanic/Latino blood lipidome reveals an additional locus and attenuated genetic prediction.

While lipids have been extensively investigated, genetic regulation of the circulating lipidome in diverse populations remains poorly understood. We conducted a lipidome-wide genome-wide association study (GWAS) of 830 lipid species in 2,287 Hispanic/Latino participants and performed predictive modeling across omics layers. We identified 7,593 genome-wide significant SNPs mapping to 208 genes. Conditional analysis disentangled the long-range linkage disequilibrium artifacts from the pleiotropic FADS1/2/3 cluster. Separately, we discovered an association at the GPLD1 locus for a circulating ceramide. Colocalization revealed shared genetic architecture with conventional lipids alongside distinct, species-specific pathways. Incorporating Native/Indigenous American expression quantitative trait loci (eQTLs) within a multi-omic framework uncovered 62 likely regulatory genes missed by European-centric gene expression models. Finally, genetically regulated predictive models demonstrated performance declining from transcriptomics to proteomics to lipidomics, reflecting increased distance from gene action along the molecular cascade. Our study provides a genetic landscape of lipid metabolism in a highly burdened population and highlights the challenges in predicting lipid abundance.

Hispanic/Latino population↗

PubMed2025

Association Analysis of the Circulating Proteome With Sarcopenia-Related Traits Reveals Potential Drug Targets for Sarcopenia.

BACKGROUND: Sarcopenia severely affects the physical health of the elderly. Currently, there is no specific drug available for sarcopenia. This study aims to identify pathogenic proteins and druggable targets for sarcopenia through Mendelian randomization (MR)-based analytical framework. METHODS: A sequential stepwise screening method that includes two-sample MR, Steiger filtering test and colocalization (MRSC) was applied to identify causal proteins associated with sarcopenia-related traits. In the MR analyses, 4372 circulating proteins with valid instrumental variables (IVs) from eight proteomic genome-wide association studies were utilized as exposures, and nine sarcopenia-related traits were utilized as outcomes. IVs were classified into cis-protein quantitative trait loci (pQTLs) and trans-pQTLs based on their positions. We conducted cis-only MRSC analyses and cis + trans MRSC analyses using cis-pQTLs and cis + trans pQTLs as IVs, respectively. Post-MRSC analyses were conducted on the prioritized findings of MRSC, including annotation of protein-altering variants (PAVs), assessment of overlap between pQTLs and expression quantitative trait loci (eQTLs), protein-protein interaction (PPI) analysis, pathway enrichment analysis and annotation of drug targets. Utilizing data from the UK Biobank, we performed an observational study to explore the associations between baseline circulating protein levels and the longitudinal changes in nine sarcopenia-related traits. RESULTS: A total of 181 causal associations for 65 proteins were prioritized by the cis-only MRSC analyses and 227 associations for 91 proteins were prioritized by the cis + trans MRSC analyses. Among the prioritized proteins, the majority of them employed non-PAVs as IVs and most of their cis-pQTLs overlapped with corresponding eQTLs and exhibited consistent directionality, with only one trans-pQTL overlapping with an eQTL. The PPI network of cis-only MRSC-prioritized proteins (p = 4.04 × 10-4) and cis + trans MRSC-prioritized proteins (p = 8.76 × 10-5) showed significantly more interactions than expected. Reactome, KEGG and GO pathway enrichment analyses for cis-only MRSC-prioritized proteins identified 52, 12 and 79 enriched pathways, respectively (adjusted p < 0.05). For proteins identified by cis + trans MRSC analyses, only 15 pathways were enriched through the GO pathway enrichment analyses. In the observational study, 197 circulating proteins were identified to be associated with one or more sarcopenia-related traits (p < 0.05/2923). Among them, the significant associations of CTSB (negative association) and ASGR1 (positive association) with sarcopenia-related traits were observed to have consistent directional associations in both MR-based studies and observational studies. Drug target annotations suggested that 52 MRSC-prioritized proteins and 145 biomarkers are drug targets or druggable. CONCLUSIONS: This study identified 89 potential pathogenic proteins and 197 candidate biomarkers for sarcopenia, providing valuable clues for the development of therapeutic drugs for sarcopenia.

Humans↗

PubMed2026

Alternative tandem transcription initiation links noncoding variants to human disease through translational control.

Alternative tandem transcription initiation is a pervasive mechanism of gene regulation, yet its genetic impact on human disease remains largely unknown. Here, we systematically quantify the genetic regulation of alternative tandem transcription initiation across 25,859 samples from 49 normal human tissues and 33 tumor tissues. We identify approximately 0.4 million genetic variants associated with alternative transcription initiation in 5295 genes, with 32% operating independently of gene expression. Moreover, we discover 2238 multi-tissue alternative tandem transcription initiation outliers enriched for rare deleterious promoter and 5' UTR variants, demonstrating that both common and rare variants modulate transcription initiation. Strikingly, 74% of disease variants that colocalize with genetic variants regulating alternative transcription initiation cannot be identified through expression quantitative trait loci. Transcriptome-wide association studies identify 614 disease susceptibility genes associated with alternative transcription initiation, including known cancer drivers such as MAFF and MLLT10. Functional validation uncovers OSGEP as a breast cancer risk gene, where the alternative allele lengthens the 5' UTR and reduces protein abundance through upstream open reading frame-mediated translation repression, and suppresses breast cancer cell proliferation. Our findings establish alternative transcription initiation as a major, underappreciated mechanism associating noncoding variation with disease, providing a critical resource for interpreting disease risk loci.

Humans↗

PubMed2026

Transcriptome-wide association analysis of Alzheimer's disease: construction and clinical validation of transcriptomic risk scores.

Early identification of individuals at high risk for Alzheimer's disease (AD) is crucial for disease prevention and intervention. This study aims to develop AD-specific transcriptomic risk scores (TRSs) through multi-tissue transcriptome-wide association study (TWAS) and to evaluate its clinical utility in AD diagnosis and risk prediction. Using GWAS summary statistics combined with expression quantitative trait loci (eQTL) data from 14 tissues, a multi-tissue TWAS approach was applied to identify AD-associated genes. Peripheral blood RNA expression data from the ADNI and GEO databases were used to construct the AD-specific TRSs. The associations of TRSs with AD pathological features and cognitive function were assessed in two independent cohorts. Furthermore, the diagnostic performance, differential diagnostic capability, and risk prediction efficiency of TRSs were evaluated. The TWAS identified 131 genes significantly associated with AD. The TRSs were significantly elevated in patients with AD and mild cognitive impairment (MCI) compared to cognitively normal (CN) individuals, and showed significant correlations with AD pathological markers and cognitive performance. When combined with APOE4 status, the TRSs demonstrated robust diagnostic ability for AD and MCI. When combined with age, the TRSs showed good diagnostic performance in distinguishing AD from frontotemporal dementia (FTD) (AUC = 0.86). Additionally, the TRSs effectively predicted the risk of progression to AD in non-AD individuals (HR = 1.74). The AD-specific TRSs developed in this study shows promising clinical utility in AD diagnosis, differential diagnosis, and risk prediction, providing valuable translational medical evidence for early screening and precision prevention of Alzheimer's disease.

Humans↗

PubMed2026

A cross-sectional study of oxidative stress pathway genotypes and their interactions with environmental pollutant levels identifies associations with gene expression and lung function.

BACKGROUND: Asthma is a heterogeneous disease influenced by genetic and environmental factors. Fine particulate matter (PM2.5) exacerbates asthma, likely through oxidative stress pathways, but whether genetic variation modifies this effect remains unclear. METHODS: We analysed data on 948 adults with asthma from the Severe Asthma Research Program (SARP), linking ZIP-code-level PM2.5 exposure with whole-genome sequencing data. We tested 4337 single nucleotide polymorphisms (SNPs) in 120 oxidative stress pathway genes for gene-environment (GxE) interactions with PM2.5 on lung function (forced expiratory volume in 1 s [FEV1] % predicted) using weighted linear regression. Gene expression data from bronchial epithelial cells (n = 170) were used to assess cis-expression quantitative trait loci (eQTLs). FINDINGS: Higher PM2.5 exposure was associated with lower FEV1% predicted (β per μg/m3 = -0.7, p = 0.01). We identified 20 SNPs across seven genes (OXSR1, PXDN, TPO, LRRK2, APP, MSRA, MSRB2) with significant GxE interactions after multiple-testing correction. Five SNPs were also eQTLs, linking PM2.5-modified gene expression to lung function. Minor alleles in OXSR1 and PXDN were associated with reduced gene expression and worsened FEV1% under high PM2.5 exposure. Conversely, TPO variants were associated with higher baseline expression and lower lung function, but under increasing PM2.5 exposure, minor allele carriers showed suppressed TPO expression and improved FEV1%. INTERPRETATION: This study identified 20 SNPs in oxidative stress pathway genes that modify the effect of PM2.5 on lung function in asthma. These findings highlight the importance of integrating environmental context in genetic studies and suggest potential therapeutic targets for pollution-sensitive asthma phenotypes. FUNDING: Supported by NIH grants.

Cross-Sectional Studies↗

PubMed2026

Integrative genetic and transcriptomic analyses prioritize CDC16 as a candidate marker for gastric cancer.

BackgroundGastric cancer (GC) remains a major cause of cancer-related mortality, and biomarkers for early detection are needed.MethodsStomach and blood expression quantitative trait loci were integrated with two GC genome-wide association studies using Mendelian randomization (MR), Bayesian colocalization, and summary-data-based MR/heterogeneity in dependent instruments (SMR/HEIDI) testing. Bulk and single-cell transcriptomic analyses characterized candidate expression and lesion-associated patterns. CDC16 protein expression was evaluated by immunohistochemistry in 53 paired GC and non-neoplastic tissues, followed by paired and exploratory receiver operating characteristic analyses.ResultsMR prioritized PILRB, CDC16, and GABPB1-AS1; SMR/HEIDI provided complementary support, while colocalization for CDC16 and GABPB1-AS1 was suggestive and model-dependent. Bulk-tissue CDC16 abundance was higher in GC, but the modest TCGA-STAD tumor-normal difference (log2FC = 0.210, FDR = 0.019) was attenuated after proliferation adjustment (log2FC = -0.002, FDR = 0.987), indicating close coupling with proliferative activity. Single-cell analysis localized CDC16 predominantly to epithelial populations, and the proportion of CDC16-detectable epithelial cells increased across lesion categories (ρ = 0.735; permutation P = 0.031). CDC16 H-scores were higher in GC than in paired non-neoplastic tissues (161.15 ± 45.11 vs 102.15 ± 54.50; P < 0.001), with higher cancer-tissue scores in 41 of 53 cases. Exploratory AUC was 0.794 (95% CI, 0.704-0.874; sensitivity, 66.0%; specificity, 79.2%).ConclusionsConvergent genetic, transcriptomic, and protein-level evidence prioritizes CDC16 as a GC-associated candidate tissue marker whose expression is closely linked to proliferative activity. Prospective validation in independent cohorts, including appropriate disease controls and blood-based evaluation, is warranted.

Stomach Neoplasms↗

PubMed2024

Lipid-lowering drugs and risk of rapid renal function decline: a mendelian randomization study.

BACKGROUND: Chronic kidney disease (CKD) patients face the risk of rapid kidney function decline leading to adverse outcomes like dialysis and mortality. Lipid metabolism might contribute to acute kidney function decline in CKD patients. Here, we utilized the Mendelian Randomization approach to investigate potential causal relationships between drug target-mediated lipid phenotypes and rapid renal function decline. METHODS: In this study, we utilized two methodologies: summarized data-based Mendelian randomization (SMR) and inverse variance-weighted Mendelian randomization (IVW-MR), to approximate exposure to lipid-lowering drugs. This entailed leveraging expression quantitative trait loci (eQTL) for drug target genes and genetic variants proximal to drug target gene regions, which encode proteins associated with low-density lipoprotein (LDL) cholesterol, as identified in genome-wide association studies. The objective was to investigate causal associations with the progression of rapid kidney function decline. RESULTS: The SMR analysis revealed a potential association between high expression of PCSK9 and rapid kidney function decline (OR = 1.11, 95% CI= [1.001-1.23]; p = 0.044). Similarly, IVW-MR analysis demonstrated a negative association between LDL cholesterol mediated by HMGCR and kidney function decline (OR = 0.74, 95% CI = 0.60-0.90; p = 0.003). CONCLUSION: Genetically predicted inhibition of HMGCR is linked with the progression of kidney function decline, while genetically predicted PCSK9 inhibition is negatively associated with kidney function decline. Future research should incorporate clinical trials to validate the relevance of PCSK9 in preventing kidney function decline.

Mendelian Randomization Analysis↗

PubMed2026

Integrative Multi-Omics Analysis Identifies Thrombosis-Associated Molecular Features Linked to Germline Susceptibility and Immune Cell Communication in Gastric Cancer.

Emerging evidence indicates that coagulation-related molecular programs are associated with thrombosis, tumor progression, and molecular dysregulation in gastric cancer (GC). However, thrombosis-associated molecular features in GC and their potential links to inherited susceptibility remain insufficiently understood. Integrated analyses of transcriptomic data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were performed to identify thrombosis-associated genes and establish a machine learning-based prognostic signature. Genome-wide association study (GWAS), expression quantitative trait loci (eQTL), transcriptome-wide association study (TWAS), and Mendelian randomization (MR) analyses were conducted to investigate susceptibility-associated transcriptional programs in GC. Functional assays were used to evaluate candidate genes associated with malignant phenotypes. Single-cell RNA sequencing (scRNA-seq) and cell-cell communication analyses were further performed to characterize cell-type-specific expression patterns and potential intercellular interactions. A total of 22 differentially expressed thrombosis-associated genes were identified, and a prognostic signature comprising 14 genes was established. The signature stratified patients into high- and low-risk groups and showed prognostic performance in both the training and validation cohorts. Integrative GWAS, eQTL, and TWAS analyses identified susceptibility-associated transcriptional programs that were positively correlated with the thrombosis-associated risk score. Silencing ACTN2 and CRYAB significantly reduced GC cell migration and invasion. scRNA-seq analysis revealed relatively high CRYAB expression in neutrophils, and CellChat analysis suggested potential neutrophil-B cell interactions involving COLLAGEN-related signaling. This integrative multi-omics study identified a thrombosis-associated molecular signature linked to prognosis and germline susceptibility-associated transcriptional programs in GC. ACTN2 and CRYAB may represent candidate genes associated with GC cell migration and invasion, while single-cell analysis suggested potential immune-related communication features.

Humans↗

PubMed2026

Prioritizing Parkinson's disease risk-associated mitochondrial candidate genes via multi-omics integrative analysis.

BACKGROUND: Mitochondrial dysfunction has been implicated in Parkinson's disease (PD), but the genetically regulated mitochondrial genes associated with PD risk remain incompletely defined. METHODS: We conducted a summary-data-based genetic epidemiology study integrating summary-based Mendelian randomization (SMR), Heterogeneity in dependent instruments (HEIDI) filtering, and Bayesian colocalization to prioritize mitochondrial-related molecular features associated with PD risk. Mitochondrial-related genes were defined using MitoCarta3.0. Genetically predicted gene expression and plasma protein abundance were evaluated using expression quantitative trait loci (eQTL) data from eQTLGen and GTEx v8, and protein quantitative trait loci (pQTL) data was assessed using International Parkinson's Disease Genomics Consortium (IPDGC) as the discovery genome-wide association study (GWAS) and FinnGen as the replication dataset. Prespecified QTL analyses were interpreted using FDR correction, HEIDI filtering, and colocalization support. DNA methylation QTL analysis, mitochondrial phenotype MR, and single-nucleus RNA-seq analysis were performed as complementary analyses. RESULTS: In the primary eQTL analysis, higher genetically predicted TTC19 expression was associated with lower PD risk (OR = 0.80, 95% CI: 0.74-0.87, PPH4 = 0.80), whereas higher MALSU1 expression was associated with increased PD risk (OR = 2.21, 95% CI: 1.59-3.06, PPH4 = 0.96). Both associations survived FDR correction, passed HEIDI filtering, and showed colocalization support. GTEx whole-blood data supported the direction of the TTC19 association. No mitochondrial protein reached significance after FDR correction and colocalization filtering in the primary pQTL analysis. Complementary methylation analysis highlighted cg06270993 as an exploratory regulatory signal for MALSU1. CONCLUSIONS: This MR-colocalization study prioritizes TTC19 and MALSU1 as genetically supported mitochondrial-related candidate genes associated with PD risk. Further validation is required to define their functional roles in PD pathogenesis.

Humans↗

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↗

PubMed2025

Gut fungi are associated with human genetic variation and disease risk.

Human genetic determinants of the gut mycobiome remain uninvestigated despite decades of research highlighting tripartite relationships between gut bacteria, genetic background, and disease. Here, we present the first genome-wide association study on the number and types of human genetic loci influencing gut fungi relative abundance. We detect 148 fungi-associated variants (FAVs) across 7 chromosomes that statistically associate with 9 fungal taxa. Of these FAVs, several occur in the protein-coding genes PTPRC, ANAPC10, NAV2, and CDH13. Additional FAVs link to tissue-specific gene expression as fungi-associated expression quantitative trait loci. Notably, the relative abundance of gut yeast Kazachstania associates with genetic variation in CDH13 encoding T-cadherin, a protein linked to cardiovascular disease. Kazachstania forms a causal relationship with cardiovascular disease risk in a mendelian two-sample randomization analysis. These findings establish previously unrecognized connections between human genetics, gut fungi, and chronic disease, broadening the paradigm of human-microbe interactions in the gut to the mycobiome.

Humans↗

PubMed2024

Relationship between inflammation/immunity and epilepsy: A multi-omics mendelian randomization study integrating GWAS, eQTL, and mQTL data.

OBJECTIVES: Increasing evidence suggests that activated innate/adaptive immunity induces an inflammatory response, thereby participating in epileptogenesis. However, the biological explanation of inflammation/immunity as a potential cause for epilepsy remains largely unknown. This research aimed to determine the causal effects of inflammation/immune-related genes in epilepsy based on multi-omics mendelian randomization (MR). METHODS: We employed summary-data-based MR (SMR) approach to combine GWAS for epilepsy (12,891 cases and 312,803 control) with gene expression quantitative trait loci (cis-eQTL, 31,684 participants) and DNA methylation QTL (cis-mQTL, 1,980 participants) data. Five additional MR methods were then used for sensitivity analyses to confirm the reliability of causal associations. In addition, enrichment analysis of key genes was conducted to provide insight into the biological functions of epilepsy risk variants. RESULTS: A total of 386 inflammation/immune-related genes were selected for further analyses. Primary SMR analysis indicated that 37 DNA methylation sites and six genes regulated by them had potential causal relationship with epilepsy. MR analysis further refined the results, identifying three genes that had a causal effect on epilepsy. Notably, VEGFA (OR: 0.925; 95 % CI: 0.862-0.994) expression was negatively correlated with epilepsy risk, whereas IL16 (OR: 1.076; 95 % CI: 1.028-1.126) and HLA-DPA1 (OR: 1.041; 95 % CI: 1.009-1.074) expressions were positively associated with epilepsy risk. Functional enrichment analysis revealed that the identified genes were involved in GO-BP terms related to VEGF activation signaling and chemotaxis regulation. CONCLUSION: This analysis confirms the causal role of inflammation/immunity in epilepsy, and the identified candidate genes provide clues for drug development in clinical practice.

Humans↗

PubMed2026

Genetic evidence supports the prioritization of CD40 among prespecified immune-related candidate drug targets in myasthenia gravis.

AIM: To prioritize prespecified immune-related candidate drug targets in myasthenia gravis for further validation based on integrated genetic evidence. METHODS: We integrated drug-target Mendelian randomization (MR) using cis-expression quantitative trait loci (cis-eQTLs), protein-level MR of plasma CD40 abundance using plasma protein quantitative trait loci (pQTLs), and colocalization analyses to evaluate genetically proxied associations with overall MG, early-onset myasthenia gravis (EOMG), and late-onset myasthenia gravis (LOMG). RESULTS: In this study, CD40 showed the most consistent genetic evidence among the six prespecified targets. Effect estimates are reported as odds ratios (ORs) with 95% confidence intervals (CIs). Higher CD40 expression proxied by cis-eQTLs was associated with increased risk of overall MG (OR = 1.14, 95% CI: 1.05-1.24, Bonferroni-adjusted p = 0.022) and EOMG (OR = 1.32, 95% CI: 1.12-1.56, Bonferroni-adjusted p = 0.015). Genetically predicted higher plasma CD40 protein abundance was associated with increased overall MG risk (OR = 1.31, 95% CI: 1.08-1.57, Bonferroni-adjusted p = 0.010), whereas the protein-level MR result for EOMG was directionally consistent but not statistically significant. Colocalization analysis provided suggestive but not definitive evidence of colocalization between CD40 expression and EOMG risk. FCGRT, IL2RA, and SYK showed additional exploratory MR signals requiring further validation. CONCLUSION: CD40 showed the most consistent genetic support among the prespecified targets, supporting its prioritization for functional validation and further therapeutic investigation in MG.

CD40↗

PubMed2026

Immune cell-specific genetic architecture of Alzheimer's disease revealed by multi-omics analysis for therapeutic target discovery and prioritization.

Alzheimer's disease (AD) is a multifactorial neurodegenerative condition in which accumulating genetic and molecular evidence implicates dysregulation of peripheral immune processes in disease pathogenesis. Nevertheless, the contribution of distinct peripheral immune cell subsets and associated gene regulatory landscapes to AD risk remains incompletely defined. To address this gap, we integrated single-cell expression quantitative trait loci (sc‑eQTL) data from the OneK1K cohort with AD GWAS summary statistics. We systematically interrogated immune cell-specific genes for their contributions to AD risk by integrating genetic causal inference with Bayesian colocalization analyses, and identified 24 eGenes that passed both the MR significance threshold (P < 0.05) and the criterion for strong shared genetic signals (PP.H4 > 0.8). Notable candidates included GATS, HLA-DOB, HLA-DQA1, PM20D1, and others, with each gene demonstrating a cell-type-specific association restricted to its corresponding immune cell type, such as monocytes, CD8 + T cells, or B cells. Independent peripheral blood single-cell transcriptomic data further supported disease-associated shifts in cell-type-specific expression patterns in AD. Phenome-wide association studies (PheWAS) indicated limited associations with off-target traits, indicating a favorable safety profile for therapeutic intervention, with the exceptions of B4GALNT3, PM20D1, and CNN2. Integration of immune gene targets with pharmacological databases yielded three candidate compound, including NSC321521 (targeting HLA-DQA1), phenoxybenzamine (targeting GSTP1), and rimexolone (targeting BIN1). Among these compounds, Predicted blood-brain barrier permeability was observed only for phenoxybenzamine and rimexolone, with docking studies indicating stable interactions, such as those between NSC321521 and HLA-DQA1, phenoxybenzamine and GSTP1, and rimexolone and BIN1. This integrative approach highlights key immune‑cell‑specific genes involved in AD and proposes repurposable drugs with central nervous system potential, paving the way for more targeted immunomodulatory strategies in AD.

Humans↗