SEARCH · Search PubMed

Results for “eQTL”

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 109 records · Page 6Linked to original sources

PubMed2006

How replicable are mRNA expression QTL?

Applying quantitative trait analysis methods to genome-wide microarray-derived mRNA expression phenotypes in segregating populations is a valuable tool in the attempt to link high-level traits to their molecular causes. The massive multiple-testing issues involved in analyzing these data make the correct level of confidence to place in mRNA abundance quantitative trait loci (QTL) a difficult problem. We use a unique resource to directly test mRNA abundance QTL replicability in mice: paired recombinant inbred (RI) and F(2) data sets derived from C57BL/6J (B6) and DBA/2J (D2) inbred strains and phenotyped using the same Affymetrix arrays. We have one forebrain and one striatum data set pair. We describe QTL replication at varying stringencies in these data. For instance, 78% of mRNA expression QTL (eQTL) with genome-wide adjusted p < or = 0.0001 in RI data replicate at a genome-wide adjusted p < 0.05 or better. Replicated QTL are disproportionately putatively cis-acting, and approximately 75% have higher apparent expression levels associated with B6 genotypes, which may be partly due to probe set generation using B6 sequence. Finally, we note that while trans-acting QTL do not replicate well between data sets in general, at least one cluster of trans-acting QTL on distal Chr 1 is notably preserved between data sets.

Animals↗

PubMed2026

Esketamine multi-omic biomarker evaluation in major depressive disorder (EMBER-MDD): concept, objectives and methodologies of a non-clinical investigator-initiated study.

Treatment resistance (TR) in major depressive disorder (MDD) affects a substantial minority of patients and is hard to recognize early, delaying intensified care. The Esketamine multi-omic biomarker evaluation in MDD (EMBER-MDD) is a non-interventional, investigator-initiated, in-vitro study within the EU Psych-STRATA programme, analyzing biospecimens collected in the randomized INTENSIFY study and the mirror OBS-TR cohort after participants complete treatment. EMBER-MDD aims to discover individual-omic and integrated multi-omic (hypothesis-free) biomarkers and signatures associated with TR risk, and molecular correlates of clinical response to esketamine nasal spray versus treatment as usual (TAU). Biomaterials will derive from approximately 420 adults with MDD (estimated n = 210 esketamine; n = 210 TAU) and include whole blood, RNA-stabilized whole blood, plasma and serum, sampled at baseline and, when feasible, during and after treatment (up to ~ 5,040 aliquots stored at - 80 °C). Genomics will use baseline DNA genotyping on Illumina Infinium GSA v3.0+MD arrays; epigenomics will profile genome-wide DNA methylation across time points using MethylationEPIC v2.0; transcriptomics will employ mRNA-seq (NovaSeq X/ X Plus); and proteomics/ metabolomics will be generated using high-throughput Olink and/ or Biocrates platforms. Each layer will undergo state-of-the-art preprocessing and analyses (e.g., GWAS/ PRS, EWAS, differential expression, WGCNA, pathway and network analyses), followed by integrative strategies including QTL mapping (meQTL/ eQTL/ pQTL/ mQTL) and intermediate-fusion machine learning with nested cross-validation, explainable AI (SHAP/ LIME) and treatment-effect modelling. All outputs are research-only and will not support individual efficacy, tolerability, or clinical decision-making. The study will deliver robust biosignatures and mechanistic hypotheses to guide future validation and inform stratified, molecularly guided intervention strategies in subsequent prospective trials. Trial registration number: 2023-506617-21-00 and 2025-178-f-S.

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

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↗

PubMed2004

Molecular markers from the transcribed/expressed region of the genome in higher plants.

