The age at first calving (AFC) is an important trait to be considered in breeding programmes of dairy buffaloes, where new approaches and technologies, such as genomic selection, are constantly applied. Thus, the objective of this study was to compare the predictive ability of different genomic single-step methods using AFC information from Murrah buffaloes. From a pedigree file containing 3320 buffaloes, 2247 cows had AFC records and 553 animals were genotyped. The following models were performed: pedigree-based BLUP (PBLUP), single-step GBLUP (ssGBLUP), weighted single-step GBLUP (WssGBLUP), and single-step Bayesian regression methods (ssBR-BayesA, BayesBπ, BayesCπ, Bayes-Lasso, and BayesRR). To compare the methodologies, the accuracy and dispersion of (G)EBVs were assessed using the LR method. Accuracy estimates for the genotyped animals ranged from 0.30 (PBLUP) to 0.39 (WssGBLUP). Predictions with the traditional model (PBLUP) were very dispersed from what was expected, while BayesCπ (0.99) and WssGBLUP (1.00) obtained the lowest dispersion. The results indicate that the use of genomic information can improve the genetic gain for AFC by increasing the accuracy and reducing inflation/deflation of predictions compared to the traditional pedigree-based model. In addition, among all genomic single-step models studied, WssGBLUP and single-step BayesA were the most advantageous methods to be used in the genomic evaluation of AFC of buffaloes from this population.

The objective of the present study was to evaluate the relationship between body weight (BW) and hip width (HW) in dairy buffaloes (Bubalus bubalis). HW was measured in 215 Murrah buffaloes with a BW of 341 ± 161.6 kg, aged between three months and five years, and raised in southeastern Mexico. Linear and non-linear regressions were used to construct the prediction models. The goodness of fit of the models was evaluated using the Akaike information criterion (AIC), Bayesian information criterion (BIC), coefficient of determination (R2), mean squared error (MSE), and root MSE (RMSE). Additionally, the developed models were evaluated through internal and external cross-validation (k-folds) using independent data. The ability of the fitted models to predict the observed values was assessed based on the root mean square error of prediction (RMSEP), R2, and mean absolute error (MAE). The relationship between BW and HW showed a high correlation coefficient (r = 0.96, P < 0.001). The chosen fitted model to predict BW was: −176.33 (± 40.83***) + 8.74 (± 1.79***) × HW + 0.04 (± 0.01*) × HW2, because it presented the lowest MSE, RMSE, and AIC values, which were 1228.64, 35.05 and 1532.41, respectively. Therefore, with reasonable accuracy, the quadratic model using hip width may be suitable for predicting body weight in buffaloes.

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 dm3) 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 (R2), 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, R2, 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.

This study evaluated the effect of replacing soybean meal by a blend of ground corn and urea-ammonium sulphate (GCU-S) in the diet of lactating buffaloes on milk production and composition, digestibility, N balance and blood metabolites. Twelve multiparous dairy Murrah buffaloes (Bubalus bubalis), at 100 ± 4 d in milk and yielding 10 ± 2.5 kg/d, were randomly distributed in a triple 4 × 4 Latin square, with four different inclusions of GCU-S at U-S levels: 0 (control), 8.0, 16.4 and 24.1 g/kg dry matter (DM) total. Replacing soybean meal with GC-US had a significant depressing effect on absolute DM intake, which was still numerically evident but no longer significant when expressed on a body weight-related basis. Intakes of crude protein (CP) and N as well as N-urinary excretion were also significantly depressed. Digestibility of dry matter was improved and, as a proportion of DMI, intake of total digestible nutrients increased significantly. Perhaps as a result, feed efficiency (kg DMI required per kg 6% fat-corrected milk) was significantly improved. Efficiency of protein use for milk production was significantly improved, but not when expressed as milk protein output. N transfer into milk, as a proportion of total N intake, also increased significantly. However, GC-US inclusion had no significant effect on milk production or milk composition, nor did it affect serum metabolites, digestibilities (apart from dry matter) or N measures of balance apart from those mentioned above. Inclusion of this blend of ground corn with urea and ammonium suplate can be recommended for dairy buffalo because, although it decreases DMI and N-excretion, it improves DM digestibility and feeding efficiency whilst maintaining milk production and composition.

