SEARCH · PubMed Health

Results for “eQTL analysis”

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

Identifying causal genetic variants for high-altitude adaptation through blood eQTL analysis in plateau populations.

A substantial number of genetic variants have been associated with high-altitude adaptation (HAA), yet most of them are located in non-coding genomic regions, leaving their specific functions and underlying mechanisms largely unknown. In this study, we analyze whole-genome and transcriptome sequencing data from a self-established cohort comprising 61 native highlanders (NHs) and 164 acclimatized newcomers (ANs), identifying 6,586 cis- and 34,203 trans-expression quantitative trait loci (eQTLs), along with 130 cell type-specific eQTLs. By further combining these data with a large East Asia (~30% Tibetan) genome-wide association study (GWAS) cohort, we employ colocalization and causal inference analyses to prioritize 85 cis-eQTLs associated with HAA and identify several novel candidate causal genes, including EXOC8, which is experimentally confirmed to regulate erythroid differentiation. Additionally, network analysis of these causal genes uncovers multiple regulatory pathways, mainly involving energy metabolism, autophagy, ubiquitination and inflammation. Our study offers a comprehensive eQTL map and reveals causal chains of "variant-gene-phenotype" for HAA-related traits, which provides new insights into potential regulatory mechanisms and targets for prevention and treatment of altitude sickness.

Quantitative Trait Loci↗

PubMed2026

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

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

Journal Article↗

PubMed2026

Multi-trait GWAS identifies pleiotropic loci shared between early pregnancy bleeding and psychiatric traits.

INTRODUCTION: Early pregnancy bleeding is a common pregnancy complication, yet its genetic basis and potential links with psychiatric traits remain poorly understood. This study aimed to characterize the shared genetic architecture between early pregnancy bleeding and reproductive, psychiatric, and cardiometabolic traits. METHODS: We integrated linkage disequilibrium score regression (LDSC), local genetic correlation analysis (LAVA), and multi-trait genome-wide association analysis (MTAG). LDSC was used to estimate genome-wide genetic correlations, including sex-stratified analyses. LAVA was applied to identify genomic regions contributing to local genetic sharing. Guided by these correlation patterns, MTAG was performed to improve locus discovery, followed by cis-eQTL analysis using GTEx v8 to explore potential regulatory mechanisms. RESULTS: LDSC revealed significant positive genetic correlations between early pregnancy bleeding and reproductive traits, including endometriosis, miscarriage, and uterine fibroids. Strong positive correlations were also observed with several psychiatric disorders, including major depressive disorder, post-traumatic stress disorder, and attention deficit hyperactivity disorder. Sex-stratified analyses suggested stronger genetic correlations with emotional reactivity-related traits in females, whereas social and behavioral traits were more prominent in males. LAVA localized these shared signals to specific genomic regions and identified pleiotropic hotspots at 8q21 near RUNX1T1 and 9p21 near CDKN2A/B. MTAG identified two novel loci, 15q15.1 marked by rs45457497 and 11q13.1 marked by rs2452681. Cis-eQTL analysis showed that the lead variant at 15q15.1 regulates RMDN3 expression across multiple brain regions, while the 11q13.1 locus regulates PACS1, GAL3ST3, and SF3B2 expression in brain tissues and the pituitary. DISCUSSION: These findings position early pregnancy bleeding as a multifactorial trait shaped by shared reproductive, psychiatric, neuroendocrine, and stress-related biology. The implication of RMDN3, which encodes a mitochondrial outer membrane protein involved in ER-mitochondria tethering and calcium homeostasis, suggests a potential molecular link between neuroendocrine stress pathways, psychiatric susceptibility, and reproductive vulnerability.

RMDN3↗

PubMed2026

Tonsillar expression quantitative trait loci verify and expand genetic contributors to childhood atopic diseases.