In recent years, molecular marker technology in higher plants has witnessed a shift from the so-called random DNA markers (RDMs), developed in the past arbitrarily from genomic DNA and cDNA, to the molecular markers representing the transcriptome and the other coding sequences. These markers have been described as gene targeted markers (GTMs). Another specific class of markers includes the so-called functional markers (FMs), which are supposed to have a cause and effect relationship with the traits of interest. In this review, we first describe the development of these markers representing the transcriptome or genes per se; we then discuss the uses of these markers in some detail and finally add a note on the future directions of research and the implications of the wider application of these markers in crop improvement programmes. Using suitable examples, we describe markers of different classes derived from cDNA clones, expressed sequence tags (ESTs), gene sequences and the unique (coding) sequences obtained through methyl filtration or genome normalization (high C(0) t fraction) from gDNA libraries. While we briefly describe RFLPs, SSRs, AFLPs and SNPs developed from the transcriptome (cDNA clones and EST databases), we have discussed in more detail some of the novel markers developed from the transcriptome and specific genes. These novel markers include expressed sequence tag polymorphisms (ESTPs), conserved orthologue set (COS) markers, amplified consensus genetic markers (ACGMs), gene specific tags (GSTs), resistance gene analogues (RGAs) and exon-retrotransposon amplification polymorphism (ERAP). Uses of these markers have been discussed in some detail under the following headings: development of transcript and functional maps, estimations of genetic diversity, marker-assisted selection (MAS), candidate-gene (CG) approach and map-based cloning, genetical genomics and identification of eQTLs, study of genome organization and taxonomic and phylogenetic studies. At the end, we also append a list of websites relevant to further studies on the transcriptome. For want of space, considerable information including voluminous data in the form of 12 tables, and a long list of references cited in these tables, has been placed on the Internet as electronic supplementary material (ESM), which the readers may find useful.

Chromosome Mapping↗

PubMed2006

Differential gene expression of wheat progeny with contrasting levels of transpiration efficiency.

High water use efficiency or transpiration efficiency (TE) in wheat is a desirable physiological trait for increasing grain yield under water-limited environments. The identification of genes associated with this trait would facilitate the selection for genotypes with higher TE using molecular markers. We performed an expression profiling (microarray) analysis of approximately 16,000 unique wheat ESTs to identify genes that were differentially expressed between wheat progeny lines with contrasting TE levels from a cross between Quarrion (high TE) and Genaro 81 (low TE). We also conducted a second microarray analysis to identify genes responsive to drought stress in wheat leaves. Ninety-three genes that were differentially expressed between high and low TE progeny lines were identified. One fifth of these genes were markedly responsive to drought stress. Several potential growth-related regulatory genes, which were down-regulated by drought, were expressed at a higher level in the high TE lines than the low TE lines and are potentially associated with a biomass production component of the Quarrion-derived high TE trait. Eighteen of the TE differentially expressed genes were further analysed using quantitative RT-PCR on a separate set of plant samples from those used for microarray analysis. The expression levels of 11 of the 18 genes were positively correlated with the high TE trait, measured as carbon isotope discrimination (Delta(13)C). These data indicate that some of these TE differentially expressed genes are candidates for investigating processes that underlie the high TE trait or for use as expression quantitative trait loci (eQTLs) for TE.

Carbon Isotopes↗

PubMed2026

Mitochondrial Function-Related Genes in Sleep Disorders: A Multi-Omics Mendelian Randomization Study.

Mitochondrial dysfunction is linked to sleep disorders in previous report, but the potential roles of specific genes remain unclear. This study aimed to dissect different subtype-specific genetic associations and their underlying mechanisms. A multi-omics Summary-data-based Mendelian Randomization (SMR) approach was performed to identify potential causal links between mitochondrial function-related genes and sleep disorders. We integrated GWAS data from FinnGen database (the discovery set), independent GWAS datasets (covering different sleep-disorder subtypes and used for validation), and cis-QTLs (including mQTLs, eQTLs, and pQTLs) to perform systematic exploration. Specially, we performed targeted validation of tissue-specific effects, leveraging gene expression data from disease-relevant brain regions within the GTEx database. Our SMR analysis identified mitochondrial function-related genes potentially modulating sleep disorders across biological layers, initially identifying 102 genes at the methylation level, 48 at the gene expression level, and 6 at the protein abundance level. Integrative analysis subsequently prioritized DCXR and ACADVL and revealed their distinct, subtype-specific associations. DCXR exhibited a protective role in sleep apnea while ACADVL showed a paradoxical risk conferring role in daytime sleepiness. In addition, the analysis identified an epigenetic regulatory mechanism for DCXR in which its expression and protein levels are modulated by DNA methylation. Finally, validation in brain-hypothalamus tissue confirmed DCXR as a significant potential protective factor (OR = 0.929, 95% CI: 0.887-0.973, P_HEIDI = 0.999, FDR = 0.2449). Our findings implicate key mitochondrial genes, particularly DCXR and ACADVL, in the pathophysiology of specific sleep disorder subtypes, highlighting potential avenues for precision medicine. Clinical trial number: Not applicable.

