Comparison of non-linear models to describe the lactation curves for milk yield and composition in buffaloes (Bubalus bubalis)

In order to describe the lactation curves of milk yield (MY) and composition in buffaloes, seven non-linear mathematical equations (Wood, Dhanoa, Sikka, Nelder, Brody, Dijkstra and Rook) were used. Data were 116 117 test-day records for MY, fat (FP) and protein (PP) percentages of milk from the first three lactations of buffaloes which were collected from 893 herds in the period from 1992 to 2012 by the Animal Breeding Center of Iran. Each model was fitted to monthly production records of dairy buffaloes using the NLIN and MODEL procedures in SAS and the parameters were estimated. The models were tested for goodness of fit using adjusted coefficient of determination root means square error (RMSE), Durbin–Watson statistic and Akaike’s information criterion (AIC). The Dijkstra model provided the best fit of MY and PP of milk for the first three parities of buffaloes due to the lower values of RMSE and AIC than other models. For the first-parity buffaloes, Sikka and Brody models provided the best fit of FP, but for the second- and third-parity buffaloes, Sikka model and Brody equation provided the best fit of lactation curve for FP, respectively. The results of this study showed that the Wood and Dhanoa equations were able to estimate the time to the peak MY more accurately than the other equations. In addition, Nelder and Dijkstra equations were able to estimate the peak time at second and third parities more accurately than other equations, respectively. Brody function provided more accurate predictions of peak MY over the first three parities of buffaloes. There was generally a positive relationship between 305-day MY and persistency measures and also between peak yield and 305-day MY, calculated by different models, within each lactation in the current study. Overall, evaluation of the different equations used in the current study indicated the potential of the non-linear models for fitting monthly productive records of buffaloes.     

Modelling of lactation curves of dairy cows based on monthly test day milk yield records under inconsistent milk recording scenarios

The objective of this study was to describe the lactation curve of dairy cattle in Kenya using a suitable lactation function in order to facilitate inclusion of partial lactations in national dairy cattle evaluation and to assess the effect of data characteristics on lactation curve parameters. Six functions were fitted to test day (TD) milk yield records from six parities of Ayrshire, Guernsey, Holstein Friesian, Jersey and Sahiwal cattle. Five datasets: DS-1 (12-TD dataset with randomly missing records), DS-2 (10-TD dataset without missing records), DS-3 (10-TD dataset with randomly missing records), DS-4 (7-TD dataset, with only TD 4 to 10 records) and DS-5 (7-TD dataset, with TD 1 to 4, 6, 8 and 10 records) depicting various recording circumstances were derived to assess the effects of data characteristics on lactation curves and to assess the feasibility of reducing the number of TD samples per lactation. The fit of the functions was evaluated using adjusted R(2) and their predictive abilities were compared using mean square prediction error, percentage of squared bias and the correlation between the predicted and actual milk yield. These criteria plus the changes in the parameters of curve functions and their associated standard errors were used in determining the effects of data characteristics on lactation curves. The mechanistic functions of Dijkstra (DIJ) and Pollott (APOL), and the incomplete gamma function of Wood (WD) had the highest adjusted R(2) > 0.75. The APOL function was eliminated due to convergence failures when analysis of individual lactations within breeds was carried out. Both DIJ and WD had good predictive ability, although DIJ performed slightly better. Convergence difficulties were noted in some DIJ analysis where data were limiting. Missing records, especially at the beginning of a lactation, greatly influenced parameters a and b of the functions. It also resulted in estimates with large standard errors. Missing records in later lactation hardly affected the parameter estimates. The WD and DIJ functions showed superior fit to the data. The WD function demonstrated higher adaptability to various data characteristics than DIJ and could be used in situations where animal recording is not consistently practised and where recording of animal performance is routinely practised. DIJ function had high data requirements, which restricts it to dairy systems with consistent recording, despite easy physiological interpretation of its parameters. The number of TD per lactation could be reduced by minimising sampling frequency in the later lactation while maintaining the monthly sampling frequency in early lactation.     

Analysis of lactation shapes in extended lactations

