Modelling lactation curve for milk fat to protein ratio in Iranian buffaloes (Bubalus bubalis) using non-linear mixed models

Comparison of nonlinear functions to describe lactation curves for cumulative milk production in buffalo

The aim of this study was to examine the suitability of different growth functions (linear, sinusoidal, Gompertz, Schumacher and Richards) to fit cumulative milk production data from buffalo cows. Cumulative milk production at each day in milk was calculated from two published datasets reporting (i) fortnightly test-day milk yield records of the first lactation of Murrah buffalo that had calved during 1977-2012 and (ii) the first lactation records of Jaffarabadi buffalo collected from history-cum-pedigree registers for each quinquennium between 1991 and 2010. Each function was fitted to the lactation curves using nonlinear regression procedures. The Richards and sinusoidal equations provided the smallest root mean square error values, Akaike's and Bayesian information criteria and, therefore, the best fit for the cumulative lactation curves for milk yield. The Richards equation appeared to provide the most accurate estimate of the cumulative milk production at peak milk yield. Sinusoidal and flexible classical growth functions are appropriate to describe cumulative milk production curves and estimate lactation traits in buffalo.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

Predicting live weight using body volume formula in lactating water buffalo

Live weight (LW) is an important piece of information within production systems, as it is related to several other economic characteristics. However, in the main buffalo-producing regions in the world, it is not common to periodically weigh the animals. We develop and evaluate linear, quadratic, and allometric mathematical models to predict LW using the body volume (BV) formula in lactating water buffalo (Bubalus bubalis) reared in southeastern Mexico. The LW (391.5 ± 138.9 kg) and BV (333.62 ± 58.51 dm³) were measured in 165 lactating Murrah buffalo aged between 3 and 10 years. The goodness-of-fit of the models was evaluated using the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), coefficient of determination (R²), mean-squared error (MSE) and root MSE (RMSE). In addition, the developed models were evaluated through cross-validation (k-folds). The ability of the fitted models to predict the observed values was evaluated based on the RMSEP, R², and mean absolute error (MAE). LW and BV were significantly positively and strongly correlated (r = 0.81; P < 0.001). The quadratic model had the lowest values of MSE (2788.12) and RMSE (52.80). On the other hand, the allometric model showed the lowest values of BIC (1319.24) and AIC (1313.07). The Quadratic and allometric models had lower values of MSEP and MAE. We recommend the quadratic and allometric models to predict the LW of lactating Murrah buffalo using BV as a predictor.

A meta-analysis of the genetic contribution estimates to major indicators for ketosis in dairy cows

The present study aimed to perform a meta-analysis using the random-effects model to merge published genetic parameter estimates for major indicators of ketosis [milk concentrations of acetone (ACETm) and β-hydroxybutyrate (BHBAm), and blood concentration of β-hydroxybutyrate (BHBAb)] in dairy cows. Overall, 51 heritability estimates and 130 genetic correlations from 19 papers published between 2012 and 2022 were used in this study. The average heritability estimates for ACETm, BHBAm, and BHBAb were 0.164, 0.123, and 0.141, respectively. The genetic correlation estimates between BHBAm and milk yield (MY), milk protein percentage (PP), and body condition score (BCS) were negative and moderate (-0.252, -0.200, and - 0.314, respectively). Genetic correlation estimates between BHBAm and milk fat percentage (FP), milk fat to protein ratio (FPR), and ketosis (KET) were moderate to high (0.411, 0.512, and 0.614, respectively). The genetic correlation estimates between BHBAb and MY and FP were low and equal to 0.128 and 0.035, respectively. The genetic correlation estimates between ACETm-MY and ACETm-PP were negative and moderate (-0.374 and - 0.398, respectively). Estimates of genetic correlation between ACETm and FP, FPR, and KET were moderate to high (0.455, 0.626, and 0.876, respectively). The results of this meta-analysis indicated the existence of additive genetic variation for ketosis indicator metabolites which could be exploited in genetic selection programs to reduce ketosis in dairy cows. Moreover, the results propose that selection for lower concentrations of indicator traits could be an effective plan for indirect improvement of production and reproduction performance, and health in dairy cows.