Humans↗

PubMed2026

AKR7A3 rs1738023 association with susceptibility to female hepatocellular carcinoma and its role in AFB1 metabolism and tumor.

BACKGROUND: Hepatocellular carcinoma (HCC) is one of the most common cancer worldwide. In this study, we performed a two-stage exome-chip association analysis and found that the aldo-keto reductase family7 member A3 (AKR7A3) rs1738023 may be a potential susceptibility locus for HCC in females. We aimed to explore its role and mechanism. METHODS: The association between genotype and phenotype was analyzed through GWAS method. The expression of AKR7A3 in cancer tissue and blood analysis by qRT-PCR. The relationship of AKR7A3 and aflatoxin B1 (AFB1) was also analyzed. The effect of AKR7A3 on the biological behavior of HCC cell line was investigated on proliferation and invasion. The potential mechanism was analyzed by transcriptome analysis and western blot. RESULTS: Through genome-wide association analysis (GWAS), AKR7A3 (rs1738023), KIF2C (rs4342887), and CYP3A5 (rs6977165 and rs4646450) were found to be associated with susceptibility to hepatocellular carcinoma (HCC) in women. Further expression quantitative trait loci (eQTL) analysis showed that only AKR7A3 (rs1738023) was significantly associated with gene expression. The expression of AKR7A3 was significantly lower in HCC than adjacent non-tumorous tissues (P < 0.001). The genotype of rs1738023 was significantly associated with AKR7A3 expression (P = 0.0085). Rs1738023[C] genotype had a low AKR7A3 expression level and limited detoxification ability of AFB1. Literature data showed that AKR7A3 is involved in the metabolism of aflatoxin B1 (AFB1). Functional experimental results showed that overexpression of AKR7A3 in the normal liver cell line HL-7702 could significantly reduce AFB1-induced ROS levels and DNA adduct formation, suggesting that it plays a protective role in AFB1 metabolic detoxification. Cell function test showed that overexpression of AKR7A3 inhibit the proliferation, migration and invasion of HCC cells, and block the cell cycle. Transcriptome sequencing and KEGG pathway enrichment analysis revealed that overexpression of AKR7A3 affected the PI3K signaling pathway and led to downregulation of HIF1A and its downstream VEGFA protein expression. The validation results were confirmed in HCC cell lines Huh-7 and SUN-387. CONCLUSION: Overexpression of AKR7A3 contributes to inhibition of HCC progression and reduction of aflatoxin toxicity. AKR7A3 may serve as a potential prognostic and therapeutic target for HCC patients, although further validation is needed.

AKR7A3↗

PubMed2026

Unraveling causal links between chronic rhinosinusitis and peripheral artery diseases: insights from genetic correlations through genome-wide association studies.

OBJECTIVES: Chronic Rhinosinusitis (CRS) shares epidemiological links with Cardiovascular Diseases (CVDs), however, their shared genetic basis remains unclear. We hypothesized that pleiotropic genetic variants underlie CRS-CVDs links via distinct biological pathways. METHODS: Using large-scale GWAS data from European-ancestry individuals, we assessed global and local genetic correlations. We applied Genomic Structural Equation Modeling (Genomic SEM) to dissect shared genetic architecture, performed bidirectional Mendelian Randomization (MR) to infer causality, and conducted cis-eQTL colocalization to identify shared genetic signals. Finally, in vitro endothelial models (HUVECs) validated the functional dynamics of candidate genes under CRS-mimicking inflammatory stress. RESULTS: CRS showed significant genetic correlations with multiple CVDs. Genomic SEM revealed a latent factor structuring shared genetic risk through three pathways: artery diseases, myocardial diseases, and heart failure. Local genetic correlations identified significant local genetic correlations specifically between CRS and Peripheral Atherosclerosis (PAS)/Peripheral Artery Disease (PAD) specifically within the chr6: 31.57‒33.24 Mb locus. MR demonstrated causal effects of CRS on PAD (OR = 1.23, p = 0.022) and PAS (OR = 1.21, p = 0.011), but not vice versa. Genetically predicted HLA-DRB1, APOM, and COL11A2 expression conferred protection, while HLA-DQA2 increased risk. Crucially, in vitro validation corroborated these pathogenic trajectories, inflammatory stress significantly downregulated the protective APOM and upregulated the risk-associated HLA-DQA2 alongside pro-atherogenic VCAM-1, while HLA-DRB1 exhibited a compensatory upregulation (p < 0.05). CONCLUSION: CRS shares global genetic liability with CVDs, structured through three primary etiological pathways. Causal effects of CRS on peripheral artery diseases are mediated by immune and lipid-related genes within the chr6 locus, revealing divergent pleiotropic mechanisms. Our integrated genetic and in vitro evidence provides a mechanistic framework wherein chronic mucosal inflammation contributes to systemic endothelial vulnerability, thereby highlighting candidate targets for mechanism-directed therapy.

