Bayesian analysis of lactation curves of Holstein-Friesian cattle using a nonlinear model

A Bayesian procedure was developed for fitting Wood's incomplete Gamma function to test-day milk records of Spanish Holstein Friesian cattle. Each parameter of Wood's function was considered as a dependent variable in a submodel that accounted for systematic effects and genetic relationships among animals. Marginal posterior distributions of model parameters were obtained using Gibbs sampling. Variables of economic interest, such as 305-d yield, persistency, peak yield, and days in milk at peak day were predicted as functions of Wood's function curve parameters. Heritability estimates were 0.26, 0.32, and 0.19 for parameters of Wood's function and 0.26, 0.14, 0.26, and 0.05 for 305-d yield, persistency, peak yield, and days in milk at peak yield. These estimates indicate that it is possible to modify the shape of the lactation curve through genetic selection. Genetic correlations between parameters of Wood's curve and the aforementioned functions of these parameters suggest that selection for 305-d milk yield would result in higher and later peak yield, but only a slight improvement in persistency is expected.     

Effects of continuous milking during the dry period or once daily milking in the first 4 weeks of lactation on metabolism and productivity of dairy cows

The objective was to compare the effects of 3 management systems in high-yielding dairy cows on metabolic profiles and milk production. Thirty-six multiparous Brown Swiss cows were randomly assigned to 1 of 3 treatment groups (n = 12 cows/group): the control (C) group, in which cows were dried off 56 d before calving and milked twice daily throughout next lactation (305 d); the once daily milking (ODM) group, in which cows were dried off 56 d before calving and milked once daily for the first 4 wk of lactation and twice daily for the remaining lactation; and the continuous milking (CM) group, in which cows were milked twice daily until calving and also during the subsequent lactation. Serum glucose concentrations decreased between wk 1 and 4 exclusively in C cows. Serum concentrations of NEFA and BHBA in the first 4 wk of lactation were highest in C cows compared with ODM and CM cows. Decreased backfat thickness during early lactation and reduction of body condition score were markedly more pronounced in C cows compared with ODM and CM cows. Mean lactational milk yield of C cows [11,310 ± 601 kg of energy-corrected milk (ECM)/305 d] was approximately 16% higher compared with ODM cows (9,531 ± 477 kg of ECM/305 d) and CM cows (9,447 ± 310 kg of ECM/305 d). The lactation curve of CM cows compared with C cows was characterized by a similar time of peak yield (wk 3), a reduced peak yield, and no obvious differences in persistency. Mean percentage of milk protein was significantly higher for CM cows (3.91%) compared with C cows (3.52%). In contrast, once daily milking was accompanied by a reduced and significantly delayed peak yield (wk 8) compared with the control treatment, whereas persistency was better and milk protein (3.79%) was higher in ODM cows than in C cows. In conclusion, continuous milking and once daily milking, targeting the interval before or after calving, respectively, substantially reduced the metabolic challenge of fresh cows and improved milk protein percentage. Continuous milking and once daily milking increased milk protein percentage markedly; furthermore, once daily milking during the first 4 wk of lactation improved the lactation curve.     

Genetic relationships between persistency and reproductive performance in first-lactation Canadian holsteins

The main objective of this study was to estimate genetic relationships between lactation persistency and reproductive performance in first lactation. Relationships with day in milk at peak milk yield and estimated 305-d milk yield were also investigated. The data set contained 33,312 first-lactation Canadian Holsteins with first-parity reproductive, persistency, and productive information. Reproductive performance traits included age at first insemination, nonreturn rate at 56 d after first insemination as a virgin heifer and as a first-lactation cow, calving difficulty at first calving and calving interval between first and second calving. Lactation persistency was defined as the Wilmink b parameter for milk yield and was calculated by fitting lactation curves to test day records using a multiple-trait prediction procedure. An 8-trait genetic analysis was performed using the Variance Component Estimation package (VCE 5) via Gibbs sampling to estimate genetic parameters for all traits. Heritabilities of persistency, day in milk at peak milk yield and estimated 305-d milk yield were 0.18, 0.09 and 0.45, respectively. Heritabilities of reproduction were low and ranged from 0.03 to 0.19. The highest heritability was for age at first insemination. Heifer reproductive traits were lowly genetically correlated, whereas cow reproductive traits were moderately correlated. Heifers younger than average when first inseminated and/or conceived successfully at first insemination tended to have a more persistent first lactation. First lactation was more persistent if heifers had difficulty calving (r(g) = 0.43), or conceived successfully at first insemination in first lactation (r(g) = 0.32) or had a longer interval between first and second calving (r(g) = 0.17). Estimates of genetic correlations of reproductive performance with estimated 305-d milk yield were different in magnitude, but similar in sign to those with persistency (0.02 to 0.51).     

