画一个基因的3个可变剪切转录本示意图

Posted on by

我最近看到了一个综述，在讨论 TNFSF15 基因 ，亦称TL1A或VEGI(血管内皮细胞生长抑制因子)，为肿瘤坏死因子超家族成员之一。

  • TNFSF15在染色体上的位置是9q32，包含4个外显子区域，按照剪切位点的不同可生成3种剪切体。
  • 根据其翻译后蛋白中含有氨基酸的数目进行命名，分别为TNFSF15~174、TNFSF15~192 及TNFSF15~251。
    这个综述文章：2020,47(07):567-571. DOI:10.3971/j.issn.1000-8578.2020.19.1622
    然后我去ensembl数据库查询了一下： http://grch37.ensembl.org/Homo_sapiens/Gene/Summary?db=core;g=ENSG00000181634;r=9:117546915-117568406
    有意思的是仅仅是看到了两个转录本哦，染色体位置没有错，而且也确实是4个外显子：
    image-20210613105247088
    所以我就布置了一个学徒作业： 是画出TNFSF15基因的的3个转录本：TNFSF15~174、TNFSF15~192及TNFSF15~251的结构示意图！
    学徒给出来的代码超级简单，值得分享：

    library(ggbio) 
# hg38
library(EnsDb.Hsapiens.v86)
ensdb <- EnsDb.Hsapiens.v86
autoplot(ensdb, GeneNameFilter("TNFSF15")) + 
theme_bw()

    如下所示，也是一样的结果， 仅仅是两个转录本 ，虽然是4个外显子：
    image-20210613105433261
    那么问题来了，是数据库更新不及时吗？还是不同的数据库有不同的结果呢？

​

Comments are closed.