Humans↗

PubMed2024

Genetic variants reduced POPs-related colorectal cancer risk via altering miRNA binding affinity and m6A modification.

Exposure to persistent organic pollutants (POPs) may contribute to colorectal cancer risk, but the underlying mechanisms of crucial POPs exposure remain unclear. Hence, we systematically investigated the associations among POPs exposure, genetics and epigenetics and their effects on colorectal cancer. A case-control study was conducted in the Chinese population for detecting POPs levels. We measured the concentrations of 24 POPs in the plasma using gas chromatography-tandem mass spectrometry (GC-MS/MS) and evaluated the clinical significance of POPs by calculating the area under the receiver operating characteristic curve (AUC). To assess the associations between candidate genetic variants and colorectal cancer risk, unconditional logistic regression was used. Compared with healthy control individuals, individuals with colorectal cancer exhibited higher concentrations of the majority of POPs. Exposure to PCB153 was positively associated with colorectal cancer risk, and PCB153 demonstrated superior accuracy (AUC=0.72) for predicting colorectal cancer compared to other analytes. On PCB153-related genes, the rs67734009 C allele was significantly associated with reduced colorectal cancer risk and lower plasma levels of PCB153. Moreover, rs67734009 exhibited an expression quantitative trait locus (eQTL) effect on ESR1, of which the expression level was negatively related to PCB153 concentration. Mechanistically, the risk allele of rs67734009 increased ESR1 expression via miR-3492 binding and m6A modification. Collectively, this study sheds light on potential genetic and epigenetic mechanisms linking PCB153 exposure and colorectal cancer risk, thereby providing insight into the accurate protection against POPs exposure.

Humans↗

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

Genetic evidence for a causal relationship between melatonin metabolism and depression.

To investigate the causal relevance of melatonin metabolism, which provides the biological basis for circulating melatonin levels, to specific depression symptom subtypes, we performed a targeted systematic review of melatonin metabolism pathways in the human brain and liver. Using two-sample Mendelian randomization (MR), we assessed the causal effects of metabolism pathways and/or individual genes on major depressive disorder (MDD) and nine symptom subtypes derived from Patient Health Questionnaire-9 (PHQ-9). Instrumental variables (IVs) were expression quantitative trait loci (eQTL) for eight individual genes, one synthesis route, and three degradation routes. Results were assessed using Bayesian colocalization and phenome-wide association analyses. At the pathway-level, the genetically proxied synthesis-route signal was associated with PHQ-9 Assessment 5 (PHQ9A5, OR: 0.89, 95% CI: 0.85-0.93), but sensitivity analyses suggested this association was primarily driven by TPH1 and may reflect serotonin-related biology. In contrast, higher brain melatonin degradation raised the risk of both PHQ9A1 (OR: 1.03, 95% CI: 1.02-1.04) and PHQ9A7 (OR: 1.03, 95% CI: 1.02-1.03). Within degradation, up-regulation of the kynurenine sub-pathway increased the odds of PHQ9A3 (OR: 1.05, 95% CI: 1.02-1.07), PHQ9A4 (OR = 1.04, 95% CI: 1.02-1.06) and PHQ9A7 (OR: 1.05, 95% CI: 1.02-1.07). Gene-level analyses were largely concordant, except for SULT1A1, whose higher expression was genetically protective for PHQ9A3 but risk-increased for PHQ9A1 and PHQ9A4. Overall, these results demonstrate that melatonin metabolism exerts symptom-specific and pathway-specific causal effects on depression. A stratified view of melatonin's role may help optimize the application of exogenous melatonin supplementation.