Nonlinear modeling to describe the pattern of 15 milk protein and nonprotein compounds over lactation in dairy cows

The protein profile of milk includes several caseins, whey proteins, and nonprotein nitrogen compounds, which influence milk's value for human nutrition and its cheesemaking properties for the dairy industry. To fill in the gap in current knowledge of the patterns of these individual nitrogenous compounds throughout lactation, we tested the ability of a parametric nonlinear lactation model to describe the pattern of each N compound expressed qualitatively (as % of total milk N), quantitatively (in g/L milk), and as daily yield (in g/d). The lactation model was tested on a data set of detailed milk nitrogenous compound profiles (15 fractions—12 protein traits and 3 nonproteins—for each expression mode: 45 traits) obtained from 1,342 cows reared in 41 multibreed herds. Our model was a modified version of Wilmink's model, often used for describing milk yield during lactation because of its reliability and ease of parameter interpretation from a biological point of view. We allowed the sign of the persistency coefficient (parameter c) that explained the variation in the long-term milk component (parameter a) to be positive or negative. We also allowed the short-term milk component (parameter b) to be positive or negative, and we estimated a specific speed of adaptation parameter (parameter k) for each trait rather than assumed a value a priori, as in the original model (k = 0.05). These 4 parameters were included in a nonlinear mixed model with cow breed and parity order as fixed effects, and herd-date as random. Combinations of the positive and negative signs of the b and c parameters allowed us to identify 4 differently shaped lactation curves, all found among the patterns exhibited by the nitrogenous fractions as follows: the “zenith” curve (with a maximum peak; for milk yield and 10 other N traits), the “nadir” curve (with a minimum point; for 20 traits, including almost all those expressed in g/L of milk), the “downward” curve (continuously decreasing; for 14 traits, including almost all those in g/d), and the “upward” curve (continuously increasing; only for κ-casein, in % N). Direct estimation of the k parameters specific to each trait showed the large variability in the adaptation speed of fresh cows and greatly increased the model's flexibility. The results indicated that nonlinear parametric mathematical models can effectively describe the different and complex patterns exhibited by individual nitrogenous fractions during lactation; therefore, they could be useful tools for interpreting milk composition variations during lactation.     

Evidence that leptin genotype is associated with fertility, growth, and milk production in Holstein cows

This study investigated associations between single nucleotide polymorphisms (SNP) in the leptin, leptin receptor, and neuropeptide Y (NPY) genes with growth, milk production, and fertility traits. Holstein Friesian heifers from 19 UK dairy farms were recruited at birth, providing an initial population of 509. Animals were monitored until they either reached the end of their second lactation or were culled. Size (weight, height, length, girth) was measured at 1, 6, and 15 mo to assess growth traits. Heifer fertility was assessed by recording age at service, age at conception, age at first calving, and number of services. Cow fertility was assessed by recording days from calving to service and conception, services per conception, percentage of animals in calf at 100 d after calving and the calving interval in each lactation. Milk production was recorded as days in milk, total milk per lactation, 305-d yield, milk per day, and peak yield. Mixed model analyses revealed that leptin SNP were associated with early skeletal growth (height, A1457G; length, A59V), fertility (UASMS1, UASMS2, A1457G, A59V) and milk production (A59V). The leptin receptor SNP (T945M) was only weakly associated with milk per day and days to first service. The NPY SNP (NPY1) was associated with the prevalence of the animal being in calf 100 d after calving and 305-d milk yield in the first lactation. The association of leptin SNP with fertility traits in heifers, in addition to lactating cows, suggests that some effects on fertility are direct and not necessarily mediated via altered tissue mobilization. In accord with this, other work has shown that leptin can affect oocyte quality and early embryo development. These results support the use of leptin SNP to inform marker-assisted selection in dairy cows.     

Cow-level factors associated with subclinical hypocalcemia at calving in multiparous Jersey cows

