Deriving Lactation Yields from Test-Day Yields Adjusted for Lactation Stage, Age, Pregnancy, and Herd Test Date

Lactation records for milk, fat, and protein yields were calculated from test-day data adjusted for the effects of lactation stage, age, previous days open, days pregnant, and test-day class (herd, test date, and milking frequency). Those lactation records reflect the improved accounting of environmental effects from a test-day model and can be combined with historical lactation records. Test-day data were adjusted with existing lactation multiplicative adjustments to maintain variance characteristics. Then additive adjustments for lactation stage, age, previous days open, and days pregnant were applied. The current multiplicative adjustments for previous days open were not applied because its effect was expected to differ by lactation stage. To remove genetic differences, the estimated breeding value from the previous evaluation divided by 305 was subtracted. Effects of test-day class, and permanent environment within and across parities were estimated within herd. The effect of test-day class was subtracted from adjusted test-day yield, and the breeding value restored. Those deviations then were combined with the best prediction procedure into a lactation measure. Heritabilities and repeatabilities of lactation records that were adjusted for test-day class were higher than for current lactation records. The adjusted records should improve the accuracy of evaluations and allow the use of test-day data as well as provide for the continued use of historical data when test-day data are not available.     

Estimated heterogeneity of phenotypic variance of test-day yield with a structural variance model

First-lactation test-day milk, fat, and protein yields from New York, Wisconsin, and California herds from 1990 through 2000 were adjusted additively for age and lactation stage. A random regression model with third-order Legendre polynomials for permanent environmental and genetic effects was used. The model included a random effect with the same polynomial regressions for 2 yr of calvings within herd (herd-time effect) to provide herd-specific lactation curves that can change every 2 yr. (Co)variance components were estimated using expectation-maximization REML simultaneously with phenotypic variances that were modeled using a structural variance model. Maximum heritability for test-day milk yield was estimated to be approximately 20% around 200 to 250 d in milk; heritabilities were slightly lower for test-day fat and protein yields. Herd-time effects explained 12 to 20% of phenotypic variance and had the greatest impact at start of lactation. Variances of test-day yields increased with time, subclass size, and milking frequency. Test month had limited influence on variance. Variance increased for cows in herds with low and high milk yields and for early and late lactation stages. Repeatabilities of variances observed for a given class of herd, test-day, and milking frequency were 14 to 17% across nested variance subclasses based on lactation stage.     

Genetic parameters of milking frequency and milk production traits in Canadian Holsteins milked by an automated milking system

Twice-a-day milking is currently the most frequently used milking schedule in Canadian dairy cattle. However, with an automated milking system (AMS), dairy cows can be milked more frequently. The objective of this study was to estimate genetic parameters for milking frequency and for production traits of cows milked within an AMS. Data were 141,927 daily records of 953 primiparous Holstein cows from 14 farms in Ontario and Quebec. Most cows visited the AMS 2 (46%) or 3 (37%) times a day. A 2-trait [daily (24-h) milking frequency and daily (24-h) milk yield] random regression daily animal model and a multiple-trait (milk, fat, protein yields, somatic cell score, and milking frequency) random regression test-day animal model were used for the estimation of (co)variance components. Both models included fixed effect of herd × test-date, fixed regressions on days in milk (DIM) nested within age at calving by season of calving, and random regressions for additive genetic and permanent environmental effects. Both fixed and random regressions were fitted with fourth-order Legendre polynomials on DIM. The number of cows in the multiple-trait test-day model was smaller compared with the daily animal model. Heritabilities from the daily model for daily (24-h) milking frequency and daily (24-h) milk yield ranged between 0.02 and 0.08 and 0.14 and 0.20, respectively. Genetic correlations between daily (24-h) milk yield and daily (24-h) milking frequency were largest at the end of lactation (0.80) and smallest in mid-lactation (0.27). Heritabilities from the test-day model for test-day milking frequency, milk, fat and protein yield, and somatic cell score were 0.14, 0.26, 0.20, 0.21, and 0.20, respectively. The genetic correlation was positive between test-day milking frequency and official test-day milk, fat, and protein yields, and negative between official test-day somatic cell score and test-day milking frequency.     

