Multiple-trait Random Regression Research Articles

Genetic Selection Approach for Semen Characteristics in Thai Native Grandparent Roosters (Pradu Hang Dum) Using Random Regression Test-Day Models and Selection Indices.

The aim of this study was to analyze suitable genetic models and selection indices to estimate the genetic parameters and breeding values of native Thai roosters. A total of 3475 records of seven semen traits (mass movement, semen pH, semen color, volume, sperm viability, sperm abnormalities, and sperm concentration) from 242 Thai native grandparent roosters were analyzed. Multiple-trait random regression test-day models with five covariance functions were used to analyze the variance components, genetic parameters, and breeding values. The selection index (SI) was calculated to determine the optimal genetic value for different selection percentages. The results showed that a multiple-trait random regression test-day model with a second-order Legendre polynomial function was the most appropriate genetic model for this population. The estimated heritability values were low to moderate, ranging from 0.110 to 0.112 (mass movement), 0.040 to 0.051 (semen pH), 0.092 to 0.097 (semen color), 0.220 to 0.225 (semen volume), 0.067 to 0.083 (sperm viability), 0.086 to 0.099 (sperm abnormalities), and 0.134 to 0.138 (sperm concentration). The repeatability values exceeded the heritability values and were within the range of 0.133 to 0.688. The genetic correlations among semen traits ranged from -0.332 to 0.677, and phenotypic correlations ranged from -0.260 to 0.460. When considering heritability and genetic correlation values, semen volume, sperm concentration, and mass movement were the top three priority semen traits calculated as selection indices. Finally, the top 10% of the selection index was recommended for creating the next generation. Our findings provide useful information on genetic parameters and an appropriate selection index of semen traits for selecting the genetics of individual Thai native grandparent roosters. The heritability estimates for semen traits reported here suggest an adequate response to selection through a genetic evaluation approach. Our results indicate that it is possible to select grandparent roosters with better reproductive performance.

Estimation of genetic parameters for feed efficiency traits using random regression models in dairy cattle

Feed efficiency has become an increasingly important research topic in recent years. As feed costs rise and the environmental impacts of agriculture become more apparent, improving the efficiency with which dairy cows convert feed to milk is increasingly important. However, feed intake is expensive to measure accurately on large populations, making the inclusion of this trait in breeding programs difficult. Understanding how the genetic parameters of feed efficiency and traits related to feed efficiency vary throughout the lactation period is valuable to gain understanding into the genetic nature of feed efficiency. This study used 121,226 dry matter intake (DMI) records, 120,500 energy-corrected milk (ECM) records, and 98,975 metabolic body weight (MBW) records, collected on 7,440 first-lactation Holstein cows from 6 countries (Canada, Denmark, Germany, Spain, Switzerland, and the United States), from January 2003 to February 2022. Genetic parameters were estimated using a multiple-trait random regression model with a fourth-order Legendre polynomial for all traits. Weekly phenotypes for DMI were re-parameterized using linear regressions of DMI on ECM and MBW, creating a measure of feed efficiency that was genetically corrected for ECM and MBW, referred to as genomic residual feed intake (gRFI). Heritability (SE) estimates varied from 0.15 (0.03) to 0.29 (0.02) for DMI, 0.24 (0.01) to 0.29 (0.03) for ECM, 0.55 (0.03) to 0.83 (0.05) for MBW, and 0.12 (0.03) to 0.22 (0.06) for gRFI. In general, heritability estimates were lower in the first stage of lactation compared with the later stages of lactation. Additive genetic correlations between weeks of lactation varied, with stronger correlations between weeks of lactation that were close together. The results of this study contribute to a better understanding of the change in genetic parameters across the first lactation, providing insight into potential selection strategies to include feed efficiency in breeding programs.

Multiple-trait random regression modeling of feed efficiency in US Holsteins

