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 19 recordsLinked to original sources

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

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

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

Expanded Chromatin Accessibility Mapping Explains Genetic Variation Associated with Complex Traits in Liver.

Genome-wide association studies (GWAS) have identified thousands of loci associated with a variety of common, complex human traits. Recent efforts have focused on characterizing chromatin accessibility to discover regulatory elements that modify the expression of nearby genes, suggesting that trait associations are mediated through changes in gene regulation. Genetic variants associated with differences in chromatin accessibility, known as chromatin accessibility quantitative trait loci (caQTLs), are established contributors to gene expression differences, providing mechanistic hypotheses for signals identified by GWAS. Using the assay for transposase-accessible chromatin with sequencing (ATAC-seq), we assessed chromatin accessibility in 189 diverse human liver samples, identifying over two million accessible chromatin regions enriched for gene regulatory features and, in 175 of these samples, over 14,000 caQTLs. Focusing subsequently on liver-relevant complex traits, we obtained publicly available blood lipids GWAS data and identified 157 loci where caQTLs, expression quantitative trait loci (eQTLs), and GWAS signals colocalized. This generated specific molecular hypotheses about regulatory elements, affected genes, and, in some cases, implicated transcription factors. Finally, we enumerated the set of blood lipid trait signals that lack an obvious proposed mechanism beyond catalogs of liver caQTLs and eQTLs. After integrating 10 multi-omic QTL regulatory mechanism datasets whilst considering limitations in statistical power, we found that approximately 20% of blood lipid GWAS signals lacked a statistical link to a proposed mechanism. Our results demonstrate the value of integrating multiple genomic datasets to improve understanding of GWAS signals, while emphasizing the need for additional experimental approaches to fully characterize complex trait associations.

Journal Article↗

PubMed2026

Plasma proteins are integral to cross-tissue gene regulatory networks implicated in cardiometabolic disorders and coronary artery disease.

The plasma proteome has demonstrated promise for identifying diagnostic markers for cardiometabolic disorders (CMDs) and coronary artery disease (CAD). However, identifying the organ of origin for these biomarkers is critical for establishing biological relevance. We performed a multi-omic integrative analysis across multiple tissues from the STARNET study by profiling 974 plasma proteins in 532 CAD patients, integrating RNA sequencing (RNA-seq) data from the arterial wall, major metabolic organs, and blood. We identified 144 cis-protein quantitative trait loci in plasma, colocalizing with tissue cis-expression quantitative trait loci. Additionally, by mapping tissue mRNA "seed genes," we traced 262 plasma proteins to their source organs, primarily the liver. Crucially, we found that 851 plasma proteins are associated with the activity of cross-tissue gene regulatory networks (GRNs), including GRNs implicated in CMD and CAD development. Our findings demonstrate that plasma proteins are integral components of GRNs, with potential for developing reliable diagnostics and precise therapeutic targets. A record of this paper's transparent peer review process is included in the supplemental information.

cardiometabolic disorders↗

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

A spectral framework to map QTLs affecting joint differential networks of gene co-expression.

Studying the mechanisms underlying the genotype-phenotype association is crucial in genetics. Gene expression studies have deepened our understanding of the genotype  →  expression  →  phenotype mechanisms. However, traditional expression quantitative trait loci (eQTL) methods often overlook the critical role of gene co-expression networks in translating genotype into phenotype. This gap highlights the need for more powerful statistical methods to analyze genotype  →  network  →  phenotype mechanism. Here, we develop a network-based method, called spectral network quantitative trait loci analysis (snQTL), to map quantitative trait loci affecting gene co-expression networks. Our approach tests the association between genotypes and joint differential networks of gene co-expression via a tensor-based spectral statistics, thereby overcoming the ubiquitous multiple testing challenges in existing methods. We demonstrate the effectiveness of snQTL in the analysis of three-spined stickleback (Gasterosteus aculeatus) data. Compared to conventional methods, our method snQTL uncovers chromosomal regions affecting gene co-expression networks, including one strong candidate gene that would have been missed by traditional eQTL analyses. Our framework suggests the limitation of current approaches and offers a powerful network-based tool for functional loci discoveries.

