2023-05-22比较基因组分析大致流程

1. 下机数据的质检,组装。基因组survey可以预测基因组大小。

(三代测序长,但不准,可以用二代协助组装。还可以测HiC进行染色体级别的组装。不同的测序手段选择不同的组装软件)

二代组装策略:

小片段 insert size:170bp;500bp; 800bp
大片段 insert size:1kb;2kb;5kb;10kb;20kb;
主流组装工具
ABySS
IDBA
SOAPdenovo2 动植物基因组组装常用且高效
Velvet
SPAdes 小基因组(<100Mb)组装时的首选
Minia 速度快,内存资源最省
MaSuRCA
plantnus 适合高杂合度组装

可能见到的名词

de Bruijn Graph(DBG)算法
Overlap-Layout-Consensus(OLC)算法

三代组装策略

测序方法:Pacific Biosciences (PacBio) Single Molecule Real Time (SMRT) 和Oxford Nanopore Technology (ONT)
纠错
组装
polish提升准确度

二、三代混合组装策略

二代数据组装得到contig
三代长读长序列延长contig,提升组装效果
提升重复序列注释效果

评价指标

contig N50 使得累加后长度超过组装总长度一半的contig的长度就是N50
对组装效果的评价主要依由于据组装序列的连续性、完整性和准确性。——3C原则:
连续性(Contiguity): 得到的contig要足够的长
正确性(Correctness): 组装的contig错误率要低
完整性(Completeness):尽可能包含整个原始序列
组装评价工具:BUSCO。(conda即可安装)

2. 注释

可分为重复序列注释,基因结构注释,基因功能注释,ncRNA注释

重复序列注释

denovo:TRF, LTR-finder / LTR-harvest / LTR-STRUC,Repeatmodeler …
Homlog-based: RepeatMasker (Repbase),RepeatProteinMasker

基因结构注释

这个就很多了 也是分denovo(augustus等)和homlog(genewise等)。
然后多种方法得到的基因集用EVM软件合并。有fasta文件和gff文件就可以得到基因序列和基因蛋白序列。BUSCO检查基因集完整度。

基因功能注释

GO,KEGG等,在线即可完成。需要pep文件。

ncRNA注释

我还没做过 不清楚

3.有基因集后,可分析

WGD鉴定,有基因集就能做。种间互相比,和已知WGD次数的模式物种比。生成一个曲线图。
共线性。(MSCcanX)没有染色体可能不好比

4.树是基本,不对的话就要调整。

orthofinder可以提单拷贝同源基因(物种选择不要太多)也可以生成一个初级的树。orthofinder可以找到单拷贝同源基因(近缘物种间一般要达到200个基因左右 树才是准确性较高的)
单拷贝同源基因构树:多序列比对(muscle)---->提取保守序列(Gblocks)---->选择模型(prottest)---->使用RAxML构建系统发育树。
基因家族收缩扩张(CAFE5)在树的基础上去推算分化时间(r8s,timetree,mcmctree,iqtree)后做。

5.特色基因的分析

多看文献,找到相关基因家族(如MADS-box)在拟南芥(或者水稻)中相关的基因序列,然后blast比对,找到自己物种的基因序列,然后看数量,结构,序列构成的变化,进化选择压力。进行详细分析。

345可以同时进行。互不影响

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

推荐阅读更多精彩内容