https://www.bilibili.com/video/av85069377?p=2

大致流程

一：準(zhǔn)備文件

1. 擬南芥的基因組，CDS，蛋白泳炉，GFF以及HMM文件。

GCF_000001735.4_TAIR10.1_genomic.fna
GCF_000001735.4_TAIR10.1_genomic.gff
GCF_000001735.4_TAIR10.1_protein.faa
GCF_000001735.4_TAIR10.1_rna.fna
NB-ARC.hmm

2.利用HMM文件搜尋pep文件

#利用hmmsearch搜尋pep
$hmmsearch --cut_tc --domtblout NBS-ABC.out NBS-ARC.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa
#過(guò)濾篩選得到E-value小于1*10-20,先拿到序列號(hào)
$grep -v "#" NBS-ABC.out|awk '($7 + 0) < 1E-20'|cut -f1 -d  " "|sort -u > NBS-ARC_qua_id.txt
#再根據(jù)篩選的序列號(hào)嚎杨，從擬南芥pep文件中提取序列
$less Arabidopsis_thaliana.TAIR10.pep.all.fa | /home/spider/project/yuantao/soft/seqkit grep -f NBS-ARC_qua_id.txt > NBS-ARC_qua.fa

hmmsearch結(jié)果文件

3.利用clustalw進(jìn)行多序列比對(duì)花鹅，構(gòu)建物種特定的NB-ARC基因家族hmm模型

流程如下：

clustalw -> 1 -> NB-ARC_qua.aln ->2 ->1--> NBS-ARC_qua.aln ->NB-ARC_qua.dnd ->X

最終生成NBS-ARC_qua.aln和NBS-ARC_qua.aln文件，其中NB.aln文件就是用于構(gòu)建新的馬爾科夫模型文件枫浙。

$hmmbuild new-NB-ARC.hmm NBS-ARC_qua.aln

在利用新生成的擬南芥特異的馬爾科夫模型尋找pep文件中的NB-ARC基因家族序列刨肃。

$hmmsearch --cut_tc --domtblout NBS-ARC.second.out new-NB-ARC.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa

對(duì)結(jié)果文件篩選，E-value值小于0.01箩帚，ID取重復(fù)重復(fù)真友，最后將ID存為final_E02.NB_ids.list。

grep -v "#" NBS-ARC.second.out|awk '($7 + 0) < 1E-20' | cut -f1 -d " "|sort -u >final.NBS.list
$less Arabidopsis_thaliana.TAIR10.pep.all.fa | /home/spider/project/yuantao/soft/seqkit grep -f final.NBS.list > final_NBS-ARC_qua.fa

這個(gè)序列就可以用于后續(xù)的進(jìn)化樹(shù)紧帕，保守結(jié)構(gòu)域盔然，motif等分析。

二：通過(guò)blast鑒定基因家族成員（適用于沒(méi)有該基因家族HMM模型）

利用CDD，Pfam愈案，interProscan確定基因家族成員保守結(jié)構(gòu)域挺尾。
利用Plant CARE和PLACE預(yù)測(cè)成員順式作用原件。
利用別的物種的基因家族的蛋白序列刻帚，實(shí)例用玉米和水稻的NBS序列尋找擬南芥的NBS序列。

1.從NCBI下載NBS序列（盡可能下載多個(gè)其他物種的NBS序列）

取NCBI得protein條目下搜索NBS序列涩嚣。

2.blastp比對(duì)并篩選目標(biāo)物種中符合要求的序列崇众。

用query=ref_nbs.plant.fa去擬南芥蛋白序列中搜索

$makeblastdb -in Arabidopsis_thaliana.TAIR10.pep.all.fa -dbtype prot -out TAIR
$blastp -num_threads 16 -db TAIR -query ref_nbs.plant.fa -evalue 1e-20 -outfmt 7 -seg yes > reftoTAIR_blastp.out
#刪除#開(kāi)頭的行
$sed -i '/^#/d' reftoTAIR_blastp.out
$awk '{print$2}' reftoTAIR_blastp.out >NBS_second.ids
##最后利用sort和uniq函數(shù)對(duì)NBS_second.ids去重復(fù)，再將去重復(fù)的HMM和blast找到的id取交集航厚。
$cat NBS_second.ids|sort|uniq >uniq_NBS_second.ids
$comm -12 uniq_NBS_second.ids final.NBS.list > common.list
#最后根據(jù)交集中的ID取擬南芥的pep文件中搜索可信的NB-ARC序列
$less Arabidopsis_thaliana.TAIR10.pep.all.fa | /home/spider/project/yuantao/soft/seqkit grep -f common.list > final_all_NBS-ARC_qua.fa

拿到這些蛋白序列之后顷歌，就需要到NCBI，Pfam等網(wǎng)站上搜索是否含有NBS保守的結(jié)構(gòu)域幔睬。

3.取Pfam查找NBS結(jié)構(gòu)域（max=500條序列）

Pfam網(wǎng)站

例子

結(jié)果

4.取NCBI的的CDD

不支持文件上傳赦抖，只能復(fù)制粘貼。

CDD

搜索

5.利用interPro搜索結(jié)構(gòu)域

interpro

結(jié)果