Genetic components of heat stress for dairy cattle with multiple lactations

Data included 585,119 test-day records for milk, fat, and protein yields from the first, second, and third parities of 38,608 Holsteins in Georgia. Daily temperature-humidity indexes (THI) were available from public weather stations. Models included a repeatability test-day model with a random regression on a function of THI and a test-day random regression model using linear splines with knots at 5, 50, 200, and 305 d in milk and a function of THI. Random effects were additive genetic and permanent environmental in the repeatability model and additive genetic, permanent environmental, and herd year in the random regression model. Additionally, models included fixed effects for herd test day, calving age, milking frequency, and lactation stage. Phenotypic variance increased by 50 to 60% from the first to second parity for all yield traits with the repeatability model and by 12 to 15% from the second to third parity. General additive genetic variance increased by 25 to 35% from the first to second parity for all yield traits but decreased slightly from the second to third parity for milk and protein yields. Genetic variance for heat tolerance doubled from the first to second parity and increased by 20 to 100% from the second to third parity. Genetic correlations among general additive effects were lowest between the first and second parities (0.84 to 0.88) and were highest between the second and third parities (0.96 to 0.98). Genetic correlations among parities for the effect of heat tolerance ranged from 0.56 to 0.79. Genetic correlations between general and heat-tolerance effects across parities and yield traits ranged from −0.30 to −0.50. With the random regression model, genetic variance for heat tolerance for milk yield was approximately one-half that of the repeatability model. For milk yield, the most negative genetic correlation (approximately −0.45) between general and heat-tolerance effects was between 50 and 200 d in milk for the first parity and between 200 and 305 d in milk for the second and third parities. The genetic variance of heat tolerance increased substantially from the first to third parity. Genetic estimates of heat tolerance may be inflated with the repeatability model because of timing of lactations to avoid peak yield during hot seasons.     

Genetic effects of heat stress on milk yield of Thai Holstein crossbreds

The threshold for heat stress on milk yield of Holstein crossbreds under climatic conditions in Thailand was investigated, and genetic effects of heat stress on milk yield were estimated. Data included 400,738 test-day milk yield records for the first 3 parities from 25,609 Thai crossbred Holsteins between 1990 and 2008. Mean test-day milk yield ranged from 12.6 kg for cows with <87.5% Holstein genetics to 14.4 kg for cows with ≥93.7% Holstein genetics. Daily temperature and humidity data from 26 provincial weather stations were used to calculate a temperature-humidity index (THI). Test-day milk yield varied little with THI for first parity except above a THI of 82 for cows with ≥93.7% Holstein genetics. For third parity, test-day milk yield started to decline after a THI of 74 for cows with ≥87.5% Holstein genetics and declined more rapidly after a THI of 82. A repeatability test-day model with parities as correlated traits was used to estimate heat stress parameters; fixed effects included herd-test month-test year and breed groups, days in milk, calving age, and parity; random effects included 2 additive genetic effects, regular and heat stress, and 2 permanent environment, regular and heat stress. The threshold for effect of heat stress on test-day milk yield was set to a THI of 80. All variance component estimates increased with parity; the largest increases were found for effects associated with heat stress. In particular, genetic variance associated with heat stress quadrupled from first to third parity, whereas permanent environmental variance only doubled. However, permanent environmental variance for heat stress was at least 10 times larger than genetic variance. Genetic correlations among parities for additive effects without heat stress considered ranged from 0.88 to 0.96. Genetic correlations among parities for additive effects of heat stress ranged from 0.08 to 0.22, and genetic correlations between effects regular and heat stress effects ranged from −0.21 to −0.33 for individual parities. Effect of heat stress on Thai Holstein crossbreds increased greatly with parity and was especially large after a THI of 80 for cows with a high percentage of Holstein genetics (≥93.7%). Individual sensitivity to heat stress was more environmental than genetic for Thai Holstein crossbreds.     

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on 1/6/2010orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.     

