SEARCH · PubMed Health

Results for “Single-cell expression quantitative trait loci”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

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↗

PubMed2026

Integrative TWAS and multi-omics analyses prioritize HSPE1 as a candidate risk gene for bipolar disorder with immune cell-specific regulatory evidence.

BACKGROUND: Bipolar disorder (BD) is a severe psychiatric disorder associated with substantial disability. Although genome-wide association studies have identified multiple BD-associated loci, the underlying genes and mechanisms remain incompletely understood. METHODS: We integrated a European-ancestry BD genome-wide association dataset with cross-tissue and tissue-specific transcriptome-wide association studies (TWAS) and complementary gene-based analysis. Candidate genes were further evaluated using differential expression analysis, consensus clustering, immune infiltration analysis, machine learning, summary-data-based Mendelian randomization, Mendelian randomization using single-cell expression quantitative trait locus data, single-nucleus transcriptomics, phenome-wide association analysis, and virtual screening. RESULTS: The integrative analyses prioritized 37 candidate genes. Peripheral-blood differential-expression analysis identified 14 genes that remained significant after FDR correction, and their expression profiles separated BD samples into two expression-defined clusters. Machine-learning analysis selected UNC50, LMAN2L, LYG2, HSPE1, and KANSL3 for an exploratory classification nomogram. SMR associated genetically predicted higher HSPE1 expression with increased BD risk in two blood eQTL datasets. Cell-type-specific analyses indicated HSPE1-related associations in T-cell and natural killer cell subsets, while single-nucleus analysis descriptively showed higher HSPE1 expression in medial thalamic T cells from BD samples. PheWAS identified no genome-wide significant associations for HSPE1, whereas virtual screening identified candidate compounds with favorable predicted docking scores against the HSPE1 structure. CONCLUSION: This integrative multi-omics study identified HSPE1 as a candidate BD risk gene with immune-cell-related regulatory evidence, providing insight into BD pathogenesis and supporting functional validation.

Humans↗

PubMed2026

Potential mitochondria-associated pathogenic genes in sepsis: a multi-omics Mendelian randomization study.

BACKGROUND: Mitochondrial dysfunction has been implicated in the pathophysiology of sepsis. However, human genetic evidence linking mitochondria-related genes to sepsis susceptibility remains limited. This study aimed to identify mitochondria-related genes associated with sepsis risk using a multi-omics Mendelian randomization framework. METHODS: Summary-data-based Mendelian randomization (SMR) was applied using sepsis genome-wide association study (GWAS) summary statistics from the UK Biobank and FinnGen databases. Expression, methylation, single-cell, and protein quantitative trait loci (QTLs) were used as genetic instruments. Colocalization analyses were conducted to evaluate whether SMR associations were driven by shared genetic variants. Expression of prioritized candidate genes was further examined in clinical septic samples, and correlations with disease severity (SOFA scores) were assessed. RESULTS: SMR analysis prioritized 13 mitochondria-related genes associated with sepsis risk. Immune cell-specific eQTL analysis suggested that genetically predicted SURF1 expression in memory B cells and naïve T cells was associated with sepsis risk. Differential expression of 12 candidate genes was confirmed in septic patients by qPCR, and PPOX expression showed a negative correlation with SOFA scores. Integration of mQTL and eQTL data supported a regulatory relationship between methylation at cg06661924 and AK4 expression. Increased genetically predicted AK4 expression was associated with higher sepsis risk (OR = 1.21, 95% CI 1.02-1.42). Protein-level analysis identified DUT as a potential sepsis-associated candidate, with consistent evidence across streptococcal and pneumococcal septicemia subtypes. Subtype analyses also suggested heterogeneous genetic signals across different sepsis subtypes. CONCLUSION: This study prioritized several mitochondria-related genes associated with sepsis susceptibility based on human genetic evidence. These findings provide candidate targets for further mechanistic and translational investigation.

Humans↗

PubMed2026