Buffalo milk production has become of significant importance on the world scale, however, there are few studies involving biotechnological tools specifically for buffalo. To verify the effects caused by subclinical mastitis on the components of milk and to study the innate immune system in the udder of dairy buffaloes with subclinical mastitis, we evaluated the levels of expression of the lactoferrin (LTF), tumor necrosis factor alpha (TNF-α), interleukin-1 beta (IL-1β), interleukin-8 (IL-8), and toll-like receptors 2 (TLR-2) and 4 (TLR-4) genes in buffaloes with and without subclinical mastitis. Milk samples were collected for the determination of milk components: somatic cell score (SCS), fat, protein, lactose, total solids and solids-not-fat (SNF), as well as for RNA extraction of milk cells, complementary DNA synthesis, and expression profile quantification by quantitative real-time PCR. For gene expression, the ΔΔCt was estimated using contrasts of the target genes expression adjusted for the expression of the housekeeping genes between both groups. Linear regression analysis was performed to determine the relationship between the genes studied and the milk components. Subclinical mastitis induced changes in the fat, lactose and SNF in milk of buffaloes, and the messenger RNA abundance was upregulated for TLR-2, TLR-4, TNF-α, IL-1β and IL-8 genes in milk cells of buffaloes with subclinical mastitis, whereas the LTF gene was not differentially expressed. Results of linear regression analysis showed that TLR-2 gene expression most explains the variation in SCS, and the change in a unit of ΔCt of the TNF-α gene would result in a higher increase in SCS. The study of these immune function genes that are active in the mammary gland is important to characterize the action mechanism of the innate immunity that occurs in subclinical mastitis in dairy buffaloes and may aid the development of strategies to preserve the health of the udder.

Water buffalo (Bubalus bubalis) is important for the economy of Asia, South America and parts of Europe. Coccidiosis is an important cause of neonatal mortality in livestock, including buffalo. Of more than 12 species of Eimeria in buffalo, Eimeria bareillyi is the most pathogenic. There are uncertainties concerning its asexual and sexual development. During a previously reported outbreak of fatal enteritis associated with E. bareillyi in buffaloes in the Netherlands, sections of small intestine were re-evaluated histologically and by transmission electron microscopy (TEM) to seek details of endogenous development. Profuse asexual multiplication occurred in the jejunum and ileum. Light microscopic examination revealed that parasites divided in two (probably endodyogeny) or more organisms. There were two or more generations of morphologically different merozoites; some of these observations were confirmed by TEM. Details of gametogonic development, including oocyst wall formation are provided. Schizogonic and gametogonic development described in the present study can serve as a guide for differential diagnosis of Eimeria species in histological sections of intestines of buffaloes.

The buffalo has a seasonal reproduction activity with mating and non-mating periods occurring from late autumn to winter and from late spring to beginning of autumn, respectively. Sperm glycocalyx plays an important role in reproduction as it is the first interface between sperm and environment. Semen quality is poorer during non-mating periods, so we aimed to evaluate if there were also seasonal differences in the surface glycosylation pattern of mating period spermatozoa (MPS) compared with non-mating period spermatozoa (NMPS). The complexity of carbohydrate structures makes their analysis challenging, and recently the high-throughput microarray approach is now providing a new tool into the evaluation of cell glycosylation status. We adopted a novel procedure in which spermatozoa was spotted on microarray slides, incubated with a panel of 12 biotinylated lectins and Cy3-conjugated streptavidin, and then signal intensity was detected using a microarray scanner. Both MPS and NMPS microarrays reacted with all the lectins and revealed that the expression of (i) O-glycans with NeuNAcα2-3Galβ1,3(±NeuNAcα2-6)GalNAc, Galβ1,3GalNAc and GalNAcα1,3(l-Fucα1,2)Galβ1,3/4GlcNAcβ1 was not season dependent; (ii) O-linked glycans terminating with GalNAc, asialo N-linked glycans terminating with Galβ1,4GlcNAc, GlcNAc, as well as α1,6 and α1,2-linked fucosylated oligosaccharides was predominant in MPS; (iii) high mannose- and biantennary complex types N-glycans terminating with α2,6 sialic acids and terminal galactose were lower in MPS. Overall, this innovative cell microarray method was able to identify specific glycosylation changes that occur on buffalo bull sperm surface during the mating and non-mating periods.