Residual feed intake (RFI) and feed saved (FS) are important feed efficiency traits that have been increasingly considered in genetic improvement programs. Future sustainability of these genetic evaluations will depend upon greater flexibility to accommodate sparsely recorded dry matter intake (DMI) records on many more cows, especially from commercial environments. Recent multiple-trait random regression (MTRR) modeling developments have facilitated days in milk (DIM)-specific inferences on RFI and FS, particularly in modeling the effect of change in metabolic body weight (MBW). The MTRR analyses, using daily data on the core traits of DMI, MBW, and milk energy (MilkE), were conducted separately for 2,532 primiparous and 2,379 multiparous US Holstein cows from 50 to 200 DIM. Estimated MTRR variance components were used to derive genetic RFI and FS and DIM-specific genetic partial regressions of DMI on MBW, MilkE, and change in MBW. Estimated daily heritabilities of RFI and FS varied across lactation for both primiparous (0.05-0.07 and 0.11-0.17, respectively) and multiparous (0.03-0.13 and 0.10-0.17, respectively) cows. Genetic correlations of RFI across DIM varied (>0.05) widely compared with FS (>0.54) within either parity class. Heritability estimates based on average lactation-wise measures were substantially larger than daily heritabilities, ranging from 0.17 to 0.25 for RFI and from 0.35 to 0.41 for FS. The partial genetic regression coefficients of DMI on MBW (0.11 to 0.16 kg/kg0.75 for primiparous and 0.12 to 0.14 kg/kg0.75 for multiparous cows) and of DMI on MilkE (0.45 to 0.68 kg/Mcal for primiparous and 0.36 to 0.61 kg/Mcal for multiparous cows) also varied across lactation. In spite of the computational challenges encountered with MTRR, the model potentially facilitates an efficient strategy for harnessing more data involving a wide variety of data recording scenarios for genetic evaluations on feed efficiency.

Determining Heat Stress Effects of Multiple Genetic Traits in Tropical Dairy Cattle Using Single-Step Genomic BLUP.

Heat stress is becoming a significant problem in dairy farming, especially in tropical countries, making accurate genetic selection for heat tolerance a priority. This study investigated the effect of heat stress manifestation on genetics for milk yield, milk quality, and dairy health traits with and without genomic information using single-step genomic best linear unbiased prediction (ssGBLUP) and BLUP in Thai−Holstein crossbred cows. The dataset contained 104,150 test-day records from the first lactation of 15,380 Thai−Holstein crossbred cows. A multiple-trait random regression test-day model on a temperature−humidity index (THI) function was used to estimate the genetic parameters and genetic values. Heat stress started at a THI of 76, and the heritability estimates ranged from moderate to low. The genetic correlation between those traits and heat stress in both BLUP methods was negative. The accuracy of genomic predictions in the ssGBLUP method was higher than the BLUP method. In conclusion, heat stress negatively impacted milk production, increased the somatic cell score, and disrupted the energy balance. Therefore, in dairy cattle genetic improvement programs, heat tolerance is an important trait. The new genetic evaluation method (ssGBLUP) should replace the traditional method (BLUP) for more accurate genetic selection.

Genetic analysis of milk urea concentration and its genetic relationship with selected traits of interest in dairy cows

The aim of this study was to estimate genetic parameters of milk urea concentration (MU) and its genetic correlations with milk production traits, longevity, and functional traits in the first 3 parities in dairy cows. The edited data set consisted in 9,107,349 MU test-day records from the first 3 parities of 560,739 cows in 2,356 herds collected during the years 1994 to 2020. To estimate the genetic parameters of MU, data of 109 randomly selected herds, with a total of 770,016 MU test-day records, were used. Genetic parameters and estimated breeding values were estimated using a multiple-trait (parity) random regression model. Herd-test-day, age-year-season of calving, and days in milk classes (every 5 d as a class) were used as fixed effects, whereas effects of herd-year of calving, permanent environment, and animal were modeled using random regressions and Legendre polynomials of order 2. The average daily heritability and repeatability of MU during days in milk 5 to 365 in the first 3 parities were 0.19, 0.22, 0.20, and 0.48, 0.48, 0.47, respectively. The mean genetic correlation estimated among MU in the first 3 parities ranged from 0.96 to 0.97. The average daily estimated breeding values for MU of the selected bulls (n = 1,900) ranged from -9.09 to 7.37 mg/dL. In the last 10 yr, the genetic trend of MU has gradually increased. The genetic correlation between MU and 11 traits of interest ranged from -0.28 (milk yield) to 0.28 (somatic cell score). The findings of this study can be used as the first step for development of a routine genetic evaluation for MU and its inclusion into the genetic selection program in the Walloon Region of Belgium.

