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
|
Malheiros JM, Correia BSB, Ceribeli C, Bruscadin JJ, Diniz WJS, Banerjee P, da Silva Vieira D, Cardoso TF, Andrade BGN, Petrini J, Cardoso DR, Colnago LA, Bogusz Junior S, Mourão GB, Coutinho LL, Palhares JCP, de Medeiros SR, Berndt A, de Almeida Regitano LC. Ruminal and feces metabolites associated with feed efficiency, water intake and methane emission in Nelore bulls. Sci Rep 2023; 13:18001. [PMID: 37865691 PMCID: PMC10590413 DOI: 10.1038/s41598-023-45330-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Accepted: 10/18/2023] [Indexed: 10/23/2023] Open
Abstract
The objectives of this study were twofold: (1) to identify potential differences in the ruminal and fecal metabolite profiles of Nelore bulls under different nutritional interventions; and (2) to identify metabolites associated with cattle sustainability related-traits. We used different nutritional interventions in the feedlot: conventional (Conv; n = 26), and by-product (ByPr, n = 26). Thirty-eight ruminal fluid and 27 fecal metabolites were significantly different (P < 0.05) between the ByPr and Conv groups. Individual dry matter intake (DMI), residual feed intake (RFI), observed water intake (OWI), predicted water intake (WI), and residual water intake (RWI) phenotypes were lower (P < 0.05) in the Conv group, while the ByPr group exhibited lower methane emission (ME) (P < 0.05). Ruminal fluid dimethylamine was significantly associated (P < 0.05) with DMI, RFI, FE (feed efficiency), OWI and WI. Aspartate was associated (P < 0.05) with DMI, RFI, FE and WI. Fecal C22:1n9 was significantly associated with OWI and RWI (P < 0.05). Fatty acid C14:0 and hypoxanthine were significantly associated with DMI and RFI (P < 0.05). The results demonstrated that different nutritional interventions alter ruminal and fecal metabolites and provided new insights into the relationship of these metabolites with feed efficiency and water intake traits in Nelore bulls.
Collapse
Affiliation(s)
| | | | - Caroline Ceribeli
- Institute of Chemistry, University of São Paulo/USP, São Carlos, São Paulo, Brazil
- Department of Food Science, University of Copenhagen, Copenhagen, Denmark
| | | | - Wellison J S Diniz
- Departament of Animal Sciences, Auburn University, Auburn, AL, 36849, USA
| | - Priyanka Banerjee
- Departament of Animal Sciences, Auburn University, Auburn, AL, 36849, USA
| | | | | | - Bruno Gabriel Nascimento Andrade
- Embrapa Southeast Livestock, São Carlos, São Paulo, Brazil
- Computer Science Department, Munster Technological University, MTU/ADAPT, Cork, Ireland
| | - Juliana Petrini
- Department of Animal Science, University of São Paulo/ESALQ, Piracicaba, São Paulo, Brazil
| | | | | | | | - Gerson Barreto Mourão
- Department of Animal Science, University of São Paulo/ESALQ, Piracicaba, São Paulo, Brazil
| | - Luiz Lehmann Coutinho
- Department of Animal Science, University of São Paulo/ESALQ, Piracicaba, São Paulo, Brazil
| | | | | | | | | |
Collapse
|
3
|
Faggion S, Carnier P, Franch R, Babbucci M, Pascoli F, Dalla Rovere G, Caggiano M, Chavanne H, Toffan A, Bargelloni L. Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings. Genet Sel Evol 2023; 55:22. [PMID: 37013478 PMCID: PMC10069116 DOI: 10.1186/s12711-023-00796-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Accepted: 03/21/2023] [Indexed: 04/05/2023] Open
Abstract
BACKGROUND The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers. RESULTS Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively). CONCLUSIONS The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.
Collapse
Affiliation(s)
- Sara Faggion
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy.
| | - Paolo Carnier
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy
| | - Rafaella Franch
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy
| | - Massimiliano Babbucci
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy
| | - Francesco Pascoli
- Division of Comparative Biomedical Sciences, OIE Reference Centre for Viral Encephalopathy and Retinopathy, Istituto Zooprofilattico Sperimentale delle Venezie (IZSVe), Padova, Italy
| | - Giulia Dalla Rovere
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy
| | - Massimo Caggiano
- Panittica Italia Società Agricola S.R.L., Strada del Procaccio, 72016, Torre Canne di Fasano, Italy
| | - Hervé Chavanne
- Panittica Italia Società Agricola S.R.L., Strada del Procaccio, 72016, Torre Canne di Fasano, Italy
| | - Anna Toffan
- Division of Comparative Biomedical Sciences, OIE Reference Centre for Viral Encephalopathy and Retinopathy, Istituto Zooprofilattico Sperimentale delle Venezie (IZSVe), Padova, Italy
| | - Luca Bargelloni
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università, 16, 35020, Legnaro, PD, Italy
| |
Collapse
|
4
|
Brunes LC, de Faria CU, Magnabosco CU, Lobo RB, Peripolli E, Aguilar I, Baldi F. Genomic prediction ability and genetic parameters for residual feed intake calculated using different approaches and their associations with growth, reproductive, and carcass traits in Nellore cattle. J Appl Genet 2023; 64:159-167. [PMID: 36376720 DOI: 10.1007/s13353-022-00734-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Revised: 09/03/2022] [Accepted: 10/26/2022] [Indexed: 11/16/2022]
Abstract
This study aimed to estimate prediction ability and genetic parameters for residual feed intake (RFI) calculated using a regression equation for each test (RFItest) and for the whole population (RFIpop) in Nellore beef cattle. It also aimed to evaluate the correlations between RFIpop and RFItest with growth, reproductive, and carcass traits. Genotypic and phenotypic records from 8354 animals were used. An analysis of variance (ANOVA) was performed to verify the adequacy of the regression equations applied to estimate the RFItest and RFIpop. The (co)variance components were obtained using the single-step genomic best linear unbiased prediction under single and two-trait animal model analyses. The genetic and phenotypic correlations between RFItest and RFIpop with dry matter intake, frame, growth, reproduction, and carcass-related traits were evaluated. The prediction ability and bias were estimated to compare the RFItest and RFIpop genomic breeding values (GEBV). The RFIpop ANOVA showed a higher significance level (p < 0.0001) than did the RFItest for the fixed effects. The RFIpop displayed higher additive genetic variance estimated than the RFItest, although the RFIpop and RFItest displayed similar heritabilities. Overall, the RFItest showed higher residual correlations with growth, reproductive, and carcass traits, while the RFIpop displayed higher genetic correlations with such traits. The GEBV for the RFItest was slightly biased than GEBV RFIpop. The approach to calculate the RFI influenced the decomposition and estimation of variance components and genomic prediction for RFI. The application of RFIpop would be more appropriate for genetic evaluation purpose to adjust or correct for non-genetic effects and to decrease the prediction bias for RFI.
Collapse
Affiliation(s)
| | | | | | | | - Elisa Peripolli
- Departament of Animal Science, College of Agricultural and Veterinary Sciences, Sao Paulo State University (UNESP), Jaboticabal, 14884-900, Brazil
| | - Ignacio Aguilar
- Instituto Nacional de Investigación Agropecuaria (INIA), 11500, Montevideo, Uruguay
| | - Fernando Baldi
- Departament of Animal Science, College of Agricultural and Veterinary Sciences, Sao Paulo State University (UNESP), Jaboticabal, 14884-900, Brazil
| |
Collapse
|
5
|
Carrara ER, Peixoto MGCD, da Silva AA, Bruneli FAT, Ventura HT, Zadra LEF, Josahkian LA, Veroneze R, Lopes PS. Genomic prediction in Brazilian Guzerá cattle: application of a single-step approach to productive and reproductive traits. Trop Anim Health Prod 2023; 55:48. [PMID: 36705782 DOI: 10.1007/s11250-023-03484-9] [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: 09/13/2022] [Accepted: 01/23/2023] [Indexed: 01/28/2023]
Abstract
This study aimed to investigate the feasibility of genomic prediction for productive and reproductive traits in Guzerá cattle using single-step genomic best linear unbiased prediction (ssGBLUP). Evaluations included the 305-day cumulative yields (first lactation, in kg) of milk, lactose, protein, fat, and total solids; adjusted body weight (kg) at the ages of 450, 365, and 210 days; and age at first calving (in days), from a database containing 197,283 measurements from Guzerá males and females born between 1954 and 2018. The pedigree included 433,823 animals spanning up to 14 overlapping generations. A total of 1618 animals were genotyped. The analyses were performed using ssGBLUP and traditional BLUP methods. Predictive ability and bias were accessed using cross-validation: predictive ability was similar between the methods and ranged from 0.27 to 0.47 for the genomic-based model and from 0.30 to 0.45 for the pedigree-based model; the bias was also similar between the methods, ranging from 0.88 to 1.35 in the genomic-based model and from 0.96 to 1.41 in the pedigree-based model. The individual accuracies of breeding values were evidently increased in the genomic evaluation, with values ranging from 0.41 to 0.56 in the genomic-based model and from 0.26 to 0.54 in the pedigree-based model. Even based on a small number of genotyped animals and a small database for some traits, the results suggest that ssGBLUP is feasible and may be applied to national genetic evaluation of the breed to increase the accuracy of breeding values without greatly impacting predictive ability and bias.
Collapse
Affiliation(s)
- Eula Regina Carrara
- Department of Animal Science, Federal University of Viçosa, Viçosa, Minas Gerais, Brazil.
| | | | - Alessandra Alves da Silva
- Department of Agricultural Sciences, School of Agricultural and Veterinarian Sciences, São Paulo State University, Jaboticabal, São Paulo, Brazil
| | | | | | - Lenira El Faro Zadra
- Brazilian Center for the Genetic Improvement of Guzerá, Belo Horizonte, Minas Gerais, Brazil
| | | | - Renata Veroneze
- Department of Animal Science, Federal University of Viçosa, Viçosa, Minas Gerais, Brazil
| | - Paulo Sávio Lopes
- Department of Animal Science, Federal University of Viçosa, Viçosa, Minas Gerais, Brazil
| |
Collapse
|
6
|
Ding C, Weng Y, Byram TD, Bartlett BD, Raley EM. Post hoc experimental designs improve genetic trial analyses: A case study of cherrybark oak (Quercus pagoda Raf.) genetic evaluation in the western Gulf region, USA. PLoS One 2023; 18:e0285150. [PMID: 37172062 PMCID: PMC10180598 DOI: 10.1371/journal.pone.0285150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 04/18/2023] [Indexed: 05/14/2023] Open
Abstract
Oaks (Quercus spp.) are widespread hardwood trees in the Northern Hemisphere and of high ecological, economic, and social values. Optimal experimental design of genetic trials is essential for accurate estimates of genetic parameters and improving the genetic merit of breeding stock. Here, we evaluate the use of post hoc row-column factors combined with spatial adjustment to improve genetic analyses of parents and individual trees in field progeny tests of plantation hardwoods, using cherrybark oak (Quercus pagoda Raf.) as an example. For tree height, post hoc incomplete blocking reduced ~14% more of the within-block environmental variance compared to the randomized complete block design (RCBD) model. Incomplete blocking also improved the heritability estimates for height by 7% to 14% compared to the original RCBD model. No clinal trend for growth breeding values was identified due to provenances. Our approach warrants the initial selection for height as early as age ~10 based on its moderate narrow-sense heritability of 0.2; however, diameter and volume need longer evaluation times. The post hoc incomplete blocking is more efficient and promising to improve the genetic analysis of Q. pagoda to minimize the environmental heterogeneity influences. Adjusting competition and spatial effects, including the distance principal components and autoregressive residual structure notably improves the model fit based on the observed reductions in AICs and BICs. Employing our approach is promising for hardwood genetic improvement in the southern USA.
Collapse
Affiliation(s)
- Chen Ding
- Western Gulf Forest Tree Improvement Program, Texas A&M Forest Service, Texas A&M University System, College Station, Texas, United States of America
- College of Forestry, Wildlife and Environment, Auburn University, Auburn, Alabama, United States of America
| | - Yuhui Weng
- Arthur Temple College of Forestry and Agriculture, Stephen F. Austin State University, Nacogdoches, Texas, United States of America
| | - Tom D Byram
- Western Gulf Forest Tree Improvement Program, Texas A&M Forest Service, Texas A&M University System, College Station, Texas, United States of America
| | - Benjamin D Bartlett
- Western Gulf Forest Tree Improvement Program, Texas A&M Forest Service, Texas A&M University System, College Station, Texas, United States of America
| | - Earl M Raley
- Western Gulf Forest Tree Improvement Program, Texas A&M Forest Service, Texas A&M University System, College Station, Texas, United States of America
| |
Collapse
|
7
|
Tahir MS, Porto-Neto LR, Reverter-Gomez T, Olasege BS, Sajid MR, Wockner KB, Tan AWL, Fortes MRS. Utility of multi-omics data to inform genomic prediction of heifer fertility traits. J Anim Sci 2022; 100:skac340. [PMID: 36239447 PMCID: PMC9733504 DOI: 10.1093/jas/skac340] [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: 06/18/2022] [Accepted: 10/12/2022] [Indexed: 12/15/2022] Open
Abstract
Biologically informed single nucleotide polymorphisms (SNPs) impact genomic prediction accuracy of the target traits. Our previous genomics, proteomics, and transcriptomics work identified candidate genes related to puberty and fertility in Brahman heifers. We aimed to test this biological information for capturing heritability and predicting heifer fertility traits in another breed i.e., Tropical Composite. The SNP from the identified genes including 10 kilobases (kb) region on either side were selected as biologically informed SNP set. The SNP from the rest of the Bos taurus genes including 10-kb region on either side were selected as biologically uninformed SNP set. Bovine high-density (HD) complete SNP set (628,323 SNP) was used as a control. Two populations-Tropical Composites (N = 1331) and Brahman (N = 2310)-had records for three traits: pregnancy after first mating season (PREG1, binary), first conception score (FCS, score 1 to 3), and rebreeding score (REB, score 1 to 3.5). Using the best linear unbiased prediction method, effectiveness of each SNP set to predict the traits was tested in two scenarios: a 5-fold cross-validation within Tropical Composites using biological information from Brahman studies, and application of prediction equations from one breed to the other. The accuracy of prediction was calculated as the correlation between genomic estimated breeding values and adjusted phenotypes. Results show that biologically informed SNP set estimated heritabilities not significantly better than the control HD complete SNP set in Tropical Composites; however, it captured all the observed genetic variance in PREG1 and FCS when modeled together with the biologically uninformed SNP set. In 5-fold cross-validation within Tropical Composites, the biologically informed SNP set performed marginally better (statistically insignificant) in terms of prediction accuracies (PREG1: 0.20, FCS: 0.13, and REB: 0.12) as compared to HD complete SNP set (PREG1: 0.17, FCS: 0.10, and REB: 0.11), and biologically uninformed SNP set (PREG1: 0.16, FCS: 0.10, and REB: 0.11). Across-breed use of prediction equations still remained a challenge: accuracies by all SNP sets dropped to around zero for all traits. The performance of biologically informed SNP was not significantly better than other sets in Tropical Composites. However, results indicate that biological information obtained from Brahman was successful to predict the fertility traits in Tropical Composite population.
Collapse
Affiliation(s)
- Muhammad S Tahir
- School of Chemistry and Molecular Biosciences, The University of Queensland, St. Lucia Campus, Brisbane 4072, QLD, Australia
| | - Laercio R Porto-Neto
- Commonwealth Scientific and Industrial Research Organization, St. Lucia, Brisbane 4072, QLD, Australia
| | - Toni Reverter-Gomez
- Commonwealth Scientific and Industrial Research Organization, St. Lucia, Brisbane 4072, QLD, Australia
| | - Babatunde S Olasege
- School of Chemistry and Molecular Biosciences, The University of Queensland, St. Lucia Campus, Brisbane 4072, QLD, Australia
| | - Mirza R Sajid
- Department of Statistics, University of Gujrat, 50700 Punjab, Pakistan
| | - Kimberley B Wockner
- Queensland Department of Agriculture and Fisheries, Brisbane 4072, QLD, Australia
| | - Andre W L Tan
- School of Chemistry and Molecular Biosciences, The University of Queensland, St. Lucia Campus, Brisbane 4072, QLD, Australia
| | - Marina R S Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, St. Lucia Campus, Brisbane 4072, QLD, Australia
| |
Collapse
|
8
|
Ribeiro G, Baldi F, Cesar ASM, Alexandre PA, Peripolli E, Ferraz JBS, Fukumasu H. Detection of potential functional variants based on systems-biology: the case of feed efficiency in beef cattle. BMC Genomics 2022; 23:774. [PMID: 36434498 PMCID: PMC9700932 DOI: 10.1186/s12864-022-08958-y] [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: 08/25/2021] [Accepted: 10/20/2022] [Indexed: 11/26/2022] Open
Abstract
BACKGROUND Potential functional variants (PFVs) can be defined as genetic variants responsible for a given phenotype. Ultimately, these are the best DNA markers for animal breeding and selection, especially for polygenic and complex phenotypes. Herein, we described the identification of PFVs for complex phenotypes (in this case, Feed Efficiency in beef cattle) using a systems-biology driven approach based on RNA-seq data from physiologically relevant organs. RESULTS The systems-biology coupled with deep molecular phenotyping by RNA-seq of liver, muscle, hypothalamus, pituitary, and adrenal glands of animals with high and low feed efficiency (FE) measured by residual feed intake (RFI) identified 2,000,936 uniquely variants. Among them, 9986 variants were significantly associated with FE and only 78 had a high impact on protein expression and were considered as PFVs. A set of 169 significant uniquely variants were expressed in all five organs, however, only 27 variants had a moderate impact and none of them a had high impact on protein expression. These results provide evidence of tissue-specific effects of high-impact PFVs. The PFVs were enriched (FDR < 0.05) for processing and presentation of MHC Class I and II mediated antigens, which are an important part of the adaptive immune response. The experimental validation of these PFVs was demonstrated by the increased prediction accuracy for RFI using the weighted G matrix (ssGBLUP+wG; Acc = 0.10 and b = 0.48) obtained in the ssGWAS in comparison to the unweighted G matrix (ssGBLUP; Acc = 0.29 and b = 1.10). CONCLUSION Here we identified PFVs for FE in beef cattle using a strategy based on systems-biology and deep molecular phenotyping. This approach has great potential to be used in genetic prediction programs, especially for polygenic phenotypes.
Collapse
Affiliation(s)
- Gabriela Ribeiro
- grid.11899.380000 0004 1937 0722Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, Pirassununga, Sao Paulo, 13635-900 Brazil
| | - Fernando Baldi
- grid.410543.70000 0001 2188 478XDepartment of Animal Science, São Paulo State University (UNESP), Jaboticabal, São Paulo, Brazil
| | - Aline S. M. Cesar
- grid.11899.380000 0004 1937 0722Escola Superior de Agricultura “Luiz de Queiroz”, University of Sao Paulo, Piracicaba, São Paulo, Brazil
| | - Pâmela A. Alexandre
- grid.11899.380000 0004 1937 0722Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, Pirassununga, Sao Paulo, 13635-900 Brazil ,CSIRO Agriculture & Food, 306 Carmody Rd., St. Lucia, Brisbane, QLD 4067 Australia
| | - Elisa Peripolli
- grid.11899.380000 0004 1937 0722Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, Pirassununga, Sao Paulo, 13635-900 Brazil ,grid.410543.70000 0001 2188 478XDepartment of Animal Science, São Paulo State University (UNESP), Jaboticabal, São Paulo, Brazil
| | - José B. S. Ferraz
- grid.11899.380000 0004 1937 0722Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, Pirassununga, Sao Paulo, 13635-900 Brazil
| | - Heidge Fukumasu
- grid.11899.380000 0004 1937 0722Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, Pirassununga, Sao Paulo, 13635-900 Brazil
| |
Collapse
|
9
|
|
10
|
Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers. Mamm Genome 2022; 33:629-641. [PMID: 35840822 DOI: 10.1007/s00335-022-09960-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2021] [Accepted: 06/28/2022] [Indexed: 10/17/2022]
Abstract
Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.
Collapse
|
11
|
Sepulveda BJ, Muir SK, Bolormaa S, Knight MI, Behrendt R, MacLeod IM, Pryce JE, Daetwyler HD. Eating Time as a Genetic Indicator of Methane Emissions and Feed Efficiency in Australian Maternal Composite Sheep. Front Genet 2022; 13:883520. [PMID: 35646089 PMCID: PMC9130857 DOI: 10.3389/fgene.2022.883520] [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: 02/25/2022] [Accepted: 04/25/2022] [Indexed: 11/13/2022] Open
Abstract
Previous studies have shown reduced enteric methane emissions (ME) and residual feed intake (RFI) through the application of genomic selection in ruminants. The objective of this study was to evaluate feeding behaviour traits as genetic indicators for ME and RFI in Australian Maternal Composite ewes using data from an automated feed intake facility. The feeding behaviour traits evaluated were the amount of time spent eating per day (eating time; ETD; min/day) and per visit (eating time per event; ETE; min/event), daily number of events (DNE), event feed intake (EFI; g/event) and eating rate (ER; g/min). Genotypes and phenotypes of 445 ewes at three different ages (post-weaning, hogget, and adult) were used to estimate the heritability of ME, RFI, and the feeding behaviour traits using univariate genomic best linear unbiased prediction models. Multivariate models were used to estimate the correlations between these traits and within each trait at different ages. The response to selection was evaluated for ME and RFI with direct selection models and indirect models with ETE as an indicator trait, as this behaviour trait was a promising indicator based on heritability and genetic correlations. Heritabilities were between 0.12 and 0.18 for ME and RFI, and between 0.29 and 0.47 for the eating behaviour traits. In our data, selecting for more efficient animals (low RFI) would lead to higher methane emissions per day and per kg of dry matter intake. Selecting for more ETE also improves feed efficiency but results in more methane per day and per kg dry matter intake. Based on our results, ETE could be evaluated as an indicator trait for ME and RFI under an index approach that allows simultaneous selection for improvement in emissions and feed efficiency. Selecting for ETE may have a tremendous impact on the industry, as it may be easier and cheaper to obtain than feed intake and ME data. As the data were collected using individual feeding units, the findings on this research should be validated under grazing conditions.
Collapse
Affiliation(s)
- Boris J Sepulveda
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia
| | | | - Sunduimijid Bolormaa
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC, Australia
| | | | - Ralph Behrendt
- Agriculture Victoria, Hamilton Centre, Hamilton, VIC, Australia
| | - Iona M MacLeod
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC, Australia
| | - Jennie E Pryce
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia
| | - Hans D Daetwyler
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia
| |
Collapse
|
12
|
Rodriguez Neira JD, Peripolli E, de Negreiros MPM, Espigolan R, López-Correa R, Aguilar I, Lobo RB, Baldi F. Prediction ability for growth and maternal traits using SNP arrays based on different marker densities in Nellore cattle using the ssGBLUP. J Appl Genet 2022; 63:389-400. [PMID: 35133621 DOI: 10.1007/s13353-022-00685-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Revised: 01/25/2022] [Accepted: 02/02/2022] [Indexed: 11/25/2022]
Abstract
This study aimed to investigate the prediction ability for growth and maternal traits using different low-density customized SNP arrays selected by informativeness and distribution of markers across the genome employing single-step genomic BLUP (ssGBLUP). Phenotypic records for adjusted weight at 210 and 450 days of age were utilized. A total of 945 animals were genotyped with high-density chip, and 267 individuals born after 2008 were selected as validation population. We evaluated 11 scenarios using five customized density arrays (40 k, 20 k, 10 k, 5 k and 2 k) and the HD array was used as desirable scenario. The GEBV predictions and BIF (Beef Improvement Federation) accuracy were obtained with BLUPF90 family programs. Linear regression was used to evaluate the prediction ability, inflation, and bias of GEBV of each customized array. An overestimation of partial GEBVs in contrast with complete GEBVs and increase of BIF accuracy with the density arrays diminished were observed. For all traits, the prediction ability was higher as the array density increased and it was similar with customized arrays higher than 10 k SNPs. Level of inflation was lower as the density array increased of and was higher for MW210 effect. The bias was susceptible to overestimation of GEBVs when the density customized arrays decreased. These results revealed that the BIF accuracy is sensible to overestimation using low-density customized arrays while the prediction ability with least 10,000 informative SNPs obtained from the Illumina BovineHD BeadChip shows accurate and less biased predictions. Low-density customized arrays under ssGBLUP method could be feasible and cost-effective in genomic selection.
Collapse
Affiliation(s)
- Juan Diego Rodriguez Neira
- Departamento de Zootecnia, Faculdade de Ciências Agrarias e Veterinárias, Universidade Estadual Paulista (Unesp), Jaboticabal, 14884-900, Brazil.
| | - Elisa Peripolli
- Departamento de Zootecnia, Faculdade de Ciências Agrarias e Veterinárias, Universidade Estadual Paulista (Unesp), Jaboticabal, 14884-900, Brazil
| | - Maria Paula Marinho de Negreiros
- Departamento de Medicina Veterinária, Faculdade de Zootecnia e Engenharia de Alimentos, Universidade de São Paulo (Usp), Pirassununga, 13535-900, Brazil
| | - Rafael Espigolan
- Departamento de Medicina Veterinária, Faculdade de Zootecnia e Engenharia de Alimentos, Universidade de São Paulo (Usp), Pirassununga, 13535-900, Brazil
| | - Rodrigo López-Correa
- Departamento de Genética y Mejoramiento Animal, Facultad de Veterinaria, Universidad de La República, Montevideo, Uruguay
| | - Ignacio Aguilar
- Instituto Nacional de Investigación Agropecuaria (INIA), Montevideo, Uruguay
| | - Raysildo B Lobo
- Associação Nacional de Criadores e Pesquisadores (ANCP), Ribeirão Preto, Brazil
| | - Fernando Baldi
- Departamento de Zootecnia, Faculdade de Ciências Agrarias e Veterinárias, Universidade Estadual Paulista (Unesp), Jaboticabal, 14884-900, Brazil
| |
Collapse
|
13
|
Naserkheil M, Mehrban H, Lee D, Park MN. Evaluation of Genome-Enabled Prediction for Carcass Primal Cut Yields Using Single-Step Genomic Best Linear Unbiased Prediction in Hanwoo Cattle. Genes (Basel) 2021; 12:genes12121886. [PMID: 34946834 PMCID: PMC8701981 DOI: 10.3390/genes12121886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 11/16/2021] [Accepted: 11/24/2021] [Indexed: 11/16/2022] Open
Abstract
There is a growing interest worldwide in genetically selecting high-value cut carcass weights, which allows for increased profitability in the beef cattle industry. Primal cut yields have been proposed as a potential indicator of cutability and overall carcass merit, and it is worthwhile to assess the prediction accuracies of genomic selection for these traits. This study was performed to compare the prediction accuracy obtained from a conventional pedigree-based BLUP (PBLUP) and a single-step genomic BLUP (ssGBLUP) method for 10 primal cut traits-bottom round, brisket, chuck, flank, rib, shank, sirloin, striploin, tenderloin, and top round-in Hanwoo cattle with the estimators of the linear regression method. The dataset comprised 3467 phenotypic observations for the studied traits and 3745 genotyped individuals with 43,987 single-nucleotide polymorphisms. In the partial dataset, the accuracies ranged from 0.22 to 0.30 and from 0.37 to 0.54 as evaluated using the PBLUP and ssGBLUP models, respectively. The accuracies of PBLUP and ssGBLUP with the whole dataset varied from 0.45 to 0.75 (average 0.62) and from 0.52 to 0.83 (average 0.71), respectively. The results demonstrate that ssGBLUP performed better than PBLUP averaged over the 10 traits, in terms of prediction accuracy, regardless of considering a partial or whole dataset. Moreover, ssGBLUP generally showed less biased prediction and a value of dispersion closer to 1 than PBLUP across the studied traits. Thus, the ssGBLUP seems to be more suitable for improving the accuracy of predictions for primal cut yields, which can be considered a starting point in future genomic evaluation for these traits in Hanwoo breeding practice.
Collapse
Affiliation(s)
- Masoumeh Naserkheil
- Animal Breeding and Genetics Division, National Institute of Animal Science, Cheonan-si 31000, Chungcheongnam-do, Korea;
| | - Hossein Mehrban
- Department of Animal Science, Shahrekord University, Shahrekord 88186-34141, Iran;
| | - Deukmin Lee
- Department of Animal Life and Environment Sciences, Hankyong National University, Jungang-ro 327, Anseong-si 17579, Gyeonggi-do, Korea
- Correspondence: (D.L.); (M.N.P.); Tel.: +82-31-670-5091 (D.L.); +82-41-580-3355 (M.N.P.)
| | - Mi Na Park
- Animal Breeding and Genetics Division, National Institute of Animal Science, Cheonan-si 31000, Chungcheongnam-do, Korea;
- Correspondence: (D.L.); (M.N.P.); Tel.: +82-31-670-5091 (D.L.); +82-41-580-3355 (M.N.P.)
| |
Collapse
|
14
|
Buaban S, Lengnudum K, Boonkum W, Phakdeedindan P. Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model. J Dairy Sci 2021; 105:468-494. [PMID: 34756438 DOI: 10.3168/jds.2020-19826] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2020] [Accepted: 08/24/2021] [Indexed: 12/26/2022]
Abstract
Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.
Collapse
Affiliation(s)
- S Buaban
- Bureau of Animal Husbandry and Genetic Improvement, Department of Livestock Development, Pathum Thani 12000, Thailand
| | - K Lengnudum
- Bureau of Biotechnology in Livestock Production, Department of Livestock Development, Pathum Thani 12000, Thailand
| | - W Boonkum
- Department of Animal Science, Faculty of Agriculture, Khon Kaen University, Khon Kaen 40002, Thailand
| | - P Phakdeedindan
- Department of Animal Husbandry, Faculty of Veterinary Science, Chulalongkorn University, Bangkok 10330, Thailand; Genomics and Precision Dentistry Research Unit, Department of Physiology, Faculty of Dentistry, Chulalongkorn University, Bangkok 10330, Thailand.
| |
Collapse
|
15
|
Zhu S, Guo T, Yuan C, Liu J, Li J, Han M, Zhao H, Wu Y, Sun W, Wang X, Wang T, Liu J, Tiambo CK, Yue Y, Yang B. Evaluation of Bayesian alphabet and GBLUP based on different marker density for genomic prediction in Alpine Merino sheep. G3 (BETHESDA, MD.) 2021; 11:6310012. [PMID: 34849779 PMCID: PMC8527494 DOI: 10.1093/g3journal/jkab206] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 06/01/2021] [Indexed: 01/20/2023]
Abstract
The marker density, the heritability level of trait and the statistical models adopted are critical to the accuracy of genomic prediction (GP) or selection (GS). If the potential of GP is to be fully utilized to optimize the effect of breeding and selection, in addition to incorporating the above factors into simulated data for analysis, it is essential to incorporate these factors into real data for understanding their impact on GP accuracy, more clearly and intuitively. Herein, we studied the GP of six wool traits of sheep by two different models, including Bayesian Alphabet (BayesA, BayesB, BayesCπ, and Bayesian LASSO) and genomic best linear unbiased prediction (GBLUP). We adopted fivefold cross-validation to perform the accuracy evaluation based on the genotyping data of Alpine Merino sheep (n = 821). The main aim was to study the influence and interaction of different models and marker densities on GP accuracy. The GP accuracy of the six traits was found to be between 0.28 and 0.60, as demonstrated by the cross-validation results. We showed that the accuracy of GP could be improved by increasing the marker density, which is closely related to the model adopted and the heritability level of the trait. Moreover, based on two different marker densities, it was derived that the prediction effect of GBLUP model for traits with low heritability was better; while with the increase of heritability level, the advantage of Bayesian Alphabet would be more obvious, therefore, different models of GP are appropriate in different traits. These findings indicated the significance of applying appropriate models for GP which would assist in further exploring the optimization of GP.
Collapse
Affiliation(s)
- Shaohua Zhu
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Tingting Guo
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Chao Yuan
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Jianbin Liu
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Jianye Li
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Mei Han
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Hongchang Zhao
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Yi Wu
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Weibo Sun
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.,Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Xijun Wang
- Gansu Provincial Sheep Breeding Technology Extension Station, Sunan 734400, China
| | - Tianxiang Wang
- Gansu Provincial Sheep Breeding Technology Extension Station, Sunan 734400, China
| | - Jigang Liu
- Gansu Provincial Sheep Breeding Technology Extension Station, Sunan 734400, China
| | - Christian Keambou Tiambo
- Centre for Tropical Livestock Genetics and Health (CTLGH), International Livestock Research Institute, Nairobi 00100, Kenya
| | - Yaojing Yue
- Sheep Breeding Engineering Technology Center, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Bohui Yang
- Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| |
Collapse
|
16
|
Buaban S, Prempree S, Sumreddee P, Duangjinda M, Masuda Y. Genomic prediction of milk-production traits and somatic cell score using single-step genomic best linear unbiased predictor with random regression test-day model in Thai dairy cattle. J Dairy Sci 2021; 104:12713-12723. [PMID: 34538484 DOI: 10.3168/jds.2021-20263] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2021] [Accepted: 08/04/2021] [Indexed: 12/15/2022]
Abstract
Cow genotypes are expected to improve the accuracy of genomic estimated breeding values (GEBV) for young bulls in relatively small populations such as Thai Holstein-Friesian crossbred dairy cattle in Thailand. The objective of this study was to investigate the effect of cow genotypes on the predictive ability and individual accuracies of GEBV for young dairy bulls in Thailand. Test-day data included milk yield (n = 170,666), milk component traits (fat yield, protein yield, total solids yield, fat percentage, protein percentage, and total solids percentage; n = 160,526), and somatic cell score (n = 82,378) from 23,201, 82,378, and 13,737 (for milk yield, milk component traits, and SCS, respectively) cows calving between 1993 and 2017, respectively. Pedigree information included 51,128; 48,834; and 32,743 animals for milk yield, milk component traits, and somatic cell score, respectively. Additionally, 876, 868, and 632 pedigreed animals (for milk yield, milk component traits, and SCS, respectively) were genotyped (152 bulls and 724 cows), respectively, using Illumina Bovine SNP50 BeadChip. We cut off the data in the last 6 yr, and the validation animals were defined as genotyped bulls with no daughters in the truncated set. We calculated GEBV using a single-step random regression test-day model (SS-RR-TDM), in comparison with estimated breed value (EBV) based on the pedigree-based model used as the official method in Thailand (RR-TDM). Individual accuracies of GEBV were obtained by inverting the coefficient matrix of the mixed model equations, whereas validation accuracies were measured by the Pearson correlation between deregressed EBV from the full data set and (G)EBV predicted with the reduced data set. When only bull genotypes were used, on average, SS-RR-TDM increased individual accuracies by 0.22 and validation accuracies by 0.07, compared with RR-TDM. With cow genotypes, the additional increase was 0.02 for individual accuracies and 0.06 for validation accuracies. The inflation of GEBV tended to be reduced using cow genotypes. Genomic evaluation by SS-RR-TDM is feasible to select young bulls for the longitudinal traits in Thai dairy cattle, and the accuracy of selection is expected to be increased with more genotypes. Genomic selection using the SS-RR-TDM should be implemented in the routine genetic evaluation of the Thai dairy cattle population. The genetic evaluation should consider including genotypes of both sires and cows.
Collapse
Affiliation(s)
- S Buaban
- The Bureau of Biotechnology in Livestock Production, Department of Livestock Development, Pathum Thani 12000, Thailand
| | - S Prempree
- The Bureau of Biotechnology in Livestock Production, Department of Livestock Development, Pathum Thani 12000, Thailand
| | - P Sumreddee
- The Bureau of Biotechnology in Livestock Production, Department of Livestock Development, Pathum Thani 12000, Thailand
| | - M Duangjinda
- Department of Animal Science, Khon Kaen University, Meaung, Khon Kaen 40002, Thailand.
| | - Y Masuda
- Department of Animal and Dairy Science, University of Georgia, Athens 30602
| |
Collapse
|
17
|
Tonussi RL, Londoño-Gil M, de Oliveira Silva RM, Magalhães AFB, Amorim ST, Kluska S, Espigolan R, Peripolli E, Pereira ASC, Lôbo RB, Aguilar I, Lourenço DAL, Baldi F. Accuracy of genomic breeding values and predictive ability for postweaning liveweight and age at first calving in a Nellore cattle population with missing sire information. Trop Anim Health Prod 2021; 53:432. [PMID: 34373940 DOI: 10.1007/s11250-021-02879-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Accepted: 07/30/2021] [Indexed: 11/30/2022]
Abstract
The multiple sire system (MSS) is a common mating scheme in extensive beef production systems. However, MSS does not allow paternity identification and lead to inaccurate genetic predictions. The objective of this study was to investigate the implementation of single-step genomic BLUP (ssGBLUP) in different scenarios of uncertain paternity in the evaluation for 450-day adjusted liveweight (W450) and age at first calving (AFC) in a Nellore cattle population. To estimate the variance components using BLUP and ssGBLUP, the relationship matrix (A) with different proportions of animals with missing sires (MS) (scenarios 0, 25, 50, 75, and 100% of MS) was created. The genotyped animals with MS were randomly chosen, and ten replicates were performed for each scenario and trait. Five groups of animals were evaluated in each scenario: PHE, all animals with phenotypic records in the population; SIR, proven sires; GEN, genotyped animals; YNG, young animals without phenotypes and progeny; and YNGEN, young genotyped animals. The additive genetic variance decreased for both traits as the proportion of MS increased in the population when using the regular REML. When using the ssGBLUP, accuracies ranged from 0.13 to 0.47 for W450 and from 0.10 to 0.25 for AFC. For both traits, the prediction ability of the direct genomic value (DGV) decreased as the percentage of MS increased. These results emphasize that indirect prediction via DGV of young animals is more accurate when the SNP effects are derived from ssGBLUP with a reference population with known sires. The ssGBLUP could be applied in situations of uncertain paternity, especially when selecting young animals. This methodology is shown to be accurate, mainly in scenarios with a high percentage of MS.
Collapse
Affiliation(s)
- Rafael Lara Tonussi
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | - Marisol Londoño-Gil
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil.
| | | | - Ana Fabrícia Braga Magalhães
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | - Sabrina Thaise Amorim
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | - Sabrina Kluska
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | - Rafael Espigolan
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | - Elisa Peripolli
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| | | | - Raysildo Barbosa Lôbo
- Associação Nacional de Criadores E Pesquisadores (ANCP), Ribeirão Preto, SP, CEP 14020-230, Brazil
| | - Ignácio Aguilar
- Instituto Nacional de Pesquisa Agropecuária (INIA), CEP 90200, Las Brujas, Uruguay
| | | | - Fernando Baldi
- Grupo de Melhoramento Animal, Faculdade de Ciências Agrárias E Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, CEP 14884-900, Brazil
| |
Collapse
|
18
|
Salek Ardestani S, Jafarikia M, Sargolzaei M, Sullivan B, Miar Y. Genomic Prediction of Average Daily Gain, Back-Fat Thickness, and Loin Muscle Depth Using Different Genomic Tools in Canadian Swine Populations. Front Genet 2021; 12:665344. [PMID: 34149806 PMCID: PMC8209496 DOI: 10.3389/fgene.2021.665344] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2021] [Accepted: 04/15/2021] [Indexed: 12/12/2022] Open
Abstract
Improvement of prediction accuracy of estimated breeding values (EBVs) can lead to increased profitability for swine breeding companies. This study was performed to compare the accuracy of different popular genomic prediction methods and traditional best linear unbiased prediction (BLUP) for future performance of back-fat thickness (BFT), average daily gain (ADG), and loin muscle depth (LMD) in Canadian Duroc, Landrace, and Yorkshire swine breeds. In this study, 17,019 pigs were genotyped using Illumina 60K and Affymetrix 50K panels. After quality control and imputation steps, a total of 41,304, 48,580, and 49,102 single-nucleotide polymorphisms remained for Duroc (n = 6,649), Landrace (n = 5,362), and Yorkshire (n = 5,008) breeds, respectively. The breeding values of animals in the validation groups (n = 392–774) were predicted before performance test using BLUP, BayesC, BayesCπ, genomic BLUP (GBLUP), and single-step GBLUP (ssGBLUP) methods. The prediction accuracies were obtained using the correlation between the predicted breeding values and their deregressed EBVs (dEBVs) after performance test. The genomic prediction methods showed higher prediction accuracies than traditional BLUP for all scenarios. Although the accuracies of genomic prediction methods were not significantly (P > 0.05) different, ssGBLUP was the most accurate method for Duroc-ADG, Duroc-LMD, Landrace-BFT, Landrace-ADG, and Yorkshire-BFT scenarios, and BayesCπ was the most accurate method for Duroc-BFT, Landrace-LMD, and Yorkshire-ADG scenarios. Furthermore, BayesCπ method was the least biased method for Duroc-LMD, Landrace-BFT, Landrace-ADG, Yorkshire-BFT, and Yorkshire-ADG scenarios. Our findings can be beneficial for accelerating the genetic progress of BFT, ADG, and LMD in Canadian swine populations by selecting more accurate and unbiased genomic prediction methods.
Collapse
Affiliation(s)
| | - Mohsen Jafarikia
- Canadian Centre for Swine Improvement, Ottawa, ON, Canada.,Centre for Genetic Improvement of Livestock (CGIL), Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, ON, Canada.,Select Sires Inc., Plain City, OH, United States
| | - Brian Sullivan
- Canadian Centre for Swine Improvement, Ottawa, ON, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| |
Collapse
|
19
|
Brunes LC, Baldi F, Lopes FB, Lobo RB, Espigolan R, Costa MFO, Magnabosco CU. Selection criteria for feed efficiency-related traits and their association with growth, reproductive and carcass traits in Nelore cattle. ANIMAL PRODUCTION SCIENCE 2021. [DOI: 10.1071/an20487] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Context
Livestock feed costs have a higher impact on the profitability of beef production systems and are directly related to feed efficiency. However, these traits are hard and have high costs to measure, reducing the availability of phenotypic records and reliability of genetic evaluations. Thus, the use of genomic information can increase the robustness of genetic studies that address them.
Aims
The aim of the present study was to estimate genetic parameters for feed efficiency, growth, reproductive and carcass traits in Nelore cattle and the correlated response among them, using genomic information.
Methods
Residual feed intake (RFI), dry-matter intake, feed conversion ratio, feed efficiency (FE), residual average daily gain (RG), residual feed intake and average daily gain (RIG), birthweight, weight at 120, 240, 365 and 450 days of age, scrotal circumference at 365 and 450 days of age, rib-eye area, backfat thickness and rump fat thickness were evaluated. The genetic parameters were estimated using the single-step genomic best linear unbiased prediction approach.
Key results
The FE-related traits showed low to moderate heritability ranging from 0.07 to 0.23. Feed efficiency-related traits showed low genetic correlations with reproductive (–0.24 to 0.27), carcass (–0.17 to 0.27) and growth (–0.19 to 0.24) traits, except for growth with dry-matter intake (0.32–0.56) and weight at 365 days of age with FE (–0.40).
Conclusions
The selection to improve growth, reproductive and carcass traits would not change RFI, RG and RIG. The choice of the most adequate selection criterion depends on the production system, that is, RFI might be used for low-input beef cattle systems, and RIG would be used for more intensive and without-any-dietary-restrictions beef cattle systems.
Implications
The estimates of heritability and genetic correlations suggest that genetic selection for feed efficiency using RFI, RG and RIG in Nellore cattle leads to higher genetic gain than does that using FE and feed conversion ratio without affecting other profitability traits.
Collapse
|
20
|
Brunes LC, Baldi F, Lopes FB, Narciso MG, Lobo RB, Espigolan R, Costa MFO, Magnabosco CU. Genomic prediction ability for feed efficiency traits using different models and pseudo-phenotypes under several validation strategies in Nelore cattle. Animal 2020; 15:100085. [PMID: 33573965 DOI: 10.1016/j.animal.2020.100085] [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: 03/27/2020] [Revised: 09/09/2020] [Accepted: 09/15/2020] [Indexed: 10/22/2022] Open
Abstract
There is a growing interest to improve feed efficiency (FE) traits in cattle. The genomic selection was proposed to improve these traits since they are difficult and expensive to measure. Up to date, there are scarce studies about the implementation of genomic selection for FE traits in indicine cattle under different scenarios of pseudo-phenotypes, models, and validation strategies on a commercial large scale. Thus, the aim was to evaluate the feasibility of genomic selection implementation for FE traits in Nelore cattle applying different models and pseudo-phenotypes under validation strategies. Phenotypic and genotypic information from 4 329 and 3 467 animals were used, respectively, which were tested for residual feed intake, DM intake, feed efficiency, feed conversion ratio, residual BW gain, and residual intake and BW gain. Six prediction methods were used: single-step genomic best linear unbiased prediction, Bayes A, Bayes B, Bayes Cπ, Bayesian least absolute shrinkage and selection operator (BLASSO), and Bayes R. Phenotypes adjusted for fixed effects (Y*), estimated breeding value (EBV), and EBV deregressed (DEBV) were used as pseudo-phenotypes. The validation approaches used were: (1) random: the data was randomly divided into ten subsets and the validation was done in each subset at a time; (2) age: the partition into training and testing sets was based on year of birth and testing animals were born after 2016; and (3) EBV accuracy: the data was split into two groups, being animals with accuracy above 0.45 the training set; and below 0.45 the validation set. In the analyses that used the Y* as pseudo-phenotype, prediction ability (PA) was obtained by dividing the correlation between pseudo-phenotype and genomic EBV (GEBV) by the square root of the heritability of the trait. When EBV and DEBV were used as the pseudo-phenotype, the simple correlation of this quantity with the GEBV was considered as PA. The prediction methods show similar results for PA and bias. The random cross-validation presented higher PA (0.17) than EBV accuracy (0.14) and age (0.13). The PA was higher for Y* than for EBV and DEBV (30.0 and 34.3%, respectively). Random validation presented the highest PA, being indicated for use in populations composed mainly of young animals and traits with few generations of data recording. For high heritability traits, the validation can be done by age, enabling the prediction of the next-generation genetic merit. These results would support breeders to identify genomic approaches that are more viable for genomic prediction for FE-related traits.
Collapse
Affiliation(s)
- L C Brunes
- Animal Science Department, Goiás Federal University, 74690-900 Goiânia, GO, Brazil; Embrapa Rice and Beans, GO-462, km 12, 75375-000 Santo Antônio de Goiás, GO, Brazil.
| | - F Baldi
- Animal Science Department, São Paulo State University - Júlio de Mesquita Filho (UNESP), Prof. Paulo Donato Castelane, 14884-900 Jaboticabal, SP, Brazil
| | - F B Lopes
- Cobb-Vantress, Inc., 72761 Siloam Springs, AR, USA
| | - M G Narciso
- Embrapa Rice and Beans, GO-462, km 12, 75375-000 Santo Antônio de Goiás, GO, Brazil
| | - R B Lobo
- National Association of Breeders and Researchers, 14020-230 Ribeirão Preto, Brazil
| | - R Espigolan
- Department of Veterinary Medicine, Faculty of Animal Science and Food Engineering, University of Sao Paulo, 13635-900 Pirassununga, SP, Brazil
| | - M F O Costa
- Embrapa Rice and Beans, GO-462, km 12, 75375-000 Santo Antônio de Goiás, GO, Brazil
| | - C U Magnabosco
- Embrapa Cerrados, BR-020, 18 Sobradinho, 70770-901 Brasilia, DF, Brazil
| |
Collapse
|
21
|
Banos G, Lindsay V, Desta TT, Bettridge J, Sanchez-Molano E, Vallejo-Trujillo A, Matika O, Dessie T, Wigley P, Christley RM, Kaiser P, Hanotte O, Psifidi A. Integrating Genetic and Genomic Analyses of Combined Health Data Across Ecotypes to Improve Disease Resistance in Indigenous African Chickens. Front Genet 2020; 11:543890. [PMID: 33193617 PMCID: PMC7581896 DOI: 10.3389/fgene.2020.543890] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Accepted: 09/04/2020] [Indexed: 12/30/2022] Open
Abstract
Poultry play an important role in the agriculture of many African countries. The majority of chickens in sub-Saharan Africa are indigenous, raised in villages under semi-scavenging conditions. Vaccinations and biosecurity measures rarely apply, and infectious diseases remain a major cause of mortality and reduced productivity. Genomic selection for disease resistance offers a potentially sustainable solution but this requires sufficient numbers of individual birds with genomic and phenotypic data, which is often a challenge to collect in the small populations of indigenous chicken ecotypes. The use of information across-ecotypes presents an attractive possibility to increase the relevant numbers and the accuracy of genomic selection. In this study, we performed a joint analysis of two distinct Ethiopian indigenous chicken ecotypes to investigate the genomic architecture of important health and productivity traits and explore the feasibility of conducting genomic selection across-ecotype. Phenotypic traits considered were antibody response to Infectious Bursal Disease (IBDV), Marek's Disease (MDV), Fowl Cholera (PM) and Fowl Typhoid (SG), resistance to Eimeria and cestode parasitism, and productivity [body weight and body condition score (BCS)]. Combined data from the two chicken ecotypes, Horro (n = 384) and Jarso (n = 376), were jointly analyzed for genetic parameter estimation, genome-wide association studies (GWAS), genomic breeding value (GEBVs) calculation, genomic predictions, whole-genome sequencing (WGS), and pathways analyses. Estimates of across-ecotype heritability were significant and moderate in magnitude (0.22-0.47) for all traits except for SG and BCS. GWAS identified several significant genomic associations with health and productivity traits. The WGS analysis revealed putative candidate genes and mutations for IBDV (TOLLIP, ANGPTL5, BCL9, THEMIS2), MDV (GRM7), SG (MAP3K21), Eimeria (TOM1L1) and cestodes (TNFAIP1, ATG9A, NOS2) parasitism, which warrant further investigation. Reliability of GEBVs increased compared to within-ecotype calculations but accuracy of genomic prediction did not, probably because the genetic distance between the two ecotypes offset the benefit from increased sample size. However, for some traits genomic prediction was only feasible in across-ecotype analysis. Our results generally underpin the potential of genomic selection to enhance health and productivity across-ecotypes. Future studies should establish the required minimum sample size and genetic similarity between ecotypes to ensure accurate joint genomic selection.
Collapse
Affiliation(s)
- Georgios Banos
- The Roslin Institute, The University of Edinburgh, Edinburgh, United Kingdom
- Scotland’s Rural College, Edinburgh, United Kingdom
- Centre for Tropical Livestock Genetics and Health, Edinburgh, United Kingdom
| | - Victoria Lindsay
- Royal Veterinary College, University of London, London, United Kingdom
| | - Takele T. Desta
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Judy Bettridge
- Institute of Infection and Global Health, University of Liverpool, Liverpool, United Kingdom
- LiveGene – Centre for Tropical Livestock Genetics and Health, International Livestock Research Institute, Addis Ababa, Ethiopia
- Natural Resources Institute, University of Greenwich, London, United Kingdom
| | | | | | - Oswald Matika
- The Roslin Institute, The University of Edinburgh, Edinburgh, United Kingdom
| | - Tadelle Dessie
- LiveGene – Centre for Tropical Livestock Genetics and Health, International Livestock Research Institute, Addis Ababa, Ethiopia
| | - Paul Wigley
- Institute of Infection and Global Health, University of Liverpool, Liverpool, United Kingdom
| | - Robert M. Christley
- Institute of Infection and Global Health, University of Liverpool, Liverpool, United Kingdom
| | - Peter Kaiser
- The Roslin Institute, The University of Edinburgh, Edinburgh, United Kingdom
| | - Olivier Hanotte
- Centre for Tropical Livestock Genetics and Health, Edinburgh, United Kingdom
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
- LiveGene – Centre for Tropical Livestock Genetics and Health, International Livestock Research Institute, Addis Ababa, Ethiopia
| | - Androniki Psifidi
- The Roslin Institute, The University of Edinburgh, Edinburgh, United Kingdom
- Centre for Tropical Livestock Genetics and Health, Edinburgh, United Kingdom
- Royal Veterinary College, University of London, London, United Kingdom
| |
Collapse
|
22
|
Stolpovsky YA, Piskunov AK, Svishcheva GR. Genomic Selection. I: Latest Trends and Possible Ways of Development. RUSS J GENET+ 2020. [DOI: 10.1134/s1022795420090148] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
23
|
Vieira Ventura R, Fonseca E Silva F, Manuel Yáñez J, Brito LF. Opportunities and challenges of phenomics applied to livestock and aquaculture breeding in South America. Anim Front 2020; 10:45-52. [PMID: 32368412 PMCID: PMC7189274 DOI: 10.1093/af/vfaa008] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Ricardo Vieira Ventura
- Department of Animal Nutrition and Production, Faculty of Veterinary Medicine and Animal Science, University of São Paulo (FMVZ/USP), Pirassununga, SP, Brazil
| | | | - José Manuel Yáñez
- Faculty of Veterinary and Animal Sciences, University of Chile, Santa Rosa, La Pintana, Santiago, Chile
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN
| |
Collapse
|
24
|
Ramos PVB, E Silva FF, da Silva LOC, Santiago GG, Menezes GRDO, Soriano Viana JM, Torres Júnior RAA, Gondo A, Brito LF. Genomic evaluation for novel stayability traits in Nellore cattle. Reprod Domest Anim 2020; 55:266-273. [PMID: 31880841 DOI: 10.1111/rda.13612] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Accepted: 12/19/2019] [Indexed: 02/01/2023]
Abstract
Cow stayability plays a major role on the overall profitability of the beef cattle industry, as it is directly related to reproductive efficiency and cow's longevity. Stayability (STAY63) is usually defined as the ability of the cow to calve at least three times until 76 months of age. This is a late-measured and lowly heritable trait, which consequently constrains genetic progress per time unit. Thus, the use of genomic information associated with novel stayability traits measured earlier in life will likely result in higher prediction accuracy and faster genetic progress for cow longevity. In this study, we aimed to compare pedigree-based and single-step GBLUP (ssGBLUP) methods as well as to estimate genetic correlations between the proposed stayability traits: STAY42, STAY53 and STAY64, which are measured at 52, 64 and 76 months of cow's age, considering at least 2, 3 and 4 calving, respectively. ssGBLUP yielded the highest prediction accuracy for all traits. The heritability estimates for STAY42, STAY53, STAY63 and STAY64 were 0.090, 0.151, 0.152 and 0.143, respectively. The genetic correlations between traits ranged from 0.899 (STAY42 and STAY53) to 0.985 (STAY53 and STAY63). The high genetic correlation between STAY42 and STAY53 suggests that besides being related to cow longevity, STAY53 is also associated with the early-stage reproductive efficiency. Thus, STAY53 is recommended as a suitable selection criterion for reproductive efficiency due to its higher heritability, favourable genetic correlation with other traits, and measured earlier in life, compared with the conventional stayability trait, that is STAY63.
Collapse
Affiliation(s)
| | | | | | - Gustavo Garcia Santiago
- Faculty of Veterinary Medicine and Animal Science, Federal University of Mato Grosso do Sul, Campo Grande, Brazil
| | | | | | | | | | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, USA
| |
Collapse
|
25
|
Feitosa FLB, Pereira ASC, Amorim ST, Peripolli E, Silva RMDO, Braz CU, Ferrinho AM, Schenkel FS, Brito LF, Espigolan R, de Albuquerque LG, Baldi F. Comparison between haplotype-based and individual snp-based genomic predictions for beef fatty acid profile in Nelore cattle. J Anim Breed Genet 2019; 137:468-476. [PMID: 31867831 DOI: 10.1111/jbg.12463] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 11/18/2019] [Accepted: 11/24/2019] [Indexed: 11/28/2022]
Abstract
The aim of this study was to evaluate the genomic predictions using the single-step genomic best linear unbiased predictor (ssGBLUP) method based on SNPs and haplotype markers associated with beef fatty acids (FAs) profile in Nelore cattle. The data set contained records from 963 Nelore bulls finished in feedlot (±90 days) and slaughtered with approximately 24 months of age. Meat samples from the Longissimus dorsi muscle were taken for FAs profile measurement. FAs were quantified by gas chromatography using a SP-2560 capillary column. Animals were genotyped with the high-density SNP panel (BovineHD BeadChip assay) containing 777,962 markers. SNPs with a minor allele frequency and a call rate lower than 0.05 and 0.90, respectively, monomorphic, located on sex chromosomes, and with unknown position were removed from the data set. After genomic quality control, a total of 469,981 SNPs and 892 samples were available for subsequent analyses. Missing genotypes were imputed and phased using the FImpute software. Haplotype blocks were defined based on linkage disequilibrium using the Haploview software. The model to estimate variance components and genetic parameters and to predict the genomic values included the random genetic additive effects, fixed effects of the contemporary group and the age at slaughter as a linear covariate. Accuracies using the haplotype-based approach ranged from 0.07 to 0.31, and those SNP-based ranged from 0.06 to 0.33. Regression coefficients ranged from 0.07 to 0.74 and from 0.08 to 1.45 using the haplotype- and SNP-based approaches, respectively. Despite the low to moderate accuracies for the genomic values, it is possible to obtain genetic progress trough selection using genomic information based either on SNPs or haplotype markers. The SNP-based approach allows less biased genomic evaluations, and it is more feasible when taking into account the computational and operational cost underlying the haplotypes inference.
Collapse
Affiliation(s)
- Fabieli Loise Braga Feitosa
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| | - Angélica Simone Cravo Pereira
- Faculdade de Zootecnia e Engenharia de Alimentos, Departamento de Nutrição e Produção Animal, Universidade de São Paulo, Pirassununga, Brazil
| | - Sabrina Thaise Amorim
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| | - Elisa Peripolli
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| | | | - Camila Urbano Braz
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| | - Adrielle Matias Ferrinho
- Faculdade de Zootecnia e Engenharia de Alimentos, Departamento de Nutrição e Produção Animal, Universidade de São Paulo, Pirassununga, Brazil
| | | | | | - Rafael Espigolan
- Faculdade de Zootecnia e Engenharia de Alimentos, Departamento de Medicina Veterinária, Universidade de São Paulo, Pirassununga, Brazil
| | - Lucia Galvão de Albuquerque
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| | - Fernando Baldi
- Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Universidade Estadual Paulista Júlio de Mesquita Filho, UNESP, Jaboticabal, Brazil
| |
Collapse
|
26
|
de las Heras-Saldana S, Clark SA, Duijvesteijn N, Gondro C, van der Werf JHJ, Chen Y. Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle. BMC Genomics 2019; 20:939. [PMID: 31810463 PMCID: PMC6898931 DOI: 10.1186/s12864-019-6270-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Accepted: 11/07/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions. RESULTS The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15. CONCLUSIONS Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.
Collapse
Affiliation(s)
| | - Samuel A. Clark
- School of Environmental and Rural Science, University of New England, Armidale, NSW Australia
| | - Naomi Duijvesteijn
- School of Environmental and Rural Science, University of New England, Armidale, NSW Australia
| | - Cedric Gondro
- School of Environmental and Rural Science, University of New England, Armidale, NSW Australia
- Department of Animal Science, Michigan State University, East Lansing, MI USA
| | | | - Yizhou Chen
- Department of Primary Industries, Elizabeth Macarthur Agricultural Institute, Menangle, NSW Australia
| |
Collapse
|
27
|
Sollero BP, Howard JT, Spangler ML. The impact of reducing the frequency of animals genotyped at higher density on imputation and prediction accuracies using ssGBLUP1. J Anim Sci 2019; 97:2780-2792. [PMID: 31115442 DOI: 10.1093/jas/skz147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Accepted: 04/25/2019] [Indexed: 11/12/2022] Open
Abstract
The largest gains in accuracy in a genomic selection program come from genotyping young selection candidates who have not yet produced progeny and who might, or might not, have a phenotypic record recorded. To reduce genotyping costs and to allow for an increased amount of genomic data to be available in a population, young selection candidates may be genotyped with low-density (LD) panels and imputed to a higher density. However, to ensure that a reasonable imputation accuracy persists overtime, some parent animals originally genotyped at LD must be re-genotyped at a higher density. This study investigated the long-term impact of selectively re-genotyping parents with a medium-density (MD) SNP panel on the accuracy of imputation and on the genetic predictions using ssGBLUP in a simulated beef cattle population. Assuming a moderately heritable trait (0.25) and a population undergoing selection, the simulation generated sequence data for a founder population (100 male and 500 female individuals) and 9,000 neutral markers, considered as the MD panel. All selection candidates from generation 8 to 15 were genotyped with LD panels corresponding to a density of 0.5% (LD_0.5), 2% (LD_2), and 5% (LD_5) of the MD. Re-genotyping scenarios chose parents at random or based on EBV and ranged from 10% of male parents to re-genotyping all male and female parents with MD. Ranges in average imputation accuracy at generation 15 were 0.567 to 0.936, 0.795 to 0.985, and 0.931 to 0.995 for the LD_0.5, LD_2, and LD_5, respectively, and the average EBV accuracies ranged from 0.453 to 0.735, 0.631 to 0.784, and 0.748 to 0.807 for LD_0.5, LD_2, and LD_5, respectively. Re-genotyping parents based on their EBV resulted in higher imputation and EBV accuracies compared to selecting parents at random and these values increased with the size of LD panels. Differences between re-genotyping scenarios decreased when the density of the LD panel increased, suggesting fewer animals needed to be re-genotyped to achieve higher accuracies. In general, imputation and EBV accuracies were greater when more parents were re-genotyped, independent of the proportion of males and females. In practice, the relationship between the density of the LD panel used and the target panel must be considered to determine the number (proportion) of animals that would need to be re-genotyped to enable sufficient imputation accuracy.
Collapse
|
28
|
Karimi K, Sargolzaei M, Plastow GS, Wang Z, Miar Y. Opportunities for genomic selection in American mink: A simulation study. PLoS One 2019; 14:e0213873. [PMID: 30870528 PMCID: PMC6417779 DOI: 10.1371/journal.pone.0213873] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Accepted: 03/01/2019] [Indexed: 12/25/2022] Open
Abstract
Genomic selection can be considered as an effective tool for developing breeding programs in American mink. However, the genetic gains for economically important traits can be influenced by the accuracy of genomic predictions. The objective of this study was to investigate the prediction accuracies of traditional best linear unbiased prediction (BLUP), multi-step genomic BLUP (GBLUP) and single-step GBLUP (ssGBLUP) methods in American mink using simulated data with different levels of heritability, marker density, training set (TS) sizes and selection designs based on either phenotypic performance or estimated breeding values (EBVs). Under EBV selection design, the accuracy of BLUP predictions was increased by 38% and 44% for h2 = 0.10, 27% and 29% for h2 = 0.20, and 5.8% and 6% for h2 = 0.50 using GBLUP and ssGBLUP methods, respectively. Under phenotypic selection design, the accuracies of prediction by ssGBLUP method were 11.8% and 15.4% higher than those obtained by GBLUP for heritability of 0.10 and 0.20, respectively. However, the efficiency of ssGBLUP and GBLUP was not influenced by selection design at higher level of heritability (h2 = 0.50). Furthermore, higher selection intensity increased the bias of predictions in both pedigree-based and genomic evaluations. Regardless of selection design, TS sizes for GBLUP and ssGBLUP methods should be at least 3000 to achieve more accuracy than using BLUP for heritability of 0.50 and marker density of 10k and 50k. Overall, more accurate predictions were obtained using ssGBLUP method particularly for lowly heritable traits and low density of markers. Our results indicated that TS sizes should be optimized in accordance with heritability level, marker density, selection design and prediction method for genomic selection in American mink. The results provided an initial framework for designing genomic selection in mink breeding programs.
Collapse
Affiliation(s)
- Karim Karimi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, Ontario, Canada
- Select Sires Inc., Plain City, Ohio, United States of America
| | - Graham Stuart Plastow
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
| | - Zhiquan Wang
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| |
Collapse
|
29
|
Genomic selection for meat quality traits in Nelore cattle. Meat Sci 2019; 148:32-37. [DOI: 10.1016/j.meatsci.2018.09.010] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Revised: 09/11/2018] [Accepted: 09/17/2018] [Indexed: 11/24/2022]
|
30
|
Mrode R, Ojango JMK, Okeyo AM, Mwacharo JM. Genomic Selection and Use of Molecular Tools in Breeding Programs for Indigenous and Crossbred Cattle in Developing Countries: Current Status and Future Prospects. Front Genet 2019; 9:694. [PMID: 30687382 PMCID: PMC6334160 DOI: 10.3389/fgene.2018.00694] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2018] [Accepted: 12/11/2018] [Indexed: 11/23/2022] Open
Abstract
Genomic selection (GS) has resulted in rapid rates of genetic gains especially in dairy cattle in developed countries resulting in a higher proportion of genomically proven young bulls being used in breeding. This success has been undergirded by well-established conventional genetic evaluation systems. Here, the status of GS in terms of the structure of the reference and validation populations, response variables, genomic prediction models, validation methods, and imputation efficiency in breeding programs of developing countries, where smallholder systems predominate and the basic components for conventional breeding are mostly lacking is examined. Also, the application of genomic tools and identification of genome-wide signatures of selection is reviewed. The studies on genomic prediction in developing countries are mostly in dairy and beef cattle usually with small reference populations (500-3,000 animals) and are mostly cows. The input variables tended to be pre-corrected phenotypic records and the small reference populations has made implementation of various Bayesian methods feasible in addition to GBLUP. Multi-trait single-step has been used to incorporate genomic information from foreign bulls, thus GS in developing countries would benefit from collaborations with developed countries, as many dairy sires used are from developed countries where they may have been genotyped and phenotyped. Cross validation approaches have been implemented in most studies resulting in accuracies of 0.20-0.60. Genotyping animals with a mixture of HD and LD chips, followed by imputation to the HD have been implemented with imputation accuracies of 0.74-0.99 reported. This increases the prospects of reducing genotyping costs and hence the cost-effectiveness of GS. Next-generation sequencing and associated technologies have allowed the determination of breed composition, parent verification, genome diversity, and genome-wide selection sweeps. This information can be incorporated into breeding programs aiming to utilize GS. Cost-effective GS in beef cattle in developing countries may involve usage of reproductive technologies (AI and in-vitro fertilization) to efficiently propagate superior genetics from the genomics pipeline. For dairy cattle, sexed semen of genomically proven young bulls could substantially improve profitability thus increase prospects of small holder farmers buying-in into genomic breeding programs.
Collapse
Affiliation(s)
- Raphael Mrode
- Animal Biosciences, International Livestock Research Institute, Nairobi, Kenya
- Animal and Veterinary Science, Scotland Rural College, Edinburgh, United Kingdom
| | - Julie M. K Ojango
- Animal Biosciences, International Livestock Research Institute, Nairobi, Kenya
| | - A. M. Okeyo
- Animal Biosciences, International Livestock Research Institute, Nairobi, Kenya
| | - Joram M. Mwacharo
- Small Ruminant Genomics, International Centre for Agricultural Research in the Dry Areas (ICARDA), Addis Ababa, Ethiopia
| |
Collapse
|
31
|
Genomic prediction ability for beef fatty acid profile in Nelore cattle using different pseudo-phenotypes. J Appl Genet 2018; 59:493-501. [DOI: 10.1007/s13353-018-0470-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2018] [Revised: 08/28/2018] [Accepted: 09/17/2018] [Indexed: 11/26/2022]
|
32
|
Genetic correlations and heritability estimates for dry matter intake, weight gain and feed efficiency of Nellore cattle in feedlot. Livest Sci 2018. [DOI: 10.1016/j.livsci.2018.06.013] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
33
|
Hay EH, Roberts A. Genome-wide association study for carcass traits in a composite beef cattle breed. Livest Sci 2018. [DOI: 10.1016/j.livsci.2018.04.018] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
34
|
Genomic predictions combining SNP markers and copy number variations in Nellore cattle. BMC Genomics 2018; 19:441. [PMID: 29871610 PMCID: PMC5989480 DOI: 10.1186/s12864-018-4787-6] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Accepted: 05/14/2018] [Indexed: 11/26/2022] Open
Abstract
Background Due to the advancement in high throughput technology, single nucleotide polymorphism (SNP) is routinely being incorporated along with phenotypic information into genetic evaluation. However, this approach often cannot achieve high accuracy for some complex traits. It is possible that SNP markers are not sufficient to predict these traits due to the missing heritability caused by other genetic variations such as microsatellite and copy number variation (CNV), which have been shown to affect disease and complex traits in humans and other species. Results In this study, CNVs were included in a SNP based genomic selection framework. A Nellore cattle dataset consisting of 2230 animals genotyped on BovineHD SNP array was used, and 9 weight and carcass traits were analyzed. A total of six models were implemented and compared based on their prediction accuracy. For comparison, three models including only SNPs were implemented: 1) BayesA model, 2) Bayesian mixture model (BayesB), and 3) a GBLUP model without polygenic effects. The other three models incorporating both SNP and CNV included 4) a Bayesian model similar to BayesA (BayesA+CNV), 5) a Bayesian mixture model (BayesB+CNV), and 6) GBLUP with CNVs modeled as a covariable (GBLUP+CNV). Prediction accuracies were assessed based on Pearson’s correlation between de-regressed EBVs (dEBVs) and direct genomic values (DGVs) in the validation dataset. For BayesA, BayesB and GBLUP, accuracy ranged from 0.12 to 0.62 across the nine traits. A minimal increase in prediction accuracy for some traits was noticed when including CNVs in the model (BayesA+CNV, BayesB+CNV, GBLUP+CNV). Conclusions This study presents the first genomic prediction study integrating CNVs and SNPs in livestock. Combining CNV and SNP marker information proved to be beneficial for genomic prediction of some traits in Nellore cattle. Electronic supplementary material The online version of this article (10.1186/s12864-018-4787-6) contains supplementary material, which is available to authorized users.
Collapse
|
35
|
Cardoso DF, de Albuquerque LG, Reimer C, Qanbari S, Erbe M, do Nascimento AV, Venturini GC, Scalez DCB, Baldi F, de Camargo GMF, Mercadante MEZ, do Santos Gonçalves Cyrillo JN, Simianer H, Tonhati H. Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle. Genet Sel Evol 2018; 50:22. [PMID: 29720080 PMCID: PMC5930444 DOI: 10.1186/s12711-018-0381-2] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Accepted: 02/20/2018] [Indexed: 12/11/2022] Open
Abstract
Background This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection. Results Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (FST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4. Conclusions The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight. Electronic supplementary material The online version of this article (10.1186/s12711-018-0381-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Diercles F Cardoso
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.
| | - Lucia Galvão de Albuquerque
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.,National Counsel of Technological and Scientific Development (CNPq), Brasília, DF, Brazil
| | - Christian Reimer
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Goettingen, Germany
| | - Saber Qanbari
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Goettingen, Germany
| | - Malena Erbe
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Goettingen, Germany.,Institute for Animal Breeding, Bavarian State Research Center for Agriculture, Grub, Germany
| | - André V do Nascimento
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | - Guilherme C Venturini
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | - Daiane C Becker Scalez
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | - Fernando Baldi
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.,National Counsel of Technological and Scientific Development (CNPq), Brasília, DF, Brazil
| | - Gregório M Ferreira de Camargo
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | - Maria E Zerlotti Mercadante
- National Counsel of Technological and Scientific Development (CNPq), Brasília, DF, Brazil.,APTA Beef Cattle Center, Institute of Animal Science, Sertãozinho, SP, Brazil
| | | | - Henner Simianer
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Goettingen, Germany
| | - Humberto Tonhati
- Department of Animal Science, Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.,National Counsel of Technological and Scientific Development (CNPq), Brasília, DF, Brazil
| |
Collapse
|
36
|
Zhang C, Kemp RA, Stothard P, Wang Z, Boddicker N, Krivushin K, Dekkers J, Plastow G. Genomic evaluation of feed efficiency component traits in Duroc pigs using 80K, 650K and whole-genome sequence variants. Genet Sel Evol 2018; 50:14. [PMID: 29625549 PMCID: PMC5889553 DOI: 10.1186/s12711-018-0387-9] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2017] [Accepted: 03/27/2018] [Indexed: 12/30/2022] Open
Abstract
BACKGROUND Increasing marker density was proposed to have potential to improve the accuracy of genomic prediction for quantitative traits; whole-sequence data is expected to give the best accuracy of prediction, since all causal mutations that underlie a trait are expected to be included. However, in cattle and chicken, this assumption is not supported by empirical studies. Our objective was to compare the accuracy of genomic prediction of feed efficiency component traits in Duroc pigs using single nucleotide polymorphism (SNP) panels of 80K, imputed 650K, and whole-genome sequence variants using GBLUP, BayesB and BayesRC methods, with the ultimate purpose to determine the optimal method to increase genetic gain for feed efficiency in pigs. RESULTS Phenotypes of average daily feed intake (ADFI), average daily gain (ADG), ultrasound backfat depth (FAT), and loin muscle depth (LMD) were available for 1363 Duroc boars from a commercial breeding program. Genotype imputation accuracies reached 92.1% from 80K to 650K and 85.6% from 650K to whole-genome sequence variants. Average accuracies across methods and marker densities of genomic prediction of ADFI, FAT, LMD and ADG were 0.40, 0.65, 0.30 and 0.15, respectively. For ADFI and FAT, BayesB outperformed GBLUP, but increasing marker density had little advantage for genomic prediction. For ADG and LMD, GBLUP outperformed BayesB, while BayesRC based on whole-genome sequence data gave the best accuracies and reached up to 0.35 for LMD and 0.25 for ADG. CONCLUSIONS Use of genomic information was beneficial for prediction of ADFI and FAT but not for that of ADG and LMD compared to pedigree-based estimates. BayesB based on 80K SNPs gave the best genomic prediction accuracy for ADFI and FAT, while BayesRC based on whole-genome sequence data performed best for ADG and LMD. We suggest that these differences between traits in the effect of marker density and method on accuracy of genomic prediction are mainly due to the underlying genetic architecture of the traits.
Collapse
Affiliation(s)
- Chunyan Zhang
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, T6G 2R3, Canada
| | | | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, T6G 2R3, Canada
| | - Zhiquan Wang
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, T6G 2R3, Canada
| | | | - Kirill Krivushin
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, T6G 2R3, Canada
| | - Jack Dekkers
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Graham Plastow
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, T6G 2R3, Canada.
| |
Collapse
|
37
|
Yoshida GM, Carvalheiro R, Rodríguez FH, Lhorente JP, Yáñez JM. Single-step genomic evaluation improves accuracy of breeding value predictions for resistance to infectious pancreatic necrosis virus in rainbow trout. Genomics 2018; 111:127-132. [PMID: 29357303 DOI: 10.1016/j.ygeno.2018.01.008] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Revised: 12/30/2017] [Accepted: 01/14/2018] [Indexed: 11/26/2022]
Abstract
The aim of this study was to compare the accuracy of breeding values (EBVs) predicted using the traditional pedigree based Best Linear Unbiased Prediction (PBLUP) and the single-step genomic Best Linear Unbiased Prediction (ssGBLUP) for resistance against infectious pancreatic necrosis virus (IPNV) in rainbow trout. A total of 2278 animals were challenged against IPNV and 768 individuals were genotyped using a 57 K single nucleotide polymorphism array for rainbow trout. Accuracies for both methods were assessed using five-fold cross-validation. The heritabilities were higher for PBLUP compared to ssGBLUP. The ssGBLUP accuracies outperformed PBLUP in 7 and 11% for days to death and binary survival, respectively. The ssGBLUP could be an alternative approach to improve the accuracy of breeding values for resistance against infectious pancreatic necrosis virus in rainbow trout, using information from genotyped and non-genotyped animals.
Collapse
Affiliation(s)
- Grazyella M Yoshida
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Av Santa Rosa 11735, La Pintana, Santiago 8820808, Chile; Animal Science Department, Universidade Estadual Paulista "Júlio de Mesquita Filho" (UNESP), Faculdade de Ciências Agrárias e Veterinárias (FCAV), Campus Jaboticabal, Via de Acesso Prof. Paulo Donato Castellane, 14884-900 Jaboticabal, Brazil
| | - Roberto Carvalheiro
- Animal Science Department, Universidade Estadual Paulista "Júlio de Mesquita Filho" (UNESP), Faculdade de Ciências Agrárias e Veterinárias (FCAV), Campus Jaboticabal, Via de Acesso Prof. Paulo Donato Castellane, 14884-900 Jaboticabal, Brazil
| | - Francisco H Rodríguez
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Av Santa Rosa 11735, La Pintana, Santiago 8820808, Chile; Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional del Altiplano, Av. Floral 1153, Puno, Perú
| | | | - José M Yáñez
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Av Santa Rosa 11735, La Pintana, Santiago 8820808, Chile; Aquainnovo, Cardonal S/N, Puerto Montt, Chile; Núcleo Milenio INVASAL, Concepción, Chile.
| |
Collapse
|
38
|
Evaluation of the potential use of a meta-population for genomic selection in autochthonous beef cattle populations. Animal 2017; 12:1350-1357. [PMID: 29094666 DOI: 10.1017/s175173111700283x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
This study investigated the potential application of genomic selection under a multi-breed scheme in the Spanish autochthonous beef cattle populations using a simulation study that replicates the structure of linkage disequilibrium obtained from a sample of 25 triplets of sire/dam/offspring per population and using the BovineHD Beadchip. Purebred and combined reference sets were used for the genomic evaluation and several scenarios of different genetic architecture of the trait were investigated. The single-breed evaluations yielded the highest within-breed accuracies. Across breed accuracies were found low but positive on average confirming the genetic connectedness between the populations. If the same genotyping effort is split in several populations, the accuracies were lower when compared with single-breed evaluation, but showed a small advantage over small-sized purebred reference sets over the accuracies of subsequent generations. Besides, the genetic architecture of the trait did not show any relevant effect on the accuracy with the exception of rare variants, which yielded slightly lower results and higher loss of predictive ability over the generations.
Collapse
|
39
|
Guo P, Zhu B, Xu L, Niu H, Wang Z, Guan L, Liang Y, Ni H, Guo Y, Chen Y, Zhang L, Gao X, Gao H, Li J. Genomic prediction with parallel computing for slaughter traits in Chinese Simmental beef cattle using high-density genotypes. PLoS One 2017; 12:e0179885. [PMID: 28723906 PMCID: PMC5516975 DOI: 10.1371/journal.pone.0179885] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2017] [Accepted: 06/06/2017] [Indexed: 11/18/2022] Open
Abstract
Genomic selection has been widely used for complex quantitative trait in farm animals. Estimations of breeding values for slaughter traits are most important to beef cattle industry, and it is worthwhile to investigate prediction accuracies of genomic selection for these traits. In this study, we assessed genomic predictive abilities for average daily gain weight (ADG), live weight (LW), carcass weight (CW), dressing percentage (DP), lean meat percentage (LMP) and retail meat weight (RMW) using Illumina Bovine 770K SNP Beadchip in Chinese Simmental cattle. To evaluate the abilities of prediction, marker effects were estimated using genomic BLUP (GBLUP) and three parallel Bayesian models, including multiple chains parallel BayesA, BayesB and BayesCπ (PBayesA, PBayesB and PBayesCπ). Training set and validation set were divided by random allocation, and the predictive accuracies were evaluated using 5-fold cross validations. We found the accuracies of genomic predictions ranged from 0.195±0.084 (GBLUP for LMP) to 0.424±0.147 (PBayesB for CW). The average accuracies across traits were 0.327±0.085 (GBLUP), 0.335±0.063 (PBayesA), 0.347±0.093 (PBayesB) and 0.334±0.077 (PBayesCπ), respectively. Notably, parallel Bayesian models were more accurate than GBLUP across six traits. Our study suggested that genomic selections with multiple chains parallel Bayesian models are feasible for slaughter traits in Chinese Simmental cattle. The estimations of direct genomic breeding values using parallel Bayesian methods can offer important insights into improving prediction accuracy at young ages and may also help to identify superior candidates in breeding programs.
Collapse
Affiliation(s)
- Peng Guo
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- College of Computer and Information Engineering, Tianjin Agricultural University, Tianjin, China
| | - Bo Zhu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Lingyang Xu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- * E-mail: (LYX); (JYL)
| | - Hong Niu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Zezhao Wang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Long Guan
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Yonghu Liang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Hemin Ni
- Animal Science and Technology College, Beijing University of Agriculture, Beijing, China
| | - Yong Guo
- Animal Science and Technology College, Beijing University of Agriculture, Beijing, China
| | - Yan Chen
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Lupei Zhang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Huijiang Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- * E-mail: (LYX); (JYL)
| |
Collapse
|
40
|
Chiaia HLJ, Peripoli E, Silva RMDO, Aboujaoude C, Feitosa FLB, Lemos MVAD, Berton MP, Olivieri BF, Espigolan R, Tonussi RL, Gordo DGM, Bresolin T, Magalhães AFB, Júnior GAF, Albuquerque LGD, Oliveira HND, Furlan JDJM, Ferrinho AM, Mueller LF, Tonhati H, Pereira ASC, Baldi F. Genomic prediction for beef fatty acid profile in Nellore cattle. Meat Sci 2017; 128:60-67. [PMID: 28214693 DOI: 10.1016/j.meatsci.2017.02.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2016] [Revised: 01/14/2017] [Accepted: 02/07/2017] [Indexed: 12/25/2022]
Abstract
The objective of this study was to compare SNP-BLUP, BayesCπ, BayesC and Bayesian Lasso methodologies to predict the direct genomic value for saturated, monounsaturated, and polyunsaturated fatty acid profile, omega 3 and 6 in the Longissimus thoracis muscle of Nellore cattle finished in feedlot. A total of 963 Nellore bulls with phenotype for fatty acid profiles, were genotyped using the Illumina BovineHD BeadChip (Illumina, San Diego, CA) with 777,962 SNP. The predictive ability was evaluated using cross validation. To compare the methodologies, the correlation between DGV and pseudo-phenotypes was calculated. The accuracy varied from -0.40 to 0.62. Our results indicate that none of the methods excelled in terms of accuracy, however, the SNP-BLUP method allows obtaining less biased genomic evaluations, thereby; this method is more feasible when taking into account the analyses' operating cost. Despite the lowest bias observed for EBV, the adjusted phenotype is the preferred pseudophenotype considering the genomic prediction accuracies regarding the context of the present study.
Collapse
Affiliation(s)
| | - Elisa Peripoli
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | - Carolyn Aboujaoude
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | | | - Mariana Piatto Berton
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | - Rafael Espigolan
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | - Rafael Lara Tonussi
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | - Tiago Bresolin
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | | | | | | | | | | | - Lenise Freitas Mueller
- Faculdade de Zootecnia e Engenharia de Alimentos, USP, Pirassununga, SP, 13635-900, Brazil
| | - Humberto Tonhati
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| | | | - Fernando Baldi
- Faculdade de Ciências Agrárias e Veterinárias, UNESP, Jaboticabal, SP, 14884-000, Brazil
| |
Collapse
|
41
|
Olivieri BF, Mercadante MEZ, Cyrillo JNDSG, Branco RH, Bonilha SFM, de Albuquerque LG, Silva RMDO, Baldi F. Genomic Regions Associated with Feed Efficiency Indicator Traits in an Experimental Nellore Cattle Population. PLoS One 2016; 11:e0164390. [PMID: 27760167 PMCID: PMC5070821 DOI: 10.1371/journal.pone.0164390] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Accepted: 09/23/2016] [Indexed: 01/18/2023] Open
Abstract
The objective of this study was to identify genomic regions and metabolic pathways associated with dry matter intake, average daily gain, feed efficiency and residual feed intake in an experimental Nellore cattle population. The high-density SNP chip (Illumina High-Density Bovine BeadChip, 777k) was used to genotype the animals. The SNP markers effects and their variances were estimated using the single-step genome wide association method. The (co)variance components were estimated by Bayesian inference. The chromosome segments that are responsible for more than 1.0% of additive genetic variance were selected to explore and determine possible quantitative trait loci. The bovine genome Map Viewer was used to identify genes. In total, 51 genomic regions were identified for all analyzed traits. The heritability estimated for feed efficiency was low magnitude (0.13±0.06). For average daily gain, dry matter intake and residual feed intake, heritability was moderate to high (0.43±0.05; 0.47±0.05, 0.18±0.05, respectively). A total of 8, 17, 14 and 12 windows that are responsible for more than 1% of the additive genetic variance for dry matter intake, average daily gain, feed efficiency and residual feed intake, respectively, were identified. Candidate genes GOLIM4, RFX6, CACNG7, CACNG6, CAPN8, CAPN2, AKT2, GPRC6A, and GPR45 were associated with feed efficiency traits. It was expected that the response to selection would be higher for residual feed intake than for feed efficiency. Genomic regions harboring possible QTL for feed efficiency indicator traits were identified. Candidate genes identified are involved in energy use, metabolism protein, ion transport, transmembrane transport, the olfactory system, the immune system, secretion and cellular activity. The identification of these regions and their respective candidate genes should contribute to the formation of a genetic basis in Nellore cattle for feed efficiency indicator traits, and these results would support the selection for these traits.
Collapse
Affiliation(s)
- Bianca Ferreira Olivieri
- Universidade Estadual Paulista, Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Via de acesso Prof. Paulo Donato Castellane, s/no, CEP 14884-900 Jaboticabal, SP, Brazil
| | - Maria Eugênia Zerlotti Mercadante
- Instituto de Zootecnia, Centro Avançado de Pesquisa Tecnológica do Agronegócio de Bovinos de Corte, Rodovia Carlos Tonanni, km 94, CEP 14.174-000, Sertãozinho, SP, Brazil
| | | | - Renata Helena Branco
- Instituto de Zootecnia, Centro Avançado de Pesquisa Tecnológica do Agronegócio de Bovinos de Corte, Rodovia Carlos Tonanni, km 94, CEP 14.174-000, Sertãozinho, SP, Brazil
| | - Sarah Figueiredo Martins Bonilha
- Instituto de Zootecnia, Centro Avançado de Pesquisa Tecnológica do Agronegócio de Bovinos de Corte, Rodovia Carlos Tonanni, km 94, CEP 14.174-000, Sertãozinho, SP, Brazil
| | - Lucia Galvão de Albuquerque
- Universidade Estadual Paulista, Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Via de acesso Prof. Paulo Donato Castellane, s/no, CEP 14884-900 Jaboticabal, SP, Brazil
| | - Rafael Medeiros de Oliveira Silva
- Universidade Estadual Paulista, Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Via de acesso Prof. Paulo Donato Castellane, s/no, CEP 14884-900 Jaboticabal, SP, Brazil
| | - Fernando Baldi
- Universidade Estadual Paulista, Faculdade de Ciências Agrárias e Veterinárias, Departamento de Zootecnia, Via de acesso Prof. Paulo Donato Castellane, s/no, CEP 14884-900 Jaboticabal, SP, Brazil
| |
Collapse
|