Genetic parameters for first and second lactation milk yields of Polish black and white cattle with random regression test-day models

Single- and two-trait random regression models were applied to estimate variance components of test-day records of milk, fat, and protein yields in the first and second lactation of Polish Black and White cattle. The model included fixed herd test-day effect, three covariates to describe lactation curve nested within age-season classes, and random regressions for additive genetic and permanent environmental effects. In two-parity models, each parity was treated as a separate trait. For milk and the two-parity model, heritabilities were in the range of 0.14 to 0.19 throughout first lactation and 0.10 to 0.16 throughout second lactation. For fat, heritabilities were within 0.11 to 0.16 and 0.11 to 0.22 throughout first and second lactations, respectively. For protein in the two-parity model, heritabilities were within 0.10 to 0.15 throughout most of first lactation and within 0.06 to 0.15 throughout the most of second lactation. For milk, genetic correlations between the first and second parities were 0.6 at the beginning of the lactation, rising to 0.9 in the middle, and 0.8 at the end of the lactation. For fat, the corresponding correlations were 0.6, 0.8, and 0.7, respectively, and for protein were 0.6, 0.8, and 0.8, respectively. Heritability estimates for all traits were flatter for the two-parity model. Relatively smooth genetic and permanent environmental variances with the two-parity model indicated that large swings of heritabilities could be artifacts of single-trait random regression models. High correlations between most of test day records across lactations suggested that a repeatability model could be considered as an alternative to a multiple-trait model to analyze multiple parities.     

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.     

Experience with a test-day model

The Canadian Test-Day Model is a 12-trait random regression animal model in which traits are milk, fat, and protein test-day yields, and somatic cell scores on test days within each of first three lactations. Test-day records from later lactations are not used. Random regressions (genetic and permanent environmental) were based on Wilmink's three parameter function that includes an intercept, regression on days in milk, and regression on an exponential function to the power -0.05 times days in milk. The model was applied to over 22 million test-day records of over 1.4 million cows in seven dairy breeds for cows first calving since 1988. A theoretical comparison of test-day model to 305-d complete lactation animal model is given. Each animal in an analysis receives 36 additive genetic solutions (12 traits by three regression coefficients), and these are combined to give one estimated breeding value (EBV) for each of milk, fat, and protein yields, average daily somatic cell score and milk yield persistency (for bulls only). Correlation of yield EBV with previous 305-d lactation model EBV for bulls was 0.97 and for cows was 0.93 (Holsteins). A question is whether EBV for yield traits for each lactation should be combined into one overall EBV, and if so, what method to combine them. Implementation required development of new methods for approximation of reliabilities of EBV, inclusion of cows without test day records in analysis, but which were still alive and had progeny with test-day records, adjustments for heterogeneous herd-test date variances, and international comparisons. Efforts to inform the dairy industry about changes in EBV due to the model and recovering information needed to explain changes in specific animals' EBV are significant challenges. The Canadian dairy industry will require a year or more to become comfortable with the test-day model and to realize the impact it could have on selection decisions.     

Adjustment for heterogeneous covariance due to herd milk yield by transformation of test-day random regressions

A method of accounting for differences in covariance components of test-day milk records was developed based on transformation of regressions for random effects. Preliminary analysis indicated that genetic and nongenetic covariance structures differed by herd milk yield. Differences were found for phenotypic covariances and also for genetic, permanent environmental, and herd-time covariances. Heritabilities for test-day milk yield tended to be lower at the end and especially at the start of lactation; they also were higher (maximum of ∼25%) for high-yield herds and lower (maximum of 15%) for low-yield herds. Permanent environmental variances were on average 10% lower in high-yield herds. Relative herd-time variances were ∼10% at start of lactation and then began to decrease regardless of herd yield; high-yield herds increased in midlactation followed by another decrease, and medium-yield herds increased at the end of lactation. Regressors for random regression effects were transformed to adjust for heterogeneity of test-day yield covariances. Some animal reranking occurred because of this transformation of genetic and permanent environmental effects. When genetic correlations between environments were allowed to differ from 1, some additional animal re-ranking occurred. Correlations of variances of genetic and permanent-environmental regression solutions within herd, test-day, and milking frequency class with class mean milk yields were reduced with adjustment for heterogeneous covariance. The method suggests a number of innovative solutions to issues related to heterogeneous covariance structures, such as adjusted estimates in multibreed evaluation.     

