SEARCH · Search PubMed

Results for “eQTL analysis”

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 73 records · Page 4Linked to original sources

PubMed2026

The bioinformatics approach to identifying pathogenic variants for colorectal cancer (CRC).

Colorectal cancer (CRC) is the third most prevalent cancer globally, accounting for 9.6% of newly diagnosed cases and 9.3% of cancer-related deaths. It develops from the uncontrolled proliferation of glandular cells in the colon and rectum and is categorized into three primary types: sporadic, hereditary, and colitis-associated. While genetic susceptibility is a key factor in CRC pathogenesis, identifying high-impact pathogenic variants remains a significant challenge. This study integrates bioinformatics and population genetics approaches to identify CRC-associated single-nucleotide polymorphisms (SNPs) with potential clinical significance. CRC-associated SNPs were extracted from the Genome-Wide Association Studies (GWAS) Catalog, functionally annotated via HaploReg, and validated via Ensembl. In addition, expression quantitative trait locus (eQTL) data from the GTEx database were used to assess the effects of these variants on gene expression across human tissues. Our analysis identified three high-priority SNPs (rs9379084, rs3184504, and rs11557154) associated with the RREB1, ATXN2, SH2B3, and DCAF12 genes, which exhibited marked allele frequency differences among populations. These findings suggest potential biomarkers for CRC risk assessment and highlight the importance of genetic screening across diverse populations.

Bioinformatics↗

PubMed2026

Multi-omics Mendelian Randomization Prioritizes Neutrophil Extracellular Trap-related Genes Associated with Atrial Fibrillation Risk.

BACKGROUND: Neutrophil extracellular traps (NETs) participate in thrombosis, inflammation, and cardiovascular remodeling, yet whether NET-related genes (NRGs) are associated with atrial fibrillation (AF) risk across multiple molecular layers remains unclear. This study used a multiomics Mendelian randomization framework to prioritize NRGs supported by methylation, expression, and protein quantitative trait loci (QTL) data. METHODS: Genome-wide significant cis instruments (P < 5 × 10-8) were obtained for 90 methylation QTLs (mQTLs), 100 expression QTLs (eQTLs), and 38 protein QTLs (pQTLs) mapped to 137 literature- curated NRG entries. Summary-data-based Mendelian randomization (SMR) coupled with the heterogeneity in dependent instruments (HEIDI) test was applied using whole-blood mQTL data (n = 1,980), eQTLGen blood eQTL data (n = 31,684), and deCODE plasma pQTL data (n = 35,559). AF outcome data were obtained from a meta-analysis including 60,620 cases and 970,216 controls of European ancestry. RESULTS: At the methylation level, 21 CpG-feature associations across 13 genes remained significant after HEIDI filtering and false discovery rate (FDR) correction. Expression-level analysis identified eight significant gene-AF associations, whereas protein-level analysis identified seven significant features representing five unique proteins. Cross-omics integration prioritized C3, MAPK3, and STAT3 as Tier 1 genes, CTSC, LPAR3, and THBD as Tier 2 genes, and fourteen additional genes as Tier 3 candidates. C3 showed risk-increasing protein-level associations together with multiple significant CpG signals, whereas MAPK3 and STAT3 showed directionally protective expression/protein or methylation/protein patterns. DISCUSSION: The cross-omics convergence on C3, MAPK3, and STAT3 is consistent with complement activation, immune-fibrotic signaling, and cytokine-regulatory pathways implicated in AF biology, but the findings should be interpreted as genetic prioritization rather than definitive intervention-ready causality. CpG-level heterogeneity at the C3 locus and the blood/plasma origin of the QTL resources further support a cautious interpretation. Modest colocalization support and the unresolved possibility of pQTL sample overlap further support this cautious, hypothesis-generating interpretation. CONCLUSION: Multi-omics SMR prioritizes C3, MAPK3, and STAT3 as the most consistently supported NET-related genes associated with AF risk. These findings provide a framework for atrialtissue replication and mechanistic validation of NET-related pathways in AF.

