Selection for Test-Day Milk Yield and Thermotolerance in Brazilian Holstein Cattle

Unravelling the impact of heat stress on daughter pregnancy rate in Mehsana buffalo through innovative breeding interventions

Heat stress profoundly affects the reproductive success of buffaloes, which are vital for the dairy industry due to their unique anatomical and physiological characteristics, necessitating careful evaluation under such conditions. Hence, this guided our search for quantifying heat stress' impact on Mehsana buffaloes using the best THI model and evaluating sires' performance. Fertility records (days open converted to daughter pregnancy rate) were collected in the span of over 24 years, w.e.f. 1989 to 2012. Finally, 3070 records of first lactation cows, daughters of 117 sires from DURDA, Gujarat, India, were used in the analysis. Meteorological data were retrieved from IMD, Pune, to understand the relationship between daughter pregnancy rate (DPR) and heat stress indicators. Several heat stress models were compared based on R2, adjusted R2, AIC, and BIC values, and the impact of heat stress was quantified. The year was classified into different heat stress zones, viz., Non heat stress zone (NHSZ), Heat stress zone (HSZ), and critical heat stress zone (CHSZ), drawing from the findings of DPR and THI. The THI 4th model based on dry and wet bulb temperature was identified as the best-fit model, and DPR significantly changed (P < 0.01) by 1.14% per unit change in THI value based on the month of calving. The average EBVs of the sires for DPR were found to be 20.78% (NHSZ), 38.09% (HSZ), and 39.08% (CHSZ) using BLUP-SM and 20.78% (NHSZ), 37.30% (HSZ), and 38.87% (HSZ) using BLUP-AM. Subsequently, the optimum sire for each of the zones was prioritized. It is noteworthy that bulls that performed better in NHSZ did not perform as well in HSZ and CHSZ, and vice versa. This supports the possibility of evaluating bulls independently in each heat stress zone.

Clustering of countries based on dairy productivity characteristics of Holstein cattle for breeding material selection

Background and Aim

The aim of any breeding process is to create a herd based on certain parameters that reflect an ideal animal vision. Targeted herding involves selecting the source of breeding material to be imported from another country. Therefore, there is a problem in selecting a breeding material importer to rapidly form a uterine canopy with the required properties. The purpose of this study was to evaluate a set of predictive milk productivity traits in Holstein cattle across countries.

Materials and Methods

This research was based on records of 819,358 recorded animals from 28 countries born after January 1, 2018, from open databases. We used the Euclidean metric to construct dendrograms characterizing the similarity of countries according to the complex milk productivity traits of the daughters of bulls. The Ward method was used to minimize intracluster variance when forming clusters and constructing the corresponding diagrams. Principal component analysis was used to reduce dimensionality and eliminate the effect of multicollinearity. The principal components were selected using the Kaiser-Harris criteria.

Results

A ranking of multidimensional complex milk productivity traits in different countries over the past 5 years was performed. A group of leading countries led by the USA was established according to the studied indicators, and the possible reasons for such a division into groups were described.

Conclusion

The pressure of purposeful artificial selection prevails in comparison with the pressure of natural selection concerning milk productivity traits in a certain group of countries, which allows specialists to choose suppliers when buying breeding animals and materials. The findings are based solely on data from recorded animals, which may not represent the entire breed population within each country, especially in regions where record-keeping may be inconsistent. It is expected that further studies will include regional data from large enterprises not part of Interbull, with mandatory verification and validation. An important element of such work is seen as the ability to compare the milk productivity of populations from different countries using a different scale, as well as studying the differentiation of countries by other selection traits of dairy.

Impact of heat stress on genetic evaluation of oocyte and embryo production in Gir dairy cattle

