Random regression models using different functions to model milk flow in dairy cows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fluoride concentrations of milk, infant formulae, and soy-based products commercially available in Brazil

To assess the fluoride (F) content in commercially available milk formulae in Brazil and to estimate the F intake in children from this source in the first year of life. Samples of cow's milk (n = 51), infant formulae (n = 15), powdered milk (n = 13), and soy-based products (n = 4) purchased in Araçatuba (Brazil) had their F content measured using an ion-specific electrode, after hexamethyldisiloxane-facilitated diffusion. Powdered milk and infant formulae were reconstituted with deionized water, while ready-to-drink products were analyzed without any dilution. Using average infant body masses and suggested volumes of formula consumption for infants 1-12 months of age, possible F ingestion per body mass was estimated. Data were analyzed by descriptive analysis. Mean F content ranged from 0.02 to 2.52 mg/L in all samples. None of the cow's milk provided F intake higher than 0.07 mg/kg. However, two infant formulae, one powdered milk, and one soy-milk led to a daily F intake above the suggested threshold for fluorosis when reconstituted with deionized water. Assuming reconstitution of products with tap water at 0.7 ppm F, two infant formulae, five powdered milk, and four soymilks led to daily F intake ranging from 0.108 to 0.851 mg/kg. The results suggest that the consumption of some brands of infant formulae, powdered milk, and soy-based milk in the first year of age could increase the risk of dental fluorosis, reinforcing the need for periodic surveillance of the F content of foods and beverages typically consumed by young children.

Short communication: Principal components and factor analytic models for test-day milk yield in Brazilian Holstein cattle

A total of 46,089 individual monthly test-day (TD) milk yields (10 test-days), from 7,331 complete first lactations of Holstein cattle were analyzed. A standard multivariate analysis (MV), reduced rank analyses fitting the first 2, 3, and 4 genetic principal components (PC2, PC3, PC4), and analyses that fitted a factor analytic structure considering 2, 3, and 4 factors (FAS2, FAS3, FAS4), were carried out. The models included the random animal genetic effect and fixed effects of the contemporary groups (herd-year-month of test-day), age of cow (linear and quadratic effects), and days in milk (linear effect). The residual covariance matrix was assumed to have full rank. Moreover, 2 random regression models were applied. Variance components were estimated by restricted maximum likelihood method. The heritability estimates ranged from 0.11 to 0.24. The genetic correlation estimates between TD obtained with the PC2 model were higher than those obtained with the MV model, especially on adjacent test-days at the end of lactation close to unity. The results indicate that for the data considered in this study, only 2 principal components are required to summarize the bulk of genetic variation among the 10 traits.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

Validation of models for predicting milk urea nitrogen concentrations, estimating dry matter intake by the NRC (2001)

The objective of this study was to validate three different models for predicting milk urea nitrogen using field conditions, attempting to evaluate the nutritional adequacy diets for dairy cows and prediction of nitrogen excreted to the environment. Observations (4,749) from 855 cows were used. Milk yield, body weight (BW), days in milk and parity were recorded on the milk sampling days. Milk was sampled monthly, for analysis of milk urea nitrogen (MUN), fat, protein, lactose and total solids concentration and somatic cells count. Individual dry matter intake was estimated using the NRC (2001). The three models studied were derived from a first one to predict urinary nitrogen (UN). Model 1 was MUN = UN/12.54, model 2 was MUN = UN/17.6 and model 3 was MUN = UN/(0.0259 × BW), adjusted by body weight effect. To evaluate models, they were tested for accuracy, precision and robustness. Despite being more accurate (mean bias = 0.94 mg/dL), model 2 was less precise (residual error = 4.50 mg/dL) than model 3 (mean bias = 1.41 and residual error = 4.11 mg/dL), while model 1 was the least accurate (mean bias = 6.94 mg/dL) and the least precise (residual error = 5.40 mg/dL). They were not robust, because they were influenced by almost all the variables studied. The three models for predicting milk urea nitrogen were different with respect to accuracy, precision and robustness.

