Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and Bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records. © 2013 American Dairy Science Association.

Genetic parameters for test-day milk yield, 305-day milk yield, and lactation length in Guzerat cows

Milk production in tropical environments requires the use of crossbreeding systems including breeds well adapted to harsh conditions, but with lower productivities when compared to specialized breeds. Besides the genetic improvement for milk production, lactation lengths also need to be studied for most of these breeds. Accordingly, genetic parameters were estimated for 305-day cumulative milk yield (MY305), test-day milk yield (TDMY), and lactation length (LL) using information from the first lactations of 2816 Guzerat cows selected for milk production in 28 herds in Brazil. Contemporary groups were defined as herd, year and season of the test for TDMY, and as herd, year and season of calving for MY305 and LL. Variance components were estimated with the restricted maximum likelihood method under a multi-trait animal model. Heritabilities estimated for TDMY ranged from 0.16 to 0.24, and were 0.24 and 0.12 for MY305 and LL, respectively. Genetic correlations were high and positive, ranging from 0.51 to 0.99 among TDMY records, from 0.81 to 0.98 between each TDMY and MY305, and from 0.71 to 0.94 between each TDMY and LL. Genetic parameters obtained in this study indicated the possibility of using test-day records for the prediction of breeding values for milk yield in this population of the Guzerat breed. The use of TDMY as selection criteria would result in indirect gains in MY305 and LL. However, the highest response to selection for MY305 would be obtained by direct selection for this trait. © 2012 Elsevier B.V.

Effect of Psychrotrophic Growth on the Milk Fat Fraction at Different Temperatures of Storage

The use of cooling, without using adequate hygienic practices in primary milk production, allows for the growth of psychrotrophic microorganisms that produce the thermoresistant lipases that give milk a rancid flavor. This study aimed to verify how the variation in temperature influences the lipolytic metabolism of the psychrotrophic organisms. Samples of raw milk were collected and submitted to laboratorial analysis as follows: psychrotrophic bacteria count, lipolytic bacteria count, and free fatty acids dosage. Each sample was divided into 3 aliquots and then incubated at 4, 8, and 12 °C, respectively. For each temperature, analyses were repeated after 12, 24, and 48 h of storage. Despite the psychrotrophs growth increase, according to temperature rise, the lipolytic metabolism was not consistent and presented the lower index at 8 °C, suggesting an intensification of the proteolytic compensatory activity at this temperature. © 2013 Institute of Food Technologists®.

Occurrence of mycobacteria in bovine milk samples from both individual and collective bulk tanks at farms and informal markets in the southeast region of Sao Paulo, Brazil

Background: Mycobacterium spp. is one of the most important species of zoonotic pathogens that can be transmitted from cattle to humans. The presence of these opportunistic, pathogenic bacteria in bovine milk has emerged as a public-health concern, especially among individuals who consume raw milk and related dairy products. To address this concern, the Brazilian control and eradication program focusing on bovine tuberculosis, was established in 2001. However, bovine tuberculosis continues to afflict approximately 1,3 percent of the cattle in Brazil. In the present study, 300 samples of milk from bovine herds, obtained from both individual and collective bulk tanks and informal points of sale, were cultured on Löwenstein-Jensen and Stonebrink media. Polymerase chain reaction (PCR)-based tests and restriction-enzyme pattern analysis were then performed on the colonies exhibiting phenotypes suggestive of Mycobacterium spp., which were characterized as acid-fast bacilli.Results: Of the 300 bovine milk samples that were processed, 24 were positively identified as Mycobacterium spp.Molecular identification detected 15 unique mycobacterial species: Mycobacterium bovis, M. gordonae, M. fortuitum, M. intracellulare, M. flavescens, M. duvalii, M. haemophilum, M. immunogenum, M. lentiflavum, M. mucogenicum, M. novocastrense, M. parafortuitum, M. smegmatis, M. terrae and M. vaccae. The isolation of bacteria from the various locations occurred in the following proportions: 9 percent of the individual bulk-tank samples, 7 percent of the collective bulk-tank samples and 8 percent of the informal-trade samples. No statistically significant difference was observed between the presence of Mycobacterium spp. in the three types of samples collected, the milk production profiles, the presence of veterinary assistance and the reported concerns about bovine tuberculosis prevention in the herds.Conclusion: The microbiological cultures associated with PCR-based identification tests are possible tools for the investigation of the presence of Mycobacterium spp. in milk samples. Using these methods, we found that the Brazilian population may be regularly exposed to mycobacteria by consuming raw bovine milk and related dairy products. These evidences reinforces the need to optimize quality programs of dairy products, to intensify the sanitary inspection of these products and the necessity of further studies on the presence of Mycobacterium spp. in milk and milk-based products. © 2013 Franco et al.; licensee BioMed Central Ltd.

Polymorphisms in the ghrelin gene and their associations with milk yield and quality in water buffaloes

