SEARCH · Search PubMed

Results for “gsMap”

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.

7 recordsLinked to original sources

PubMed2026

Spatially Contextualized Integrative Genomics Highlights Neuronal and Glial Regulatory Programs in Low Back Pain.

PURPOSE: Low back pain (LBP) is a heterogeneous pain condition with a measurable genetic contribution, but the genes, brain cell types, and spatial tissue contexts through which inherited risk is expressed remain unclear. We aimed to define cell-type-specific and spatially contextualized genetic mechanisms underlying LBP. METHODS: FinnGen R12 LBP GWAS summary statistics (42,521 cases and 353,224 controls) were integrated with brain single-nuclei eQTL data across eight major brain cell classes. We evaluated genome-wide polygenic signal using LDSC, prioritized genes using MAGMA and PoPS, and performed brain cell-type-specific eQTL-anchored Mendelian randomization, primarily based on single-instrument Wald ratio estimates, followed by Bayesian colocalization. Spatial genetic mapping was conducted using gsMap in an E16.5 murine embryonic atlas and two adult human lumbar spinal cord Visium sections. Selected candidates were assessed by RT-qPCR in neuronal-like and astroglial-like inflammatory cell models. RESULTS: LDSC supported interpretable polygenic signal for LBP. MAGMA and PoPS showed partial gene-level convergence, with TCF4 and TMEFF2 supported by both approaches. Across 1641 tested gene-cell type exposures, significant eQTL-anchored MR associations were concentrated in excitatory neurons, oligodendrocytes, inhibitory neurons, and astrocytes. Integrated eQTL-anchored MR, colocalization, and gene-prioritization evidence highlighted CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate, but not strong, colocalization support. gsMap localized LBP-associated enrichment to neuroaxis-related embryonic regions, including brain, spinal cord, sympathetic nerve, and dorsal root ganglion, and to neuronal-like niches in adult lumbar spinal cord. RT-qPCR showed model-dependent expression changes, with QPRT and LGI4 preferentially responsive in neuronal-like SH-SY5Y cells and GMPPB and DPYSL5 responsive in astroglial-like U251 cells. CONCLUSION: These findings support neuronal and glial regulatory programs as plausible contributors to LBP genetic susceptibility and highlight CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate colocalization support. The results provide a spatially contextualized framework for candidate prioritization in LBP, while emphasizing the need for larger cell-type-specific eQTL resources and functional validation before therapeutic or mechanistic conclusions can be drawn.

Mendelian randomization↗

PubMed2026

Genetic evidence links hypertension to accelerated brain aging.

Hypertension affects one-third of adults and is a major comorbidity of neurocognitive disorders. The causal relationship, shared genetic architecture, and upstream mechanisms linking hypertension to brain aging remain unclear. Hypertension GWAS datasets from MVP and FinnGen R12 were meta-analyzed as the exposure, and a European-ancestry brain age gap (BAG) GWAS derived from the UK Biobank and LIFE-Adult cohorts was used as the outcome. MR and GSMR assessed causality. LDSC, HDL, and S-LDSC estimated genetic correlation. Four TWAS methods (MAGMA, FUSION, JTI-PrediXcan, FOCUS) mapped associations to genes, followed by SMR for causal validation and PoPS for prioritization. GSMAP with spatial transcriptomics characterized regional and cell-type enrichment. Hypertension and brain aging were genetically correlated, and MR and GSMR analyses suggested a causal effect of hypertension on increased brain age gap. TWAS identified 15 shared Hypertension-BAG genes, 10 supported by SMR. PoPS prioritized TRIM47 as the core gene. Shared signals were enriched in meninges, fiber tracts, cortical layer 1, and CA1 stratum lacunosum/radiatum, with cell-type enrichment in meninges, smooth muscle cells, oligodendrocytes, and astrocyte subtypes. Hypertension is genetically correlated with, and shows evidence of a causal effect on, accelerated brain aging. TRIM47 is a core gene bridging hypertension and BAG. GSMAP-based spatial enrichment provides a hypothesis-generating framework for understanding vascular, meningeal, and myelin-related pathways linking hypertension to increased brain age gap.

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↗