In order to describe the temporal evolution of milk yield (MY) and composition in extended lactations, 21 658 lactations of Italian Holstein cows were analyzed. Six empirical mathematical models currently used to fit 305 standard lactations (Wood, Wilmink, Legendre, Ali and Schaeffer, quadratic and cubic splines) and one function developed specifically for extended lactations (a modification of the Dijkstra model) were tested to identify a suitable function for describing patterns until 1000 days in milk (DIM). Comparison was performed on individual patterns and on average curves grouped according to parity (primiparous and multiparous) and lactation length (standard ⩽305 days, and extended from 600 to 1000 days). For average patterns, polynomial models showed better fitting performances when compared with the three or four parameters models. However, LEG and spline regression, showed poor prediction ability at the extremes of the lactation trajectory. The Ali and Schaeffer polynomial and Dijkstra function were effective in modelling average curves for MY and protein percentage, whereas a reduced fitting ability was observed for fat percentage and somatic cell score. When individual patterns were fitted, polynomial models outperformed nonlinear functions. No detectable differences were observed between standard and extended patterns in the initial phase of lactation, with similar values of peak production and time at peak. A considerable difference in persistency was observed between 200 and 305 DIM. Such a difference resulted in an estimated difference between standard and extended cycle of about 7 and 9 kg/day for daily yield at 305 DIM and of 463 and 677 kg of cumulated milk production at 305 DIM for the first- and second-parity groups, respectively. For first and later lactation animals, peak yield estimates were nearly 31 and 38 kg, respectively, and occurred at around 65 and 40 days. The asymptotic level of production was around 9 kg for multiparous cows, whereas the estimate was negative for first parity.     

Genome-wide associations for investigating time-dependent genetic effects for milk production traits in dairy cattle

Phenotypic variation in milk production traits has been described over the course of a lactation as well as between different parities. The objective of this study was to investigate whether variation in production is affected by different loci across lactations. A genome-wide association study (GWAS) using a 50-k SNP chip was conducted in 152 divergent German Holstein Friesian cows to test for association with milk production traits over different lactations. The first four lactations were analysed regarding milk yield, fat, protein, lactose, milk urea nitrogen yield and content as well as somatic cell score. Two approaches were used: (i) Wilmink curve parameters were used to assess the genetic effects over the course of a lactation and (ii) test-day yield deviations (YD) were used as a normative approach for a GWAS. The significant effects were largest for markers affecting curve parameters for which there was a statistical power <0.8 of detection even in this small design. While significant markers for YDs were detected in this study, the power to detect effects of a similar magnitude was only 0.11, suggesting that many loci may have been missed with this approach in the present design. Furthermore, all significant effects were specific for a single lactation, leading to the conclusion that the variance explained by a certain locus changes from lactation to lactation. We confirm the common evidence that most production traits vary in the degree of persistency after the peak as a result of genetic influence.     

Genetic parameters of a biological lactation model: early lactation and secretion rate traits of dairy heifers

Early lactation parameters are difficult to estimate from commercial dairy records due to the small number of records available before the peak of production. A biological model of lactation was used with weekly milk records from a single Holstein herd to estimate these early lactation parameters and the secretion rate of milk from the average cell throughout lactation. A genetic analysis of the lactation curve parameters, calculated curve characteristics and secretion rate traits was undertaken. Early lactation traits were found to have little genetic variation and effectively zero heritability. Secretion rate traits for milk, protein, lactose and water were all moderately heritable and highly genetically correlated (>0.87) but fat secretion rate had lower genetic correlations with the other secretion rates. A similar pattern of correlations was seen between total lactation yield traits for fat, protein, lactose and water. The genetic correlations between the lactation curve traits and the secretion rate traits were calculated. Total milk yield, peak yield and maximum secretion potential were all highly correlated with milk, lactose and water secretion rates but less so with fat and protein secretion rates. In particular, fat secretion rate had a moderate to low genetic correlation with these lactation curve traits. Persistency of lactation was highly correlated with fat and protein secretion rates, more persistent lactations being associated with lower rates of secretion of these milk components. Similar levels of heritability were found, where trait genetic parameters were directly equivalent to those derived from the same dataset by random regression methods. However, by using a biological model of lactation to analyse lactation traits new insights into the biology of lactation are possible and ways to select cows on a range of lactation traits may be achieved.     

The Effect of Calf Gender on Milk Production in Seasonal Calving Cows and Its Impact on Genetic Evaluations