BACKGROUND: The spectrum of causal variants, mechanisms, and immunologic gene networks that influence pediatric atopic traits is not completely understood. Human genetic variation associated with transcript abundance (expression quantitative trait loci [eQTLs]) can help to advance our understanding, yet prior work has focused on profiling immune cell populations collected from peripheral blood primarily in adult populations, leaving tissue-resident lymphocytes collected from children uncharacterized. OBJECTIVE: We sought to characterize gene expression of 4 populations of tonsil-derived immune cell types collected from pediatric patients. METHODS: We collected naive B, germinal center B, naive T, and T follicular helper cells from the discarded tonsils of 103 children across development (age range 1-19). Following genotyping and RNA sequencing of samples, we performed differential expression and eQTL analysis, then statistically linked eQTL signals to relevant atopic traits via colocalization. RESULTS: We found differentially expressed genes across cell types and identified 13,393 expression genes (eGenes) (1,793 eGenes not previously reported in similar datasets) influenced by 27,603 eQTLs (5,199 eQTLs not previously reported). We linked eQTLs to associations identified in pediatric and adult asthma and atopy traits, nominating 78 eGenes including TRAF3, ZBTB10, and JAZF1 in disease-relevant cell types. CONCLUSIONS: Our freely available resource exemplifies the importance of discovery in native tissues and across human development.

Expression quantitative trait locus↗

PubMed2026

Identifying Single-Cell Expression Quantitative Trait Loci Using a Bootstrap Penalized Hurdle Model.

BACKGROUND: Expression quantitative trait loci (eQTL) analysis links genetic variants to gene expression levels, helping to uncover how genetic variation contributes to gene regulation. While traditional eQTL analyses rely on bulk RNA-seq data, recent advances in single-cell RNA sequencing (scRNA-seq) have made it possible to detect cell-type-specific eQTLs. However, the inherent sparsity and heterogeneity of scRNA-seq data present major challenges for standard modeling approaches. METHODS: In this paper, we propose a novel statistical framework, Bootstrap Penalized Hurdle regression model (BPHurdle), designed specifically for scRNA-seq data. BPHurdle employs a hurdle modeling framework, where a logistic component accounts for the excess zeros in single-cell expression data, and a Poisson component jointly evaluates the effects of multiple SNPs on positive gene expression levels. RESULTS: Through simulation studies, we show that BPHurdle achieves high accuracy and robustness in identifying regulatory variants. We further demonstrate its utility on a real dataset through a case study focusing on a subset of differentially expressed genes, where it successfully identifies reliable cell-type-specific eQTLs. CONCLUSIONS: Overall, BPHurdle offers an advanced and flexible approach for single-cell eQTL mapping, providing deeper insight into the genetic regulation of gene expression at cellular resolution.

Quantitative Trait Loci↗

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

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

Integrative genomic and transcriptomic analysis of hypertension in a Taiwanese population.

OBJECTIVES: Hypertension is highly prevalent in Asian populations and represents a major cardiovascular risk factor. However, most genome-wide association studies (GWASs) and transcriptome-wide association studies (TWASs) have focused primarily on Caucasian cohorts. This study aimed to identify genetic loci and gene expression signatures associated with hypertension in an Asian population. METHODS: We analyzed 10 739 hypertensive patients and 49 668 controls from the Taiwan Biobank, testing 4 512 191 genome-wide single nucleotide polymorphisms (SNPs). Integrated GWAS, TWAS, and expression quantitative trait locus (eQTL) analyses were conducted to characterize genetic risk. Additionally, a polygenic risk score (PRS) was constructed using a split-sample design to evaluate genetic risk stratification. RESULTS: We identified 14 loci significantly associated with hypertension, including a novel locus at 5p13.1. eQTL analysis linked this locus to DAB2 expression in whole blood. TWAS detected 55 hypertension-associated genes, with 20 (36%) overlapping GWAS loci. Several novel genes outside GWAS loci, including FBXL15, KCNIP2, and CRIP3, were highly significant and implicated in vascular biology and hypertension mechanisms. PRS analysis effectively differentiated hypertension risk, with individuals in the top 10% showing a > 3.5-fold increased risk compared to the bottom 10%. CONCLUSIONS: Our findings provide new insights into the genetic and transcriptomic landscape of hypertension in Asians. The identification of novel loci and genes advances understanding of disease biology and may guide precision medicine approaches for risk prediction and therapeutic development.

Female↗

PubMed2025

