生物信息学在新基因全长cDNA电子克隆中的应用

生物技术通报ＢＩＯＴＥＣＨＮＯＬＯＧＹＢＵＬＬＥＴＩＮ・技术与方法・２００７年第４期收稿日期：２００７－０１－３１作者简介：胡皝（１９８４－），男，在读硕士，研究方向：生物信息学随着人类基因组测序工作的基本完成，人类进入到了后基因组时代，基因组学的研究从结构基因组学过渡到了功能基因组学［１］，即从“是什么”过渡到“为什么”的研究。

然而，全基因组序列的解读，并不能使人类对编码基因这一层次有更明确的认识。

因此，ｃＤＮＡ的测序成为人们了解编码基因结构与功能的关键所在。

要理解新基因的结构和功能，仅有不完整的ｃＤＮＡ片段是不够的。

全长ｃＤＮＡ的获得是基因克隆的重要内容，也是目前基因组研究中的一个重要方面。

目前获取基因全长ｃＤＮＡ序列较常见的方法有：ｃＤＮＡ文库筛选法［２］、快速ｃＤＮＡ末端扩增法［３］和电子克隆法［４］等。

电子克隆法是近年来基于表达序列标签（ｅｘｐｒｅｓｓｅｄｓｅｑｕｅｎｃｅｔａｇ，ＥＳＴ）和基因组数据库发展起来的基因克隆新型技术，其利用生物信息学知识和计算机技术对ＥＳＴ或基因组数据库中进行同源性比较分析、整理拼接出新基因的编码序列，确认完整后根据序列设计引物进行ＲＴ－ＰＣＲ验证获得全长基因。

具有效率高、成本低、对实验条件要求低等特点［５］。

现以新基因全长ｃＤＮＡ电子克隆与分析的步骤为顺序，就生物信息学在其间的应用作一简单介绍。

１新基因全长ｃＤＮＡ电子克隆的方法及生物信息学在其中的应用１．１基于ＥＳＴ数据库的电子克隆ＥＳＴ是从ｃＤＮＡ克隆中随机挑选出来进行一次性测序的结果，一般长约２００ｂｐ￣５００ｂｐ，通常作为基因的标志。

近年来ＥＳＴ数据库容量扩增迅速，基于ＥＳＴ数据库由一个已知的基因利用生物信息学的方法进行功能基因的电子克隆已经成为目前最常用的基因克隆手段，许多新基因就是通过ＥＳＴ序列的拼接发现的［６，７］。

基于ＥＳＴ数据库的电子克隆大致步骤如下：第一步，选择其他物种尤其是亲缘关系较近的物种某生物信息学在新基因全长ｃＤＮＡ电子克隆中的应用胡皝萧浪涛（湖南农业大学生物科学技术学院，长沙４１０１２８）摘要：新基因全长ｃＤＮＡ序列的获得常常是生物学工作者面临的难题，电子克隆是利用生物信息学手段得到新基因全长ｃＤＮＡ序列的新方法。

介绍了电子克隆的方法及其生物信息学在其间的具体应用，并概述了一些生物信息学在序列分析中的应用。

