[发明专利]基于生物信息的选育良种系统及其算法

说明书

技术领域

本发明涉及生物技术领域，具体涉及了一种基于生物信息的选育良种系统及其算法。

背景技术

在农业生物(作物、畜禽和水产生物)的产业链中,良种选育占有至关重要的地位。传统的育种方法大多基于遗传统计原理并结合生产实际,利用候选个体本身和(或)亲属的性状记录估计得到其育种值，并利用育种值为指标通过进行多世代的人工选育以实现品种培育和性状改良。然而这种方法往往存在着育种周期长,人力物力消耗大等缺点。如法国的“大约克”种猪就花了三十多年时间才选育成功。在乳牛育种领域通过传统后裔测定方式验证一头候选种公牛需要5-6年时间，大约投入5万美元。

基于候选个体的遗传性状由其所携带的基因来决定的原理，通过利用候选个体的基因信息在基因层面来直接评估而非仅仅利用系谱信息和记录性状的表型数据来考量候选个体的遗传性状的优劣，将能够大大缩短育种时间及成本投入。利用候选个体的基因信息的育种方法主要包括两种。第一种方法称为标记辅助选择(Marker Assisted Selection,MAS)方法，其试图通过部分基因组信息，即确定控制表型性状的基因信息来实现对目标性状的选择。该方法被有效利用的前提条件是与控制表型性状的基因相关的数量性状位点(QTL)能够确定并被精确定位。然而在实际育种应用中，决定某一经济性状的往往有很多基因，而每个基因的效应都不明显，从而造成QTL检测和定位的困难。即便存在有对经济性状贡献较大的主效基因，在构成表型性状的所有遗传变异中应用于MAS的遗传标记只能捕获主效基因所带来的那部分变异，而小效应累加起来所带来的变异却被忽视了。所以MAS方法只能有限地应用在经济性状的控制基因明确且控制基因的效应显著的性状选择中。

为了捕获构成表型的所有遗传变异，需要在候选个体整个基因组水平上检测影响目标性状的所有QTL来进行选择利用，即在候选个体全基因组范围内进行标记辅助选择，这就是第二种方法：本发明所涉及的全基因组选择方法。该方法利用个体整个基因组高密度遗传标记图谱以保证影响目标性状的每一个QTL都与高密度全基因组标记图谱中的至少一个遗传标记处于连锁不平衡(Linkage Disequilibrium，LD)状态。通过这种连锁不平衡信息来实现对所有影响目标性状的QTL的捕捉，再通过LDMIP方法来实现良种选育。该方法比目前全基因组选择的另一常用方法GBLUP有更高的准确率。

全基因组遗传标记还可以提供系谱中个体间更加精确的亲缘遗传关系，即连锁分析(Linkage Analysis，LA)信息。利用个体及其父母的全基因组信息可以计算个体等位基因从父母处遗传的遗传几率，可以用来示踪等位基因在系谱中一代一代的传递。这种连锁分析信息对于例如检测和确定与疾病相关的功能基因位置方面非常重要:利用遗传标记可以示踪附近对疾病有效应的基因。这样如果标记附近有功能基因，并且样本足够大，携带不同标记等位基因的个体的表现就会显著不同。在品种选育中，人们除了追求高产优质等特性之外，能够有效抵御疾病的抗病能力也往往是重要的选育特性。通过连锁分析信息便可以检测在一个家系中等位基因与疾病的传递是否相关，有助于选育抗病能力强的个体。

发明内容

为解决上述技术问题，我们提出了一种基于生物信息的选育良种系统及其算法，其目的：准确预测以实现高效率特异性育种。

为达到上述目的，本发明的技术方案如下：

一种基于生物信息的选育良种系统，具体实施步骤如下：

一、采用全基因组选择方法,首先获取参考群体和候选群体个体的基因组数据，其中，采用SNP遗传标记来获取基因型数据；

二、获取系谱文件并进行预处理，向上追溯父母系谱信息至少5代及以上，系谱文件需要包含个体编号、父亲编号和母亲编号至少三段信息，在系谱文件中筛选具有基因型数据的基因分型个体，向上追溯父母系谱信息至少5代及以上以获得GA系谱；通过基因分型个体的基因型数据信息对系谱文件中父母及后代亲缘关系进行亲子鉴定以核查GA系谱文件，当根据基因型数据获得的亲子鉴定结果与原始GA系谱文件所显示的亲子关系不一致时，应以基因型数据获得的亲子鉴定结果为基准重新编排原始GA系谱文件；

三、在完成基因型数据的预处理并得到了GA系谱文件后，通过LDMIP方法来完成，具体的LDMIP方法，如下:

A、GA系谱中未基因分型的祖先个体，推知其基因型数据；

B、根据基因型数据和GA系谱信息，计算个体等位基因从父母处遗传的遗传几率；