Atrial fibrillation↗

PubMed2026

Network-based stratification of allele-specific expression reveals patient subgroups in Huntington's disease.

MOTIVATION: Huntington's disease (HD) exhibits substantial variability in age of onset and disease progression that is not fully explained by CAG repeat length alone. Part of this residual variation is heritable, implicating additional genetic mechanisms. cis-regulatory variation, genetic variants that alter transcription and splicing of nearby genes, represents one such mechanism that can be quantified through allele-specific expression (ASE) analysis. However, methods for integrating ASE profiles into patient stratification frameworks remain underdeveloped, particularly for rare diseases with small cohorts and sparse data. RESULTS: We adapt a network-based stratification algorithm, originally developed for somatic tumour mutations, to ASE data. By propagating gene-level ASE imbalance profiles through a protein-protein interaction network, we stratified 20 HD patients into three distinct biological patient subgroups. Differential gene expression analysis highlights neuroinflammatory pathways, including microglial activation, immune cell activation, and cytokine regulation, as key sources of inter-patient heterogeneity, while differential ASE analysis implicates proteasomal and ubiquitin-dependent protein catabolic processes, immune activation, and central nervous system development. Intersection of differentially imbalanced and expressed genes identified FAM181B as a candidate gene with potential eQTL-mediated regulation, supported by independent cis-eQTL evidence for rs3780 in the caudate and putamen, the primary HD-affected striatal regions. FAM181B encodes a nuclear protein expressed in neural tissues acting as an interactor of the Hippo pathway TEAD transcription factors, implicating transcriptional regulatory variation as a potential contributor to molecular heterogeneity between patient subgroups. Differences in cortical and striatal neuropathological scores between clusters, even when adjusted for CAG repeat length, provide clinical support for the biological relevance of the identified subgroups. AVAILABILITY: All analysis code, Docker containers, and conda environments are available at https://github.com/macsbio/HD-ASE-NBS.

Huntington Disease↗

PubMed2026

Integrative Multi-Omics Analysis Identifies Thrombosis-Associated Molecular Features Linked to Germline Susceptibility and Immune Cell Communication in Gastric Cancer.

Emerging evidence indicates that coagulation-related molecular programs are associated with thrombosis, tumor progression, and molecular dysregulation in gastric cancer (GC). However, thrombosis-associated molecular features in GC and their potential links to inherited susceptibility remain insufficiently understood. Integrated analyses of transcriptomic data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were performed to identify thrombosis-associated genes and establish a machine learning-based prognostic signature. Genome-wide association study (GWAS), expression quantitative trait loci (eQTL), transcriptome-wide association study (TWAS), and Mendelian randomization (MR) analyses were conducted to investigate susceptibility-associated transcriptional programs in GC. Functional assays were used to evaluate candidate genes associated with malignant phenotypes. Single-cell RNA sequencing (scRNA-seq) and cell-cell communication analyses were further performed to characterize cell-type-specific expression patterns and potential intercellular interactions. A total of 22 differentially expressed thrombosis-associated genes were identified, and a prognostic signature comprising 14 genes was established. The signature stratified patients into high- and low-risk groups and showed prognostic performance in both the training and validation cohorts. Integrative GWAS, eQTL, and TWAS analyses identified susceptibility-associated transcriptional programs that were positively correlated with the thrombosis-associated risk score. Silencing ACTN2 and CRYAB significantly reduced GC cell migration and invasion. scRNA-seq analysis revealed relatively high CRYAB expression in neutrophils, and CellChat analysis suggested potential neutrophil-B cell interactions involving COLLAGEN-related signaling. This integrative multi-omics study identified a thrombosis-associated molecular signature linked to prognosis and germline susceptibility-associated transcriptional programs in GC. ACTN2 and CRYAB may represent candidate genes associated with GC cell migration and invasion, while single-cell analysis suggested potential immune-related communication features.

Humans↗

PubMed2004

