Milk coagulation properties (MCP) are traditionally expressed using rennet coagulation time (RCT), time to curd firmness (CF) of 20mm (k20), and CF 30min after enzyme addition (a30) values, all of ...which are single-point measures taken from the output of computerized renneting meters, such as the Formagraph. Thus, traditional MCP use only some of the available information. Moreover, because of the worldwide spreading of breeds such as the Holstein-Friesian, characterized by late-coagulating milk, it happens often that some samples do not coagulate at all, that a30 is strongly and negatively related to RCT, and that k20 is not measurable. The aim of the present work was to model CF as a function of time (CFt, mm) over a 30-min interval. The model tested was CFt=CFP×1−e−kCF×(t−RCT), where CFP (mm) is the potential asymptotical CF at an infinite time, kCF (min−1) is the curd firming rate constant, and RCT is measured inminutes. The CFt model was initially applied to data of milk of each of 105 Brown Swiss cows from 7 herds, each sampled once (trial 1). Four samples did not coagulate within 30min. Eighty-seven of the 101 individual equations obtained fit the CF data of milk samples very well, even though the samples differed in composition, and were produced by cows of different ages and days in milk, reared on different farms (coefficient of determination >0.99; average residual standard deviation=0.21mm). Samples with a very late RCT (slowly coagulating samples) yielded so few observational data points that curve parameters could not be precisely estimated. The repeatability of CFt equation parameters was estimated using data obtained from 5 replicates of each of 2 samples of bulk milk from 5 Holstein-Friesian cows analyzed every day for 5 consecutive days (trial 2). Repeatability of RCT was better than that of the other 2 parameters. Moreover, traditional MCP values (RCT, a30, and k20) can be obtained from the individual CFt equations, using all available information. The MCP estimated from equations were very similar to the single-point measures yielded by the computerized renneting meter (coefficient of determination >0.97), but repeatability was slightly better. The model allowed the estimation of k20 for samples with a very late coagulation or with very slow curd firming. Finally, the 3 novel parameters used to assess different milk samples were less interdependent than are the traditional measures, and their practical and scientific utility requires further study.
The role of each step of cheese and ricotta making in development of flavor of cheese and other dairy products is not yet well known. The objectives of this study were to characterize volatile ...organic compounds (VOC) in cheese and ricotta making with bulk milk from cows grazing in a highland area and to evaluate their evolution in the various dairy products and by-products obtained during the production processes. A group of 148 cows was grazed day and night on pasture from June to September. A total of 7 cheese-making sessions were carried out using the bulk milk collected every 2 wk during summer pasturing according to the artisanal procedure used for Malga cheese production. All milks, products, and by-products were sampled, and the VOC content of milk, cream, whey, ricotta, scotta (residual liquid), fresh cheeses, and cheeses ripened for 6 and 12 mo was determined by solid-phase microextraction gas chromatography–mass spectrometry. Forty-nine compounds were identified belonging to the following chemical families: alcohols (13), aldehydes (9), esters (8), free fatty acids (6), ketones (5), lactones (2), sulfurs (2), terpenes (2), phenol (1), and benzene (1). The results showed that the amounts of VOC in the various dairy products differed significantly. Comparisons between the VOC of 4 types of milk (whole evening, skim evening, whole morning, mixed in the vat) showed that the skimming process had the greatest effect, with about half of all the VOC analyzed affected, followed by time of milking (evening milking vs. morning milking) and mixing (skim evening milk mixed with whole morning milk). In general, among fresh products, cream had higher contents of fatty acids, sulfurs, and terpene volatile compounds than fresh cheese and ricotta, whereas ricotta showed a very high VOC amount compared with fresh cheese, probably due to its high processing temperature. The effects of the progressive nutrient depletion in milk during processing were investigated by comparing the amounts of VOC in vat milk, whey, and scotta. Although milk contained greater amounts of nutrients, whey and especially scotta had higher concentrations of VOC, with the exception of esters, sulfurs, terpenes, and phenolic compounds, as a result of physicochemical and microbial modifications during processing. Finally, the effect of ripening was tested by comparing the VOC of fresh and ripened cheeses (6 and 12 mo), revealing that VOC release increased dramatically during the first semester and further with increasing the ripening period to 1 yr. In particular, some alcohols (butan-2-ol), aldehydes (2-methylpropanal, hexanal, and heptanal), esters (ethyl butanoate and ethyl hexanoate), fatty acids (acetic, butanoic, and hexanoic acids), and ketones (butan-2-one, pentan-2-one, and heptan-2-one) showed a very large increase. In conclusion, according to the artisanal milk processing carried out for Malga cheese production, the quantity of VOC was shown to increase about 3 times during cheese making (from milk in vat to fresh cheese plus whey), almost 4 times during ricotta making (from whey to ricotta plus scotta), and about 16 times during 1 yr of ripening of cheese.
Milk coagulation properties (MCP) are conventionally measured using computerized renneting meters, mechanical or optical devices that record curd firmness over time (CFt). The traditional MCP are ...rennet coagulation time (RCT, min), curd firmness (a30, mm), and curd-firming time (k20, min). The milk of different ruminant species varies in terms of CFt pattern. Milk from Holstein-Friesian and some Scandinavian cattle breeds yields higher proportions of noncoagulating samples, samples with longer RCT and lower a30, and samples for which k20 is not estimable, than does milk from Brown Swiss, Simmental, and other local Alpine breeds. The amount, proportion, and genetic variants (especially κ-casein) of milk protein fractions strongly influence MCP and explain variable proportions of the observed differences among breeds and among individuals of the same breed. In addition, other major genes have been shown to affect MCP. Individual repeatability of MCP is high, whereas any herd effect is low; thus, the improvement of MCP should be based principally on selection. Exploitable additive genetic variation in MCP exists and has been assessed using different breeds in various countries. Several models have been formulated that either handle noncoagulating samples or not. The heritability of MCP is similar to that of other milk quality traits and is higher than the heritability of milk yield. Rennet coagulation time and a30 are highly correlated, both phenotypically and genetically. This means that the use of a30 data does not add valuable information to that obtainable from RCT; both traits are genetically correlated mainly with milk acidity. Moreover, a30 is correlated with casein content. The major limitations of traditional MCP can be overcome by prolonging the observation period and by using a novel CFt modeling, which uses all available information provided by computerized renneting meters and allows the estimation of RCT, the potential asymptotic curd firmness, the curd-firming rate, and the syneresis rate. Direct measurements of RCT obtained from both mechanical and optical devices show similar heritabilities and exhibit high phenotypic and genetic correlations. Moreover, mid-infrared reflectance spectroscopy can predict MCP. The heritabilities of predicted MCP are higher than those of measured MCP, and the 2 sets of values are strongly correlated. Therefore, mid-infrared reflectance spectroscopy is a reliable and cheap method whereby MCP can be improved at the population level; this is because such spectra are already routinely acquired from the milk of cows enrolled in milk recording schemes.
The aim of this study was to investigate associations between pathogen-specific cases of subclinical mastitis and milk yield, quality, protein composition, and cheese-making traits. Forty-one ...multibreed herds were selected for the study, and composite milk samples were collected from 1,508 cows belonging to 3 specialized dairy breeds (Holstein Friesian, Brown Swiss, and Jersey) and 3 dual-purpose breeds of Alpine origin (Simmental, Rendena, and Grey Alpine). Milk composition i.e., fat, protein, casein, lactose, pH, urea, and somatic cell count (SCC) was analyzed, and separation of protein fractions was performed by reversed-phase high performance liquid chromatography. Eleven coagulation traits were measured: 5 traditional milk coagulation properties time from rennet addition to milk gelation (RCT, min), curd-firming rate as the time to a curd firmness (CF) of 20 mm (k20, min), and CF at 30, 45, and 60 min from rennet addition (a30, a45, and a60, mm), and 6 new curd firming and syneresis traits potential asymptotical CF at an infinite time (CFP, mm), curd-firming instant rate constant (kCF, % × min−1), curd syneresis instant rate constant (kSR, % × min−1), modeled RCT (RCTeq, min), maximum CF value (CFmax, mm), and time at CFmax (tmax, min). We also measured 3 cheese yield traits, expressing the weights of total fresh curd (%CYCURD), dry matter (%CYSOLIDS), and water (%CYWATER) in the curd as percentages of the weight of the processed milk, and 4 nutrient recovery traits (RECPROTEIN, RECFAT, RECSOLIDS, and RECENERGY), representing the percentage ratio between each nutrient in the curd and milk. Milk samples with SCC > 100,000 cells/mL were subjected to bacteriological examination. All samples were divided into 7 clusters of udder health (UH) status: healthy (cows with milk SCC < 100,000 cells/mL and uncultured); culture-negative samples with low, medium, or high SCC; and culture-positive samples divided into contagious, environmental, and opportunistic intramammary infection (IMI). Data were analyzed using a linear mixed model. Significant variations in the casein to protein ratio and lactose content were observed in all culture-positive samples and in culture-negative samples with medium to high SCC compared to normal milk. No differences were observed among contagious, environmental, and opportunistic pathogens, suggesting an effect of inflammation rather than infection. The greatest impairment in milk quantity and composition, clotting ability, and cheese production was observed in the 2 UH status groups with the highest milk SCC (i.e., contagious IMI and culture-negative samples with high SCC), revealing a discrepancy between the bacteriological results and inflammatory status, and thus confirming the importance of SCC as an indicator of udder health and milk quality.
Fourier-transform infrared (FTIR) spectra are used to predict the fat, protein, casein, and lactose contents of milk. These estimates are currently used to predict the individual estimated breeding ...values of animals. The objective of the present study was to estimate the genetic variation and heritabilities of the milk transmittance spectrum at each individual FTIR wave. Milk was sampled once per cow from a total of 1,064 Italian Brown Swiss cows from 30 herds, sired by 50 artificial insemination sires. The FTIR spectra of all samples were collected within 3h of sampling from 25mL of milk. The obtained spectral range comprised wavenumbers 5,000 to 930×cm−1, corresponding to wavelengths 2.00 to 10.76μm and frequencies from 149.9 to 27.9THz, for a total of 1,056 waves. These were acquired using a MilkoScan FT120 FTIR interferometer (Foss Electric A/S, Hillerød, Denmark). Each spectral data point was treated as a single trait and analyzed using an animal model REML method. The results indicated that the transmittance of the bovine milk FTIR spectrum was heritable for most individual waves in the wavenumber interval from 5,000 to 930×cm−1. Moreover, the transmittance of contiguous FTIR waves was much more highly correlated in terms of the average value and phenotypic variation, compared with genetic variation. In the present study, we characterized 5 regions of the FTIR spectrum that were relevant to the analysis of milk; 2 regions, one in the transition area between the short-wavelength infrared (SWIR) and mid-wavelength infrared (MWIR) divisions of the electromagnetic spectrum (SWIR-MWIR region) and another very short region in the MWIR division (MWIR-2 region), were characterized by very high phenotypic variability in the transmittance of individual milk samples within each wave. This was caused by the absorption peaks of water, which can mask the effects of other important milk components. These regions also showed high genetic variability in transmittance, and the heritability estimates of individual waves were generally very low (with some exceptions). The 3 other identified regions contained many transmittance peaks that represented important chemical bonds; these showed much lower phenotypic and genetic variability in terms of individual waves, but relatively higher and less variable heritability estimates. Among them, the SWIR region (near-infrared) showed a peculiar cyclic pattern of the heritability coefficients of transmittance, the MWIR-1 region was particularly important for the estimation of fat, and the MWIR-LWIR region (also known also as the “fingerprint region”) had 3 areas of relatively high heritability. In summary, we found that the transmittance data from the FTIR spectra of milk have genetic variability that may prove useful for the direct genetic improvement of dairy species, rather than only through indirect phenotypic predictions of individual milk quality and technological traits.
The list of standard abbreviations for JDS is available at adsa.org/jds-abbreviations-24. Nonstandard abbreviations are available in the Notes.
This study aimed to set up a life cycle assessment ...(LCA) approach at the level of individual animals to assess the effects of a 3-breed crossbreeding program on the environmental impact of cows. The study involved 564 cows, 279 purebred Holstein Friesian (HO) and 285 crossbred cows (CR), that originated from a 3-breed crossbreeding program based on the rotational use of Viking Red, Montebèliarde, and HO sires and kept in 2 dairy herds of northern Italy (224 and 340 cows/herd, respectively). The reference unit of the LCA model was the lifetime of cows, from the birth to culling or death. Data were collected at different levels: individual animal-based data referred to the whole life (birth, calving, dry, cull or death dates, and milk production); individual test-date collection of body measures and BCS, used to predict BW and to estimate energy requirements; common farm-based data concerning herd management (diets composition, and materials used). Data were used to compute DMI, milk and milk components production, gross income (GI), and income over feed costs (IOFC) pertaining to the lifespan of cows. An individual LCA-derived approach was set up to compute global warming potential (GWP), acidification and eutrophication potential (AP and EP, respectively), and land occupation (LO), which have been associated with different functional units (cow in her whole life or per day of life; kilograms of milk fat plus protein, and GI and IOFC in euros produced in the herd life). Data were analyzed using a generalized linear model including the fixed effects of genetic group (CR vs. HO), farm, and their interaction (genetic group × farm). Compared with HO, CR cows completed more lactations (+12%), had earlier first calving (−2 wk), yielded more fat plus protein in milk both in the lifespan (+8%) and per day of life (+4%). Concerning the environmental impact, when compared with HO herd mates, CR cows had nominal greater emissions per cow in the whole life, similar emissions per day of life and ∼3% lower GWP, AP, and EP per kilogram of fat plus protein yielded in lifespan. Income over feed costs per unit of emission tended to be ∼4% greater in CR compared with HO cows. Also, the use of land tended to be lower in CR compared with HO in most indicators considered. In conclusion, LCA could be adapted to represent individual animals. Moreover, managing dairy cows according to a 3-breed rotational crossbreeding scheme may be regarded as a strategy that can contribute to mitigate the emissions and to improve the environmental impact of dairy operations.
Culled dairy cows represent a considerable source of meat production, but their carcasses may vary greatly in quality because of the wide variation in the age, stage of lactation, breed, body ...condition, and other characteristics of the cows at slaughter. However, the effect of crossbreeding on the value of culled cows has so far received little investigation. The aim of this observational study was to compare a range of carcass attributes of cull cows from 3-breed rotational crossbreeding using Viking Red, Montbéliarde (MO), and Holstein (HO) bulls with those of HO purebred cows. Data on 1,814 dairy cows were collected. Cows were reared together in one herd and slaughtered in 4 slaughterhouses. The carcass weight, fleshiness, and fatness scores, the total value, and the price (€/kg) of each cow carcass were recorded. The culling of a few cows in the sample (n = 86) was classified by the farm manager as “urgent” following a diagnosis of injury or sickness, and this information was recorded. Carcass traits were analyzed with a mixed model which included the fixed effects of parity, days in milk, genetic group (purebred HO, 787 cows, and crossbred cows, classified according to the breed of sire within crossbreds, with 309, 428, and 290 cows sired by Viking Red, MO, and HO bulls, respectively), and interactions, and the random effects of month × year of the date of slaughter, and slaughterhouse. Logistic regression was used to investigate the association of parity, days in milk and purebred or crossbred origin with unplanned, “urgent” culling compared with regular culling. Average carcass weight across genetic groups was 297 ± 65 kg, average price €2.03 ± 0.53/kg, and average value €631 ± 269. Compared with HO, crossbred carcasses were 7 to 12% heavier depending on the breed of sire, were graded + 0.12 to + 0.28 units higher for fleshiness and + 0.26 to + 0.30 units higher for fatness, and fetched an 8 to 11% higher price. As a consequence, compared with purebred HO, carcasses from crossbreds had 15 to 24% higher value (€84 to €133 more per cow), with crossbred cows sired by MO showing the greatest values. Moreover, compared with the HO cows, the crossbred cows had a 37% lower risk of being urgently removed from the herd, which raises welfare concerns and may reduce the salvage value of cull cows. Because cull cows represent a supplemental source of income for dairy farmers, the greater overall value of crossbred cull cows should be taken into account in evaluating the economic effectiveness of crossbreeding schemes.
Crossbreeding is a strategy to counter the declining fertility, resilience, and longevity of purebred Holstein (Ho) cows. However, little is known of the effects of long-term systematic rotational ...crossbreeding on milk technological properties and cheese yield (CY). In this study, we compared the individual milk composition, milk coagulation properties (MCP), and CY of 468 purebred Ho and 648 crossbred (CR) cows obtained from two 3-breed rotational crossbreeding systems using Viking Red (VR), Montbéliarde (Mo), and Ho sires over 4 generations. Individual milk samples were collected once from 1,116 primiparous and multiparous cows kept in 2 dairy herds, raised for the production of Grana Padano (high milk yield, total mixed ration based on corn silage) and Parmigiano Reggiano (moderate milk yield, only dry feeds) cheeses. In both herds, a 3-breed rotational mating system was used in which Ho cows were first inseminated with VR, whereas Mo and Ho semen was used in the subsequent generations. In one herd, the sequence Mo-VR-Ho was also used. Individual milk samples were analyzed for milk composition, single-point MCP, and parameters for modeling curd firming over time, whereas CY and milk nutrient recovery in the curd were assessed through a laboratory cheese-making procedure. Compared with Ho, CR cows produced 5.8% less milk, which had comparable fat but greater protein and casein contents and lower lactose contents and somatic cell scores. Milk from CR cows tended to reach a curd firmness of 20 mm more quickly and exhibited greater curd firmness at 45 and 60 min from rennet addition. Holstein and CR cows yielded milk with similar CY and recovery in the curd traits. The milk fat content, somatic cell scores, curd firmness traits, and CY of CR cows relative to the Ho cows differed in the 2 herds, and the favorable effects on the CR cows were more evident in the herd with the greatest milk yield and the worst MCP traits. Crossbred cows of the 4 generations performed similarly, with the exception of the better MCP of the milk from first-generation CR cows. The 2 rotational systems using different sire-breed sequences also performed similarly. In summary, both rotational crossbreeding programs exhibited some advantage over the Ho purebred breeding system in terms of milk composition and MCP but not CY. Future research is needed to investigate the interactions between crossbreeding schemes and dairy systems.
The aim of this study was to elucidate the relationships between various cheesemaking-related traits, namely the well-known traditional milk coagulation properties (MCP), the new curd firming and ...syneresis traits, the cheese yield, and the curd nutrient recoveries or whey losses (all measured at the individual level). Data were obtained from 1,167 Brown Swiss cows reared in 85 herds. A 2-L milk sample was collected once from each animal and assessed for 10 phenotypes related to changes in curd firmness (CF) over time, plus 7 cheesemaking traits. The CF-related traits included 4 traditional single-point lactodynamographic properties rennet coagulation time (RCT, min); time to a CF of 20mm, min; and the CF 30 and 45min after rennet addition (a30 and a45, respectively), 4 parameters used to model the 360 CF data recorded over time for each milk sample the potential asymptotic CF at infinite time (CFP, mm); the CF instant rate constant, % × min−1; the syneresis instant rate constant, % × min−1; and the RCT obtained from modeling individual samples, and 2 traits calculated from individual equations the maximum CF(CFmax, mm); and the time at CFmax, min. The cheesemaking traits included 3 cheese yield traits (weights of the fresh curd, curd solids and curd moisture as percent of the weights of the processed milk) and 4milk nutrient recoveries in the curd (calculated as the percent ratios between a given nutrient in the curd versus that in the processed milk). Bayesian methodology-based multivariate analyses were used to estimate the phenotypic, additive genetic, herd/date, and residual relationships between the aforementioned traits, whereas statistical inferences were based on the marginal posterior distributions of the parameters of concern. The a45, CFP, and CFmax traits were genetically associated with all of the percent cheese yield traits (the additive genetic correlations varied from 0.752 to 0.855 for a45; 0.496 to 0.583 for CFP; and 0.750 to 0.801 for CFmax) and the nutrient recovery traits (additive genetic correlations varied from 0.296 to 0.901 for a45; 0.428 to 0.697 for CFP; and 0.412 to 0.941 for CFmax). Moreover, the nutrient recoveries for fat, solids, and energy exhibited large additive genetic correlations with the other coagulation and curd firming traits. In particular, recovery of protein and fat were found to be powerful instruments for understanding the relationships between milk technological properties and cheese quantity or quality. We observed only weak genetic relationships with the milk quality and MCP traits, suggesting that the highly heritable trait of protein recovery should perhaps be included as a genetic index when seeking to improve cheesemaking efficiency at the population level. In contrast, we found that fat recovery exhibited moderate genetic variation and could be improved through the CF over time traits, especially using those recorded during the late phase of the curd firming process. Moreover, our results demonstrated that the traditional MCP have limited relevance for predicting individual cheese yield. Therefore, their use for this purpose in the dairy industry and breeding programs seems questionable.
In the present study we estimated the genetic parameters of enteric methane emissions (EME) traits predicted from milk fatty acid profile (FA) and those of their predictors in 1,091 Brown Swiss cows ...reared on 85 farms in order to assess the potential of using EME-related phenotypes in selective breeding. Univariate and bivariate genetic models were fitted in a Bayesian framework. The means of the marginal posterior distribution of intra-herd heritability ranged from 0.12 for estimated methane production (g/d/cow) to 0.24 for estimated methane yield (g/kg dry matter intake DMI), with intermediate values for estimated methane intensity, increasingly higher when expressed per kg of corrected milk (0.13), fresh cheese (0.16), or cheese solids (0.20). Regarding the correlations, the milk quality traits and percentage cheese yields were generally moderately correlated with the estimated EME traits, and were variable in terms of sign. Daily milk and cheese yield traits were, as expected, all highly positively correlated with estimated daily methane production. In contrast, they were negatively correlated with estimated methane yield and intensity, the estimates being large in the case of phenotypic and herd correlations, and low in the case of additive genetic and residual correlations. With the exception of the negative correlations with daily methane production, EME traits exhibited trivial correlations with body size and BCS of cows, which, in turn, were negatively correlated with milk yield. Although the results should be validated on a larger population and different breeds, our study demonstrate the presence of additive genetic variation of EME traits, which could be exploited in breeding programmes for the improvement in both milk production and the ecological footprint of dairy farming.
Highlights
Enteric methane emissions (EME) of dairy cows can be estimated on the basis of milk fatty acid profile.
EME exhibited exploitable genetic variation.
Genetic selection could be preferentially based on predicted methane intensity per kg of milk, or per kg of cheese in countries where milk production is used mainly for cheese-making.