Identifying and selecting genotypes tolerant to heat stress might improve reproductive traits in dairy cattle, including oocyte and embryo production. The temperature-humidity index (THI) was used, via random regression models, to investigate the impact of heat stress on genetic parameters and breeding values of oocyte and embryo production in Gir dairy cattle. We evaluated records of total oocytes (TO), viable oocytes (VO), cleaved embryos (CE), and viable embryos (VE) from dairy Gir donors. Twenty-four models were tested, considering age at ovum pick-up (AOPU) and THI means as a regressor in the genetic evaluation. We computed THI in eight periods, from 0 to 112 days before ovum pick-up, which were adjusted by different orders of Legendre polynomials (second, third, and fourth). The best-fit model according to Akaike's information criterion (AIC) and Model Posterior Probabilities (MPP) considered Legendre polynomials of third order and THI means of 112 days for TO, fourth order and 56 days for VO, second order and 28 days for CE, and second order and 42 days for VE, respectively. The heritability (h2) estimates across AOPU and THI scales ranged from 0.34 to 0.62 for TO, 0.31 to 0.58 for VO, 0.26 to 0.39 for CE, and 0.15 to 0.26 for VE, respectively. The fraction of the phenotypic variance explained by the permanent environment in different AOPU and THI scales ranged from 0.03 to 0.25 for TO, 0.05 to 0.26 for VO, 0.09 to 0.36 for CE, and 0.15 to 0.27 for VE, respectively. Spearman's rank correlation between the estimated breeding values in different AOPU and THI scale from the top 5% sires and females ranged from 0.18 to 0.90 for TO, 0.31 to 0.95 for VO, 0.14 to 0.85 for CE, and 0.47 to 0.94 for VE, respectively. The h2 estimates for all evaluated traits varied from moderate to high magnitude across AOPU and THI scales, indicating that genetic selection can result in rapid genetic progress for the evaluated traits. There was a reranking among the best animals in different AOPU and THI. It is possible to select dairy Gir cattle tolerant to heat stress to improve oocyte and embryo production.

How does reproduction account for dairy farm sustainability?

Sustainability - the new hype of the 21st century has brought discomfort for the government and society. Sustainable agriculture is essential to face our most concerning challenges: climate change, food security, and the environmental footprint, all of which add to consumers' opinions and choices. Improvements in reproductive indexes can enhance animal production and efficiency, guaranteeing profit and sustainability. Estrus detection, artificial insemination (AI), embryo transfer (ET), estrus synchronization (ES), and multiple ovulations are some strategies used to improve animal reproduction. This review highlights how reproductive strategies and genetic selection can contribute to sustainable ruminant production. Improved reproductive indices can reduce the number of nonproductive cows in the herd, reducing methane emissions and land use for production while preserving natural resources.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

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.

Genes and models for estimating genetic parameters for heat tolerance in dairy cattle

Dairy cattle are highly susceptible to heat stress. Heat stress causes a decline in milk yield, reduced dry matter intake, reduced fertility rates, and alteration of physiological traits (e.g., respiration rate, rectal temperature, heart rates, pulse rates, panting score, sweating rates, and drooling score) and other biomarkers (oxidative heat stress biomarkers and stress response genes). Considering the significant effect of global warming on dairy cattle farming, coupled with the aim to reduce income losses of dairy cattle farmers and improve production under hot environment, there is a need to develop heat tolerant dairy cattle that can grow, reproduce and produce milk reasonably under the changing global climate and increasing temperature. The identification of heat tolerant dairy cattle is an alternative strategy for breeding thermotolerant dairy cattle for changing climatic conditions. This review synthesizes information pertaining to quantitative genetic models that have been applied to estimate genetic parameters for heat tolerance and relationship between measures of heat tolerance and production and reproductive performance traits in dairy cattle. Moreover, the review identified the genes that have been shown to influence heat tolerance in dairy cattle and evaluated the possibility of using them in genomic selection programmes. Combining genomics information with environmental, physiological, and production parameters information is a crucial strategy to understand the mechanisms of heat tolerance while breeding heat tolerant dairy cattle adapted to future climatic conditions. Thus, selection for thermotolerant dairy cattle is feasible.

Changes in the pattern of heat waves and the impacts on Holstein cows in a subtropical region