Role of IFIT1 and IFIT3 in systemic lupus erythematosus: modeling a diagnosis and exploring immune regulation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disorder characterized by multi-organ involvement and a protracted clinical course. Current diagnostic strategies, which rely heavily on clinical symptoms and serology, are often insufficient for early detection. Therefore, highly accurate diagnostic biomarkers are urgently needed to facilitate early intervention and optimize personalized treatment strategies. D atasets GSE61635 and GSE135779 were integrated to identify differentially expressed genes. Weighted gene co-expression network analysis (WGCNA) was performed to isolate the module with the strongest clinical relevance. Mendelian randomization and single‑cell RNA‑seq were used to identify key disease‑relevant genes. A diagnostic model was then constructed, and gene set variation analysis (GSVA), along with gene set enrichment analysis (GSEA), was conducted to elucidate the underlying molecular pathways. IFIT1 and IFIT3 were identified as 2 core genes highly expressed in monocytes and T cells of SLE patients. Functional enrichment analysis revealed that these genes were enriched in immune-related pathways, metabolic pathways related to inflammation and genomic stability. The diagnostic model showed good accuracy, with an area under the curve (AUC) of 0.974 on the training set and 0.912 on the validation set. IFIT1 and IFIT3 represent promising biomarkers for diagnosing SLE and appear to mediate key immune and metabolic disturbances. Furthermore, the developed model serves as an accurate and reliable instrument for early diagnosis and personalized therapy. Large-scale clinical studies are warranted to further validate these findings and evaluate their clinical application.

Humans↗

PubMed2026

Integrative multi-omics reveals a fibroblast-centered, ZFHX3-prioritized regulatory framework linking sick sinus syndrome and atrial fibrillation.

OBJECTIVE: To define shared genetic and multi-scale mechanisms underlying comorbidity between sick sinus syndrome (SSS) and atrial fibrillation (AF). METHODS: We integrated genome-wide association study (GWAS) summary statistics for SSS and AF with Genotype-Tissue Expression (GTEx) expression and splicing quantitative trait loci (eQTL/sQTL), atrial single-cell and spatial transcriptomics, and epigenomics. We identified trait-relevant tissues and pathways, prioritized shared cell types, quantified genome-wide and local genetic sharing, detected joint loci by cross-trait meta-analysis, and linked loci to regulatory programs via colocalization and cell-prioritized co-expression networks. RESULTS: Both traits showed strongest enrichment in cardiac tissue, especially Heart Atrial Appendage. Fibroblasts from the left atrial appendage were consistently prioritized as the key shared cell population. SSS and AF displayed significant positive genome-wide genetic correlation, with multiple locally shared regions, including six major loci. Cross-trait meta-analysis identified eight joint-phenotype SNPs implicating four susceptibility genes. ZFHX3 was the leading tissue-cell-gene candidate, acting as a hub in fibroblast co-expression modules and colocalizing with cardiac regulatory signals. CONCLUSION: Shared liability for SSS and AF is highly tissue- and cell-specific, converging on regulatory networks in atrial appendage fibroblasts, with ZFHX3 serving as a central mechanistic and biomarker node.

Humans↗

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↗

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↗

PubMed2026

Deciphering novel targets in salivary gland pleomorphic adenoma by integrating plasma proteomics and parotid transcriptomics analyses.

BACKGROUND/PURPOSE: Pleomorphic adenoma (PA) is the most common salivary gland benign tumor, with its molecular drivers elusive due to a lack of experimental models. This study aimed to decipher novel targets in PA by systematically integrating plasma protein quantitative trait loci (pQTL)-based Mendelian randomization (MR) with multi-omics profiling of parotid gland tissues. MATERIALS AND METHODS: We performed two-sample MR using 5450 plasma pQTLs and genome-wide association study summary for benign or broader salivary gland diseases from FinnGen consortium. Bulk RNA-sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq) comparing PA to normal tissue were used for transcriptomic validation. Immunohistochemistry (IHC) was applied for protein-level validation in human PA, adenoid cystic carcinoma (ACC), and murine inflammatory lesions. RESULTS: MR identified 12 plasma proteins associated with benign salivary gland tumor risk. Transmembrane serine protease 6 (TMPRSS6) was the only protein significantly risk-increasing for both benign and broader salivary gland diseases. Strikingly, mitogen-activated protein kinase kinase 4 (MAP2K4) showed opposite MR effects between benign and all-lesion outcomes. Bulk RNA-seq showed limited concordance with MR findings, while scRNA-seq revealed a unique plastic epithelium and partially validated candidates at cellular resolution. Critically, IHC confirmed MAP2K4 protein overexpression specifically in human PA, but not in ACC or inflammatory lesions, while TMPRSS6 was downregulated in established pathologies despite its genetic risk association. CONCLUSION: By integrating plasma proteome-based causal inference with parotid tissue multi-omics, this study unveils MAP2K4 as a potential PA-specific driver. This integrative framework provides novel, context-specific targets for further functional investigation in salivary gland tumorigenesis.

