基因预测总结
- 格式:docx
- 大小:19.79 KB
- 文档页数:7
基因预测总结
1、基因预测
对于真菌来说有四个ab initio预测软件:
GlimmerHMM,SNAP,Genearkes,augustus 以及同源预测(homology)。四个软件中:GeneMarkes是通过隐马模型工作的,但是它不需要参考物种,是自身训练的,不需要参考序列,当处理一个新物种,没有理想的或者较近缘的已测序物种时可以采用这种方法。Augustus,GlimmerHMM,SNAP都需要参考训练集的。
总流程:
perl /nas/MG01/FUNGUS/PGAP/FGAP.pl [options]
Genome.fa
Options
--all run all analysis for Fungi
--cutlen cut the scaffolds longer than this
--predict select the method to predict genes:
augustus,genemarkes,snap,glimmerhmm or homology
--prepara set the parament for augustus,snap,homology
--repeat set repeat method, defalut: repbase-proteinmasker-trf
--ncRNA set ncRNA type, default: tRNA-rRNA-miRNA-sRNA-snRNA
--rRNA_ref set Reference for rRNA, if null rRNA will be
predicted by rRNAmmer
--function set dbs for gene function annotaion,default:
nr-swissprot-trembl-cog-kegg-iprscan
--lib set the lib for synteny analysis and gene family analysis,
needed
--synteny synteny analysis
--family Gene Family analysi --species species tree, default, created by lib information
--category category file, default, created by lib information
--cpu set the cpu number to use in parallel, default 20 for
qsub and 5 for multi --run set the parallel type, qsub, or multi,
default=qsub
--outdir set the result directory, default="."
--prefix set a prefix name for results
--help output help information to screen
分步流程程序路径:/nas/MG01/FUNGUS/PGAP/gene-prediction/bin/gene-predict.pl
perl gene-predict.pl [options]
--glimmer run glimmer by self training
--genemark run genemark by self training
--shape set the shape of prokaryote DNA,
circular,linear,partial, default=partial --glimmerhmm run
glimmerhmm and give a glimmerhmm parameter directory
--snap run snap and give a snap parameter file
--genemarkes run genemarkes by self traning
--augustus run augustus and set species
--homology predict genes based on proteins on a homology
species
--genemarkM run genemarkM for mata gene prediction
--metagene run metagene for meta gene prediction
--metageneA run metageneA for meta gene prediction
--cpu set the cpu number to use in parallel, default=3
--run set the parallel type, qsub, or multi, default=qsub
--prefix set gene id prefix
--outdir set the result directory, default="./"
--verbose output running progress information to screen
--help output help information to screen 1.1Genemarkes预测:
Self-training algorithm GeneMark-ES
a) splits input sequence at such "NN...N" strings
b) runs gene finding GeneMark.hmm on contigs
c) maps back predictions to original super-contig sequence
As a result, incomplete gene structures can be predicted inside
super-contig sequences.
Script:perl ./gene-predict.pl --genemarkes
GeneMarkES 输出结果为./genemark_hmm.gtf
1.2 Homology预测
Homology(同源预测)是通过基因组序列和参考蛋白集进行比对来确定基因位置的,预测的结果特点是基因数目少,但是准确率很高。通过genewise软件预测,用此方法需要提供近缘物种的氨基酸序列,由于同源预测的结果比真实的少,所以建议大家再用此方法时不要只找一个参考物种,多找几个近缘物种,尽量能多预测出一些基因来。
perl ./gene-predict.pl --homology ../input/CNGH2S.fa.pep -prefix RHO
*.pep 序列为参考序列的蛋白序列。
1. 3 SNAP预测
SNAP也是通过隐马尔科夫模型工作的,自身是没有现成的训练集的,需要参考物种进
行训练集的构建。如果要使用这个两个软件,必须寻找一个参考物种,得到它的基因组序列和基因位置信息的gff文件,自己来构建训练集来进行预测
SNAP基因预测包括两步:
Step 1:参考物种训练集的构建。训练结果为*.hmm文件。
perl /nas/MG01/FUNGUS/PGAP/gene-prediction/bin/train_snap.pl --name gram
为参考基因组序列
为基因位置信息的gff文件 Step 2 :
perl ./gene-predict.pl --snap .. /SNAP/train/gram.hmm --prefix RHO
1. 4 Augustus预测
Augustus运用隐马尔科夫模型,隐马模型ab initi o预测。基本要素包括两个状态(观察状态、隐含状态)和三个概率(初始概率、转移概率和两态对应概率)。模型在DNA序列和基因结构上定义了一个概率分布,采用的是维特比算法。Augustu s自带了真菌的48个训练集,大家可以选择较近缘的物种来进行预测,当自带的训练集中没有和目标基因组近缘或者亲缘关系都特别远不太理想时,这时可以自己寻找一个参考物种,得到它的基因组序列和基因位置信息的gff文件,自己来构建训练集来进行预测。
第一种情况:如果自带的训练集中有和目标基因组近缘,可以用Augustus自带的训练集进行预测。
perl ./gene-predict.pl –augustus
saccharomyces_cerevisiae_S288C
说明:saccharomyces_cerevisiae_S288C 为augustus自带46个训练集其中的一个标签。Augustus训练集路径为:/nas/MG01/liule/bin/software/augustus.2.5/config/。
第二种情况:如果自带的训练集中没有和目标基因组近缘或者亲缘关系都特别远不太理想时,可以自己寻找一个参考物种,得到它的基因组序列和基因位置信息的gff文件,自己来构建训练集来进行预测。自己构建的训练结果会添加到Augustus自带的训练集中以—species
命名。
perl
/nas/MG01/liule/bin/software/augustus.2.5/scripts/autoAugTrain.pl --genome=/path/ref.seq --trainingset=/path/ref.gff --species =Curvularia_lunata –useexisting
--genome :近源物种的参考序列的绝对路径
--trainingset :参考序列的gff 文件的绝对路径 --species :物种名字,自己命名。
--useexisting 如果它们存在species,使用和改变目前的配置和参数文件。
1.1.5 GlimmerHMM预测
GlimmerHMM,自身是没有现成的训练集的,需要参考物种进行训练集的构建。如果要使用这个软件,必须寻找一个近源性非常高的参考物种,得到它的基因组序列和基因位置信息的gff文件,自己来构建训练集来进行预测。
Step1 :参考物种训练集的构建。
Script :/nas/MG01/FUNGUS/PGAP/software/GlimmerHMM/train/trainGlimmerHMM
-d train
参考序列得基因组文件,FASTA格式
is a file with the exon coordinates relative to the sequences
contained in the ; different genes are separated by a
blank line; I am assuming a format like below:
seq1 5 15
seq1 20 34
seq1 50 48