PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “expression quantitative trait loci”

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 37 records · Page 2Linked to original sources

Integrating genotypic and expression data in a segregating mouse population to identify 5-lipoxygenase as a susceptibility gene for obesity and bone traits.

Forward genetic approaches to identify genes involved in complex traits such as common human diseases have met with limited success. Fine mapping of linkage regions and validation of positional candidates are time-consuming and not always successful. Here we detail a hybrid procedure to map loci involved in complex traits that leverages the strengths of forward and reverse genetic approaches. By integrating genotypic and expression data in a segregating mouse population, we show how clusters of expression quantitative trait loci linking to regions of the genome accurately reflect the underlying perturbation to the transcriptional network induced by DNA variations in genes that control the complex traits. By matching patterns of gene expression in a segregating population with expression responses induced by single-gene perturbation experiments, we show how genes controlling clusters of expression and clinical quantitative trait loci can be mapped directly. We demonstrate the utility of this approach by identifying 5-lipoxygenase as underlying previously identified quantitative trait loci in an F(2) cross between strains C57BL/6J and DBA/2J and showing that it has pleiotropic effects on body fat, lipid levels and bone density.

Animals↗

Global eQTL mapping reveals the complex genetic architecture of transcript-level variation in Arabidopsis.

The genetic architecture of transcript-level variation is largely unknown. The genetic determinants of transcript-level variation were characterized in a recombinant inbred line (RIL) population (n = 211) of Arabidopsis thaliana using whole-genome microarray analysis and expression quantitative trait loci (eQTL) mapping of transcript levels as expression traits (e-traits). Genetic control of transcription was highly complex: one-third of the quantitatively controlled transcripts/e-traits were regulated by cis-eQTL, and many trans-eQTL mapped to hotspots that regulated hundreds to thousands of e-traits. Several thousand eQTL of large phenotypic effect were detected, but almost all (93%) of the 36,871 eQTL were associated with small phenotypic effects (R(2) < 0.3). Many transcripts/e-traits were controlled by multiple eQTL with opposite allelic effects and exhibited higher heritability in the RILs than their parents, suggesting nonadditive genetic variation. To our knowledge, this is the first large-scale global eQTL study in a relatively large plant mapping population. It reveals that the genetic control of transcript level is highly variable and multifaceted and that this complexity may be a general characteristic of eukaryotes.

Arabidopsis↗

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↗

Expressionview: visualization of quantitative trait loci and gene-expression data in Ensembl.

We present here a software tool for combined visualization of gene-expression data and quantitative trait loci (QTL). The application is implemented as an extension to the Ensembl project and caters for a direct transition from microarray experiments of gene or protein expression levels to the genomic context of individual genes and QTL. It supports the visualization of gene clusters and the selection of functional candidate genes in the context of research on complex traits.

Animals↗

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

Genetics of gene expression surveyed in maize, mouse and man.

Treating messenger RNA transcript abundances as quantitative traits and mapping gene expression quantitative trait loci for these traits has been pursued in gene-specific ways. Transcript abundances often serve as a surrogate for classical quantitative traits in that the levels of expression are significantly correlated with the classical traits across members of a segregating population. The correlation structure between transcript abundances and classical traits has been used to identify susceptibility loci for complex diseases such as diabetes and allergic asthma. One study recently completed the first comprehensive dissection of transcriptional regulation in budding yeast, giving a detailed glimpse of a genome-wide survey of the genetics of gene expression. Unlike classical quantitative traits, which often represent gross clinical measurements that may be far removed from the biological processes giving rise to them, the genetic linkages associated with transcript abundance affords a closer look at cellular biochemical processes. Here we describe comprehensive genetic screens of mouse, plant and human transcriptomes by considering gene expression values as quantitative traits. We identify a gene expression pattern strongly associated with obesity in a murine cross, and observe two distinct obesity subtypes. Furthermore, we find that these obesity subtypes are under the control of different loci.

Animals↗

Mapping determinants of gene expression plasticity by genetical genomics in C. elegans.

