软件介绍: BLUPF90的无敌和寂寞

前同桌说, 要带粉丝了, 一个换一个, 我就把题目改了一下, 原来的题目是《遗传育种软件三剑客之一：BLUPF90》，但我还是喜欢现在这个题目《BLUPF90的无敌和寂寞》。

BLUPF90号称遗传评估领域的三剑客之一(另外两个分别是ASREML和DMU), 其应该是最先开发了支持基因组数据的模块, 并且里面的基因组数据过滤筛选非常友好。除了上述优势外，BLUPF90还具有以下功能：

• SSGWAS：之前的GWAS都是基于G矩阵或SNP数据, 现在可以同时利用系谱和SNP信息；
• PCA分析：之前的PCA都是基于G矩阵或者SNP数据, 现在可以利用H矩阵来进行群体划分；
• 亲子鉴定：BLUPF90的seekparents模块在进行系谱矫正、亲子鉴定分析时, 速度非常快!
• 重复样本和异常样本检查：BLUPF90的模块在检查样本可能存在的重复或者异常现象时使用方便；
• H逆矩阵的输出与再利用：BLUPF90输出的H逆矩阵，在DMU和ASREML中同样适用。

1. ABLUP VS SSGBLUP

传统ABLUP与SSGBLUP的区别在于，原来的A逆矩阵变为了H逆矩阵。

1.1 传统的动物模型计算BLUP值

根据系谱计算A矩阵，然后使用Henderson方程组计算BLUP（EBV值）。

1.2 SSGBLUP计算BLUP值

2. BLUPF90进行SSGBLUP分析

• RENUMF90对数据进行预处理
• BLUPF90进行SSGBLUP的计算（已知方差组分时）
• AIREMLF90进行SSGBLUP方差组分估算（未知方差组分时）
• PREGSF90进行H逆矩阵的计算

2.1 renumf90预处理数据

2.2 blupf90计算H矩阵以及育种值

2.3 pregsf90计算H矩阵

主要作用:

• 根据基因组数据构建G矩阵
• 根据系谱构建A矩阵
• 根据公式构建H逆矩阵

3. 构建H矩阵的参数设置

4. 基因组数据的筛选

5. 亲子鉴定的作用

6. G矩阵结果数据挖掘

6.1 检测异常个体

G矩阵中, 某些个体对角线有较高的值。这说明有可能这个个体不是群体内的个体, 可能来源于其它群体或者家系, 或者call rate值较低。

6.2 检测重复样本

如果某两个个体的亲缘关系大于0.9, 则表明这两个个体可能是重复样本。

6.3 G矩阵和A22矩阵的关系

G矩阵和A22矩阵是相同个体构建的G矩阵和A矩阵，因此二者应该具有很高的相似性。如果对角线和非对角线相似度较低, 这说明出现了一些问题，需要引起重视。可能是测序个体ID错误, 也可能是数据量较少导致等等。

6.4 根据基因组数据进行PCA分析

7. 构建A22矩阵时的高效方法

可以看出, 构建A22矩阵时, 使用57,000系谱数据, 6,500测序个体, Tabular用了311s, 内存占用12G, 而Colleau method用了45s, 内存占用322Mb。因此使用Colleau方法更合适。

8. BLUPF90构建H逆矩阵输出

如果想要使用DMU, ASREML或者WOMBAT利用BLUPF90构建好的H逆矩阵，需要输出Original ID的形式。然后转化为DMU和ASREML的格式即可。