1
|
Mota LFM, Arikawa LM, Santos SWB, Fernandes Júnior GA, Alves AAC, Rosa GJM, Mercadante MEZ, Cyrillo JNSG, Carvalheiro R, Albuquerque LG. Benchmarking machine learning and parametric methods for genomic prediction of feed efficiency-related traits in Nellore cattle. Sci Rep 2024; 14:6404. [PMID: 38493207 PMCID: PMC10944497 DOI: 10.1038/s41598-024-57234-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Accepted: 03/15/2024] [Indexed: 03/18/2024] Open
Abstract
Genomic selection (GS) offers a promising opportunity for selecting more efficient animals to use consumed energy for maintenance and growth functions, impacting profitability and environmental sustainability. Here, we compared the prediction accuracy of multi-layer neural network (MLNN) and support vector regression (SVR) against single-trait (STGBLUP), multi-trait genomic best linear unbiased prediction (MTGBLUP), and Bayesian regression (BayesA, BayesB, BayesC, BRR, and BLasso) for feed efficiency (FE) traits. FE-related traits were measured in 1156 Nellore cattle from an experimental breeding program genotyped for ~ 300 K markers after quality control. Prediction accuracy (Acc) was evaluated using a forward validation splitting the dataset based on birth year, considering the phenotypes adjusted for the fixed effects and covariates as pseudo-phenotypes. The MLNN and SVR approaches were trained by randomly splitting the training population into fivefold to select the best hyperparameters. The results show that the machine learning methods (MLNN and SVR) and MTGBLUP outperformed STGBLUP and the Bayesian regression approaches, increasing the Acc by approximately 8.9%, 14.6%, and 13.7% using MLNN, SVR, and MTGBLUP, respectively. Acc for SVR and MTGBLUP were slightly different, ranging from 0.62 to 0.69 and 0.62 to 0.68, respectively, with empirically unbiased for both models (0.97 and 1.09). Our results indicated that SVR and MTGBLUBP approaches were more accurate in predicting FE-related traits than Bayesian regression and STGBLUP and seemed competitive for GS of complex phenotypes with various degrees of inheritance.
Collapse
Affiliation(s)
- Lucio F M Mota
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil.
| | - Leonardo M Arikawa
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil
| | - Samuel W B Santos
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil
| | - Gerardo A Fernandes Júnior
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil
| | - Anderson A C Alves
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil
| | - Guilherme J M Rosa
- Department of Animal and Dairy Sciences, University of Wisconsin, Madison, WI, 53706, USA
| | - Maria E Z Mercadante
- Institute of Animal Science, Beef Cattle Research Center, Sertãozinho, SP, 14174-000, Brazil
- National Council for Science and Technological Development, Brasilia, DF, 71605-001, Brazil
| | - Joslaine N S G Cyrillo
- Institute of Animal Science, Beef Cattle Research Center, Sertãozinho, SP, 14174-000, Brazil
| | - Roberto Carvalheiro
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil
- National Council for Science and Technological Development, Brasilia, DF, 71605-001, Brazil
| | - Lucia G Albuquerque
- School of Agricultural and Veterinarian Sciences, São Paulo State University (UNESP), Jaboticabal, SP, 14884-900, Brazil.
- National Council for Science and Technological Development, Brasilia, DF, 71605-001, Brazil.
| |
Collapse
|
2
|
Cai W, Hu J, Fan W, Xu Y, Tang J, Xie M, Zhang Y, Guo Z, Zhou Z, Hou S. Strategies to improve genomic predictions for 35 duck carcass traits in an F 2 population. J Anim Sci Biotechnol 2023; 14:74. [PMID: 37147656 PMCID: PMC10163724 DOI: 10.1186/s40104-023-00875-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 04/02/2023] [Indexed: 05/07/2023] Open
Abstract
BACKGROUND Carcass traits are crucial for broiler ducks, but carcass traits can only be measured postmortem. Genomic selection (GS) is an effective approach in animal breeding to improve selection and reduce costs. However, the performance of genomic prediction in duck carcass traits remains largely unknown. RESULTS In this study, we estimated the genetic parameters, performed GS using different models and marker densities, and compared the estimation performance between GS and conventional BLUP on 35 carcass traits in an F2 population of ducks. Most of the cut weight traits and intestine length traits were estimated to be high and moderate heritabilities, respectively, while the heritabilities of percentage slaughter traits were dynamic. The reliability of genome prediction using GBLUP increased by an average of 0.06 compared to the conventional BLUP method. The Permutation studies revealed that 50K markers had achieved ideal prediction reliability, while 3K markers still achieved 90.7% predictive capability would further reduce the cost for duck carcass traits. The genomic relationship matrix normalized by our true variance method instead of the widely used [Formula: see text] could achieve an increase in prediction reliability in most traits. We detected most of the bayesian models had a better performance, especially for BayesN. Compared to GBLUP, BayesN can further improve the predictive reliability with an average of 0.06 for duck carcass traits. CONCLUSION This study demonstrates genomic selection for duck carcass traits is promising. The genomic prediction can be further improved by modifying the genomic relationship matrix using our proposed true variance method and several Bayesian models. Permutation study provides a theoretical basis for the fact that low-density arrays can be used to reduce genotype costs in duck genome selection.
Collapse
Affiliation(s)
- Wentao Cai
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Jian Hu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
- Shandong New Hope Liuhe Group Co., Ltd., Qingdao, 266108, China
| | - Wenlei Fan
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
- College of Animal Science and Technology, Qingdao Agricultural University, Qingdao, 266109, China
| | - Yaxi Xu
- College of Animal Science and Technology, Beijing University of Agriculture, Beijing, 102206, China
| | - Jing Tang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Ming Xie
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Yunsheng Zhang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Zhanbao Guo
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Zhengkui Zhou
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Shuisheng Hou
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| |
Collapse
|
3
|
Ogawa S, Taniguchi Y, Watanabe T, Iwaisaki H. Fitting Genomic Prediction Models with Different Marker Effects among Prefectures to Carcass Traits in Japanese Black Cattle. Genes (Basel) 2022; 14:24. [PMID: 36672767 PMCID: PMC9859149 DOI: 10.3390/genes14010024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 12/16/2022] [Accepted: 12/20/2022] [Indexed: 12/25/2022] Open
Abstract
We fitted statistical models, which assumed single-nucleotide polymorphism (SNP) marker effects differing across the fattened steers marketed into different prefectures, to the records for cold carcass weight (CW) and marbling score (MS) of 1036, 733, and 279 Japanese Black fattened steers marketed into Tottori, Hiroshima, and Hyogo prefectures in Japan, respectively. Genotype data on 33,059 SNPs was used. Five models that assume only common SNP effects to all the steers (model 1), common effects plus SNP effects differing between the steers marketed into Hyogo prefecture and others (model 2), only the SNP effects differing between Hyogo steers and others (model 3), common effects plus SNP effects specific to each prefecture (model 4), and only the effects specific to each prefecture (model 5) were exploited. For both traits, slightly lower values of residual variance than that of model 1 were estimated when fitting all other models. Estimated genetic correlation among the prefectures in models 2 and 4 ranged to 0.53 to 0.71, all <0.8. These results might support that the SNP effects differ among the prefectures to some degree, although we discussed the necessity of careful consideration to interpret the current results.
Collapse
Affiliation(s)
- Shinichiro Ogawa
- Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
- Division of Meat Animal and Poultry Research, Institute of Livestock and Grassland Science, Tsukuba 305-0901, Japan
| | - Yukio Taniguchi
- Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
| | - Toshio Watanabe
- National Livestock Breeding Center, Fukushima 961-8511, Japan
- Maebashi Institute of Animal Science, Livestock Improvement Association of Japan, Inc., Maebashi 371-0121, Japan
| | - Hiroaki Iwaisaki
- Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
- Sado Island Center for Ecological Sustainability, Niigata University, Niigata 952-0103, Japan
| |
Collapse
|
4
|
Nantongo JS, Potts BM, Klápště J, Graham NJ, Dungey HS, Fitzgerald H, O'Reilly-Wapstra JM. Genomic selection for resistance to mammalian bark stripping and associated chemical compounds in radiata pine. G3 (BETHESDA, MD.) 2022; 12:jkac245. [PMID: 36218439 PMCID: PMC9635650 DOI: 10.1093/g3journal/jkac245] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 08/29/2022] [Indexed: 07/28/2023]
Abstract
The integration of genomic data into genetic evaluations can facilitate the rapid selection of superior genotypes and accelerate the breeding cycle in trees. In this study, 390 trees from 74 control-pollinated families were genotyped using a 36K Axiom SNP array. A total of 15,624 high-quality SNPs were used to develop genomic prediction models for mammalian bark stripping, tree height, and selected primary and secondary chemical compounds in the bark. Genetic parameters from different genomic prediction methods-single-trait best linear unbiased prediction based on a marker-based relationship matrix (genomic best linear unbiased prediction), multitrait single-step genomic best linear unbiased prediction, which integrated the marker-based and pedigree-based relationship matrices (single-step genomic best linear unbiased prediction) and the single-trait generalized ridge regression-were compared to equivalent single- or multitrait pedigree-based approaches (ABLUP). The influence of the statistical distribution of data on the genetic parameters was assessed. Results indicated that the heritability estimates were increased nearly 2-fold with genomic models compared to the equivalent pedigree-based models. Predictive accuracy of the single-step genomic best linear unbiased prediction was higher than the ABLUP for most traits. Allowing for heterogeneity in marker effects through the use of generalized ridge regression did not markedly improve predictive ability over genomic best linear unbiased prediction, arguing that most of the chemical traits are modulated by many genes with small effects. Overall, the traits with low pedigree-based heritability benefited more from genomic models compared to the traits with high pedigree-based heritability. There was no evidence that data skewness or the presence of outliers affected the genomic or pedigree-based genetic estimates.
Collapse
Affiliation(s)
- Judith S Nantongo
- Corresponding author: National Agricultural Research Organization, P.O Box 1752, Mukono, Uganda.
| | - Brad M Potts
- School of Natural Sciences, University of Tasmania, Hobart, TAS 7001, Australia
- ARC Training Centre for Forest Value, Hobart, TAS 7001, Australia
| | - Jaroslav Klápště
- Scion (New Zealand Forest Research Institute Ltd.), Rotorua 3046, New Zealand
| | - Natalie J Graham
- Scion (New Zealand Forest Research Institute Ltd.), Rotorua 3046, New Zealand
| | - Heidi S Dungey
- Scion (New Zealand Forest Research Institute Ltd.), Rotorua 3046, New Zealand
| | - Hugh Fitzgerald
- School of Natural Sciences, University of Tasmania, Hobart, TAS 7001, Australia
| | - Julianne M O'Reilly-Wapstra
- School of Natural Sciences, University of Tasmania, Hobart, TAS 7001, Australia
- ARC Training Centre for Forest Value, Hobart, TAS 7001, Australia
| |
Collapse
|
5
|
Fathoni A, Boonkum W, Chankitisakul V, Duangjinda M. An Appropriate Genetic Approach for Improving Reproductive Traits in Crossbred Thai-Holstein Cattle under Heat Stress Conditions. Vet Sci 2022; 9:vetsci9040163. [PMID: 35448661 PMCID: PMC9031002 DOI: 10.3390/vetsci9040163] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Revised: 03/19/2022] [Accepted: 03/26/2022] [Indexed: 01/16/2023] Open
Abstract
Thailand is a tropical country affected by global climate change and has high temperatures and humidity that cause heat stress in livestock. A temperature−humidity index (THI) is required to assess and evaluate heat stress levels in livestock. One of the livestock types in Thailand experiencing heat stress due to extreme climate change is crossbred dairy cattle. Genetic evaluations of heat tolerance in dairy cattle have been carried out for reproductive traits. Heritability values for reproductive traits are generally low (<0.10) because environmental factors heavily influence them. Consequently, genetic improvement for these traits would be slow compared to production traits. Positive and negative genetic correlations were found between reproductive traits and reproductive traits and yield traits. Several selection methods for reproductive traits have been introduced, i.e., the traditional method, marker-assisted selection (MAS), and genomic selection (GS). GS is the most promising technique and provides accurate results with a high genetic gain. Single-step genomic BLUP (ssGBLUP) has higher accuracy than the multi-step equivalent for fertility traits or low-heritability traits.
Collapse
Affiliation(s)
- Akhmad Fathoni
- Department of Animal Science, Faculty of Agriculture, Khon Kaen University, Khon Kaen 40002, Thailand; (A.F.); (W.B.); (V.C.)
- Department of Animal Breeding and Reproduction, Faculty of Animal Science, Universitas Gadjah Mada, Yogyakarta 55281, Indonesia
| | - Wuttigrai Boonkum
- Department of Animal Science, Faculty of Agriculture, Khon Kaen University, Khon Kaen 40002, Thailand; (A.F.); (W.B.); (V.C.)
- Network Center for Animal Breeding and OMICS Research, Khon Kaen University, Khon Kaen 40002, Thailand
| | - Vibuntita Chankitisakul
- Department of Animal Science, Faculty of Agriculture, Khon Kaen University, Khon Kaen 40002, Thailand; (A.F.); (W.B.); (V.C.)
- Network Center for Animal Breeding and OMICS Research, Khon Kaen University, Khon Kaen 40002, Thailand
| | - Monchai Duangjinda
- Department of Animal Science, Faculty of Agriculture, Khon Kaen University, Khon Kaen 40002, Thailand; (A.F.); (W.B.); (V.C.)
- Network Center for Animal Breeding and OMICS Research, Khon Kaen University, Khon Kaen 40002, Thailand
- Correspondence: ; Tel.: +66-81-872-4207
| |
Collapse
|
6
|
Mehrban H, Naserkheil M, Lee D, Ibáñez-Escriche N. Multi-Trait Single-Step GBLUP Improves Accuracy of Genomic Prediction for Carcass Traits Using Yearling Weight and Ultrasound Traits in Hanwoo. Front Genet 2021; 12:692356. [PMID: 34394186 PMCID: PMC8363309 DOI: 10.3389/fgene.2021.692356] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 07/02/2021] [Indexed: 11/25/2022] Open
Abstract
There has been a growing interest in the genetic improvement of carcass traits as an important and primary breeding goal in the beef cattle industry over the last few decades. The use of correlated traits and molecular information can aid in obtaining more accurate estimates of breeding values. This study aimed to assess the improvement in the accuracy of genetic predictions for carcass traits by using ultrasound measurements and yearling weight along with genomic information in Hanwoo beef cattle by comparing four evaluation models using the estimators of the recently developed linear regression method. We compared the performance of single-trait pedigree best linear unbiased prediction [ST-BLUP and single-step genomic (ST-ssGBLUP)], as well as multi-trait (MT-BLUP and MT-ssGBLUP) models for the studied traits at birth and yearling date of steers. The data comprised of 15,796 phenotypic records for yearling weight and ultrasound traits as well as 5,622 records for carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score), resulting in 43,949 single-nucleotide polymorphisms from 4,284 steers and 2,332 bulls. Our results indicated that averaged across all traits, the accuracy of ssGBLUP models (0.52) was higher than that of pedigree-based BLUP (0.34), regardless of the use of single- or multi-trait models. On average, the accuracy of prediction can be further improved by implementing yearling weight and ultrasound data in the MT-ssGBLUP model (0.56) for the corresponding carcass traits compared to the ST-ssGBLUP model (0.49). Moreover, this study has shown the impact of genomic information and correlated traits on predictions at the yearling date (0.61) using MT-ssGBLUP models, which was advantageous compared to predictions at birth date (0.51) in terms of accuracy. Thus, using genomic information and high genetically correlated traits in the multi-trait model is a promising approach for practical genomic selection in Hanwoo cattle, especially for traits that are difficult to measure.
Collapse
Affiliation(s)
- Hossein Mehrban
- Department of Animal Science, Shahrekord University, Shahrekord, Iran
| | - Masoumeh Naserkheil
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran.,Department of Animal Life and Environment Sciences, Hankyong National University, Gyeonggi-do, South Korea
| | - Deukhwan Lee
- Department of Animal Life and Environment Sciences, Hankyong National University, Gyeonggi-do, South Korea
| | - Noelia Ibáñez-Escriche
- Institute for Animal Science and Technology, Universitat Politècnica de València, València, Spain
| |
Collapse
|
7
|
Martinez-Castillero M, Then C, Altarriba J, Srihi H, López-Carbonell D, Díaz C, Martinez P, Hermida M, Varona L. Detection of Genomic Regions with Pleiotropic Effects for Growth and Carcass Quality Traits in the Rubia Gallega Cattle Breed. Animals (Basel) 2021; 11:ani11061682. [PMID: 34200089 PMCID: PMC8227173 DOI: 10.3390/ani11061682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Revised: 05/25/2021] [Accepted: 06/02/2021] [Indexed: 11/16/2022] Open
Abstract
Simple Summary The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. We have developed a ssGWAS by backsolving the SNP effects after implementing a ssGBLUP. The results showed an apparent heterogeneity of the additive genetic variance across the genome. Some of the genomic regions explaining the most of this additive variance were shared across traits, indicating the presence of pleiotropic effects, which were reflected in their genetic correlations. Abstract The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. In recent years, genomic evaluation has been implemented by using a ssGBLUP (single-step Genomic Best Linear Unbiased Prediction). This procedure can reparameterized to perform ssGWAS (single-step Genome Wide Association Studies) by backsolving the SNP (single nucleotide polymorphisms) effects. Therefore, the objective of this study was to identify genomic regions associated with the genetic variability in growth and carcass quality traits. We implemented a ssGBLUP by using a database that included records for Birth Weight (BW-327,350 records-), Weaning Weight (WW-83,818-), Cold Carcass Weight (CCW-91,621-), Fatness (FAT-91,475-) and Conformation (CON-91,609-). The pedigree included 464,373 individuals, 2449 of which were genotyped. After a process of filtering, we ended up using 43,211 SNP markers. We used the GBLUP and SNPBLUP model equivalences to obtain the effects of the SNPs and then calculated the percentage of variance explained by the regions of the genome between 1 Mb. We identified 7 regions of the genome for CCW; 8 regions for BW, WW, FAT and 9 regions for CON, which explained the percentage of variance above 0.5%. Furthermore, a number of the genome regions had pleiotropic effects, located at: BTA1 (131–132 Mb), BTA2 (1–11 Mb), BTA3 (32–33 Mb), BTA6 (36–38 Mb), BTA16 (24–26 Mb), and BTA 21 (56–57 Mb). These regions contain, amongst others, the following candidate genes: NCK1, MSTN, KCNA3, LCORL, NCAPG, and RIN3.
Collapse
Affiliation(s)
- Maria Martinez-Castillero
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
- Correspondence:
| | - Carlos Then
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
| | - Juan Altarriba
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
| | - Houssemeddine Srihi
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
| | - David López-Carbonell
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
| | - Clara Díaz
- Instituto Nacional de Investigación y Tecnología Agraria (INIA), 28040 Madrid, Spain;
| | - Paulino Martinez
- Facultad de Veterinaria, Universidad de Santiago de Compostela, 27002 Lugo, Spain; (P.M.); (M.H.)
| | - Miguel Hermida
- Facultad de Veterinaria, Universidad de Santiago de Compostela, 27002 Lugo, Spain; (P.M.); (M.H.)
| | - Luis Varona
- Instituto Agroalimentario de Aragón (IA2), Universidad de Zaragoza, 50013 Zaragoza, Spain; (C.T.); (J.A.); (H.S.); (D.L.-C.); (L.V.)
| |
Collapse
|
8
|
Karaman E, Su G, Croue I, Lund MS. Genomic prediction using a reference population of multiple pure breeds and admixed individuals. Genet Sel Evol 2021; 53:46. [PMID: 34058971 PMCID: PMC8168010 DOI: 10.1186/s12711-021-00637-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 05/11/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND In dairy cattle populations in which crossbreeding has been used, animals show some level of diversity in their origins. In rotational crossbreeding, for instance, crossbred dams are mated with purebred sires from different pure breeds, and the genetic composition of crossbred animals is an admixture of the breeds included in the rotation. How to use the data of such individuals in genomic evaluations is still an open question. In this study, we aimed at providing methodologies for the use of data from crossbred individuals with an admixed genetic background together with data from multiple pure breeds, for the purpose of genomic evaluations for both purebred and crossbred animals. A three-breed rotational crossbreeding system was mimicked using simulations based on animals genotyped with the 50 K single nucleotide polymorphism (SNP) chip. RESULTS For purebred populations, within-breed genomic predictions generally led to higher accuracies than those from multi-breed predictions using combined data of pure breeds. Adding admixed population's (MIX) data to the combined pure breed data considering MIX as a different breed led to higher accuracies. When prediction models were able to account for breed origin of alleles, accuracies were generally higher than those from combining all available data, depending on the correlation of quantitative trait loci (QTL) effects between the breeds. Accuracies varied when using SNP effects from any of the pure breeds to predict the breeding values of MIX. Using those breed-specific SNP effects that were estimated separately in each pure breed, while accounting for breed origin of alleles for the selection candidates of MIX, generally improved the accuracies. Models that are able to accommodate MIX data with the breed origin of alleles approach generally led to higher accuracies than models without breed origin of alleles, depending on the correlation of QTL effects between the breeds. CONCLUSIONS Combining all available data, pure breeds' and admixed population's data, in a multi-breed reference population is beneficial for the estimation of breeding values for pure breeds with a small reference population. For MIX, such an approach can lead to higher accuracies than considering breed origin of alleles for the selection candidates, and using breed-specific SNP effects estimated separately in each pure breed. Including MIX data in the reference population of multiple breeds by considering the breed origin of alleles, accuracies can be further improved. Our findings are relevant for breeding programs in which crossbreeding is systematically applied, and also for populations that involve different subpopulations and between which exchange of genetic material is routine practice.
Collapse
Affiliation(s)
- Emre Karaman
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark.
| | - Guosheng Su
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark
| | | | - Mogens S Lund
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark
| |
Collapse
|
9
|
Mehrban H, Naserkheil M, Lee DH, Cho C, Choi T, Park M, Ibáñez-Escriche N. Genomic Prediction Using Alternative Strategies of Weighted Single-Step Genomic BLUP for Yearling Weight and Carcass Traits in Hanwoo Beef Cattle. Genes (Basel) 2021; 12:genes12020266. [PMID: 33673102 PMCID: PMC7917987 DOI: 10.3390/genes12020266] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Revised: 02/08/2021] [Accepted: 02/11/2021] [Indexed: 01/20/2023] Open
Abstract
The weighted single-step genomic best linear unbiased prediction (GBLUP) method has been proposed to exploit information from genotyped and non-genotyped relatives, allowing the use of weights for single-nucleotide polymorphism in the construction of the genomic relationship matrix. The purpose of this study was to investigate the accuracy of genetic prediction using the following single-trait best linear unbiased prediction methods in Hanwoo beef cattle: pedigree-based (PBLUP), un-weighted (ssGBLUP), and weighted (WssGBLUP) single-step genomic methods. We also assessed the impact of alternative single and window weighting methods according to their effects on the traits of interest. The data was comprised of 15,796 phenotypic records for yearling weight (YW) and 5622 records for carcass traits (backfat thickness: BFT, carcass weight: CW, eye muscle area: EMA, and marbling score: MS). Also, the genotypic data included 6616 animals for YW and 5134 for carcass traits on the 43,950 single-nucleotide polymorphisms. The ssGBLUP showed significant improvement in genomic prediction accuracy for carcass traits (71%) and yearling weight (99%) compared to the pedigree-based method. The window weighting procedures performed better than single SNP weighting for CW (11%), EMA (11%), MS (3%), and YW (6%), whereas no gain in accuracy was observed for BFT. Besides, the improvement in accuracy between window WssGBLUP and the un-weighted method was low for BFT and MS, while for CW, EMA, and YW resulted in a gain of 22%, 15%, and 20%, respectively, which indicates the presence of relevant quantitative trait loci for these traits. These findings indicate that WssGBLUP is an appropriate method for traits with a large quantitative trait loci effect.
Collapse
Affiliation(s)
- Hossein Mehrban
- Department of Animal Science, Shahrekord University, Shahrekord 88186-34141, Iran;
| | - Masoumeh Naserkheil
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj 77871-31587, Iran;
- Department of Animal Life and Environment Sciences, Hankyong National University, Jungang-ro 327, Anseong-si 17579, Gyeonggi-do, Korea
| | - Deuk Hwan Lee
- Department of Animal Life and Environment Sciences, Hankyong National University, Jungang-ro 327, Anseong-si 17579, Gyeonggi-do, Korea
- Correspondence: ; Tel.: +82-316-705-091
| | - Chungil Cho
- Hanwoo Genetic Improvement Center, NongHyup Agribusiness Group Inc., Seosan 31948, Korea;
| | - Taejeong Choi
- Animal Breeding and Genetics Division, National Institute of Animal Science, Rural Development Administration, Cheonan 31000, Korea; (T.C.); (M.P.)
| | - Mina Park
- Animal Breeding and Genetics Division, National Institute of Animal Science, Rural Development Administration, Cheonan 31000, Korea; (T.C.); (M.P.)
| | - Noelia Ibáñez-Escriche
- Institute for Animal Science and Technology, Universitat Politècnica de València, 46022 València, Spain;
| |
Collapse
|
10
|
Lopes FB, Baldi F, Passafaro TL, Brunes LC, Costa MFO, Eifert EC, Narciso MG, Rosa GJM, Lobo RB, Magnabosco CU. Genome-enabled prediction of meat and carcass traits using Bayesian regression, single-step genomic best linear unbiased prediction and blending methods in Nelore cattle. Animal 2020; 15:100006. [PMID: 33516009 DOI: 10.1016/j.animal.2020.100006] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2020] [Revised: 06/05/2020] [Accepted: 06/08/2020] [Indexed: 10/22/2022] Open
Abstract
Several methods have been used for genome-enabled prediction (or genomic selection) of complex traits, for example, multiple regression models describing a target trait with a linear function of a set of genetic markers. Genomic selection studies have been focused mostly on single-trait analyses. However, most profitability traits are genetically correlated, and an increase in prediction accuracy of genomic breeding values for genetically correlated traits is expected when using multiple-trait models. Thus, this study was carried out to assess the accuracy of genomic prediction for carcass and meat quality traits in Nelore cattle, using single- and multiple-trait approaches. The study considered 15 780, 15 784, 15 742 and 526 records of rib eye area (REA, cm2), back fat thickness (BF, mm), rump fat (RF, mm) and Warner-Bratzler shear force (WBSF, kg), respectively, in Nelore cattle, from the Nelore Brazil Breeding Program. Animals were genotyped with a low-density single nucleotide polymorphism (SNP) panel and subsequently imputed to arrays with 54 and 777 k SNPs. Four Bayesian specifications of genomic regression models, namely, Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression; blending methods, BLUP; and single-step genomic best linear unbiased prediction (ssGBLUP) methods were compared in terms of prediction accuracy using a fivefold cross-validation. Estimates of heritability ranged from 0.20 to 0.35 and from 0.21 to 0.46 for RF and WBSF on single- and multiple-trait analyses, respectively. Prediction accuracies for REA, BF, RF and WBSF were all similar using the different specifications of regression models. In addition, this study has shown the impact of genomic information upon genetic evaluations in beef cattle using the multiple-trait model, which was also advantageous compared to the single-trait model because it accounted for the selection process using multiple traits at the same time. The advantage of multi-trait analyses is attributed to the consideration of correlations and genetic influences between the traits, in addition to the non-random association of alleles.
Collapse
Affiliation(s)
- F B Lopes
- Department of Animal Science, São Paulo State University - Júlio de Mesquita Filho (UNESP), Access way Prof. Paulo Donato Castelane, Jaboticabal, SP 14884-900, Brazil; Embrapa Cerrados, BR-020, 18, Sobradinho, Brasilia, DF 70770-901, Brazil.
| | - F Baldi
- Department of Animal Science, São Paulo State University - Júlio de Mesquita Filho (UNESP), Access way Prof. Paulo Donato Castelane, Jaboticabal, SP 14884-900, Brazil
| | - T L Passafaro
- Department of Animal Sciences, University of Wisconsin-Madison, Madison, WI 53705, USA
| | - L C Brunes
- Department of Animal Science, Federal University of Goiás, Goiânia, GO 75345-000, Brazil
| | - M F O Costa
- Embrapa Rice and Beans, GO-462, km 12, Santo Antônio de Goiás, GO 75375-000, Brazil
| | - E C Eifert
- Embrapa Cerrados, BR-020, 18, Sobradinho, Brasilia, DF 70770-901, Brazil
| | - M G Narciso
- Embrapa Rice and Beans, GO-462, km 12, Santo Antônio de Goiás, GO 75375-000, Brazil
| | - G J M Rosa
- Department of Animal Sciences, University of Wisconsin-Madison, Madison, WI 53705, USA; Department of Biostatistics and Medical Informatics, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - R B Lobo
- National Association of Breeders and Researchers, Ribeirão Preto SP 14020-230, Brazil
| | - C U Magnabosco
- Embrapa Cerrados, BR-020, 18, Sobradinho, Brasilia, DF 70770-901, Brazil
| |
Collapse
|
11
|
Alvarenga AB, Veroneze R, Oliveira HR, Marques DBD, Lopes PS, Silva FF, Brito LF. Comparing Alternative Single-Step GBLUP Approaches and Training Population Designs for Genomic Evaluation of Crossbred Animals. Front Genet 2020; 11:263. [PMID: 32328083 PMCID: PMC7162606 DOI: 10.3389/fgene.2020.00263] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Accepted: 03/05/2020] [Indexed: 02/06/2023] Open
Abstract
As crossbreeding is extensively used in some livestock species, we aimed to evaluate the performance of single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) methods to predict Genomic Estimated Breeding Values (GEBVs) of crossbred animals. Different training population scenarios were evaluated: (SC1) ssGBLUP based on a single-trait model considering purebred and crossbred animals in a joint training population; (SC2) ssGBLUP based on a multiple-trait model to enable considering phenotypes recorded in purebred and crossbred training animals as different traits; (SC3) WssGBLUP based on a single-trait model considering purebred and crossbred animals jointly in the training population (both populations were used for SNP weights' estimation); (SC4) WssGBLUP based on a single-trait model considering only purebred animals in the training population (crossbred population only used for SNP weights' estimation); (SC5) WssGBLUP based on a single-trait model and the training population characterized by purebred animals (purebred population used for SNP weights' estimation). A complex trait was simulated assuming alternative genetic architectures. Different scaling factors to blend the inverse of the genomic (G -1) and pedigree (A 22 - 1 ) relationship matrices were also tested. The predictive performance of each scenario was evaluated based on the validation accuracy and regression coefficient. The genetic correlations across simulated populations in the different scenarios ranged from moderate to high (0.71-0.99). The scenario mimicking a completely polygenic trait (h Q T L 2 = 0) yielded the lowest validation accuracy (0.12; for SC3 and SC4). The simulated scenarios assuming 4,500 QTLs affecting the trait andh Q T L 2 = h 2 resulted in the greatest GEBV accuracies (0.47; for SC1 and SC2). The regression coefficients ranged from 0.28 (for SC3 assuming polygenic effect) to 1.27 (for SC2 considering 4,500 QTLs). In general, SC3 and SC5 resulted in inflated GEBVs, whereas other scenarios yielded deflated GEBVs. The scaling factors used to combine G -1 andA 22 - 1 had a small influence on the validation accuracies, but a greater effect on the regression coefficients. Due to the complexity of multiple-trait models and WssGBLUP analyses, and a similar predictive performance across the methods evaluated, SC1 is recommended for genomic evaluation in crossbred populations with similar genetic structures [moderate-to-high (0.71-0.99) genetic correlations between purebred and crossbred populations].
Collapse
Affiliation(s)
- Amanda B. Alvarenga
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- Department of Animal Science, Federal University of Viçosa, Viçosa, Brazil
| | - Renata Veroneze
- Department of Animal Science, Federal University of Viçosa, Viçosa, Brazil
| | - Hinayah R. Oliveira
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | | | - Paulo S. Lopes
- Department of Animal Science, Federal University of Viçosa, Viçosa, Brazil
| | - Fabyano F. Silva
- Department of Animal Science, Federal University of Viçosa, Viçosa, Brazil
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| |
Collapse
|