Quantitative genetic variation: a post-modern view.

It has become commonplace to map individual quantitative trait loci (QTL) in experimental organisms; the means (line-crosses and dense maps of markers) and motivation (the close relationship between continuous physiological traits and common, complex diseases) are self-evident. Progress in mapping human QTL has been more gradual, an inevitable consequence of genetic mapping in a natural population setting. The common objective of these studies has been to understand the molecular mechanisms underlying individual QTL. Recent theoretical and practical advances shift this focus to a more comprehensive or genomic perspective on quantitative variation. Fisher's infinitesimal model of adaptive evolution, which satisfied quantitative geneticists for over 50 years, has been modified in the light of data from QTL mapping experiments in plants and animals. The resulting exponential model provides a pleasing empirical fit to the distribution of QTL effect sizes, predicts that a large amount of quantitative variation will be explained by a limited number of genes and suggests a new mathematical framework for linkage mapping. Molecular analysis of QTL suggests that coding variants (e.g. allozymes) underlie a fraction of quantitative variation and that variants that affect gene expression (expression QTL, eQTL) have a substantial role. This is supported by genomic experiments that combine expression profiling with classical genetic mapping approaches to reveal a remarkable wealth of quantitative heritable variation in the transcriptome and that cis-and trans-acting regulatory factors are organized in networks reflecting pleiotropy. It is hoped that these advances will enhance our understanding of the genetic basis of complex inherited diseases.

Animals↗

PubMed2026

Predicting gene-specific regulation with transcriptomic and epigenetic single-cell data.

MOTIVATION: Analysis of single cell ATAC-seq and RNA-seq data has allowed to gain unprecedented insights into gene regulation by allowing to define cell type-specific regulatory regions and their effects on gene expression. While powerful, such analysis is challenging due to the inherent sparsity of single cell data. RESULTS: We present a new approach, MetaFR, to learn gene-specific models that link open-chromatin variation from scATAC-seq data to gene expression from scRNA-seq. Using efficient regression trees, we illustrate that accurate expression prediction models can be learned on the single-cell or meta-cell level. Validation was done using fine-mapped eQTLs. Meta-cell models were found to outperform single-cell models for most genes. Comparison to the SOTA method SCARlink revealed advantages of MetaFR in terms of runtime and prediction performance. MetaFR thus allows time-efficient analysis and obtains reliable models of gene expression prediction, which can be used to study gene regulation in any organism for which scRNA-seq and scATAC-seq data is available. AVAILABILITY AND IMPLEMENTATION: MetaFR is available under https://github.com/SchulzLab/MetaFR.

Single-Cell Analysis↗

PubMed2026

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

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

CEBPZ↗

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↗

PubMed2024

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

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

Humans↗

PubMed2006

Differential gene expression of wheat progeny with contrasting levels of transpiration efficiency.

High water use efficiency or transpiration efficiency (TE) in wheat is a desirable physiological trait for increasing grain yield under water-limited environments. The identification of genes associated with this trait would facilitate the selection for genotypes with higher TE using molecular markers. We performed an expression profiling (microarray) analysis of approximately 16,000 unique wheat ESTs to identify genes that were differentially expressed between wheat progeny lines with contrasting TE levels from a cross between Quarrion (high TE) and Genaro 81 (low TE). We also conducted a second microarray analysis to identify genes responsive to drought stress in wheat leaves. Ninety-three genes that were differentially expressed between high and low TE progeny lines were identified. One fifth of these genes were markedly responsive to drought stress. Several potential growth-related regulatory genes, which were down-regulated by drought, were expressed at a higher level in the high TE lines than the low TE lines and are potentially associated with a biomass production component of the Quarrion-derived high TE trait. Eighteen of the TE differentially expressed genes were further analysed using quantitative RT-PCR on a separate set of plant samples from those used for microarray analysis. The expression levels of 11 of the 18 genes were positively correlated with the high TE trait, measured as carbon isotope discrimination (Delta(13)C). These data indicate that some of these TE differentially expressed genes are candidates for investigating processes that underlie the high TE trait or for use as expression quantitative trait loci (eQTLs) for TE.

