基因组预测和注释软件怎样检查基因组注释结果的可靠性

基因组注释主要包括四个研究方向：重复序列的识别；非编码RNA的预测；基因结构预测和基因功能注释。我们将分别对这四个领域进行阐述。

1：重复序列的识别。

重复序列的研究背景和意义：重复序列可分为串联重复序列（Tendam repeat）和散在重复序列(Interpersed repeat)两大类。其中串联重复序列包括有微卫星序列，小卫星序列等等；散在重复序列又称转座子元件，包括以DNA-DNA方式转座的DNA转座子和反转录转座子(retrotransposon)。常见的反转录转座子类别有LTR,LINE和SINE等。

重复序列识别的发展现状：目前，识别重复序列和转座子的方法为序列比对和从头预测两类。序列比对方法一般采用Repeatmasker软件，识别与已知重复序列相似的序列，并对其进行分类。常用Repbase重复序列数据库。从头预测方法则是利用重复序列或转座子自身的序列或结构特征构建从头预测算法或软件对序列进行识别。从头预测方法的优点在于能够根据转座子元件自身的结构特征进行预测，不依赖于已有的转座子数据库，能够发现未知的转座子元件。常见的从头预测方法有Recon，Piler，Repeatscout,LTR-finder，ReAS等等。

重复序列识别的研究内容：获得组装好的基因组序列后，我们首先预测基因组中的重复序列和转座子元件。一方面，我们采用RepeatScout、LTR-finder、Tendem Repeat Finder、Repeatmoderler、Piler等从头预测软件预测重复序列。为了获得从头预测方法得到的重复序列的类别信息，我们把这些序列与Repbase数据库比对，将能够归类的重复序列进行分类。另一方面，我们利用Repeatmasker识别与已知重复序列相似的重复序列或蛋白质序列。通过构建Repbase数据库在DNA水平和蛋白质水平的重复序列，Repeatmasker能够分别识别在DNA水平和蛋白质水平重复的序列，提高了识别率。

重复序列识别的关键技术难点：

1）：第二代测序技术测基因组，有成本低、速度快等优点。但是由于目前产生的读长（reads）较短。由于基因组序列采用kmer算法进行组装，高度相似的重复序列可能会被压缩到一起，影响对后续的重复序列识别。

2）：某些高度重复的序列用现有的组装方法难以组装出来，成为未组装reads（unassembled reads）。有必要同时分析未组装reads以得到更为完整的重复序列分布图。之前，华大已开发了ReAS软件，专门用于识别未组装reads中的重复序列。但该软件目前只能处理传统测序技术(如sanger测序)生成的较长片段的reads，需要进一步改进方可用于分析第二代测序技术得到的reads。同时，未组装的短片段reads重复度更高，识别其重复区域具有较大难度。

重复序列识别的研究方向：

1）：整合现有的重复序列预测方法，对组装好的基因组序列进行分析。

2）：综合考虑并结合短序列组装策略，校正重复序列识别的结果。

3）：开发识别未组装reads重复序列的算法和流程并构建一致性序列。

2：非编码RNA序列的预测。

非编码RNA预测的研究背景和意义：非编码RNA，指的是不被翻译成蛋白质的RNA，如tRNA, rRNA等，这些RNA不被翻译成蛋白质，但是具有重要的生物学功能。miRNA结合其靶向基因的mRNA序列结合，将mRNA降解或抑制其翻译成蛋白质，具有沉默基因的功能。tRNA (转运RNA)携带氨基酸进入核糖体，使之在mRNA指导下合成蛋白质。rRNA(核糖体RNA)与蛋白质结合形成核糖体，其功能是作为mRNA的支架，提供mRNA翻译成蛋白质的场所。snRNA（小核RNA）主要参与RNA前体的加工过程，是RNA剪切体的主要成分。