This study aimed to evaluate the change in the air temperature and the impacts of heat waves using Climate Change Indexes on the physiological and productive responses of lactating Holstein cows. Daily data of maximum and minimum air temperature for 1981-2021 were used. Heat waves were determined using six Climate Change Indexes. Individual data on respiratory rate, rectal temperature, and milk yield were collected in the summers of 2018, 2019, and 2021. The temperature trend analysis showed a significant (p < 0.0001) increase in maximum temperature, minimum temperature, and days in a heat wave. All six indexes increased significantly (p > 0.01). The increase in warm nights (> 20 °C) and the hottest days (> 35 °C) was the highest since 2010. Heat waves were classified into short (< 5 days) and long (> 5 days) of greater (> 36 °C) or lesser (< 36 °C) intensity. During the long and short heat waves of greater intensity, the respiratory rate increased (p < 0.05) until the fourth day. On the other hand, rectal temperature was higher (p < 0.05) from the fourth day until the end of the long heat waves. Therefore, the decrease in milk yield was significantly greater from the fourth or fifth day onwards. Finally, the evaluation method based on indexes was efficient to demonstrate the negative effects on physiological parameters and milk yield and can be indicated to evaluate heat stress in lactating cows.

Heteroscedastic Reaction Norm Models Improve the Assessment of Genotype by Environment Interaction for Growth, Reproductive, and Visual Score Traits in Nellore Cattle

The assessment of the presence of genotype by environment interaction (GxE) in beef cattle is very important in tropical countries with diverse climatic conditions and production systems. The present study aimed to assess the presence of GxE by using different reaction norm models for eleven traits related to growth, reproduction, and visual score in Nellore cattle. We studied five reaction norm models (RNM), fitting a linear model considering homoscedastic residual variance (RNM_homo), and four models considering heteroskedasticity, being linear (RNM_hete), quadratic (RNM_quad), linear spline (RNM_l-l), and quadratic spline (RNM_q-q). There was the presence of GxE for age at first calving (AFC), scrotal circumference (SC), weaning to yearling weight gain (WYG), and yearling weight (YW). The best models were RNM_l-l for YW and RNM_q-q for AFC, SC, and WYG. The heritability estimates for RNM_l-l ranged from 0.07 to 0.20, 0.42 to 0.61, 0.24 to 0.42, and 0.47 to 0.63 for AFC, SC, WYG, and YW, respectively. The heteroskedasticity in reaction norm models improves the assessment of the presence of GxE for YW, WYG, AFC, and SC. Additionally, the trajectories of reaction norms for these traits seem to be affected by a non-linear component, and selecting robust animals for these traits is an alternative to increase production and reduce environmental sensitivity.

Genotype by Environment Interaction and Selection Response for Milk Yield Traits and Conformation in a Local Cattle Breed Using a Reaction Norm Approach

Local breeds are often reared in various environmental conditions (EC), suggesting that genotype by environment interaction (GxE) could influence genetic progress. This study aimed at investigating GxE and response to selection (R) in Rendena cattle under diverse EC. Traits included milk, fat, and protein yields, fat and protein percentage, and somatic cell score, three-factor scores and 24 linear type traits. The traits belonged to 11,085 cows (615 sires). Variance components were estimated in a two-step reaction norm model (RNM). A single trait animal model was run to obtain the solutions of herd-EC effect, then included in a random regression sire model. A multivariate response to selection (R) in different EC was computed for traits under selection including beef traits from a performance test. GxE accounted on average for 10% of phenotypic variance, and an average rank correlation of over 0.97 was found between bull estimated breeding values (EBVs) by either including or not including GxE, with changing top ranks. For various traits, significantly greater genetic components and R were observed in plain farms, loose housing rearing system, feeding total mixed ration, and without summer pasture. Conversely, for beef traits, a greater R was found for mountain farms, loose housing, hay-based feeding and summer pasture.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Direct-Contact Air, Gravel, Ground Heat Exchanger in Air Treatment Systems for Cowshed Air Conditioning

This paper describes the analysis of the possibility of use of the direct-contact air, gravel, ground heat exchanger (acronym GAHE), patented at the Wroclaw University of Science and Technology, as a means of improving microclimate parameters in dairy cows’ barns. Different possibilities of introducing GAHE to the standard mechanical ventilation system of cowsheds have been proposed and investigated. Based on literature data, the required air parameters in the barns of dairy cows were determined and discussed. Computer simulations were carried out and the results obtained were compared to the baseline model. Year-round changes in microclimate parameters, especially air temperature, relative humidity, and THI index were investigated. The benefits of GAHE use were indicated. The possible increase in the minimum air volume of ventilation during the winter season and the decrease in the maximum values of this parameter in the summer were presented. Indications were made of the systems where the application of GAHE could be the most beneficial. A further research path has been proposed.