SEARCH · Search PubMed

Results for “Expression quantitative trait loci”

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 55 records · Page 3Linked to original sources

PubMed2026

Immune cell-specific genetic architecture of Alzheimer's disease revealed by multi-omics analysis for therapeutic target discovery and prioritization.

Alzheimer's disease (AD) is a multifactorial neurodegenerative condition in which accumulating genetic and molecular evidence implicates dysregulation of peripheral immune processes in disease pathogenesis. Nevertheless, the contribution of distinct peripheral immune cell subsets and associated gene regulatory landscapes to AD risk remains incompletely defined. To address this gap, we integrated single-cell expression quantitative trait loci (sc‑eQTL) data from the OneK1K cohort with AD GWAS summary statistics. We systematically interrogated immune cell-specific genes for their contributions to AD risk by integrating genetic causal inference with Bayesian colocalization analyses, and identified 24 eGenes that passed both the MR significance threshold (P < 0.05) and the criterion for strong shared genetic signals (PP.H4 > 0.8). Notable candidates included GATS, HLA-DOB, HLA-DQA1, PM20D1, and others, with each gene demonstrating a cell-type-specific association restricted to its corresponding immune cell type, such as monocytes, CD8 + T cells, or B cells. Independent peripheral blood single-cell transcriptomic data further supported disease-associated shifts in cell-type-specific expression patterns in AD. Phenome-wide association studies (PheWAS) indicated limited associations with off-target traits, indicating a favorable safety profile for therapeutic intervention, with the exceptions of B4GALNT3, PM20D1, and CNN2. Integration of immune gene targets with pharmacological databases yielded three candidate compound, including NSC321521 (targeting HLA-DQA1), phenoxybenzamine (targeting GSTP1), and rimexolone (targeting BIN1). Among these compounds, Predicted blood-brain barrier permeability was observed only for phenoxybenzamine and rimexolone, with docking studies indicating stable interactions, such as those between NSC321521 and HLA-DQA1, phenoxybenzamine and GSTP1, and rimexolone and BIN1. This integrative approach highlights key immune‑cell‑specific genes involved in AD and proposes repurposable drugs with central nervous system potential, paving the way for more targeted immunomodulatory strategies in AD.

Humans↗

PubMed2026

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 association between GLP-1R expression and cardiovascular-kidney-metabolic-related diseases in non-diabetic and non-obese population: evidence triangulation using Mendelian randomization, observational and polygenic score association analysis.

BACKGROUND: Glucagon-like peptide-1 receptor (GLP-1R) agonists are emerging as promising therapies for cardiovascular-kidney-metabolic (CKM) related diseases in individuals with type 2 diabetes mellitus (T2DM) or obesity. But their effects in non-obese and non-diabetic individuals are unclear. This study triangulates evidence using Mendelian randomization (MR), polygenic scores (PGS) and observational analyses to estimate the associations of GLP-1R expression with chronic kidney disease (CKD), heart failure (HF) and metabolic dysfunction-associated steatotic liver disease (MASLD). METHODS: For the MR analysis, instruments mimicking GLP-1R expression were identified using pancreas-specific cis-expression quantitative trait loci from GTEx (N ≤ 305). MR-Robust method was used as the primary MR approach. PGS and observational analyses were performed both in non-diabetic and non-obese individuals separately. A genome-wide association study (GWAS) for MASLD (14,231 cases and 348,091 controls) was performed in the general population using data from UK Biobank. RESULTS: GLP-1R expression showed robust effects on CKD (odds ratio [OR] 0.96, 95%CI 0.95 to 0.97, q = 1.7 × 10- 10 ), HF (OR = 0.96, 95%CI 0.94 to 0.97, q = 2.5 × 10- 8) and MASLD (OR = 0.96, 95%CI 0.93 to 0.98, q = 1.3 × 10- 3) in the general population. Consistent results were observed in validation analyses. Furthermore, PGS and observational analyses among non-T2DM and non-obese individuals found little evidence to support its association with CKD, HF or MASLD. GWAS analysis identified eight conditionally independent variants associated with MASLD, in which rs563199662 was a new signal located at TFPI region. CONCLUSIONS: This study provides multilayered evidence for GLP-1R expression in mitigating CKD, HF and MASLD risks in the general population, while de-prioritized its effect on CKM-related diseases in non-obese and non-diabetic individuals. Further clinical trials are needed to validate the effects of GLP-1R agonists in relative health population.

Humans↗

PubMed2026