The objective of our study was to identify cow-level factors associated with subclinical hypocalcemia at calving (SCH) in multiparous Jersey cows. A total of 598 Jersey and 218 Jersey × Holstein crossbreed cows from 2 commercial dairy herds were enrolled in a retrospective cohort study. Blood samples to determine total Ca concentration were collected from the coccygeal vessels at 3 h 19 min (±2 h 33 min) after calving. We used 2 serum Ca concentration thresholds to define SCH: <2.00 mmol/L (SCH-2.00) and <2.12 mmol/L (SCH-2.12). We evaluated the association of cow-level factors with SCH with multivariable Poisson regression models. Variables evaluated for association with SCH were herd; parity (2, 3, and ≥4); breed; previous lactation length and 305-d mature-equivalent milk yield; previous lactation first test milk yield and last test somatic cell count; lengths of calving interval, gestation, dry, and close-up periods; body condition and locomotion scores at calving; calving ease; and calf sex for singletons. We categorized continuous variables into quartiles (≤25th percentile, interquartile range and ≥75th percentile). The prevalence of SCH among Jersey cows was 40 (SCH-2.00) and 64% (SCH-2.12). Jersey cows of higher parity had greater risk of SCH-2.00 and SCH-2.12. The risk of SCH-2.12 was higher after birthing male calves. We also found a tendency for previous lactation length and previous lactation 305-d mature-equivalent milk yield effect to affect risk of SCH-2.12. The risk of SCH-2.12 was lower for cows that had a previous lactation length shorter than the 25th percentile compared with cows that had a previous lactation length within the interquartile range. The risk of SCH-2.12 was higher for cows that had a previous lactation 305-d mature-equivalent milk yield below the 25th percentile compared with cows that had a previous lactation 305-d mature-equivalent milk yield above the 75th percentile. Also, Jersey × Holstein crossbreed was associated with increased risk of SCH-2.00. In the multivariable analysis, we observed no association between SCH and previous lactation first test milk yield; last test somatic cell count; lengths of calving interval, gestation, dry, and close-up periods; body condition and locomotion scores at calving; and calving ease. Our study identified parity, breed, calf sex, previous lactation length, and previous lactation 305-d mature-equivalent milk yield as cow-level factors associated with SCH in multiparous Jersey cows.     

Predicting Milk Yield of Holstein Cows from 306 to 395 Days in Milk

Prediction equations were determined to estimate daily milk yield from 306 to 395 d in milk for forecasting herd milk sales from Holstein cows in lactation >305 d. Data were test day milk weights for 65,322 primiparous and 119,220 pluriparous lactations of > 305 d from the Southern US. A forecast model was developed using same lactation 305 d milk yield (in classes of 500 kg increments) that gave similar predicted daily yields as models utilizing last sample milk weight information. This model has the advantage of early forecasting of later milk using projected 305-d yields.Reduced forecast models ignoring days pregnant, yield class, or both accounted for 95, 68, and 59%, and 91, 67, and 56% as much variation in daily milk as the full model for the primiparous and pluriparous cows. Percentage of 305-d milk yielded in mo 11, 12, and 13, depending on 305-d yield class, ranged from 7.1 to 7.0%, 6.2 to 6.0%, and 5.4 to 5.0%, and 5.4 to 5.0%, 4.3 to 3.9%, and 3.3 to 3.0% for first parity and pluriparous cows calving in winter and 125 d open. Cows not calving in winter or with more than 125 d open yielded more milk in extended lactation. These percentages are larger than generally assumed in studies of days open, thus indicating that cost of days open may have been overestimated.     

Relationships among milk yield,body condition,cow weight,and reproduction in spring-calved Holstein-Friesians