Gender of the calf whose birth initiates lactation could influence whole lactation milk yield of the dam due to hormonal influences on mammary gland development, or through calf gender effects on gestation length. Fetal gender could influence late lactation yields because cows become pregnant at peak lactation. The effects of calf gender sequences in parities 1–3 were assessed by separately fitting animal models to datasets from New Zealand comprising 274 000 Holstein Friesian and 85 000 Jersey cows, decreasing to 12 000 and 4 000 cows by parity 3. The lactation initiated by the birth of a female rather than a male calf was associated with a 0.33–1.1% (p≤0.05) higher milk yield. Female calf gender had carryover effects associated with higher milk yield in second lactations for Holstein Friesians (0.24%; p = 0.01) and third lactations for Jerseys (1.1%; p = 0.01). Cows giving birth to bull calves have 2 day longer gestations, which reduces lactation length in seasonal calving herds. Adding a covariate for lactation length to the animal model eroded some of these calf gender effects, such that calving a female led to higher milk yield only for second lactation Holstein Friesians (1.6%; p = 0.002). The interval centering method generates lower estimates of whole lactation yield when Wood’s lactation curves are shifted to the right by 2 days for male calves and this explained the higher yield in female calves when differences in lactation length were considered. Correlations of estimated breeding values between models including or excluding calf gender sequence were 1.00 for bulls or cows. Calf gender primarily influences milk yield through increased gestation length of male calves, and bias associated with the interval centering method used to estimate whole lactation milk yields. Including information on calf gender is unlikely to have an effect on selection response in New Zealand dairy cattle.     

Evaluation of reference lactation length in Chios dairy sheep

Definition and establishment of a fixed reference lactation length could provide useful tools for on-farm comparison of ewes and flock management as well as genetic evaluations for the breeding programme. The objectives of this study were to (i) evaluate different reference lactation lengths for the Chios dairy sheep and (ii) define the most suitable reference length for the breed. A total of 260 042 test-day milk records from 24 474 ewes in 130 flocks collected between 2003 and 2014 were used; 15 different lactation lengths were evaluated ranging from 120 to 260 days, defined at 10-day intervals as reference for the Chios sheep. The evaluation criteria included: (a) heritability and repeatability of milk yield in each reference lactation, (b) genetic correlation of reference lactation milk yield with actual lactation milk yield and yield at first test-day record and (c) correlated response in reference lactation milk yield from selection based on first test-day milk yield. The latter emulates genetic gains achieved for milk yield based on early lactation selection. Heritability and repeatability estimates of reference lactation milk yield and genetic correlation with actual lactation yield favoured long reference lactations (180 to 230 days). On the contrary, correlation with first test-day record milk yield was higher for short lactations (120 to 170 days). Moreover, selection on first test-day record milk yield would lead to a correlated response in reference yield in 220 days equal to 85% of the highest estimate achieved in the maximum reference length of 260 days (190 days when only considering first lactation milk yield). Based on the results of the present study, an overall reference lactation length for the Chios breed of 220 days post-lambing and a first lactation reference length of 190 days post-lambing are recommended.     

Effects of stage of pregnancy on variance components,daily milk yields and 305-day milk yield in Holstein cows,as estimated by using a test-day model

Pregnancy and calving are elements indispensable for dairy production, but the daily milk yield of cows decline as pregnancy progresses, especially during the late stages. Therefore, the effect of stage of pregnancy on daily milk yield must be clarified to accurately estimate the breeding values and lifetime productivity of cows. To improve the genetic evaluation model for daily milk yield and determine the effect of the timing of pregnancy on productivity, we used a test-day model to assess the effects of stage of pregnancy on variance component estimates, daily milk yields and 305-day milk yield during the first three lactations of Holstein cows. Data were 10 646 333 test-day records for the first lactation; 8 222 661 records for the second; and 5 513 039 records for the third. The data were analyzed within each lactation by using three single-trait random regression animal models: one model that did not account for the stage of pregnancy effect and two models that did. The effect of stage of pregnancy on test-day milk yield was included in the model by applying a regression on days pregnant or fitting a separate lactation curve for each days open (days from calving to pregnancy) class (eight levels). Stage of pregnancy did not affect the heritability estimates of daily milk yield, although the additive genetic and permanent environmental variances in late lactation were decreased by accounting for the stage of pregnancy effect. The effects of days pregnant on daily milk yield during late lactation were larger in the second and third lactations than in the first lactation. The rates of reduction of the 305-day milk yield of cows that conceived fewer than 90 days after the second or third calving were significantly (P<0.05) greater than that after the first calving. Therefore, we conclude that differences between the negative effects of early pregnancy in the first, compared with later, lactations should be included when determining the optimal number of days open to maximize lifetime productivity in dairy cows.     