PubMed2026

Shared genetic basis and spatial cellular atlas of psoriasis and metabolic syndrome.

BACKGROUND: Psoriasis (PS) and metabolic syndrome (MetS) frequently co-occur. Characterizing their shared genetic architecture and spatially enriched cellular populations may clarify the context of their co-occurrence and generate hypotheses for functional validation. METHODS: We integrated genome-wide association study (GWAS) summary statistics for PS, MetS, and five related components with spatially resolved single-cell transcriptomic data. Global and local genetic correlations were assessed using linkage disequilibrium score regression, genetic covariance analysis, high-definition likelihood, and local analysis of variant association. A bivariate causal mixture model quantified polygenic overlap. Conditional/conjunctional false discovery rate and composite-null pleiotropy analyses identified shared susceptibility loci. Finally, gsMap evaluated trait-associated enrichment across annotated embryonic tissues at single-cell resolution. RESULTS: Genetic approaches identified significant genome-wide correlations and polygenic sharing between PS, MetS, and its components. Local and cross-trait analyses identified region-specific signals and cross-validated shared loci. gsMap revealed trait-specific tissue enrichment. PS showed the strongest enrichment in the epidermis (pCauchy = 1.0573 × 10  - ⁴), adipose tissue (pCauchy = 1.5366 × 10 - ⁴), and liver (pCauchy = 1.0167 × 10 - ³). Across MetS, FBG, HDL-C, hypertension, and TG, enriched regions mainly involved the liver, adipose tissue, and epidermis. WC enrichment was predominantly observed in adipose tissue (pCauchy = 1.7823 × 10 - ⁴), with no significant liver or epidermal enrichment. CONCLUSION: Integrating GWAS with single-cell transcriptomic and spatial information characterized shared genetic architecture between PS and MetS-related phenotypes and their spatial enrichment patterns. These findings provide a framework for generating testable hypotheses about comorbidity biology and guiding future functional and clinical validation.

Psoriasis↗

PubMed2026

Shared genetic architecture of smoking dependence and Crohn's disease: A cross-trait analysis of GWAS summary statistics.

INTRODUCTION: Smoking dependence (SD) and Crohn's disease (CD) are epidemiologically associated, but whether this relationship reflects shared genetic susceptibility remains unclear. METHODS: We conducted a cross-trait genetic analysis of SD and CD using publicly available genome-wide association study (GWAS) summary statistics from European-ancestry populations. Genome-wide genetic correlation was estimated using linkage disequilibrium score regression (LDSC) and high-definition likelihood (HDL). Pleiotropic variants were identified using PLACO and mapped to genomic loci using FUMA. Regional signal sharing was assessed by Bayesian colocalization. Functional analyses included stratified LDSC, Multi-marker Analysis of GenoMic Annotation (MAGMA), GTEx tissue analysis, and Metascape. Expression-linked candidate genes were prioritized using expression quantitative trait locus (eQTL)-based summary-data-based Mendelian randomization (SMR) with heterogeneity in dependent instruments (HEIDI) testing. Genetically informed spatial mapping of cells for complex traits (gsMap) was used for spatial mapping. RESULTS: SD and CD showed positive genetic correlation by LDSC (rg=0.2090, p=0.0008) and HDL (rg=0.3817, p=0.00106). PLACO identified 81 genome-wide significant pleiotropic SNPs, which were mapped by FUMA to three loci at 1p31.3, 5p13.1, and 12q12, represented by rs11209031, rs1395152, and rs17467116, respectively. MAGMA identified 22 FDR-significant genes, four of which remained Bonferroni significant: LRRK2, TNFRSF6B, ZGPAT, and RP4-583P15.15. Cross-trait tissue analysis showed significant enrichment of the shared genetic signal in whole blood and small intestine, while gene-set analysis highlighted inflammatory response (pbon=1.86×10-5) and T-helper 17 cell differentiation (pbon=7.37×10-4). SMR/HEIDI analysis further prioritized RPS6KB1 as a shared expression-linked candidate. Spatial mapping revealed a prominent signal in the embryonic gastrointestinal tract and gene-specific regional patterns involving LRRK2 and SLC2A13 in the adult mouse brain. CONCLUSIONS: SD and CD showed measurable shared genetic susceptibility, with convergent evidence from pleiotropic loci, immune-inflammatory pathway enrichment, tissue-level associations, and spatial transcriptomic mapping.