Recent genetical genomics studies have provided intimate views on gene regulatory networks. Gene expression variations between genetically different individuals have been mapped to the causal regulatory regions, termed expression quantitative trait loci. Whether the environment-induced plastic response of gene expression also shows heritable difference has not yet been studied. Here we show that differential expression induced by temperatures of 16 degrees C and 24 degrees C has a strong genetic component in Caenorhabditis elegans recombinant inbred strains derived from a cross between strains CB4856 (Hawaii) and N2 (Bristol). No less than 59% of 308 trans-acting genes showed a significant eQTL-by-environment interaction, here termed plasticity quantitative trait loci. In contrast, only 8% of an estimated 188 cis-acting genes showed such interaction. This indicates that heritable differences in plastic responses of gene expression are largely regulated in trans. This regulation is spread over many different regulators. However, for one group of trans-genes we found prominent evidence for a common master regulator: a transband of 66 coregulated genes appeared at 24 degrees C. Our results suggest widespread genetic variation of differential expression responses to environmental impacts and demonstrate the potential of genetical genomics for mapping the molecular determinants of phenotypic plasticity.

Animals↗

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↗

Design of microarray experiments for genetical genomics studies.

Microarray experiments have been used recently in genetical genomics studies, as an additional tool to understand the genetic mechanisms governing variation in complex traits, such as for estimating heritabilities of mRNA transcript abundances, for mapping expression quantitative trait loci, and for inferring regulatory networks controlling gene expression. Several articles on the design of microarray experiments discuss situations in which treatment effects are assumed fixed and without any structure. In the case of two-color microarray platforms, several authors have studied reference and circular designs. Here, we discuss the optimal design of microarray experiments whose goals refer to specific genetic questions. Some examples are used to illustrate the choice of a design for comparing fixed, structured treatments, such as genotypic groups. Experiments targeting single genes or chromosomic regions (such as with transgene research) or multiple epistatic loci (such as within a selective phenotyping context) are discussed. In addition, microarray experiments in which treatments refer to families or to subjects (within family structures or complex pedigrees) are presented. In these cases treatments are more appropriately considered to be random effects, with specific covariance structures, in which the genetic goals relate to the estimation of genetic variances and the heritability of transcriptional abundances.

Animals↗

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↗

Integrated gene expression profiling and linkage analysis in the rat.

The combined application of genome-wide expression profiling from microarray experiments with genetic linkage analysis enables the mapping of expression quantitative trait loci (eQTLs) which are primary control points for gene expression across the genome. This approach allows for the dissection of primary and secondary genetic determinants of gene expression. The cis-acting eQTLs in practice are easier to investigate than the trans-regulated eQTLs because they are under simpler genetic control and are likely to be due to sequence variants within the gene itself or its neighboring regulatory elements. These genes are therefore candidates both for variation in gene expression and for contributions to whole-body phenotypes, particularly when these are located within known and relevant physiologic QTLs. Multiple trans-acting eQTLs tend to cluster to the same genetic location, implying shared regulatory control mechanisms that may be amenable to network analysis to identify gene clusters within the same metabolic pathway. Such clusters may ultimately underlie development of individual complex, whole-body phenotypes. The combined expression and linkage approach has been applied successfully in several mammalian species, including the rat which has specific features that demonstrate its value as a model for studying complex traits.

Animals↗

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↗

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↗

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↗

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↗

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&#x2009;=&#x2009;1.02, 95%CI&#x2009;=&#x2009;1.01-1.04, FDR adjusted P&#x2009;=&#x2009;1.69&#x2009;&#xd7;&#x2009;10-2; OR in putamen&#x2009;=&#x2009;1.02, 95% CI&#x2009;=&#x2009;1.01-1.03, PFDR&#x2009;=&#x2009;3.37&#x2009;&#xd7;&#x2009;10-2, OR in nucleus accumbens&#x2009;=&#x2009;1.02, 95% CI&#x2009;=&#x2009;1.01-1.04, PFDR&#x2009;=&#x2009;3.37&#x2009;&#xd7;&#x2009;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↗

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↗

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&#x2009;<&#x2009;0.001). The genotype of rs1738023 was significantly associated with AKR7A3 expression (P&#x2009;=&#x2009;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↗