Follicular Dynamics during Estrous Cycle of Pubertal, Mature and Postpartum Crossbred (Nili Ravi × Jianghan) Buffaloes

The follicular dynamics is used as a reliable indicator for reproductive management in livestock. However, the follicular dynamics (follicle wave emergence, estrus cycle length, diameter of dominant follicle, follicular growth and atretic phases) during the estrous cycle of crossbred (Nili Ravi-Jianghan) buffalo is still unexplored. Therefore, the present study aimed to observe the follicular dynamics in estrous cycle of crossbred buffaloes at different physiological stages (pubertal; n = 28, sexual mature; n = 22 and postpartum; n = 18). In the present study, the follicular dynamics were ultrasonically examined at 12 h intervals throughout an estrous cycle during the breeding season. The results indicate that about 86.76% (59/68) crossbred buffaloes, irrespective of physiological stage, exhibited two follicular waves in estrous cycle with an average estrus cycle length was 20.7 ± 0.4 days. The estrus cycle length was significantly shorter (p < 0.05) in pubertal buffaloes (19.4 ± 0.4 days) compared with sexual mature (21.5 ± 0.3 days) and postpartum (21.9 ± 0.4 days) buffaloes. The first follicular wave emerged on same day during one- (pubertal vs. postpartum), two- (pubertal vs. mature vs. postpartum) or three-wave (mature vs. postpartum) estrous cycle buffaloes. The maximum diameter of dominant follicle (DF) in pubertal, sexually mature and postpartum crossbred buffaloes was 9.6 ± 2.0 mm, 10.6 ± 0.5 mm and 12.6 ± 0.7 mm with growth rate of 1.08 ± 0.04 mm/day, 0.92 ± 0.04 mm/day, and 0.9 ± 0.07 mm/day, respectively. In conclusion, similar to other buffalo breeds, Nili Ravi-Jianghan crossbred buffaloes showed the two-wave follicular pattern dominantly with an average duration of ~20 days estrous cycle. The observed follicular dynamics can be used as a reliable indicator for synchronization and fixed-time artificial insemination (FTAI) programs to improve the fertility of crossbred buffaloes.

Genetic Analysis of Persistency for Milk Fat Yield in Iranian Buffaloes (Bubalus bubalis)

This study aimed to estimate heritabilities and genetic trends for different persistency measures for milk fat yield and their genetic correlations with 270-day milk yield in Iranian buffaloes. The records of test-day milk fat yield belonging to the first three lactations of buffaloes within 523 herds consisting of 43,818 records were got from the Animal Breeding Center and Promotion of Animal Products of Iran from 1996 to 2012. To fit the lactation curves based on a random regression test-day model, different orders of Legendre polynomial (LP) functions were selected. Three persistency measures were altered according to the specific condition of the lactation curve in buffaloes: (1) The average of estimated breeding values (EBVs) for test day fat yield from day 226 to day 270 as a deviation from the average of EBVs from day 44 to day 62 (PM₁), (2) A summation of contribution for each day from day 53 to day 247 as a deviation from day 248 (PM₂), and (3) The difference between EBVs for day 257 and day 80 (PM₃). The estimates of heritability for PM₁, PM₂, and PM₃ ranged from 0.20 to 0.48, from 0.36 to 0.47, and from 0.19 to 0.35 over the first three lactations, respectively. The estimate of genetic trends for different persistency measures of milk fat yield was not significant over the lactations (P > 0.05). Genetic correlation estimates between various measures of persistency were generally high over the first three lactations. Also, genetic correlations estimates between persistency measures and 270-day milk yield were mostly low and varied from 0.00 to 0.24 (between PM₁ and 270-day milk yield), from −0.19 to 0.13 (between PM₂ and 270-day milk yield), and from −0.02 to 0.00 (between PM₁ and 270-day milk yield) over the first three lactations, respectively. Persistency measures that showed low genetic correlations with milk fat yield were considered the most suitable measures in selection schemes. Besides, medium to high heritability estimates for different persistency measures for milk fat yield indicated that relevant genetic variations detected for these characters could be regarded in outlining later genetic improvement programs of Iranian buffaloes.

