Production,partial cash flows and greenhouse gas emissions of simulated dairy herds with extended lactations

The transition period is the most critical period in the lactation cycle of dairy cows. Extended lactations reduce the frequency of transition periods, the number of calves and the related labour for farmers. This study aimed to assess the impact of 2 and 4 months extended lactations on milk yield and net partial cash flow (NPCF) at herd level, and on greenhouse gas (GHG) emissions per unit of fat- and protein-corrected milk (FPCM), using a stochastic simulation model. The model simulated individual lactations for 100 herds of 100 cows with a baseline lactation length (BL), and for 100 herds with lactations extended by 2 or 4 months for all cows (All+2 and All+4), or for heifers only (H+2 and H+4). Baseline lactation length herds produced 887 t (SD: 13) milk/year. The NPCF, based on revenues for milk, surplus calves and culled cows, and costs for feed, artificial insemination, calving management and rearing of youngstock, was k€174 (SD: 4)/BL herd per year. Extended lactations reduced milk yield of the herd by 4.1% for All+2, 6.9% for All+4, 1.1% for H+2 and 2.2% for H+4, and reduced the NPCF per herd per year by k€7 for All+2, k€12 for All+4, k€2 for H+2 and k€4 for H+4 compared with BL herds. Extended lactations increased GHG emissions in CO₂-equivalents per t FPCM by 1.0% for All+2, by 1.7% for All+4, by 0.2% for H+2 and by 0.4% for H+4, but this could be compensated by an increase in lifespan of dairy cows. Subsequently, production level and lactation persistency were increased to assess the importance of these aspects for the impact of extended lactations. The increase in production level and lactation persistency increased milk production of BL herds by 30%. Moreover, reductions in milk yield for All+2 and All+4 compared with BL herds were only 0.7% and 1.1% per year, and milk yield in H+2 and H+4 herds was similar to BL herds. The resulting NPCF was equal to BL for All+2 and All+4 and increased by k€1 for H+2 and H+4 due to lower costs for insemination and calving management. Moreover, GHG emissions per t FPCM were equal to BL herds or reduced (0% to −0.3%) when lactations were extended. We concluded that, depending on lactation persistency, extending lactations of dairy cows can have a positive or negative impact on the NPCF and GHG emissions of milk production.     

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.     

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.     

Milk yield persistency in Brazilian Gyr cattle based on a random regression model

With the objective of evaluating measures of milk yield persistency, 27,000 test-day milk yield records from 3362 first lactations of Brazilian Gyr cows that calved between 1990 and 2007 were analyzed with a random regression model. Random, additive genetic and permanent environmental effects were modeled using Legendre polynomials of order 4 and 5, respectively. Residual variance was modeled using five classes. The average lactation curve was modeled using a fourth-order Legendre polynomial. Heritability estimates for measures of persistency ranged from 0.10 to 0.25. Genetic correlations between measures of persistency and 305-day milk yield (Y305) ranged from -0.52 to 0.03. At high selection intensities for persistency measures and Y305, few animals were selected in common. As the selection intensity for the two traits decreased, a higher percentage of animals were selected in common. The average predicted breeding values for Y305 according to year of birth of the cows had a substantial annual genetic gain. In contrast, no improvement in the average persistency breeding value was observed. We conclude that selection for total milk yield during lactation does not identify bulls or cows that are genetically superior in terms of milk yield persistency. A measure of persistency represented by the sum of deviations of estimated breeding value for days 31 to 280 in relation to estimated breeding value for day 30 should be preferred in genetic evaluations of this trait in the Gyr breed, since this measure showed a medium heritability and a genetic correlation with 305-day milk yield close to zero. In addition, this measure is more adequate at the time of peak lactation, which occurs between days 25 and 30 after calving in this breed.     

Does increasing milk yield per cow reduce greenhouse gas emissions? A system approach