Gene expression profiling↗

PubMed2026

Endocrine-disrupting chemical-induced gene networks confer coronary heart disease risk revealed by causal inference and single-cell analyses.

BACKGROUND: Endocrine-disrupting chemicals (EDCs) are linked to coronary heart disease (CHD), but underlying mechanisms remain unclear. We aimed to identify EDC-related genes and evaluate their causal roles in CHD. METHODS: We curated EDC-related genes from a compound-gene interaction database and integrated them with CHD genome-wide association study (GWAS) summary statistics and tissue-specific expression quantitative trait loci (eQTL) data. Two-sample Mendelian randomization (MR) and Bayesian colocalization were applied to infer causality. Functional enrichment, single-cell RNA sequencing of human coronary arteries, and EDC-gene networks were further analyzed. RESULTS: After FDR correction, 39 genes were significantly associated with CHD risk via MR. Four genes-ZNF827, FCHO1, IPO9 (protective), and RPL13 (risk-increasing)-showed strong colocalization (PPH4 > 0.9). Pathway and single-cell analyses of coronary artery tissue indicated that vascular and immune pathways mediate these effects. An interaction network highlighted associations between specific EDCs and candidate genes implicated in CHD susceptibility. CONCLUSION: This integrative genomic study provides evidence that EDCs influence CHD susceptibility through distinct gene networks, revealing potential mechanisms and molecular targets for prevention and therapy.

Humans↗

PubMed2026

Mapping the Immune cell-specific gene regulatory network in bipolar disorder: A framework from scTWMR to exploratory drug-target annotation.

BACKGROUND: Although the involvement of the immune system in the genetic susceptibility of bipolar disorder (BD) is widely acknowledged, the causal relationship between gene expression in specific immune cell subtypes and BD requires systematic elucidation. METHODS: We implemented an analytical framework integrating single-cell transcriptome-wide Mendelian randomization (scTWMR) with colocalization analysis. This approach utilized cis-expression quantitative trait loci (cis-eQTLs) derived from 14 distinct immune cell types as instrumental variables to interrogate BD genome-wide association study (GWAS) summary statistics (comprising 41,917 cases and 371,549 controls). Subsequent investigations encompassed functional enrichment analysis, protein-protein interaction (PPI) network construction, phenome-wide association study (PheWAS), and performed an exploratory drug-target annotation. RESULTS: Our analysis identified 33 gene-immune cell associations. Colocalization analysis provided robust evidence (PPH4 > 90%) for shared causal variants implicating the MAD1L1, APOM, and NFKBIL1 loci. Significantly enriched biological pathways included cell cycle regulation, circadian rhythm entrainment, and neuroinflammation. The PPI network revealed a core regulatory module centered on histone-encoding and immune-related genes. Exploratory drug-target annotation nominated compounds for further investigation for compounds targeting APOM, TMEM258, and NFKBIL1. CONCLUSION: This study systematically delineates a genetically supported regulatory network of immune cell-specific gene expression in BD, predominantly implicating CD8⁺ effector T cells, plasma cells, and B cells. The findings corroborate established pathological pathways while uncovering novel cell type-specific therapeutic targets, thereby providing a genetic framework for prioritizing candidate targets for future investigation.

Bipolar disorder↗

PubMed2026

