Metabolomics Approach Explore Diagnostic Biomarkers and Metabolic Changes in Heat-Stressed Dairy Cows

Heat Stress Effects on Physiological and Milk Yield Traits of Lactating Holstein Friesian Crossbreds Reared in Tanga Region, Tanzania

Global warming caused by climate change is a challenge for dairy farming, especially in sub-Saharan countries. Under high temperatures and relative humidity, lactating dairy cows suffer from heat stress. The objective of this study was to investigate the effects and relationship of heat stress (HS) measured by the temperature-humidity index (THI) regarding the physiological parameters and milk yield and composition of lactating Holstein Friesian crossbred dairy cows reared in the humid coastal region of Tanzania. A total of 29 lactating Holstein Friesian x Zebu crossbred dairy cows with 50% (HF50) and 75% (HF75) Holstein Friesian gene levels in the second and third months of lactation were used. The breed composition of Holstein Friesians was determined based on the animal recording system used at the Tanzania Livestock Research Institute (TALIRI), Tanga. The data collected included the daily temperature, relative humidity, daily milk yield, and physiological parameters (core body temperature, rectal temperature, respiratory rate, and panting score). THI was calculated using the equation of the National Research Council. The THI values were categorized into three classes, i.e., low THI (76-78), moderate THI (79-81), and high THI (82-84). The effects of THI on the physiological parameters and milk yield and composition were assessed. The effects of the genotype, the parity, the lactation month, and the interaction of these parameters with THI on the milk yield, milk composition, and physiological parameters were also investigated. The results show that THI and its interaction with genotypes, parity, and the lactation month had a highly significant effect on all parameters. THI influenced (p ˂ 0.05) the average daily milk yield and milk fat %, protein %, lactose %, and solids-not-fat %. As the THI increased from moderate to high levels, the average daily milk yield declined from 3.49 ± 0.04 to 3.43 ± 0.05 L/day, while the fat % increased from 2.66 ± 0.05% to 3.04 ± 0.06% and the protein decreased from 3.15 ± 0.02% to 3.13 ± 0.03%. No decline in lactose % was observed, while the solid-not-fat % declined from 8.56 ± 0.08% to 8.55 ± 0.10% as the THI values increased from moderate to high. Also, the THI influenced physiological parameters (p ˂ 0.05). The core body temperature (CBT), rectal temperature (RT), respiratory rate (RR) and panting score (PS) increased from 35.60 ± 0.01 to 36.00 ± 0.01 °C, 38.03 ± 0.02 to 38.30 ± 0.02 °C, 62.53 ± 0.29 to 72.35 ± 0.28 breaths/min, and 1.35 ± 0.01 to 1.47 ± 0.09, respectively, as the THI increased from low to high. The THI showed a weak positive correlation with the average daily milk yield and fat percentage, whereas the protein, lactose, and solids-not-fat percentages showed negative relationships with THI (p ≤ 0.05). CBT, RT, RR, and PS showed positive relationships (p ≤ 0.05) with THI. These negative relationships indicate that there is an antagonistic correlation between sensitivity to HS and the level of production. It is concluded that the THI, the genotype, the parity, and the lactation month, along with their interactions with THI, significantly influenced the milk yield, milk composition, and physiological parameters of lactating Holstein Friesian dairy crosses at THI thresholds ranging from 77 to 84.

Effects of heat stress on feed intake, milk yield, milk composition, and feed efficiency in dairy cows: A meta-analysis