Quantitative Trait Loci↗

PubMed2026

Integrative cross-tissue transcriptome-wide association and metabolomic analysis reveals novel genetic risk loci for aortic aneurysm.

BACKGROUND: Aortic aneurysm (AA) is a life-threatening cardiovascular condition with a strong genetic component, however, its molecular mechanisms remain poorly understood. Although genome-wide association studies (GWAS) have identified numerous risk loci, most prior studies have investigated genetic and metabolic factors separately, leaving the causal pathways from genetic variants to disease largely unexplored. METHODS: We established an integrative framework combining cross-tissue transcriptome-wide association studies (TWAS) with metabolomic mediation analysis. First, we integrated GWAS data from FinnGen R12 with multi-tissue expression quantitative trait loci (eQTL) data from Genotype-Tissue Expression Project (GTEx) V8, then performed cross-tissue TWAS using the Unified Test for MOlecular SignaTures (UTMOST) and single-tissue validation with the Functional Summary-based Imputation (FUSION) to prioritize susceptibility genes. Second, we applied Mendelian randomization (MR), colocalization, and Fine-mapping Of CaUsal gene Sets (FOCUS) to assess causality and identify high-confidence genes. Third, we performed metabolite mediation analysis to uncover metabolic pathways linking genetic variants to disease risk. Finally, we validated key findings in mouse models of thoracic aortic aneurysm (TAA) and abdominal aortic aneurysm (AAA) using Quantitative Real-Time Reverse Transcription Polymerase Chain Reaction (RT-qPCR) and Western blotting. RESULTS: We identified multiple novel susceptibility genes for AA and its subtypes. Key genes included ADH family members (ADH1A, ADH1B, ADH4, ADH6) and ZNF827, which showed cross-subtype associations with strong colocalization evidence in vascular tissues. Metabolite mediation analysis revealed significant pathways involving N-acetylphenylalanine and methionine sulfoxide. Functional enrichment revealed distinct biological mechanisms: AA and AAA were primarily associated with metabolic pathways, whereas TAA-related genes were enriched in developmental and contractile processes. PheWAS indicated no significant off-target associations. Critically, experimental validation in mouse models confirmed significant upregulation of ZNF827 in TAA and ADH6 in AAA at both mRNA and protein levels, corroborating the genetic predictions. CONCLUSION: This integrated cross-omics analysis identifies novel genetic loci and, crucially, uncovers specific nutrient-related metabolic pathways that mediate genetic risk. These findings provide a mechanistic basis for future nutritional and metabolic intervention studies in AA and its subtypes.

MAGMA↗

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)↗

PubMed2025

Prioritization of causal genes from genome-wide association studies by Bayesian data integration across loci.

MOTIVATION: Genome-wide association studies (GWAS) have identified genetic variants, usually single-nucleotide polymorphisms (SNPs), associated with human traits, including disease and disease risk. These variants (or causal variants in linkage disequilibrium with them) usually affect the regulation or function of a nearby gene. A GWAS locus can span many genes, however, and prioritizing which gene or genes in a locus are most likely to be causal remains a challenge. Better prioritization and prediction of causal genes could reveal disease mechanisms and suggest interventions. RESULTS: We describe a new Bayesian method, termed SigNet for significance networks, that combines information both within and across loci to identify the most likely causal gene at each locus. The SigNet method builds on existing methods that focus on individual loci with evidence from gene distance and expression quantitative trait loci (eQTL) by sharing information across loci using protein-protein and gene regulatory interaction network data. In an application to cardiac electrophysiology with 226 GWAS loci, only 46 (20%) have within-locus evidence from Mendelian genes, protein-coding changes, or colocalization with eQTL signals. At the remaining 180 loci lacking functional information, SigNet selects 56 genes other than the minimum distance gene, equal to 31% of the information-poor loci and 25% of the GWAS loci overall. Assessment by pathway enrichment demonstrates improved performance by SigNet. Review of individual loci shows literature evidence for genes selected by SigNet, including PMP22 as a novel causal gene candidate.

