SEARCH · Search PubMed

Results for “eQTL analysis”

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 91 records · Page 5Linked to original sources

PubMed2005

Genetic regulation of gene expression during shoot development in Arabidopsis.

The genetic control of gene expression during shoot development in Arabidopsis thaliana was analyzed by combining quantitative trait loci (QTL) and microarray analysis. Using oligonucleotide array data from 30 recombinant inbred lines derived from a cross of Columbia and Landsberg erecta ecotypes, the Arabidopsis genome was scanned for marker-by-gene linkages or so-called expression QTL (eQTL). Single-feature polymorphisms (SFPs) associated with sequence disparities between ecotypes were purged from the data. SFPs may alter the hybridization efficiency between cDNAs from one ecotype with probes of another ecotype. In genome scans, five eQTL hot spots were found with significant marker-by-gene linkages. Two of the hot spots coincided with classical QTL conditioning shoot regeneration, suggesting that some of the heritable gene expression changes observed in this study are related to differences in shoot regeneration efficiency between ecotypes. Some of the most significant eQTL, particularly those at the shoot regeneration QTL sites, tended to show cis-chromosomal linkages in that the target genes were located at or near markers to which their expression was linked. However, many linkages of lesser significance showed expected "trans-effects," whereby a marker affects the expression of a target gene located elsewhere on the genome. Some of these eQTL were significantly linked to numerous genes throughout the genome, suggesting the occurrence of large groups of coregulated genes controlled by single markers.

Arabidopsis↗

PubMed2026

Genetically predicted CXCL16 expression is associated with Parkinson's disease risk and peripheral immune cell dysregulation: a two-sample mendelian randomization study.

BACKGROUND: Parkinson's disease (PD) is a progressive neurodegenerative disorder with limited disease-modifying therapies. PANoptosis, an integrated form of programmed cell death involving apoptosis, pyroptosis, and necroptosis, has been implicated in neuroinflammation-related neurodegeneration. However, the roles of PANoptosis-related genes in PD remain unclear. METHODS: We performed two-sample Mendelian randomization (MR) using cis-eQTL instruments from the eQTLGen Consortium for 30 PANoptosis-related genes, with PD GWAS data from Nalls et al. 2019 as the outcome. Instrumental variables were selected using a hierarchical strategy, with genome-wide significant cis-eQTLs as primary instruments and a relaxed threshold applied only for genes with fewer than three independent SNPs. Sensitivity analyses included MR-Egger, weighted median, MR-PRESSO, MR-RAPS, and leave-one-out analyses. SMR/HEIDI testing and two-step MR mediation using 731 peripheral immune traits were also performed. RESULTS: Genetically predicted higher CXCL16 expression was associated with increased PD risk (OR = 1.115, 95% CI 1.060-1.173, p = 2.4 × 10-5), while higher FADD expression was associated with reduced PD risk (OR = 0.861, 95% CI 0.790-0.939, p = 7.1 × 10-4). CASP1 and IFI27 were nominally significant and considered exploratory. Sensitivity analyses were directionally consistent, although MR-Egger estimates were imprecise. SMR/HEIDI supported CXCL16. Exploratory mediation analysis identified 63/66 candidate immune mediators after FDR correction. CONCLUSION: These findings provide MR-based genetic evidence linking CXCL16 expression to PD risk, with exploratory mediation through peripheral immune phenotypes. The CXCL16-immune cell-PD axis warrants further experimental validation.

Humans↗

PubMed2026

SNPannotator: automated functional annotation of genetic variants and linked proxies.

