BPGA分析細(xì)菌泛基因組

導(dǎo)讀

BPGA是Bacterial Pan Genome Analysis tool的簡寫盆偿,16年發(fā)表的工具淮野,17年最后一版更新痊远,內(nèi)置KEGG COG數(shù)據(jù)(老了)伸头,依賴usearch(32bit 可免費(fèi)用),速度很快陆盘,其他一般普筹,win linux均支持,可做參考隘马。

文獻(xiàn):BPGA- an ultra-fast pan-genome analysis pipeline. sci rep 2016
引用:293

1 下載太防,解壓,獲取依賴usearch gnuplot酸员,配置蜒车,啟動(dòng)

官網(wǎng):https://iicb.res.in/bpga/index.html,下載幔嗦,解壓酿愧,BPGA


usearch官網(wǎng):http://www.drive5.com/usearch/download.html
下載,解壓邀泉,重命名為usearch.exe嬉挡,移動(dòng)到BPGA bin文件夾,

根據(jù)BPGA User Guide呼渣,下載安裝gnuplot棘伴。

啟動(dòng)BPGA進(jìn)行初始化,正常啟動(dòng)屁置,

2 泛基因組分析 -- 默認(rèn)

準(zhǔn)備【1】> 蛋白文件【4】> 選擇文件 > 默認(rèn)分析【2】> usearch聚類 > 50%一致性 > 等待焊夸。。蓝角。

一大堆結(jié)果文件阱穗,然后,

exclusively absent genes/proteins:
orthologous families that contain genes from all genomes except one specific genome

這里列出的是每個(gè)基因組的基因分類使鹅,全部加和是遠(yuǎn)高于泛基因組基因數(shù)的揪阶。Supporting_files/pan_default.txt給出了泛基因組基因數(shù),如下患朱。不僅如此Sequences中的代表性序列的加和也是泛基因組基因數(shù)鲁僚。

泛基因組和核心基因組增長趨勢:

各基因組基因家族數(shù):

新基因數(shù)(與某一基因組相比?裁厅?冰沙?):

3 高級(jí)分析

完成后一大堆結(jié)果,
泛基因組和核心基因組执虹,又來拓挥??袋励?:

image.png

系統(tǒng)發(fā)生樹 -- 泛基因組 & 核心基因組:

KEGG注釋分類:

COG注釋分類:

實(shí)戰(zhàn):Linux中使用BPGA

獲取Linux版BPGA侥啤,獲取Linux版usearch到BPGA bin文件夾
啟動(dòng)

./BPGA-Version-1.3

基礎(chǔ)pangenome分析:
1 INPUT PREPARATION FOR CLUSTERING
2 Use any Protein Fasta files
3 enter full path to the Directory where *.fasta
4 DEFAULT PAN GENOME ANALYSIS
5 Use USEARCH Clustring Algorithm (Ultra-fast)
6 Choose Sequence Identity Cut-off for Clustering: 0.8

其他過程同window版本当叭,其實(shí)也就是輸入文件指定略有不同,似乎如此盖灸。

節(jié)點(diǎn)132G內(nèi)存蚁鳖,使用4G足以,大數(shù)據(jù)更加耗內(nèi)存

高級(jí)分析 - 進(jìn)化分析:
1 Neighbour Joining Tree (NJ):pan phylogeny
2 MLST based core phylogeny
3 Neighbour Joining Tree (NJ): core gene phylogeny

默認(rèn)僅獲得pan phylogenetic nwk赁炎,在此建樹則有core phylogenetic nwk

結(jié)果整理

out="result"
mkdir $out
mv gi_name $out
mv INPUT_all.seq $out
mv list $out
mv Results $out
mv Sequences $out
mv Supporting_files $out

更多閱讀:
BPGA - 一款泛基因組分析軟件

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末才睹,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子甘邀,更是在濱河造成了極大的恐慌,老刑警劉巖垮庐,帶你破解...
    沈念sama閱讀 216,372評(píng)論 6 498
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件松邪,死亡現(xiàn)場離奇詭異,居然都是意外死亡哨查,警方通過查閱死者的電腦和手機(jī)逗抑,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,368評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來寒亥,“玉大人邮府,你說我怎么就攤上這事「绒龋” “怎么了褂傀?”我有些...
    開封第一講書人閱讀 162,415評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長加勤。 經(jīng)常有香客問我仙辟,道長,這世上最難降的妖魔是什么鳄梅? 我笑而不...
    開封第一講書人閱讀 58,157評(píng)論 1 292
  • 正文 為了忘掉前任叠国,我火速辦了婚禮,結(jié)果婚禮上戴尸,老公的妹妹穿的比我還像新娘粟焊。我一直安慰自己,他們只是感情好孙蒙,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,171評(píng)論 6 388
  • 文/花漫 我一把揭開白布项棠。 她就那樣靜靜地躺著,像睡著了一般马篮。 火紅的嫁衣襯著肌膚如雪沾乘。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,125評(píng)論 1 297
  • 那天浑测,我揣著相機(jī)與錄音翅阵,去河邊找鬼歪玲。 笑死,一個(gè)胖子當(dāng)著我的面吹牛掷匠,可吹牛的內(nèi)容都是我干的滥崩。 我是一名探鬼主播,決...
    沈念sama閱讀 40,028評(píng)論 3 417
  • 文/蒼蘭香墨 我猛地睜開眼讹语,長吁一口氣:“原來是場噩夢(mèng)啊……” “哼钙皮!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起顽决,我...
    開封第一講書人閱讀 38,887評(píng)論 0 274
  • 序言:老撾萬榮一對(duì)情侶失蹤短条,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后才菠,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體茸时,經(jīng)...
    沈念sama閱讀 45,310評(píng)論 1 310
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,533評(píng)論 2 332
  • 正文 我和宋清朗相戀三年赋访,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了可都。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,690評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡蚓耽,死狀恐怖渠牲,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情步悠,我是刑警寧澤签杈,帶...
    沈念sama閱讀 35,411評(píng)論 5 343
  • 正文 年R本政府宣布,位于F島的核電站贤徒,受9級(jí)特大地震影響芹壕,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜接奈,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,004評(píng)論 3 325
  • 文/蒙蒙 一踢涌、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧序宦,春花似錦睁壁、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,659評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至秕噪,卻和暖如春钳降,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背腌巾。 一陣腳步聲響...
    開封第一講書人閱讀 32,812評(píng)論 1 268
  • 我被黑心中介騙來泰國打工遂填, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留铲觉,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,693評(píng)論 2 368
  • 正文 我出身青樓吓坚,卻偏偏與公主長得像撵幽,于是被迫代替她去往敵國和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子礁击,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,577評(píng)論 2 353

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