R語(yǔ)言里做基因組共線性可視化R包~GENESPACE

論文

GENESPACE: syntenic pan-genome annotations for eukaryotes

https://www.biorxiv.org/content/10.1101/2022.03.09.483468v1

還沒(méi)有發(fā)表

github主頁(yè)

https://github.com/jtlovell/GENESPACE

詳細(xì)介紹

https://htmlpreview.github.io/?https://github.com/jtlovell/GENESPACE/blob/master/doc/genespaceOverview.html

windows系統(tǒng)還不能用 只能在MacOS或者在Linux系統(tǒng)下使用,我試試在linux下使用

首先安裝orthofinder

conda install -c bioconda orthofinder 

安裝MCScanX

https://github.com/wyp1125/MCScanX

git clone https://github.com/wyp1125/MCScanX.git
cd MCScanX
make
image.png

這里出現(xiàn)了三個(gè)error,但是也出現(xiàn)了三個(gè)可執(zhí)行程序撮胧,試了一下可以運(yùn)行叔锐,不知道后面會(huì)不會(huì)有影響

image.png

安裝依賴的R包

conda install r-data.table r-dbscan r-R.utils r-devtools
conda install bioconductor-Biostrings bioconductor-rtracklayer

安裝GENESPAE

# 啟動(dòng)R radian
devtools::install_github("jtlovell/GENESPACE", upgrade = F)

運(yùn)行示例數(shù)據(jù)

library(GENESPACE)
runwd<-file.path("./testGenespace/")
make_exampleDataDir(writeDir = runwd) ## 這一步會(huì)下載示例數(shù)據(jù)

gids<-c("human","chimp","rhesus") 
gpar<-init_genespace(genomeIDs = gids,speciesIDs = gids,versionIDs = gids,ploidy = rep(1,3),wd = runwd,gffString = "gff",pepString = "pep",path2orthofinder = "orthofinder",path2mcscanx = "/home/myan/scratch/apps/mingyan/Biotools/MCScanX",path2diamond = "diamond",diamondMode = "fast",orthofinderMethod = "fast",rawGenomeDir = file.path(runwd,"rawGenomes")) 

parse_annotations(gsParam = gpar,gffEntryType = "gene",gffIdColumn ="locus",gffStripText = "locus=",headerEntryIndex = 1,headerSep = " ",headerStripText = "locus=") 
# 上面這行代碼沒(méi)有看懂是在干啥

gpar<-run_orthofinder(gsParam = gpar)  

## 運(yùn)行這行代碼出現(xiàn)警告信息
Warning message:
In system2(gsParam$paths$orthofinderCall, com, stdout = TRUE, stderr = TRUE) :
  running command ''orthofinder' -b ./testGenespace//orthofinder -t 4 -a 1 -X -og 2>&1' had status 120 and error message 'Interrupted system call'
## 不知道時(shí)候?qū)罄m(xù)有影響 有可能是 runwd<-file.path("./testGenespace/") 這行代碼最后多了一個(gè)斜線 重新運(yùn)行了一遍沒(méi)有問(wèn)題了

gpar<-synteny(gsParam = gpar)

## 畫(huà)圖展示

pdf(file="abc.pdf",width = 10,height = 8)
plot_riparianHits(gpar)
dev.off()
image.png

畫(huà)圖更多的參數(shù)

pdf(file="abc.pdf",width = 9.6,height = 4)
plot_riparianHits(gpar, refGenome = "chimp",invertTheseChrs = data.frame(genome = "rhesus", chr = 2),genomeIDs = c("chimp", "human", "rhesus"),labelTheseGenomes = c("chimp", "rhesus"),gapProp = .001,refChrCols = c("#BC4F43", "#F67243"),blackBg = FALSE,returnSourceData = T, verbose = F)
dev.off()
image.png

還可以自定義感興趣的區(qū)域

regs <- data.frame(genome = c("human", "human", "chimp", "rhesus"),chr = c(3, 3, 4, 5),start = c(0, 50e6, 0, 60e6),end = c(10e6, 70e6, 50e6, 90e6),cols = c("pink", "gold", "cyan", "dodgerblue"))
 pdf(file = "abc2.pdf",width = 9.6,height = 4)
plot_riparianHits(gpar, onlyTheseRegions = regs,blackBg = FALSE)
dev.off()

image.png