The objectives of this study were to analyse buffalo milk fat composition, to verify the activity of Delta(9)-desaturase enzyme in the mammary gland, as well as to estimate additive genetic variances for milk, fat and protein yield, and milk cis-9,trans-11 conjugated linoleic acid percentage (cis-9,trans-11 CLA%). A total of 3929 lactation milk yields (MY) records from 2130 buffaloes and 1598 lactation fat (FY) and protein (PY) yield records from 914 buffaloes were analysed. For cis-9,trans-11 CLA%percentage, a total of 661 milk samples from 225 buffaloes, daughters of 8 sires, belonging to 4 herds and calving in 2003 and 2004, were used. The genetic parameters and variance components were estimated by Restricted Maximum Likelihood applying an animal model. The fixed effects considered in the model were: contemporary group (herd, year, calving season) and age at calving (linear and quadratic effects) and lactation length (linear and quadratic effects) as covariables. Additive genetic and permanent environment effects were considered as random. The MY, FY, PY and CLA% means were 1482±355 kg, 90·1±24·6 kg, 56·9±15·2 kg and 0·69±0·16%, respectively. Heritability estimates for MY, FY, PY and CLA% were 0·28±0·05, 0·26±0·11, 0·25±0·11 and 0·35±0·14, respectively. There is enough additive genetic variation for buffalo milk, protein and fat yield to improve these traits through selection. The cis-9,trans-11 CLA% can be enhanced by selection in buffaloes and will contribute to improving human health. The activity and efficiency of Delta(9)-desaturase in the mammary was measured and confirmed.

In order to contribute to the breeding programmes of Asian water buffalo, the aim of this study was to analyse the influence of genetic effects in the stayability of Murrah dairy buffaloes. The stayability trait (ST) was defined as the female's ability to stay in the herd for one (ST1), two (ST2), three (ST3), four (ST4), five (ST5) or six years (ST6) after the first calving. The same trait was also considered as continuous and was designated stayability in days up to one (STD1), two (STD2), three (STD3), four (STD4), five (STD5) or six years (STD6) after the first calving. Data from 1016 females reared in nine herds located in the State of São Paulo, Brazil, were analysed. Statistical models included the additive genetic effect of the animal and the fixed effects of the buffalo breeding herd, birth year and birth season. Additive effects for ST were estimated by approximate restricted maximum likelihood using a threshold model, while for STD, the additive effects were estimated by restricted maximum likelihood. Heritability estimates were lower for ST, except for ST1, (0·11±0·07, 0·17±0·06, 0·23±0·06, 0·16±0·08, 0·14±0·09 and 0·16±0·10 for ST1, ST2, ST3, ST4, ST5 and ST6, respectively) when compared with STD (0·05±0·06, 0·18±0·08, 0·40±0·10, 0·49±0·11, 0·41±0·11 and 0·30±0·13, for STD1, STD2, STD3, STD4, STD5 and STD6, respectively). Considering the values of heritability and owing to the serial nature of STD to a specific age, selection for STD3 should have a favourable influence on STD to other ages.

In order to contribute to the genetic breeding programs of buffaloes, this study aimed to determine the influence of environmental effects on the stayability (ST) of dairy female Murrah buffalo in the herd. Data from 1016 buffaloes were used. ST was defined as the ability of the female to remain in the herd for 1, 2, 3, 4, 5 or 6 years after the first calving. Environmental effects were studied by survival analysis, adjusted to the fixed effects of farm, year and season of birth, class of first-lactation milk yield and age at first calving. The data were analyzed using the LIFEREG procedure of the SAS program that fits parametric models to failure time data (culling or ST = 0), and estimates parameters by maximum likelihood estimation. Breeding farm, year of birth and first-lactation milk yield significantly influenced (P < 0.0001) the ST to the specific ages (1 to 6 years after the first calving). Buffaloes that were older at first calving presented higher probabilities of being culled 1 year after the first calving, without any effect on culling at older ages. Buffaloes with a higher milk yield at first calving presented a lower culling probability and remained for a longer period of time in the herd. The effects of breeding farm, year of birth and first-lactation milk yield should be included in models used for the analysis of ST in buffaloes.

We censused wild buffalo Bubalus arnee in Koshi Tappu Wildlife Reserve, Nepal, in March 2004 using methods employed in earlier surveys, and estimated a population of 159 animals. Since the last census in 2000 the management situation has deteriorated. Guard posts have been evacuated due to the Maoist insurgency. The entire Army Battalion usually posted in the Reserve was at Headquarters at the time of this study and therefore there were no patrols over most of the Reserve, and much human encroachment. Mortality from flooding and road deaths, and possibly poaching for meat, were evident, and males suffered more mortality overall than females. Despite these threats the population had increased since 2000, albeit at a lower rate than previously. Active management, including interventions within the Reserve and a translocation of some individuals to Chitwan National Park, are recommended.