Short communication: Effect of preadjusting test-day yields for stage of pregnancy on variance component estimation in Canadian Ayrshires

Preadjustment of phenotypic records is an alternative to accounting for the effect of pregnancy within the genetic evaluation model. Variance components used in the Canadian Test-Day Model may need to be re-estimated after preadjusting for pregnancy. The objective of this study was to assess the effect of preadjusting test-day yields on variance components and estimated breeding values using a random regression test-day model in a random sample of Ayrshire cows. A random sample of 981 Canadian Ayrshire cows from 18 complete herds (average of 54.5 cows/herd) was analyzed. Two data sets were created using the same animals, one with unadjusted milk, fat, and protein yields, and one data set with test-day records adjusted for pregnancy effects. Pregnancy effect estimates from a previous study were used for additive preadjustment of records. Variance components were estimated using both data sets. Results were very similar between the 2 data sets for all estimated genetic parameters (heritabilities, genetic, and permanent environmental correlations). The relative squared differences were very small: 0.05% for heritabilities, 0.20% for genetic correlations, and 0.18% for permanent environmental correlations. Furthermore, paired Student's t-tests showed that the differences between the genetic parameters of data sets adjusted and unadjusted for pregnancy effect were not significantly different from 0. Results from this study show that preadjusting data for pregnancy did not yield changes in covariance component estimates, thus suggesting that preadjusting test-day records could be a feasible solution to account for pregnancy in the Canadian Test-Day Model without changing the current model. Estimated breeding values (EBV) were calculated for both data sets to observe the impact of preadjusting for pregnancy. Overall, the largest changes in EBV seen when preadjusting for pregnancy (compared with unadjusted records) occurred for nonpregnant elite cows, whose EBV declined. Preadjusting for pregnancy before genetic evaluations improves the estimation of breeding values by adding the negative impact of pregnancy back onto pregnant cow test-day records, causing an increase in their production EBV.     

Genetic parameters for a multiple-trait multiple-lactation random regression test-day model in Italian Holsteins

The objectives of this study were to estimate variance components for test-day milk, fat, and protein yields and average daily SCS in 3 subsets of Italian Holsteins using a multiple-trait, multiple-lactation random regression test-day animal model and to determine whether a genetic heterogeneous variance adjustment was necessary. Data were test-day yields of milk, fat, and protein and SCS (on a log₂ scale) from the first 3 lactations of Italian Holsteins collected from 1992 to 2002. The 3 subsets of data included 1) a random sample of Holsteins from all herds in Italy, 2) a random sample of Holsteins from herds using a minimum of 75% foreign sires, and 3) a random sample of Holsteins from herds using a maximum of 25% foreign sires. Estimations of variances and covariances for this model were achieved by Bayesian methods using the Gibbs sampler. Estimated 305-d genetic, permanent environmental, and residual variance was higher in herds using a minimum of 75% foreign sires compared with herds using a maximum of 25% foreign sires. Estimated average daily heritability of milk, fat, and protein yields did not differ among subsets. Heritability of SCS in the first lactation differed slightly among subsets and was estimated to be the highest in herds with a maximum of 25% foreign sire use (0.19 ± 0.01). Genetic correlations across lactations for milk, fat, and protein yields were similar among subsets. Genetic correlations across lactations for SCS were 0.03 to 0.08 higher in herds using a minimum of 75% or a maximum of 25% foreign sires, compared with herds randomly sampled from the entire population. Results indicate that adjustment for heterogeneous variance at the genetic level based on the percentage of foreign sire use should not be necessary with a multiple-trait random regression test-day animal model in Italy.     

Genetic parameters for test-day yield of milk, fat and protein in buffaloes estimated by random regression models

