Accuracy of genomic selection for age at puberty in a multi-breed population of tropically adapted beef cattle

A zero altered Poisson random forest model for genomic-enabled prediction

In genomic selection choosing the statistical machine learning model is of paramount importance. In this paper, we present an application of a zero altered random forest model with two versions (ZAP_RF and ZAPC_RF) to deal with excess zeros in count response variables. The proposed model was compared with the conventional random forest (RF) model and with the conventional Generalized Poisson Ridge regression (GPR) using two real datasets, and we found that, in terms of prediction performance, the proposed zero inflated random forest model outperformed the conventional RF and GPR models.

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

A Multivariate Poisson Deep Learning Model for Genomic Prediction of Count Data

The paradigm called genomic selection (GS) is a revolutionary way of developing new plants and animals. This is a predictive methodology, since it uses learning methods to perform its task. Unfortunately, there is no universal model that can be used for all types of predictions; for this reason, specific methodologies are required for each type of output (response variables). Since there is a lack of efficient methodologies for multivariate count data outcomes, in this paper, a multivariate Poisson deep neural network (MPDN) model is proposed for the genomic prediction of various count outcomes simultaneously. The MPDN model uses the minus log-likelihood of a Poisson distribution as a loss function, in hidden layers for capturing nonlinear patterns using the rectified linear unit (RELU) activation function and, in the output layer, the exponential activation function was used for producing outputs on the same scale of counts. The proposed MPDN model was compared to conventional generalized Poisson regression models and univariate Poisson deep learning models in two experimental data sets of count data. We found that the proposed MPDL outperformed univariate Poisson deep neural network models, but did not outperform, in terms of prediction, the univariate generalized Poisson regression models. All deep learning models were implemented in Tensorflow as back-end and Keras as front-end, which allows implementing these models on moderate and large data sets, which is a significant advantage over previous GS models for multivariate count data.

Enhancing Genetic Gain through Genomic Selection: From Livestock to Plants

Although long-term genetic gain has been achieved through increasing use of modern breeding methods and technologies, the rate of genetic gain needs to be accelerated to meet humanity's demand for agricultural products. In this regard, genomic selection (GS) has been considered most promising for genetic improvement of the complex traits controlled by many genes each with minor effects. Livestock scientists pioneered GS application largely due to livestock's significantly higher individual values and the greater reduction in generation interval that can be achieved in GS. Large-scale application of GS in plants can be achieved by refining field management to improve heritability estimation and prediction accuracy and developing optimum GS models with the consideration of genotype-by-environment interaction and non-additive effects, along with significant cost reduction. Moreover, it would be more effective to integrate GS with other breeding tools and platforms for accelerating the breeding process and thereby further enhancing genetic gain. In addition, establishing an open-source breeding network and developing transdisciplinary approaches would be essential in enhancing breeding efficiency for small- and medium-sized enterprises and agricultural research systems in developing countries. New strategies centered on GS for enhancing genetic gain need to be developed.

Multivariate genomic predictions for age at puberty in tropically adapted beef heifers

Heifers that have an earlier age at puberty often have greater lifetime productivity. Age at puberty is moderately heritable so selection should effectively reduce the number of days to puberty, and improve heifer productivity and profitability as a result. However, recording age at puberty is intensive, requiring repeat ovarian scanning to determine age at first corpus luteum (AGECL). Genomic selection has been proposed as a strategy to select for earlier age at puberty; however, large reference populations of cows with AGECL records and genotypes would be required to generate accurate GEBV for this trait. Reproductive maturity score (RMS) is a proxy trait for age at puberty for implementation in northern Australia beef herds, where large scale recording of AGECL is not feasible. RMS assigns a score of 0 to 5 from a single ovarian scan to describe ovarian maturity at ~600 d. Here we use multivariate genomic prediction to evaluate the value of a large RMS data set to improve accuracy of GEBV for age at puberty (AGECL). There were 882 Brahman and 990 Tropical Composite heifers with AGECL phenotypes, and an independent set of 974 Brahman, 1,798 Santa Gertrudis, and 910 Droughtmaster heifers with RMS phenotypes. All animals had 728,785 real or imputed SNP genotypes. The correlation of AGECL and RMS (h2 = 0.23) was estimated as -0.83 using the genomic information. This result also demonstrates that using genomic information it is possible to estimate genetic correlations between traits collected on different animals in different herds, with minimal or unknown pedigree linkage between them. Inclusion of heifers with RMS in the multi-trait model improved the accuracy of genomic evaluations for AGECL. Accuracy of RMS GEBV generally did not improve by adding heifers with AGECL phenotypes into the reference population. These results suggest that RMS and AGECL may be used together in a multi-trait prediction model to increase the accuracy of prediction for age at puberty in tropically adapted beef cattle.

Towards multi-breed genomic evaluations for female fertility of tropical beef cattle

