H. Tonhati

GENETIC PARAMETER ESTIMATES FOR BUFFALO MILK YIELD, MILK QUALITY AND MOZZARELLA PRODUCTION AND BAYESIAN INFERENCE ANALYSIS OF THEIR RELATIONSHIPS

Buffalo milk has excellent physical and chemical qualities as a consequence of the high percentage of constituents. This milk property is desirable for the dairy industry because it facilitates manufacture of mozzarella cheese. We estimated genetic parameters for milk yield, milk fat and protein and their effects on mozzarella cheese production using Bayesian inference. Using information from 4907 lactation records of buffaloes, genetic and non-genetic parameters were estimated for accumulated 305-day milk yield (MY), milk fat (%F) and protein (%P) percentages and mozzarella production per lactation (MP). The (co)variance components were obtained by Bayesian inference using a multiple trait model, which included as fixed effects contemporary group, milking number and buffalo age at calving as covariables (linear and quadratic), along with the additive genetic, permanent environmental and residual random effects. Mean a posteriori heritability distributions for MY, %F, %P, and MP were 0.25, 0.30, 0.38, and 0.23, respectively. The genetic correlation estimates between MY with %P and %F were negative and moderate. Positive genetic correlation estimates varying from 0.19 (%P/MP) to 0.95 (MY/MP) were obtained among the traits. Milk yield, milk components, and mozzarella production in Murrah buffaloes have enough genetic variation for selection purposes. We conclude that selection to increase milk yield would be effective in improving mozzarella production.

Association between single-nucleotide polymorphisms and milk production traits in buffalo

The aim of this study was to identify single-nucleotide polymorphisms (SNPs) in buffaloes associated with milk yield and content, in addition to somatic cell scores based on the cross-species transferability of SNPs from cattle to buffalo. A total of 15,745 SNPs were analyzed, of which 1562 showed 1% significance and 4742 with 5% significance, which were associated for all traits studied. After application of Bonferronis correction for multiple tests of the traits analyzed, we found 2 significant SNPs placed on cattle chromosomes BTA15 and BTA20, which are homologous to buffalo chromosomes BBU16 and BBU19, respectively. In this genome association study, we found several significant SNPs affecting buffalo milk production and quality. Furthermore, the use of the high-density bovine BeadChip was suitable for genomic analysis in buffaloes. Although extensive chromosome arm homology was described between cattle and buffalo, the exact chromosomal position of SNP markers associated with these economically important traits in buffalo can be determined only through buffalo genome sequencing.

Genetic parameters estimates for milk, fat and protein yield analyzed by test day models for Murrah buffaloes in Brazil

Abstract The objectives of this study were to estimate genetic parameters for test-day milk, fat and protein yields, in Murrah buffaloes. In this study 4,757 complete lactations of Murrah buffaloes were analyzed. The (co) variance components were estimated by restricted maximum likelihood using MTDFREML software. The bi-trait animal test-day models included genetic additive direct and permanent environment effects, as random effects, and the fixed effects of contemporary group (herds-year-month of control) and age of the cow at calving as linear and quadratic covariable. The heritability estimate at first control was 0.19, increased until the third control (0.24), decreasing thereafter, reaching the lowest value at the ninth control (0.09). The highest heritability estimates for fat and protein yield were 0.23 (first control) and 0.33 (third control), respectively. For milk yield, genetic and phenotypic correlation estimates ranged from 0.37 to 0.99 and from 0.52 to 0.94, respectively. Genetic correlations were higher than phenotypic ones. For fat and protein yields, genetic correlation estimates ranged from 0.42 to 0.97.

Polymorphism in the A2M gene associated with high-quality milk in Murrah buffaloes (Bubalus bubalis)

The study of genes associated with host defense mechanisms, such as the A2M gene, plays a critical role in preventing diseases that reduce milk yield and its constituents. The aim of this study was to identify polymorphisms in the A2M gene in Murrah buffaloes (Bubalus bubalis), and investigate their associations with milk yield, fat and protein production, fat and protein percentages, and somatic cell count. Hair follicle samples of 136 animals were collected for DNA extraction, and polymorphisms were identified by polymerase chain reactions and sequencing. Statistical analyses were performed to ascertain the allelic and genotypic frequencies, the Hardy-Weinberg equilibrium, and association analysis was conducted between the polymorphisms and the traits studied. Comparative analysis between buffalo and bovine sequences revealed seven nucleotide substitutions. Alignments among the buffalo sequences identified three single nucleotide polymorphisms (SNPs), including one in exon 29, g.241A>G, which was used in subsequent statistical analyses. A Hardy-Weinberg test indicated that this SNP was in equilibrium in this population, and was significantly associated (P < 0.05) with fat production and fat and protein percentages. Therefore, this SNP can be used as a molecular marker for these traits.

Multivariate analysis of test-day and total milk yield in goats

The objective of this study was to estimate genetic parameters for 305-day cumulative milk yield (MY305) and its association with test-day milk yield (TDMY) in Saanen and Alpine goats in order to provide information that allows the use of TDMY as selection criteria. This was done using standard multi-trait and reduced rank models. Data from 1157 lactations, including the first three kiddings, and 5435 test-day records from 683 Saanen and 449 Alpine goats were used. MY305 was analyzed together with TDMY by multi-trait analysis, from the first to tenth test-day, using records of the first three lactations as repeated measures. Three multi-trait models were used: a standard (SM) and two reduced rank models that fitted the first two (PC2) and three (PC3) genetic principal components. Akaike and Schwarz Bayesian information criteria were used to compare models. Heritability for TDMY estimated with the SM ranged from 0.20 to 0.66, whereas the range calculated from the PC2 model was 0.16 to 0.63. Genetic correlations between TDMY and MY305 were positive and moderate to high, ranging from 0.56 to 0.98 when estimated with the SM, and 0.91 to 1.00 when estimated with the PC2. The standard multi-trait model produced estimates that were more accurate than the reduced rank models. Although the SM provided the worst fit according to the two model selection criteria, it was the best in this dataset.