Mastitis diagnostics: quantitative PCR for Staphylococcus aureus genotype B in bulk tank milk

A novel real-time quantitative PCR assay for detecting the pathogenic and contagious Staphylococcus aureus genotype B (GTB) in bulk tank milk was developed and evaluated. The detection of this pathogen in bulk tank milk would greatly facilitate its control, as it is responsible for great economic loss in Swiss dairy herds. The assay is based on the simultaneous detection of 3 GTB-typical target sequences, including 2 enterotoxin genes and a polymorphism within the leucotoxin E gene. A variety of mastitis-associated bacteria was used to validate the assays, resulting in an analytical specificity of 100% and high repeatability. The analytical sensitivity in milk was 40 cfu/mL. An exponential association between simulated cow prevalence and quantitative PCR result was observed. An initial field study revealed 1 GTB-positive herd among the 33 studied herds. This novel assay for bulk tank milk analysis is suitable for routine purposes and is expected to be an effective tool for minimizing Staph. aureus GTB in Swiss dairy herds.

Suitability of low-dosage oxytocin treatment to induce milk ejection in dairy cows

Chronic use of high oxytocin (OT) dosages can cause a reduced response to endogenous OT. In this study the OT dosages used in the milking practice of 82 dairy cow farms were recorded. The OT dosages per cow used were high, especially when injected i.m. (23+/-2 IU) compared with i.v. (7+/-1 IU). In addition, the minimum OT dosages needed to obtain normal milk removal in cows with disturbed milk ejection were investigated. Seventeen cows routinely treated with OT during milking (group T) and 17 cows without previous OT treatment were used (group C). After cessation of spontaneous milk flow, both T and C groups were injected i.v. with a low dosage of OT (0.2 or 0.5 IU/cow). The time from injection until cessation of the OT-induced milk flow was recorded (response phase). The response phase and the amounts of removed milk by effect of the OT injection increased with increasing OT dosage. Values for 0.2 and 0.5 IU/cow of OT injected i.v. were (response phase and amount of milk removed) 198+/-27 and 302+/-18s and 3.4+/-0.7 kg and 6.5+/-1.3 kg, respectively, for the C group, and 157+/-15 and 221+/-16s and 3.2+/-0.5 and 5.5+/-1.0 kg, respectively, for the T group. Within 20 min of the OT injection, plasma concentrations returned to basal levels. The threshold OT concentration at cessation of milk flow after injection of 0.2 or 0.5 IU/cow of OT was calculated based on the OT plasma half-life. The threshold increased with increasing dosages of OT and was higher in group T (8+/-1 and 14+/-1 pg/mL for 0.2 and 0.5 IU/cow, respectively) than in group C (7+/-1 and 11+/-1 pg/mL for 0.2 and 0.5 IU/cow, respectively). In conclusion, desensitization of the udder toward OT occurs when the udder is exposed to elevated OT plasma concentrations, both short-term during the actual milking and long-term due to chronic high-dosage OT treatment. However, low-dosage OT treatments to induce normal milk removal can minimize the observed side effects.

Proteomic and peptidomic study of proteolysis in quarter milk after infusion with lipoteichoic acid from Staphylococcus aureus

Mastitic milk is associated with increased bovine protease activity, such as that from plasmin and somatic cell enzymes, which cause proteolysis of the caseins and may reduce cheese yield and quality. The aim of this work was to characterize the peptide profile resulting from proteolysis in a model mastitis system and to identify the proteases responsible. One quarter of each of 2 cows (A and B) was infused with lipoteichoic acid from Staphylococcus aureus. The somatic cell counts of the infused quarters reached a peak 6h after infusion, whereas plasmin activity of those quarters also increased, reaching a peak after 48 and 12h for cow A and B, respectively. Urea-polyacrylamide gel electrophoretograms of milk samples of cow A and B obtained at different time points after infusion and incubated for up to 7 d showed almost full hydrolysis of beta- and alpha(S1)-casein during incubation of milk samples at peak somatic cell counts, with that of beta-casein being faster than that of alpha(S1)-casein. Two-dimensional gel electrophoretograms of milk 6h after infusion with the toxin confirmed hydrolysis of beta- and alpha(S1)-casein and the appearance of lower-molecular-weight products. Peptides were subsequently separated by reversed-phase HPLC and handmade nanoscale C(18) columns, and identified by matrix-assisted laser desorption/ionization time-of-flight tandem mass spectrometry. Twenty different peptides were identified and shown to originate from alpha(s1)- and beta-casein. Plasmin, cathepsin B and D, elastase, and amino- and carboxypeptidases were suggested as possible responsible proteases based on the peptide cleavage sites. The presumptive activity of amino- and carboxypeptidases is surprising and may indicate the activity of cathepsin H, which has not been reported in milk previously.