New insights into genetic comorbidity mechanisms: type 2 diabetes and primary open-angle glaucoma.

AIMS: To investigate the shared genetic mechanisms between type 2 diabetes (T2D) and primary open-angle glaucoma (POAG). Using large-scale genome-wide association study (GWAS) data, we performed single nucleotide polymorphism (SNP) level analysis to detect pleiotropic variants and loci, paired eQTL mapping analysis and gene-level analysis to identify candidate pleiotropic genes. In addition, Mendelian randomisation (MR) analysis was performed to assess causal associations. MATERIALS AND METHODS: We used POAG GWAS data from Finngen (9565 cases and 430 250 controls) and T2D GWAS data from 55 555 European ancestry samples. We used Linkage Disequilibrium SCore (LDSC) regression to assess the genetic association between T2D and POAG and further used PLeiotropic Analysis under the COmposite null hypothesis (PLACO) to identify shared genetic variants between paired traits. Finally, we further used MR analysis to explore the causal association between T2D and POAG at the genetic level. RESULTS: The LDSC results and MR analysis revealed that the T2D effect was significantly higher than that of the POAG (OR=1.09, 95% CI 1.03 to 1.14, p=1.50×10-3). The PLACO property analysis determined that the T2D sum POAG shared 178 individual SNPs, separate localisation of 79 individual causes. The five most popular choices are based on the effectiveness of CCND2, SVEP1, ST6GAL1, TCF7L2 and HMGA2. expression quantitative trait loci mapping further revealed 36 genes with regulatory roles in optic nerve-related brain tissues. Functional enrichment analyses indicated that these pleiotropic genes are involved in neurodevelopmental, neuroprotective and metabolic pathways, with tissue-specific enrichment observed in neural, pancreatic, adipose and retinal tissues. It is possible to present the main comorbid mechanisms of T2D and POAG. CONCLUSIONS: Our study provides new insights into the aetiology and pathogenesis of T2D and POAG at the genetic level.

Humans↗

PubMed2026

Mitochondria-Related Pathogenic Genes in Paediatric Asthma: A Multi-Omics Mendelian Randomization Study.

Mitochondrial dysfunction is implicated in asthma pathogenesis, but causal roles of mitochondrial-related genes in paediatric asthma remain unclear. We performed a multi-omics Mendelian randomization study integrating GWAS data from paediatric asthma cohorts with blood-based methylation quantitative trait loci (mQTLs), expression QTLs (eQTLs) and protein QTLs (pQTLs) datasets. Causal inference was assessed using Summary-data-based Mendelian Randomization (SMR) and HEIDI testing, complemented by colocalization analysis. Findings were validated in independent cohorts and evaluated for tissue specificity using GTEx. Functional enrichment and protein-protein interaction (PPI) network analyses were conducted. SMR analysis identified 80 methylation sites spanning 54 genes, 26 gene expressions, and three proteins significantly associated with paediatric asthma. Colocalization analysis confirmed strong evidence for 10 methylation sites (7 genes), the STX17 eQTL (PP.H4 = 0.98) and the UNG pQTL (PP.H4 = 0.84). Tissue-specific eQTL validation replicated the STX17 association. Multi-omics integration associated ALAS1 (cg13241645, cg15698299) and TXNRD1 (cg09884423) with asthma at both methylation and expression levels, with colocalization supporting both ALAS1 associations. Furthermore, integrated mQTL-eQTL analysis suggests that DNA methylation potentially regulates ALAS1 and TXNRD1 expression. Functional enrichment and network analyses revealed that these candidate genes converge on mitochondrial metabolic pathways and identified seven hub genes with potential regulatory significance (SDHB, MFN2, GLDC, PHB2, TXNRD1, ATP5MC1 and PHB). This study provides multi-omics evidence supporting a causal role for mitochondrial-related genes, particularly ALAS1 and TXNRD1, in paediatric asthma, offering new insights into pathogenesis and potential therapeutic targets.

Humans↗

PubMed2026

Integrated bioinformatics analyses for GSDMB in carcinogenesis and progression of bladder cancer.