Ghrelin is a gastrointestinal hormone that acts in releasing growth hormone and influences the body general metabolism. It has been proposed as a candidate gene for traits such as growth, carcass quality, and milk production of livestock because it influences feed intake. In this context, the aim of this study was to verify the existence of polymorphisms in the ghrelin gene and their associations with milk, fat and protein yield, and percentage in water buffaloes (Bubalus bubalis). A group of 240 animals was studied. Five primer pairs were used and 11 single nucleotide polymorphisms (SNP) were found in the ghrelin gene by sequencing. The animals were genotyped for 8 SNP by PCR-RFLP. The SNP g.960G>A and g.778C>T were associated with fat yield and the SNP g.905T>C was associated with fat yield and percentage and protein percentage. These SNP are located in intronic regions of DNA and may be in noncoding RNA sites or affect transcriptional efciency. The ghrelin gene in buffaloes influences milk fat and protein synthesis. The polymorphisms observed can be used as molecular markers to assist selection. © 2013 American Dairy Science Association.

Effect of shade and water sprinkling on physiological responses and milk yields of Holstein cows in a semi-arid region

The effects of shading and water sprinkling on physiological responses and milk production were studied in Holstein cows managed in an equatorial semi-arid region. Five cows were observed for 36 days during May and June 2009. We studied the effects of sun, shade, shade plus single sprinkling, shade plus sprinkling every 2. h and shade plus sprinkling every hour on physiological responses, i.e., rectal temperature, coat surface temperature and respiratory rate, and on milk production. Environmental variables, i.e., air temperature, black globe temperature, mean radiant temperature, relative humidity and wind speed, were recorded every 30. min throughout the observation period. The statistical analysis utilized a 5×5 replicated latin square, and the results showed reduced values in physiological variables in cows that were protected from short wave radiation and that received water sprinkling; milk yield increased by approximately 3.5. kg in such animals relative to those not receiving such treatment. However, the milk yield for the morning milking presented similar mean values regardless of whether the animals received shading and sprinkling benefits, which indicates that frequency of sprinkling and amount of time (permanence) in the shade were critical for the effect on milk yield. © 2013 Elsevier B.V.

Comparison of factor-analytic and reduced rank models for test-day milk yield in Gyr dairy cattle (Bos indicus)

We analyzed 46,161 monthly test-day records of milk production from 7453 first lactations of crossbred dairy Gyr (Bos indicus) x Holstein cows. The following seven models were compared: standard multivariate model (M10), three reduced rank models fitting the first 2, 3, or 4 genetic principal components, and three models considering a 2-, 3-, or 4-factor structure for the genetic covariance matrix. Full rank residual covariance matrices were considered for all models. The model fitting the first two principal components (PC2) was the best according to the model selection criteria. Similar phenotypic, genetic, and residual variances were obtained with models M10 and PC2. The heritability estimates ranged from 0.14 to 0.21 and from 0.13 to 0.21 for models M10 and PC2, respectively. The genetic correlations obtained with model PC2 were slightly higher than those estimated with model M10. PC2 markedly reduced the number of parameters estimated and the time spent to reach convergence. We concluded that two principal components are sufficient to model the structure of genetic covariances between test-day milk yields. © FUNPEC-RP.

Genetic Parameters for Productive and Reproductive Traits for Milk Buffalo in Brazil

Knowing the genetic parameters of productive and reproductive traits in milking buffaloes is essential for planning and implementing of a program genetic selection. In Brazil, this information is still scarce. The objective of this study was to verify the existence of genetic variability in milk yield of buffaloes and their constituents, and reproductive traits for the possibility of application of the selection. A total of 9,318 lactations records from 3,061 cows were used to estimate heritabilities for milk yield (MY), fat percentage (%F), protein percentage (%P), length of lactation (LL), age of first calving (AFC) and calving interval (CI) and the genetic correlations among traits MY, %F and %P. The (co) variance components were estimated using multiple-trait analysis by Bayesian inference method, applying an animal model, through Gibbs sampling. The model included the fixed effects of contemporary groups (herd-year and calving season), number of milking (2 levels), and age of cow at calving as (co) variable (quadratic and linear effect). The additive genetic, permanent environmental, and residual effects were included as random effects in the model. Estimated heritability values for MY, % F, % P, LL, AFC and CI were 0.24, 0.34, 0.40, 0.09, 0.16 and 0.05, respectively. The genetic correlation estimates among MY and % F, MY and % P and % F and % P were -0.29, -0.18 and 0.25, respectively. The production of milk and its constituents showed enough genetic variation to respond to a selection program. Negative estimates of genetic correlations between milk production and its components suggest that selection entails a reduction in the other.

Physiological and production response of dairy goats bred in a tropical climate

