Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Modeling genetic components of hatch of fertile in broiler breeders

Reproductive efficiency such as fertility and hatch of fertile (HoF) are of economic importance and concern to breeding companies becaue of their effects on chick output. Similar to other traits of economic importance in poultry breeding, the rate of response for HoF is largely dependent on the use of an appropriate model for evaluating the trait. Therefore, the objectives of this study were to estimate genetic parameters from cumulative, repeatability, fixed regression, random regression, and multitrait models for HoF from a pure-line broiler breeder. The data available for this study consisted of weekly HoF records from 11,729 hens with a total pedigree record of 38,260. Estimates of heritability from the various models ranged from 0.04 to 0.22 with the highest estimate obtained from the cumulative model and the lowest from the repeatability model. Responses to selection estimated for the different models ranged from 0.03 to 0.08% gain per year of the phenotypic mean. In general, the cumulative and the repeatability models underestimated response to selection. The multitrait and random regression models gave similar results for response to selection at 0.08 percentage change in phenotypic mean. In conclusion, the cumulative model is not optimal for modeling HoF, and likewise, the repeatability model. The random regression and multitrait models should be considered instead as they offered a higher response to selection. However, if a multitrait analysis is to be considered, it is recommended to split up the production period in such a way as to avoid computational constraints due to overparameterization.

Estimation of genetic parameters for test-day milk yield in Girolando cows using a random regression model

ABSTRACT The objective of this study was to estimate the components of variance and genetic parameters of test-day milk yield in first lactation Girolando cows, using a random regression model. A total of 126,892 test-day milk yield (TDMY) records of 15,351 first-parity Holstein, Gyr, and Girolando breed cows were used, obtained from the Associação Brasileira dos Criadores de Girolando. To estimate the components of (co) variance, the additive genetic functions and permanent environmental covariance were estimated by random regression in three functions: Wilmink, Legendre Polynomials (third order) and Linear spline Polynomials (three knots). The Legendre polynomial function showed better fit quality. The genetic and permanent environment variances for TDMY ranged from 2.67 to 5.14 and from 9.31 to 12.04, respectively. Heritability estimates gradually increased from the beginning (0.13) to mid-lactation (0.19). The genetic correlations between the days of the control ranged from 0.37 to 1.00. The correlations of permanent environment followed the same trend as genetic correlations. The use of Legendre polynomials via random regression model can be considered as a good tool for estimating genetic parameters for test-day milk yield records.

Genetic parameters for test-day milk yield, lactation persistency, and fertility in low-, medium-, and high-production systems in Kenya

Genetic parameters for test-day milk yield, lactation persistency, and age at first calving (as a fertility trait) were estimated for the first 4 lactations in multiple-breed dairy cows in low-, medium-, and high-production systems in Kenya. Data included 223,285 test-day milk yield records from 11,450 cows calving from 1990 to 2015 in 148 herds. A multivariate random regression model was used to estimate variance and covariance components. The fixed effects in the model included herd, year, and test month, and age as a covariate. The lactation profile over days in milk (DIM) was fitted as a cubic smoothing spline. Random effects included herd, year, and test month interaction effects, genetic group effects, and additive genetic and permanent environmental effects modeled with a cubic Legendre polynomial function. The residual variance was heterogeneous with 11 classes. Consequently, the variance components were varied over the lactation and with the production system. The estimated heritability for milk yield was lower in the low-production system (0.04-0.48) than in the medium- (0.22-0.59) and high-production (0.21-0 60) systems. The genetic correlations estimated between different DIM within lactations decreased as the time interval increased, becoming negative between the ends of the lactations in the low- and medium-production systems. Low (0.05) to medium (0.60) genetic correlations were estimated among first lactation test-day milk yields across the 3 production systems. Genetic correlations between the first lactation test-day milk yield and age at first calving ranged from 0.27 to 0.49, 0 to 0.81, and -0.08 to 0.27 in the low-, medium-, and high-production systems, respectively. Medium to high heritabilities (0.17-0.44) were estimated for persistency, with moderate to high (0.30-0.87) genetic correlations between 305-d milk yield and persistency. This indicates that genetic improvement in persistency would lead to increased milk yield. The low to medium genetic correlations between test-day milk yield between production systems indicate that sires may be re-ranked between production systems. Therefore, we conclude that sires should be selected based on a genetic evaluation within the target production system.

Genotype by environment interaction for somatic cell score in Holstein cattle of southern Brazil via reaction norms

Objective

The objective of this study was to evaluate the genetic behavior of a population of Holstein cattle in response to the variation of environmental temperature by analyzing the effects of genotype by environment interaction (GEI) through reaction norms for the somatic cell score (SCS).

Methods

Data was collected for 67,206 primiparous cows from the database of the Paraná Holstein Breeders Association in Brazil, with the aim of evaluating the temperature effect, considered as an environmental variable, distinguished under six gradients, with the variation range found being 17°C to 19.5°C, over the region. A reaction norm model was adopted utilizing the fourth order under the Legendre polynomials, using the mixed models of analysis by the restricted maximum likelihood method by the WOMBAT software. Additionally, the genetic behavior of the 15 most representative bulls was assessed, in response to the changes in the temperature gradient.

Results

A mean score of 2.66 and a heritability variation from 0.17 to 0.23 was found in the regional temperature increase. The correlation between the environmental gradients proved to be higher than 0.80. Distinctive genetic behaviors were observed according to the increase in regional temperature, with an observed increase of up to 0.258 in the breeding values of some animals, as well as a reduction in the breeding of up to 0.793, with occasional reclassifications being observed as the temperature increased.

Conclusion

Non-relevant GEI for SCS were observed in Holstein cattle herds of southern Brazil. Thus, the inclusion of the temperature effect in the model of genetic evaluation of SCS for the southern Brazilian Holstein breed is not required.