Escherichia coli producing CMY-2 β-lactamase in bovine mastitis milk

An Escherichia coli isolate producing the CMY-2 β-lactamase was found in the milk of a cow with recurrent subclinical mastitis. The isolate was resistant to the antibiotics commonly used for intramammary mastitis treatment, such as penicillins, cephalosporins, β-lactam/β-lactamase inhibitor combinations, aminoglycosides, tetracyclines, and sulfonamides. This is the first report of a plasmid-mediated AmpC-producing Enterobacteriaceae in bovine milk.

Milk cortisol concentration in automatic milking systems compared with auto-tandem milking parlors

Milk cortisol concentration was determined under routine management conditions on 4 farms with an auto-tandem milking parlor and 8 farms with 1 of 2 automatic milking systems (AMS). One of the AMS was a partially forced (AMSp) system, and the other was a free cow traffic (AMSf) system. Milk samples were collected for all the cows on a given farm (20 to 54 cows) for at least 1 d. Behavioral observations were made during the milking process for a subset of 16 to 20 cows per farm. Milk cortisol concentration was evaluated by milking system, time of day, behavior during milking, daily milk yield, and somatic cell count using linear mixed-effects models. Milk cortisol did not differ between systems (AMSp: 1.15 +/- 0.07; AMSf: 1.02 +/- 0.12; auto-tandem parlor: 1.01 +/- 0.16 nmol/L). Cortisol concentrations were lower in evening than in morning milkings (1.01 +/- 0.12 vs. 1.24 +/- 0.13 nmol/L). The daily periodicity of cortisol concentration was characterized by an early morning peak and a late afternoon elevation in AMSp. A bimodal pattern was not evident in AMSf. Finally, milk cortisol decreased by a factor of 0.915 in milking parlors, by 0.998 in AMSp, and increased by a factor of 1.161 in AMSf for each unit of ln(somatic cell count/1,000). We conclude that milking cows in milking parlors or AMS does not result in relevant stress differences as measured by milk cortisol concentrations. The biological relevance of the difference regarding the daily periodicity of milk cortisol concentrations observed between the AMSp and AMSf needs further investigation.

Cell population, viability, and some key immunomodulatory molecules in different milk somatic cell samples in dairy cows

Immune cells in the milk are most important in combating pathogens that invade the mammary gland. This study investigated the immune competence and viability of somatic milk cells that are already resident in milk and udders free of infection. Cells were studied in freshly removed milk to simulate conditions in the udder. Effects of incubation, cell preparation, and immunological stimulation with 0.5 mug/ml lipopolysaccharide (LPS) from Escherichia coli were analysed. Viability and differential counts of milk cells between high and low somatic cell count (SCC) quarters, and cisternal and alveolar milk with and without LPS stimulation were compared. Incubation and preparation of cells caused a cell loss which further increased with time independently of SCC and milk fraction. The viability of these cells was stable until 3 h post incubation and decreased until 6 h. Cell populations differed between both investigations, but did not change during the course of the experiment. mRNA expression of immune and apoptosis factors of the cells, measured by qPCR, did not change substantially: mRNA expression of caspase 3, Toll like receptor 4, and GM-CSF did not change, whereas the expression of the death receptor Fas/APO-1 (CD95), lactoferrin and lysozyme was decreased at 6 h. Cyclooxygenase-2 and TNF-alpha mRNA expression were decreased after 6 h of LPS treatment. In comparison with other studies in vivo or in vitro (in cell culture), in this study where cells are studied ex vivo (removed from the udder but kept in their natural environment, the milk) resident milk cells seem to be more vulnerable, less viable, less able to respond to stimulation, and thus less immune competent compared with cells that have freshly migrated from blood into milk after pathogen stimulation. The cell viability and differential cell count differed between high- and low-SCC milk and between cisternal and alveolar milk depending on the individual cow. In conclusion, the results support the view that for a most effective defence against invading pathogens the mammary gland is reliant on the recruitment of fresh immune cells from the blood.

