TCGA微卫星不稳定灶(MSI)数据处理

基本概念

微卫星(Microsatellite):基因组中的一类短串联重复DNA序列,一般由1-6个核苷酸组成,呈串联重复排列。由于其核心重复单元重复次数差异,微卫星具有群体多态性。

微卫星不稳定性(Microsatellite Instability,MSI):是指由于在DNA复制时插入或缺失突变引起的MS序列长度改变的现象,常由错配修复(MMR)功能缺陷引起。MSI现象于1993年被Jacobs等人在结直肠癌中首次发现,与癌症发生有关,可用于癌症检测。

产生原因:MSI的发生是由于肿瘤组织的DNA错配修复出现功能性缺陷导致。

这是对MSI泛癌的研究
Landscape of Microsatellite Instability Across 39 Cancer Types - PMC (nih.gov)
需要了解的知识:
1.计算MSI分数的工具:MANTIS,默认阈值0.4,高于阈值为MSI-H,低于阈值为MSS(无明显的MSI出现)。
2.最早再结直肠癌种发现,是预后良好的标志,MSI结直肠癌5年生存率要显著高于MSS结直肠癌,MSI-H结直肠癌比MSS结直肠癌有更好的预后。

数据获取

对于数据分析,首先要得到数据

之前习惯性的在UCSC 数据库中下载数据,没有发现相关msi的数据信息。搜索之后,发现有用cBioPortalData package这个包进行临床数据下载,其中有MSI的数据。

实战操作

安装

BiocManager::install("cBioPortalData")
library(cBioPortalData)

获得相应的研究信息

cbio <- cBioPortal()
studies = getStudies(cbio)
head(studies$studyId)

相应的研究信息类型在cancerTypeId这个字段,可以看一下都有什么样的研究类型

table(studies$cancerTypeId)

选择一种类型,进行下载

id = "blca_tcga_pan_can_atlas_2018"
clinical = clinicalData(cbio, id)
colnames(clinical)

我下载的是膀胱癌相关数据,可以根据cancerTypeId选择合适的研究
然后,定义MSI相关的类型

df = na.omit(clinical[,c("patientId","MSI_SCORE_MANTIS")])
colnames(df)[2] = "MSI_score"
df$MSI_score = as.numeric(df$MSI_score)
k= df$MSI_score >0.4
table(k)

发现膀胱癌相关的MSI大于0.4的并不多,可能和相应的肿瘤有一定的关系,后续研究中再深入思考一下。

后记

其实,上面已经得到了相应的数据,可以根据自己的研究,再看看有什么有趣的发现。

但是始终没有从UCSC上面得到这一部分数据,不知道是自己还没找到,还是本身就没有这一部分数据。有了解的也可以指导一下。

参考文章
100篇泛癌研究文献解读之微卫星不稳定性
cBioPortal 数据库 API 使用
TCGA的微卫星不稳定性数据获取和可视化

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,324评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,303评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,192评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,555评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,569评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,566评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,927评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,583评论 0 257
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,827评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,590评论 2 320
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,669评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,365评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,941评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,928评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,159评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,880评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,399评论 2 342

推荐阅读更多精彩内容