A regulatory network underlying idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease in which genetic susceptibility interacts with epithelial, immune, and mesenchymal remodeling. Although the chromosome 11p15.5 locus contains established IPF susceptibility signals near MUC5B and TOLLIP, the broader regulatory architecture of this region remains incompletely resolved. METHODS: We integrated IPF genome-wide association study summary statistics with methylation, expression, and protein quantitative trait loci using summary-data-based Mendelian randomization (SMR). SMR-prioritized candidates were evaluated in independent transcriptomic and methylation cohorts and further contextualized using microRNA, transcription-factor, protein-interaction, machine-learning, single-cell, and spatial transcriptomic analyses. Fibrosis-associated expression patterns were assessed in a bleomycin-induced pulmonary fibrosis rat model. RESULTS: The analyses recovered the established MUC5B and TOLLIP signals and prioritized BRSK2 as a comparatively underexplored candidate supported by eQTL-based SMR and independent molecular evidence. The BRSK2 pQTL association did not pass the HEIDI test and was therefore not interpreted as convergent protein-level genetic evidence. Network analyses linked BRSK2 to cell-cycle, metabolic-stress, and senescence-related programs, while cross-cohort machine learning prioritized FOXA2, CDC25B, and NFE2 as informative network features. Single-cell and spatial analyses localized BRSK2 preferentially to fibroblast and myofibroblast compartments and to regions with greater histological fibrosis severity. In fibrotic rat lungs, BRSK2 expression increased, whereas FOXA2 and CDC25B decreased at the transcript and protein levels. CONCLUSIONS: These findings refine the molecular landscape of the chromosome 11p15.5 IPF susceptibility locus and prioritize BRSK2 as a candidate component of an IPF-associated profibrotic fibroblast state. Its causal contribution, direct regulatory relationships, and therapeutic tractability require targeted mechanistic validation.

Idiopathic Pulmonary Fibrosis↗

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↗

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

Mapping the immune-genetic architecture of Epstein-Barr virus-related phenotypes and multiple sclerosis through a single-cell genetic framework for target prioritization and pharmacologic hypothesis generation.

BACKGROUND: Multiple sclerosis (MS) is a severe neuroinflammatory disease causing substantial long-term disability. Strong epidemiologic evidence links Epstein-Barr virus (EBV) exposure with MS risk, but genetic evidence for immune target prioritization in EBV-related phenotypes remains limited. METHODS: We integrated single-cell cis-eQTL data from 14 immune cell types with GWASs of an EBV-related clinical phenotype and MS using a single-cell Mendelian randomization framework with colocalization analyses. Candidate eGenes were evaluated in independent cohorts. For multi-SNP instruments, we performed heterogeneity, pleiotropy, MR-Egger, weighted median, mode-based, and MR-PRESSO sensitivity analyses. We also conducted phenome-wide association analyses and queried DrugBank to annotate candidate compounds targeting prioritized genes. RESULTS: We prioritized 43 immune-cell-specific candidate eGenes with convergent genetic support, including 6 for the EBV-related phenotype and 37 for MS. SERPINB1 in NK cells was associated with increased risk of the EBV-related phenotype, whereas HLA-G was associated with decreased risk. For MS, APOM and MSH5 showed protective associations, while AHI1 showed cell-type-dependent, bidirectional associations across immune lineages. Colocalization and independent cohort evaluation supported these findings. Among FDR-significant multi-SNP associations, MR-Egger intercept tests did not indicate directional pleiotropy, although a small subset showed heterogeneity or MR-PRESSO signals. Phenome-wide analyses identified no significant adverse phenotypic associations among evaluable genes at the prespecified threshold. DrugBank annotation nominated sodium nitroprusside, fasudil, artenimol, and choline as hypothesis-generating compounds for experimental follow-up. CONCLUSIONS: This study provides a single-cell genetic framework for prioritizing immune-cell-specific candidate targets for EBV-related phenotypes and MS, and nominates genetically supported targets and pharmacologic hypotheses for experimental investigation.

Humans↗

PubMed2026

Immune Cell-Stratified Regulatory Contexts Associated With BMI-Related Multi-System Disease Risk: A Cell-Stratified Mendelian Randomization Study Using Single-Cell eQTL Data.

AIMS: Body mass index (BMI) is associated with multisystem disease risk, but the immune cell-specific regulatory contexts underlying BMI-related genetic associations with disease outcomes remain unclear. METHODS: We applied a cell-stratified Mendelian randomization framework integrating European-ancestry BMI GWAS data, GWAS datasets for 33 disease outcomes across five disease systems, single-cell cis-eQTL data from 28 peripheral blood immune cell types, and dynamic CD4+ T cell eQTL data. SuSiE-based colocalization was used to identify BMI-associated loci sharing causal variants with immune-cell gene expression. These variants were used as cell-stratified instruments for Mendelian randomization. RESULTS: Across 28 immune cell types, 1326 colocalized variants regulating 1426 genes were identified. In primary MR analyses, genetically predicted BMI showed Bonferroni-significant associations with 26 disease outcomes. Cell-stratified analyses identified 87 Bonferroni-significant associations across 17 disease outcomes. Cardiovascular diseases showed the broadest cell-stratified associations, followed by respiratory and metabolic diseases. CD4+ T cell regulatory contexts contributed one of the largest shares of prioritized associations, and BMI-related effects varied across CD4+ T cell activation states. Cross-disease prioritization highlighted recurrent immune feature genes, including TRAF3 and FGFR1. CONCLUSION: These findings prioritize CD4+ T cell regulatory contexts as potential immunogenetic links between BMI and multi-system disease risk, while requiring further validation in diverse populations and mechanistic models.