The test-day yields of milk, fat and protein were analysed from 1433 first lactations of buffaloes of the Murrah breed, daughters of 113 sires from 12 herds in the state of S?o Paulo, Brazil, born between 1985 and 2007. For the test-day yields, 10 monthly classes of lactation days were considered. The contemporary groups were defined as the herd-year-month of the test day. Random additive genetic, permanent environmental and residual effects were included in the model. The fixed effects considered were the contemporary group, number of milkings (1 or 2 milkings), linear and quadratic effects of the covariable cow age at calving and the mean lactation curve of the population (modelled by third-order Legendre orthogonal polynomials). The random additive genetic and permanent environmental effects were estimated by means of regression on third- to sixth-order Legendre orthogonal polynomials. The residual variances were modelled with a homogenous structure and various heterogeneous classes. According to the likelihood-ratio test, the best model for milk and fat production was that with four residual variance classes, while a third-order Legendre polynomial was best for the additive genetic effect for milk and fat yield, a fourth-order polynomial was best for the permanent environmental effect for milk production and a fifth-order polynomial was best for fat production. For protein yield, the best model was that with three residual variance classes and third- and fourth-order Legendre polynomials were best for the additive genetic and permanent environmental effects, respectively. The heritability estimates for the characteristics analysed were moderate, varying from 0·16±0·05 to 0·29±0·05 for milk yield, 0·20±0·05 to 0·30±0·08 for fat yield and 0·18±0·06 to 0·27±0·08 for protein yield. The estimates of the genetic correlations between the tests varied from 0·18±0·120 to 0·99±0·002; from 0·44±0·080 to 0·99±0·004; and from 0·41±0·080 to 0·99±0·004, for milk, fat and protein production, respectively, indicating that whatever the selection criterion used, indirect genetic gains can be expected throughout the lactation curve.     

Genotype x environment interaction for yield and somatic cell score with automatic and conventional milking systems

The objective of this study was to quantify genotype by environment interaction (G x E) between automatic milking systems (AMS) and conventional milking systems (CMS) for test-day milk, fat, and protein yield and for test-day somatic cell score (SCS) in The Netherlands. The G x E was studied in 2 ways: 1) between AMS farms and CMS farms in the same period and 2) within farms comparing the period before introduction of AMS with the period after introduction of AMS. For both sub-objectives, a separate data set was generated. Test-day records were used to be more flexible with respect to the introduction date of AMS. Multivariate, fixed regression, test-day sire models were used to estimate variance components. Genetic correlations between AMS farms and CMS farms in the same period were 0.93, >0.99, 0.98, and 0.79 for test-day milk yield, fat yield, protein yield, and SCS, respectively. Genetic correlations within farms between the period before and after introduction of AMS were lower for production traits and higher for SCS: 0.89, 0.91, 0.87, and >0.99, respectively, for test-day milk yield, fat yield, protein yield, and SCS. Heterogeneity of variance was observed between AMS and CMS in both data sets. Especially the residual variance increased with automatic milking. As a consequence, the heritability tended to be lower for automatic milking. It was concluded that effects of G x E are small between AMS and CMS. Therefore, AMS farms can select sires accurately based on national rankings.     

Comparison of random regression models with Legendre polynomials and linear splines for production traits and somatic cell score of Canadian Holstein cows

A random regression model with both random and fixed regressions fitted by Legendre polynomials of order 4 was compared with 3 alternative models fitting linear splines with 4, 5, or 6 knots. The effects common for all models were a herd-test-date effect, fixed regressions on days in milk (DIM) nested within region-age-season of calving class, and random regressions for additive genetic and permanent environmental effects. Data were test-day milk, fat and protein yields, and SCS recorded from 5 to 365 DIM during the first 3 lactations of Canadian Holstein cows. A random sample of 50 herds consisting of 96,756 test-day records was generated to estimate variance components within a Bayesian framework via Gibbs sampling. Two sets of genetic evaluations were subsequently carried out to investigate performance of the 4 models. Models were compared by graphical inspection of variance functions, goodness of fit, error of prediction of breeding values, and stability of estimated breeding values. Models with splines gave lower estimates of variances at extremes of lactations than the model with Legendre polynomials. Differences among models in goodness of fit measured by percentages of squared bias, correlations between predicted and observed records, and residual variances were small. The deviance information criterion favored the spline model with 6 knots. Smaller error of prediction and higher stability of estimated breeding values were achieved by using spline models with 5 and 6 knots compared with the model with Legendre polynomials. In general, the spline model with 6 knots had the best overall performance based upon the considered model comparison criteria.     