Genome-Wide Association Study↗

PubMed2025

Endogenous fine-mapping and prioritization of functional regulatory elements in complex genetic loci.

Most genetic loci linked to polygenic traits are in non-coding regions, with complex regulation and linkage disequilibrium (LD), complicating causal variant and gene prioritization. We used multiplexed single-cell CRISPR interference and activation perturbations to investigate cis-regulatory element (CRE) and gene expression relationships within tight LD in the endogenous chromatin context. We demonstrated the prevalence of multiple causality in perfect LD (pLD) for independent expression quantitative trait loci (eQTLs) and uncovered fine-grained genetic effects on gene expression within pLD, which are difficult to decipher using traditional eQTL fine-mapping or existing computational methods. We found that over one-third of the causal CREs lack classical epigenetic markers prior to perturbation, and we functionally validated one of these hidden regulatory mechanisms. Leveraging Multiome single-cell epigenetic and sequence perturbations, we highlighted the regulatory plasticity of the human genome. Our study will guide the exploration of missing causal mechanisms underlying molecular trait regulation and disease development.

Humans↗

PubMed2026

DNA Methylation-Mediated Regulation of TAGLN2 Expression Promotes Pulmonary Arterial Hypertension.

BACKGROUND: Succinylation, a key post-translational modification, is implicated in the metabolic reprogramming and vascular remodeling of pulmonary arterial hypertension (PAH). While epigenetic regulation, particularly DNA methylation, potentially governs succinylation-related gene expression, its causal links to PAH remain unclear. METHODS: We performed an integrative causal analysis using two-sample Mendelian randomization (MR) and summary-data-based MR (SMR) to identify succinylation-related genes that influence PAH risk. We leveraged PAH GWAS data (FinnGen) and gene expression quantitative trait loci (eQTLGen). Subsequently, methylation-mediated effect decomposition was applied using DNA methylation data (GoDMC) to explore epigenetic regulation. Experimental validation was conducted in lung tissues from a monocrotaline (MCT)-induced PAH rat model via quantitative reverse transcription polymerase chain reaction (qRT-PCR). RESULTS: Genetic analyses identified a significant causal effect of elevated Transgelin 2 (TAGLN2) expression on increased PAH risk. This effect was mediated by two specific DNA methylation sites, cg13892570 and cg16107628, which influenced PAH pathogenesis by regulating TAGLN2 transcription, with mediation proportions of 86.46 and 97.65%, respectively. Sensitivity analyses supported the robustness of these findings. Consistent with the genetic evidence, TAGLN2 mRNA was significantly upregulated in the lungs of MCT-induced PAH rats. CONCLUSIONS: This study establishes a clear epigenetic causal pathway in which DNA methylation regulates TAGLN2 expression to promote PAH. TAGLN2 is validated as a key disease driver and presents a promising target for diagnostic and therapeutic strategies in PAH.

Animals↗

PubMed2024

Genetic and functional analysis of Raynaud's syndrome implicates loci in vasculature and immunity.

Raynaud's syndrome is a dysautonomia where exposure to cold causes vasoconstriction and hypoxia, particularly in the extremities. We performed meta-analysis in four cohorts and discovered eight loci (ADRA2A, IRX1, NOS3, ACVR2A, TMEM51, PCDH10-DT, HLA, and RAB6C) where ADRA2A, ACVR2A, NOS3, TMEM51, and IRX1 co-localized with expression quantitative trait loci (eQTLs), particularly in distal arteries. CRISPR gene editing further showed that ADRA2A and NOS3 loci modified gene expression and in situ RNAscope clarified the specificity of ADRA2A in small vessels and IRX1 around small capillaries in the skin. A functional contraction assay in the cold showed lower contraction in ADRA2A-deficient and higher contraction in ADRA2A-overexpressing smooth muscle cells. Overall, our study highlights the power of genome-wide association testing with functional follow-up as a method to understand complex diseases. The results indicate temperature-dependent adrenergic signaling through ADRA2A, effects at the microvasculature by IRX1, endothelial signaling by NOS3, and immune mechanisms by the HLA locus in Raynaud's syndrome.

