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.

Lázaro SF ，Tonhati H ，Oliveira HR ，Silva AA ，Scalez DCB ，Nascimento AV ，Santos DJA ，Stefani G ，Carvalho IS ，Sandoval AF ，Brito LF ... -  《-》

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models.

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Lázaro SF ，Tonhati H ，Oliveira HR ，Silva AA ，Nascimento AV ，Santos DJA ，Stefani G ，Brito LF ... -  《-》

Random regression models using B-splines functions provide more accurate genomic breeding values for milk yield and lactation persistence in Murrah buffaloes.

There is a great worldwide demand for cheese made with buffalo milk, due to its flavour and nutritional properties. In this context, there is a need for increasing the efficiency of buffalo milk production (including lactation persistence), which can be achieved through genomic selection. The most used methods for the genetic evaluation of longitudinal data, such as milk-related traits, are based on random regression models (RRM). The choice of the best covariance functions and polynomial order for modelling the random effects is an important step to properly fit RRM. To our best knowledge, there are no studies evaluating the impact of the order and covariance function (Legendre polynomials-LEG and B-splines-BSP) used to fit RRM for genomic prediction of breeding values in dairy buffaloes. Therefore, the main objectives of this study were to estimate variance components and evaluate the performance of LEG and BSP functions of different orders on the predictive ability of genomic breeding values for the first three lactations of milk yield (MY1, MY2, and MY3) and lactation persistence (LP1, LP2, and LP3) of Brazilian Murrah. Twenty-two models for each lactation were contrasted based on goodness of fit, genetic parameter estimates, and predictive ability. Overall, the models of higher orders of LEG or BSP had a better performance based on the deviance information criterion (DIC). The daily heritability estimates ranged from 0.01 to 0.30 for MY1, 0.08 to 0.42 for MY2, and from 0.05 to 0.47 for MY3. For lactation persistence (LP), the heritability estimates ranged from 0.09 to 0.32 for LP1, from 0.15 to 0.33 for LP2, and from 0.06 to 0.32 for LP3. In general, the curves plotted for variance components and heritability estimates based on BSP models presented lower oscillation along the lactation trajectory. Similar predictive ability was observed among the models. Considering a balance between the complexity of the model, goodness of fit, and credibility of the results, RRM using quadratic B-splines functions based on four or five segments to model the systematic, additive genetic, and permanent environment curves provide better fit with no significant differences between genetic variances estimates, heritabilities, and predictive ability for the genomic evaluation of dairy buffaloes.

Silva AA ，Brito LF ，Silva DA ，Lazaro SF ，Silveira KR ，Stefani G ，Tonhati H ... -  《-》

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes.

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal.

Borquis RR ，Neto FR ，Baldi F ，Hurtado-Lugo N ，de Camargo GM ，Muñoz-Berrocal M ，Tonhati H ... -  《-》

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle.

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.

Oliveira HR ，Lourenco DAL ，Masuda Y ，Misztal I ，Tsuruta S ，Jamrozik J ，Brito LF ，Silva FF ，Cant JP ，Schenkel FS ... -  《-》