Relationships among milk production, body condition score (BCS), body weight (BW), and reproduction were studied using logistic regression on data from 6433 spring-calving Holstein-Friesian dairy cows in 74 commercial herds. Multivariate models were adjusted for herd, breeding value for milk yield, proportion of Holstein-Friesian genes, lactation number, calving period, and degree of calving assistance. Significant associations between reproductive measures and components of energy balance were identified. Higher 200-d milk protein content and higher protein-to-fat ratio at start of breeding were associated with increased likelihood of submission for breeding in the first 21 d of the breeding season (SR21). High 100-d cumulative milk yield as a proportion of estimated 305-d milk yield (low persistency) was associated with a lower likelihood of pregnancy to first service (PREG1), whereas cows reaching peak milk yields earlier tended to have higher PREG1. Cows that reached nadir milk protein content relatively late in lactation had lower PREG1. Milk yield at first service and 305-d milk protein content were positively associated with the likelihood of pregnancy after 42 d of breeding (PR42). Higher 305-d milk lactose content was associated with increased PREG1 and PR42. Mean BCS at 60 to 100 d of lactation was positively associated with both SR21 and PR42, whereas nadir BCS was positively associated with PREG1. Cows with precalving BCS > 3.0 that also lost > 0.5 BCS unit by first service had lower PR42. More BW gain for 90 d after start of breeding was associated with higher SR21 and PREG1; more BW gain for 90 d after first service was associated with higher PR42. Milk protein and lactose content, BCS, and BW changes are important tools to identify cows at risk of poor reproduction.     

Short communication: do Holstein lactations of varied lengths have different characteristics?

The calving intervals of Holstein cows in many countries have been increasing in length over recent years. Many reports of lactation characteristics refer only to a standard 305-d lactation. This paper attempts to describe the characteristics of lactations of different lengths using a large database of nationally recorded data. Three different lactation lengths were studied: 305 d (the traditional annual calving interval), 370 d (the current UK national average), and 440 d (equivalent to an 18-mo calving interval). Test-day milk yield; fat, protein, and lactose composition; and somatic cell counts were analyzed. Characteristics of each of 29,838 lactations were calculated using a biological model of lactation and were analyzed to identify which factors affected them. Maximum secretion potential was the only trait not affected by lactation length. Day of peak yield, peak yield, persistency, and total milk yield all increased with lactation length, whereas relative cell death rate and the rate of milk increase in early lactation both decreased as lactation length increased. For most weeks of lactation, the 3 lactation-length groups differed from each other in most traits, following the order of their lengths; for example, milk yield was always higher for lactations of 440 d and lower for those of 305 d. The post-peak trends in all traits were found to continue in longer lactations. Thus, daily milk production and lactose percentage continued to decrease as lactations became longer, whereas fat percentage, protein percentage, and somatic cell count continued to increase. Pregnancy was found to affect all traits, leading to an accentuation of these trends in late lactation. However, the effect of pregnancy depended on the yield at about the fourth month of pregnancy. Lactations of 305, 370, and 440 d all had different characteristics and not solely due to increasing length.     

Agreement of predicted 305-day milk yields relative to actual 305-day milk weight yields

The objectives of this study were to compare the multiple trait prediction (MTP) model estimate of 305-d lactation yield with the 305-d daily milk yield data from on-farm automated meters and software and to examine the accuracy of electronic identification (ID). Twenty-four-hour milk and component yields are calculated by using milk weights and samples collected 8 to 10 times/yr by Dairy Herd Improvement (DHI) organizations. Daily milk weights were collected from cows on 20 Canadian farms that used parlor milking systems with electronic ID and that were enrolled in a regular DHI program. A total of 10,175 DHI test days from 1,103 cows with complete 305-d lactation yields were entered into the MTP model, and lactation yields were predicted. Test days were grouped into first, second, and third and greater lactations and within each lactation group, days in milk were categorized in 3 stages (5 to 60, 61 to 120, and 120 to 305 d in milk) for a total of 9 classes. Agreement analysis was used to compare the 305-d sum of daily milk to the MTP 305-d lactation yield predictions by using inputs from test days throughout the lactations. Results indicated that the MTP model overestimated lactation yields across all parity groups, ranging from 310 to 1,552 kg in parity 1, 640 to 2,000 kg in parity 2, and 567 to 1,476 kg in parity 3 and greater. A preliminary examination of electronic ID accuracy was conducted on 4 farms. Two electronic ID systems were examined for cow ID accuracy by verifying the ID number appearing in the parlor with the corresponding ear tag number. There were no ID errors on 3 of 4 farms tested and only a very small number of errors (3/80) on the fourth farm, indicating that the electronic ID systems used in milking parlors identify cows accurately.     

Genotype x environment interaction for milk yield in Holsteins using Luxembourg and Tunisian populations