Application of single-step genomic evaluation using multiple-trait random regression test-day models in dairy cattle

Test-day traits are important for genetic evaluation in dairy cattle and are better modeled by multiple-trait random regression models (RRM). The reliability and bias of genomic estimated breeding values (GEBV) predicted using multiple-trait RRM via single-step genomic best linear unbiased prediction (ssGBLUP) were investigated in the 3 major dairy cattle breeds in Canada (i.e., Ayrshire, Holstein, and Jersey). Individual additive genomic random regression coefficients for the test-day traits were predicted using 2 multiple-trait RRM: (1) one for milk, fat, and protein yields in the first, second, and third lactations, and (2) one for somatic cell score in the first, second, and third lactations. The predicted coefficients were used to derive GEBV for each lactation day and, subsequently, the daily GEBV were compared with traditional daily parent averages obtained by BLUP. To ensure compatibility between pedigree and genomic information for genotyped animals, different scaling factors for combining the inverse of genomic (G-1) and pedigree (A-122) relationship matrices were tested. In addition, the inclusion of only genotypes from animals with accurate breeding values (defined in preliminary analysis) was compared with the inclusion of all available genotypes in the analyzes. The ssGBLUP model led to considerably larger validation reliabilities than the BLUP model without genomic information. In general, scaling factors used to combine the G-1 and A-122 matrices had small influence on the validation reliabilities. However, a greater effect was observed in the inflation of GEBV. Less inflated GEBV were obtained by the ssGBLUP compared with the parent average from traditional BLUP when using optimal scaling factors to combine the G-1 and A-122 matrices. Similar results were observed when including either all available genotypes or only genotypes from animals with accurate breeding values. These findings indicate that ssGBLUP using multiple-trait RRM increases reliability and reduces bias of breeding values of young animals when compared with parent average from traditional BLUP in the Canadian Ayrshire, Holstein, and Jersey breeds.

Estimates of genetic trends for daily milk yield and somatic cell score of primiparous Holstein cattle in Iran

ABSTRACT Genetic and permanent environmental correlations between milk yield and somatic cell score (SCS) were estimated. The data set used in this analyse contained a total of 204,438 test-day (TD) observations of milk yield and SCS from 22,903 primiparous daughters of 554 dairy Holstein sires. A multiple-trait random regression TD model was employed in the analyses using the Gibbs sampling method. Obtained results showed that daily heritabilities were highest in the second half of lactation and varied from 0.086 to 0.250 for milk yield and 0.041 to 0.086 for SCS. Undesirable genetic relationships with milk yield were found for SCS with genetic correlations ranging from −0.093 to 0.319 that the highest was estimated between SCS days in milk (DIM) 305 and milk DIM 5. Based on these results, we can conclude that the cows with high yield potential in the initial stage of lactation had the highest genetic potential for SCS at the end of lactation. This emphasizes the need for including udder health traits in the breeding goal. However results of this investigation, based on genetic trends, indicate that more emphasis is paid to increase milk yield of Holstein cows in Iran with no or a little emphasis on SCC.

Short communication: Single-step genomic evaluation of milk production traits using multiple-trait random regression model in Chinese Holsteins