An encyclopedia of human enhancer-gene regulatory interactions.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the effect of human genetic variation on disease1-6. Here we create and evaluate a resource of more than 92 million enhancer-gene regulatory interactions across 1,458 biosamples covering 369 cell types and tissues, by integrating predictive models, chromatin states, three-dimensional contacts and large-scale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,356 element-gene pairs measured in CRISPR perturbation experiments, more than 30,000 fine-mapped expression quantitative trait loci and 569 fine-mapped genome-wide association study (GWAS) variants linked to a probable causal gene. Using this framework, we develop ENCODE-rE2G, a predictive model achieving state-of-the-art performance across several prediction tasks, demonstrating that iterative perturbations and supervised machine learning can build increasingly accurate predictive models of enhancer regulation. Using ENCODE-rE2G, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, revealing global properties of enhancer networks, identifying differences in regulatory complexity across genes and improving analyses linking noncoding variants to target genes and cell types for common complex diseases. By interpreting the model, we find that beyond enhancer activity and three-dimensional enhancer-promoter contacts, additional features that guide enhancer-promoter communication include promoter class and enhancer-enhancer synergy. These genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

Humans↗

PubMed2026

Exploring genomic regions regulating the liver transcriptome and energy homeostasis in pigs.

In pigs, energy homeostasis has an impact on meat quality and health. In a Duroc pig population, 30 quantitative trait locus (QTL) regions associated with fatty acid (FA) composition in adipose tissue, plasma, liver and muscle were previously identified. Mapping of expression quantitative trait locus (eQTL) regions will provide a molecular hypothesis for genotype-phenotype interactions and may allow the identification of shared causal variants, key to increasing our understanding of the genetic regulation of FA composition and energy homeostasis. However, gene expression is impacted by environmental factors, while individual-level allelic imbalance (AI) can be more reliable and can be surveyed via allelic-specific expression (ASE) analysis. Furthermore, treatment of ASE as a quantitative trait allows the identification of allele-specific expression quantitative trait loci (aseQTLs), which are variants whose heterozygosity is linked to the AI of a nearby single-nucleotide polymorphism (SNP), pointing to regulatory elements. In this study, liver was selected as a key metabolic hub with an important role in the regulation of energy homeostasis, and 310 liver RNA sequencing samples were analysed using a combination of (1) eQTL mapping, (2) ASE analysis, and (3) aseQTL mapping methods. A total of 2 188 eQTL regions were identified, mostly cis-eQTL regions (73.17%). ASE analysis reported 1 964 ASE SNPs, associated with 633 genes. Finally, aseQTL mapping reported 64 172 aseQTL, associated with the AI of 31 genes. Colocalisation analysis combined with ASE analysis showed that the expression of FADS1 and FADS2 genes is associated with the polyunsaturated FA composition in several tissues, where microRNA regulation may be present. Finally, in the DGAT2 gene, annotated in a QTL region associated with multiple FAs in adipose tissue, ASE revealed allelic imbalance in the 3' untranslated region (UTR) of this gene. Allelic differential expression can be caused by a 13-bp insertion affecting messenger RNA stability, previously described, exemplifying how allelic imbalance is caused by a post-transcriptional regulatory mechanism undetectable by eQTL mapping. Furthermore, aseQTLs were associated with this gene, linked to a previously identified copy number variant not yet associated with DGAT2 expression. These results demonstrate how ASE analysis and aseQTL mapping can complement eQTL mapping, as they resolved a complex region affected by allelic heterogeneity, a main confounding effect of QTL mapping. In conclusion, the combination of eQTL mapping, ASE analysis and aseQTL mapping allowed the characterisation of the regulation of liver gene expression, improving our understanding of the genetic determinism of energy homeostasis.

Allele-specific expression↗

PubMed2026

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

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

CEBPZ↗

PubMed2026

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↗

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

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

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

Mendelian randomization↗

PubMed2025

Cross-Phenotype Genome-Wide Association Study on the Shared Genetic Susceptibility to Systemic Sclerosis and Primary Biliary Cholangitis.