Test-day (TD) milk yield records of first-lactation Holstein cows in Luxembourg and Tunisia were analyzed using within-and between-country random regression TD models. Edited data used for within-country analysis included 661,453 and 281,913 TD records in Luxembourg and Tunisia, respectively. The joint data included 730,810 TD records of 87,734 cows and 231 common sires. Both data sets covered calving years 1995 to 2006. Fourth-order Legendre polynomials for random effects and a Gibbs sampling method were used to estimate variance components of lactation curve parameters in separate and joint analyses. Genetic variances of the first 3 coefficients from Luxembourg data were 46 to 69% larger than corresponding estimates from the Tunisian data. Inversely, the Tunisian permanent environment variances for the same coefficients were 52 to 65% larger than the Luxembourg ones. Posterior mean heritabilities of 305-d milk yield and persistency, defined as estimated breeding values (EBV) at 280 days in milk-EBV at 80 days in milk, from between-country analysis were 0.42 and 0.12 and 0.19 and 0.08 in Luxembourg and Tunisia, respectively. Heritability estimates for the same traits from within-country analyses, mainly from the Tunisian data, were lower than those from the joint analysis. Genetic correlations for 305-d milk yield and persistency between countries were 0.60 and 0.36. Product moment and rank correlations between EBV of common sires for 305-d milk yield and persistency from within-country analyses were 0.38 and 0.41 and 0.27 and 0.26, respectively. Differences between genetic variances found in both countries reflect different milk production levels. Moreover, low genetic and rank correlations suggest different ranking of sires in the 2 environments, which implies the existence of a genotype × environment interaction for milk yield in Holsteins.     

Phenotypic effects of calving ease on the subsequent fertility and milk production of dam and calf in UK Holstein-Friesian heifers

The effect of calving ease on the fertility and production performance of both dam and calf was studied in approximately 50,000 and 10,000 UK Holstein-Friesian heifers and heifer calves, respectively. The first objective of this study was to estimate the effect of a difficult calving on the subsequent first-lactation milk production by estimating lactation curves using cubic splines. This methodology allows the estimation of daily milk, protein, and fat yields following calvings of differing degrees of difficulty. Losses in milk yield after a difficult calving have been quantified previously; however, estimates are generally restricted to the accumulated yields at specific days in lactation. By fitting cubic splines, gaps (in which the shape of the lactation curve can be merely guessed) between estimations were avoided. The second objective of this study was to estimate the effect of a difficult birth on the subsequent performance of the calf as an adult animal. Even though the calving process is known to involve cooperation between dam and calf, the effect of a difficult calving has, until now, only been estimated for the subsequent performance of the dam. Addressing the effects of a difficult birth on the adult calf strengthens the importance of calving ease as a selection trait because it suggests that the benefit of genetic improvement may currently be underestimated. The effect of calving ease on the subsequent reproductive performance of dam and calf was analyzed using linear regression and with calving ease score fitted as a fixed effect. Dams with veterinary-assisted calvings required 0.7 more services to conception and 8 more days to first service and experienced a 28-d longer calving interval in first lactation compared with dams that were not assisted at calving. Effects of calving ease on the reproductive performance of the adult calf in first lactation were not detected. Losses in milk yield of the dam were significant between d 9 to 90 in milk subsequent to a veterinary-assisted calving, creating a loss of approximately 2 kg of milk per day, compared with a nonassisted calving. Calves being born with difficulties showed a significant reduction in milk yield in first lactation, demonstrating the lifelong effect of a difficult birth. Compared with nonassisted calves, veterinary-assisted calves showed a loss of 710 kg in accumulated 305-d milk yield, which was significant from 129 to 261 d in milk. This suggests that from birth to production, physiological effects of a bad calving are not negated. Results furthermore suggest a beneficial effect of farmer assistance at calving on the milk yield of both dam and calf, when moderate difficulties occurred.     

Effect of dry period length on milk yield over multiple lactations