The effect of lactation length on greenhouse gas emissions from the national dairy herd

Many governments have signed up to greenhouse gas emission (GHGE) reduction programmes under their national climate change obligations. Recently, it has been suggested that the use of extended lactations in dairy herds could result in reduced GHGE. Dairy GHGE were modelled on a national basis and the model was used to compare emissions from lactations of three different lengths (305, 370 and 440 days), and a current ‘base’ scenario on the basis of maintaining current milk production levels. In addition to comparing GHGE from the average ‘National Herd’ under these scenarios, results were used to investigate how accounting for lactations of different lengths might alter the estimation of emissions calculated from the National Inventory methodology currently recommended by Intergovernmental Panel on Climate Change. Data for the three lactation length scenarios were derived from nationally recorded dairy performance information and used in the GHGE model. Long lactations required fewer milking cows and replacements to maintain current milk yield levels than short ones, but GHGEs were found to rise from 1214 t of CO₂ equivalent (CE)/farm per year for lactations of 305 days to 1371 t CE/farm per year for 440-day lactations. This apparent anomaly can be explained by the less efficient milk production (kg milk produced per kg cow weight) found in later lactation, a more pronounced effect in longer lactations. The sensitivity of the model to changes in replacement rate, persistency and level of milk yield was investigated. Changes in the replacement rate from 25% to 20% and in persistency by −10% to +20% resulted in very small changes in GHGE. Differences in GHGE due to the level of milk yield were much more dramatic with animals in the top 10% for yield, producing about 25% less GHGE/year than the average animal. National Inventory results were investigated using a more realistic spread of lactation lengths than recommended for such calculations using emissions calculated in the first part of the study. Current UK emission calculations based on the National Inventory were 329 Gg of methane per year from the dairy herd. Using the national distribution of lactation lengths, this was found to be an underestimate by about 10%. This work showed that the current rise in lactation length or a move towards calving every 18 months would increase GHGE by 7% to 14% compared with the current scenario, assuming the same milk yield in all models. Increased milk yield would have a much greater effect on reducing GHGE than changes to lactation length, replacement rate or persistency. National Inventory methodology appears to underestimate GHGE when the distribution of lactation lengths is considered and may need revising to provide more realistic figures.     

Genetic and phenotypic parameters of lactations longer than 305 days (extended lactations)

Test-day milk yield and somatic cell count data over extended lactation (lactation to 540-600 days) were analysed considering part lactations as different traits and fitting random regression (RR) models. Data on Australian Jersey and Holstein Friesian (HF) were used to demonstrate the shape of the lactation curve and data on HF were used for genetic study. Test-day data from about 100 000 cows that calved between 1998 and 2005 were used for this study. In all analyses, a sire model was used.When part lactations were considered as different traits, protein yield early in the lactation (e.g. first 2 months) had a genetic correlation of about 0.8 with protein yield produced after 300 days of lactation. Genetic correlations between lactation stages that are adjacent to each other were high (0.9 or more) within parity. Across parities, genetic correlations were high for both protein and milk yield if they are within the same stage of lactation. Phenotypic correlations were lower than genetic correlations. Heritability of milk-yield traits estimated from the RR model varied from 0.15 at the beginning of the lactation to as high as 0.37 by the 4th month of lactation. All genetic correlations between different days in milk were positive, with the highest correlations between adjacent days in milk and decreasing correlations with increasing time-span. The pattern of genetic correlations between milk yield in the second 300 days (301 to 600 days of lactation) do not markedly differ from the pattern in the first 300 days of lactation. The lowest estimated genetic correlation was 0.15 between milk yield on days 45 and 525 of lactation. The result from this study shows that progeny of bulls with high estimated breeding values for yield traits and those that produce at a relatively high level in the first few months are the most likely candidates for use in herds favouring extended lactations.     