Raynaud Disease↗

PubMed2025

Genomic loci and molecular genetic mechanisms for hidradenitis suppurativa.

BACKGROUND: Hidradenitis suppurativa (HS) is a common, chronic and debilitating inflammatory disease that most commonly affects intertriginous skin. Despite its high heritability, the genetic underpinnings of HS remain poorly understood. OBJECTIVES: To identify genetic signals associated with HS, determine genetic relationships with other diseases and investigate potential molecular genetic mechanisms. METHODS: We performed a genome-wide association meta-analysis of six studies, totalling 4540 patients with HS and > 1 million control participants, and identified genetic correlations with other common diseases. We integrated the HS data with expression quantitative trait loci from 10 trait-relevant tissues, epigenomic and transcriptomic data from human scalp, differential expression data from HS lesions vs. adjacent skin and mesenchymal Hi-C chromatin looping data. To identify functional noncoding variants, we performed transcriptional reporter assays for signals near KLF5 and SOX9. RESULTS: We identified 11 significant HS signals across 7 loci: 4 corresponded to previously reported associations, 4 represented novel signals within known loci and 3 were signals in newly implicated loci. We identified significant genetic correlations between HS and other inflammatory conditions, particularly inflammatory bowel disease, rheumatoid arthritis, type 2 diabetes mellitus and asthma. We prioritized candidate genes for the 11 signals. The risk allele at KLF5 exhibited 10-fold greater transcriptional activity than the nonrisk allele, while risk alleles at SOX9 showed significantly reduced transcriptional activity. CONCLUSIONS: Our results provide insights into potential genetic mechanisms underlying HS and suggest potential therapeutic targets for this challenging condition.

Humans↗

PubMed2026

Identification of Critical Genes for Recurrent Aphthous Ulcer by Transcriptome Data Analysis and Mendelian Randomization.

PURPOSE: Recurrent aphthous ulcer (RAU) is a common oral mucosal disorder with a poorly understood etiology, significantly affecting patients' quality of life. This study aims to investigate critical genes linked to RAU and explore their biological mechanisms using transcriptomic data and Mendelian randomization (MR) analysis. MATERIALS AND METHODS: RAU-related gene expression data from the GEO database (GSE37265) were analyzed to identify differentially expressed genes (DEGs). A two-sample MR approach was used to assess the causal impact of expression quantitative trait loci (eQTL) on RAU. Critical genes were identified by intersecting DEGs with significant MR findings. GO and KEGG pathway enrichment analyses were performed, along with GSEA and immune cell infiltration analysis, to investigate the functions and mechanisms of these genes in RAU. RESULTS: A total of 184 differentially expressed genes (DEGs) were identified, while 339 RAU-associated genes were screened through MR analysis. Cross-validation further identified 7 critical genes. Among these, CCR1, ERP27, HCK, MICB, and SLC2A3 showed protective associations with RAU risk, whereas CD177 and IFITM1 were positively associated with increased risk. Enrichment analysis revealed that these genes are involved in specific biological processes, including cell migration, immune response, and metabolic regulation, which are closely linked to RAU pathogenesis. CONCLUSION: This systematic study comprehensively investigates the critical causative genes underlying RAU, emphasizing the intricate relationships between immune regulation and metabolic disturbances in its pathology. These findings lay a solid foundation for the development of novel biomarkers and may inform future research on targeted therapeutic strategies for RAU.

Stomatitis, Aphthous↗

PubMed2026

The tissue-specific effects of glucose-lowering drug targets on aging mediated through DNA methylation: a multi-omics genetic study.

