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 ... -  《-》

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 ... -  《-》

Genome-wide association for milk production traits and somatic cell score in different lactation stages of Ayrshire, Holstein, and Jersey dairy cattle.

We performed genome-wide association analyses for milk, fat, and protein yields and somatic cell score based on lactation stages in the first 3 parities of Canadian Ayrshire, Holstein, and Jersey cattle. The genome-wide association analyses were performed considering 3 different lactation stages for each trait and parity: from 5 to 95, from 96 to 215, and from 216 to 305 d in milk. Effects of single nucleotide polymorphisms (SNP) for each lactation stage, trait, parity, and breed were estimated by back-solving the direct breeding values estimated using the genomic best linear unbiased predictor and single-trait random regression test-day models containing only the fixed population average curve and the random genomic curves. To identify important genomic regions related to the analyzed lactation stages, traits, parities and breeds, moving windows (SNP-by-SNP) of 20 adjacent SNP explaining more than 0.30% of total genetic variance were selected for further analyses of candidate genes. A lower number of genomic windows with a relatively higher proportion of the explained genetic variance was found in the Holstein breed compared with the Ayrshire and Jersey breeds. Genomic regions associated with the analyzed traits were located on 12, 8, and 15 chromosomes for the Ayrshire, Holstein, and Jersey breeds, respectively. Especially for the Holstein breed, many of the identified candidate genes supported previous reports in the literature. However, well-known genes with major effects on milk production traits (e.g., diacylglycerol O-acyltransferase 1) showed contrasting results among lactation stages, traits, and parities of different breeds. Therefore, our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the analyzed traits across breeds, parities, and lactation stages. Further functional studies are needed to validate our findings in independent populations.

Oliveira HR ，Cant JP ，Brito LF ，Feitosa FLB ，Chud TCS ，Fonseca PAS ，Jamrozik J ，Silva FF ，Lourenco DAL ，Schenkel FS ... -  《-》

Genomic prediction of lactation curves for milk, fat, protein, and somatic cell score in Holstein cattle.

Application of random regression models (RRM) in a 2-step genomic prediction might be a feasible way to select young animals based on the complete pattern of the lactation curve. In this context, the prediction reliability and bias of genomic estimated breeding value (GEBV) for milk, fat, and protein yields and somatic cell score over days in milk (DIM) using a 2-step genomic approach were investigated. In addition, the effect of including cows in the training and validation populations was investigated. Estimated breeding values for each DIM (from 5 to 305 d) from the first 3 lactations of Holstein animals were deregressed and used as pseudophenotypes in the second step. Individual additive genomic random regression coefficients for each trait were predicted using RRM and genomic best linear unbiased prediction and further used to derive GEBV for each DIM. Theoretical reliabilities of GEBV obtained by the RRM were slightly higher than theoretical reliabilities obtained by the accumulated yield up to 305 d (P305). However, validation reliabilities estimated for GEBV using P305 were higher than for GEBV using RRM. For all traits, higher theoretical and validation reliabilities were estimated when incorporating genomic information. Less biased GEBV estimates were found when using RRM compared with P305, and different validation reliability and bias patterns for GEBV over time were observed across traits and lactations. Including cows in the training population increased the theoretical reliabilities and bias of GEBV; nonetheless, the inclusion of cows in the validation population does not seem to affect the regression coefficients and the theoretical reliabilities. In summary, the use of RRM in 2-step genomic prediction produced fairly accurate GEBV over the entire lactation curve for all analyzed traits. Thus, selecting young animals based on the pattern of lactation curves seems to be a feasible alternative in genomic selection of Holstein cattle for milk production traits.

Oliveira HR ，Brito LF ，Silva FF ，Lourenco DAL ，Jamrozik J ，Schenkel FS ... -  《-》

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 ... -  《-》