OBJECTIVE: An increased risk of primary biliary cholangitis (PBC) has been reported in patients with systemic sclerosis (SSc). Our study aims to investigate the shared genetic susceptibility between the two disorders and to define candidate causal genes using cross-phenotype genome-wide association study (GWAS) meta-analysis. METHODS: We performed cross-phenotype GWAS meta-analysis and Bayesian colocalization analysis for patients with SSc and patients with PBC. We performed both genome-wide and locus-based analysis, including tissue and pathway enrichment analyses, fine-mapping, Bayesian colocalization analyses with expression quantitative trait loci and protein quantitative trait loci (pQTL) datasets, and phenome-wide association studies. Finally, we used an integrative approach to prioritize candidate causal genes from the novel loci. RESULTS: We detected a strong genetic correlation between SSc and PBC (global genetic correlation = 0.84, P = 1.7 × 10-6). In the cross-phenotype GWAS meta-analysis, we identified 44 nonhuman leukocyte antigens loci that reached genome-wide significance (P < 5 × 10-8). Evidence of shared causal variants between patients with SSc and patients with PBC was found for nine loci, five of which were novel. Integrating multiple sources of evidence, we prioritized CD40, ERAP1, PLD4, SPPL3, and CCDC113 as novel candidate causal genes. The CD40 risk locus colocalized with trans-pQTLs of multiple plasma proteins involved in B cell function. CONCLUSION: Our study supports a strong shared genetic susceptibility between SSc and PBC. Using cross-phenotype analyses, we have prioritized several novel candidate causal genes and pathways for these disorders.

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↗

PubMed2024

Dissecting the association between blood pressure traits, hypertension, antihypertensive medications and epilepsy: A Mendelian randomization study.

BACKGROUND: Observational studies suggest that hypertension and epilepsy have a high co-occurrence, and antihypertensive medications may have impacts on the prevention and treatment of epilepsy. However, the directionality of causation between them is elusive. METHOD: By leveraging genome-wide association studies (GWAS) summary data of each trait, we firstly performed bidirectional univariate Mendelian randomization (UVMR) to assess the strength and direction of the associations between pairs of traits, then multivariate MR (MVMR) was conducted to adjust for potential confounders in causalities. Cochran's Q statistics, leave-one-out analysis, MR-Egger regression and MR-Pleiotropy Residual Sum and Outlier methods (MR-PRESSO) were employed to evaluate the robustness of the results. Drug target MR was proceeded to assess the association between five classes of first-line antihypertensive medications and epilepsy. Specifically, single nucleotide polymorphisms (SNPs) extracted from GWAS data on systolic blood pressure (SBP)/diastolic blood pressure (DBP), along with expression quantitative trait loci (eQTL) were utilized as proxies for antihypertensive medications, respectively. RESULTS: Forward UVMR results provided evidence that genetically predicted blood pressure traits and hypertension have causal effects on epilepsy, while reverse UVMR indicated no causal impacts of epilepsy on blood pressure traits or hypertension. The sensitivity analysis results were robust. The causalities between DBP, hypertension and epilepsy remained remarkable after adjustment by MVMR. Inverse-variance-weighted MR (IVW-MR) yielded evidence of positive association only between Beta-Blockers target genes based on DBP GWAS screening and epilepsy. Summary-data-based MR (SMR) identified a positive correlation between Beta-Blockers target gene ADRA1D and epilepsy risk. CONCLUSIONS: Hypertension has a causal effect on epilepsy and managing DBP in patients with hypertension through Beta-Blockers may help prevent epilepsy.

Humans↗

PubMed2017

Functional mapping and annotation of genetic associations with FUMA.

A main challenge in genome-wide association studies (GWAS) is to pinpoint possible causal variants. Results from GWAS typically do not directly translate into causal variants because the majority of hits are in non-coding or intergenic regions, and the presence of linkage disequilibrium leads to effects being statistically spread out across multiple variants. Post-GWAS annotation facilitates the selection of most likely causal variant(s). Multiple resources are available for post-GWAS annotation, yet these can be time consuming and do not provide integrated visual aids for data interpretation. We, therefore, develop FUMA: an integrative web-based platform using information from multiple biological resources to facilitate functional annotation of GWAS results, gene prioritization and interactive visualization. FUMA accommodates positional, expression quantitative trait loci (eQTL) and chromatin interaction mappings, and provides gene-based, pathway and tissue enrichment results. FUMA results directly aid in generating hypotheses that are testable in functional experiments aimed at proving causal relations.

Chromatin↗

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↗

PubMed2025

