SEARCH · PubMed Health

Results for “eQTL”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

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

Saturating the eQTL map in Drosophila: Genome-wide patterns of cis and trans regulation of transcriptional variation in outbred populations.

Most genetic polymorphisms associated with complex traits are found in non-coding regions of the genome. Characterizing their effect presents a formidable challenge, and expression quantitative trait locus (eQTLs) mapping has been a key approach to do so. As comprehensive eQTL maps are available only for a few species, here we developed the Drosophila outbred synthetic population (Dros-OSP) and used it to characterize the landscape of transcriptional regulation in Drosophila melanogaster. We collected head and body transcriptomes and genomes from 1,286 outbred flies and mapped local and distant eQTLs for 98% of the genes. We characterized the network organization of the transcriptome across tissues and described the properties of local and distal eQTLs in terms of genetic diversity, heritability, connectivity, and pleiotropy. These results provide new insights into the genetic basis of transcriptional regulation in the fruit fly and offer a new mapping resource that will expand the possibilities currently available for the Drosophila community.

Animals↗

PubMed2026

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

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

Quantitative Trait Loci↗

PubMed2026

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

Identification of novel type 1 and type 2 diabetes genes by co-localization of human islet eQTL and GWAS variants with colocRedRibbon.

Over 1,000 genetic variants have been associated with diabetes by genome-wide association studies (GWASs), but for most, their functional impact is unknown; only 7% alter gene expression in pancreatic islets in expression quantitative trait locus (eQTL) studies. To fill this gap, we developed a co-localization pipeline, colocRedRibbon, that prefilters eQTLs by the direction of effect on gene expression and shortlists overlapping eQTL and GWAS variants prior to co-localization. Applying colocRedRibbon to recent diabetes and glycemic trait GWASs, we identified 292 co-localizing gene regions, including 24 co-localizations for type 1 diabetes and 268 for type 2 diabetes and glycemic traits, representing a 4-fold increase. A low-frequency type 2 diabetes protective variant increases islet MYO5C expression, and a type 1 diabetes protective variant increases FUT2 expression. These novel co-localizations advance the understanding of diabetes genetics and its impact on human islet biology. colocRedRibbon has broad applicability to co-localize GWASs and various QTLs.

Humans↗

PubMed2026

Immune Cell-Stratified Regulatory Contexts Associated With BMI-Related Multi-System Disease Risk: A Cell-Stratified Mendelian Randomization Study Using Single-Cell eQTL Data.

AIMS: Body mass index (BMI) is associated with multisystem disease risk, but the immune cell-specific regulatory contexts underlying BMI-related genetic associations with disease outcomes remain unclear. METHODS: We applied a cell-stratified Mendelian randomization framework integrating European-ancestry BMI GWAS data, GWAS datasets for 33 disease outcomes across five disease systems, single-cell cis-eQTL data from 28 peripheral blood immune cell types, and dynamic CD4+ T cell eQTL data. SuSiE-based colocalization was used to identify BMI-associated loci sharing causal variants with immune-cell gene expression. These variants were used as cell-stratified instruments for Mendelian randomization. RESULTS: Across 28 immune cell types, 1326 colocalized variants regulating 1426 genes were identified. In primary MR analyses, genetically predicted BMI showed Bonferroni-significant associations with 26 disease outcomes. Cell-stratified analyses identified 87 Bonferroni-significant associations across 17 disease outcomes. Cardiovascular diseases showed the broadest cell-stratified associations, followed by respiratory and metabolic diseases. CD4+ T cell regulatory contexts contributed one of the largest shares of prioritized associations, and BMI-related effects varied across CD4+ T cell activation states. Cross-disease prioritization highlighted recurrent immune feature genes, including TRAF3 and FGFR1. CONCLUSION: These findings prioritize CD4+ T cell regulatory contexts as potential immunogenetic links between BMI and multi-system disease risk, while requiring further validation in diverse populations and mechanistic models.

Humans↗

PubMed2026

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

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

Journal Article↗

PubMed2026

Cell-type-specific genetic associations in Lewy body dementia identified using single-cell eQTL-based Mendelian randomization.

