TBtools基因家族分析详细教程(3)基因家族成员的进化分析2

  • 基因-共线性的定义与常见算法原理
  • 物种内的共线性分析
    文件准备(物种比对到自身的.blast文件,物种基因信息文件.gff文件),运行MCScanX,输出collinear和tandem文件
  • 基因家族成员的来源分析(如何复制得到)
  • 不同物种之间的共线性分析
  • 共线性分析结果可视化

1 共线性分析:与同线性的联系

用途:

  • 识别直系同源gene
  • 蛋白编码基因注释
  • 发现进化事件

2物种内的共线性分析

3基因家族来源分析

4不同物种之间的共线性分析

共线性分析

数据文件下载genome.fa,gff3,protein.fa

2数据文件格式转换(TBtools)

3共线性分析

4解读文本输出结果

-----开始----

1 下载菠萝,水稻,拟南芥,香蕉的基因组和注释文件

  • 通过TBtools由上述文件得到CDS和protein文件(前面已讲)
  • 把菠萝蛋白比对到自身(用时相对较长)得到blast结果文件

2 获得所有基因的位置信息


如下



下面可以把刚才得到的blast结果文件简化,也可以不做,做的话,下面

3 菠萝自身的比对的结果如下

image.png

可视化

先得到串联重复序列的link文件
上面得到的.tandem文件用excel打开并进行分列,另存为txt文件

image.png

结果
GRAS基因家族在染色体上的位置并显示串联重复序列

可以看到有串联重复序列
再把pineapple2pineapple.blast.tab.collinearity文件转换为link文件
获得基因间关系的link文件

结果如下
image.png

可以看出和视频中不一样,因为我和作者用的不是同一个基因家族
对于比对到自身的(单个基因组)的还可以做其他的
设置

结果

也可以选择性展示


circle gene view.300dpi.jpg

4不同物种之间的共线性分析

分析菠萝与水稻之间的共线性区块

  • 需要菠萝的所有蛋白序列比对到水稻的所有蛋白序列
  • 两个基因组的所有基因的位置关系

按前述步骤分别得到水稻的CDS和protein,方法不再赘述
这里需要说明的是,视频中CDS的总序列数为66338,我下载了几个水稻品种包括reference geonome均不是,就用以下信息吧


image.png

取最长的可变剪切本,以使下一步分析更加准确



接下来,这一步看具体请看需要不需要做

################################

下面开始两个基因组比较

开始菠萝和水稻比对

image.png

然后水稻比对菠萝

都比对完之后,开始merge两个比对后的blast文件


image.png

同样gff文件也要merge

image.png

然后



得到pineapple_rice.collinearity文件
然后


mutiple synteney plotter

mutiple synteney plotter

新建一个multiple文件夹

接下来在做菠萝和香蕉的比对
步骤按上面
提取cds,pro(考虑可变剪切,可以选择最大长度可变剪切序列),然后互相比对得到blast结果
上面个gff和blast结果分别merge,就可以比对了

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

推荐阅读更多精彩内容

  • MCScanX是检测基因共线性和进化分析的常用工具之一,2012发表至今引用数200+,作者之一的唐海宝老师是国内...
    rapunzel0103阅读 70,308评论 21 102
  • 写在前面:非常感谢CJ开发出如此强大的TBtools工具,作者详细讲解了基因家族的分析过程和意义,课程购买地址,讲...
    Y大宽阅读 160,585评论 18 390
  • 基因组组装完成后,或者是完成了草图,就不可避免遇到一个问题,需要对基因组序列进行注释。注释之前首先得构建基因模型,...
    xuzhougeng阅读 50,311评论 11 184
  • 到了25岁,满世界的被催婚,好像再不结婚,就是不断的在贬值,再过一两年,都可以放在超市的特价栏,和稀烂的苹果香蕉一...
    木枭阅读 202评论 0 0
  • 你问我身在何方,为何独自飞翔? 我想长路漫漫,谁能陪我永远流浪? 轻风抚过我的手掌,雨水吻干我的泪痕。 清寂艰苦的...
    陌雅安阅读 155评论 3 5