BACKGROUND: Emerging evidence suggests that pyroptosis influences the development of various diseases. Gasdermin B (GSDMB), an intracellular protein that executes pyroptosis, has recently attracted attention for its potential role in tumor biology. However, its specific function in bladder cancer (BLCA) remains unclear. Therefore, this study aimed to investigate the potential role of GSDMB in the carcinogenesis and prognosis of BLCA patients. METHODS: Mendelian randomization (MR) studies were conducted to examine relationships between the expression of GSDMB and BLCA with expression quantitative trait loci (eQTL) data. Then, GSDMB mRNA expression data and clinical characteristics of BLCA patients were retrieved from The Cancer Genome Atlas (TCGA) database. Cox regression was used to explore the relationship between GSDMB mRNA expression and patients' survival. Additionally, the correlation between GSDMB and the immune microenvironment, tumor mutational burden (TMB), tumor microenvironment (TME), and drug sensitivity in BLCA was examined. RESULTS: According to MR analysis based on eQTLs, GSDMB mRNA expression has positive causal effects on bladder carcinogenesis and the need for bladder surgery (P<0.05). The analyses of TCGA demonstrated an increased expression of GSDMB in BLCA tissues, correlating with improved patient survival. Additionally, elevated GSDMB mRNA expression was identified as an independent protective prognostic factor for BLCA, and it was associated with immune cell infiltration, TMB, TME score, and drug sensitivity. CONCLUSIONS: Elevated mRNA expression of GSDMB has a causal link to a higher risk of BLCA and the likelihood of bladder surgery, but also indicates a better prognosis. Thus, GSDMB exhibits dual effects and might serve as a potential biomarker for predicting onset and progression of BLCA. Nevertheless, further investigation of pathogenesis and mechanisms underlying GSDMB is warranted.

Bladder cancer (BLCA)↗

PubMed2026

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

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

Humans↗

PubMed2025

Association of genetically proxied cancer-targeted drugs with cardiovascular diseases through Mendelian randomization analysis.

BACKGROUND: Cancer-targeted therapies are progressively pivotal in oncological care. Observational studies underscore the emergence of cancer therapy-related cardiovascular toxicity (CTR-CVT), impacting patient outcomes. We aimed to investigate the causal relationship between different types of cancer-targeted therapies and cardiovascular disease (CVD) outcomes through a two-sample Mendelian randomization (MR) study. METHODS: This genome-wide association study was conducted using a two-sample Mendelian randomization framework. Genetic instruments for drug target gene expression were extracted from the eQTLGen consortium (31684 individuals, 37 cohorts). Genome-wide association study (GWAS) summary statistics for 19 cardiovascular diseases were derived from the FinnGen database. Primary analysis was carried out using the summary-data-based MR (SMR) method, with sensitivity analysis for validation. Colocalization analysis identifies shared causal variants between exposure eQTLs and CVD-associated single-nucleotide polymorphisms (SNPs). RESULTS: Among the 39 drug target genes, 8 were identified with detectable cis-eQTLs and were subsequently validated through positive control analysis for further investigation. In the SMR and sensitivity analyses, genetically proxied VEGFA inhibition showed significantly strong association with stroke (odds ratio [OR] = 1.17, 95% confidence interval [CI] = 1.09-1.26, p = 1.33 × 10- 5). Additionally, the inhibition of FGFR1, FLT1, and MAP2K2 exhibited suggestive association with corresponding cardiovascular disease outcomes. Nevertheless, only VEGFA expression and stroke shared a causal variant (93.6%), whereas FGFR1, MAP2K2, and FLT1 did not share causal variants with corresponding cardiovascular diseases in the colocalization analysis. CONCLUSIONS: This genetic association study revealed evidence supporting the genetic association between the use of VEGFA inhibitors and increased stroke risk, highlighting the need for enhanced pharmacovigilance. These findings underscore the delicate balance between cardiovascular toxicity risk and the benefits of cancer-targeted therapy.

Humans↗

PubMed2025

Possible linking and treatment between Parkinson's disease and inflammatory bowel disease: a study of Mendelian randomization based on gut-brain axis.