Carbon Isotopes↗

PubMed2025

Complex Genetics and Regulatory Drivers of Hypermobile Ehlers-Danlos Syndrome: Insights from Genome-Wide Association Study Meta-analysis.

BACKGROUND: Hypermobile Ehlers-Danlos syndrome (hEDS) is the most common subtype of EDS, a group of heritable connective tissue disorders. Clinically, hEDS is defined by generalized joint hypermobility and chronic musculoskeletal pain, but its impact extends beyond the musculoskeletal system. Affected individuals frequently experience autonomic, gastrointestinal, immune, and neuropsychiatric involvement, highlighting both the multisystemic nature of the condition and challenges of diagnosis. In contrast to other EDS subtypes with defined genetic causes, the molecular basis of hEDS has remained elusive. METHODS: We conducted a genome-wide association study (GWAS) of hEDS across three case controls studies, including 1,815 cases and 5,008 ancestry-matched controls. Fixed-effects meta-analysis of 6.2 million variants was complemented with LDAK gene-based association testing, transcriptome-wide association studies, and integrative annotation across multiple tissues and cell types including eQTLs, enhancer marks and open chromatin accessibility profiles, supported by luciferase assays on one candidate variant. LD-score genetic correlations were assessed between hEDS and 19 frequently reported comorbid conditions. RESULTS: Two loci reached genome-wide significance, including a regulatory region near the atypical chemokine receptor 3 gene (ACKR3) on chromosome 2. Functional annotation supports ACKR3 risk alleles colocalize with eQTLs in tibial nerve, alter enhancer activity, and generate a de novo AHR transcription factor regulatory site, implicating neuroimmune and pain signaling pathways. Gene-based and transcriptome-wide analyses identified common variants in a locus containing multiple candidates, including SLC39A13, a zinc transporter critical for connective tissue development previously implicated in a rare form of EDS, and PSMC3, a gene involved in central nervous system development. LD-score regression revealed significant genetic correlations between hEDS and joint hypermobility, myalgic encephalomyelitis/chronic fatigue syndrome, fibromyalgia, depression, anxiety, autism spectrum disorder, migraine, and gastrointestinal diseases. CONCLUSIONS: These results establish the first evidence of common variant contributions to hEDS, supporting a complex, multisystem model involving neuroimmune-stromal dysregulation. Our findings add novel indications to hEDS pathogenesis and provide solid foundations for future molecular definition and therapeutic discovery.

Genome-wide association study↗

PubMed2026

Mitochondrial Function-Related Genes in Sleep Disorders: A Multi-Omics Mendelian Randomization Study.

Mitochondrial dysfunction is linked to sleep disorders in previous report, but the potential roles of specific genes remain unclear. This study aimed to dissect different subtype-specific genetic associations and their underlying mechanisms. A multi-omics Summary-data-based Mendelian Randomization (SMR) approach was performed to identify potential causal links between mitochondrial function-related genes and sleep disorders. We integrated GWAS data from FinnGen database (the discovery set), independent GWAS datasets (covering different sleep-disorder subtypes and used for validation), and cis-QTLs (including mQTLs, eQTLs, and pQTLs) to perform systematic exploration. Specially, we performed targeted validation of tissue-specific effects, leveraging gene expression data from disease-relevant brain regions within the GTEx database. Our SMR analysis identified mitochondrial function-related genes potentially modulating sleep disorders across biological layers, initially identifying 102 genes at the methylation level, 48 at the gene expression level, and 6 at the protein abundance level. Integrative analysis subsequently prioritized DCXR and ACADVL and revealed their distinct, subtype-specific associations. DCXR exhibited a protective role in sleep apnea while ACADVL showed a paradoxical risk conferring role in daytime sleepiness. In addition, the analysis identified an epigenetic regulatory mechanism for DCXR in which its expression and protein levels are modulated by DNA methylation. Finally, validation in brain-hypothalamus tissue confirmed DCXR as a significant potential protective factor (OR = 0.929, 95% CI: 0.887-0.973, P_HEIDI = 0.999, FDR = 0.2449). Our findings implicate key mitochondrial genes, particularly DCXR and ACADVL, in the pathophysiology of specific sleep disorder subtypes, highlighting potential avenues for precision medicine. Clinical trial number: Not applicable.