Milk production and composition in reindeer (Rangifer tarandus): effect of lactational stage

Milk yield and composition of major milk constituents were measured in captive, nursing reindeer. Registration of milk production was performed during two successive lactations (2001 and 2002). The milk yield was significantly affected by week of lactation (P<0.001) and by individual (P<0.001). The lactation curve had an asymmetrical peak 3 weeks postpartum and the milk yield at peak lactation was 983 g/day (range 595-1239). The length of lactation varied from 24 to 26 weeks and average total milk production was 99.5 kg. From peak lactation the milk production decreased linearly (P<0.001) until milk production was terminated. Mean values for content of major milk constituents were 15.5% fat, 9.9% protein and 2.5% lactose. The content of fat and protein increased markedly with the lactation stage (P<0.001), while lactose showed a slight decrease (P<0.001). The milk composition was significantly affected by stage of lactation (P<0.001). There was a marginally significant decrease in protein:fat ratio (P=0.06) as protein was substituted by fat with stage of lactation. The caloric value of the milk averaged 8.7 kJ/g and increased significantly with the stage of lactation (P<0.001). The overall increase in milk gross energy content during lactation was 67.6%. The energy output averaged 7996 kJ/day at peak lactation and decreased significantly during the course of lactation (P=0.002).     

Environmental effects on lactation curves included in a test-day model genetic evaluation

A large number of environmental factors affect the daily milk production of a cow. Lactation curves included in the French test-day model are modelled as a function of days in milk with semi-parametric curves (splines). The proper modelling of environmental effects in the test-day analysis was investigated using test-day records collected from the first three lactations of French Montbéliarde cows from 1988 to 2005. Four lactation-curve effects describing calving month, length of dry period, age at calving and gestation defined within parity-class were fitted. The shape of lactation curves did not depend on year of calving, which can be modelled as a constant over the whole lactation. To reduce computational requirements and time, data were pre-adjusted in a first step for fixed effects with no year interaction, and then used for genetic evaluation. Correlations for each lactation between 305-day estimates of genetic and permanent environment effects computed using pre-adjustment factors obtained at a 4-year interval were virtually one. The use of a two-step procedure had a very limited impact on the estimates of genetic and permanent environment effects. The minimum correlations with values estimated with a one-step procedure were 0.9984 and 0.9974, respectively. The knowledge of systematic environmental effects affecting the cow daily yield through lactation curves offers interesting perspectives to predict future daily milk production.     

Genetic modelling of test day records in dairy sheep using orthogonal Legendre polynomials

Test day milk yields of three lactations in Sfakia sheep were analyzed fitting a random regression (RR) model, regressing on orthogonal polynomials of the stage of the lactation period, i.e. days in milk. Univariate (UV) and multivariate (MV) analyses were also performed for four stages of the lactation period, represented by average days in milk, i.e. 15, 45, 70 and 105 days, to compare estimates obtained from RR models with estimates from UV and MV analyses. The total number of test day records were 790, 1314 and 1041 obtained from 214, 342 and 303 ewes in the first, second and third lactation, respectively. Error variances and covariances between regression coefficients were estimated by restricted maximum likelihood. Models were compared using likelihood ratio tests (LRTs). Log likelihoods were not significantly reduced when the rank of the orthogonal Legendre polynomials (LPs) of lactation stage was reduced from 4 to 2 and homogenous variances for lactation stages within lactations were considered. Mean weighted heritability estimates with RR models were 0.19, 0.09 and 0.08 for first, second and third lactation, respectively. The respective estimates obtained from UV analyses were 0.14, 0.12 and 0.08, respectively. Mean permanent environmental variance, as a proportion of the total, was high at all stages and lactations ranging from 0.54 to 0.71. Within lactations, genetic and permanent environmental correlations between lactation stages were in the range from 0.36 to 0.99 and 0.76 to 0.99, respectively. Genetic parameters for additive genetic and permanent environmental effects obtained from RR models were different from those obtained from UV and MV analyses.     

Relationships between fertility, peak milk yields and lactational persistency in dairy cows

