[发明专利]一种利用二代和三代转录组测序数据的基因组注释方法

说明书

本发明涉及基因组注释技术领域，提供一种利用二代和三代转录组测序数据的基因组注释方法。本发明的方法包括下述步骤：步骤1：将三代全长转录组测序序列比对到目标基因组，得到各编码基因的初始结构信息；步骤2：将二代转录组测序短序列比对到目标基因组，从比对文件中提取内含子剪切位点信息；步骤3：结合各编码基因的初始结构信息和内含子剪切位点信息，获得各编码基因的最终结构信息。本发明既能充分利用三代转录组无须拼接即可得到全长转录本序列从而准确性高的优点，又能充分考虑二代转录组测序数据能够提供大量的内含子剪切位点证据的优点，大大提高了基因组注释的准确性和效率。

技术领域

本发明涉及基因组注释技术领域，涉及一种真核生物全基因组编码基因结构的注释方法，特别是涉及一种利用二代和三代转录组测序数据的基因组注释方法。

背景技术

一般全基因组序列组装完毕之后，需要对其蛋白编码基因的结构进行预测。预测通常综合采用三种策略：从头预测、基于近缘物种的序列同源性进行预测以及基于转录组数据的预测。由于转录组是一个物种所表达基因的最直接的证据，因此基于转录组数据的预测结果被认为是最可信的，在整合所有策略的预测结果得到最终基因集的时候该策略所给的权重一般也是最高的。

基于转录组数据预测编码基因的方法中，比较有代表性的有PASA、tophat+cufflinks、hisat2+stringtie等。二代和三代转录组测序数据都可以作为PASA的输入。三代全长转录组因为无须组装即可得到转录本的全长序列，具有准确可靠的优点。通常三代数据可以直接输入给PASA，二代数据则需先用组装软件如Trinity组装后才能输入。而tophat+cufflinks和hisat2+stringtie则比较相似，都是基于二代测序的预测策略，前者可以将测序的短序列比对到基因组上，而后者可以将比对结果信息提取出来从而达到预测编码基因的目的。这两种策略的优点是能充分利用二代测序数据提供的信息，但不能直接利用三代测序数据。PASA的缺点则是不能充分利用二代测序数据所能提供的大量内含子位置信息，且速度较慢。

可见，现有的基因组注释方法各有优缺点，但都无法同时利用二代和三代测序的优点，对基因组预测的精度和效率都不够高。而且随着测序技术的进步，使得每个物种可以利用的二代和三代转录组数据越来越多，如何充分利用这些数据得到更为可靠的预测结果成为迫切需要解决的难题。

发明内容

针对现有技术存在的问题，本发明提供一种利用二代和三代转录组测序数据的基因组注释方法，能够同时利用二代和三代转录组测序数据的优点，大大提高基因组注释的准确性和效率。

本发明的技术方案为：

一种利用二代和三代转录组测序数据的基因组注释方法，其特征在于，包括下述步骤：

步骤1：将三代全长转录组测序序列比对到目标基因组，得到各编码基因的初始结构信息；

步骤2：将二代转录组测序短序列比对到目标基因组，从比对文件中提取内含子剪切位点信息；

步骤3：结合各编码基因的初始结构信息和内含子剪切位点信息，获得各编码基因的最终结构信息。

进一步的，所述步骤1包括下述步骤：

步骤1.1：预测三代全长转录组测序序列的蛋白编码序列：对三代全长转录组测序序列进行分析得到全长非嵌合序列，预测全长非嵌合序列可能的蛋白编码序列；

步骤1.2：将蛋白编码序列比对到目标基因组：将预测得到的所有蛋白编码序列比对到目标基因组，获得比对结果文件；

步骤1.3：提取外显子序列并重新比对回目标基因组：根据比对结果文件及基因组文件提取各外显子序列，将各外显子序列比对回目标基因组，得到各编码基因的初始结构信息。

本发明的有益效果为：