Humans↗

PubMed2026

GFPT1 as a cross-ancestry validated target for degenerative spinal disease: genetic association in a Chinese cohort and functional characterization in zebrafish.

Degenerative spinal disease (DSD), including spinal stenosis and spondylosis, lacks effective pharmacological treatment. To identify druggable targets and assess cross-ancestry applicability, we integrate multi-omics analyses using Summary-data-based Mendelian Randomization (SMR), colocalization, and two-sample Mendelian randomization with European whole-blood, peripheral-blood, and CSF eQTL/pQTL datasets, followed by whole-genome sequencing (WGS) validation in a Chinese cohort. We identify 7 genes/proteins associated with spinal stenosis and 5 with spondylosis, with GFPT1, GPX1, and SERPINA1 shared by both. Two-sample MR further supports the causal associations of these targets with DSD. Phenome-wide MR prioritization selects GFPT1 and GPX1 as favorable candidates with no predicted adverse effects and potential beneficial effects on hypertension. In the Chinese cohort (67 lumbar spinal stenosis patients and 100 controls), WGS identifies 4 GFPT1 cis-eQTL loci (rs13016371, rs35392088, rs12997521, and rs13019789) associated with lumbar spinal stenosis risk; all risk alleles are linked to increased GFPT1 expression, and all 24 variant carriers show L4/L5 stenosis on imaging. Druggability analysis identifies IOX1 as the sole preclinical-stage compound targeting GFPT1, and molecular docking supports robust binding to GFPT1 (- 6.39 kcal/mol). Functional assays show that IOX1 directly inhibits GFPT1 enzymatic activity and induces fructose-6-phosphate accumulation. In zebrafish, IOX1 significantly rescues GFPT1-induced degenerative phenotypes. These findings establish GFPT1 as a cross-ancestry validated therapeutic target for DSD and nominate IOX1 as a promising disease-modifying candidate.

Animals↗

PubMed2025

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

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

Humans↗

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↗

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↗

PubMed2006

Population and systems genetics analyses of cortisol in pigs divergently selected for stress.

This study presents a systems genetic analysis on the physiology of cortisol in mice and pigs with an aim to show the potential of a comprehensive computational approach to quickly identify candidate genes and avoid a costly whole-genome quantitative trait locus (QTL) mapping. Population genetics analyses were performed on measurements of cortisol from a pig selection experiment. Expression QTL were mapped and gene networks were built using gene expressions for Crhr1 (corticotrophin-releasing hormone receptor) gene and single nucleotide polymorphisms from public mouse data. Results from mouse data were used to infer potential candidate regulatory genes involved in pig cortisol regulation, using a comparative or translational systems genetics approach. The pig data used were from a 10-yr divergent genetic selection experiment, providing data on 417 individuals. Population genetics analysis showed that cortisol is highly genetically determined with heritabilities of 0.40-0.70. Furthermore, a major gene with an additive effect of 86 ng/ml is segregating. Genetical-genomics investigations revealed two trans-acting eQTL for Crhr1 gene expression on chromosomes 2 and 13. Candidate gene search under trans-eQTL peaks yielded 63 genes for Crhr1 expression phenotypes. Functional links for Crhr1 genes with other genes/proteins in the gene network using mouse data were shown for the first 10 statistically significant genes involved. Results show translational or comparative systems genetics approaches reduce costs and time in large-scale genetics and "-omics" investigations. This is the first study to report a strong genetic basis for cortisol physiology using a systems approach.

Animals↗

PubMed2026

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

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

Mendelian randomization↗