Humans↗

PubMed2026

Integrated single-cell transcriptomics, Mendelian randomization, and machine learning identify CEBPZ as an immune-related biomarker in oral lichen planus.

BACKGROUND: Oral lichen planus (OLP) is a chronic, immune-mediated oral mucosal disease with complex pathophysiology and potential for malignant transformation. Understanding its molecular basis is critical for the development of precise diagnostic and therapeutic strategies. OBJECTIVES: We aimed to identify key immune-related biomarkers and characterize cellular dynamics in OLP, with a particular focus on the role of CEBPZ in disease pathogenesis. MATERIAL AND METHODS: We analyzed single-cell RNA sequencing (scRNA-seq) data from OLP lamina propria samples (GSE211630) to identify disease-specific T-cell subpopulations using high-dimensional weighted gene co-expression network analysis (hdWGCNA) for oxidative stress-related gene modules.-data-based Mendelian randomization (SMR) integrated FinnGen genome-wide association study (GWAS; 342,499 Europeans) data with Genotype-Tissue Expression (GTEx) expression quantitative trait loci (eQTL) data to identify causal genes. Machine learning (ML) models (least absolute shrinkage and selection operator (LASSO) and convolutional neural network (CNN)) were developed using bulk RNA-seq datasets (GSE52130 and GSE38616) for diagnostic purposes. RESULTS: We identified OLP-specific T-cell populations (clusters 0, 3, 5, 7, 13, and 15) with enhanced migration inhibition factor (MIF) pathway signaling toward B cells and monocytes. Two oxidative stress-associated modules contained hub genes, including CEBPZ. Summary-data-based Mendelian randomization analysis identified 231 OLP-associated genes, with CEBPZ uniquely intersecting LASSO-selected markers (odds ratio (OR) = 1.057, 95% confidence interval (95% CI) = 1.013-1.102, p = 0.010). Machine learning models achieved area under the curve (AUC) values ranging from 0.653 to 0.745, with the CNN model reaching a validation accuracy of 0.735. CEBPZ showed elevated expression in OLP T cells and correlated with enhanced MIF-(CD74+CXCR4) signaling. CONCLUSIONS: This integrative approach identifies CEBPZ as a pivotal biomarker linking genetic susceptibility, oxidative stress, and immune dysregulation in OLP. Our diagnostic models offer promising tools for OLP management.

CEBPZ↗

PubMed2026

Microglial PICALM: A novel genetic driver and therapeutic target in vascular dementia.

BACKGROUND: Vascular dementia (VaD) lacks well-defined genetic mechanisms. Cell-type-specific effects of GWAS loci remain unexplored. METHODS: We integrated single‑cell eQTL data (183 donors, eight cell types) with VaD GWAS (3624 cases, 475,484 controls) using Mendelian randomization and Bayesian colocalization, replicated in an independent cohort (2074 cases, 456,366 controls). Subtype, snRNA‑seq, cell‑cell communication, PheWAS, expression profiling, and drug prediction with BBB permeability assessment were performed. RESULTS: Microglial PICALM was the only robustly replicated signal (OR = 0.8334, p = 5.3 × 10⁻⁴; colocalization PP.H4 > 0.75). The effect was strongest in multiple infarctions dementia (OR = 0.7746). Exploratory snRNA-seq analysis (4 VaD vs. 4 controls; GSE282111) provided supporting evidence for microglial PICALM enrichment and downregulation (p < 0.001). PICALM‑high microglia showed enhanced neurovascular‑ and phagocytosis‑related communication (e.g., SPP1, GAS6, GRN). PheWAS revealed no pleiotropy. In silico drug repurposing prioritised three FDA-approved BBB-penetrant compounds (disopyramide, benzocaine, amantadine) as candidates warranting further mechanistic validation. CONCLUSIONS: Microglial PICALM is identified as a likely genetic determinant of VaD, especially in the multiple infarctions subtype. Upregulating PICALM may be associated with a neuroprotective microglial phenotype, highlighting PICALM as a candidate therapeutic target warranting further experimental validation.