非编码RNA预测的发展现状：由于ncRNA种类繁多，特征各异，缺少编码蛋白质的基因所具有的典型特征，现有的ncRNA预测软件一般专注于搜索单一种类的ncRNA，如tRNAScan-SE 搜索tRNA、snoScan 搜索带C/D盒的snoRNAs、SnoGps 搜索带H/ACA 盒的snoRNAs、mirScan 搜索microRNA等等。Sanger实验室开发了Infernal软件，建立了1600多个RNA家族，并对每个家族建立了一致性二级结构和协方差模型，形成了Rfam数据库。采用Rfam数据库中的每个RNA的协方差模型，结合Infernal软件可以预测出已有RNA家族的新成员。Rfam/Infernal方法应用广泛，可以预测各种RNA家族成员，但是特异性较差。我们建议：如果有更好的专门预测某一类非编码RNA的软件，那么采用该软件进行预测；否则，使用Rfam/Infernal流程。

非编码RNA预测的研究内容：利用Rfam家族的协方差模型，我们采用Rfam自带的Infernal软件预测miRNA和snRNA序列。由于rRNA的保守性很强，为此我们用序列比对已知的rRNA序列，识别基因组中的rRNA序列。tRNAscan-SE工具中综合了多个识别和分析程序，通过分析启动子元件的保守序列模式、tRNA二级结构的分析、转录控制元件分析和除去绝大多数假阳性的筛选过程，据称能识别99%的真tRNA基因。

非编码RNA预测中拟解决的关键技术难点：

识别非编码RNA的假基因：基因组中很多序列由非编码RNA基因复制而来，与非编码RNA基因序列相似，但不具有非编码RNA的功能。目前我们采用的非编码RNA序列的预测方法都是基于序列比对和结构预测，不能够很好的去除这类非编码RNA的假基因。针对这个问题，我们考虑结合RNA表达信息如RNA-seq数据进行筛选。

非编码RNA预测的研究方向：

1）：专门检测小片段RNA序列的方法现在已经得到广泛应用，利用小片段RNA序列数据进行非编码RNA的预测是我们的重要研究方向。

2）：开发miRNA靶向基因预测流程：miRNA通过调控其靶向基因的mRNA稳定性或翻译来控制生命活动的进程。预测miRNA靶向基因能够给我们研究miRNA功能带来提示。由于miRNA在动物和植物中对靶向基因的调控机制差别较大，我们建议对动物和植物分别建立靶向基因预测流程，提高预测准确度。

3：基因结构预测。

基因结构预测的研究背景和意义：通过基因结构预测，我们能够获得基因组详细的基因分布和结构信息，也将为功能注释和进化分析工作提供重要的原料。基因结构预测包括预测基因组中的基因位点、开放性阅读框架（ORF）、翻译起始位点和终止位点、内含子和外显子区域、启动子、可变剪切位点以及蛋白质编码序列等等。

基因结构预测的发展现状：原核生物基因的各种信号位点（如启动子和终止子信号位点）特异性较强且容易识别，因此相应的基因预测方法已经基本成熟。Glimmer是应用最为广泛的原核生物基因结构预测软件，准确度高。而真核生物的基因预测工作的难度则大为增加。首先，真核生物中的启动子和终止子等信号位点更为复杂，难以识别。其次，真核生物中广泛存在可变剪切现象，使外显子和内含子的定位更为困难。因此，预测真核生物的基因结构需要运用更为复杂的算法，常用的有隐马尔科夫模型等。常用的软件有Genscan、SNAP、GeneMark、Twinscan等。

基因结构预测的研究内容：基因结构预测主要通过序列比对结合从头预测方法进行。序列比对方法采用blat和pasa等比对方法，将基因组序列与外部数据进行比对，以找到可能的基因位置信息。常用的数据包括物种自身或其近缘物种的蛋白质序列、EST序列、全长cDNA序列、unigene序列等等。这种方法对数据的依赖性很高，并且在选择数据的同时要充分考虑到物种之间的亲缘关系和进化距离。基因从头预测方法则是通过搜索基因组中的重要信号位点进行的。常用的软件有Genscan、SNAP、Augustus、Glimmer、GlimmerHMM等等。同时采用多种方法进行基因预测将产生众多结果，因此最后需要对结果进行整合以得到基因的一致性序列。常用软件有Glean，EVM等。

基因结构预测中拟解决的关键技术难点：

目前，真核生物的基因结构预测方法仍有较大改进空间，主要面临以下的技术难点。