The aim of this work was to determine the adaptability of Saanen and A1/2Saanen x A1/2Anglo-Nubian (A1/2SA1/2AN) goats bred in tropical climates. The study included 30 goats, 15 Saanen and 15 A1/2SA1/2AN. The data was collected during the rainy and dry seasons. During the whole experimental period, the environment variables were recorded, as well as rectal temperature (RT), superficial temperature (ST), respiratory rate (RR) and heart rate (HR) and milk production (MP). The adaptability coefficient (AC) was calculated for both genotypes. The averages were evaluated by ANOVA at 5 % probability. There was a genotype and period of year effect, as well as the interaction genotype x period of year. Pearson's simple correlation analysis was then carried out between milk production and physiological and environment variables. There was a statistical difference (p < 0.05) between the seasons for RT, ST and RR. RT, RR and HR were lower for A1/2SA1/2AN than Saanen goats, regardless of the season. MP was greater in the dry season (p < 0.05) (2.52 A +/- 0.50 kg/day for A1/2SA1/2AN and 2.41 A +/- 0.38 kg/day for Saanen) than the rainy season (2.17 A +/- 0.27 kg/day for A1/2SA1/2AN and 2.28 A +/- 0.53 kg/day for Saanen). The MP correlations were very significant (p < 0.05), however low and negative, where it was higher when correlated with RR in Saanen goats. Under the conditions of the present study, it is concluded that the goats were influenced by climatic factors, where the rainy period was more likely to cause thermal stress in the animals.

A comparative study of production performance and animal health practices in organic and conventional dairy systems

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Predicting breeding values for milk yield of Guzera (Bos indicus) cows using random regression models

Given the importance of Guzera breeding programs for milk production in the tropics, the objective of this study was to compare alternative random regression models for estimation of genetic parameters and prediction of breeding values. Test-day milk yields records (TDR) were collected monthly, in a maximum of 10 measurements. The database included 20,524 records of first lactation from 2816 Guzera cows. TDR data were analyzed by random regression models (RRM) considering additive genetic, permanent environmental and residual effects as random and the effects of contemporary group (CG), calving age as a covariate (linear and quadratic effects) and mean lactation curve as fixed. The genetic additive and permanent environmental effects were modeled by RRM using Wilmink, All and Schaeffer and cubic B-spline functions as well as Legendre polynomials. Residual variances were considered as heterogeneous classes, grouped differently according to the model used. Multi-trait analysis using finite-dimensional models (FDM) for testday milk records (TDR) and a single-trait model for 305-days milk yields (default) using the restricted maximum likelihood method were also carried out as further comparisons. Through the statistical criteria adopted, the best RRM was the one that used the cubic B-spline function with five random regression coefficients for the genetic additive and permanent environmental effects. However, the models using the Ali and Schaeffer function or Legendre polynomials with second and fifth order for, respectively, the additive genetic and permanent environmental effects can be adopted, as little variation was observed in the genetic parameter estimates compared to those estimated by models using the B-spline function. Therefore, due to the lower complexity in the (co)variance estimations, the model using Legendre polynomials represented the best option for the genetic evaluation of the Guzera lactation records. An increase of 3.6% in the accuracy of the estimated breeding values was verified when using RRM. The ranks of animals were very close whatever the RRM for the data set used to predict breeding values. Considering P305, results indicated only small to medium difference in the animals' ranking based on breeding values predicted by the conventional model or by RRM. Therefore, the sum of all the RRM-predicted breeding values along the lactation period (RRM305) can be used as a selection criterion for 305-day milk production. (c) 2014 Elsevier B.V. All rights reserved.

Evaluation of the effect of supplementing fermented milk with quinoa flour on probiotic activity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Random regression models using different functions to model milk flow in dairy cows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Occurrence of heat waves in Triangulo Mineiro and Alto Paranaiba and their effect on milk yield and feed intake in cattle

The aim of this study was to estimate milk production and food consumption during the occurrence of heat waves in the Triangulo Mineiro and Alto Paranaiba, MG by means of bioclimatic zoning based on the Temperature and Humidity Index (THI). Therefore a history of heat wave occurrence between the years 2000-2010 was compiled. The decline in milk production (DMP) and reduced food consumption (RFC) were simulated in cities where periods of heat waves were identified. Frutal and Ituiutaba had the highest rate of heat wave occurrence per year. The DMP and RFC showed bioclimatic differences between the cities of Uberaba, Ituiutaba and Frutal. The cities with the best bioclimatic conditions were Sacramento and Patrocinio, as they presented a THI classified outside of the emergency range, with a night THI of below 76 and without heat waves. Therefore, the occurrence of heat waves can impair food intake and decrease milk production, thereby most effectively demonstrating the effects of thermal stress on dairy cows in the Triangulo Mineiro and Alto Paranaiba, MG region.

Short communication: Variable number of tandem repeat polymorphisms in DGAT1 gene of buffaloes (Bubalus bubalis) is associated with milk constituents

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)