The energy expenditure of 2 Holstein cow strains in an organic grazing system

Until recently, measurements of energy expenditure (EE; herein defined as heat production) in respiration chambers did not account for the extra energy requirements of grazing dairy cows on pasture. As energy is first limiting in most pasture-based milk production systems, its efficient use is important. Therefore, the aim of the present study was to compare EE, which can be affected by differences in body weight (BW), body composition, grazing behavior, physical activity, and milk production level, in 2 Holstein cow strains. Twelve Swiss Holstein-Friesian (HCH; 616 kg of BW) and 12 New Zealand Holstein-Friesian (HNZ; 570 kg of BW) cows in the third stage of lactation were paired according to their stage of lactation and kept in a rotational, full-time grazing system without concentrate supplementation. After adaption, the daily milk yield, grass intake using the alkane double-indicator technique, nutrient digestibility, physical activity, and grazing behavior recorded by an automatic jaw movement recorder were investigated over 7d. Using the (13)C bicarbonate dilution technique in combination with an automatic blood sampling system, EE based on measured carbon dioxide production was determined in 1 cow pair per day between 0800 to 1400 h. The HCH were heavier and had a lower body condition score compared with HNZ, but the difference in BW was smaller compared with former studies. Milk production, grass intake, and nutrient digestibility did not differ between the 2 cow strains, but HCH grazed for a longer time during the 6-h measurement period and performed more grazing mastication compared with the HNZ. No difference was found between the 2 cow strains with regard to EE (291 ± 15.6 kJ) per kilogram of metabolic BW, mainly due to a high between-animal variation in EE. As efficiency and energy use are important in sustainable, pasture-based, organic milk production systems, the determining factors for EE, such as methodology, genetics, physical activity, grazing behavior, and pasture quality, should be investigated and quantified in more detail in future studies.

Milk production during the colostral period is not related to the later lactational performance in dairy cows

In dairy cows, milk yield increases rapidly after parturition until a peak at around wk 6 of lactation. However, the description of the shape of the lactation curve is commonly based on weekly average milk yields. For a more detailed analysis of the milk production curve from the very beginning of lactation including the colostral period and the effect of colostrum yield on further lactational performance, the first 10 milkings after parturition, daily milk yields from d 1 to 28 of lactation, and the cumulative milk production on d 100 to 305 of lactation were investigated in 17 primiparous and 39 multiparous cows milked twice daily. Milk yield at the first milking after parturition (colostrum) ranged from 1.3 to 20.7kg (Δ=19.4kg) in multiparous and from 1.8 to 10.9kg in primiparous animals (Δ=9.1kg). At the tenth milking, milk production ranged from 9.2 to 21.5kg (Δ=12.3kg) in multiparous and from 7.0 to 15.2kg (Δ=8.2kg) in primiparous animals. Immediately after parturition, daily milk production increased rapidly, but after approximately 1wk in lactation, the slope of the daily milk production curve flattened and continued more linear. A nonlinear regression equation was used to determine this timely change, which occurred earlier in primiparous (d 6.9±0.3) than in multiparous cows (d 8.2±0.2). The correlation between the amount of first colostrum and milk production during further lactation decreased already from 0.47 on d 5 to 0.32 on d 14. In multiparous cows, the correlation between total milk production of the previous 305d standard lactation and the amount of first colostrum was not significant (correlation=0.29), whereas the correlation with the daily production increased from 0.45 on d 5 to 0.69 on d 14. However, in primiparous animals, correlations between first-colostrum yield and daily milk yields up to d 28 of lactation were not significant, possibly due to the smaller sample size compared with multiparous animals. First-colostrum yield and cumulative milk production of 100, 200, and 305 lactation days were not significantly correlated in multiparous and primiparous cows. In conclusion, the milk production during the first few milkings is widely independent from the overall production level of a cow. Potentially, genetic selection toward lower milk yield during the very first days after parturition at a simultaneously high lactational performance may be a tool to ensure sufficient colostrum quality and to reduce the metabolic load around parturition.