关键词：生物信息学电子克隆全长ｃＤＮＡＡｐｐｌｉｃａｔｉｏｎｏｆＢｉｏｉｎｆｏｒｍａｔｉｃｓｉｎＦｕｌｌ－ｌｅｎｇｔｈｃＤＮＡＳｅｑｕｅｎｃｅｉｎＳｉｌｉｃｏｎＣｌｏｎｉｎｇｏｆＮｏｖｅｌＧｅｎｅｓＨｕＨｕａｎｇＸｉａｏＬａｎｇｔａｏ（ＣｏｌｌｅｇｅｏｆＢｉｏｓｃｉｅｎｃｅａｎｄＢｉｏｔｅｃｈｎｏｌｏｇｙ，ＨｕｎａｎＡｇｒｉｃｕｌｔｕｒａｌＵｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ４１０１２８）Ａｂｓｔｒａｃｔ：Ｔｏｏｂｔａｉｎｔｈｅｆｕｌｌ－ｌｅｎｇｔｈｃＤＮＡｓｅｑｕｅｎｃｅｏｆａｎｏｖｅｌｇｅｎｅｉｓａｐｒｏｂｌｅｍｆｏｒｒｅｓｅａｒｃｈｅｒｓ．Ｉｎｓｉｌｉｃｏｃｌｏｎｉｎｇｗａｓａｎｅｗｗａｙｏｆｏｂｔａｉｎｔｈｅｆｕｌｌ－ｌｅｎｇｔｈｃＤＮＡｓｅｑｕｅｎｃｅｂｙｂｉｏｉｎｆｏｍａｔｉｃｓ．Ｔｈｉｓａｒｔｉｃｌｅｉｎｔｒｏｄｕｃｅｓｔｈｅｗａｙｏｆｉｎｓｉｌｉｃｏｃｌｏｎｉｎｇａｎｄｔｈｅｂｉｏｉｎｆｏｒｍａｔｉｃｓａｐｐｌｉｃａｔｉｏｎｉｎｉｔ．Ｉｔｗａｓａｌｓｏｓｕｍｍａｒｉｚａｔｉｏｎｏｆｔｈｅｂｉｏｉｎｆｏｒｍａｔｉｃｓａｐｐｌｉｃａｔｉｏｎｉｎｔｈｅｓｅｑｕｅｎｃｅａｎａｌｙｓｅ．Ｋｅｙｗｏｒｄｓ：ＢｉｏｉｎｆｏｒｍａｔｉｃｓＩｎｓｉｌｉｃｏｎｃｌｏｎｉｎｇＦｕｌｌ－ｌｅｎｇｔｈｃＤＮＡ生物技术通报ＢｉｏｔｅｃｈｎｏｌｏｇｙＢｕｌｌｅｔｉｎ２００７年第４期基因全长ｃＤＮＡ序列或ＥＳＴ序列为查询探针或者以该物种某基因ＥＳＴ为查询探针，搜索ＥＳＴ数据库进行Ｂｌａｓｔ比对，得到许多ＥＳＴ序列，从中寻找感兴趣的ＥＳＴ（标准的选择与可预计的同源基因的同源程度有关。

通常为：同源长度≥１００ｂｐ，同源性５０％以上，８５％以下）。

第二步，把感兴趣的ＥＳＴ基于ＧｅｎＢａｎｋ中的非冗余数据库进行Ｂｌａｓｔ分析，判断其是否是已知基因的一部分，筛选出新颖的ＥＳＴ。

第三步，将筛选出的ＥＳＴ在该物种的ＥＳＴ数据库中进行搜索，找到部分重叠的ＥＳＴ进行拼接，经严格聚类分析，尽量避免含有旁系同源基因，拼接后产生的序列重叠群），相当于实验中的一部分ｃＤＮＡ步移工作。

第四步，以新获得的重叠群为新的查询探针，继续搜索ＥＳＴ数据库，直到没有新的ＥＳＴ可供拼接为止。

将拼接得到的序列对非冗余数据库进行搜索，以证明这是一个全新的序列。

这种策略也存在一定的局限性，许多拷贝数较低的基因很难涵盖在ＥＳＴ数据库中，这些基因只能通过分析基因组序列才能被发现。

ＥＳＴ序列的拼接是电子克隆中非常重要的环节，用于ＥＳＴ序列的拼接的软件有很多，表１列出了一些比较常用的拼接软件［８］，使用者可按具体情况选择不同的拼接软件以得到最好的结果。

另外，还可以将序列提交到ＮＣＢＩ的ＵｎｉＧｅｎｅ［９］数据库上。

ＮＣＢＩ的Ｕｎｉｇｅｎｅ系统是ＧｅｎＢａｎｋ中的序列另外分离出来形成一个非冗余的基因簇。

数据库中除包含已确定的基因以外，还包括数以万计的表达序列标签，每个簇包含惟一的非冗余的基因序列，表达的组织类型和基因图谱位点。

