GEO数据挖掘(一)--简单快速下载GEO数据

之前我们已经和大家分享过如何上传GEO数据,同样,我们也可以从公共数据库下载其他研究者的数据进行分析,探寻自己感兴趣的研究方向。这种通过对现有公共数据挖掘数据进行分析的研究方式目前正是生物信息分析的热门领域。GEO数据库储存着海量的二代测序数据(特别是转录数据)和芯片数据,就像一座“藏宝山”,我们可以通过挖掘前人研究数据,提取相关研究结果,构建自身研究方向的基础架构,再辅以“湿实验”验证或针对性的其他组学测序分析,最后解决自身研究方向的复杂问题。那么接下来小诺将为大家带来全套的GEO数据挖掘系列教程,包括五步:

搜寻数据集及下载数据 差异分析--->差异表达基因 --->五大数据库注释--->蛋白互作等网络和通路注释

今天先为大家分享第一步“数据下载”,以便后续挖掘感兴趣的基因及构建通路。

首先我们可以需要先寻找贴合我们研究方向的论文所使用的数据集,下载里面的GPL文件和表达矩阵“series matrix”做基因的表达分析;再依据包含样本生存数据的临床特征数据集,去做生存性分析。

在进入GEO数据库官网前,需要先理解一下GEO的数据编号含义:

一篇文章可能包含至少一个GSE数据集,一个GSE数据集里面可能包含至少一个GSM样本。多个研究的GSM样本根据研究目的会整合为一个GDS,不过GDS运用的很少。而每个数据集都有着对应的芯片平台,就是GPL。

通过网页下载:

首先,登录GEO官网 https://www.ncbi.nlm.nih.gov/geo/,在右侧的搜索框输入一个GSE号,我们以gse21933为例,在搜索栏输入后点search进行检索。

检索后会进入Accession Display界面,里面包含着这个GSE数据的基本信息,如:标题、物种、研究概要、作者、样本描述、测序平台等等,当然还有我们最需要的原始数据。

上图可知道,这组数据包含的是肺癌与正常组织的基因表达图谱,我们想从这组数据中寻找差异基因,需要三个文件:原始文件、表型文件、注释文件。

1.原始文件,这里面储存的是每个样本中各个基因的表达量页面底部会提供原始数据,如图所示,点击http下载即可,文件是tar格式,下载下来以后需要解压缩。

2.表型文件,该文件备注每个样本是属于正常组还是癌症组,想要比较肿瘤和正常样本的区别,我们需要知道每组里面都是样本类型;

这里存储的是样本基因表达量信息;Series Matrix File,(基因表达量矩阵)

3.注释文件,因为原始数据处理得到的差异基因是以探针号表示的,我们需要知道这些探针号代表的基因是什么,这就需要用到注释文件。

有了这三类数据后,就相当于获得了 “烹饪”的原材料,之后便可以对这组数据进行“煎炒炸炖”,依据自己的需求和研究方向进行个性化的数据挖掘。

通过R包下载:

通常来说,数据挖掘之所以叫“挖掘”是因为需要在海量的数据内,通过检索收集才可以整合出可用信息从而满足我们研究目标。因此这里将介绍一项用于GEO数据下载的利器GEOquery,它是由Davis开发的一款针对下载GEO数据库的R包,目前在R开源网站bio-conductor内,运用它可以简单高效的下载所需的GEO数据。

如何使用GEOquery下载:

调用:

getGEO函数会加载GSE的matrix文件,默认会下载其注释探针信息,并对表达矩阵中的探针予以注释,但往往注释文件比较大,会出现parse保存的问题,所以一般建议把注释关掉了:getGPL=F,然后在后续分析步骤里进行手动注释。

我们下载了这些数据,就走完了第一步 “下载数据”,大家如果想了解更多后续个性化分析,请持续关注诺禾致源“GEO数据挖掘”系列文章,每周一个实用干货带您了解上手生信分析。

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

推荐阅读更多精彩内容