Heat stress compromises dairy production by decreasing feed intake and milk yield, and it may also alter milk composition and feed efficiency. However, little information is available for evaluating such effects across different levels of heat stress and cows enrolled in heat stress studies. The objectives of this study were to evaluate the effects of heat stress on dry matter intake (DMI), energy-corrected milk (ECM), milk composition, and feed efficiency (kg ECM/kg DMI) and to investigate the relationship between such effects and heat stress intervention and animal characteristics by using meta-analytical approaches. Data from 31 studies (34 trials) fulfilled the inclusion criteria and were used for analysis. Results showed that heat stress decreased DMI, ECM, and milk protein concentration, but did not alter milk fat concentration or feed efficiency. Meta-regression confirmed that such reductions in DMI and ECM were significantly associated with increasing temperature-humidity index (THI). Over the period of heat stress, for each unit increase in THI, DMI and ECM decreased by 4.13% and 3.25%, respectively, in mid-lactation cows. Regression models further revealed the existence of a strong interaction between THI and lactation stage, which partially explained the large heterogeneity in effect sizes of DMI and ECM. The results indicated a need for more research on the relationship between the effect of heat stress and animal characteristics. This study calls for the implementation of mitigation strategies in heat-stressed herds due to the substantial decrease in productivity.

Plasma and milk metabolomics revealed changes in amino acid metabolism in Holstein dairy cows under heat stress

Our understanding of metabolic alterations triggered by heat stress is incomplete, which limits the designing of nutritional strategies to mitigate negative productive and health effects. Thus, this study aimed to explore the metabolic responses of heat-stressed dairy cows to dietary supplementation with vitamin D3/Ca and vitamin E/Se. Twelve multiparous Holstein cows were enrolled in a split-plot Latin square design with two distinct vitamin E/Se supplementation levels, either at a low (ESe-, n = 6, 11.1 IU/kg vitamin E and 0.55 mg/kg Se) or a high dose (ESe+, n = 6 223 IU/kg vitamin E and 1.8 mg/kg Se) as the main plot. Treatment subplots, arranged in a replicated 3 × 3 Latin square design, comprised heat challenge (Temperature Humidity Index, THI: 72.0-82.0) supplemented with different levels of vitamin D3/Ca: either low (HS/DCa-, 1 012 IU/kg and 0.73%, respectively) or high (HS/DCa+, 3 764 IU/kg and 0.97%, respectively), and a pair-fed control group in thermoneutrality (THI = 61.0-64.0) receiving the low dose of vitamin D3/Ca (TN). The liquid chromatography-mass spectrometry-based metabolome profile was determined in blood plasma and milk sampled at the beginning (day 0) and end (day 14) of each experimental period. The results were analyzed for the effect of (1) TN vs. HS/ESe-/DCa-, and (2) the vitamin E/Se and vitamin D3/Ca supplementation. No group or group × day effects were detected in the plasma metabolome (false discovery rate, FDR > 0.05), except for triglyceride 52:2 being higher (FDR = 0.03) on day 0 than 14. Taurine, creatinine and butyryl-carnitine showed group × day interactions in the milk metabolome (FDR ≤ 0.05) as creatinine (+22%) and butyryl-carnitine (+190%) were increased (P < 0.01) on day 14, and taurine was decreased (-65%, P < 0.01) on day 14 in the heat stress (HS) cows, compared with day 0. Most compounds were unaffected by vitamin E/Se or vitamin D3/Ca supplementation level or their interaction (FDR > 0.05) in plasma and milk, except for milk alanine which was lower (-69%, FDR = 0.03) in the E/Se+ groups, compared with E/Se-. Our results indicated that HS triggered more prominent changes in the milk than in the plasma metabolome, with consistent results in milk suggesting increased muscle catabolism, as reflected by increased creatinine, alanine and citrulline levels. Supplementing with high levels of vitamin E/Se or vitamin D3/Ca or their combination did not appear to affect the metabolic remodeling triggered by HS.

Differential Responses of Physiological Parameters, Production Traits, and Blood Metabolic Profiling between First- and Second-Parity Holstein Cows in the Comparison of Spring versus Summer Seasons