Peak milk yield, lactational persistency and conception rates were studied using 5928 lactation records of high milk-producing cows at three California dairies. Log-linear analysis was used to study relationships between peak milk yield, lactational persistency, dairy of origin, lactation number and conception rates in 3850 completed lactations. Cows with peak milk yields greater than the median (38.2 kg milk per day) were less likely to have conceived in one or two breedings than cows with peak milk yields lower than or equal to the median. Cows with a higher than median (0.755) lactational persistency were less likely to have conceived in one or two breedings than cows with a lactational persistency lower than or equal to the median. Dairy of origin had a significant effect on the probability of conceiving in one or two breedings. Cows in the first lactation were more likely than those in subsequent lactations to conceive in one or two breedings. This retrospective study demonstrated that subfertility is associated with high peak lactational yields in high milk-producing California cows.     

Estimates of correlations among yield traits and somatic cell score with different models to adjust for bovine somatotropin effects on Holstein dairy cows

Records of Holstein cows from the Dairy Records Processing Center at Raleigh, NC were edited to obtain three data sets: 65,720 first, 50,694 second, and 65,445 later lactations. Correlations among yield traits and somatic cell score were estimated with three different models: 1) bovine somatotropin (bST) administration ignored, 2) bST administration as a fixed effect and 3) administration of bST as part of the contemporary group (herd-year-month-bST). Heritability estimates ranged from 0.13 to 0.17 for milk, 0.12 to 0.20 for fat, 0.14 to 0.16 for protein yields, and 0.08 to 0.09 for somatic cell score. Estimates were less for later than first lactations. Estimates of genetic correlations among yields ranged from 0.35 to 0.85 with no important differences between estimates with the 3 models. Estimates for lactation 2 agreed with estimates for lactation 1. Estimates of genetic correlations for later lactations were generally greater than for lactations 1 and 2 except between milk and protein yields. Estimates of genetic correlations between yields and somatic cell score were mostly negative or small (-0.45 to 0.11). Estimates of environmental correlations among yield traits were similar with all models (0.77 to 0.97). Estimates of environmental correlations between yields and somatic cell score were negative (-0.22 to -0.14). Estimates of phenotypic correlations among yield traits ranged from 0.70 to 0.95. Estimates of phenotypic correlations between yields and somatic cell score were small and negative. For all three data sets and all traits, no important differences in estimates of genetic parameters were found for the two models that adjusted for bST and the model that did not.     

Data transformation for rank reduction in multi-trait MACE model for international bull comparison

Since many countries use multiple lactation random regression test day models in national evaluations for milk production traits, a random regression multiple across-country evaluation (MACE) model permitting a variable number of correlated traits per country should be used in international dairy evaluations. In order to reduce the number of within country traits for international comparison, three different MACE models were implemented based on German daughter yield deviation data and compared to the random regression MACE. The multiple lactation MACE model analysed daughter yield deviations on a lactation basis reducing the rank from nine random regression coefficients to three lactations. The lactation breeding values were very accurate for old bulls, but not for the youngest bulls with daughters with short lactations. The other two models applied principal component analysis as the dimension reduction technique: one based on eigenvalues of a genetic correlation matrix and the other on eigenvalues of a combined lactation matrix. The first one showed that German data can be transformed from nine traits to five eigenfunctions without losing much accuracy in any of the estimated random regression coefficients. The second one allowed performing rank reductions to three eigenfunctions without having the problem of young bulls with daughters with short lactations.     

Relation of Milk Production Loss to Milk Somatic Cell Count

Milk production loss was studied in relation to increased somatic cell count (SCC). Available data were weekly test-day milk yields and SCC (in 1,000 cells/ml), and mastitis incidences. In total, 18,131 records from 274 cows were used. Production loss was determined for test-day kg milk, kg protein, and kg energy-corrected milk. Least-squares analysis of variance was used to estimate the direct effect of Log10(SCC) on production. The recorded measures of production were first corrected for fixed effects, with adjustment factors estimated from a healthy data-set. The average daily milk yield was 19.7 kg/day in first lactation and 22.0 in later lactations. The geometric mean of SCC was 63.1 in first lactation and 107.2 in later lactations. The incidence of clinical mastitis treated by a veterinarian was 19.8% of the lactations-at-risk. Linear relationships were found between the production parameters and Log10(SCC). Quadratic and cubic effects were evaluated, but were found to contribute little to the overall fit of the models. The individual milk yield loss was 1.29 kg/day for each unit increase in Log10(SCC) for cows in first lactation. Milk yield decreased by 2.04 kg/day per unit Log10(SCC) for older cows. Corresponding values for protein yield were 0.042 and 0.067 kg/day for first and later lactations, respectively.     