Shortening or omitting the dry period (DP) can improve the energy balance of dairy cows in early lactation through a decrease in milk yield after calving. Little is known about the effect of a short or no DP on milk yield over multiple lactations. Our objectives were (1) to assess the effect of DP length over multiple lactations on milk yield, and (2) to assess if the prediction of milk yield in response to DP length could be improved by including individual cow characteristics before calving. Lactation data (2007 to 2015) of 16 Dutch dairy farms that apply no or short DP were used to compute cumulative milk yield in the 60 d before calving (additional yield) and in the 305 d after calving (305-d yield), and the mean daily yield over the interval from 60 d before calving to 60 d before next calving (effective lactation yield). The DP categories were no (0 to 2 wk), short (3 to 5 wk), standard (6 to 8 wk), and long (9 to 12 wk). The effect of current DP and previous DP on yields was analyzed with mixed models (n = 1,420 lactations). The highest effective lactation yield of fat- and protein-corrected milk (FPCM) was observed for cows with a standard current DP (27.6 kg per day); a daily decrease was observed of 0.6 kg for a long DP, 1.0 kg for a short DP, and 2.0 kg for no DP. Previous DP did not significantly affect the effective lactation yield. Thus, cows can be managed with short or no DP over consecutive lactations without a change in quantity of milk losses. Cows that received no DP for consecutive lactations had a lower additional yield before calving (?172 kg of FPCM), but a higher 305-d yield (+560 kg of FPCM), compared with cows that received no DP for the first time. This could lessen the improvement of the energy balance in early lactation when no DP is applied a second time compared with the first time. For the second objective, a basic model was explored to predict effective lactation yield based on parity, DP length, and first-parity 305-d yield (n = 2,866 lactations). The basic model was subsequently extended with data about recent yield, days open, and somatic cell count. Extending the model reduced the error of individual predictions by only 6%. Therefore, the basic model seems sufficient to predict the effect of DP length on effective lactation yield. Other individual cow characteristics can still be relevant, however, to make a practical and tailored decision about DP length.     

Genetic analysis of milk production traits of polish black and white cattle using large-scale random regression test-day models

Genetic parameters for milk, fat, and protein yield and persistency in the first 3 lactations of Polish Black and White cattle were estimated. A multiple-lactation model was applied with random herd-test-day effect, fixed regressions for herd-year and age-season of calving, and random regressions for the additive genetic and permanent environmental effects. Three data sets with slightly different edits on minimal number of days in milk and the size of herd-year class were used. Each subset included more than 0.5 million test-day records and more than 58,000 cows. Estimates of covariance components and genetic parameters for each trait were obtained by Bayesian methods using the Gibbs sampler. Due to the large size and a good structure of the data, no differences in estimates were found when additional criteria for record selection were applied. More than 95% of the genetic variance for all traits and lactations was explained by the first 2 principal components, which were associated with the mean yield and lactation persistency. Heritabilities of 305-d milk yield in the first 3 lactations (0.18, 0.16, 0.17) were lower than those for fat (0.12, 0.11, 0.12) and protein (0.13, 0.14, 0.15). Estimates of daily heritabilities increased in general with days in milk for all traits and lactations, with no apparent abnormalities at the beginning or end of lactation. Genetic correlations between yields in different lactations ranged from 0.74 (fat yield in lactations 1 and 3) to 0.89 (milk yield in lactations 2 and 3). Persistency of lactation was defined as the linear regression coefficient of the lactation curve. Heritability of persistency increased with lactation number for all traits and genetic correlations between persistency in different lactations were smaller than those for 305d yield. Persistency was not genetically correlated with the total yield in lactation.     

Effect of lactation number, year, and season of initiation of lactation on milk yield of cows hormonally induced into lactation and treated with recombinant bovine somatotropin

Records representing data from 1,500 barren Holstein cows over an 8-yr period from a large commercial dairy farm in northern Mexico were analyzed to determine the effects of lactation number and season and year of initiation of lactation on milk production of cows induced hormonally into lactation and treated with recombinant bovine somatotropin (rbST) throughout lactation. Peak and 305-d milk yields were also assessed as predictors of total milk yield in cows induced into lactation. A significant quadratic relationship was found between 305-d milk yield and number of lactation [7,607 ± 145 and 9,548 ± 181 kg for first- and ≥6-lactation cows, respectively; mean ± standard error of the mean (SEM)] with the highest production occurring in the fifth lactation. Total milk yields of cows with ≤2 lactations were approximately 4,500 kg less than milk yields of adult cows (the overall average ± standard milk yield was 13,544 ± 5,491 kg per lactation and the average lactation length was 454 ± 154 d). Moreover, 305-d milk production was depressed in cows induced into lactation in spring (8,804 ± 153 kg; mean ± SEM) and summer (8,724 ± 163 kg) than in fall (9,079 ± 151 kg) and winter (9,085 ± 143 kg). Partial regression coefficients for 305-d milk yield and peak milk yield indicated an increment of 157 kg of milk per lactation per 1-kg increase in peak milk yield (r² = 0.69). Neither peak milk yield (r² = 0.18) nor 305-d milk yield (r² = 0.29) was accurate for predicting total milk yield per lactation. Year, parity, and season effects had significant influence on milk yield of cows induced into lactation and treated with rbST throughout lactation, and peak milk yield can assist in the prediction of 305-d milk yield but not total milk yield. This study also showed that hormonal induction of lactation in barren high-yielding cows is a reliable, practical, and affordable technique in countries where rbST treatment and prolonged steroid administration of dairy cows are legally permitted.     