Association among blood pressure, antihypertensive drugs, and amyotrophic lateral sclerosis.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) is a fatal and incurable neurodegenerative disease. The impacts of antihypertensive drugs and blood pressure (BP) on ALS are currently debatable. OBJECTIVE: To evaluate the causal relationship involving antihypertensive drugs, BP, and ALS through a Mendelian randomization (MR) analysis. METHODS: The causal relationship between BP and ALS was evaluated by a bidirectional two-sample MR analysis. Then, a sensitivity analysis was performed using a secondary BP genome-wide association study. The drug-target MR was employed to evaluate the impact of antihypertensive drugs on ALS. Furthermore, we used cis-expression quantitative trait loci (cis-eQTLs) data from brain tissue and blood to validate the positive results by a summary-based MR method. RESULTS: We found that an increment in systolic BP (SBP) could elevate the risk of ALS (inverse-variance weighted [IVW] odds ratio [OR] = 1.003; 95% confidence interval [95%CI]: 1.001-1.006; per 10-mmHg increment) and ALS might be protected by angiotensin-converting enzyme inhibitors (ACEIs; OR = 0.970; 95%CI: 0.956-0.984; p = 1.96 × 10-5; per 10-mmHg decrement). A causal relationship was not observed between diastolic BP and other antihypertensive drugs in ALS. CONCLUSION: In the present study, genetic support for elevated SBP serves as a risk factor for ALS. Besides, ACEIs hold promise as a candidate for ALS.

Humans↗

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↗

PubMed2025

EP300-mediated lactylation leads to ulcerative colitis via CD86-positive plasmacytoid dendritic cells: A Mendelian randomization and mediation analysis.

This study explores the potential mechanism between lactylation and ulcerative colitis (UC) using two-sample Mendelian randomization and multi-omics analysis. This study employed expression quantitative trait loci and protein quantitative trait loci as exposures, with UC from the Finnish database as the outcome, to conduct Mendelian randomization analysis on lactylation-related target genes, aiming to investigate the causal relationships between these exposures and the outcome. Sensitivity and pleiotropy tests, combined with colocalization analysis, are performed to identify the best target genes and ensure the robustness of the results. Finally, immune cells are included for mediation analysis between lactylation and UC to explore potential mechanisms of action. Through Mendelian randomization analysis combined with sensitivity and pleiotropy tests, 2 lactylation target genes were found to have a significant causal relationship with UC. Subsequent colocalization analysis confirmed EP300 as a potential gene target. After including immune cells in the mediation analysis, it was discovered that there is a potential mechanism involving EP300, CD86+ plasmacytoid dendritic cells (pDCs), and UC. There is a significant causal relationship between lactylation and UC. Furthermore, the lactylation-modified gene EP300 may lead to UC occurrence by regulating CD86+ pDCs.

Humans↗

PubMed2026

Deciphering miRNA-mediated genetic architecture of immune cell subsets in hypertrophic scars and keloids: A 2-step Mendelian randomization study unveiling causal associations.

This study aimed to investigate the potential causal roles of specific circulating microRNAs (miRNAs) and immune cell subsets in the pathogenesis of hypertrophic scars and keloids using a 2-step Mendelian randomization framework. We employed a 2-sample Mendelian randomization approach to evaluate the causal relationships between miRNAs, immune cell genotypes, and scar phenotypes. The analysis integrated miRNA expression quantitative trait loci, immune cell genome-wide association studies, and scar datasets. A 2-step mediation analysis was conducted to assess the indirect effects of miRNAs on scars through immune cell genotypes, using inverse variance weighted methods and complementary sensitivity analyses to ensure robustness. Our analysis identified significant associations between specific miRNAs and scar phenotypes. Notably, miR-6887-5p exhibited a total effect on keloid formation risk (β = 0.324, 95% confidence interval [CI]: 0.073-0.576) and a direct effect (β = 0.283, 95% CI: 0.027, 0.538), with a marginally significant mediation effect through B-cell activating factor receptor on CD20- CD38- B cells (β = 0.042, 95% CI: -0.001, 0.084, P = .047). For hypertrophic scars, miR-345-5p demonstrated a significant total effect (β = -0.501, 95% CI: -0.903, -0.099) and direct effect (β = -0.469, 95% CI: -0.872, -0.066), with a significant mediation effect through CD28+ CD45RA- CD8dim T cell percentage (β = -0.032, 95% CI: -0.062, -0.002, P = .034). miR-4801 showed a significant total effect (β = -0.246, 95% CI: -0.429, -0.064) and direct effect (β = -0.218, 95% CI: -0.402, -0.033), with a marginally significant mediation effect through T cell absolute count (β = -0.028, 95% CI: -0.057, -0.000, P = .043). These findings highlight the interplay between miRNAs and immune cell subsets in scar pathogenesis. This study provides preliminary evidence for the causal roles of specific miRNAs and immune cell subsets in scar formation, emphasizing the potential of miRNA-immune cell axes as therapeutic targets. While the identified associations offer important insights into the molecular mechanisms of scar heterogeneity, further validation through mechanistic studies and clinical trials is necessary to translate these genetic insights into clinical interventions.

Humans↗