SUMMARY: Genome-wide association studies (GWASs) have identified thousands of genetic variants associated with complex traits and diseases. However, explaining the mechanisms underlying phenotypic variation remains challenging. Here, we introduce SNPannotator, an automated post-GWAS analysis software package designed to streamline the interpretation of GWAS findings. Our pipeline implements a multi-step process that identifies proxy variants in high linkage disequilibrium (LD) with associated lead variants, then queries comprehensive resources (including Ensembl, the GTEx Portal, the eQTL Catalog, and STRING DB) for genomic position, deleteriousness, regulatory annotations, clinical significance, trait associations, expression (eQTLs) and splicing quantitative trait loci (sQTLs), and functional enrichment analyses and compiles the results into user-friendly reports. This package is implemented in the R programming language and includes auxiliary functions for variant lookup and LD exploration. SNPannotator provides a practical framework for efficiently deriving biologically meaningful insights from GWAS data and for assisting researchers in prioritizing candidate variants for functional validation. AVAILABILITY AND IMPLEMENTATION: The SNPannotator package is available from the Comprehensive R Archive Network (CRAN) at https://cran.r-project.org/web/packages/SNPannotator. The development version and tutorial is available on GitHub (https://github.com/omicslaboratory/SNPannotator). The online version of the package is available at https://omicslab.org/snpannotator.

Software↗

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

Genetic dissection of cardiac iron regulation using transcriptome network analysis and systems genetics in BXD mice.

Cardiac iron homeostasis is essential for myocardial energy metabolism and contractile function, yet the genetic and molecular mechanisms governing iron levels within the heart remain poorly understood. We used a systems genetics approach to dissect the transcriptional regulation of cardiac iron homeostasis. Myocardial iron level varies substantially across BXD strains (40-112 μg/g) and is under heritable genetic control (H2 = 0.38). Elevated cardiac iron is associated with reduced ventricular mass, increased ventricular ectopy, and prolonged atrioventricular conduction in the BXD population. Weighted gene co-expression network analysis of the BXD heart transcriptome identified a co-expression module that was significantly and negatively correlated with cardiac iron levels in both young and old BXD mice and enriched for pathways related to metabolic regulation, cyclic AMP (cAMP) signaling, circadian entrainment, and cardiovascular physiology. The module showed substantial overlap with a curated cardiac iron gene set, and cross-species enrichment analysis confirmed its conservation in human cardiomyopathy differentially expressed genes (enrichment ratio = 1.49; false discovery rate [FDR] = 0.0342). Quantitative trait locus (QTL) mapping of the first principal component of the overlapping module iron genes (n = 38), corroborated by individual gene mapping, identified trans-eQTL hotspots on multiple chromosomes, implicating Fcho2, Gcc2, and Rmdn1 as candidate upstream regulators operating through sequential steps of intracellular iron trafficking. Together, these findings establish a systems-level map of cardiac iron gene regulation, identify candidate genetic regulators, and provide a molecular framework linking disruption of iron-related transcriptional networks to structural and electrical cardiac dysfunction with implications for iron-related heart diseases.

BXD mouse population↗

PubMed2026

Integrative Multi-Omics Mendelian Randomization Analysis Identifies NIT2 as a Potential Metabolic Risk Gene in Hepatocellular Carcinoma.

BACKGROUND: Metabolic pathways are crucial in hepatocellular carcinoma (HCC) pathogenesis, but causal metabolic genes remain unclear. This study used Summary data-based Mendelian Randomization (SMR) and colocalization to identify metabolism-related genetic loci influencing HCC risk. METHODS: Differentially expressed genes in hepatic malignancy phenotype versus normal tissues from TCGA and GTEx were analyzed. Metabolism-related candidates were examined via SMR and colocalization using multi-omics data: methylation (mQTL), expression (eQTL), and protein (pQTL) quantitative trait loci. RESULTS: Multi-omics integration identified NIT2 as a key metabolic regulator for HCC. The cg13016775 locus of NIT2 was associated with elevated HCC risk at gene (OR = 1.618, 95% CI: 1.199-2.182) and protein (OR = 4.432, 95% CI: 1.783-11.018) levels. Colocalization supported a shared causal variant (PPH4 > 0.6), linking NIT2 to hepatocarcinogenesis via metabolic regulation. CONCLUSIONS: This study provides multi-omics evidence for NIT2 as a potential causal gene in HCC, enhancing understanding of metabolic contributions to HCC pathogenesis and highlighting integrative genomics for uncovering causal relationships.

Carcinoma, Hepatocellular↗

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↗

PubMed2025

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

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

Drug Repositioning↗

PubMed2026

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

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

Hispanic/Latino population↗

PubMed2026

Integrative multi-omics quantitative trait loci prioritize CASP7 as a candidate protective gene for cataract.