Modeling extended lactations of dairy cows

Nine mathematical models were compared for their ability to predict daily milk yields (n = 294,986) in standard 305-d and extended lactations of dairy cows of Costa Rica. Lactations were classified by parity (first and later), lactation length (9 to 10, 11 to 12, 13 to 14, 15 to 16, and 16 to 17 mo), and calving to conception interval (1 to 2, 3 to 4, 5 to 6, 7 to 8, and 9 to 10 mo). Of the nine models, the diphasic model and lactation persistency model resulted in the best goodness of fit as measured by adjusted coefficient of determination, residual standard deviation, and Durbin-Watson coefficient. All other models showed less accuracy and positively correlated residuals. In extended lactations, models were also fitted using only test-day records before 305 d, which resulted in a different ranking. The diphasic model showed the best prediction of milk yield in standard and extended lactations. We concluded that the diphasic model provided accurate estimates of milk yield for standard and extended lactations. Interpretation of parameters deserves further attention because of the large variation observed. As expected, the calving to conception interval was found to have a negative effect on milk yield for cows with a standard lactation length. In extended lactations, these negative effects of pregnancy on milk yield were not observed.     

Genetic parameters of feed intake, production, body weight, body condition score, and selected type traits of Holstein cows in commercial tie-stall barns

The objectives of this study were to determine the feasibility of measuring feed intake in commercial tie-stall dairies and infer genetic parameters of feed intake, yield, somatic cell score, milk urea nitrogen, body weight (BW), body condition score (BCS), and linear type traits of Holstein cows. Feed intake, BW, and BCS were measured on 970 cows in 11 Pennsylvania tie-stall herds. Historical test-day data from these cows and 739 herdmates who were contemporaries during earlier lactations were also included. Feed intake was measured by researchers once per month over a 24-h period within 7 d of 6 consecutive Dairy Herd Information test days. Feed samples from each farm were collected monthly on the same day that feed intake was measured and were used to calculate intakes of dry matter, crude protein, and net energy of lactation. Test-day records were analyzed with multiple-trait animal models, and 305-d fat-corrected milk yield, dry matter intake, crude protein intake, net energy of lactation intake, average BW, and average BCS were derived from the test-day models. The 305-d traits were also analyzed with multiple-trait animal models that included a prediction of 40-wk dry matter intake derived from National Research Council equations. Heritability estimates for 305-d intake of dry matter, crude protein, and net energy of lactation ranged from 0.15 to 0.18. Genetic correlations of predicted dry matter intake with 305-d dry matter, crude protein, and net energy of lactation intake were 0.84, 0.90, and 0.94, respectively. Genetic correlations among the 3 intake traits and fat-corrected milk yield, BW, and stature were moderate to high (0.52 to 0.63). Results indicate that feed intake measured in commercial tie-stalls once per month has sufficient accuracy to enable genetic research. High-producing and larger cows were genetically inclined to have higher feed intake. The genetic correlation between observed and predicted intakes was less than unity, indicating potential variation in feed efficiency.     

Genome-wide association studies for genetic effects that change during lactation in dairy cattle

