Integrating whole genome and transcriptome sequencing to characterize the genetic architecture of isoform variation
- Chunyu Liu
- Roby Joehanes
- Jiantao Ma
- Jiuyong Xie
- Jian Yang
- Mengyao Wang
- Tianxiao Huan
- Shih-Jen Hwang
- Jia Wen
- Quan Sun
- Cumhur Y. Demirkale
- Nancy L. Heard-Costa
- Peter Orchard
- April P. Carson
- Jeffrey W. Haessler
- Laura M. Raffield
- Alex P. Reiner
- Nora Franceschini
- Paul L. Auer
- Charles Kooperberg
- Yun Li
- George O’Connor
- Joanne M. Murabito
- Peter Munson
- Daniel Levy
2025-11-22
We present a whole-blood isoform ratio QTL (irQTL) resource by analyzing genome-wide isoform-to-gene expression ratios using sequencing data. In Framingham Heart Study (FHS, n = 2622) discovery, we identify over 1.1 million cis -irQTLs (minor allele frequency [MAF] ≥ 0.01, ±1 Mb of 10,883 isoform transcripts, P < 5 × 10 −8 ) across 4,971 genes. Among 11,425 sentinel cis -irQTLs, 72% replicate ( P < 1 × 10 −4 ) in the Women’s Health Initiative (WHI; n = 2005). Notably, 20% of cis -irQTLs have no significant association with overall gene expression, indicating isoform-specific regulation. These variants are enriched at splice donor/acceptor sites and genome-wide association study loci ( P < 1 × 10 −10 ). We also identify 1870 sentinel trans -irQTLs (MAF ≥ 0.01, P < 1.5 × 10 −13 ) for 1,084 isoforms across 590 genes, and 2327 rare cis -irQTLs (0.003 < MAF < 0.01) for 2467 isoforms of 1428 genes in FHS, with external replication rates of 61% and 41% in WHI, respectively. We highlight rs12898397 in ULK3 , which alters splice site usage and reduces expression of a full-length isoform. Mendelian randomization supports a causal role between this isoform shift and reduced diastolic blood pressure. These findings highlight the power of irQTL mapping to uncover transcript-specific regulatory mechanisms underlying complex traits.