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.

Heat stress on breeding value prediction for milk yield and composition of a Brazilian Holstein cattle population

Due to the high milk production of Holstein cows, many countries have chosen to import semen to improve local dairy herds. This strategy would be more effective if this semen was used in the same environment conditions in which the bulls were selected. If the effect of genotype by environment (G × E) interaction is not considered, the estimated breeding values (EBVs) may vary, potentially reducing the selection response. We evaluate the impact of heat stress on selection for milk yield and composition of Holstein cows using random regression models. To verify the interference of heat stress in milk yield (MY) and composition traits (fat, protein, total saturated, and total unsaturated fatty acids content in milk), temperature-humidity index (THI) on test-day milk records was used. The threshold value to divide the environments using test-day information from Brazilian Holstein cows was 72 units of THI, i.e., < 72 represented no heat stress and > 72 represented heat stress. Legendre polynomials of second-order (Leg 2) model and two lactation points (33 and 122 DIM) were used to estimate heritabilities and EBVs for five important dairy traits. The heritabilities of milk components and fatty acids were low (0.09-0.29), regardless of lactation period and degree of heat stress, with the exception of protein content (0.30-0.35). Fat content was the only milk component that was reduced according to the degree of heat stress and lactation period. The EBVs tended to decrease in heat stress conditions, thus animals with high genetic potential demonstrated evidence of G × E interaction. However, acclimatization of dairy cows to heat stress in the farm production systems may have been responsible for the low differences among genetic parameters and EBVs with and without heat stress found in this study.

An Integrative Genomic Prediction Approach for Predicting Buffalo Milk Traits by Incorporating Related Cattle QTLs

Background: The 90K Axiom Buffalo SNP Array is expected to improve and speed up various genomic analyses for the buffalo (Bubalus bubalis). Genomic prediction is an effective approach in animal breeding to improve selection and reduce costs. As buffalo genome research is lagging behind that of the cow and production records are also limited, genomic prediction performance will be relatively poor. To improve the genomic prediction in buffalo, we introduced a new approach (pGBLUP) for genomic prediction of six buffalo milk traits by incorporating QTL information from the cattle milk traits in order to help improve the prediction performance for buffalo. Results: In simulations, the pGBLUP could outperform BayesR and the GBLUP if the prior biological information (i.e., the known causal loci) was appropriate; otherwise, it performed slightly worse than BayesR and equal to or better than the GBLUP. In real data, the heritability of the buffalo genomic region corresponding to the cattle milk trait QTLs was enriched (fold of enrichment > 1) in four buffalo milk traits (FY270, MY270, PY270, and PM) when the EBV was used as the response variable. The DEBV as the response variable yielded more reliable genomic predictions than the traditional EBV, as has been shown by previous research. The performance of the three approaches (GBLUP, BayesR, and pGBLUP) did not vary greatly in this study, probably due to the limited sample size, incomplete prior biological information, and less artificial selection in buffalo. Conclusions: To our knowledge, this study is the first to apply genomic prediction to buffalo by incorporating prior biological information. The genomic prediction of buffalo traits can be further improved with a larger sample size, higher-density SNP chips, and more precise prior biological information.

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.

Integrating High-Throughput Phenotyping and Statistical Genomic Methods to Genetically Improve Longitudinal Traits in Crops

The rapid development of remote sensing in agronomic research allows the dynamic nature of longitudinal traits to be adequately described, which may enhance the genetic improvement of crop efficiency. For traits such as light interception, biomass accumulation, and responses to stressors, the data generated by the various high-throughput phenotyping (HTP) methods requires adequate statistical techniques to evaluate phenotypic records throughout time. As a consequence, information about plant functioning and activation of genes, as well as the interaction of gene networks at different stages of plant development and in response to environmental stimulus can be exploited. In this review, we outline the current analytical approaches in quantitative genetics that are applied to longitudinal traits in crops throughout development, describe the advantages and pitfalls of each approach, and indicate future research directions and opportunities.

Genetic correlations and trends for traits of economic importance in dairy buffalo