Environmental sensitivity for milk yield in Luxembourg and Tunisian Holsteins by herd management level

Milk production data of Luxembourg and Tunisian Holstein cows were analyzed using herd management (HM) level. Herds in each country were clustered into high, medium, and low HM levels based on solutions of herd-test-date and herd-year of calving effects from national evaluations. Data from both populations included 730,810 test-day (TD) milk yield records from 87,734 first-lactation cows. A multi-trait, random regression TD model was used to estimate (co)variance components for milk yield within and across country HM levels. Additive genetic and permanent environmental variances of TD milk yields varied with management level in Tunisia and Luxembourg. Additive variances were smaller across HM levels in Tunisia than in Luxembourg, whereas permanent environmental variances were larger in Tunisian HM levels. Highest heritability estimates of 305-d milk yield (0.41 and 0.21) were found in high HM levels, whereas lowest estimates (0.31 and 0.12, respectively) were associated with low HM levels in both countries. Genetic correlations among Luxembourg HM levels were >0.96, whereas those among Tunisian HM levels were below 0.80. Respective rank orders of sires ranged from 0.73 to 0.83 across Luxembourg environments and from 0.33 to 0.42 across Tunisian HM levels indicating high re-ranking of sires in Tunisia and only a scaling effect in Luxembourg. Across-country environment analysis showed that estimates of genetic variance in the high, medium, and low classes of Tunisian environments were 45, 69, and 81% lower, respectively, than the estimate found in the high Luxembourg HM level. Genetic correlations among 305-d milk yields in Tunisian and Luxembourg HM environments ranged from 0.39 to 0.79. The largest estimated genetic correlation was found between the medium Luxembourg and high Tunisian HM levels. Rank correlations for common sires’ estimated breeding values among HM environments were low and ranged from 0.19 to 0.39, implying the existence of genotype by environment interaction. These results indicate that daughters of superior sires in Luxembourg have their genetic expression for milk production limited under Tunisian environments. Milk production of cows in the medium and low Luxembourg environments were good predictors of that of their paternal half-sisters in the high Tunisian HM level. Breeding decisions in low-input Tunisian environment should utilize semen from sires with daughters in similar production environments rather than semen of bulls proven in higher management levels.     

Within-herd effects of age at test day and lactation stage on test-day yields

Variance ratios were estimated for random within-herd effects of age at test day and lactation stage, on test-day yield and somatic cell score to determine whether including these effects would improve the accuracy of estimation. Test-day data starting with 1990 calvings for the entire US Jersey population and Holsteins from California, Pennsylvania, Wisconsin, and Texas were analyzed. Test-day yields were adjusted for across-herd effects using solutions from a regional analysis. Estimates of the relative variance (fraction of total variance) due to within-herd age effects were small, indicating that regional adjustments for age were adequate. The relative variances for within-herd lactation stage were large enough to indicate that accuracy of genetic evaluations could be improved by including herd stage effects in the model for milk, fat, and protein, but not for somatic cell score. Because the within-herd lactation stage effect is assumed to be random, the effect is regressed toward the regional effects for small herds, but in large herds, lactation curves become herd specific. Model comparisons demonstrated the greater explanatory power of the model with a within-herd-stage effect as prediction error standard deviations were greater for the model without this effect. The benefit of the within-herd-stage effects was confirmed in a random regression model by comparing variance components from models with and without random within-herd regressions and through log-likelihood ratio tests.     

Short communication: best prediction of 305-day lactation yields with regional and seasonal effects

