如何使用SPDE进行基因家族分析---蛋白结构域、启动子元件、保守结构域

在进化树建完后,需要对序列的结构进行分析。这里需要区分两个概念,即蛋白结构域与保守结构域,它们两个不是一回事。从上一章,通过pfam方法同学们可以获得蛋白结构域的信息(即结构域的名字,从哪儿开始到哪儿结束),而保守结构域,它指的是在你所输入的这些序列中保守的序列是哪些。至于这些保守的序列到底能不能到蛋白结构域(即可以预测有功能的结构域)的那个层次还不好说。

保守结构域通过meme网站(MEME - Submission form (meme-suite.org)):


输入DNA或者蛋白文件

选择要检测出多少保守结构域


选择text格式:


接下来开始绘图:

需要基因id以及每个基因的长度(即输入到meme网站的那些~):


输入文件即是用于meme分析的那些

需要格式化meme的结果同时也需要这些基因的ID, ID的获取:



格式化meme:


这些文件都是自动完成,大家只需要点点点就好:

绘图:

点击draw

这样就完成了保守结构域的绘制。且每单击一次default colors, 颜色都会自动配置一次,然后再度点击draw即可完成颜色的更改。其他参数参照之前的说明。上一张图吧:


在上一章讲到的pfam的结果可用于提取蛋白结构域。建议同学们:如果你的文件里,基因名字过于复杂,就是含有额外信息,建议大家利用格式化的功能将这些额外去除掉,因为这些额外信息在后续的分析中基本不会用到但如果保留它们则会对后续的一些分析产生干扰。


格式化pfam文件

这里要注意一个问题,如果所研究的家族有多个结构域,需要依次对它们的pfam结果进行格式化。格式化的结果类似于这样:


格式化的结果

之后将该结果合并在一起(就是直接复制粘贴)以及上一步的基因ID+长度的那个文件输入作图模块中进行作图即可。

再就是启动子元件的展示。需要提取启动子序列,功能在这里:


执行此功能,同学们需要注意一个问题,就是你蛋白序列文件中的ID可能与gff3文件里的ID有出入,需要稍微调整,例如:


如上图

通过比较可以发现蛋白文件中的多了“.p” 而gff文件中的多了“.v1.ABR2.1”,因此需要进行替换(就是一定要保证蛋白文件中的ID与gff文件中mRNA那一行对应的ID是一样的),替换功能在这里:


当然,同学们如果不习惯用SPDE也可以用其他软件例如word替换,这里设置的替换功能,真正的用处是在执行多个关键词同时替换的时候,这种单个关键词的替换,完全展现不出威力~。替换完成后,如下:


这个文件将作为提取启动子的关键词。对基因组序列文件进行格式化:

提取的启动子如下:

将启动子放到plantcare(PlantCARE, a database of plant promoters and their cis-acting regulatory elements (ugent.be))进行启动子序列的分析。


这里有个点儿需要注意:一个基因的启动子上有很多元件,如果将这些元件都展示在启动子上会将图整的跟食品包装袋后的条码似的,这并不利于我们展示重点想要展示的东西,因此,这里需要同学们根据自己的实验目的进行删减,例如,如果同学们想要展示的是胁迫条件下的情况,则可以将一些无关元件(如TATA等)删除,只保留重点先要展示的,这样也使得我们的研究结果一目了然。具体操作步骤与之前相同:


结果:


启动子元件

接下来是展示内含子和外显子的分布情况,当然还是要先处理文件,需要的是基因ID和GFF文件:


得到的文件如下:


之后画图

嗯,大致是这样子的吧

同学们在做的时候如果基因数量多,大可考虑将这些基因分成几组展示~

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

推荐阅读更多精彩内容