BACKGROUND: Lewy body dementia (LBD) is a complex neurodegenerative disorder marked by α-synuclein aggregation and dual impairment of cognitive and motor function.While genome-wide association studies have identified risk loci, the cellular mechanisms linking genetic variation to disease susceptibility remain largely unexplored. METHODS: We performed single-cell transcriptome-wide Mendelian randomization using brain cell-type-specific eQTLs across eight major cell types. Genetic associations were evaluated using inverse-variance weighted models, followed by Bayesian colocalization analysis. Replication was performed in independent stratified LBD cohorts based on APOE ε4 carrier status. Phenome-wide association analysis was included as a supplementary, descriptive assessment of cross-trait associations. RESULTS: Expression of ANKRD65 in excitatory neurons was significantly associated with reduced LBD risk (odds ratio = 0.65, 95 % CI: 0.52-0.81, p = 0.00013). This association passed a false discovery rate of 0.1 and showed strong evidence of colocalization (posterior probability = 0.93). Effect direction was consistent across APOE ε4+ and ε4- LBD subgroups in independent cohorts. No genome-wide significant associations were observed with non-neurological traits in the phenome-wide analysis. CONCLUSIONS: Our findings identify a genetically supported, cell-type-resolved association between ANKRD65 expression in excitatory neurons and LBD risk. This study demonstrates the value of integrating cell-resolved transcriptomic regulation with genetic inference to pinpoint functionally relevant targets in neurodegenerative diseases.

Humans↗

PubMed2024

eQTLs identify regulatory networks and drivers of variation in the individual response to sepsis.

Sepsis is a clinical syndrome of life-threatening organ dysfunction caused by a dysregulated response to infection, for which disease heterogeneity is a major obstacle to developing targeted treatments. We have previously identified gene-expression-based patient subgroups (sepsis response signatures [SRS]) informative for outcome and underlying pathophysiology. Here, we aimed to investigate the role of genetic variation in determining the host transcriptomic response and to delineate regulatory networks underlying SRS. Using genotyping and RNA-sequencing data on 638 adult sepsis patients, we report 16,049 independent expression (eQTLs) and 32 co-expression module (modQTLs) quantitative trait loci in this disease context. We identified significant interactions between SRS and genotype for 1,578 SNP-gene pairs and combined transcription factor (TF) binding site information (SNP2TFBS) and predicted regulon activity (DoRothEA) to identify candidate upstream regulators. Overall, these approaches identified putative mechanistic links between host genetic variation, cell subtypes, and the individual transcriptomic response to infection.

Humans↗

PubMed2026

Cell type resolved MR based on brain single cell eQTLs corroborated by single cell RNA sequencing uncovers neuroimmune and vascular programs in intracerebral hemorrhage.

BACKGROUND: Intracerebral hemorrhage (ICH) lacks effective neuroprotective therapies. We integrated cell type–resolved genetic inference with single-cell profiling to map putative causal programs and multicellular circuitry relevant to ICH. METHODS: Cis-eQTLs from eight human brain cell types were used as instruments for two-sample Mendelian randomization (MR), with an ICH meta-analysis from large biobanks and a stroke consortium as the outcome. Instruments were LD-pruned and restricted to strong variants (F > 10). Inverse-variance weighting (IVW) was the primary estimator, supported by robustness methods, heterogeneity/pleiotropy diagnostics, and false discovery rate control. Experimental validation used mouse collagenase ICH single-cell RNA-seq at 24 h (n = 3 sham; n = 3 ICH) with Seurat integration, composition testing, Slingshot pseudotime, and CellChat. An independent mouse cohort underwent qRT–PCR for selected genes. RESULTS: The ICH meta-analysis showed acceptable genomic control, supporting downstream MR. We identified 524 nominal gene–cell type associations, with a glia-weighted signal landscape. Enrichment implicated autophagy/mitophagy, antigen processing, cytoskeletal and vesicular trafficking, endothelial matrix–adhesion programs, ferroptosis, and myelin stress pathways. In mouse scRNA-seq, disease-associated microglia expanded with reciprocal loss of homeostatic microglia and increased neutrophils and T cells. Prioritized genes showed directional concordance; qRT–PCR confirmed ARPC3 and EIF2AK2 upregulation and TBCK and SPECC1 downregulation in ICH versus sham. Pseudotime supported a shift toward disease-associated microglial states, and CellChat indicated increased network interaction strength with microglia and endothelium as hubs. CONCLUSIONS: Cell type–specific MR combined with single-cell validation highlights neuroimmune and neurovascular programs in ICH and links genetic signals to state transitions and inferred intercellular communication.

