用R繪制兩個物種間染色體共線性
2021-12-26 20:16 作者:Albert_Xin | 我要投稿
一、加載R包
1. 安裝R包[1]
2. 加載R包:
二、 準備數(shù)據(jù)
a. 染色體信息數(shù)據(jù)見圖1,第一列數(shù)據(jù)為染色體名稱,第二列為染色體開始位置,第三列為染色體結(jié)束為止,第三列為染色體填充顏色,第五列為物種名稱,第六列是物種名稱在圖片中顯示的大小,第七列是物種名稱的顏色,支持RGB和十六進制的顏色。

圖1 染色體信息數(shù)據(jù)
b.共線性數(shù)據(jù),具體格式見圖 2,第一列為圖一中物種排名第一的數(shù)據(jù)如此處為Grape,且第一列為染色體ID,跟隨其后的兩行為所在染色體為止的信息,第二個物種就是圖一中排序在后的物種的染色體ID,跟隨兩行為染色體位置信息,最后一列為填充的顏色,根據(jù)你感興趣的基因,設置彩色,并將其排序到最后以方便在圖片中顯示。

圖 2 共線性數(shù)據(jù)
三、運行R包
ideogram包中karyotype為染色體信息數(shù)據(jù),synteny為共線性數(shù)據(jù),output為svg圖片輸出路徑,結(jié)果如圖 3。

[1] Hao Z, Lv D, Ge Y, et al. RIdeogram: drawing SVG graphics to visualize and map genome-wide data on the idiograms[J]. PeerJ Computer Science, 2020, 6: e251.
標簽: