Evaluation of a multi-line broiler chicken population using a single-step genomic evaluation procedure

Effects on prediction of analysing a multi-line chicken population as one line were evaluated. Body weight records were provided by Cobb-Vantress for two lines of broiler chickens. Phenotypic records for 183 695 and 164 149 broilers and genotypic records for 3195 and 3001 broilers were available for each line. Lines were combined to create a multi-line population and analysed using a single-step procedure combining the additive relationship matrix and the genomic relationship matrix (G). G was scaled using allele frequencies from each line, the multi-line population, or 0.5. When allele frequencies were calculated from each line, distributions of diagonal elements were bimodal. When allele frequencies were calculated from the multi-line population, the distribution of diagonal elements had one peak. When allele frequency 0.5 was used, the distribution was bimodal. Genomic estimated breeding values (GEBVs) were predicted using each allele frequency. GEBVs differed with allele frequency but had ≥ 0.99 correlations with GEBVs predicted with correct allele frequencies. Means of each line and differences in mean between the lines differed based on allele frequencies. Assumed allele frequencies have little impact on ranking within line but larger impact on ranking across lines. G may be used to evaluate multiple populations simultaneously but must be adjusted to obtain properly scaled estimates when population structure is unknown.     

Investigating the persistence of accuracy of genomic predictions over time in broilers

Accuracy of genomic predictions is an important component of the selection response. The objectives of this research were: 1) to investigate trends for prediction accuracies over time in a broiler population of accumulated phenotypes, genotypes, and pedigrees and 2) to test if data from distant generations are useful to maintain prediction accuracies in selection candidates. The data contained 820K phenotypes for a growth trait (GT), 200K for two feed efficiency traits (FE1 and FE2), and 42K for a carcass yield trait (CY). The pedigree included 1,252,619 birds hatched over 7 years, of which 154,318 from the last 4 years were genotyped. Training populations were constructed adding 1 year of data sequentially, persistency of accuracy over time was evaluated using predictions from birds hatched in the three generations following or in the years after the training populations. In the first generation, before genotypes became available for the training populations (first 3 years of data), accuracies remained almost stable with successive additions of phenotypes and pedigree to the accumulated dataset. The inclusion of 1 year of genotypes in addition to 4 years of phenotypes and pedigree in the training population led to increases in accuracy of 54% for GT, 76% for FE1, 110% for CY, and 38% for FE2; on average, 74% of the increase was due to genomics. Prediction accuracies declined faster without than with genomic information in the training populations. When genotypes were unavailable, the average decline in prediction accuracy across traits was 41% from the first to the second generation of validation, and 51% from the second to the third generation of validation. When genotypes were available, the average decline across traits was 14% from the first to the second generation of validation, and 3% from the second to the third generation of validation. Prediction accuracies in the last three generations were the same when the training population included 5 or 2 years of data, and a decrease of ~7% was observed when the training population included only 1 year of data. Training sets including genomic information provided an increase in accuracy and persistence of genomic predictions compared with training sets without genomic data. The two most recent years of pedigree, phenotypic, and genomic data were sufficient to maintain prediction accuracies in selection candidates. Similar conclusions were obtained using validation populations per year.     

Genetic components of heat stress in finishing pigs: development of a heat load function

The objective of this study was to quantify the effect of heat stress during the life of a pig on its final weight, as a first step toward a genetic evaluation for heat tolerance. Data included carcass weights of 23,556 crossbred pigs [Duroc x (Landrace x Large White)] raised on 2 farms in North Carolina and slaughtered from May 2005 through December 2006. Weather data were available from a nearby weather station. Lifetime of a pig was assumed to be partitioned into 2 periods. During an initial period, the effect of heat stress was assumed to be negligible or compensated for later. During the second period ending in slaughtering, the ADG was assumed to be affected linearly by heat load. Weekly heat load was calculated as degrees of average temperature-humidity index in excess of a threshold (18 degrees C). The total heat load (H) was the sum of weekly heat loads during the second period. During the months of January to May H was 0; H reached a peak in September. The final BW during the peak of heat stress decreased about 6 kg compared with BW during months of non-heat stress. Weekly and monthly averages of carcass weight generally moved similarly to H. However, there were large fluctuations unrelated to H; the fluctuations were different on the 2 farms. The model included the effects of farm-year of slaughter, sex, age at slaughter, and H, where age at slaughter and H were linear regressions. In analyses, the threshold was varied from 16 to 20 degrees C, and the second period was varied from 8 to 16 wk. The greatest R(2) (10.4%) was at the threshold of temperature-humidity index = 18 degrees C for a period of 10 wk. Varying the threshold and the length of time reduced R(2) less than 1%. Least squares means of year-month and year-week of carcass weight were calculated using a model with the fixed effects farm-year-month or farm-year-week of slaughter, sex, and age at slaughter (linear covariate), and the random effect of birth litter. Changes in BW of finisher pigs due to heat stress can be quantified by H during the last 10 wk of the life of the pig.     

Evaluation of methods for computing approximate accuracies of predicted breeding values in maternal random regression models for growth traits in beef cattle

The objective of this study was to determine the suitability of 2 methods for computing approximate accuracies of predicted breeding values, in which accuracy was defined as the squared correlation between the predicted and true breeding value, when modeling growth traits in beef cattle using random regression (RR) models. The first method (Strabel et al., S-M-B) was designed for use with multitrait models; thus, its use with RR models requires the clustering of measurements into different traits. The second method (Tier and Meyer, T-M) was more general, because it accounted for random coefficients other than zeros and ones and thus it could be used directly when fitting RR models. To investigate the performance of both methods, their results were compared with the true accuracies using a balanced simulated data set. The largest difference between approximate and true average accuracies for direct effects was observed at 205 d when S-M-B was used (4.6% males and 8.8% females). With regard to maternal effects, the largest differences in average accuracies were observed at 205 d in males when S-M-B was used (31.8%) and at the same age in females but when using T-M (33.3%). In general, bias increased for direct effect accuracies in males at the tails of the accuracy range, but for females and for maternal effect accuracies in both sexes, bias increased as accuracy increased. When a population was simulated to create large numbers of progeny for base females that did not have individual records, much greater errors were observed in the regression of approximate values on the true ones. When both approximate methods were compared using a real beef cattle data set, a good agreement was observed, particularly for direct effect accuracies in sires [i.e., at 205 d, the regressions were 0.98 (direct) and 0.95 (maternal) with r(2) over 0.99]. The largest discrepancies for sires between the methods were observed at 205 d for direct (2.7%) and maternal (16.3%) effect accuracies. For dams, the largest differences between methods were also observed at 205 d, 9.3% (direct), and 15.2% (maternal). The differences between methods for nonparent cattle were greater than for dams for maternal effect accuracies but intermediate between sires and dams for direct effect accuracies. In spite of the less biased results provided by T-M, its use could be problematic when employed in evaluations of large populations due to its greater memory and computation requirements (e.g., 170 and 478% more than S-M-B for a population of 11 million).     

Modeling genetic differences of combined broiler chicken populations in single-step GBLUP

The introduction of animals from a different environment or population is a common practice in commercial livestock populations. In this study, we modeled the inclusion of a group of external birds into a local broiler chicken population for the purpose of genomic evaluations. The pedigree was composed of 242,413 birds and genotypes were available for 107,216 birds. A five-trait model that included one growth, two yield, and two efficiency traits was used for the analyses. The strategies to model the introduction of external birds were to include a fixed effect representing the origin of parents and to use unknown parent groups (UPG) or metafounders (MF). Genomic estimated breeding values (GEBV) were obtained with single-step GBLUP using the Algorithm for Proven and Young. Bias, dispersion, and accuracy of GEBV for the validation birds, that is, from the most recent generation, were computed. The bias and dispersion were estimated with the linear regression (LR) method,whereas accuracy was estimated by the LR method and predictive ability. When fixed UPG were fit without estimated inbreeding, the model did not converge. In contrast, models with fixed UPG and estimated inbreeding or random UPG converged and resulted in similar GEBV. The inclusion of an extra fixed effect in the model made the GEBV unbiased and reduced the inflation. Genomic predictions with MF were slightly biased and inflated due to the unbalanced number of observations assigned to each metafounder. When combining local and external populations, the greatest accuracy can be obtained by adding an extra fixed effect to account for the origin of parents plus UPG with estimated inbreeding or random UPG. To estimate the accuracy, the LR method is more consistent among scenarios, whereas the predictive ability greatly depends on the model specification.     

The Possible Involvement of Salsolinol and Hypothalamic Prolactin in the Central Regulatory Processes in Ewes During Lactation

Salsolinol, a dopamine‐related compound and prolactin‐producing cells were found in the ovine hypothalamus. This study was designed to test the hypothesis that salsolinol, acting from the CNS level, is able to stimulate pituitary prolactin release as well as prolactin mRNA expression in the anterior pituitary cells (AP) and in the mediobasal hypothalamus (MBH) in lactating ewes. The intracerebroventricular infusions of salsolinol in two doses, total of 50 ng or 5 μg, were performed in a series of five 10‐min infusions at 20‐min intervals. All infusions were made from 12:30 to 15:00 and the pre‐infusion period was from 10:00 to 12.30 h. The prolactin concentration in plasma samples, collected every 10 min, was determined by radioimmunoassay; prolactin mRNA expression in AP and MBH tissues was determined by real‐time PCR. The obtained results showed that salsolinol infused at the higher dose significantly (p < 0.001) increased plasma prolactin concentration in lactating ewes, when compared with the concentration noted before the infusion and with that in lactating controls. In lactating ewes, the relative levels of prolactin mRNA expression in the AP and MBH were up to twofold and fivefold higher respectively than in non‐lactating ewes (p < 0.05). In our experimental design, salsolinol did not significantly affect the ongoing process of prolactin gene expression in these tissues. We conclude that in ewes, salsolinol may be involved, at least, in the process of stimulation of prolactin release during lactation and that hypothalamic prolactin plays an important role in the central mechanisms of adaptation to lactation.     