Animals↗

PubMed2024

Relationship between inflammation/immunity and epilepsy: A multi-omics mendelian randomization study integrating GWAS, eQTL, and mQTL data.

OBJECTIVES: Increasing evidence suggests that activated innate/adaptive immunity induces an inflammatory response, thereby participating in epileptogenesis. However, the biological explanation of inflammation/immunity as a potential cause for epilepsy remains largely unknown. This research aimed to determine the causal effects of inflammation/immune-related genes in epilepsy based on multi-omics mendelian randomization (MR). METHODS: We employed summary-data-based MR (SMR) approach to combine GWAS for epilepsy (12,891 cases and 312,803 control) with gene expression quantitative trait loci (cis-eQTL, 31,684 participants) and DNA methylation QTL (cis-mQTL, 1,980 participants) data. Five additional MR methods were then used for sensitivity analyses to confirm the reliability of causal associations. In addition, enrichment analysis of key genes was conducted to provide insight into the biological functions of epilepsy risk variants. RESULTS: A total of 386 inflammation/immune-related genes were selected for further analyses. Primary SMR analysis indicated that 37 DNA methylation sites and six genes regulated by them had potential causal relationship with epilepsy. MR analysis further refined the results, identifying three genes that had a causal effect on epilepsy. Notably, VEGFA (OR: 0.925; 95 % CI: 0.862-0.994) expression was negatively correlated with epilepsy risk, whereas IL16 (OR: 1.076; 95 % CI: 1.028-1.126) and HLA-DPA1 (OR: 1.041; 95 % CI: 1.009-1.074) expressions were positively associated with epilepsy risk. Functional enrichment analysis revealed that the identified genes were involved in GO-BP terms related to VEGF activation signaling and chemotaxis regulation. CONCLUSION: This analysis confirms the causal role of inflammation/immunity in epilepsy, and the identified candidate genes provide clues for drug development in clinical practice.

Humans↗

PubMed2025

Multi-omics Mendelian randomization integrating RNA-seq, eQTL and pQTL data revealed CPXM1 as a potential drug target for osteoporosis.

Osteoporosis, a prevalent skeletal disorder characterized by decreased bone mineral density and increased fracture risk, continues to be a major global health concern. Traditional treatments for osteoporosis have limited efficacy and safety profiles, highlighting the need for novel therapeutic targets. This study integrates multi-omics data, including RNA-seq, expression quantitative trait loci (eQTL), and protein quantitative trait loci (pQTL) data, through Mendelian randomization (MR) to identify potential drug targets for osteoporosis. By leveraging bidirectional two-sample MR analysis, we identified CPXM1 (Carboxypeptidase X, M14 family member 1) as a novel gene that is causally linked to osteoporosis risk. Through transcriptomic and proteomic validation, we demonstrate that CPXM1 was upregulated in aged bone tissues and osteoporotic conditions in both human and murine models. Gene set enrichment analysis (GSEA) revealed significant dysregulation of bone homeostasis pathways, including increased extracellular matrix degradation and suppression of osteoblast differentiation in aged mice. Furthermore, phenome-wide association studies (PheWAS) confirmed minimal off-target effects of CPXM1, reinforcing its potential as a therapeutic target. Finally, computational drug repurposing predicted several promising drug candidates, including Doxorubicin, 5-Fluorouracil, and 2-Methylcholine, which may target CPXM1 pathways for osteoporosis treatment. These findings highlight CPXM1 as a potential biomarker and therapeutic target, offering new avenues for osteoporosis therapy.

Osteoporosis↗

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

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

Integration of Genome-Wide Association Studies With Single-Cell and Bulk Expression Quantitative Trait Locus to Identify Stroke Susceptibility Genes.