The objectives of this study were to evaluate the prediction performance of the single-step genomic BLUP method using a multi-trait random regression model in genomic evaluation for milk production traits of Chinese Holsteins, and investigate how parameters w, τ, and ω used in the construction of the combined relationship matrix (H) affected prediction accuracy and bias. A total of 2.8 million test-day records from 0.2 million cows were available for milk, protein, and fat yields. Pedigree information included 0.3 million animals and 7,577 of them were genotyped with medium-density single nucleotide polymorphism marker panels. Genotypes were imputed into Geneseek Genomic Profiler HD (GeneSeek, Lincoln, NE) including 77K markers. A reduced data set for evaluating models was extracted from the full data set by removing test-day records from the last 4 yr. Bull and cow validation populations were constructed for each trait. We evaluated the prediction performance of the multiple-trait multiple-lactation random regression single-step genomic BLUP (RR-ssGBLUP) models with different values of parameters w, τ, and ω in the H matrix, taking consideration of inbreeding. We compared RR-ssGBLUP with the multiple-trait multiple-lactation random regression model based on pedigree and genomic BLUP. De-regressed proofs for 305-d milk, protein, and fat yields averaged over 3 lactations, which were calculated from the full data set, were used for posteriori validations. The results showed that RR-ssGBLUP was feasible for implementation in breeding practice, and its prediction performance was superior to the other 2 methods in the comparison, including prediction accuracy and unbiasedness. For bulls, RR-ssGBLUP models with w0.05τ2.0ω1.0,w0.05τ2.5ω1.0, and w0.1τ1.6ω0.4 achieved the best performance for milk, protein, and fat yields, respectively. For cows, the RR-ssGBLUP with w0.2τ1.6ω0.4 performed the best for all 3 traits. The H matrix constructed with larger τ and smaller ω gave better convergence in solving mixed model equations. Among different RR-ssGBLUP models, the differences in validation accuracy were small. However, the regression coefficient indicating prediction bias varied substantially. The increase of w and τ, and decrease of ω, led to an increase in the regression coefficient. The results demonstrated RR-ssGBLUP is a good alternative to the multi-step approach, but the optimal choice of parameters should be found via preliminary validation study to achieve the best performance.

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (-0.62 to -0.59) and with protein yield (-0.32 to -0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.

Genetic correlations between the cumulative pseudo-survival rate, milk yield, and somatic cell score during lactation in Holstein cattle in Japan using a random regression model

Trends in genetic correlations between longevity, milk yield, and somatic cell score (SCS) during lactation in cows are difficult to trace. In this study, changes in the genetic correlations between milk yield, SCS, and cumulative pseudo-survival rate (PSR) during lactation were examined, and the effect of milk yield and SCS information on the reliability of estimated breeding value (EBV) of PSR were determined. Test day milk yield, SCS, and PSR records were obtained for Holstein cows in Japan from 2004 to 2013. A random subset of the data was used for the analysis (825 herds, 205,383 cows). This data set was randomly divided into 5 subsets (162-168 herds, 83,389-95,854 cows), and genetic parameters were estimated in each subset independently. Data were analyzed using multiple-trait random regression animal models including either the residual effect for the whole lactation period (H0), the residual effects for 5 lactation stages (H5), or both of these residual effects (HD). Milk yield heritability increased until 310 to 351 d in milk (DIM) and SCS heritability increased until 330 to 344 DIM. Heritability estimates for PSR increased with DIM from 0.00 to 0.05. The genetic correlation between milk yield and SCS increased negatively to under -0.60 at 455 DIM. The genetic correlation between milk yield and PSR increased until 342 to 355 DIM (0.53-0.57). The genetic correlation between the SCS and PSR was -0.82 to -0.83 at around 180 DIM, and decreased to -0.65 to -0.71 at 455 DIM. The reliability of EBV of PSR for sires with 30 or more recorded daughters was 0.17 to 0.45 when the effects of correlated traits were ignored. The maximum reliability of EBV was observed at 257 (H0) or 322 (HD) DIM. When the correlations of PSR with milk yield and SCS were considered, the reliabilities of PSR estimates increased to 0.31-0.76. The genetic parameter estimates of H5 were the same as those for HD. The rank correlation coefficients of the EBV of PSR between H0 and H5 or HD were greater than 0.9. Additionally, the reliabilities of EBV of PSR of H0 were similar to those for H5 and HD. Therefore, the genetic parameter estimates in H0 were not substantially different from those in H5 and HD. When milk yield and SCS, which were genetically correlated with PSR, were used, the reliability of PSR increased. Estimates of the genetic correlations between PSR and milk yield and between PSR and SCS are useful for management and breeding decisions to extend the herd life of cows.

Bayesian Models combining Legendre and B-spline polynomials for genetic analysis of multiple lactations in Gyr cattle

We aimed with this study to combine Legendre polynomials (LEG) and linear B-splines (BSP) to describe simultaneously the first and second lactation of Gyr dairy cattle under a multiple-trait random regression models (MTRRM) framework. Additionally we proposed the application of self-organizing map to define the classes of residual variances under these models. A total of 26,438 and 23,892 milk yield test-day records were used, respectively, for the first and second lactations of 3253 Gyr cows. Two preliminary MTRRM analyses considering 10 residual classes were performed: the first one was based on LEG for systematic and random effects for both lactations; and the second one was based on BSP. Three classes were defined by using a self-organizing map: from 6 to 35; 36–185 and 186–305 days in milk. After definition of residual variance classes, a total of 16 MTRRM combining LEG and BSP were compared. The MTRRM based on BSP to describe the systematic effects of the first and second lactation, BSP to describe the random effects of the first lactation and LEG to describe the random effects of the second lactation (BSP-BSP-BSP-LEG) outperformed all other models. From the BSP-BSP-BSP-LEG model, heritability estimates for milk yield over time ranged from 0.1107 to 0.2902, and from 0.2036 to 0.3967, for the first and second lactation, respectively. In general, additive genetic correlation estimates between days in milk within each lactation and between lactations had medium magnitude (mean of genetic correlations were 0.6630, 0.6226 and 0.4749 for the first, second and between both lactations, respectively). We concluded that combining different functions under a MTRRM framework is a feasible alternative for genetic modeling of lactation curves in Gyr dairy cattle.

Genetic correlation between growth and reproductive performance of beef females depends on environment

In tropical production systems, beef cattle are raised in highly heterogeneous environments. Heterogeneity is, therefore, expected to exist in the (co)variance components for traits of economic interest in different production environments. The main objective of the present study was to estimate genetic correlations between growth traits and reproductive performance of beef females, depending on the environment. The present study was conducted in the tropical region of Brazil, applying a multiple-trait random regression animal model to field records of heifer pregnancy (HP), hip height, bodyweight at ~18 months of age (BW18) and postweaning weight gain (PWG) from 20 893 Nelore females. As evidence of genotype by environment interaction (G × E), heterogeneity of genetic variance across environments was observed mainly for HP, PWG and BW18. Moreover, the estimates of genetic correlation within these traits reached values lower than unity on the environmental gradient. The genetic correlation among growth traits tended to be stronger in favourable environments, a fact that should favour correlated responses under these conditions. In contrast, the genetic correlations between growth traits and HP tended to become weaker and even exhibited little evidence of antagonism in more favourable environments. On the basis of these findings, selection for higher growth in extreme favourable environments should result in little or no damage to HP as a correlated response. All these results lead us to believe that the G × E is an important factor to be considered in genetic evaluations of beef cattle raised in tropical environments.

Multiple-trait structured antedependence model to study the relationship between litter size and birth weight in pigs and rabbits

BackgroundSome genetic studies need to take into account correlations between traits that are repeatedly measured over time. Multiple-trait random regression models are commonly used to analyze repeated traits but suffer from several major drawbacks. In the present study, we developed a multiple-trait extension of the structured antedependence model (SAD) to overcome this issue and validated its usefulness by modeling the association between litter size (LS) and average birth weight (ABW) over parities in pigs and rabbits.MethodsThe single-trait SAD model assumes that a random effect at time t_{j} can be explained by the previous values of the random effect (i.e. at previous times). The proposed multiple-trait extension of the SAD model consists in adding a cross-antedependence parameter to the single-trait SAD model. This model can be easily fitted using ASReml and the OWN Fortran program that we have developed. In comparison with the random regression model, we used our multiple-trait SAD model to analyze the LS and ABW of 4345 litters from 1817 Large White sows and 8706 litters from 2286 L-1777 does over a maximum of five successive parities.ResultsFor both species, the multiple-trait SAD fitted the data better than the random regression model. The difference between AIC of the two models (AIC_random regression-AIC_SAD) were equal to 7 and 227 for pigs and rabbits, respectively. A similar pattern of heritability and correlation estimates was obtained for both species. Heritabilities were lower for LS (ranging from 0.09 to 0.29) than for ABW (ranging from 0.23 to 0.39). The general trend was a decrease of the genetic correlation for a given trait between more distant parities. Estimates of genetic correlations between LS and ABW were negative and ranged from −0.03 to −0.52 across parities. No correlation was observed between the permanent environmental effects, except between the permanent environmental effects of LS and ABW of the same parity, for which the estimate of the correlation was strongly negative (ranging from −0.57 to −0.67).ConclusionsWe demonstrated that application of our multiple-trait SAD model is feasible for studying several traits with repeated measurements and showed that it provided a better fit to the data than the random regression model.

Combining different functions to describe milk, fat, and protein yield in goats using Bayesian multiple-trait random regression models.

We proposed multiple-trait random regression models (MTRRM) combining different functions to describe milk yield (MY) and fat (FP) and protein (PP) percentage in dairy goat genetic evaluation by using Bayesian inference. A total of 3,856 MY, FP, and PP test-day records, measured between 2000 and 2014, from 535 first lactations of Saanen and Alpine goats, including their cross, were used in this study. The initial analyses were performed using the following single-trait random regression models (STRRM): third- and fifth-order Legendre polynomials (Leg3 and Leg5), linear B-splines with 3 and 5 knots, the Ali and Schaeffer function (Ali), and Wilmink function. Heterogeneity of residual variances was modeled considering 3 classes. After the selection of the best STRRM to describe each trait on the basis of the deviance information criterion (DIC) and posterior model probabilities (PMP), the functions were combined to compose the MTRRM. All combined MTRRM presented lower DIC values and higher PMP, showing the superiority of these models when compared to other MTRRM based only on the same function assumed for all traits. Among the combined MTRRM, those considering Ali to describe MY and PP and Leg5 to describe FP (Ali_Leg5_Ali model) presented the best fit. From the Ali_Leg5_Ali model, heritability estimates over time for MY, FP. and PP ranged from 0.25 to 0.54, 0.27 to 0.48, and 0.35 to 0.51, respectively. Genetic correlation between MY and FP, MY and PP, and FP and PP ranged from -0.58 to 0.03, -0.46 to 0.12, and 0.37 to 0.64, respectively. We concluded that combining different functions under a MTRRM approach can be a plausible alternative for joint genetic evaluation of milk yield and milk constituents in goats.

Genetic analysis of fat-to-protein ratio, milk yield and somatic cell score of Holstein cows in Japan in the first three lactations by using a random regression model.

We estimated the genetic parameters of fat-to-protein ratio (FPR) and the genetic correlations between FPR and milk yield or somatic cell score in the first three lactations in dairy cows. Data included 3,079,517 test-day records of 201,138 Holstein cows in Japan from 2006 to 2011. Genetic parameters were estimated with a multiple-trait random regression model in which the records within and between parities were treated as separate traits. The phenotypic values of FPR increased soon after parturition and peaked at 10 to 20 days in milk, then decreased slowly in mid- and late lactation. Heritability estimates for FPR yielded moderate values. Genetic correlations of FPR among parities were low in early lactation. Genetic correlations between FPR and milk yield were positive and low in early lactation, but only in the first lactation. Genetic correlations between FPR and somatic cell score were positive in early lactation and decreased to become negative in mid- to late lactation. By using these results for genetic evaluation it should be possible to improve energy balance in dairy cows.

Modeling lactation curves and estimation of genetic parameters in Holstein cows using multiple-trait random regression models.

We compared the goodness of fit of three mathematical functions (including: Legendre polynomials, Lidauer-Mäntysaari function and Wilmink function) for describing the lactation curve of primiparous Iranian Holstein cows by using multiple-trait random regression models (MT-RRM). Lactational submodels provided the largest daily additive genetic (AG) and permanent environmental (PE) variance estimates at the end and at the onset of lactation, respectively, as well as low genetic correlations between peripheral test-day records. For all models, heritability estimates were highest at the end of lactation (245 to 305 days) and ranged from 0.05 to 0.26, 0.03 to 0.12 and 0.04 to 0.24 for milk, fat and protein yields, respectively. Generally, the genetic correlations between traits depend on how far apart they are or whether they are on the same day in any two traits. On average, genetic correlations between milk and fat were the lowest and those between fat and protein were intermediate, while those between milk and protein were the highest. Results from all criteria (Akaike's and Schwarz's Bayesian information criterion, and -2*logarithm of the likelihood function) suggested that a model with 2 and 5 coefficients of Legendre polynomials for AG and PE effects, respectively, was the most adequate for fitting the data.

Estimation of genetic parameters for daily milk yields of primiparous Iranian Holstein cows

Abstract. Applying a multiple trait random regression (MT-RR) in national level and for whole test day records of a country is a great advance in animal breeding context. Having reliable (co) variance components is a critical step in applying multiple traits genetic evaluation especially in developing countries. Genetic parameters of milk, fat and protein yields were estimated for Iranian Holstein dairy cows. Data included 276 692 test day (TD) production traits records collected of 30 705 primiparous cows belonging to 619 sires. An animal multi-trait random regression model was employed in the analyses using the restricted maximum likelihood (REML) method. The model included herd-test-date, age-season of calving (by applying a fixed regression for each subclass of this effect) and year of calving as fixed effects and random regression (RR) coefficients for additive genetic (AG) and permanent environmental (PE) effects. Obtained results showed that daily heritabilities ranged from 0.10 to 0.21 for milk, from 0.05 to 0.08 for fat and from 0.08 to 0.18 for protein yield. Estimated heritability for 305-d milk, fat and protein yields were 0.25, 0.20 and 0.25, respectively. Correlations between individual test day records within traits were high for adjacent tests (nearly 1) and decreased as the interval between tests increased. Correlations between yields of milk, fat and protein on a given test day are also high and greater during late lactation than during early or mid-lactation. Genetic correlations between 305-d yield traits ranged from 0.75 to 0.92. The largest genetic correlation, as well as permanent environmental correlation, was observed between milk and protein yield.

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21–0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal.

Genetic and environmental relationships between body condition score and milk production traits in Canadian Holsteins

The objective of this research was to estimate genetic parameters of first-lactation body condition score (BCS), milk yield, fat percentage (Fat%), protein percentage (Prot%), somatic cell score (SCS), milk urea nitrogen (MUN), lactose percentage (Lact%), and fat to protein ratio (F:P) using multiple-trait random regression animal models. Changes in covariances between BCS and milk production traits on a daily basis have not been investigated before and could be useful for determining which BCS estimated breeding values (EBV) might be practical for selection in the future. Field staff from Valacta milk recording agency (Sainte-Anne-de-Bellevue, QC, Canada) collected BCS from Québec herds several times per cow throughout the lactation. Average daily heritabilities and genetic correlations among the various traits were similar to literature values. On an average daily basis, BCS was genetically unfavorably correlated with milk yield (i.e., increased milk yield was associated with lower body condition). The unfavorable genetic correlation between BCS and milk yield became stronger as lactation progressed, but was equivalent to zero for the first month of lactation. Favorable genetic correlations were found between BCS with Prot%, SCS, and Lact% (i.e., greater BCS was associated with greater Prot%, lower SCS, and greater Lact%). These correlations were strongest in early lactation. On an average daily basis, BCS was not genetically correlated with Fat% or MUN, but was negatively correlated with F:P. Furthermore, BCS at 5 and 50 d in milk (DIM) had the most favorable genetic correlations with milk production traits over the lactation (at 5, 50, 150, and 250 DIM). Thus, early lactation BCS EBV shows potential for selection. Regardless, this study showed that the level of association BCS has with milk production traits is not constant over the lactation. Simultaneous selection for both BCS and milk production traits should be considered, mainly due to the unfavorable genetic correlation between BCS with milk yield.

Alternative parameterizations of the multiple-trait random regression model for milk yield and somatic cell score via recursive links between phenotypes

Multiple-trait random regression models with recursive phenotypic link from somatic cell score (SCS) to milk yield on the same test day and with different restrictions on co-variances between these traits were fitted to the first-lactation Canadian Holstein data. Bayesian methods with Gibbs sampling were used to derive inferences about parameters for all models. Bayes factor indicated that the recursive model with uncorrelated environmental effects between traits was the most plausible specification in describing the data. Goodness of fit in terms of a within-trait weighted mean square error and correlation between observed and predicted data was the same for all parameterizations. All recursive models estimated similar negative causal effects from SCS to milk yield (up to -0.4 in 46-115 days in milk in lactation). Estimates of heritabilities, genetic and environmental correlations for the first two regression coefficients (overall level of a trait and lactation persistency) within both traits were similar among models. Genetic correlations between milk and SCS were dependent on the restrictions on genetic co-variances for these traits. Recursive model with uncorrelated system genetic effects between milk and SCS gave estimates of genetic correlations of the opposite sign compared with a regular multiple-trait model. Phenotypic recursion between milk and SCS seemed, however, to be the only source of environmental correlations between these two traits. Rankings of sires for total milk yield in lactation, average daily SCS and persistency for both traits were similar among models. Multiple-trait model with recursive links between milk and SCS and uncorrelated random environmental effects could be an attractive alternative for a regular multiple-trait model in terms of model parsimony and accuracy.