Melatonin↗

PubMed2026

The bioinformatics approach to identifying pathogenic variants for colorectal cancer (CRC).

Colorectal cancer (CRC) is the third most prevalent cancer globally, accounting for 9.6% of newly diagnosed cases and 9.3% of cancer-related deaths. It develops from the uncontrolled proliferation of glandular cells in the colon and rectum and is categorized into three primary types: sporadic, hereditary, and colitis-associated. While genetic susceptibility is a key factor in CRC pathogenesis, identifying high-impact pathogenic variants remains a significant challenge. This study integrates bioinformatics and population genetics approaches to identify CRC-associated single-nucleotide polymorphisms (SNPs) with potential clinical significance. CRC-associated SNPs were extracted from the Genome-Wide Association Studies (GWAS) Catalog, functionally annotated via HaploReg, and validated via Ensembl. In addition, expression quantitative trait locus (eQTL) data from the GTEx database were used to assess the effects of these variants on gene expression across human tissues. Our analysis identified three high-priority SNPs (rs9379084, rs3184504, and rs11557154) associated with the RREB1, ATXN2, SH2B3, and DCAF12 genes, which exhibited marked allele frequency differences among populations. These findings suggest potential biomarkers for CRC risk assessment and highlight the importance of genetic screening across diverse populations.

Bioinformatics↗

PubMed2026

Integrative post-GWAS analysis prioritizes immune regulatory pathways and candidate effector signals in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) has a complex polygenic architecture, but translating genome-wide association signals into biologically interpretable candidates remains challenging. We applied an integrative post-GWAS framework to refine SLE-associated loci and prioritize candidate regulatory mechanisms. METHODS: European-ancestry SLE GWAS summary statistics from FinnGen and Bentham et al. were meta-analysed, comprising 8417 cases and 354,277 controls. After quality filtering, 6,782,131 SNPs were retained. Downstream analyses included LAVA regional prioritization, Bayesian colocalization with GTEx v8 whole-blood and spleen eQTLs, independent replication in the Julià et al. Spanish cohort, pathway enrichment, bivariate LAVA cross-trait local genetic correlation, and therapeutic annotation. RESULTS: The discovery meta-analysis identified 46 genome-wide significant SLE-associated loci, including putative novel signals requiring database/literature qualification. LAVA identified 14 candidate index variants across 12 high-confidence regions, of which nine index variants were retained as the primary prioritized set based on LAVA support and/or convergent regulatory evidence. The strongest association mapped to the chr6p21.3/MHC region (rs389884), where four genes showed colocalization support, including CLIC1 in whole blood and C4A in spleen. Because the chr6p21.3/MHC rs389884 region lead variant was unavailable for replication and no suitable proxy was identified, this signal was interpreted as an emerging candidate for functional validation rather than a replicated causal signal. Seven available variants replicated with concordant effects. An exploratory Roadmap immune chromatin-state overlap analysis placed 15 of 45 non-MHC lead variants (33.3%) directly, and 34 of 45 (75.6%) within ±10 kb, in active immune enhancer/promoter states. Pathway analyses highlighted type I interferon, JAK-STAT signaling, cytokine regulation, and antigen presentation, while bivariate LAVA analyses supported shared local genetic architecture with rheumatoid arthritis, systemic sclerosis, and Sjögren syndrome. CONCLUSIONS: This integrative post-GWAS analysis refines SLE association signals into biologically interpretable candidate regions and supports interferon and JAK-STAT signaling as central genetically supported pathways in SLE.

CLIC1↗

PubMed2026

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

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

Bayesian colocalization↗

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

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↗

PubMed2005

Genetical genomics in humans and model organisms.

Genetical genomics has been proposed to map loci controlling gene-expression differences (eQTLs) that might underlie functional trait variation. We briefly review the studies in model species and conclude that, although they successfully demonstrate the utility of genetical genomics, they are too limited to unlock the full potential of this approach and some results should be interpreted with caution. We subsequently elaborate on two recent studies that use this approach in humans. The many differences between these studies complicate meaningful comparisons between them. A joint analysis of the two experiments offers some scope for more powerful genetical genomics.

Animals↗