Genetic effects on milk production traits in dairy cattle might change during lactation. However, most genome-wide association studies (GWAS) for milk production traits assume that genetic effects are constant during lactation. This assumption might lead to missing these quantitative trait loci (QTL) whose effects change during lactation. This study aimed to screen the whole genome specifically for QTL whose effects change during lactation. For this purpose, 4 different GWAS approaches were performed using test-day milk protein content records: (1) separate GWAS for specific lactation stages, (2) GWAS for estimated Wilmink lactation curve parameters, (3) a GWAS using a repeatability model where SNP effects are assumed constant during lactation, and (4) a GWAS for genotype by lactation stage interaction using a repeatability model and accounting for changing genetic effects during lactation. Separate GWAS for specific lactation stages suggested that the detection power greatly differs between lactation stages and that genetic effects of some QTL change during lactation. The GWAS for estimated Wilmink lactation curve parameters detected many chromosomal regions for Wilmink parameter a (protein content level), whereas 2 regions for Wilmink parameter b (decrease in protein content toward nadir) and no regions for Wilmink parameter c (increase in protein content after nadir) were detected. Twenty chromosomal regions were detected with effects on milk protein content; however, there was no evidence that their effects changed during lactation. For 5 chromosomal regions located on chromosomes 3, 9, 10, 14, and 27, significant evidence was observed for a genotype by lactation stage interaction and thus their effects on milk protein content changed during lactation. Three of these 5 regions were only identified using a GWAS for genotype by lactation stage interaction. Our study demonstrated that GWAS for genotype by lactation stage interaction offers new possibilities to identify QTL involved in milk protein content. The performed approaches can be applied to other milk production traits. Identification of QTL whose genetic effects change during lactation will help elucidate the genetic and biological background of milk production.     

Estimation of genetic parameters for test-day records of dairy traits in a seasonal calving system

Variance components for test-day milk yield were estimated for primiparous animals in a seasonal calving system where 80% of calvings occur within a 5-mo period in the spring. The objective was to investigate if the variance components of milk production were affected by seasonality via month in milk and test month (TM). These effects were, therefore, fitted for both the permanent environment effect and additive genetic effect. Estimates of heritability (0.14-0.27) were found to be lowest during early and late lactation for all calving months. The peak in heritability (0.22-0.27) occurred later in lactation for animals that calved toward the end of the spring calving season. Genetic variance of test day milk yield for all calving months was found to be highest at the beginning of the lactation (2.23-3.36kg(2)), with a plateau toward the middle of the lactation. Genetic variance was found to be highest (2.98-3.36kg(2)) for animals calving early in the season. Genetic correlations between corresponding stages of lactation were strongest (0.99-1.00) between consecutive calving months. Genetic correlations were slightly weaker when intervals between calving months increased; however, they remained above 0.96 in all cases. Akaike's information criterion values from models both inclusive and exclusive of TM, suggest that the model exclusive of TM is the preferred model. Estimated breeding values of lactation milk yield from bulls with 20 or more daughters, predicted for a full 305-d lactation, were used to compare the model with the standard test-day model (i.e., exclusive of TM), based on DIM. Correlations were higher than 0.995 for milk yield (and higher than 0.930 for persistency) between models inclusive and exclusive of TM, suggesting that given the straightforward approach taken in the current study, an apparent benefit for including seasonality in the evaluation of test-day milk yield was not found; however, there may be benefits of including it in the estimation of persistency.     

Prediction of Average Daughter Performance from Sire Proofs for Use in Linear Programming Sire Profit Models

When linear programming is used to rank sires on daughter profits, coefficients of input are required for representative daughters of each sire. Within herd-year-season regressions of daughter's actual 305-day milk, fat percent, days open, 2-min milk yield adjusted for total yield, total milking time, and body weight at birth and first calving on sire's proofs for production, conformation, milking speed, and nonreturn rate were computed from progeny data from 71 Canadian Holstein-Friesian artificial insemination sires, to provide these coefficients. Coefficients of determination for prediction equations were .23 to .71. A production function was fitted for milk production and days open. First lactation 305-day milk was predicted from sire's milk rating and linear and quadratic terms of size ratings. Fat percent proof and its squared term predicted daughter's fat percent. Days open was predicted by milk and fat proofs. Milk yield in the first 2-min of milking, adjusted for milking yield, was predicted by proofs for milking speed and the quadratic term for fat percent. Total milking time was obtained from 2-min yield by regression (?.3837) of total milking time on 2-min yield. Daughter's body weight at first calving was predicted from sire's size proof, as was birth weight. Milk yield was curvilinearly associated with days open. Incidence of mastitis was derived indirectly from regression (.00172) on sire's milk proof for milk.