Humans↗

PubMed2026

Cell Type-Resolved Causal Inference and Spatial Transcriptomic Integration Reveal Immune-Specific Genetic Drivers of Autoimmune and Malignant Thyroid Disease.

BACKGROUND: Thyroid diseases, including autoimmune thyroid disease (AITD) and thyroid cancer, are characterized by immune dysregulation, yet the cell type-specific genetic mechanisms underlying these conditions remain poorly understood. Most genome-wide association studies (GWAS) have relied on bulk tissue expression quantitative trait loci (eQTL), which cannot resolve the heterogeneity of immune cell populations. METHODS: We performed two-sample Mendelian randomization (MR) analyses using single-cell cis-eQTLs from 14 immune cell subtypes (OneK1K cohort) as instrumental variables against GWAS summary statistics for four thyroid outcomes: autoimmune hyperthyroidism, autoimmune hypothyroidism, thyroid cancer and autoimmune thyroiditis. Causal associations were validated through Bayesian colocalization, phenome-wide association analysis (PheWAS) and multi-layered transcriptomic validation encompassing spatial transcriptomics of AITD tissue (GSE248205), bulk RNA-seq of thyroid cancer (GSE3678) and single-cell RNA-seq of thyroid tumours (GSE250521). gsMap spatial LD score regression was applied to map disease heritability onto spatial tissue architecture. RESULTS: We identified six Bonferroni-significant causal gene-cell type pairs for autoimmune hyperthyroidism, including protective effects of ABHD16A in naïve/immature B cells (OR = 0.440), HIST1H3H in CD8 NC T cells (OR = 0.324), HMGN4 in NK recruiting cells (OR = 0.556) and ZKSCAN4 in CD8 S100B T cells (OR = 0.427), with five pairs showing strong colocalization (PP.H4 ≥ 86%). Three pairs reached significance for autoimmune hypothyroidism, including a risk association of HLA-F in CD4 NC T cells (OR = 1.139). For autoimmune thyroiditis, FAM134B/RETREG1 showed consistent suggestive protective associations across both CD4 and CD8 NC T cells (PP.H4 ≥ 90% for both), suggesting a possible involvement of ER phagy regulation in thyroiditis susceptibility. Thyroid cancer showed a suggestive association with HLA-G in classical monocytes (OR = 1.899, PP.H4 = 53%). Spatial transcriptomic validation demonstrated progressive immune infiltration from control tissue to Graves' disease to Hashimoto's thyroiditis (7.7%-15.7%, 46.1%-54.1%, respectively) and strong spatial correlation between target gene expression and corresponding cell type enrichment (e.g., plasma cell-HLA-DQB1: r = 0.491, p < 10-300). HLA-G was independently validated in thyroid cancer bulk (log2fc = 0.542, p = 9.51 × 10-3, AUC = 0.857) and single-cell datasets. PheWAS revealed no significant associations detected for the core candidates. gsMap identified significant enrichment of autoimmune hypothyroidism heritability in gastrointestinal tract, adrenal gland and adipose tissue (all Bonferroni p < 0.002). CONCLUSIONS: This study establishes a multi-scale analytical framework integrating cell type-resolved genetic inference with spatial tissue validation, revealing distinct immunogenetic architectures underlying autoimmune versus malignant thyroid disease. Protective genetic programs in autoimmune hyperthyroidism converge on chromatin remodelling (HIST1H3H, HMGN4, ZKSCAN4) and lipid metabolism (ABHD16A) across lymphocyte subsets, whereas thyroid cancer risk involves immune escape mediated by HLA-G in myeloid cells. The ER-phagy receptor RETREG1 represents a candidate pathway warranting further investigation in autoimmune thyroiditis. These findings provide genetically supported, cell type-specific therapeutic targets and demonstrate a generalizable strategy for dissecting the immune-mediated mechanisms of complex thyroid diseases.

Mendelian randomization↗