|
|
| (66 intermediate revisions by 5 users not shown) |
| Line 1: |
Line 1: |
| =EPD=
| | #REDIRECT [[:Category:Genetic Evaluation]] |
| ==Utility (compared to actual/adjusted phenotypes, ratios, disjoined marker scores, etc.) (Suggested writer: Megan Rolf)==
| | Predicting genetic merit for breeding animals is one of the oldest practices that mankind has used to improve food and fiber production. Identifying animals for [[Selection and Mating | selection and mating]] has evolved from visual appraisal to sophisticated analytical models for predicting [[Glossary#A | additive genetic]] merit of animals. Additive genetic merit is the effect of genes that are passed from parent to offspring that can be used to make genetic progress through selection. |
| ==Basic Models==
| |
| ===BLUP===
| |
| ===ssGBLUP (Suggested writer: Daniela Lourenco)===
| |
|
| |
|
| Single-step genomic BLUP (ssGBLUP) <ref> Legarra, A., I. Aguilar, and I. Misztal. 2009. A relationship matrix including full pedigree and genomic information. J. Dairy Sci. 92:4656-4663. </ref><ref> Aguilar, I., I. Misztal, D. L. Johnson, A. Legarra, S. Tsuruta, and T. J. Lawlor. 2010. Hot topic: A unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. Journal of Dairy Science 93: 743-752. </ref> is a method developed to enable the inclusion of marker genotypes into the well-known BLUP machinery. The idea of ssGBLUP came from the fact that only a small portion of the animals in the pedigree is genotyped. In this way, one approach to account for all animals (i.e., genotyped and non-genotyped) in the evaluation would be to combine pedigree and genomic relationships and use this as the covariance structure in the BLUP mixed model equations. Thus, ssGBLUP uses marker information to construct genomic relationships.
| | In North America, the standard for identifying genetic merit of breeding animals is [[Expected Progeny Difference | expected progeny differences (EPDs)]]. |
| | With very few ''ad hoc'' exceptions, EPDs are produced for North American beef cattle using models based on [[Best Linear Unbiased Prediction]]. Consequently, [[BIF recommends the use of EPD]] when available. |
|
| |
|
| Legarra et al. (2009)<ref> Legarra, A., I. Aguilar, and I. Misztal. 2009. A relationship matrix including full pedigree and genomic information. J. Dairy Sci. 92:4656-4663. </ref> stated that genomic evaluations would be simpler if genomic relationships were available for all animals in the model. Then, their idea was to look at the pedigree relationship as a priori relationship and at the genomic relationship as the observed relationship. Based on that, they showed the genomic information could be extended (i.e., imputed) to non-genotyped animals. This means that in ssGBLUP pedigree relationships for non-genotyped animals are enhanced by the genomic information of their relatives. The relationship matrix that combines information for genotyped and non-genotyped animals is represented by ''H'':
| | While not all [[Economically Relevant Traits | economically relevant traits]] in all situations and in all North American breed registries have EPDs available, the number of [[Traits | traits and trait components]] that have EPDs has increased dramatically. |
| | Nearly all the major North American beef cattle breed organizations have migrated to weekly genetic evaluations, eliminating the need for [[Expected Progeny Difference#Interim EPDs | interim EPDs]]. |
|
| |
|
| <center>
| | Most of the improvements in the technologies used in genetic evaluation have been motivated by an opportunity to increase [[Accuracy | accuracy of prediction]] and reduce [[Prediction Bias | bias]]. For example, the advent of [[Genotyping | genomic information]] to enhance the [[Accuracy | accuracy]] of prediction has resulted in EPDs for most traits being produced using either [[Single-step Genomic BLUP]] or [[Single-step Hybrid Marker Effects Models]]. The BIF has developed an extensive set of recommendations for the inclusion of [[Genomic Evaluation Guidelines | genomic data in genetic evaluations]]. |
| <math>
| |
| H=
| |
| \begin{bmatrix}
| |
| A_{11}+A_{12}A_{22}^{-1}(G-A_{22})A_{22}^{-1}A_{21} & A_{12}A_{22}^{-1}G \\
| |
| GA_{22}^{-1}A_{21} & G
| |
| \end{bmatrix}
| |
| </math>
| |
| </center>
| |
|
| |
|
| Where the subscripts 1 and 2 refer to non-genotyped and genotyped animals, respectively. ''A'' is the pedigree relationship matrix and ''G'' is the genomic relationship matrix computed based on markers. If ''M'' is a matrix of marker genotypes centered for allele frequency (''p'') and has the dimension of number of animals by number of SNP (n), ''G'' is computed as <ref> VanRaden, P. M. 2008. Efficient methods to compute genomic predictions. J. Dairy Sci. 91:4414-4423. </ref>:
| | In commercial cattle production, EPDs for [[Economically Relevant Traits | economically relevant traits]] should be combined with appropriate selection tools such as [[Selection Index | selection indices]] to make optimal genetic progress toward achieving [[Breeding Objectives | breeding objectives]]. It must be remembered that EPDs are just tools to make selection decisions to make genetic progress and manage certain genetic risks. |
| <center>
| |
| <math>
| |
| G=\frac{MM^'}{2\sum_{i=1}^n p_i(1-p_i)}
| |
| </math>
| |
| </center>
| |
|
| |
|
| | In some special situations in seedstock production breeders may need to make selection decisions using EPDs that are not [[Economically Relevant Traits | economically relevant traits]] in commercial settings in order to enhance the marketability of their breed or breeding animals. For example, if a breed has a perceived defect that is limiting that breed organizations' members from expanding their market for selling germplasm, then selection to improve that characteristic should be included in the seedstock breeder's [[Breeding Objectives | breeding objectives]]. |
|
| |
|
| Although ''H'' is very complicated, ''H''<sup>''-1''</sup> is quite simple <ref> Aguilar, I., I. Misztal, D. L. Johnson, A. Legarra, S. Tsuruta, and T. J. Lawlor. 2010. Hot topic: A unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. Journal of Dairy Science 93: 743-752. </ref>:
| | Critical to genetic evaluation is having high-quality estimates of [[Variance Components | variance components]]. Knowing the heritabilities and correlations of the traits and performing [[Multiple Trait Evaluation | Multiple-Trait Evaluation]] enhances the accuracy of prediction and reduces [[Prediction Bias | bias]] from effects such as incomplete reporting. Equally critical is understanding the [[Connectedness | connectedness]] of the data in a particular data set. Disconnected data can lead to invalid comparisons. |
| <center>
| |
| <math>
| |
| H^{-1}=A^{-1}+
| |
| \begin{bmatrix}
| |
| 0 & 0 \\
| |
| 0 & G^{-1}- A_{22}^{-1}
| |
| \end{bmatrix}
| |
| </math>
| |
| </center>
| |
|
| |
| If we replace A^{-1} by H^{-1} in the BLUP mixed model equations, we have ssGBLUP<ref> Aguilar, I., I. Misztal, D. L. Johnson, A. Legarra, S. Tsuruta, and T. J. Lawlor. 2010. Hot topic: A unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. Journal of Dairy Science 93: 743-752. </ref>:
| |
| <center>
| |
| <math>
| |
| \begin{bmatrix}
| |
| {X^'X} & {X^'Z} \\
| |
| {Z^'X} & {Z^'Z}+ H^{-1}{\lambda}
| |
| \end{bmatrix}
| |
| \begin{bmatrix}
| |
| {b} \\
| |
| {u}
| |
| \end{bmatrix}
| |
| =
| |
| \begin{bmatrix}
| |
| {X^'y} \\
| |
| {Z^'y}
| |
| \end{bmatrix}
| |
| </math>
| |
| </center>
| |
| Where ''b'' and ''u'' are vectors of fixed effects and breeding values, respectively; ''X'' and ''Z'' are incidence matrices for the effects in ''b'' and ''u''; ''y'' is a vector of phenotypes, and λ is the ratio of residual to additive genetic variance.
| |
| | |
| As a combined relationship is used in ssGBLUP, the output for each animal is automatically a genomic EBV, and the mixed model equations above can be simplistically represented as:
| |
| | |
| <center>
| |
| [[File:figure1_ssGBLUP.jpg|250px]] | |
| </center>
| |
| | |
| The genomic EPD is then calculated as:
| |
| <center>
| |
| <math>
| |
| {genomic\,EPD}=\frac{genomic\,EBV}{2}
| |
| </math>
| |
| </center>
| |
| When the subject is genetic evaluation, one of the most common questions is “What is the main difference among ssGBLUP, BLUP, and genomic BLUP (GBLUP)?” In a nutshell, ssGBLUP uses phenotypes, pedigree, and genotypes for both genotyped and non-genotyped animals, whereas BLUP uses phenotypes and pedigree for all animals and GBLUP uses phenotypes and genotypes only for genotyped animals.
| |
| | |
| In the US, ssGBLUP has been used for genomic evaluation of beef and dairy cattle, pigs, chickens, and fish. Regarding to beef cattle, Angus Genetics Inc. runs ssGBLUP evaluations for American Angus and Charolais, Canadian Angus, Red Angus, and Charolais, and Maine Anjou. Moreover, Livestock Genetic Services (A Neogen Company) runs ssGBLUP evaluations for Santa Gertrudis. For more information about ssGBLUP for beef cattle evaluation check Lourenco et al. (2015)<ref> Lourenco, D. A. L., S. Tsuruta, B. O. Fragomeni, Y. Masuda, I. Aguilar, A. Legarra, J. K. Bertrand, T. Amen, L. Wang, D. W. Moser, and I. Misztal. 2015. Genetic evaluation using single-step genomic best linear unbiased predictor in American Angus. Journal of Animal Science 93: 2653-2662.</ref>, Lourenco et al. (2017)<ref> Lourenco, D.A.L., J.K. Bertrand, H.L. Bradford, S. Miller, and I. Misztal. 2017. The promise of genomics for beef improvement. BIF Meeting (http://www.bifconference.com/bif2017/proceedings/01-lourenco.pdf) </ref>, and Misztal & Lourenco (2018) <ref> Misztal, I. and D. Lourenco. 2018. Current research in unweighted and weighted ssGBLUP. In Proc. Beef Improvement Federation 11th genetic prediction workshop 11:6-13. </ref>.
| |
| | |
| <em>ssGBLUP for large genotyped populations</em> <p></p>
| |
| Running ssGBLUP evaluations for large genotyped populations can be a huge computational challenge. This is because the construction of ''H''<sup>''-1''</sup> requires the construction and inversion of ''G''. Matrix inversion has a cubic computational cost, requiring a large amount of memory. As an example, inverting ''G'' for 100,000 animals requires about 300Gb of memory and takes over 2 hours.
| |
| | |
| The algorithm for proven and young (APY) was proposed by Misztal et al. (2014) <ref> Misztal, I., A. Legarra, and I. Aguilar. 2014. Using recursion to compute the inverse of the genomic relationship matrix. J. Dairy Sci. 97: 3943–3952. </ref> to overcome this computing limitation of ssGBLUP, and was based on Henderson’s algorithm to construct ''A''<sup>''-1''</sup> <ref> Henderson, C.R. 1976. A simple method for computing the inverse of a numerator relationship matrix used in the prediction of breeding values. Biometrics 32:69-83.</ref>. In APY, ''G''<sup>''-1''</sup> is constructed directly, avoiding the matrix inversion step. In this algorithm, genotyped animals are split into two groups: core and non-core. Breeding values of non-core animals are then calculated as functions of breeding values of core animals and the genomic relationships between core and non-core. If the number of genotyped animals surpasses 100,000, using APY ''G''<sup>''-1''</sup> in ssGBLUP is highly recommended. A simple graphic representation of APY ''G''<sup>''-1''</sup> is:
| |
| <center>
| |
| [[File:figure2_APY.jpg | 250px]] | |
| </center>
| |
| | |
| | |
| <em>Marker effects in ssGBLUP</em><p></p>
| |
| Although the marker information in ssGBLUP is used to construct genomic relationships, it is possible to calculate SNP effects once we obtain genomic EBV (Wang et al., 2012<ref> Wang, H., I. Misztal, I. Aguilar, A. Legarra, and W. M. Muir. 2012. Genome-wide association mapping including phenotypes from relatives without genotypes. Genet. Res. 94(2):73-83.</ref>, Lourenco et al., 2015<ref> Lourenco, D. A. L., S. Tsuruta, B. O. Fragomeni, Y. Masuda, I. Aguilar, A. Legarra, J. K. Bertrand, T. Amen, L. Wang, D. W. Moser, and I. Misztal. 2015. Genetic evaluation using single-step genomic best linear unbiased predictor in American Angus. Journal of Animal Science 93: 2653-2662.</ref>). Marker effects can be then used to calculate predictions based only on marker genotypes for young genotyped animals that are not yet or will never make into an official evaluation.
| |
| | |
| ===[[Single-step Hybrid Marker Effects Models]] (Suggested writer: Bruce Golden)===
| |
| | |
| ==Interim Calculations==
| |
| =Bias=
| |
| ==(in)complete reporting / contemporary groups / preferential treatment (Suggested writer: Bob Weaber==
| |
| =[[Accuracy]] (Suggested writer: Matt Spangler)=
| |
| ==meaning of accuracy==
| |
| ==what impacts accuracy==
| |
| ==different definitions of accuracy (true, BIF, reliability)==
| |
| | |
| =Variance components (Suggested writer: Steve Kachman)=
| |
| ==Impact on EPD, accuracy, genetic gain (Suggested writer: Steve Kachman)==
| |
| ==[[Heterogeneous variance]]==
| |
| | |
| =Connectivity (Suggested writer: Ron Lewis)=
| |
| ==Measures of (Suggested writer: Ron Lewis)==
| |
| ==Impact on GE== (Suggested Writer: Ron Lewis)
| |
| =Current GE=
| |
| ==How each breed (organization) is modeling each trait (Suggested writers: Steve Miller, Lauren Hyde, AHA)==
| |
Predicting genetic merit for breeding animals is one of the oldest practices that mankind has used to improve food and fiber production. Identifying animals for selection and mating has evolved from visual appraisal to sophisticated analytical models for predicting additive genetic merit of animals. Additive genetic merit is the effect of genes that are passed from parent to offspring that can be used to make genetic progress through selection.
In North America, the standard for identifying genetic merit of breeding animals is expected progeny differences (EPDs).
With very few ad hoc exceptions, EPDs are produced for North American beef cattle using models based on Best Linear Unbiased Prediction. Consequently, BIF recommends the use of EPD when available.
While not all economically relevant traits in all situations and in all North American breed registries have EPDs available, the number of traits and trait components that have EPDs has increased dramatically.
Nearly all the major North American beef cattle breed organizations have migrated to weekly genetic evaluations, eliminating the need for interim EPDs.
Most of the improvements in the technologies used in genetic evaluation have been motivated by an opportunity to increase accuracy of prediction and reduce bias. For example, the advent of genomic information to enhance the accuracy of prediction has resulted in EPDs for most traits being produced using either Single-step Genomic BLUP or Single-step Hybrid Marker Effects Models. The BIF has developed an extensive set of recommendations for the inclusion of genomic data in genetic evaluations.
In commercial cattle production, EPDs for economically relevant traits should be combined with appropriate selection tools such as selection indices to make optimal genetic progress toward achieving breeding objectives. It must be remembered that EPDs are just tools to make selection decisions to make genetic progress and manage certain genetic risks.
In some special situations in seedstock production breeders may need to make selection decisions using EPDs that are not economically relevant traits in commercial settings in order to enhance the marketability of their breed or breeding animals. For example, if a breed has a perceived defect that is limiting that breed organizations' members from expanding their market for selling germplasm, then selection to improve that characteristic should be included in the seedstock breeder's breeding objectives.
Critical to genetic evaluation is having high-quality estimates of variance components. Knowing the heritabilities and correlations of the traits and performing Multiple-Trait Evaluation enhances the accuracy of prediction and reduces bias from effects such as incomplete reporting. Equally critical is understanding the connectedness of the data in a particular data set. Disconnected data can lead to invalid comparisons.