1）：如何利用现有的数据和算法，更好地识别基因的可变性剪切位点。

2）：随着测序工作的进展，许多目前研究较少的物种也将提上测序日程。大多基因结构的从头预测算法需要预先训练预测参数。现有资源和数据稀缺的物种将很难获得预测参数。

3）：克服组装错误对基因结果预测的影响

4）：建立基因结构预测的评价系统。

可变性剪切位点的预测较为困难。如何结合RNA-seq数据进行可变剪切预测将是重要的工作方向和难点。

基因结构预测的研究方向：

1）：利用RNA-seq、EST等数据校正基因结构预测结果，识别可变剪切位点。

2）：对于研究较少的物种，建议利用近缘物种的同源基因数据以训练基因结构预测软件。

3）：利用同源基因组之间的共线性信息，辅助基因结构预测。

4：基因功能注释。

基因功能注释的研究背景和意义：获得基因结构信息后，我们希望能够进一步获得基因的功能信息。基因功能注释方向包括预测基因中的模序和结构域、蛋白质的功能和所在的生物学通路等。

基因功能注释的发展现状：全基因组测序将产生大量数据，而实验方法由于成本较高，不适用于全基因组测序的后续功能分析。为此，目前普遍采用比对方法对全基因组测序的基因功能进行注释。KEGG和Gene Ontology是目前使用最为广泛的蛋白质功能数据库，分别对蛋白质的生物学通路和功能进行注释。Interpro通过整合多个记录蛋白质特征的数据库，根据蛋白质序列或结构中的特征对蛋白质进行分类。

基因功能注释的研究内容：目前，我们利用四个常用的数据库进行基因功能注释。使用的数据库有Uniprot蛋白质序列数据库、KEGG生物学通路数据库、Interpro蛋白质家族数据库和Gene Ontology基因功能注释数据库。

1）：与Uniprot蛋白质序列数据库比对，获得序列的初步信息。

2）：与KEGG数据库比对，预测蛋白质可能具有的生物学通路信息。

3）：与Interpro数据库比对将获得蛋白质的保守性序列，模序和结构域等。

4）：预测蛋白质的功能。Interpro进一步建立了与Gene Ontology的交互系统：Interpro2GO。该系统记录了每个蛋白质家族与Gene Ontology中的功能节点的对应关系，我们通过此系统便能预测蛋白质执行的生物学功能。

基因功能注释中拟解决的关键技术难点：

目前我们的功能注释工作是建立在比对的基础上，这将会带来两个比较大的问题。首先，此方法严重依赖于外部数据，对某些研究较少的物种限制很大。其次，序列相似并不表示实际生物学功能相似，考虑引入序列比对之外的方法，进一步完善基因功能注释工作。

基因功能注释的研究方向：考虑引入序列比对之外的数据（如蛋白质互作网络、基因表达谱等），利用概率模型算法进行整合，完善基因功能注释工作。

基因组注释分析主要包括哪些内容

基因组注释包括以下方面的内容：

（1）重复序列的预测。通过比对已知的重复序列数据库，找出序列中包含的重复序列，识别类型并转化为N或者X，统计各种类型重复序列的分布。

（2）编码基因的预测。通过将转录组或EST数据比对到拼接后的基因组序列上，找出编码基因位置，预测编码基因结构。或者通过专业的外显子预测软件，预测编码基因的外显子结构。

（3）小RNA基因的预测。通过比对已知的小RNA的数据库，或者通过生物信息(bioinformation)学软件预测，找出这些小RNA基因，并进行分类。

（4）调控序列和假基因的预测。

基因功能的注释，使用的数据库包括NT/NR, SwissProt/TrEMbl, InterPro, KEGG, COG, Gene ontology等，使用比对的方法，如blast，找出同源相近的基因，并注释功能。

如何使用基因注释软件dogma

