共线性分析

生信喵实验柴

发布于 2022-10-25 11:04:21

2.2K00

代码可运行

文章被收录于专栏：生信喵实验柴生信喵实验柴

运行总次数：0

代码可运行

一、共线性分析

所谓共线性，顾名思义，表示二者在一条直线上。基因组的共线性分析也类似，主要是用一种线性图的方式来比较两个或者多个基因组是否具有较好的同源性。共线性可以是核酸水平的共线性，也可以是氨基酸水平。一般氨基酸水平具有更好的同源关系，共线性也更好。

共线性这种图做起来并不难，多种比对软件都可以直接输出这种图，blast 比对可以，Mummer 和 lastz 软件也都可以。例如使用 mummer 软件包进行比对之后，使用软件包中的delta-filter 工具过滤一下 delta 格式的结果。然后直接使用 mummerplot 工具就可以直接生成共线性的比对结果了。后面一种共线性图稍微复杂一些，需要自己写程序生成。

共线性分析

二、利用 nucmer+dotPlotly 绘制共线性

#安装软件
#dotplotly
wget https://github.com/tpoorten/dotPlotly/archive/refs/heads/master.zip
#配置R包
install.packages(c("optparse", "ggplot2", "plotly"))

mkdir 50.synteny;cd 50.synteny

#利用nucmer+dotPlotly绘图
nucmer --maxmatch -l 80 -c 100 ../data/ref.fna ../data/mgh78578.fasta -p nucmer
delta-filter -r nucmer.delta >filter.delta
show-coords -c  filter.delta >filter.delta.coords
./dotPlotly-master/mummerCoordsDotPlotly.R -i filter.delta.coords -o filter.delta.plot  -m 1000 -q 300000 -k 10 -s -t -l -p 12

用Minimap2 + dotPlotly
minimap2 -x asm5 ../data/ref.fna ../data/mgh78578.fasta >minimap2.paf
./dotPlotly-master/pafCoordsDotPlotly.R -i minimap2.paf -o minimap2.plot -m 2000 -q 500000 -k 10 -s -t -l -p 12