BACKGROUND: DNA methylation plays a key role in mediating the anti-aging effects of glucose-lowering drugs. This study aims to systematically explore the potential anti-aging effects of target genes of FDA-approved glucose-lowering drugs and the underlying epigenetic mediators. METHODS: We conducted a two-sample Mendelian randomization (MR) study to investigate the putative causal relationships between the gene expression levels of glucose-lowering drug targets and 10 aging-related phenotypes, followed by a two-step MR to estimate the mediation effect of DNA methylation. Drug candidates were selected according to the latest review of clinical drug use for type 2 diabetes, and their target genes were obtained from the DGIdb. Tissue-specific cis-expression quantitative trait loci (eQTLs) from GTEx Consortium were selected as genetic instruments to proxy the expression level of drug-target genes. Glycemic phenotypes were used as positive controls to validate the instruments. The cis- and trans-methylation QTLs of Cytosine-phosphate-Guanine sites near the drug target genes were obtained from GoDMC Consortium. Additionally, we performed enrichment analyses focused on tissue specificity and aging pathways to further corroborate our findings. RESULTS: We obtained 194 target genes interacting with 36 FDA-approved anti-diabetic drugs, of which the tissue-specific eQTLs were used to proxy the drug target effects. MR showed strong evidence that nine interacting genes of six glucose-lowering drugs showed anti-aging potential on one or more aging-related phenotypes mediated by DNA methylation: EHMT2, HSPA4, IGF2BP2, IRS1, LPL, NDUFAF1, NDUFS3, SLC22A3, and TCF7L2. These genes were distributed in 17 tissues, especially in the central nervous system, suggesting a potential neural component in their anti-aging effects. For instance, expression of EHMT2 in several brain basal ganglia regions, where the gene interacted with Tolazamide, showed a protective effect on frailty (odds ratio (OR) in caudate = 1.02, 95%CI = 1.01-1.04, FDR adjusted P = 1.69 × 10-2; OR in putamen = 1.02, 95% CI = 1.01-1.03, PFDR = 3.37 × 10-2, OR in nucleus accumbens = 1.02, 95% CI = 1.01-1.04, PFDR = 3.37 × 10-2). These associations were externally validated by searching literature evidence in existing EWAS and TWAS studies, as well as evidence from enrichment analyses. CONCLUSIONS: This study prioritizes nine glucose-lowering genes as anti-aging drug targets in specific tissues and prioritizes their epigenetic regulation through DNA methylation for future drug development.

DNA Methylation↗

PubMed2026

Single-Cell Transcriptome-Wide Mendelian Randomization and Colocalization Uncover Potential Immunocytes-Related Therapeutic Targets for Obesity.

Weight-loss treatment is crucial for individuals with obesity to prevent various complications. The role of Immune cells in obesity has been recently recognized, whereas its translation into therapy requires identifying key target genes. We performed Mendelian randomization (MR) analysis to assess causal relationships between expression quantitative trait loci (eQTL) of 14 immune cells and obesity-related traits (obesity, body mass index and body fat percentage), and validated the results in colocalization analysis. For the putative causal genes identified by the MR and colocalization analyses, we conducted pathway enrichment, differential expressed gene (DEG) analysis and search of druggable evidence, and utilized a Tier system to prioritize drug targets for obesity. MR and colocalization evidence was observed for 1630 genes associated with one or more obesity-related traits, mainly expressed in CD4+ naive/central memory T cells and enriched in antigen processing and presentation pathways. Forty-one genes showed causal relationship with all three outcomes, among which 19 genes have not been reported for obesity previously. DEG analysis using single-cell RNA sequencing data of blood or adipose tissue indicated that the differential expression of UBE2Z in monocytes, ZCCHC7 in T cells, and FNBP4 in B cells between lean and obese individuals were consistent with the MR results. By searching drug-gene interaction databases, we found targeted drugs for PYGB and PRUNE1, and PYGB was the top gene ranked in the Tier system. This study provides evidence for the involvement of immune cells in obesity, and the potential cell-specific, immune-related targets for obesity treatment.

Obesity↗

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↗