Population Health Program, QIMR Berghofer Medical Research Institute, Brisbane, QLD 4006, Australia.
Faculty of Medicine, The University of Queensland, Brisbane, QLD 4006, Australia.
HGG Adv. 2023 Feb 15;4(2):100185. doi: 10.1016/j.xhgg.2023.100185. eCollection 2023 Apr 13.
Alternative splicing contributes to cancer development. Indeed, splicing analysis of cancer genome-wide association study (GWAS) risk variants has revealed likely causal variants. To systematically assess GWAS variants for splicing effects, we developed a prioritization workflow using a combination of splicing prediction tools, alternative transcript isoforms, and splicing quantitative trait locus (sQTL) annotations. Application of this workflow to candidate causal variants from 16 endometrial cancer GWAS risk loci highlighted single-nucleotide polymorphisms (SNPs) that were predicted to upregulate alternative transcripts. For two variants, sQTL data supported the predicted impact on splicing. At the 17q11.2 locus, the protective allele for rs7502834 was associated with increased splicing of an exon in a alternative transcript encoding a truncated protein in adipose tissue and is consistent with an endometrial cancer transcriptome-wide association study (TWAS) finding in adipose tissue. Notably, haploinsufficiency is protective for obesity, a well-established risk factor for endometrial cancer. At the 17q21.32 locus, the rs2278868 risk allele was predicted to upregulate a transcript that is subject to nonsense-mediated decay, concordant with a corresponding sQTL in lymphocytes. This is consistent with a TWAS finding that indicates decreased expression in blood increases endometrial cancer risk. As SKAP1 is involved in T cell immune responses, decreased expression may impact endometrial tumor immunosurveillance. In summary, our analysis has identified potentially causal endometrial cancer GWAS risk variants with plausible biological mechanisms and provides a splicing annotation workflow to aid interpretation of other GWAS datasets.
可变剪接有助于癌症的发展。事实上,对癌症全基因组关联研究 (GWAS) 风险变异的剪接分析揭示了可能的因果变异。为了系统地评估 GWAS 变异对剪接的影响,我们开发了一种使用剪接预测工具、替代转录本异构体和剪接数量性状基因座 (sQTL) 注释的组合的优先级工作流程。将该工作流程应用于 16 个子宫内膜癌 GWAS 风险位点的候选因果变异,突出显示了预测上调替代转录本的单核苷酸多态性 (SNP)。对于两个变体,sQTL 数据支持对剪接的预测影响。在 17q11.2 位点,rs7502834 的保护性等位基因与脂肪组织中编码截断蛋白的替代转录本中一个外显子的剪接增加有关,这与脂肪组织中子宫内膜癌全转录组关联研究 (TWAS) 的发现一致。值得注意的是,杂合性缺失对肥胖是一种保护作用,肥胖是子宫内膜癌的一个已确立的风险因素。在 17q21.32 位点,rs2278868 的风险等位基因被预测上调一种受无义介导的衰变调控的 转录本,这与淋巴细胞中的相应 sQTL 一致。这与 TWAS 的发现一致,即血液中 的表达降低会增加子宫内膜癌的风险。由于 SKAP1 参与 T 细胞免疫反应, 的表达降低可能会影响子宫内膜肿瘤的免疫监视。总之,我们的分析确定了具有合理生物学机制的潜在因果子宫内膜癌 GWAS 风险变异,并提供了剪接注释工作流程,以帮助解释其他 GWAS 数据集。