A case-control study to estimate the effects of acute clinical infection with the Schmallenberg virus on milk yield, fertility and veterinary costs in Swiss dairy herds.

Schmallenberg virus (SBV) was first detected in Switzerland in July 2012 and many Swiss dairy farmers reported acute clinical signs in dairy cattle during the spread of the virus until December 2012. The objectives of the present study were to investigate the effects of an acute infection with SBV on milk yield, fertility and veterinary costs in dairy farms with clinical signs of SBV infection (case farms), and to compare those farms to a matched control group of dairy farms in which cattle did not show clinical signs of SBV infection. Herd size was significantly (p<0.001) larger in case farms (33 cows, n=77) than in control farms (25 cows, n=84). Within case herds, 14.8% (median) of the cows showed acute clinical signs. Managers from case farms indicated to have observed a higher abortion rate during the year with SBV (6.5%) than in the previous year (3.7%). Analysis of fertility parameters based on veterinary bills and data from the breeding associations showed no significant differences between case and control farms. The general veterinary costs per cow from July to December 2012 were significantly higher (p=0.02) in case (CHF 19.80; EUR 16.50) than in control farms (CHF 15.90; EUR 13.25). No differences in milk yield were found between groups, but there was a significant decrease in milk production in case farms in the second half year in 2012 compared to the same period in 2011 (p<0.001) and 2013 (p=0.009). The average daily milk yield per cow (both groups together) was +0.73kg higher (p=0.03) in the second half year 2011 and +0.52kg (p=0.12) in the second half year 2013 compared to the same half year 2012. Fifty-seven percent of the cows with acute clinical signs (n=461) were treated by a veterinarian. The average calculated loss after SBV infection for a standardized farm was CHF 1606 (EUR 1338), which can be considered as low at the national level, but the losses were subject to great fluctuations between farms, so that individual farms could have very high losses (>CHF 10,000, EUR 8333).

Relationships of Efficiency to Reproductive Disorders in Danish Milk Production: A Stochastic Frontier Analysis

Relationships of various reproductive disorders and milk production performance of Danish dairy farms were investigated. A stochastic frontier production function was estimated using data collected in 1998 from 514 Danish dairy farms. Measures of farm-level milk production efficiency relative to this production frontier were obtained, and relationships between milk production efficiency and the incidence risk of reproductive disorders were examined. There were moderate positive relationships between milk production efficiency and retained placenta, induction of estrus, uterine infections, ovarian cysts, and induction of birth. Inclusion of reproductive management variables showed that these moderate relationships disappeared, but directions of coefficients for almost all those variables remained the same. Dystocia showed a weak negative correlation with milk production efficiency. Farms that were mainly managed by young farmers had the highest average efficiency scores. The estimated milk losses due to inefficiency averaged 1142, 488, and 256 kg of energy-corrected milk per cow, respectively, for low-, medium-, and high-efficiency herds. It is concluded that the availability of younger cows, which enabled farmers to replace cows with reproductive disorders, contributed to high cow productivity in efficient farms. Thus, a high replacement rate more than compensates for the possible negative effect of reproductive disorders. The use of frontier production and efficiency/ inefficiency functions to analyze herd data may enable dairy advisors to identify inefficient herds and to simulate the effect of alternative management procedures on the individual herd's efficiency.