Milk yield per cow has continuously increased in many countries over the last few decades. In addition to potential economic advantages, this is often considered an important strategy to decrease greenhouse gas (GHG) emissions per kg of milk produced. However, it should be considered that milk and beef production systems are closely interlinked, as fattening of surplus calves from dairy farming and culled dairy cows play an important role in beef production in many countries. The main objective of this study was to quantify the effect of increasing milk yield per cow on GHG emissions and on other side effects. Two scenarios were modelled: constant milk production at the farm level and decreasing beef production (as co-product; Scenario 1); and both milk and beef production kept constant by compensating the decline in beef production with beef from suckler cow production (Scenario 2). Model calculations considered two types of production unit (PU): dairy cow PU and suckler cow PU. A dairy cow PU comprises not only milk output from the dairy cow, but also beef output from culled cows and the fattening system for surplus calves. The modelled dairy cow PU differed in milk yield per cow per year (6000, 8000 and 10 000 kg) and breed. Scenario 1 resulted in lower GHG emissions with increasing milk yield per cow. However, when milk and beef outputs were kept constant (Scenario 2), GHG emissions remained approximately constant with increasing milk yield from 6000 to 8000 kg/cow per year, whereas further increases in milk yield (10 000 kg milk/cow per year) resulted in slightly higher (8%) total GHG emissions. Within Scenario 2, two different allocation methods to handle co-products (surplus calves and beef from culled cows) from dairy cow production were evaluated. Results showed that using the 'economic allocation method', GHG emissions per kg milk decreased with increasing milk yield per cow per year, from 1.06 kg CO2 equivalents (CO2eq) to 0.89 kg CO2eq for the 6000 and 10 000 kg yielding dairy cow, respectively. However, emissions per kg of beef increased from 10.75 kg CO2eq to 16.24 kg CO2eq due to the inclusion of suckler cows. This study shows that the environmental impact (GHG emissions) of increasing milk yield per cow in dairy farming differs, depending upon the considered system boundaries, handling and value of co-products and the assumed ratio of milk to beef demand to be satisfied.     

Accurate mathematical models to describe the lactation curve of Lacaune dairy sheep under intensive management

Although the intensive production system of Lacaune dairy sheep is the only profitable method for producers outside of the French Roquefort area, little is known about this type of systems. This study evaluated yield records of 3677 Lacaune sheep under intensive management between 2005 and 2010 in order to describe the lactation curve of this breed and to investigate the suitability of different mathematical functions for modeling this curve. A total of 7873 complete lactations during a 40-week lactation period corresponding to 201 281 pieces of weekly yield data were used. First, five mathematical functions were evaluated on the basis of the residual mean square, determination coefficient, Durbin Watson and Runs Test values. The two better models were found to be Pollott Additive and fractional polynomial (FP). In the second part of the study, the milk yield, peak of milk yield, day of peak and persistency of the lactations were calculated with Pollot Additive and FP models and compared with the real data. The results indicate that both models gave an extremely accurate fit to Lacaune lactation curves in order to predict milk yields (P = 0.871), with the FP model being the best choice to provide a good fit to an extensive amount of real data and applicable on farm without specific statistical software. On the other hand, the interpretation of the parameters of the Pollott Additive function helps to understand the biology of the udder of the Lacaune sheep. The characteristics of the Lacaune lactation curve and milk yield are affected by lactation number and length. The lactation curves obtained in the present study allow the early identification of ewes with low milk yield potential, which will help to optimize farm profitability.     

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.     

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.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

Effects of continuous lactation and short dry periods on mammary function and animal health

The dry period is required to facilitate cell turnover in the bovine mammary gland in order to optimize milk yield in the next lactation. Traditionally, an 8-week dry period has been a standard management practice for dairy cows based on retrospective analyses of milk yields following various dry period lengths. However, as milk production per cow has increased, transitioning cows from the nonlactating state to peak milk yield has grown more problematic. This has prompted new studies on dry period requirements for dairy cows. These studies indicate a clear parity effect on dry period requirement. First parity animals require a 60-day dry period, whereas lactations following later parities demonstrate no negative impact with 30-day dry period or even eliminating the dry period when somatotropin (ST) is also used to maintain milk yields. Shortened dry periods in first parity animals were associated with reduced mammary cell turnover during the dry period and early lactation and increased numbers of senescent cells and reduced functionality of lactating alveolar mammary cells postpartum. Use of ST and increased milking frequency postpartum reduced the impact of shortened dry periods. The majority of new intramammary infections occur during the dry period and persist into the following lactation. There is therefore the possibility of altering mastitis incidence by modifying or eliminating the dry period in older parity animals. As the composition of mammary secretions including immunoglobulins may be reduced when the dry period is reduced or eliminated, there is the possibility that the immune status of cows during the peripartum period is influenced by the length of the dry period.     

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.     

Factors influencing heifer survival and fertility on commercial dairy farms