構(gòu)建泛基因組組

pg <- pangenome(gpar)

輸出一個(gè)文件 results/human_pangenomeDB.txt.gz

打開(kāi)這個(gè)文件,部分結(jié)果如下

image.png

這個(gè)結(jié)果怎么看暫時(shí)沒(méi)看懂

幫助文檔里寫(xiě)道

This is the source data that can be manipulated programatically to extract your regions of interest. Future GENESPACE releases will have auxilary functions that let the user access the pan-genome by rules (e.g. contains these genes, in these regions etc.). For now, we’ll leave this work to scripting by the user.

接下來(lái)就是研究研究如何準(zhǔn)備自己的數(shù)據(jù)

歡迎大家關(guān)注我的公眾號(hào)

小明的數(shù)據(jù)分析筆記本

小明的數(shù)據(jù)分析筆記本 公眾號(hào) 主要分享:1置侍、R語(yǔ)言和python做數(shù)據(jù)分析和數(shù)據(jù)可視化的簡(jiǎn)單小例子冲簿;2晤柄、園藝植物相關(guān)轉(zhuǎn)錄組學(xué)擦剑、基因組學(xué)、群體遺傳學(xué)文獻(xiàn)閱讀筆記芥颈;3惠勒、生物信息學(xué)入門學(xué)習(xí)資料及自己的學(xué)習(xí)筆記!

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末爬坑,一起剝皮案震驚了整個(gè)濱河市纠屋,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌妇垢,老刑警劉巖巾遭,帶你破解...
    沈念sama閱讀 212,718評(píng)論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異闯估,居然都是意外死亡灼舍,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,683評(píng)論 3 385
  • 文/潘曉璐 我一進(jìn)店門涨薪,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)骑素,“玉大人,你說(shuō)我怎么就攤上這事刚夺∠壮螅” “怎么了末捣?”我有些...
    開(kāi)封第一講書(shū)人閱讀 158,207評(píng)論 0 348
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)创橄。 經(jīng)常有香客問(wèn)我箩做,道長(zhǎng),這世上最難降的妖魔是什么妥畏? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 56,755評(píng)論 1 284
  • 正文 為了忘掉前任邦邦,我火速辦了婚禮,結(jié)果婚禮上醉蚁,老公的妹妹穿的比我還像新娘燃辖。我一直安慰自己,他們只是感情好网棍,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,862評(píng)論 6 386
  • 文/花漫 我一把揭開(kāi)白布黔龟。 她就那樣靜靜地躺著,像睡著了一般滥玷。 火紅的嫁衣襯著肌膚如雪氏身。 梳的紋絲不亂的頭發(fā)上,一...
    開(kāi)封第一講書(shū)人閱讀 50,050評(píng)論 1 291
  • 那天惑畴,我揣著相機(jī)與錄音观谦,去河邊找鬼。 笑死桨菜,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的捉偏。 我是一名探鬼主播倒得,決...
    沈念sama閱讀 39,136評(píng)論 3 410
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼夭禽!你這毒婦竟也來(lái)了霞掺?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 37,882評(píng)論 0 268
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤讹躯,失蹤者是張志新(化名)和其女友劉穎菩彬,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體潮梯,經(jīng)...
    沈念sama閱讀 44,330評(píng)論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡骗灶,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,651評(píng)論 2 327
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了秉馏。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片耙旦。...
    茶點(diǎn)故事閱讀 38,789評(píng)論 1 341
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖萝究,靈堂內(nèi)的尸體忽然破棺而出免都,到底是詐尸還是另有隱情锉罐,我是刑警寧澤,帶...
    沈念sama閱讀 34,477評(píng)論 4 333
  • 正文 年R本政府宣布绕娘,位于F島的核電站脓规,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏险领。R本人自食惡果不足惜侨舆,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 40,135評(píng)論 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望舷暮。 院中可真熱鬧态罪,春花似錦、人聲如沸下面。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,864評(píng)論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)沥割。三九已至耗啦,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間机杜,已是汗流浹背帜讲。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 32,099評(píng)論 1 267
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留椒拗,地道東北人似将。 一個(gè)月前我還...
    沈念sama閱讀 46,598評(píng)論 2 362
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像蚀苛,于是被迫代替她去往敵國(guó)和親在验。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,697評(píng)論 2 351

推薦閱讀更多精彩內(nèi)容