BACKGROUND: Previous studies have integrated genome-wide association studies with expression quantitative trait locus (eQTL) data from bulk tissues to identify stroke susceptibility genes. However, eQTL data exhibit high cell-type specificity, and genetic variants may have distinct effects across stroke subtypes. METHODS: We applied the summary-data-based Mendelian randomization (MR) method to integrate eQTL data from 7 brain cell types with genome-wide association studies data for 5 stroke phenotypes (stroke, ischemic stroke, cardioembolic stroke, large artery stroke, and small vessel stroke). Results were compared with summary-data-based MR using eQTL data from 49 tissues in the Genotype-Tissue Expression project. Robustness of significant single-cell summary-data-based MR associations was assessed via MR and colocalization analyses. Further evaluations included single-cell RNA-seq differential expression, protein-protein interaction, druggability, and phenome-wide association studies. RESULTS: Single-cell summary-data-based MR identified many novel significant genes not detected using bulk tissue eQTL data. Validated associations revealed 2 stroke risk genes (LRCH1, ICA1L), 3 stroke protective genes (AHI1, LYRM9, CENPQ), 2 large artery stroke risk genes (LIPA, ELL), and 1 ischemic stroke protective gene (CENPQ). Single-cell RNA-seq showed significantly increased LIPA expression in mouse stroke samples compared with controls. Protein-protein interaction and druggability analyses, along with phenome-wide association studies, prioritized LIPA and LRCH1 as potential therapeutic targets for stroke while indicating possible adverse effects. CONCLUSIONS: Integrating single-cell eQTL with stroke-subtype genome-wide association studies uncovers novel cell-type-specific causal genes and highlights promising therapeutic targets, advancing understanding of stroke pathogenesis.

Genome-Wide Association Study↗

PubMed2026

Spatially Contextualized Integrative Genomics Highlights Neuronal and Glial Regulatory Programs in Low Back Pain.

PURPOSE: Low back pain (LBP) is a heterogeneous pain condition with a measurable genetic contribution, but the genes, brain cell types, and spatial tissue contexts through which inherited risk is expressed remain unclear. We aimed to define cell-type-specific and spatially contextualized genetic mechanisms underlying LBP. METHODS: FinnGen R12 LBP GWAS summary statistics (42,521 cases and 353,224 controls) were integrated with brain single-nuclei eQTL data across eight major brain cell classes. We evaluated genome-wide polygenic signal using LDSC, prioritized genes using MAGMA and PoPS, and performed brain cell-type-specific eQTL-anchored Mendelian randomization, primarily based on single-instrument Wald ratio estimates, followed by Bayesian colocalization. Spatial genetic mapping was conducted using gsMap in an E16.5 murine embryonic atlas and two adult human lumbar spinal cord Visium sections. Selected candidates were assessed by RT-qPCR in neuronal-like and astroglial-like inflammatory cell models. RESULTS: LDSC supported interpretable polygenic signal for LBP. MAGMA and PoPS showed partial gene-level convergence, with TCF4 and TMEFF2 supported by both approaches. Across 1641 tested gene-cell type exposures, significant eQTL-anchored MR associations were concentrated in excitatory neurons, oligodendrocytes, inhibitory neurons, and astrocytes. Integrated eQTL-anchored MR, colocalization, and gene-prioritization evidence highlighted CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate, but not strong, colocalization support. gsMap localized LBP-associated enrichment to neuroaxis-related embryonic regions, including brain, spinal cord, sympathetic nerve, and dorsal root ganglion, and to neuronal-like niches in adult lumbar spinal cord. RT-qPCR showed model-dependent expression changes, with QPRT and LGI4 preferentially responsive in neuronal-like SH-SY5Y cells and GMPPB and DPYSL5 responsive in astroglial-like U251 cells. CONCLUSION: These findings support neuronal and glial regulatory programs as plausible contributors to LBP genetic susceptibility and highlight CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate colocalization support. The results provide a spatially contextualized framework for candidate prioritization in LBP, while emphasizing the need for larger cell-type-specific eQTL resources and functional validation before therapeutic or mechanistic conclusions can be drawn.

Mendelian randomization↗

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

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

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

Humans↗