Heat stress (HS) negatively influences cows' welfare and productivity. Therefore, a better understanding of the physiological and molecular mechanisms of HS responses from multiple parities is paramount for the development of effective management and breeding strategies. In comparison with first-parity cows in the spring (Spring-1), first-parity cows in the summer (Summer-1) had a significantly higher rectal temperature (RT), respiration rate (RR), drooling score (DS), and daily activity (DA), while lower (P < 0.05) daily rumination (DR), seven-day average milk yield (7AMY), milk yield on sampling day (MY_S), milk yield on test day (MY_T), and lactose percentage (LP) were observed. When comparing the spring (Spring-2) and summer (Summer-2) of the second-parity cows, significant differences were also found in RT, RR, DS, DA, and DR (P < 0.05), corresponding to similar trends with the first parity while having smaller changes. Moreover, significantly negative impacts on performance traits were only observed on fat percentage (FP) and LP. These results showed that there were different biological responses between first- and second-parity Holstein cows. Further, 18 and 17 metabolites were involved in the seasonal response of first- and second-parity cows, respectively. Nine differential metabolites were shared between the two parities, and pathway analyses suggested that cows had an inhibited tricarboxylic acid cycle, increased utilization of lipolysis, and a dysregulated gut microbiome during the summer. The metabolites identified exclusively for each parity highlighted the differences in microbial response and host amino acid metabolism between two parities in response to HS. Moreover, glucose, ethanol, and citrate were identified as potential biomarkers for distinguishing individuals between Spring-1 and Summer-1. Ethanol and acetone were better predictors for distinguishing individuals between Spring-2 and Summer-2. Taken together, the present study demonstrated the impact of naturally induced HS on physiological parameters, production traits, and the blood metabolome of Holstein cows. There are different biological responses and regulation mechanisms between first- and second-parity Holstein cows.

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.

Genetic Markers Associated with Milk Production and Thermotolerance in Holstein Dairy Cows Managed in a Heat-Stressed Environment

Dairy production in Holstein cows in a semiarid environment is challenging due to heat stress. Under such conditions, genetic selection for heat tolerance appears to be a useful strategy. The objective was to validate molecular markers associated with milk production and thermotolerance traits in Holstein cows managed in a hot and humid environment. Lactating cows (n = 300) exposed to a heat stress environment were genotyped using a medium-density array including 53,218 SNPs. A genome-wide association study (GWAS) detected six SNPs associated with total milk yield (MY305) that surpassed multiple testing (p < 1.14 × 10^-6). These SNPs were further validated in 216 Holstein cows from two independent populations that were genotyped using the TaqMan bi-allelic discrimination method and qPCR. In these cows, only the SNPs rs8193046, rs43410971, and rs382039214, within the genes TLR4, GRM8, and SMAD3, respectively, were associated (p < 0.05) with MY305, rectal temperature (RT), and respiratory rate. Interestingly, these variables improved as the number of favorable genotypes of the SNPs increased from 0 to 3. In addition, a regression analysis detected RT as a significant predictor (R² = 0.362) for MY305 in cows with >1 favorable genotype, suggesting this close relationship was influenced by genetic markers. In conclusion, SNPs in the genes TLR4, GRM8, and SMAD3 appear to be involved in the molecular mechanism that regulates milk production in cows under heat-stressed conditions. These SNPs are proposed as thermotolerance genetic markers for a selection program to improve the milk performance of lactating Holstein cows managed in a semiarid environment.

Metabolomic and transcriptomic study to understand changes in metabolic and immune responses in steers under heat stress

Heat stress (HS) damages livestock by adversely affecting physiological and immunological functions. However, fundamental understanding of the metabolic and immunological mechanisms in animals under HS remains elusive, particularly in steers. To understand the changes on metabolic and immune responses in steers under HS condition, we performed RNA-sequencing and proton nuclear magnetic resonance spectroscopy-based metabolomics on HS-free (THI value: 64.92 ± 0.56) and HS-exposed (THI value: 79.13 ± 0.56) Jersey steer (n = 8, body weight: 559.67 ± 32.72 kg). This study clarifies the metabolic changes in 3 biofluids (rumen fluid, serum, and urine) and the immune responses observed in the peripheral blood mononuclear cells of HS-exposed steers. This integrated approach allowed the discovery of HS-sensitive metabolic and immunological pathways. The metabolomic analysis indicated that HS-exposed steers showed potential HS biomarkers such as isocitrate, formate, creatine, and riboflavin (P < 0.05). Among them, there were several integrative metabolic pathways between rumen fluid and serum. Furthermore, HS altered mRNA expression and immune-related signaling pathways. A meta-analysis revealed that HS decreased riboflavin metabolism and the expression of glyoxylate and dicarboxylate metabolism-related genes. Moreover, metabolic pathways, such as the hypoxia-inducible factor-1 signaling pathway, were downregulated in immune cells by HS (P < 0.05). These findings, along with the datasets of pathways and phenotypic differences as potential biomarkers in steers, can support more in-depth research to elucidate the inter-related metabolic and immunological pathways. This would help suggest new strategies to ameliorate the effects of HS, including disease susceptibility and metabolic disorders, in Jersey steers.