Crohn's disease↗

PubMed2026

Genetic and epigenetic underpinnings of biological aging: a multi-omics study integrating Mendelian randomization, spatial transcriptomics, and drug target discovery.

Inflammaging represents a hallmark of biological aging, yet the causal inflammatory mediators driving multi-dimensional epigenetic aging and their effector genes remain poorly characterized at the genetic level. We developed a four-tier analytical framework integrating causal screening, multi-omics effector gene mapping, spatial transcriptomics, and drug target evaluation. Two-sample Mendelian randomization (MR) of 91 circulating inflammatory proteins against six aging phenotypes identified IL-12B, IFNG, and IL-2 as the most robust pro-aging mediators with consistent effects across independent outcomes. Using multi-omics summary-based MR (SMR) as the core analytical engine, we integrated four-layer whole-blood molecular QTL resources eQTL (eQTLGen, n = 31,684), sQTL (GTEx, n = 755), pQTL (INTERVAL + SCALLOP, n = 34,232), and mQTL (McRae et al., n = 1,980) - with GWAS summary statistics for four epigenetic age acceleration measures. At a stringent threshold (P_SMR < 1×10⁻¹²), seven high-confidence effector genes were identified: NHLRC1, TPMT, SELP, and RIPPLY3 for IEAA; ZNF373A and PLDN for HannumAA; and EDARADD for PhenoAA. The chromosome 6p21 NHLRC1-TPMT locus, overwhelmingly driven by methylation QTL signals (-log₁₀P = 26.06), emerged as the dominant genetic node of epigenetic aging. Spatial projection via gsMap onto a mouse E16.5 embryo atlas (121,767 cells) revealed preferential enrichment in smooth muscle and lung, with EDARADD showing marked specificity in mucosal epithelium. Cross-database drug target mining classified TPMT and SELP as repurposable known targets and NHLRC1 as a high-priority novel druggable candidate. This study provides multi-omics convergent causal evidence for inflammation-driven epigenetic aging and delivers genetically anchored targets for precision anti-aging intervention.

Aging↗

PubMed2026

Deciphering the Genetic Underpinnings of Liver Cirrhosis-Heart Failure Comorbidity Through Multi-Omics: CRIM1 as a Key Endothelial Mediator.

The co-occurrence of liver cirrhosis (LC) and heart failure (HF) poses considerable clinical challenges, yet the cellular and molecular determinants of this comorbidity remain poorly characterized. To address this, we developed an integrative multi-omics pipeline encompassing GWAS meta-analysis, gsMap-based spatial transcriptomic projection, GeneEnrich functional annotation, single-cell atlas construction, seismicGWAS and ECLIPSER cell-type scoring, eCAVIAR and fastenloc colocalization, hdWGCNA network inference, scTenifoldKnk in silico gene perturbation, and GCTA-COJO fine-mapping. Quality-controlled meta-analysis yielded 12,347,758 and 9,256,862 variant-level associations for LC and HF, respectively. Spatial projection confirmed preferential enrichment of disease signals within embryonic hepatic and cardiac compartments. Pathway analyses disclosed that LC-linked loci were concentrated in lipid metabolic programs, whereas HF-linked loci implicated mitochondrial bioenergetics and lysosomal degradation. At the cellular level, endothelial cells emerged as the dominant HF-associated population. Convergent evidence from five orthogonal algorithms pinpointed CRIM1 as the sole robustly supported shared gene, selectively enriched in HF endothelial cells; virtual perturbation further identified LCP1 and PTPRC as downstream regulatory nodes. Fine-mapping of the chromosome 2 locus harboring rs12476437 revealed multiple statistically independent signals in the vicinity of CRIM1. Collectively, these findings computationally prioritize the endothelial-CRIM1 axis as a previously unappreciated candidate mechanistic bridge between LC and HF requiring experimental validation.

Humans↗