Prediction of first test day milk yield using historical records in dairy cows

The transition between two lactations remains one of the most critical periods during the productive life of dairy cows. In this study, we aimed to develop a model that predicts the milk yield of dairy cows from test day milk yield data collected in the previous lactation. In the past, data routinely collected in the context of herd improvement programmes on dairy farms have been used to provide insights in the health status of animals or for genetic evaluations. Typically, only data from the current lactation is used, comparing expected (i.e., unperturbed) with realised milk yields. This approach cannot be used to monitor the transition period due to the lack of unperturbed milk yields at the start of a lactation. For multiparous cows, an opportunity lies in the use of data from the previous lactation to predict the expected production of the next one. We developed a methodology to predict the first test day milk yield after calving using information from the previous lactation. To this end, three random forest models (nextMILK_FULL, nextMILK_PH, and nextMILK_P) were trained with three different feature sets to forecast the milk yield on the first test day of the next lactation. To evaluate the added value of using a machine-learning approach against simple models based on contemporary animals or production in the previous lactation, we compared the nextMILK models with four benchmark models. The nextMILK models had an RMSE ranging from 6.08 to 6.24 kg of milk. In conclusion, the nextMILK models had a better prediction performance compared to the benchmark models. Application-wise, the proposed methodology could be part of a monitoring tool tailored towards the transition period. Future research should focus on validation of the developed methodology within such tool.     

Changes in the yield, and carbohydrate, lipid and protein content of milk during lactation in the rat

The milk yield and composition was studied during the first three lactations of a group of rats. Milk yield increased steadily throughout the three lactations, but was somewhat lower during the first than subsequent lactations. Protein concentration was similar during all three lactations and varied little with stage of lactation. In contrast the lactose concentration, which was reasonably constant for the first 8 days post partum, increased thereafter two-fold by the end of the period studied in all three lactations. However, the N-acetyl-neuraminyl lactose concentration showed somewhat reciprocal changes. Considerable variations in the triacylglycerol concentration was found during the first lactation but few changes were observed during subsequent lactations. The free fatty acid concentration was at all times low and showed no significant changes during or between lactations. At most stages of lactation in raw milk, the major fatty acids are palmitate, oleate and linoleate. However, as lactation progresses there is an increase in the proportion of medium-chain saturated fatty acids and a corresponding decrease in the proportion of long chain unsaturated fatty acids in milk fat. Clearly the composition of milk is not invariable but changes both during and between lactations. Such changes may be expected to have some influence on the metabolism of the offspring.     

Genetic analysis for mastitis resistance and milk somatic cell score in French Lacaune dairy sheep

Genetic analysis for mastitis resistance was studied from two data sets. Firstly, risk factors for different mastitis traits, i.e. culling due to clinical or chronic mastitis and subclinical mastitis predicted from somatic cell count (SCC), were explored using data from 957 first lactation Lacaune ewes of an experimental INRA flock composed of two divergent lines for milk yield. Secondly, genetic parameters for SCC were estimated from 5 272 first lactation Lacaune ewes recorded among 38 flocks, using an animal model. In the experimental flock, the frequency of culling due to clinical mastitis (5%) was lower than that of subclinical mastitis (10%) predicted from SCC. Predicted subclinical mastitis was unfavourably associated with the milk yield level. Such an antagonism was not detected for clinical mastitis, which could result, to some extent, from its low frequency or from the limited amount of data. In practice, however, selection for mastitis resistance could be limited in a first approach to selection against subclinical mastitis using SCC. The heritability estimate of SCC was 0.15 for the lactation mean trait and varied from 0.04 to 0.12 from the first to the fifth test-day. The genetic correlation between lactation SCC and milk yield was slightly positive (0.15) but showed a strong evolution during lactation, i.e. from favourable (-0.48) to antagonistic (0.27). On a lactation basis, our results suggest that selection for mastitis resistance based on SCC is feasible. Patterns for genetic parameters within first lactation, however, require further confirmation and investigation.