Estimation of maximum thermo-hygrometric index thresholds affecting milk production in Italian Brown Swiss cattle

Relationship among thermal environment, stage of lactation, body characteristics, yield and milk constituents of dairy Gyr cows managed on pasture

Our aims were to evaluate changes in body characteristics, milk yield and milk constituents as well as to determine the relationship between the thermal environment and production characteristics during the first lactation of dairy Gyr cows managed on pasture. Between 2013 and 2015, forty-five primiparous dairy Gyr cows were evaluated from prepartum to 10 months of lactation in Southeast of Brazil. Body weight, body condition score (BCS), subcutaneous fat thickness (SFT), milk yield (305 d), and milk constituents were collected monthly and progesterone was collected weekly. Additionally, we determined the temperature humidity index (THI) based on microclimate data. Overall, the cows lost body weight until six months of lactation and there was a progressive decrease in BCS, SFT, milk yield and milk lactose as the months in lactation progressed. In contrast, there was an increase in milk fat, milk protein and milk solids. The thermal environment did not pose a consistent heat challenge, nevertheless, we found a positive correlation between the average THI two days before milk collection with milk yield, fat and lactose contents, but in contrast a negative correlation was found with total solids and protein. In conclusion, the THI and months of lactation affected the yield and constituents of milk. However, more studies are necessary to understand the impacts of body characteristics and thermal environment on yield and milk constituents throughout the productive life of Gyr dairy cows.

Effect of transgenerational environmental condition on genetics parameters of Italian Brown Swiss

The aim of this study was to infer the effects of heat stress (HS) of dams during late gestation on direct and maternal genetic parameters for traits related to milk production and milk quality parameters (90,558 records) in Italian Brown Swiss cattle (12,072 cows in 617 herds). Daily average temperature-humidity indices (THI) during the last 56 d of pregnancy were calculated, using the climate data from the nearest public weather station for each herd. Heat load effects were considered as the average across the entire periods considering a thermoneutrality condition for data below the THI 60. For parameter estimation a random regression model using the second-order Legendre polynomial regression coefficient for THI considering both animal and maternal effect for heat load. Direct heritability increased sharply from THI 60 to 65, then decreased gradually up to THI ∼72, and sharply thereafter. Maternal heritability showed a different trend, with values close to 0 up until to THI 65 and slightly increasing toward extreme THI values. The study suggests a lower threshold of THI 60 for the onset of HS. Higher heritability values indicate greater selective efficiency in the THI range of 65 to 70, even if a higher standard deviation value have been detected. The effects of high THI during intrauterine life varied among traits with different heritability levels. Genetic correlations for milk, fat, and protein content at 60 THI with increasing value of environmental variable, remained constant (∼0.90) until THI >75, where they slightly decreased (∼0.85). Fat and protein yields, as well as milk and energy-corrected milk, showed correlations dropping to 0.80 around THI 67 to 68 and stabilizing between 0.75 and 0.85 at extreme THI values. Maternal component correlations dropped close to zero, with negative values for protein content at THI 65 to 70. Antagonism between direct and maternal components was stronger for intermediate THI values but less divergent for extremes. Genotype by environment interaction was observed, indicating the selection of resilient animals would be theoretically possible. In the future, the application of climate variables in selection schemes first should take into account the dimensions of the genetic correlations to be able to decide between the simple inclusion of the environmental effect in the statistical models, rather than a real parallel genetic evaluation.

Investigation of Climate Effects on the Physiological Parameters of Dairy Livestock (Cow vs. Buffalo)

Nowadays climate change is affecting the planet's biodiversity, and livestock practices must adapt themselves to improve production without affecting animal welfare. This work investigates the influence that some climatic parameters such as Environment Temperature, Relative Humidity, Thermal excursion and Temperature-Humidity Index (THI), can have on milk quantity and quality in two different dairy species (buffaloes and cows) raised on the same farm. A further aim was to understand if THI threshold used for cows could also be used for buffaloes. The climatic parameters were recorded daily through a meteorological station located inside the farm. Milk quantity (converted into ECM) and quality (Fat Percentage-FP; Protein Percentage-PP; Somatic Cell Count-SCC) were measured. Data were analyzed with Spearman's correlation index, separately for buffaloes and cows. The results indicate a greater sensitivity of cows to heat stress and a strong negative correlation of the ECM with meteorological data (p < 0.01). The results of this study may stimulate the use of integrated technologies (sensors, software) in the dairy sector, since the IoT (sensors, software) helps to enhance animal well-being and to optimize process costs, with a precision livestock farming approach.

The autumn low milk yield syndrome in Brown Swiss cows in continental climates: hypotheses and facts