In the United States, lactation yields are calculated using best prediction (BP), a method in which test-day (TD) data are compared with breed- and parity-specific herd lactation curves that do not account for differences among regions of the country or seasons of calving. Complete data from 538,090 lactations of 348,123 Holstein cows with lactation lengths between 250 and 500 d, records made in a single herd, at least 5 reported TD, and twice-daily milking were extracted from the national dairy database and used to construct regional and seasonal lactation curves. Herds were assigned to 1 of 7 regions of the country, individual lactations were assigned to 3-mo seasons of calving, and lactation curves for milk, fat, and protein yields were estimated by parity group for regions, seasons, and seasons within regions. Multiplicative pre-adjustment factors (MF) also were computed. The resulting lactation curves and MF were tested on a validation data set of 891,806 lactations from 400,000 Holstein cows sampled at random from the national dairy database. Mature-equivalent milk, fat, and protein yields were calculated using the standard and adjusted curves and MF, and differences between 305-d mature-equivalent yields were tested for significance. Yields calculated using 50-d intervals from 50 to 250 d in milk (DIM) and using all TD to 500 DIM allowed comparisons of predictions for records in progress (RIP). Differences in mature-equivalent milk ranged from 0 to 51 kg and were slightly larger for first-parity than for later parity cows. Milk and components yields did not differ significantly in any case. Correlations of yields for 50-d intervals with those using all TD were similar across analyses. Yields for RIP were slightly more accurate when adjusted for regional and seasonal differences.     

Estimation of genetic parameters for milk fat depression in dairy cattle

The objective of this study was to apply reaction norm models to milk recording data to investigate genetic variation in and environmental sensitivity of susceptibility to milk fat depression (MFD). Data comprised 556,276 test-day records of 80,493 heifers in 1043 herds. Breeding values and genetic variances for fat percentage and fat yield were estimated by applying random regression models to average herd-test-day fat percentage. Genetic and permanent environmental correlations between fat yield expressed in different environments ranged, respectively, from 0.83 to 1.00 and from 0.29 to 1.00. Genetic and permanent environmental correlations between fat percentage expressed in different environments ranged, respectively, from 0.87 to 1.00 and from -0.05 to 0.99. Two traits were defined for MFD. The first trait reflected variation of milk fat percentage of animals within lactation after correction for year-season, herd-test-day, age-at-calving, and stage-of-lactation. This trait had an estimated heritability of about 5% and a genetic correlation between the fifth and 95th percentile of the data of 0.50. The second trait reflected the deviation of an animal's fat percentage on a test-day from its expected fat percentage based on fat percentage on the first test-day. This trait had an estimated heritability of about 4% and a genetic correlation between the fifth and 95th percentile of the data of 0.43. The correlation between estimated breeding values of sires for the 2 MFD traits was -0.3. Our results suggest that genetic variation in susceptibility to MFD is present and that selection for reduced susceptibility to MFD is possible.     

Joint international evaluation of Milking Shorthorn dairy cattle for production traits

Pedigree information and test-day records for the first 3 parities of Milking Shorthorn dairy cattle from 5 countries were analyzed. After editing, the data included 1,018,528 test-day records from 68,653 cows. A multiple-lactation random regression test-day model with Legendre polynomials of order 4 and a Bayesian method were used to estimate variance components for both single and multiple-countries. Fixed effects included herd-test-day class and regressions on DIM within age at calving-parity-season of calving. Random effects included animal genetic, permanent environmental, and residual effects. Average daily heritabilities from single country analyses ranged from 0.33 to 0.47 for milk yield and from 0.37 to 0.45 for protein yield across lactations and countries. Common sires (66) and their daughters were identified for creating a connected data set for simultaneous (co)variance component estimation of milk yield across all 5 countries. Between-country genetic correlations were low, with values from 0.08 to 0.46 and standard deviations from 0.08 to 0.12. Estimated breeding values for milk were generated for each animal using the same test-day animal model. Correlations among country estimated breeding values were higher than genetic correlations. Top 100 bull lists were generated on the scale of each country, and genetic progress was assessed. Future evaluation with increased genetic ties among countries may facilitate international comparison of Milking Shorthorns.