【GS文献】基因组选择技术在农业动物育种中的应用

中国农业大学等多家单位2017年合作发表在《遗传》杂志上的综述，笔记之。

作者中还有李宁院士，不胜唏嘘。

1.概述

• GS的两大难题：基因组分型的成本，基因组育种值（genomic estimted breeding value, GEBV）的准确性。
• 基于个体的基因组估计育种值GEBV比传统基于系谱的估计育种值（estimted breeding value, EBV）准确性更高。
• GS实施示意图：

基于单点SNP标记的GEBV估计方法

• 一类基于估计等位基因效应来计算GEBV：
  
  ①最小二乘法
  
  ②岭回归-最佳线性无偏预测（ridge regression best linear unbiased prediction, RR-BLUP）
  
  ③贝叶斯方法

• 另一类基于SNP构建基因组关系矩阵（GRM/G矩阵），再用线性混合模型估计GEBV：
  
  ①基因组最佳线性无偏预测法（genomic best linear unbiased prediction, GBLUP）
  
  ②一步法GBLUP（ssBLUP）

基于单倍型的GEBV 预测模型

• 随着个体数和标记密度提升，传统基于单点SNP 信息的统计模型无法利用所有个体表型数据或所有标记位点信息。
• 基于单倍型的GEBV 预测模型：通过将已知功能的基因组学信息，包括基因结构、甲基化区域、转录因子调控结合位点、选择信号候选区域等信息，以单倍型信息整合GEBV 预测模型，从而提高准确性（这辈子都可能看不到~）。
• 其他模型：增加显性、上位性、印记效应等。

2.基因型

高密度SNP芯片

代表：Illumina的Infinium 技术和Affymetrix的Axiom 技术

简化基因组测序：

• 如RAD-seq 、GBS、2b-RAD、GGRS、CroPS、RRL、MSG等，不同建库方法，但原理都是：通过合适的限制性内切酶对基因组进行酶切，然后连接二代测序平台所需的接头序列，再选择一定范围内的连接产物进行测序。
• GBS、RAD-seq 、2b-RAD最常用，区别在于是否在接头连接后进行片段选择，以及PCR后是否进行片段选择。

基因组重测序

• GS2.0（Hickey, 2013）：基于基因组重测序的GS，区别基于芯片的GS1.0。
• GS2.0另一种策略：首先对群体中遗传贡献较大的个体进行10~30x的中高覆盖度的重测序，结合单倍型估计技术构建群体主要的单倍型库；对于大多数个体采用SNP高密芯片或简化基因组测序技术进行分型，通过基因型填充获得全基因组的遗传变异，从而用GEBV的预测。
• 高密度SNP芯片的单价从几百到上千人民币不等，主要取决于定制量、位点密度等；简化基因组测序的单价通常芯片低，主要取决于所研究物种的基因组大小和所选择的内切酶。
• 也可以选20~30个代表性的个体进行重测序，获得这
  
  个群体具有一定代表性的SNP数据集，利用这些数据设计芯片，从而进行“廉价”的大规模基因分型。
  
  
• 挑战：分型的准确性和分析时效性。准确性可通过优化测序方案和基因型填充，时效性依靠硬件加速。

3.统计模型

贝叶斯

• BayesA 假设所有SNP都有效应，且所有SNP 效应的方差服从尺度逆卡方分布的正态分布，其中自由度和尺度参数这两个参数与遗传结构直接相关，能够确定遗传结构。BayesA 使用MCMC(马尔科夫链蒙特卡洛方法)方法构建Gibbs 抽样链，在模型中对标记效应进行求解。
• BayesB 只有一小部分标记位点有效应（有效应位点比例预先设定，一般为0.05，其效应方差服从的分布同BayesA），其他大部分染色体片段效应为0。BayesB使用混合分布作为标记效应方差的先验。
• BayesC 模型中，有/无效应位点比例是未知的，需要在模型中求解得到，其他同BayesB。
• BayesCπ 假设有效应位点的效应方差相同，而BayesC 假设其效应方差是不同的。
• Bayes LASSO 假定标记效应服从Laplace（拉普拉斯）分布，即等价于方差服从指数分布的正态分布。
• 以上是目前广泛应用的贝叶斯算法，结果稳定，准确性较高，此外有其他贝叶斯模型：BayesDπ，Bayes SSVS、fBayesB、wBSR 、BAL/IAL 、emBayesR 、EBL 、BayesRS、Bayes TA等。都是在基础模型上对先验假设进行变换和模型中的参数进行优化，以期寻找最适合群体的假设模型和参数。

GBLUP

• 通过构建基因组关系矩阵（G矩阵）替换基于系谱信息构建的分子血缘关系矩阵（NRM或A矩阵），进而使用最佳线性无偏预测（BLUP）方法估计GEBV。

• VanRaden构建G矩阵：先编码{1，0，1}，减去哈温平衡时均值，得到Z矩阵，如下计算G矩阵，根据MME求解育种值：
  
  

• 与贝叶斯方法相比，GBLUP不需要先利用参考群体估计标记效应，再计算GEBV；而是可以直接将有表型及无表型个体可以放在同一个模型中，同时估计出有表型和无表型个体的GEBV及其准确性。

ssGBLUP

• 基于系谱的BLUP（A矩阵）+ 基于SNP的GBLUP（G矩阵）= HBLUP（H矩阵）

• H矩阵：
  
  
  
  

• ssGBLUP能将无基因型个体的系谱及表型数据和有基因型的个体的基因组信息结合起来，大大提高基因组选择的准确性。

4.应用

奶牛：

应用最早最成熟，重点在缩短育种周期，难在扩大参考群和基因型鉴定。

猪：

重点在提高GEBV估计的准确性（因为世代间隔在现有育种体系中已经很短）。

难点是对传统育种中选择准确性低的性状预测，如繁殖性状、屠宰性状、肉质性状、抗病性状等。

GS 对猪的生长或胴体性状具有很好的预测准确性。

在母猪的繁殖性状（遗传力低）中应用受限于参考群体的大小。

鸡：

限制和难点是个体价值远小于基因分型成本（这个在植物中更加明显）。