Nagura Yuya, Shimada Mihoko, Kuribayashi Ryoji, Ikemoto Ko, Kiyose Hiroki, Igarashi Arisa, Kaname Tadashi, Unoki Motoko, Fujimoto Akihiro
Department of Human Genetics, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan.
Department of Genome Medicine, National Centre for Child Health and Development, Tokyo, Japan.
Genome Biol. 2025 May 8;26(1):110. doi: 10.1186/s13059-025-03583-w.
Genetic variations linked to changes in gene expression are known as expression quantitative loci (eQTLs). The identification of eQTLs helps to understand the mechanisms governing gene expression. However, prior studies have primarily utilized short-read sequencing techniques, and the analysis of eQTLs on isoforms has been relatively limited.
In this study, we employ long-read sequencing technology (Oxford Nanopore) on B cells from 67 healthy Japanese individuals to explore genetic variations associated with isoform expression levels, referred to as isoform eQTLs (ieQTLs). Our analysis reveals 17,119 ieQTLs, with 70.6% remaining undetected by a gene-level analysis. Additionally, we identify ieQTLs that have significantly different effects on isoform expression levels within a gene. A functional feature analysis demonstrates a significant enrichment of ieQTLs at splice sites and specific histone marks, such as H3K36me3, H3K4me1, H3K4me3, and H3K79me2. Through an experimental validation using genome editing, we observe that a distant genomic region can modulate isoform-specific expression. Moreover, an ieQTL analysis and minigene splicing assays unveils functionally crucial variants in splicing that splicing prediction software did not assign a high prediction score. A comparison with GWAS data reveals a higher number of colocalizations between ieQTLs and GWAS findings compared to gene eQTLs.
These findings highlight the substantial contribution of ieQTLs identified through long-read analysis in our understanding of the functional implications of genetic variations and the regulatory mechanisms governing isoforms.
与基因表达变化相关的基因变异被称为表达数量性状基因座(eQTL)。eQTL的鉴定有助于理解基因表达的调控机制。然而,先前的研究主要使用短读长测序技术,对异构体上eQTL的分析相对有限。
在本研究中,我们对67名健康日本个体的B细胞采用长读长测序技术(牛津纳米孔技术),以探索与异构体表达水平相关的基因变异,即异构体eQTL(ieQTL)。我们的分析揭示了17119个ieQTL,其中70.6%在基因水平分析中未被检测到。此外,我们鉴定出对基因内异构体表达水平有显著不同影响的ieQTL。功能特征分析表明,ieQTL在剪接位点和特定组蛋白标记(如H3K36me3、H3K4me1、H3K4me3和H3K79me2)处显著富集。通过使用基因组编辑的实验验证,我们观察到一个遥远的基因组区域可以调节异构体特异性表达。此外,ieQTL分析和小基因剪接试验揭示了剪接预测软件未赋予高预测分数的剪接功能关键变体。与全基因组关联研究(GWAS)数据的比较表明,与基因eQTL相比,ieQTL与GWAS结果之间的共定位数量更多。
这些发现突出了通过长读长分析鉴定的ieQTL在我们理解基因变异的功能影响和异构体调控机制方面的重大贡献。