现在数据库中已经包括大量模式或重要生物的ＥＳＴ序列，其中人类、老鼠和水稻的序列最多。

通过Ｕｎｉｇｅｎｅ系统可以很方便地进行序列的拼接得到新基因。

表１常用的拼接软件及网址１．２基于基因组数据库的电子克隆人类基因组及其他许多模式、重要物种基因组测序工作的完成，基于基因组序列的新基因预测软件的开发为我们利用生物信息学的方法克隆新基因带来了新的策略。

近年来，许多新基因［１０］就是通过分析基因组序列发现的。

基于基因组数据库的电子克隆大致步骤如下：第一步，选择其他物种尤其是亲缘关系较近的物种某基因全长ｃＤＮＡ序列或ＥＳＴ序列为查询探针或者以该物种某基因ＥＳＴ为查询探针，基于ＧｅｎＢａｎｋ中的非冗余数据库ｎｒ进行Ｂｌａｓｔ分析，从结果中筛选出同源性较高、含外显子的该物种基因组重叠群或ＢＡＣ克隆，并通过超级链接得到其所在的基因组序列，同时根据比对的结果对基因组序列可能造成的移码测序错误进行修正。

第二步，将这些外序列根据内含子和外显子的剪接特征“ＧＵ．．．ＡＧ”，通过人工拼接，或者通过基因预测软件预测，得到可能的新基因序列。

第三步，把可能的新基因序列基于非冗余数据库做ＢＬＡＳＴ分析，检验其新颖性。

第四步，把筛选后的新基因序列提交到ｄｂＥＳＴ数据库做ＢＬＡＳＴ分析并延伸，进一步确认其真实度。

用于基因预测的软件有很多，国际上使用最广泛的有ＧｅｎＳｃａｎ、ＧｅｎｅＦｉｎｄｅｒ和ＦＧＥＮＥＳＨ等，此外还有一些针对某一物种的专用基因预测软件，如水稻基因预测软件ＲｉｃｅＨＭＭ、玉米基因预测软件ＳｐｌｉｃｅＰｒｅｄｉｃｔｏｒ和人类基因预测软件ＨｕｍＧｅｎｅ［１１］等。

由于各种预测软件依据的原理不用，得出的结果也有所差异，使用者可同时利用几种预测软件的分析结果并结合自己的经验综合考虑，确定最有可能的基因结构。

实际操作中常常将上述两种方法结合起来使用，使最后所得的序列更加准确。

首先进行ＥＳＴ序列的拼接，无法拼接后再进行基因组比较和外显子预测，以判断ＥＳＴ拼接的完整性和正确性。

１．３全长ｃＤＮＡ的判断运用以上方法得到的ｃＤＮＡ序列还不能确定９４２００７年第４期其为全长的ｃＤＮＡ序列，需要进行判断。

直接从序列上可以从如下几个方面［１２］进行判断。

５′端：（１）有同源全长基因的比较，通过与其它生物已有的对应基因末端进行Ｂｌａｓｔ来判断。

（２）无同源基因的新基因，Ｉ．判断编码框架是否完整，ａ．在开放阅读框架的第一个ＡＴＧ上游有同框架的终止密码，需要注意的是有时真正的翻译起始密码子并非是出现在ｍＲＮＡ中的第一个ＡＴＧ，在有的真核细胞中，在起始密码子ＡＴＧ的上游非编码区会有可能出现一到几个ＡＴＧ，这称为非编码的５′ＡＴＧ。

以这种５′ＡＴＧ并不是真是的起始密码子，以其开始的开放阅读框常常很快遇到终止密码子。

ｂ．无终止密码的则考虑有保守的Ｋｏｚａｋ［１３］序列；Ⅱ．判断是否自转录起始点，有资料表明，在５′帽结构后一般都有一段富含嘧啶的区域，另外如果ｃＤＮＡ５′序列与基因组序列中经Ｓ１酶切保护的部分相同，则可以确定得到的ｃＤＮＡ是全长的。