Exploration of urinary metabolite dynamicity for early detection of pregnancy in water buffaloes

Early and precise pregnancy diagnosis can reduce the calving interval by minimizing postpartum period. The present study explored the differential urinary metabolites between pregnant and non-pregnant Murrah buffaloes (Bubalus bubalis) during early gestation to identify potential pregnancy detection biomarkers. Urine samples were collected on day 0, 10, 18, 35 and 42 of gestation from the pregnant (n = 6) and on day 0, 10 and 18 post-insemination from the non-pregnant (n = 6) animals. ¹H-NMR-based untargeted metabolomics followed by multivariate analysis initially identified twenty-four differentially expressed metabolites, among them 3-Hydroxykynurenine, Anthranilate, Tyrosine and 5-Hydroxytryptophan depicted consistent trends and matched the selection criteria of potential biomarkers. Predictive ability of these individual biomarkers through ROC curve analyses yielded AUC values of 0.6–0.8. Subsequently, a logistic regression model was constructed using the most suitable metabolite combination to improve diagnostic accuracy. The combination of Anthranilate, 3-Hydroxykynurenine, and Tyrosine yielded the best AUC value of 0.804. Aromatic amino acid biosynthesis, Tryptophan metabolism, Phenylalanine and Tyrosine metabolism were identified as potential pathway modulations during early gestation. The identified biomarkers were either precursors or products of these metabolic pathways, thus justifying their relevance. The study facilitates precise non-invassive urinary metabolite-based pen-side early pregnancy diagnostics in buffaloes, eminently before 21 days post-insemination.

Improving Genomic Selection for Heat Tolerance in Dairy Cattle: Current Opportunities and Future Directions

Heat tolerance is the ability of an animal to maintain production and reproduction levels under hot and humid conditions and is now a trait of economic relevance in dairy systems worldwide because of an escalating warming climate. The Australian dairy population is one of the excellent study models for enhancing our understanding of the biology of heat tolerance because they are predominantly kept outdoors on pastures where they experience direct effects of weather elements (e.g., solar radiation). In this article, we focus on evidence from recent studies in Australia that leveraged large a dataset [∼40,000 animals with phenotypes and 15 million whole-genome sequence variants] to elucidate the genetic basis of thermal stress as a critical part of the strategy to breed cattle adapted to warmer environments. Genotype-by-environment interaction (i.e., G × E) due to temperature and humidity variation is increasing, meaning animals are becoming less adapted (i.e., more sensitive) to changing environments. There are opportunities to reverse this trend and accelerate adaptation to warming climate by 1) selecting robust or heat-resilient animals and 2) including resilience indicators in breeding goals. Candidate causal variants related to the nervous system and metabolic functions are relevant for heat tolerance and, therefore, key for improving this trait. This could include adding these variants in the custom SNP panels used for routine genomic evaluations or as the basis to design specific agonist or antagonist compounds for lowering core body temperature under heat stress conditions. Indeed, it was encouraging to see that adding prioritized functionally relevant variants into the 50k SNP panel (i.e., the industry panel used for genomic evaluation in Australia) increased the prediction accuracy of heat tolerance by up to 10% units. This gain in accuracy is critical because genetic improvement has a linear relationship with prediction accuracy. Overall, while this article used data mainly from Australia, this could benefit other countries that aim to develop breeding values for heat tolerance, considering that the warming climate is becoming a topical issue worldwide.