Developing accurate genomic evaluations of fertility for tropical beef cattle must deal with at least two major challenges (i) recording cow fertility traits in extensive production systems on large numbers of cows and (ii) the genomic evaluations should work across the breeds, crossbreds, and composites used in tropical beef production. Here, we assess accuracy of genomic evaluations for a trait which can be collected on a large scale in extensive conditions, corpus luteum score (CLscore), which is 1 if ovarian scanning indicates a heifer has cycled by 600 d and 0 if not, in a multi-breed population. A total of 3,696 heifers, including 979 Brahmans, 914 Droughtmasters, and 1,803 Santa Gertrudis in seven herds across 3-yr cohorts with CLscores, were genotyped for 24,211 SNPs. Genotypes were imputed to 728,785 SNPs. GBLUP and BayesR were used to predict GEBV. Accuracy of GEBV was evaluated with two validation strategies. In the first strategy, the last year cohort of heifers from each herd was used for validation, such that every herd had heifers in both reference and validation populations. In the second validation strategy, each herd in turn was removed in its entirety from the reference population, and was used for validation. For both validation strategies, accuracy of GEBV for single breed and multi-breed reference populations was assessed. For the first validation strategy, accuracy of GEBV ranged from 0.2 for Brahmans to 0.4 for Droughtmasters. Increasing marker density from 24K SNPs to 728K SNPs resulted in a small increase in accuracy, and including multiple-breeds in the reference did not help improve accuracy. These results suggest that provided a herd has animals in the reference population, the accuracy of the GEBV is largely determined by within herd (linkage) information. The situation was very different when entire herds were predicted in the second validation. In this case accuracy of GEBV using only 24K SNPs and only a within breed reference was close to zero for all breeds. Accuracy increased substantially when 728K SNPs, BayesR, and a multi-breed reference were used, from 0.15 for Brahmans to 0.35 for Santa Gertrudis. Given the second validation strategy is more likely to reflect the situation for many herds in tropical beef production (no animals in the reference), genomic evaluations for fertility in tropical beef cattle should be based on high-density markers (728K SNPs) and should be multi-breed.

Genomic prediction of continuous and binary fertility traits of females in a composite beef cattle breed

Reproduction efficiency is a major factor in the profitability of the beef cattle industry. Genomic selection (GS) is a promising tool that may improve the predictive accuracy and genetic gain of fertility traits. There is a wide range of traits used to measure fertility in dairy and beef cattle including continuous (days open), discrete (pregnancy status), and count (number of inseminations) responses. In this study, a joint analysis of age of puberty (AOP), age at first calving (AOC), and the heifer pregnancy status (HPS) was performed. Data used in this study consisted of records from 1,365 Composite Gene Combination (CGC; 50% Red Angus, 25% Charolais, 25% Tarentaise) first parity females born between 2002 and 2011. The pedigree file included 5,374 animals. A total of 3,902 animals were genotyped with different density SNP chips (3K to 50K SNP). Animals genotyped with low-density arrays were imputed to higher density (BovineSNP50 BeadChip) using FImpute. Data were analyzed using univariate and multivariate classical quantitative models (pedigree based) and univariate genomic approaches. For the latter, 3 different Bayesian methods (BayesA, BayesB, and BayesCπ) were implemented and compared. Estimates of heritabilities using univariate and multivariate analyses based on pedigree relationships ranged between 0.03 (for AOC) to 0.2 (AOP). Heritability of pregnancy status was 0.15 and 0.09 using the univariate and multivariate analyses, respectively. Genetic correlation between pregnancy status and the other 2 traits was low being 0.08 with age at puberty and -0.10 with age at first calving. Heritability estimates were slightly higher using genomic rather than average additive relationships. The accuracy of genomic prediction was similar across the 3 Bayesian methods with higher accuracies for age of puberty than the age at first calving likely due to the higher heritability of the former. The prediction of the binary pregnancy status measured using the area under the curve increased by 27% to 29% compared to a random classifier. Due to the small size of the data, all estimates have large posterior standard deviations and results should be interpreted with caution.

Prediction of genomic breeding values for growth, carcass and meat quality traits in a multi-breed sheep population using a HD SNP chip

Background

New Zealand has some unique Terminal Sire composite sheep breeds, which were developed in the last three decades to meet commercial needs. These composite breeds were developed based on crossing various Terminal Sire and Maternal breeds and, therefore, present high genetic diversity compared to other sheep breeds. Their breeding programs are focused on improving carcass and meat quality traits. There is an interest from the industry to implement genomic selection in this population to increase the rates of genetic gain. Therefore, the main objectives of this study were to determine the accuracy of predicted genomic breeding values for various growth, carcass and meat quality traits using a HD SNP chip and to evaluate alternative genomic relationship matrices, validation designs and genomic prediction scenarios. A large multi-breed population (n = 14,845) was genotyped with the HD SNP chip (600 K) and phenotypes were collected for a variety of traits.

Results

The average observed accuracies (± SD) for traits measured in the live animal, carcass, and, meat quality traits ranged from 0.18 ± 0.07 to 0.33 ± 0.10, 0.28 ± 0.09 to 0.55 ± 0.05 and 0.21 ± 0.07 to 0.36 ± 0.08, respectively, depending on the scenario/method used in the genomic predictions. When accounting for population stratification by adjusting for 2, 4 or 6 principal components (PCs) the observed accuracies of molecular breeding values (mBVs) decreased or kept constant for all traits. The mBVs observed accuracies when fitting both G and A matrices were similar to fitting only G matrix. The lowest accuracies were observed for k-means cross-validation and forward validation performed within each k-means cluster.

Conclusions

The accuracies observed in this study support the feasibility of genomic selection for growth, carcass and meat quality traits in New Zealand Terminal Sire breeds using the Ovine HD SNP chip. There was a clear advantage on using a mixed training population instead of performing analyzes per genomic clusters. In order to perform genomic predictions per breed group, genotyping more animals is recommended to increase the size of the training population within each group and the genetic relationship between training and validation populations. The different scenarios evaluated in this study will help geneticists and breeders to make wiser decisions in their breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0476-8) contains supplementary material, which is available to authorized users.