基因组注释系统是MGAP 的核心,整合了许多常用的基因识别和蛋白质功能预测软件,包括GeneMarks、IPRsearch、BLASTPGP 和FASTA3 等,以及多个数据库,如非冗余蛋白质序列数据库(Non redundant , NR) 、已知三维空间结构的蛋白质序列数据库(PDBSeq) 、国际蛋白质资源信息系统( InterPro) [6] 和直系同源蛋白质家族数据库(Cluster of orthologousgroups ,COG) 等,编写了相应的模块进行自动操作,并把每一步注释结果导入数据库中。MGAP 整合的一般模块,可以被其他任何一种微生物基因组直接使用。不同实验室可根据实际研究需要,增加相应模块或数据,如蓝细菌Anabaena sp. strain PCC 7120 的蛋白质序列库等。

基因组注释不完整可以做全转录组吗

组学omics,研究的是整体.按照分析目标不同主要分为基因组学,转录组学,蛋白质组学,代谢组学.基因组学研究的主要是基因组DNA,使用方法目前以二代测序为主,将基因组拆成小片段后再用生物信息学算法进行迭代组装.当然这仅仅是第一步,随后还有繁琐的基因注释等数据分析工作.转录组学研究的是某个时间点的mRNA总和,可以用芯片,也可以用测序.芯片是用已知的基因探针,测序则有可能发现新的mRNA,蛋白组学针对的是全体蛋白,组要以2D-Gel和质谱为主,分为top-down和bottom-up分析方法.理念和基因组类似,将蛋白用特定的物料化学手段分解成小肽段,在通过质量反推蛋白序列,最后进行搜索,标识已知未知的蛋白序列.代谢组分析的代谢产物,是大分子和小分子的混合物,主要也是用液相和质谱.总而言之,这些技术都想从全局找变量,都是一种top-down的研究方法,原因很简单：避免‘只缘身在此山中’的尴尬.但因为技术局限,都各有缺点,尤其是转录组和蛋白组数据,基本上颠覆了以前一直认为的mRNA水平能代表蛋白水平的观念,因为这两组数据的重合度太低.所以目前很多研究都开始使用交叉验证方法.

在基因组中有多个注释到某一基因怎么区分

1、把EST用BLAT比对到基因组序列上，挑最好的match。

2、下载同版本的基因组注释文件。

3、比较1和2中的基因组位置关系，并找出来未被基因组注释的EST。

这时候应该还剩下上千条，其中绝大部分都是蛋白编码基因的反义RNA，当然一小部分是基因间区的非编码RNA。

4.、将affimetrix的probe sequence用blast比对到这些非编码RNA中。建议先找找关心的事情是否有tiling-array data。

另外要注意非编码RNA和它可能有关的基因一起分析，比如附近基因，或者互补的基因。

做芯片一般不去找表达或非表达，而往往是在找差异表达。

比如某非编码RNA和它附近的基因共同在肿瘤细胞中特异表达。

如何从宏基因组序列中预测出特定的基因序列

把你的参考基因多下载一些不同的，然后随便揪个mapping软件跑跑不就好了？

另一个技术路线是：先组装，然后再预测cds，然后注释cds。好处是组装了之后数据量会少非常多，坏处是组装得未必对。

最重要的是：这根本不是你的真实问题。你的真实问题是：你不会搞高通量测序分析，但却必须要做。

refseq和ucsc人类基因组注释文件哪个质量更高

UCSC Genome Browser是由University of California Santa Cruz (UCSC) 创立和维护的，该站点包含有人类、小鼠和大鼠等多个物种的基因组草图，并提供一系列的网页分析工具。站点用户可以通过它可靠和迅速地浏览基因组的任何一部分，并且同时可以得到与该部分有关的基因组注释信息，如已知基因，预测基因，表达序列标签，信使RNA，CpG岛，克隆组装间隙和重叠，染色体带型，小鼠同源性等。用户也可以因为教育或科研目的加上他们自己的注释信息。UCSC Genome Browser目前应用相当广泛，比如Ensembl 就是使用它的人类基因组序列草图为基础的。

如何选择合适的基因组结构注释文件

在线的基因注释资源对于基因组数据的分析来说必不可少注释就是对代码的解释和说明。目的是为了让别人和自己很容易看懂。基因组，Genome，一般的定义是单倍体细胞中的全套染色体为一个基因组，或是单倍体细胞中的全部基因为一个基因组。

转载请注明出处51数据库 » 基因组预测和注释软件怎样检查基因组注释结果的可靠性