３′端：（１）有同源全长基因的比较，方法同５′端；（２）编码框架的下游有终止密码；（３）有一个以上的ｐｏｌｙＡ加尾信号；（４）无明显加尾信号的则也有ｐｏｌｙＡ尾。

同源全长基因的比较可以用Ｂｌａｓｔ比对或多重序列比对软件来实现，首先搜索到其他物种（以相似性比较高的物种为宜）该基因全长ｃＤＮＡ序列，再将这些序列做Ｂｌａｓｔ或多重比对，ＣｌｕｓｔａｌＷ［１４］是目前使用最广泛的多重比对软件，使用者可以将序列提交到ｈｔｔｐ：／／ｗｗｗ．ｅｂｉ．ａｃ．ｕｋ／ｃｌｕｓｔａｌｗ／进行在线分析，或者下载该软件到本地进行比对。

确定得到的ｃＤＮＡ序列为全长的ｃＤＮＡ序列还只是在计算机上的“虚拟克隆”，最终还必须通过ＲＴ－ＲＣＲ、序列测定和Ｎｏｒｔｈｅｒｎ杂交等方法进行实验验证，以保证序列的准确性。

２生物信息学在新基因全长ｃＤＮＡ电子克隆和序列分析中的应用生物信息学对新基因序列的分析也有很大的指导作用，可以对该序列及其所编码的蛋白质进行结构和功能的详细分析，能为使用者节约大量的时间。

表２列出了一些常用的分析软件及网址。

表２常用的分析软件及网址３讨论与展望综上主要介绍了基于Ｉｎｔｅｒｎｅｔ网上生物信息资源新基因全长ｃＤＮＡ的电子克隆策略，对于那些需要经常进行序列分析，或者分析规模较大的实验室，我们可以构建本地的生物信息学数据分析平台，把一些重复性的、可程序化的过程直接交由计算机完成，用户的主要精力就可转移到对分析结果进行后续分析及实验设计上，能节约大量的人力和物力，具体的构建方法可以查阅有关此方面的一系列论著［１５，１６］。

生物信息技术的应用使新基因全长ｃＤＮＡ克隆和分析的方法不断更新，朝着快速、经济、准确的方向发展，但鉴于生物大分子结构和功能的复杂性，许多分析软件的输出结果存在较大的偏差，因此利用生物信息学进行的＂虚拟＂克隆的结果尚需回到实验室进行验证。

但是，这种分析方法为实验研究提供了重要的线索，对随后的研究起到了“事半功倍”的作用，避免走弯路，极大地提高了工作效率。

可以相信随着基因组序列信息的日益丰富，计算方法和数据库的不断完善，生物信息学将在基因全长ｃＤＮＡ克隆和分析中扮演更加重要的角色。

参考文献１ＳｔｕａｒｔＫ，Ｋｉｍ．ＣｕｒｒｅｎｔＢｉｏｌｏｇｙ，２００１，１１：８５￣８７．２ＣａｉＹ，ＧａｏＹ，ｅｔａｌ．ＭｏｌＲｅｐｒｏｄＤｅｖ，２００２，６１（１）：１２６￣１３４．３ＳｃｈａｅｆｅｒＢＣ．ＡｎａｌｙｔｉｃａｌＢｉｏｃｈｅｍｉｓｔｒｙ，１９９５，２２７：２５５￣２７３．４ＳａｎｄｒｏＢａｎｆｉ，ＡｌｅｓｓａｎｄｒｏＧｕｆｆａｎｔｉ．ＴｒｅｎｄｓｉｎＧｅｎｅｔｉｃｓ，１９９８，１４：８０￣８１．胡皝等：生物信息学在新基因全长ｃＤＮＡ电子克隆中的应用９５生物技术通报ＢｉｏｔｅｃｈｎｏｌｏｇｙＢｕｌｌｅｔｉｎ２００７年第４期人类基因功能研究和疾病基因治疗开辟了一个革命性的新领域［２３］。