Cataracts are the leading cause of vision loss worldwide. Despite surgery being the only effective treatment, its economic burden highlights the necessity of exploring the pathogenesis of cataracts. In this study, we analyzed 4 large-scale GWAS (genome-wide association study) datasets for cataracts and performed SMR analysis along with heterogeneity in dependent instruments (HEIDI) testing to explore the effects of methylation, expression, and protein QTLs on cataracts. We further validated shared genetic variants through COLOC analysis. Additionally, we searched datasets related to cataracts from the Gene Expression Omnibus (GEO) database for differentially expressed genes (DEGs) and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses. By integrating summary-based Mendelian randomization (SMR) results with bioinformatics findings, CASP7 showed a consistent protective-direction association with cataract risk (mQTL: OR [95% CI] = 0.959 [0.941-0.977], FDR-adjusted P = .039; eQTL: OR [95% CI] = 0.897 [0.860-0.937], FDR-adjusted P = .0046; pQTL: OR [95% CI] = 0.597 [0.483-0.738], FDR-adjusted P = .00083). GEO-based analyses provided transcriptomic support for CASP7 involvement in cataract-related lens biology. These findings prioritize CASP7 as a genetically supported candidate protective gene associated with cataract risk. Because this study is based on public summary-level and transcriptomic datasets, the results should be interpreted cautiously and require functional validation in human lens-relevant systems.

Quantitative Trait Loci↗

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↗

PubMed2006

Review of microarray experimental design strategies for genetical genomics studies.

Genetical genomics approaches provide a powerful tool for studying the genetic mechanisms governing variation in complex traits. By combining information on phenotypic traits, pedigree structure, molecular markers, and gene expression, such studies can be used for estimating heritability of mRNA transcript abundances, for mapping expression quantitative trait loci (eQTL), and for inferring regulatory gene networks. Microarray experiments, however, can be extremely costly and time consuming, which may limit sample sizes and statistical power. Thus it is crucial to optimize experimental designs by carefully choosing the subjects to be assayed, within a selective profiling approach, and by cautiously controlling systematic factors affecting the system. Also, a rigorous strategy should be used for allocating mRNA samples across assay batches, slides, and dye labeling, so that effects of interest are not confounded with nuisance factors. In this presentation, we review some selective profiling strategies for genetical genomics studies, including the selection of individuals for increased genetic dissimilarity and for a higher number of recombination events. Efficient designs for studying epistasis are also discussed, as well as experiments for inferring heritability of transcriptional levels. It is shown that solving an optimal design problem generally requires a numerical implementation and that the optimality criteria should be intimately related to the goals of the experiment, such as the estimation of additive, dominance, and interacting effects, localizing putative eQTL, or inferring genetic and environmental variance components associated with transcriptional abundances.

Animals↗

PubMed2024

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

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

Mendelian Randomization Analysis↗

PubMed2026

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↗

PubMed2006

Tissue-specific expression and regulation of sexually dimorphic genes in mice.

We report a comprehensive analysis of gene expression differences between sexes in multiple somatic tissues of 334 mice derived from an intercross between inbred mouse strains C57BL/6J and C3H/HeJ. The analysis of a large number of individuals provided the power to detect relatively small differences in expression between sexes, and the use of an intercross allowed analysis of the genetic control of sexually dimorphic gene expression. Microarray analysis of 23,574 transcripts revealed that the extent of sexual dimorphism in gene expression was much greater than previously recognized. Thus, thousands of genes showed sexual dimorphism in liver, adipose, and muscle, and hundreds of genes were sexually dimorphic in brain. These genes exhibited highly tissue-specific patterns of expression and were enriched for distinct pathways represented in the Gene Ontology database. They also showed evidence of chromosomal enrichment, not only on the sex chromosomes, but also on several autosomes. Genetic analyses provided evidence of the global regulation of subsets of the sexually dimorphic genes, as the transcript levels of a large number of these genes were controlled by several expression quantitative trait loci (eQTL) hotspots that exhibited tissue-specific control. Moreover, many tissue-specific transcription factor binding sites were found to be enriched in the sexually dimorphic genes.

Animals↗

PubMed2025

Genome-Wide Association Analyses of HPV16 and HPV18 Seropositivity Identify Susceptibility Loci for Cervical Cancer.