BACKGROUND: Mounting evidence suggests that Parkinson's disease (PD) and inflammatory bowel disease (IBD) are closely associated and becoming global health burdens. However, the causal relationships and common pathogeneses between them are uncertain. Furthermore, they are uncurable. Thus, we aimed to identify the causal relationships and novel therapeutic targets shared between them based on their common pathophysiological mechanisms in gut-brain-axis (GBA). METHODS: A meta-analysis on bidirectional Mendelian randomization (MR) utilizing various datasets was performed to estimate their causal relationship. Then, pleiotropic analysis under the composite null hypothesis (PLACO) with functional mapping combined with annotation of genetic associations (FUMA) analysis were conducted to identify pleiotropic genes. Next, blood, brain and intestine expression quantitative trait locus (eQTL) were taken to perform drug-target MR finding common causal genes in two diseases. Colocalization analysis ensured the eQTLs of corresponding gene colocalized with disease. Enrichment analysis and protein‒protein interaction (PPI) network were done to explore common pathogenesis pathways. Genes passed all analysis were regarded as drug targets. RESULTS: Our MR meta-analysis revealed the bidirectional causal relationship between diseases, with combined ORs for PD on IBD, CD, UC (1.050 [95% CI 1.014-1.086], 1.044 [95% CI 0.995-1.095], 1.063 [95% CI 1.016-1.120]); for IBD, CD, UC on PD (1.003 [95% CI 0.973-1.034], 1.035 [95% CI 1.004-1.067], 1.008 [95% CI 0.977-1.040]). Overall, 277, 216 and 201 genes were identified as pleiotropic genes between PD and IBD, CD, UC. Total of 733 genes were classified as tier 3 (found in only one tissue) druggable targets, 57 as tier 2 (found in two tissues, 51 protein-coding genes) and 9 as tier 3 (found in three tissues). Among 60 protein-coding druggable targets over tier 2, 18 overlapped with pleiotropic genes and enriched in mitochondria, antigen presentation, processing and immune cell regulation pathways. Three druggable genes (LRRK2, RAB29 and HLA-DQA2) passed colocalization analysis. LRRK2 and RAB29 were reported to be pleiotropic genes, and RAB29 and HLA-DQA2 were reported for the first time as potential drug targets. CONCLUSIONS: This study established a reliable causal relationship, possible shared drug targets and common pathogenesis pathways of two diseases, which had important implications for intervention and treatment of two diseases simultaneously.

Humans↗

PubMed2026

Genome-wide cis-expression Quantitative Trait Loci (eQTL) and transcriptomic signals reveal distinct molecular regulation across correlated feed efficiency traits.

INTRODUCTION: Feed efficiency (FE) is a complex trait which determines livestock production profitability, yet the molecular mechanisms behind it remain unclear. This study investigated the blood transcriptomic profile of lambs, alongside genotype data with the aim to uncover the genetic basis of FE traits such as absolute dry matter intake (DMIabsolute), DMI adjusted for body size (DMIadjusted), average daily live weight gain (ADG), and residual feed intake (RFI). MATERIALS AND METHODS: Bulk RNA-Seq and genotype data were analysed using three complementary approaches: differential gene expression (DGE) analysis, weighted gene co-expression network analysis (WGCNA), and cis-expression Quantitative Trait Loci (cis-eQTL) mapping. These methods were used independently to identify genes and regulatory networks associated with FE traits and to investigate evidence supporting multi-trait candidate gene selection. RESULTS: DGE analysis revealed 2, 24, 85 and 4 differentially expressed genes for DMIabsolute, DMIadjusted, ADG, and RFI (Padjusted < 0.05), functionally enriched in sensory perception, ATP-dependent chromatin remodeling, Notch signaling and immune response pathways. 9 gene modules significantly associated with the FE traits (P ≤ 0.05) with correlations ranging from r = -0.56 to 0.49, were identified using WGCNA. Single nucleotide polymorphism (SNP)-level cis-eQTL analysis identified 93 eSNPs associated with 74 genes (false discovery rate (FDR) < 0.05), while permutation-derived gene level analysis identified 280 eGenes (FDR < 0.2, empirical P < 0.03). Across the three analyses, applying thresholds of DGE (Padjusted < 0.05), WGCNA (correlation, P ≤ 0.05), and cis-eQTL gene-level significance (empirical P < 0.05), multiple overlapping genes were identified including DNMT3A, KANSL1, NCOR1 for DMIadjusted, ACOX2, FANCF, CIMIP2B, LOC101115106, ARMH2, LOC132657496 for ADG, and LOC114114576 for RFI representing regulators of variations in FE. DISCUSSION: The integration of DGE, WGCNA, and cis-eQTL analyses identified key genes and regulatory mechanisms associated with variation in FE traits. These results highlight that integrated multi-trait candidate gene identification approaches can reveal key genes that lower feed intake while maintaining animal growth, supporting breeding strategies aimed at improving efficiency and long-term economic sustainability in sheep.