Heat Stress Induces Shifts in the Rumen Bacteria and Metabolome of Buffalo

Exposure to the stress (HS) negatively affects physiology, performance, reproduction and welfare of buffalo. However, the mechanisms by which HS negatively affects rumen bacteria and its associated metabolism in buffalo are not well known yet. This study aimed to gain insight into the adaption of bacteria and the complexity of the metabolome in the rumen of six buffalo during HS using 16S rDNA and gas chromatography metabolomics analyses. HS increased respiratory rate (p < 0.05) and skin temperature (p < 0.01), and it decreased the content of acetic acid (p < 0.05) and butyric acid (p < 0.05) in the rumen. Omics sequencing revealed that the relative abundances of Lachnospirales, Lachnospiraceae, Lachnospiraceae_NK3A20_group and Clostridia_UCG-014 were significantly (p < 0.01) higher under HS than non-heat stress conditions. Several bacteria at different levels, such as Lactobacillales, Streptococcus, Leuconostocaceae and Leissella, were significantly (p < 0.05) more abundant in the rumen of the non-heat stress than HS condition. Thirty-two significantly different metabolites closely related to HS were identified (p < 0.05). Metabolic pathway analysis revealed four key pathways: D-Alanine metabolism; Lysine degradation, Tropane; piperidine and pyridine alkaloid biosynthesis; and Galactose metabolism. In summary, HS may negatively affected rumen fermentation efficiency and changed the composition of rumen community and metabolic function.

Metabolomics Analysis Across Multiple Biofluids Reveals the Metabolic Responses of Lactating Holstein Dairy Cows to Fermented Soybean Meal Replacement

This experiment was performed to reveal the metabolic responses of dairy cows to the replacement of soybean meal (SBM) with fermented soybean meal (FSBM). Twenty-four lactating Chinese Holstein dairy cattle were assigned to either the SBM group [the basal total mixed ration (TMR) diet containing 5.77% SBM] or the FSBM group (the experimental TMR diet containing 5.55% FSBM), in a completely randomized design. The entire period of this trial consisted of 14 days for the adjustment and 40 days for data and sample collection, and sampling for rumen liquid, blood, milk, and urine was conducted on the 34th and 54th day, respectively. When SBM was completely replaced by FSBM, the levels of several medium-chain FA in milk (i.e., C13:0, C14:1, and C16:0) rose significantly (p < 0.05), while the concentrations of a few milk long-chain FA (i.e., C17:0, C18:0, C18:1n9c, and C20:0) declined significantly (p < 0.05). Besides, the densities of urea nitrogen and lactic acid were significantly (p < 0.05) higher, while the glucose concentration was significantly (p < 0.05) lower in the blood of the FSBM-fed cows than in the SBM-fed cows. Based on the metabolomics analysis simultaneously targeting the rumen liquid, plasma, milk, and urine, it was noticed that substituting FSBM for SBM altered the metabolic profiles of all the four biofluids. According to the identified significantly different metabolites, 3 and 2 amino acid-relevant metabolic pathways were identified as the significantly different pathways between the two treatments in the rumen fluid and urine, respectively. Furthermore, glycine, serine, and threonine metabolism, valine, leucine, and isoleucine biosynthesis, and cysteine and methionine metabolism were the three key integrated different pathways identified in this study. Results mainly implied that the FSBM replacement could enhance nitrogen utilization and possibly influence the inflammatory reactions and antioxidative functions of dairy cattle. The differential metabolites and relevant pathways discovered in this experiment could serve as biomarkers for the alterations in protein feed and nitrogen utilization efficiency of dairy cows, and further investigations are needed to elucidate the definite roles and correlations of the differential metabolites and pathways.

Plasma metabolomic profiles as affected by diet and stress in Spanish goats