Infection by high-risk human papillomavirus is known to exacerbate cervical cancer development. The host immune response is crucial in disease regression. Large-scale genetic association studies for cervical cancer have identified few susceptibility variants, mainly at the human leukocyte antigen locus on chromosome 6. We hypothesized that the host immune response modifies cervical cancer risk and performed three genome-wide association analyses for HPV16, HPV18 and HPV16/18 seropositivity in 7814, 7924, and 7924 samples from the UK Biobank, followed by validation genotyping in the German Cervigen case-control series of cervical cancer and dysplasia. In GWAS analyses, we identified two loci associated with HPV16 seropositivity (6p21.32 and 15q26.2), two loci associated with HPV18 seropositivity (5q31.2 and 14q24.3), and one locus for HPV16 and/or HPV18 seropositivity (at 6p21.32). MAGMA gene-based analysis identified HLA-DQA1 and HLA-DQB1 as genome-wide significant (GWS) genes. In validation genotyping, the genome-wide significant lead variant at 6p21.32, rs9272293 associated with overall cervical disease (OR = 0.86, p = 0.004, 95% CI = 0.78-0.95, n = 3710) and HPV16 positive invasive cancer (OR = 0.73, p = 0.005, 95% CI = 0.59-0.91, n = 1431). This variant was found to be a robust eQTL for HLA-DRB1, HLA-DQB1-AS1, C4B, HLA-DRB5, HLA-DRB6, HLA-DQB1, and HLA-DPB1 in a series of cervical epithelial tissue samples. We additionally genotyped twenty-four HPV seropositivity variants below the GWS threshold out of which eleven variants were found to be associated with cervical disease in our cohort, suggesting that further seropositivity variants may determine cervical disease outcome. Our study identifies novel genomic risk loci that associate with HPV type-specific cervical cancer and dysplasia risk and provides evidence for candidate genes at one of the risk loci.

Humans↗

PubMed2026

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

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

Humans↗

PubMed2025

From genetic causality to druggable targets: A multiomics framework identifies ZSCAN16 in gout pathogenesis.

ObjectiveGout is a prevalent form of inflammatory arthritis in which many patients respond suboptimally to current therapies. Drug development is hampered by a lack of genetically validated targets, leading to high clinical trial attrition. This study aimed to systematically identify and prioritize novel, druggable targets for gout via a multilayered genetic and functional genomics approach.MethodsWe performed two-sample Mendelian randomization (MR) using cis-expression quantitative trait locus (cis-eQTL) data and dual independent gout genome-wide association study (GWAS) cohorts (openGWAS and FinnGen). The candidate genes were subjected to a rigorous validation pipeline including Bayesian colocalization, phenome-wide association studies (PheWASs) to assess pleiotropy and on-target safety, and single-cell RNA sequencing (scRNA-seq) to delineate the cellular context. Molecular docking was used to evaluate the structural druggability of prioritized targets.ResultsMR analysis revealed 15 genes causally associated with gout. Colocalization analysis (PPH4 > 0.8) prioritized two targets: ZSCAN16 (risk-increasing, OR = 1.04, 95% CI [1.02-1.06]) and TRIM10 (protective, OR = 0.96, 95% CI [0.94-0.98]). Crucially, PheWAS revealed that ZSCAN16 is highly specific to gout, whereas TRIM10 exhibited extensive pleiotropy with hematological and cardiometabolic traits, indicating significant safety risks. Single-cell analysis provided orthogonal validation, demonstrating flare-specific upregulation of ZSCAN16 in cytotoxic T/NK cells. Molecular docking confirmed ZSCAN16 as a structurally druggable target, showing high-affinity binding with known compounds (e.g. digoxin, binding energy = -9.6 kcal/mol).ConclusionsOur study identifies ZSCAN16 as a high-potential, druggable therapeutic target for gout, highlighting its genetic influence on specific immune cell activities during acute flares. Conversely, TRIM10 was deprioritized owing to substantial pleiotropic liabilities and poor chemical tractability. These findings suggest that ZSCAN16 could play a crucial role in the pathogenesis of gout and may provide a valuable lead for future drug discovery efforts.

Humans↗