average daily gain (ADG)↗

PubMed2025

Integrative omics of the genetic basis for wheat WUE and drought resilience reveal the function of TaMYB7-A1.

Improving wheat drought resilience and water use efficiency (WUE) is critical for sustaining productivity under increasing water scarcity. Here, we integrate genome-wide association study (GWAS), expression quantitative trait locus (eQTL) mapping, population-transcriptome analysis, and summary-data-based mendelian randomization (SMR), followed by functional validation using indexed EMS mutants and transgenic lines, to systematically identify key WUE regulators. GWAS across water conditions in 228 accessions identifies 73 quantitative trait loci (QTLs) for WUE-traits. Transcriptome profiling of 110 diverse accessions reveals 28 drought-responsive modules. eQTL mapping uncovers 146,966 regulatory variants, including condition-specific hotspots associated with key drought-related pathways. Integrative analysis underscores 85 high-confidence candidate genes, notably TaMYB7-A1. Overexpression of TaMYB7-A1 enhances photosynthesis, WUE, root development, and grain yield under drought condition by activating TaPIP2;2-B1 (water transport), TaRD20-D1 (stomatal regulation), and TaABCB4-B1 (root growth), reflecting reduced water loss and improved physiological resilience. Our study presents a comprehensive regulatory map and robust targets for wheat drought adaptation and resilient cultivar breeding.

Triticum↗

PubMed2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations.

MOTIVATION: Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. RESULTS: We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. AVAILABILITY AND IMPLEMENTATION: It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Quantitative Trait Loci↗

PubMed2025

Association of AGER genetic variants with chronic obstructive pulmonary disease susceptibility in Southern Chinese Han populations.

OBJECTIVE: Chronic obstructive pulmonary disease (COPD) remains a leading cause of disability and mortality among elderly populations. Studies indicate that AGER plays a critical regulatory role in the pathogenesis of respiratory disorders. However, the genetic variations in AGER to COPD susceptibility remain incompletely understood. This study employs a case-control design to investigate associations between AGER genetic variants and COPD risk in the Southern Chinese Han population. METHODS: This study enrolled 270 COPD patients and 271 healthy controls. AGER single-nucleotide polymorphisms (SNPs) were analysed using the MassARRAY iPLEX platform. Logistic regression models evaluated associations between AGER polymorphisms and COPD susceptibility, with false discovery rate (FDR) correction applied to mitigate multiple testing errors. SNP-SNP interactions were investigated through multifactor dimensionality reduction (MDR) analysis. Expression quantitative trait locus (eQTL) data from the GTEx database were further analysed to assess regulatory relationships between SNPs and AGER gene expression levels. RESULTS: This study showed that rs3134941 (G allele, OR = 0.21, 95% CI = 0.10-0.41, p (FDR) = 0.001) and rs3131300 (G allele, OR = 0.32, 95% CI = 0.20-0.49, p (FDR) = 0.0001) were significantly associated with a reduced susceptibility to COPD. MDR indicated that rs3131300 was the optimal predictive model for COPD risk. Additionally, initial mechanistic investigations utilizing the GTEx database identify rs3134941 (C > G) and rs3131300 (A > G) as significant expression quantitative trait loci for AGER mRNA in cell-cultured fibroblasts and whole blood. CONCLUSION: Our study demonstrated that AGER genetic variants might play a protective role in the progression of COPD.

Aged↗