Milk somatic cell count and its relationship with milk yield and quality traits in Italian water buffaloes

In Southern Italy, buffalo (Bubalus bubalis) milk is mostly intended for the manufacture of Mozzarella di Bufala Campana Protected Denomination of Origin (PDO) cheese. Despite the economic boost of the last 2 decades, the buffalo farming system should be improved to maximize the efficiency of the dairy industry, improve yield and quality of milk and cheese, and work toward better animal welfare. Milk somatic cell count (SCC) is used worldwide as an indicator of udder health in individual milk and is useful for monitoring farm hygiene in bulk milk. Mastitis data are currently not available on a large scale in Italy; thus, SCC is essential for identifying animals with suspected udder infection and inflammation. Moreover, high milk SCC is associated with altered composition and acidity, and poor technological properties of milk. However, payment systems of the PDO area are based simply on the delivered volume of milk rather than on quality characteristics. Hence, currently there are no penalties for elevated SCC in bulk milk in the Italian buffalo dairy industry. In addition, SCC for buffalo milk is not mentioned by either the European Community regulations or the PDO protocol, evidencing a lack of rules for the maximum SCC limit. To provide a phenotypic characterization of SCC at the population level and to improve knowledge on buffalo milk quality, 876,299 test-day records of 70,156 buffaloes reared in the PDO area were analyzed. Data revealed that around 11% of herd-test-dates (≥5 animals sampled each) showed average milk SCC ≥400,000 cells/mL (i.e., above the threshold fixed by the European Community for bovine milk). This suggests that there is room to improve SCC at both the farm and individual level. Within first parity, more than 28 and 15% of lactations had average SCC ≥200,000 and ≥300,000 cells/mL, respectively. Both percentages increased with parity and were 39 and 25% in sixth parity, respectively. Supporting this, the proportion of lactations with average SCC ≥500,000 cells/mL increased from 6% in first parity to 12% in sixth parity. Milk yield and SCC were negatively correlated with each other, especially when SCC level was high. An ANOVA was carried out on test-day record milk yield and composition traits, with fixed effects of parity, lactation stage, class of somatic cell score (n = 6), month of calving, and their interactions; buffalo, herd-test-date, and residual were considered random effects. Significantly lower milk yield and lactose percentage were estimated in progressively higher classes of somatic cell score, whereas no significant differences were observed for fat and protein percentages. This is the first attempt to investigate milk SCC in a large data set of Italian dairy buffaloes. These findings may be helpful for defining reliable and effective SCC thresholds to be adopted whenever specific penalties for high SCC are included in milk payment systems. Finally, these results could be used in mastitis monitoring plans aiming to reduce SCC and udder issues at both the individual and farm levels in the Italian buffalo population.

A review on breeding and genetic strategies in Iranian buffaloes (Bubalus bubalis)

The aim of current study was to review breeding progress and update information on genetic strategies in Iranian buffaloes. Iranian buffalo is one of the vital domestic animals throughout north, north-west, south and south-west of Iran with measurable characteristics both in milk and meat production. The species plays an important role in rural economy of the country due to its unique characteristics such as resistance to diseases and parasites, having long productive lifespan and showing higher capability of consuming low-quality forage. In Iran, total production of milk and meat devoted to buffaloes are 293,000 and 24,700 tons, respectively. Selection activities and milk yield recording are carrying out by the central government through the Animal Breeding Centre of Iran. The main breeding activities of Iranian buffaloes included the estimation of genetic parameters and genetic trends for performance traits using different models and methods, estimation of economic values and selection criteria and analysis of population structure. Incorporating different aspects of dairy buffalo management together with improved housing, nutrition, breeding and milking, is known to produce significant improvements in buffalo production. Therefore, identifying genetic potential of Iranian buffaloes, selection of superior breeds, improving nutritional management and reproduction and developing the education and increasing the skills of practical breeders can be useful in order to enhance the performance and profitability of Iranian buffaloes.