Extensive research has been conducted globally on the impact of heat stress (HS) on animal health and milk production in dairy cows. In this article, we examine the possible reasons for the decrease in milk production in Brown Swiss (BS) cows during the autumn season, known as the autumn low milk yield syndrome (ALMYS). This condition has been extensively studied in high-yielding Holstein Friesian (HF) cattle and has also been observed in BS cows with a daily milk yield of around 30 kg. Our hypothesis is that the drop in milk yield and the increased prevalence of mastitis in autumn, as found in our recent studies, may be a long-term consequence of summer HS. We re-evaluate our previous findings in light of the possible manifestation of an HS-related form of ALMYS in BS cows. As milk yield, mastitis spread, and reproductive function of cows are interrelated and have seasonal dependence, we examine the consistency of our hypothesis with existing data. The significant drop in milk yield in BS cows in autumn (by 2.0-3.2 kg), as well as the threshold of milk yield decrease (temperature-humidity index of 70.7), may point in favour of the manifestation of ALMYS in BS cows, similar to HF cows. Only the percentage effect of seasonal factor (59.4%; p < 0.05) on milk yield of BS cows was significant. HS-related ALMYS provides a robust conceptual framework for diverse sets of productive and animal health data in BS cows, similar to observations in high-yielding HF cattle. However, the limitations associated with the lack of additional data (e.g. immunological indicators) suggest the need for further research to confirm ALMYS in BS breed.

A system dynamics approach to model heat stress accumulation in dairy cows during a heatwave event

