Genetic variance in micro-environmental sensitivity for milk and milk quality in Walloon Holstein cattle

Animals that are robust to environmental changes are desirable in the current dairy industry. Genetic differences in micro-environmental sensitivity can be studied through heterogeneity of residual variance between animals. However, residual variance between animals is usually assumed to be homogeneous in traditional genetic evaluations. The aim of this study was to investigate genetic heterogeneity of residual variance by estimating variance components in residual variance for milk yield, somatic cell score, contents in milk (g/dL) of 2 groups of milk fatty acids (i.e., saturated and unsaturated fatty acids), and the content in milk of one individual fatty acid (i.e., oleic acid, C18:1 cis-9), for first-parity Holstein cows in the Walloon Region of Belgium. A total of 146,027 test-day records from 26,887 cows in 747 herds were available. All cows had at least 3 records and a known sire. These sires had at least 10 cows with records and each herd × test-day had at least 5 cows. The 5 traits were analyzed separately based on fixed lactation curve and random regression test-day models for the mean. Estimation of variance components was performed by running iteratively expectation maximization-REML algorithm by the implementation of double hierarchical generalized linear models. Based on fixed lactation curve test-day mean models, heritability for residual variances ranged between 1.01 × 10⁻³ and 4.17 × 10⁻³ for all traits. The genetic standard deviation in residual variance (i.e., approximately the genetic coefficient of variation of residual variance) ranged between 0.12 and 0.17. Therefore, some genetic variance in micro-environmental sensitivity existed in the Walloon Holstein dairy cattle for the 5 studied traits. The standard deviations due to herd × test-day and permanent environment in residual variance ranged between 0.36 and 0.45 for herd × test-day effect and between 0.55 and 0.97 for permanent environmental effect. Therefore, nongenetic effects also contributed substantially to micro-environmental sensitivity. Addition of random regressions to the mean model did not reduce heterogeneity in residual variance and that genetic heterogeneity of residual variance was not simply an effect of an incomplete mean model.     

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (?0.62 to ?0.59) and with protein yield (?0.32 to ?0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.     

Variation in fatty acid contents of milk and milk fat within and across breeds

The aim of this research was to study the potential for selection of cows with a higher nutritional quality of milk fat by studying the differences in fatty acid profiles within and across the following breeds: Dual Purpose Belgian Blue, Holstein-Friesian, Jersey, Montbeliarde, and non-Holstein Meuse-Rhine-Yssel type Red and White. Six hundred milk samples from 275 animals were taken from 7 herds. Several types of fatty acids in milk and milk fat were quantified using mid-infrared spectrometry and previously obtained calibration equations. Statistical analyses were made using a mixed linear model with a random animal effect. The variance components were estimated by using REML. Results showed breed differences for the fatty acid profile. The repeatability estimate obtained in the present study may suggest the existence of moderate additive genetic variance for the fatty acid profile within each breed. Results also indicated variation for each analyzed milk component in the whole cow population studied. Genetic improvement of the nutritional quality of milk fat based on fatty acid profiles might be possible, and further research and development are warranted.     

Effects of fatty acid supplements on milk yield and energy balance of lactating dairy cows

Saturated and unsaturated fatty acid supplements (FS) were evaluated for effects on yield of milk and milk components, concentration of milk components including milk fatty acid profile, and energy balance. Eight ruminally and duodenally cannulated cows and 8 noncannulated cows were used in a replicated 4 × 4 Latin square design experiment with 21-d periods. Treatments were control and a linear substitution of 2.5% fatty acids from saturated FS (SAT; prilled, hydrogenated free fatty acids) for partially unsaturated FS (UNS; calcium soaps of long-chain fatty acids). The SAT treatment did not change milk fat concentration, but UNS linearly decreased milk fat in cannulated cows and tended to decrease milk fat in noncannulated cows compared with control. Milk fat depression with UNS corresponded to increased concentrations of trans-10, cis-12 conjugated linoleic acid and trans C18:1 fatty acids in milk. Milk fat profile was similar for SAT and control, but UNS decreased concentration of short- and medium-chain FA. Digestible energy intake tended to decrease linearly with increasing unsaturated FS in cannulated and noncannulated cows. Increasing unsaturated FS linearly increased empty body weight and net energy gain in cannulated cows, whereas increasing saturated FS linearly increased plasma insulin. Efficiency of conversion of digestible energy to milk tended to decrease linearly with increasing unsaturated FS for cannulated cows only. Addition of SAT provided little benefit to production and energy balance, whereas UNS decreased energy intake and milk energy yield.     

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.     

Terpenes and fatty acid profiles of milk fat and “Bitto” cheese as affected by transhumance of cows on different mountain pastures

The evolution of fatty acid (FA) and terpenoid profiles was studied in milk (n = 20) and “Bitto” (n = 3), a protected designation of origin cheese produced in a restricted Italian alpine area. Milk came from 25 Italian Brown cows successively grazing pastures at 1400, 2100 and 2200 m during transhumance in June–September 2006. The fat matter was analyzed for FAs and terpenes by means of gas chromatography and purge & trap/gas chromatography–mass spectrometry, respectively. FA composition of milk fat varied significantly (p < 0.0001) in relation to contents of conjugated linoleic acid (CLA), stearic, linoleic and trans-vaccenic acids. Similar monoterpene profiles characterized milk fat from cows grazing the different pastures and the highest amount of terpenes was measured in milk coming from cows grazing at 1400 m. High levels of δ3-carene in milk fat were likely related to the important presence of Ligusticum mutellina in the pasture. Only negligible amounts of sesquiterpenes were detected in milk fat whereas they were the most abundant class in fodder. Both FA and terpene profiles of ripened (70 days) cheeses resembled those of the original milks. Overall, results confirm the influence of the botanical composition of mountain pastures both in enhancing the ruminal synthesis of CLA and in modifying the FA and terpenoid profiles of milk and “Bitto” cheese. Nevertheless, neither the FA nor the terpenoid profiles revealed here can be considered as “unique” to “Bitto” cheese and, for this reason, they can hardly be assumed to be biomarkers for defining a specific relationship among grazing area, milk and “Bitto” cheese. They better represent the chemical fingerprint of the cow feeding, adopted in mountain areas.     

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

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

Butter composition and texture from cows with different milk fatty acid compositions fed fish oil or roasted soybeans

Changing the milk fatty acid composition can improve the nutritional and physical properties of dairy products and their acceptability to consumers. A more healthful milk fatty acid composition can be achieved by altering the cow's diet, for example, by feeding supplemental fish oil (FO) or roasted soybeans (RSB), or by selecting cows with a more unsaturated milk fatty acid composition. We examined whether feeding supplemental FO or RSB to cows that had a more unsaturated milk fatty acid composition acted additively to produce butter with improved fatty acid composition and texture. Using a 3 × 3 Latin square design with 2 replications, we fed diets to multiparous Holstein cows (60 to 200 DIM) chosen for producing either more or less unsaturated milk fatty acid composition (n = 6 for each group) for three 3-wk periods. The control diet contained 3.7% crude fat and the 2 experimental diets contained, on a dry matter basis, 0.8% of additional lipids in the form of 0.9% of FO or 5% of RSB. The milk, collected in the third week of feeding, was used to make butter, which was analyzed for its fatty acid composition and physical properties. Dry matter intake, milk yield, and milk composition were not significantly affected by cow diet or by cow selection. Cows that produced a more unsaturated and healthful milk fat prior to the feeding study, according to a “health-promoting index” [HPI = (sum of % of unsaturated fatty acids)/ (%12:0 + 4 × %14:0 + %16:0)], maintained a higher HPI in their butter during the feeding study than did cows with a low HPI. Milk from cows fed supplemental FO or RSB yielded more unsaturated butters with a higher HPI. This butter also was softer when the cows were fed RSB. Feeding RSB to cows chosen for their high milk HPI yielded the most unsaturated butter with the highest HPI and softest texture. Thus, selecting cows with a more health-promoting milk fatty acid composition and feeding supplemental RSB can be used in combination to produce butter that has a consumer-friendly texture and a healthful fatty acid profile.     

Use of herd solutions from a random regression test-day model for diagnostic dairy herd management

In a random regression test-day model, environmental effects in addition to individual animal factors can be included and analyzed. Moreover, instead of herd-year classification of the management groups, the herd-test-day classification within the model better accounts for month-to-month short-term environmental variation in production and somatic cell count (SCC) traits. The herd management levels of milk yield (milk deviation from whole-country mean, kilograms/day), protein and fat concentration (protein and fat deviation, %), and SCC (SCC deviation, 1,000 cells/mL) are used in the dairy herd management Web application “Maitoisa” (in English, “Milky”). This management tool helps to recognize several management problems. For recognition of systematic patterns and single unusual test-days, a monthly time-trend analysis was developed to smooth the random fluctuations and display the yearly production pattern. In addition to analyzing single test-day deviations from the mean, modeled herd solutions assist users in identifying repeated phenomena and enable the forecasting of the management pattern for the subsequent year. The solutions are displayed either as tables or graphs plotted by calendar months. In addition to management effects of the farmer's own herd, he or she can request country or region percentiles to be displayed in the graphs. The Web service has been offered to farmers and dairy advisors since 2001, and it has proved to be a powerful tool for herd monitoring and planning.     

Relationships between milk culture results and milk yield in Norwegian dairy cattle

Associations between test-day milk yield and positive milk cultures for Staphylococcus aureus, Streptococcus spp., and other mastitis pathogens or a negative milk culture for mastitis pathogens were assessed in quarter milk samples from randomly sampled cows selected without regard to current or previous udder health status. Staphylococcus aureus was dichotomized according to sparse (≤1,500 cfu/mL of milk) or rich (>1,500 cfu/mL of milk) growth of the bacteria. Quarter milk samples were obtained on 1 to 4 occasions from 2,740 cows in 354 Norwegian dairy herds, resulting in a total of 3,430 samplings. Measures of test-day milk yield were obtained monthly and related to 3,547 microbiological diagnoses at the cow level. Mixed model linear regression models incorporating an autoregressive covariance structure accounting for repeated test-day milk yields within cow and random effects at the herd and sample level were used to quantify the effect of positive milk cultures on test-day milk yields. Identical models were run separately for first-parity, second-parity, and third-parity or older cows. Fixed effects were days in milk, the natural logarithm of days in milk, sparse and rich growth of Staph. aureus (1/0), Streptococcus spp. (1/0), other mastitis pathogens (1/0), calving season, time of test-day milk yields relative to time of microbiological diagnosis (test day relative to time of diagnosis), and the interaction terms between microbiological diagnosis and test day relative to time of diagnosis. The models were run with the logarithmically transformed composite milk somatic cell count excluded and included. Rich growth of Staph. aureus was associated with decreased production levels in first-parity cows. An interaction between rich growth of Staph. aureus and test day relative to time of diagnosis also predicted a decline in milk production in third-parity or older cows. Interaction between sparse growth of Staph. aureus and test day relative to time of diagnosis predicted declining test-day milk yields in first-parity cows. Sparse growth of Staph. aureus was associated with high milk yields in third-parity or older cows after including the logarithmically transformed composite milk somatic cell count in the model, which illustrates that lower production levels are related to elevated somatic cell counts in high-producing cows. The same association with test-day milk yield was found among Streptococcus spp.-positive pluriparous cows.     

Subclinical hypocalcemia,plasma biochemical parameters,lipid metabolism,postpartum disease,and fertility in postparturient dairy cows

A study was conducted to evaluate the potential association between Ca status at calving and postpartum energy balance, liver lipid infiltration, disease occurrence, milk yield and quality parameters, and fertility in Holstein cows. One hundred cows were assigned to 1 of 2 groups based on whole-blood ionized Ca concentration ([iCa]) on the day of calving [d 0; hypocalcemic [iCa] <1.0 mmol/L (n = 51); normocalcemic [iCa] ≥1.0 mmol/L (n = 49)]. Cows were blocked based on calving date and parity. Blood samples were collected approximately 14 d from expected calving date (d −14), the day of calving (d 0), and on d 3, 7, 14, 21, and 35 postpartum for measurement of plasma nonesterified fatty acid, iCa, total Ca, glucose, and total and direct bilirubin concentrations, and plasma aspartate aminotransferase and gamma glutamyl transferase activities. Liver biopsies were obtained from a subset of cows on d 0, 7, and 35 for quantification of lipid content. Milk samples were collected on d 3, 7, 14, 21, and 35 postpartum for measurement of somatic cell count and percentages of protein, fat, and solids-not-fat. Data for peak test-day milk yield, services per conception, and days open were obtained from Dairy Herd Improvement Association herd records. Disease occurrence was determined based on herd treatment records. Hypocalcemic cows had significantly higher nonesterified fatty acids on d 0. Hypocalcemic cows also had significantly more lipid in hepatocytes on d 7 and 35 postpartum. However, no statistically significant differences were observed between groups for plasma aspartate aminotransferase and gamma glutamyl transferase activities or total and direct bilirubin concentrations. Milk protein percentage was lower in hypocalcemic cows on d 21 and 35. However other milk quality variables (somatic cell count, milk fat percentage, and solids-not-fat) and milk yield variables (peak test-day milk yield and 305-d mature-equivalent 4% fat-corrected milk yield) did not differ between groups. No differences were observed between groups in the occurrence of clinical mastitis, ketosis, displaced abomasum, dystocia, retained placenta, metritis, or fertility measures (percentage cycling at 50–60 d postpartum, services per conception, or days open). These data suggest that early lactation fatty acid metabolism differs between cows with subclinical hypocalcemia and their normocalcemic counterparts.     

Effect of days to conception in the previous lactation on the risk of death and live culling around calving

The objective was to quantify the effect of days to conception (days open) in the previous parity on the risk of death and live culling (excluding the disposal reasons “dairy purposes” and “death”) in the proximity of a subsequent calving in Holstein cows enrolled in the Dairy Herd Improvement program. After edits, 2,075,834 observations of cows calving between 2001 and 2007 in herds located in 36 US states primarily east of the Mississippi river were available. The period at risk included the time between 14 d before expected calving to 60 d after calving. Days open were categorized in 6 periods: 0 to 45 d, 46 to 90 d, 91 to 150 d, 151 to 210 d, 211 to 300 d, and 301 to 600 d after calving in the previous parity. Other variables of interest included parity, length of the dry period before calving, relative last test-day milk yield before dry off before calving, season of calving, and the cow's relative 305-d mature equivalent (305ME) milk yield before calving. Control variables were calving year, herd size, and relative herd 305ME milk yield. All variables were categorized. Generalized mixed models were used and included herd as random effect. Interactions between days open and the variables of interest were included in the models. Least squares means for the risk of death between −14 and 60 d after calving were 2.5, 2.5, 2.9, 3.6, 4.4, and 5.8% for increasing categories of days open. Similarly, for the same categories, least squares means for the risk of live culling were 4.5, 5.0, 5.4, 6.1, 6.9, and 8.1%. The effect of days open on the risk of death was slightly greater for third-parity cows, for long dry periods, for cows with low test-day milk yield before dry off before calving, for spring calvings, and for intermediate relative 305ME milk yield before calving. Similar trends for these interactions were found for the risk of live culling, except for relative 305ME milk yield before calving, where the effect of days open was more pronounced for high-producing cows. In conclusion, increased days open in the previous parity were associated with a greater risk of death and live culling around calving.     

Phenotypic and genetic variability of production traits and milk fatty acid contents across days in milk for Walloon Holstein first-parity cows

The objective of this study was to assess the phenotypic and genetic variability of production traits and milk fatty acid (FA) contents throughout lactation. Genetic parameters for milk, fat, and protein yields, fat and protein contents, and 19 groups and individual FA contents in milk were estimated for first-parity Holstein cows in the Walloon Region of Belgium using single-trait, test-day animal models and random regressions. Data included 130,285 records from 26,166 cows in 531 herds. Heritabilities indicated that de novo synthesized FA were under stronger genetic control than FA originating from the diet and from body fat mobilization. Estimates for saturated short- and medium-chain individual FA ranged from 0.35 for C4:0 to 0.44 for C8:0, whereas those for monounsaturated long-chain individual FA were lower (around 0.18). Moreover, de novo synthesized FA were more heritable in mid to late lactation. Approximate daily genetic correlations among traits were calculated as correlations between daily breeding values for days in milk between 5 and 305. Averaged daily genetic correlations between milk yield and FA contents did not vary strongly among FA (around −0.35) but they varied strongly across days in milk, especially in early lactation. Results indicate that cows selected for high milk yield in early lactation would have lower de novo synthesized FA contents in milk but a slightly higher content of C18:1 cis-9, indicating that such cows might mobilize body fat reserves. Genetic correlations among FA emphasized the combination of FA according to their origin: contents in milk of de novo FA were highly correlated with each other (from 0.64 to 0.99). Results also showed that genetic correlations between C18:1 cis-9 and other FA varied strongly during the first 100 d in milk and reinforced the statement that the release of long-chain FA inhibits FA synthesis in the mammary gland while the cow is in negative energy balance. Finally, results showed that the FA profile in milk changed during the lactation phenotypically and genetically, emphasizing the relationship between the physiological status of cow and milk composition.     

Flaxseed supplementation improves fatty acid profile of cow milk

The objective of the study was to determine the effects of adding flaxseed or fish oil to the diet on the milk fatty acid profile of cows. The experiment was conducted in the summer of 2006 and involved 24 Friesian cows that were divided into 3 groups of 8 animals according to different type of fat supplementation: a traditional diet with no fat supplementation, a diet supplemented with whole flaxseed, and a diet supplemented with fish oil. Results suggested that whole flaxseed supplementation positively affects the milk fatty acid profile during summer. In particular, milk from cows receiving flaxseed supplementation showed a decrease in saturated fatty acid, an increase in monounsaturated fatty acid, and, together with the milk from fish oil-supplemented cows, an increase in polyunsaturated fatty acid content compared with milk from control cows. As expected, both fish oil and flaxseed supplementation increased the content of n-3 polyunsaturated fatty acids in milk fat. The increased dietary intake of C18:3 in flaxseed-supplemented cows resulted in increased levels of milk C18:1 trans-11 and increased conjugated linoleic acid C18:2 cis-9,trans-11 by Δ⁹-desaturase activity. Milk from flaxseed-supplemented cows together with the high conjugated linoleic acid content was characterized by low atherogenic and thrombogenic indices, suggesting that its use has less detrimental effects concerning the atherosclerosis and coronary thrombosis risk associated with the consumption of milk and dairy products. In conclusion, flaxseed supplementation improves composition and nutritional properties of milk from cows milked during times of high ambient temperature.     

Evaluation of cooling strategies for pumping of milk - impact of fatty acid composition on free fatty acid levels

Cooling strategies for pumping of raw milk were evaluated. Milk was pumped for 450 s at 31 degrees C, or pumped after cooling to 4 degrees C and subsequently subjected to various incubation times. Two types of milk were used; i.e. milk from cows fed a diet high in saturated fat supplements resulting in significantly larger milk fat globules than the other type of milk which comes from cows fed a low-fat diet that stimulates high de novo fat synthesis. The content of liquid fat was determined by low-field 1H NMR, which showed that milk from cows given the saturated fat diet also contained less liquid fat at both 4 degrees and 31 degrees C than the other type of milk. This can be ascribed to the differences in the fatty acid composition of the milk as a result of the fatty acid composition of the diets. After pumping of the milk at 31 degrees C, measurement of fat globule size distribution revealed a significant coalescence of milk fat globules in the milk obtained from the saturated fat diet due to pumping. Pumping at 4 degrees C or pumping the other type of milk did not result in coalescence of milk fat globules. Formation of free fatty acids increased significantly in both types of milk by pumping at 31 degrees C. Cooling the milk to 4 degrees C immediately before pumping inhibited an increased content of free fatty acids. However, when the milk was incubated at 4 degrees C for 60 min after cooling and then subjected to pumping, a significant increase in the formation of free fatty acids was observed in both types of milk. It is suggested that this increase in free fatty acids is caused by transition of polymorphic crystal forms or higher level of attached lipoprotein lipases to the milk fat globule before pumping.     

Modeling milk urea of Walloon dairy cows in management perspectives

The aim of this study was to develop an adapted random regression test-day model for milk urea (MU) and to study the possibility of using predictions and solutions given by the model for management purposes. Data included 607,416 MU test-day records of first-lactation cows from 632 dairy herds in the Walloon Region of Belgium. Several advanced features were used. First, to detect the herd influence, the classical herd × test-day effect was split into 3 new effects: a fixed herd × year effect, a fixed herd × month-period effect, and a random herd test-day effect. A fixed time period regression was added in the model to take into account the yearly oscillations of MU on a population scale. Moreover, first autoregressive processes were introduced and allowed us to consider the link between successive test-day records. The variance component estimation indicated that large variance was associated with the random herd × test-day effect (48% of the total variance), suggesting the strong influence of herd management on the MU level. The heritability estimate was 0.13. By comparing observed and predicted MU levels at both the individual and herd levels, target ranges for MU concentrations were defined to take into account features of each cow and each herd. At the cow level, an MU record was considered as deviant if it was <200 or >400 mg/L (target range used in the field) and if the prediction error was >50 mg/L (indicating a significant deviation from the expected level). Approximately 7.5% of the MU records collected between June 2007 and May 2008 were beyond these thresholds. This combination allowed for the detection of potentially suspicious cows. At the herd level, the expected MU level was considered as the sum of the solutions for specific herd effects. A herd was considered as deviant from its target range when the prediction error was greater than the standard deviation of MU averaged by herd test day. Results showed that 6.7% of the herd test-day MU levels between June 2007 and May 2008 were considered deviant. These deviations seemed to occur more often during the grazing period. Although theoretical considerations developed in this study should be validated in the field, this research showed the potential use of a test-day model for analyzing functional traits to advise dairy farmers.     

Genetic parameters of saturated and monounsaturated fatty acid content and the ratio of saturated to unsaturated fatty acids in bovine milk

Fatty acid composition influences the nutritional quality of milk and the technological properties of butter. Using a prediction of fatty acid (FA) contents by mid-infrared (MIR) spectrometry, a large amount of data concerning the FA profile in bovine milk was collected. The large number of records permitted consideration of more complex models than those used in previous studies. The aim of the current study was to estimate the effects of season and stage of lactation as well as genetic parameters of saturated (SAT) and monounsaturated (MONO) fatty acid contents in bovine milk and milk fat, and the ratio of SAT to unsaturated fatty acids (UNSAT) that reflect the hardness of butter (SAT:UNSAT), using 7 multiple-trait, random-regression test-day models. The relationship between these FA traits with common production traits was also studied. The data set contained 100,841 test-day records of 11,626 Holstein primiparous cows. The seasonal effect was studied based on unadjusted means. These results confirmed that milk fat produced during spring and summer had greater UNSAT content compared with winter (63.13 vs. 68.94% of SAT in fat, on average). The effect of stage of lactation on FA profile was studied using the same methodology. Holstein cows in early first lactation produced a lower content of SAT in their milk fat. Variance components were estimated using a Bayesian method via Gibbs sampling. Heritability of SAT in milk (0.42) was greater than heritability of SAT in milk fat (0.24). Estimates of heritability for MONO were also different in milk and fat (0.14 vs. 0.27). Heritability of SAT:UNSAT was moderate (0.27). For all of these traits, the heritability estimates and the genetic and phenotypic correlations varied through the lactation.     

Production performance and milk composition of dairy cows fed whole or ground flaxseed with or without monensin

Eight multiparous Holstein cows averaging 570 ± 43 kg of body weight and 60 ± 20 d in milk were used in a double Latin square design with four 21-d experimental periods to determine the effects of feeding ground or whole flaxseed with or without monensin supplementation (0.02% on a dry matter basis) on milk production and composition, feed intake, digestion, blood composition, and fatty acid profile of milk. Intake of dry matter was similar among treatments. Cows fed whole flaxseed had higher digestibility of acid detergent fiber but lower digestibilities of crude protein and ether extract than those fed ground flaxseed; monensin had no effect on digestibility. Milk production tended to be greater for cows fed ground flaxseed (22.8 kg/d) compared with those fed whole flaxseed (21.4 kg/d). Processing of flaxseed had no effect on 4% fat-corrected milk yield and milk protein and lactose concentrations. Monensin supplementation had no effect on milk production but decreased 4% fat-corrected milk yield as a result of a decrease in milk fat concentration. Feeding ground compared with whole flaxseed decreased concentrations of 16:0, 17:0, and cis6-20:4 and increased those of cis6-18:2, cis9, trans11-18:2, and cis3-18:3 in milk fat. As a result, there was a decrease in concentrations of medium-chain and saturated fatty acids and a trend for higher concentrations of long-chain fatty acids in milk fat when feeding ground compared with whole flaxseed. Monensin supplementation increased concentrations of cis9 and trans11-18:2 and decreased concentrations of saturated fatty acids in milk fat. There was an interaction between flaxseed processing and monensin supplementation, with higher milk fat concentration of trans11-18:1 for cows fed ground flaxseed with monensin than for those fed the other diets. Flaxseed processing and monensin supplementation successfully modified the fatty acid composition of milk fat that might favor nutritional value for consumers.     

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.