Estimation of genetic parameters for mature weight in Angus cattle

The aim of this study was to estimate genetic parameters for BW of Angus cattle up to 5 yr of age and to discuss options for including mature weight (MW) in their genetic evaluation. Data were obtained from the American Angus Association. Only records from herds with at least 500 animals and with >10% of animals with BW at ≥ 2 yr of age were considered. Traits were weaning weight (WW, n = 81,525), yearling weight (YW, n = 62,721), and BW measured from 2 to 5 yr of age (MW2, n = 15,927; MW3, n = 12,404; MW4, n = 9,805; MW5, n = 7,546). Genetic parameters were estimated using an AIREML algorithm with a multiple-trait animal model. Fixed effects were contemporary group and departure of the actual age from standard age (205, 365, 730, 1,095, 1,460, and 1,825 d of age for WW, YW, MW2, MW3, MW4, and MW5, respectively). Random effects were animal direct additive genetic, maternal additive genetic, maternal permanent environment, and residual. Estimates of direct genetic variances (kg(2)) were 298 ± 71.8, 563 ± 15.1, 925 ± 52.1, 1,221 ± 65.8, 1,406 ± 80.4, and 1,402 ± 66.9; maternal genetic variances were 167 ± 4.8, 153 ± 6.1, 123 ± 9.1, 136 ± 12.25, 167 ± 18.0, and 110 ± 14.0; maternal permanent environment variances were 124 ± 2.9, 120 ± 4.3, 61 ± 7.5, 69 ± 11.9, 103 ± 15.9, and 134 ± 35.2; and residual variances were 258 ± 3.8, 608 ± 8.6, 829 ± 34.2, 1,016 ± 38.8, 1,017 ± 52.1, and 1,202 ± 63.22 for WW, YW, MW2, MW3, MW4, and MW5, respectively. The direct genetic correlation between WW and YW was 0.84 ± 0.14 and between WW and MW ranged from 0.66 ± 0.06 (WW and MW4) to 0.72 ± 0.11 (WW and MW2). Direct genetic correlations ranged from 0.77 ± 0.08 (YW and MW5) to 0.85 ± 0.07 (YW and MW2) between YW and MW, and they were ≥ 0.95 among MW2, MW3, MW4, and MW5. Maternal genetic correlations between WW and YW and MW ranged from 0.52 ± 0.05 (WW and MW4) to 0.95 ± 0.07 (WW and YW), and among MW they ranged from 0.54 ± 0.14 (MW4 and MW5) to 0.94 ± 0.07 (MW2 and MW3). Genetic correlations suggest that a genetic evaluation for MW may be MW2-based and that including BW from older ages could be accomplished by adjusting records to the scale of MW2.     

Is single-step genomic REML with the algorithm for proven and young more computationally efficient when less generations of data are present?

Efficient computing techniques allow the estimation of variance components for virtually any traditional dataset. When genomic information is available, variance components can be estimated using genomic REML (GREML). If only a portion of the animals have genotypes, single-step GREML (ssGREML) is the method of choice. The genomic relationship matrix (G) used in both cases is dense, limiting computations depending on the number of genotyped animals. The algorithm for proven and young (APY) can be used to create a sparse inverse of G (GAPY~-1) with close to linear memory and computing requirements. In ssGREML, the inverse of the realized relationship matrix (H⁻¹) also includes the inverse of the pedigree relationship matrix, which can be dense with a long pedigree, but sparser with short. The main purpose of this study was to investigate whether costs of ssGREML can be reduced using APY with truncated pedigree and phenotypes. We also investigated the impact of truncation on variance components estimation when different numbers of core animals are used in APY. Simulations included 150K animals from 10 generations, with selection. Phenotypes (h² = 0.3) were available for all animals in generations 1–9. A total of 30K animals in generations 8 and 9, and 15K validation animals in generation 10 were genotyped for 52,890 SNP. Average information REML and ssGREML with G⁻¹ and GAPY~-1 using 1K, 5K, 9K, and 14K core animals were compared. Variance components are impacted when the core group in APY represents the number of eigenvalues explaining a small fraction of the total variation in G. The most time-consuming operation was the inversion of G, with more than 50% of the total time. Next, numerical factorization consumed nearly 30% of the total computing time. On average, a 7% decrease in the computing time for ordering was observed by removing each generation of data. APY can be successfully applied to create the inverse of the genomic relationship matrix used in ssGREML for estimating variance components. To ensure reliable variance component estimation, it is important to use a core size that corresponds to the number of largest eigenvalues explaining around 98% of total variation in G. When APY is used, pedigrees can be truncated to increase the sparsity of H and slightly reduce computing time for ordering and symbolic factorization, with no impact on the estimates.