Climate change is expected to increase the number of heat wave events, leading to prolonged exposures to severe heat stress (HS) and the corresponding adverse effects on dairy cattle productivity. Modelling dairy cattle productivity under HS conditions is complicated because it requires comprehending the complexity, non-linearity, dynamicity, and delays in animal response. In this paper, we applied the System Dynamics methodology to understand the dynamics of animal response and system delays of observed milk yield (MY) in dairy cows under HS. Data on MY and temperature-humidity index were collected from a dairy cattle farm. Model development involved: (i) articulation of the problem, identification of the feedback mechanisms, and development of the dynamic hypothesis through a causal loop diagram; (ii) formulation of the quantitative model through a stock-and-flow structure; (iii) calibration of the model parameters; and (iv) analysis of results for individual cows. The model was successively evaluated with 20 cows in the case study farm, and the relevant parameters of their HS response were quantified with calibration. According to the evaluation of the results, the proposed model structure was able to capture the effect of HS for 11 cows with high accuracy with mean absolute percent error <5%, concordance correlation coefficient >0.6, and R2 > 0.6, except for two cows (ID #13 and #20) with R2 less than 0.6, implying that the rest of the nine animals do not exhibit heat-sensitive behaviour for the defined parameter space. The presented HS model considered non-linear feedback mechanisms as an attempt to help farmers and decision makers quantify the animal response to HS, predict MY under HS conditions, and distinguish the heat-sensitive cows from heat-tolerant cows at the farm level.

Marine Microalgae as a Nutritive Tool to Mitigate Ruminal Greenhouse Gas Production: In Vitro Fermentation Characteristics of Fresh and Ensiled Maize (Zea mays L.) Forage

The aim of the present study was to evaluate the effects of marine microalgae (Dunaliella salina) as a food additive on biogas (BG), methane (CH4), carbon monoxide (CO), and hydrogen sulfide (H2S) production kinetics, as well as in in vitro rumen fermentation and the CH4 conversion efficiency of different genotypes of maize (Zea mays L.) and states of forage. The treatments were characterized by the forage of five maize genotypes (Amarillo, Montesa, Olotillo, Tampiqueño, and Tuxpeño), two states of forage (fresh and ensiled), and the addition of 3% (on DM basis) of microalgae (with and without). The parameters (b = asymptotic production, c = production rate, and Lag = delay phase before gas production) of the production of BG, CH4, CO, and H2S showed an effect (p < 0.05) of the genotype, the state of the forage, the addition of the microalgae, or some of its interactions, except for the time in the CO delay phase (p > 0.05). Moreover, the addition of microalgae decreased (p < 0.05) the production of BG, CH4, and H2S in most of the genotypes and stages of the forage, but the production of CO increased (p < 0.05). In the case of fermentation characteristics, the microalgae increased (p < 0.05) the pH, DMD, SCFA, and ME in most genotypes and forage states. With the addition of the microalgae, the fresh forage from Olotillo obtained the highest pH (p < 0.05), and the ensiled from Amarillo, the highest (p < 0.05) DMD, SCFA, and ME. However, the ensiled forage produced more (p < 0.05) CH4 per unit of SFCA, ME, and OM, and the microalgae increased it (p < 0.05) even more, and the fresh forage from Amarillo presented the highest (p < 0.05) quantity of CH4 per unit of product. In conclusion, the D. salina microalga showed a potential to reduce the production of BG, CH4, and H2S in maize forage, but its effect depended on the chemical composition of the genotype and the state of the forage. Despite the above, the energy value of the forage (fresh and ensiled) improved, the DMD increased, and in some cases, SCFA and ME also increased, all without compromising CH4 conversion efficiency.

Different Conditions during Confinement in Pasture-Based Systems and Feeding Systems Affect the Fatty Acid Profile in the Milk and Cheese of Holstein Dairy Cows

The diet of dairy cows influences the fatty acid (FA) profiles of their milk and cheese, but how these are affected by different conditions during confinement in a mixed system (MS:grazing + total mixed ration:TMR) is not known. The aim of this study was to compare the FAs of the milk and cheese from MS in a compost-bedded pack barns (CB-GRZ) versus an outdoor soil-bedded pen (OD-GRZ) during confinement, and with a confinement system (100%TMR) in a compost-bedded pack barns (CB-TMR). Individual milk samples (n = 12 cows/group), cheese, and pooled milk (MilkP) samples were collected. The saturated FA percentages in the milk and the omega 6/omega 3 ratio in the MilkP and cheese were greater for the CB-TMR (p < 0.0001), while the unsaturated and monounsaturated FA percentages in the milk were lower for the CB-TMR than the MS (p < 0.001). The milk n-3, C18:3, and conjugated linoleic acid percentages were lower for the CB-TMR than the MS (p < 0.001). The milk n-3 and C18:3 were higher for the CB-GRZ than the OD-GRZ (p < 0.01), but no differences were observed between the MS in the MilkP and cheese. In conclusion, CB-GRZ cows during confinement produced better quality milk compared to OD-GRZ cows. However, the FA profiles of the milk, MilkP, and cheese were affected to a greater extent by the feeding management than by the conditions during confinement.

A transgenerational study on the effect of great-granddam birth month on granddaughter EBV for production traits in Italian Simmental cattle

Heat tolerance is a key feature of resilient animals. Offspring of animals that suffer environmental stress during pregnancy could show physiological, morphological, and metabolic modifications. This is due to a dynamic reprogramming of the epigenetics of the mammalian genome that occurs in the early life cycle. Thus, the aim of this study was to investigate the extent of the transgenerational effect of heat stress during the pregnancy of Italian Simmental cows. The effects of dam and granddam birth months (as indicator of pregnancy period) on their daughter and granddaughter estimated breeding values (EBV) for some dairy traits as well as of the temperature-humidity index (THI) during the pregnancy were tested. A total of 128,437 EBV (milk, fat, and protein yields, and somatic cell score) were provided by the Italian Association of Simmental Breeders. The best birth months (of both dam and granddam) for milk yield and protein yield were May and June, whereas the worst were January and March. Great-granddam pregnancies developed during the winter and spring seasons positively affected the EBV for milk and protein yields of their great-granddaughters; in contrast, pregnancies during summer and autumn had negative effects. These findings were confirmed by the effects of maximum and minimum THI in different parts of the great-granddam pregnancy on the performances of their great-granddaughters. Thus, a negative effect of high temperatures during the pregnancy of female ancestors was observed. Results of the present study suggest a transgenerational epigenetic inheritance in Italian Simmental cattle due to environmental stressors.

Genotype by environment interaction due to heat stress in Brown Swiss cattle

Due to its geographical position and a highly variable orography, Italy is characterized by several climatic areas and thus, by many different dairy cow farming systems. Brown Swiss cattle, in this context, are a very appreciated genetic resource for their adaptability and low metabolic requirement. The significant heterogeneity in farming systems may consist of genotype by environment (G × E) interactions with neglected changes in animals' rank position. The objective of this study was to investigate G × E for heat tolerance in Brown Swiss cattle for several production traits (milk, fat, and protein yield in kilograms; fat, protein, and cheese yield in percentage) and 2 derivate traits (fat-corrected milk and energy-corrected milk). We used the daily maximum temperature-humidity index (THI) range, calculated according to weather stations' data from 2008 to 2018 in Italy, and 202,776 test-day records from 23,396 Brown Swiss cows from 639 herds. Two different methodologies were applied to estimate the effect of the environmental variable (THI) on genetic parameters: (1) the reaction norm model, which uses a continuous random covariate to estimate the animal additive effect, and (2) the multitrait model, which splits each production pattern as a distinct and correlated trait according to the first (a thermal comfort condition), third (a moderate heat stress condition), and fifth (a severe heat stress condition) mean THI value quintile. The results from the reaction norm model showed a descending trend of the additive genetic effect until THI reached the value of 80. Then we recorded an increase with high extreme THI values (THI 90). Permanent environmental variance at increasing THI values revealed an opposite trend: The plot of heritability and the ratio of animal permanent environmental variance to phenotypic variance showed that when the environmental condition worsens, the additive genetic and permanent environmental component for production traits play a growing role. The negative additive genetic correlation between slope and linear random coefficient indicates no linear relationship between the production traits or under heat stress conditions, except for milk yield and protein yield. In tridimensional wireframe plots, the extreme margin decreases until a minimum of ∼0.90 of genetic correlation in the ECM trait, showing that the magnitude of G × E interaction is greater than the other traits. Genetic correlation values in Brown Swiss suggest the possibility of moderate changes in animals' estimated breeding value in heat stress conditions. Results indicated a moderate G × E interaction but significant variability in sire response related to their production level.

Estimation of genotype by temperature-humidity index interactions on milk production and udder health traits in Montbeliarde cows

BACKGROUND

Heat stress negatively influences cattle welfare, health and productivity. To cope with the forecasted increases in temperature and heat waves frequency, identifying high-producing animals that are tolerant to heat is of capital importance to maintain milk production. This study, based on the joint analysis of on-farm performance and weather data, had two objectives: (1) to determine the response in production performances (milk, fat and protein yields, fat and protein contents) and udder health (somatic cell score) to temperature-humidity index (THI) variations in Montbeliarde cows, and (2) to estimate the interactions between genotype and THI, to enable the identification of the most adapted animals for facing the expected increases in temperature.

RESULTS

Test-day records from first and second lactations from 2016 to 2020 were associated with the average THI during the three days before the test-day record. In total, 446,717 test-day records from 55,650 cows in first lactation and 457,516 test-day records from 58,229 cows in second lactation were analysed. The optimal THI was below 55 (i.e. ~ 12-13 °C) for all traits. Individual responses to THI were estimated by random regression models, which also included individual responses to days in milk. Regardless of the stage of lactation, genetic correlations along the THI gradient were above 0.80, which suggests that genotype-by-THI interactions were weak for production and udder health traits. Nevertheless, a variability in the individual slope of decay could be highlighted at high THI. The genetic correlation between production level at moderate THI and the slope at high THI was negative, while for somatic cell score, it was positive, indicating that heat stress amplifies the susceptibility to mastitis.

CONCLUSIONS

The optimal THI for French Montbeliarde cows is below 55 for production and udder health traits. Genetic-by-THI interactions are weak in French Montbeliarde cows for production and udder health traits, but not all animals react in the same way to high temperatures. Even if there is little room for improvement, using a heat tolerance index in cattle selection would be relevant to anticipate the expected increases in temperature. Further investigations are needed to interpret this variability on production traits. However, the current selection for mastitis resistance seems appropriate to adapt cattle to rising temperatures.

Effect of Heat Waves on Some Italian Brown Swiss Dairy Cows' Production Patterns

Climate change is impacting worldwide efficiency and welfare standards in livestock production systems. Considering the sensibility to heat stress reported for different milk production patterns in Italian Brown Swiss, this study aims to evaluate the effect of heat waves (HWs)of different lengths on some milk production traits (fat-corrected milk, energy-corrected milk, protein and fat yield, protein percentage, cheese production at 24 h, and cheese yield). A 10-year dataset (2009–2018), containing 202,776 test-day records from 23,296 Brown Swiss cows, was used. The dataset was merged both with the daily maximum temperature–humidity index (THI) recorded by weather stations and with the daily maximum THI threshold for each trait in Italian Brown Swiss cows. The study considered 4 different HWs according to their length: 2, 3, 4, and 5 consecutive days before the test-day over the weighted THI threshold. Milk production traits were determined as the difference in losses compared to those after only 1 day before the test-day over the weighted maximum THI. All traits showed to be affected by HWs. Particularly, protein percentage losses increased from −0.047% to −0.070% after 2 consecutive days over the daily THI threshold, reaching −0.10% to −0.14% after 5 days (p &lt; 0.01), showing a worsening trend with the increasing length of HWs. First parity cows showed to be more sensitive to HWs than other parity classes, recording greater losses after shorter HWs, compared to multiparous cows, for protein yield and, consequently, for cheese production at 24 h. This suggests a less efficient metabolic response to heat stress and exposure time in primiparous, compared to multiparous cows, probably due to their incomplete growth process that overlaps milk production, making it more difficult for them to dissipate heat. Although actions to mitigate heat stress are always needed in livestock, this study points out that often time exposure to warm periods worsens milk production traits in Brown Swiss cows.