Context The planning and execution of selection programs requires estimates of the genetic correlations between traits. As genetic change is achieved for a given trait, it is important to consider possible genetic changes for other traits. Understanding the magnitude and direction of genetic correlations can assist in selection decisions. Aims The aim of the present study was to estimate the genetic correlations of reproductive traits with productive traits and with percentages of fat and protein in the milk of dairy buffalo. Additionally, genetic trends were estimated for the traits under study over the years. Methods Data from 11530 complete lactations of 3431 female buffalo were used. The following traits were analysed: milk, fat and protein yields; percentages of fat and protein; age at first calving; and calving interval. The (co)variance components were estimated by Bayesian inference in multi-trait analyses, considering a linear animal model. To calculate the genetic trends, the average annual genetic values were regressed on the year of birth. Key results The means of genetic correlations estimated between reproductive (age at first calving and calving interval) and productive (milk, fat and protein yields) traits were positive, but of moderate to low magnitude. The association between the reproductive and milk quality (fat and protein percentages) traits were negative and of low magnitude. Genetic trends for the productive traits were positive (5.25 ± 0.63, 0.15 ± 0.034 and 0.09 ± 0.038 kg/year for milk, fat and protein yields respectively). Genetic trends for the reproductive traits of age at first calving and calving interval increased by 0.47 ± 0.09 and 0.48 ± 0.10 days/year respectively. In terms of milk quality, however, the percentages of fat and protein decreased by 0.016 ± 0.003 and 0.011 ± 0.001%/year respectively. Conclusions Genetic gains in productive traits may elevate the number of days at first calving and extend the calving interval, in addition to leading to the production of milk of lower quality. Implications The use of a multi-trait selection index is an alternative, as it combines information from different sources, such that an optimal selection criterion can be achieved over time by virtue of its emphasis on appropriate weighting for all traits.

Genome-wide association study applied to type traits related to milk yield in water buffaloes (Bubalus bubalis)

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

Invited review: Advances and applications of random regression models: From quantitative genetics to genomics

An important goal in animal breeding is to improve longitudinal traits; that is, traits recorded multiple times during an individual's lifetime or physiological cycle. Longitudinal traits were first genetically evaluated based on accumulated phenotypic expression, phenotypic expression at specific time points, or repeatability models. Until now, the genetic evaluation of longitudinal traits has mainly focused on using random regression models (RRM). Random regression models enable fitting random genetic and environmental effects over time, which results in higher accuracy of estimated breeding values compared with other statistical approaches. In addition, RRM provide insights about temporal variation of biological processes and the physiological implications underlying the studied traits. Despite the fact that genomic information has substantially contributed to increase the rates of genetic progress for a variety of economically important traits in several livestock species, less attention has been given to longitudinal traits in recent years. However, including genomic information to evaluate longitudinal traits using RRM is a feasible alternative to yield more accurate selection and culling decisions, because selection of young animals may be based on the complete pattern of the production curve with higher accuracy compared with the use of traditional parent average (i.e., without genomic information). Moreover, RRM can be used to estimate SNP effects over time in genome-wide association studies. Thus, by analyzing marker associations over time, regions with higher effects at specific points in time are more likely to be identified. Despite the advances in applications of RRM in genetic evaluations, more research is needed to successfully combine RRM and genomic information. Future research should provide a better understanding of the temporal variation of biological processes and their physiological implications underlying the longitudinal traits.

How to improve breeding value prediction for feed conversion ratio in the case of incomplete longitudinal body weights

With the development of automatic self-feeders, repeated measurements of feed intake are becoming easier in an increasing number of species. However, the corresponding BW are not always recorded, and these missing values complicate the longitudinal analysis of the feed conversion ratio (FCR). Our aim was to evaluate the impact of missing BW data on estimations of the genetic parameters of FCR and ways to improve the estimations. On the basis of the missing BW profile in French Large White pigs (male pigs weighed weekly, females and castrated males weighed monthly), we compared 2 different ways of predicting missing BW, 1 using a Gompertz model and 1 using a linear interpolation. For the first part of the study, we used 17,398 weekly records of BW and feed intake recorded over 16 consecutive weeks in 1,222 growing male pigs. We performed a simulation study on this data set to mimic missing BW values according to the pattern of weekly proportions of incomplete BW data in females and castrated males. The FCR was then computed for each week using observed data (obser_FCR), data with missing BW (miss_FCR), data with BW predicted using a Gompertz model (Gomp_FCR), and data with BW predicted by linear interpolation (interp_FCR). Heritability (h) was estimated, and the EBV was predicted for each repeated FCR using a random regression model. In the second part of the study, the full data set (males with their complete BW records, castrated males and females with missing BW) was analyzed using the same methods (miss_FCR, Gomp_FCR, and interp_FCR). Results of the simulation study showed that h were overestimated in the case of missing BW and that predicting BW using a linear interpolation provided a more accurate estimation of h and of EBV than a Gompertz model. Over 100 simulations, the correlation between obser_EBV and interp_EBV, Gomp_EBV, and miss_EBV was 0.93 ± 0.02, 0.91 ± 0.01, and 0.79 ± 0.04, respectively. The heritabilities obtained with the full data set were quite similar for miss_FCR, Gomp_FCR, and interp_FCR. In conclusion, when the proportion of missing BW is high, genetic parameters of FCR are not well estimated. In French Large White pigs, in the growing period extending from d 65 to 168, prediction of missing BW using a Gompertz growth model slightly improved the estimations, but the linear interpolation improved the estimation to a greater extent. This result is due to the linear rather than sigmoidal increase in BW over the study period.