The average dairy cow survives only three lactations, reducing the availability of replacement heifers. Prenatal losses occur due to early embryonic mortality (about 40%), later embryo loss (up to 20% in high-yielding herds) or abortion (about 5%). A recent survey of 19 UK herds showed that 7.9% of calves were born dead and 3.4% died within 1 month. During the rearing phase, 6.7% of animals were lost before reaching first service at 15 months due to disease or accident and another 2.3% failed to conceive. Many potential replacements therefore never enter the milking herd. This severely limits opportunities for on-farm selection of breeding cows in addition to presenting a welfare issue and causing economic loss. The most profitable animals once lactation is reached combine good milk production with a regular calving pattern. Some aspects of performance are related to age at first calving (AFC), which in turn is influenced by heifer growth rates. Poorly growing animals required more services to conceive, calved later and subsequently performed badly. Optimum fertility and maximum yield in the first lactation were associated with an AFC of 24 to 25 months. However, heifers calving at 22 to 23 months performed best in terms of total milk yield and survival over the first 5 years, partly because good heifer fertility was associated with better fertility later. We have investigated some possible juvenile predictors of future performance. Low-birth-weight calves were more likely to come from either primiparous mothers or older dams (3+ lactations) with higher peak milk yields, suggesting that the uterine environment may limit prenatal calf growth due to competition for nutrients with maternal growth or milk production. Linear trait classification scores for frame size show genetic correlations with longevity. The skeletal measures of height and crown rump length in 1-month-old calves was correlated to subsequent stature, and frame size was correlated to weight at 15 months. It may thus be possible to predict performance from simple size measurements as juveniles. Neither endogenous nor stimulated growth hormone (GH) release in 6-month-old calves were related to milk yield in the first three lactations, but size of a stimulated GH peak was positively related to milk energy values in the first lactation. Cows with delayed ovulation (>45 days) in the first lactation had a higher GH pulse amplitude and lower IGF-I as a juvenile. Cows that partition excess energy into milk in their first lactation may suffer reduced longevity.     

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.     

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.     

Effect of different feeding strategies on lactation performance of Holstein and Normande dairy cows

The dairy farming systems of Western Europe are based on a simple feeding system composed of grazed and preserved grass, maize silage and concentrates in variable proportions. There is, nevertheless, a great diversity of feeding strategies between dairy farms. Over 5 years, we studied the direct and delayed effects of four feeding strategies on the lactation and reproduction performances of Holstein and Normande dairy cows. The four feeding strategies (denoted Hh, Hl, Lh and Ll) correspond to two total mixed rations applied in winter from calving to turnout (maize silage with 30% concentrate or grass silage with 15% concentrate), which were subsequently crossed with two levels of concentrate supplementation at grazing to 210 days. Each year, 72 dairy cows managed in grouped winter calving were assigned to the four strategies. Finally, the results of 325 lactations and 295 inseminated cows were analysed. The four strategies resulted in considerable variation in nutrient intake and, in particular, in differences in concentrates consumed, with values of 1407, 1026, 773 and 392 kg dry matter per cow for strategies Hh, Lh, Hl and Ll, respectively. Total milk production (7567, 7015, 6720 and 6238 kg per cow for treatments Hh, Lh, Hl and Ll, respectively), milk fat content (39.0, 37.1, 40.3 and 38.5 g/kg, respectively), milk protein content (33.0, 31.8, 33.1 and 31.6 g/kg, respectively), and the character of the lactation and body condition curves were all highly sensitive to the strategies applied. While no significant interaction was detected on total lactation yield, the Holstein cows reacted more dramatically to each dietary change at each period, compared with the Normande cows. Winter feeding did not affect the production of milk at pasture whereas, at pasture, the milk from the cows of the H groups in winter was higher in milk fat and protein content. Reproduction performance was unaffected by feeding strategy. The Holstein cows, well fed and producing the most milk (Hh and Hb), had the lowest rate of success at first artificial inseminations (21.5%). The dual-purpose Normande cows had a pregnancy rate 10 points higher than Holstein cows. This comparison of strongly contrasting feeding strategies confirms the immediate reactivity of dairy cows (in terms of milk performance and body condition) to variations of nutritive intake throughout lactation, with a weak carryover effect from feeding levels early in lactation. In contrast, reproduction performance was less sensitive to variation in nutrient supply.     

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.     

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.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.