The effects of high-condensed tannin (CT) diet combined with preslaughter stress have not been studied at the metabolome level in goats. This study was conducted to determine the effects of feeding sericea lespedeza (SL; Lespedeza cuneata), a high-CT legume, and transportation stress on plasma metabolome in goats. Uncastrated male Spanish goats (age = 8 months; BW = 26.0 ± 0.48 kg) were either fed ground 'Serala' SL hay (SER), bermudagrass (Cynodon dactylon) hay (BG), or bermudagrass hay-dewormed goats (BG-DW; Control) at 75% of intake, with a corn-based supplementation (25%) for 8 weeks (n = 12/Diet). At the end of the trial, goats were subjected to one of two stress treatments (ST): transported for 90 min to impose stress (TS) or held in pens (NTS) before slaughtering, in two replicates. Live and carcass weights, and blood samples were collected at 0, 30, 60 and 90 min of transportation or holding time (Time). The data were analyzed using MIXED procedures in SAS and metabolomics data were analyzed using the R software. When measured after ST, SER group had the lowest body weight (P < 0.05) among the three diet groups. Carcass weights were high in the BG-DW, low in SER, and intermediate in BG group. Plasma creatine concentrations decreased over Time (P < 0.01) in the TS goats in all diet groups. Meat crude protein percentages were higher (P < 0.05) in SER (22.5 ± 0.22) and BG-DW (22.3 ± 0.22) groups compared to the BG group (21.6 ± 0.22). At the metabolome level, SER group had the lowest (P < 0.05) glycine, alanine, threonine, taurine, trans-hydroxyproline, methionine, and histidine concentrations and highest (P < 0.01) lysine and citrulline concentrations among the Diet groups. Butyric acid, concentration was higher (P < 0.05) in the SER group compared to BG group. Eight medium- and long-chained acylcarnitines were higher (P < 0.05) in the BG-DW group than SER or BG groups. In general, amino acid levels decreased and acylcarnitine increased with Time (P < 0.05) in all groups. Sericea diet can be beneficial in enhancing stress coping abilities in goats due to elevated butyrate, lysine, and citrulline levels; however, SER resulted in lower energy level in goats compared to BG or BG-DW groups. Fatty acid metabolism is the main energy pathway in all groups during prolonged stress. Inclusion of certain varieties of SL in the diet must be carefully controlled to prevent possible negative effect.

LC-MS/MS Based Metabolomics Reveal Candidate Biomarkers and Metabolic Changes in Different Buffalo Species

Consumers have shown more and more interest in high-quality and healthy dairy products and buffalo milk is commercially more viable than other milks in producing superior dairy products due to its higher contents of fat, crude protein, and total solids. Metabolomics is one of the most powerful strategies in molecular mechanism research however, little study has been focused on the milk metabolites in different buffalo species. Therefore, the aim of this study was to explore the underlying molecular mechanism of the fatty synthesis and candidate biomarkers by analyzing the metabolomic profiles. Milk of three groups of buffaloes, including 10 Mediterranean, 12 Murrah, and 10 crossbred buffaloes (Murrah × local swamp buffalo), were collected and UPLC-Q-Orbitrap HRMS was used to obtain the metabolomic profiles. Results showed that milk fatty acid in Mediterranean buffalo was significantly higher than Murrah buffalo and crossbred buffalo. A total of 1837/726 metabolites was identified in both positive and negative electrospray ionization (ESI±) mode, including 19 significantly different metabolites between Mediterranean and Murrah buffalo, and 18 different metabolites between Mediterranean and crossbred buffalo. We found 11 of the different metabolites were both significantly different between Mediterranean vs. Murrah group and Mediterranean vs crossbred group, indicating that they can be used as candidate biomarkers of Mediterranean buffalo milk. Further analysis found that the different metabolites were mainly enriched in fat synthesis related pathways such as fatty acid biosynthesis, unsaturated fatty acid biosynthesis, and linoleic acid metabolism, indicating that the priority of different pathways affected the milk fat content in different buffalo species. These specific metabolites may be used as biomarkers in the identification of milk quality and molecular breeding of high milk fat buffalo.