0%

物种内(一个物种)与物种间(两个物种)共线性分析

物种内(一个物种)与物种间(两个物种)共线性分析

使用的相同的原理

1 准备好注释文件 gff 和“基因组”文件 fasta

2 使用一键 MCScanX,得到一大堆文件

(时间慢,有快捷方法)

3从结果中找出几个文件

.collinearity

.gff(是比对得到的文件,不是最开始你准备的gff)

.ctl (仅在物种间比对时产生,物种内不会产生)

4 物种间分析,直接导入上述文件即可

5 物种内分析

准备文件:

①打开 fasta stats,统计各染色体长度,用excel提取,另存为Chrlen.txt

②删去用不到的contig

③打开file merger for MCScanX,输入第三步得到的.collinearity,改变mode为collinear,键入输出文件名为GenePair.tab.txt

④打开File Transformat for MicroSynteny Viewer,输入GenePair.tab.txt和第三步得到的.gff,键入输出名为Link_Region.tab.txt

补充基因家族成员:

⑤打开text block extractor,输入GenePair.tab.txt,键入输出名为GenePair.family.tab.txt,输入你的基因家族成员的ID号

⑥回到File Transformat for MicroSynteny Viewer,输入GenePair.family.tab.txt,gff文件(第三步结果中的),键入输出名为Link_Region_family.tab.txt,用excel打开,在最后输入颜色的RGB值R如红色255,0,0

⑦将改好的excel全选,插入到Link_Region.tab.txt的最前边,另存为Link_Region_family.tab.red.txt

补充标签:

⑧打开table row manipulator,导入第三步得到的gff,选中基因ID所在列,导入ID号,键入输出名为family_genpos.tab.txt

⑨打开adavanced circos,键入文件

Chrlen.txt,

family_genpos.tab.txt,

Link_Region_family.tab.red.txt

修改与导出图片:

⑩修改各部分大小,字体,颜色,间距,并导出图像为pdf,svg(可在AI内修改)

谢谢你们对我的支持!

欢迎关注我的其它发布渠道