Genome-wide association study for stayability at different calvings in Nellore beef cattle.

BACKGROUNDING: Stayability, which may be defined as the probability of a cow remaining in the herd until a reference age or at a specific number of calvings, is usually measured late in the animal's life. Thus, if used as selection criteria, it will increase the generation interval and consequently might decrease the annual genetic gain. Measuring stayability at an earlier age could be a reasonable strategy to avoid this problem. In this sense, a better understanding of the genetic architecture of this trait at different ages and/or at different calvings is important. This study was conducted to identify possible regions with major effects on stayability measured considering different numbers of calvings in Nellore cattle as well as pathways that can be involved in its expression throughout the female's productive life. RESULTS: The top 10 most important SNP windows explained, on average, 17.60% of the genetic additive variance for stayability, varying between 13.70% (at the eighth calving) and 21% (at the fifth calving). These SNP windows were located on 17 chromosomes (1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 18, 19, 20, 27, and 28), and they harbored a total of 176 annotated genes. The functional analyses of these genes, in general, indicate that the expression of stayability from the second to the sixth calving is mainly affected by genetic factors related to reproductive performance, and nervous and immune systems. At the seventh and eighth calvings, genes and pathways related to animal health, such as density bone and cancer, might be more relevant. CONCLUSION: Our results indicate that part of the target genomic regions in selecting for stayability at earlier ages (from the 2th to the 6th calving) would be different than selecting for this trait at later ages (7th and 8th calvings). While the expression of stayability at earlier ages appeared to be more influenced by genetic factors linked to reproductive performance together with an overall health/immunity, at later ages genetic factors related to an overall animal health gain relevance. These results support that selecting for stayability at earlier ages (perhaps at the second calving) could be applied, having practical implications in breeding programs since it could drastically reduce the generation interval, accelerating the genetic progress.

Genetic Perspectives on Feed Event, Meal and Feed Efficiency Traits in Bos taurus indicus Beef Cattle.

Electronic feeders record feeding behaviour as feed events by tracking the animal's in-out visits to the feeder. Another way to measure feeding behaviour is based on meals. However, the two approaches provide different outcomes. The objectives of this study were to estimate genetic parameters (heritabilities and genetic and phenotypic correlations) for feed event and meal traits, and their genetic and phenotypic correlations with feed efficiency traits in Nellore cattle. The present study analysed six feed event traits (DMIFE: dry matter intake per feed event, FED: feed event duration, TBFE: time between feed events, FTd: feeding time per day, FEd: feed events per day, and FR: feeding rate), six meal traits (DMIME: DMI per meal, MED: meal duration, TBME: time between meals, MC: meal criterion, MTd: meal time per day, and MEd: meals per day), and three feed efficiency traits (ADG: average daily gain, DMI, and RFI: residual feed intake). The traits were measured in feed efficiency tests of Nellore cattle (age = 280 ± 41 days and body weight = 258 ± 47 kg at enrolment). The MC was calculated for each animal and ranged from 1.70 to 64.0 min, i.e., any pair of feed events separated by less than the MC value was considered part of the same meal. The heritabilities and correlations were estimated by fitting univariate and bivariate animal models, respectively, using single-step genomic BLUP. The highest heritabilities for feed event traits were 0.35 ± 0.06 (FED), 0.39 ± 0.06 (FTd), and 0.50 ± 0.05 (FTd), and for meal traits were 0.31 ± 0.06 (MED) and 0.45 ± 0.06 (MTd). The genetic correlation between feed event traits and meal traits were weak. FR, FED, and FTd had moderate genetic correlations with RFI (-0.56 ± 0.11, 0.44 ± 0.11, 0.60 ± 0.08, respectively). These results indicate that more efficient animals spent less time at the feeder per feed event and per day, and eat faster compared to less efficient animals. In conclusion, feed event and meal traits must be treated as distinct groups of traits since the genetic and phenotypic correlations were, in general, weak to moderate. Among feed event versus meal traits, feed event traits are more favourable to explain the genetic relationships of feeding behaviour with feed efficiency-related traits.

Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle.

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by ß-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CVr) and herd-out (CVh). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R2 value from 5% to 75% for CVr and 2% to 139% for CVh, whereas the stacking ensemble increased the R2 value from 4% to 70% for CVr and 4% to 150% for CVh. Considering the best model, with the CVr scenario, good prediction accuracies were obtained for glucose (R2 = 0.81), urea (R2 = 0.73), albumin (R2 = 0.75), total reactive oxygen metabolites (R2 = 0.79), total thiol groups (R2 = 0.76), ceruloplasmin (R2 = 0.74), total proteins (R2 = 0.81), globulins (R2 = 0.87), and Na (R2 = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

Infrared spectroscopy coupled with machine learning algorithms for predicting the detailed milk mineral profile in dairy cattle.

Milk minerals are not only essential components for human health, but they can be informative for milk quality and cow's health. Herein, we investigated the feasibility of Fourier Transformed mid Infrared (FTIR) spectroscopy for the prediction of a detailed panel of 17 macro, trace, and environmental elements in bovine milk, using partial least squares regression (PLS) and machine learning approaches. The automatic machine learning significantly outperformed the PLS regression in terms of prediction performances of the mineral elements. For macrominerals, the R2 ranged from 0.59 to 0.78. Promising predictability was achieved for Cu and B (R2 = 0.66 and 0.74, respectively) and more moderate ones for Fe, Mn, Zn, and Al (R2 from 0.48 to 0.58). These results provide a reliable basis for a rapid and cost-effective quantification of these traits, serving as a resource for dairy farmers seeking to enhance the quality of milk production and optimize cheese properties.

Structural equation models to infer relationships between energy-related blood metabolites and milk daily energy output in Holstein cows.

During lactation, high-yielding cows experience metabolic disturbances due to milk production. Metabolic monitoring offers valuable insights into how cows manage these challenges throughout the lactation period, making it a topic of considerable interest to breeders. In this study, we used Bayesian networks to uncover potential dependencies among various energy-related blood metabolites, i.e., glucose, urea, beta-hydroxybutyrate (BHB), non-esterified fatty acids (NEFA), cholesterol (CHOL), and daily milk energy output (dMEO) in 1,254 Holstein cows. The inferred causal structure was then incorporated into structural equation models (SEM) to estimate heritabilities and additive genetic correlations among these phenotypes using both pedigree and genotypes from a 100k chip. Dependencies among traits were determined using the Hill-Climbing algorithm, implemented with the posterior distribution of the residuals obtained from the standard multiple-trait model. These identified relationships were then used to construct the SEM, considering both direct and indirect relationships. The relevant dependencies and path coefficients obtained, expressed in units of measurement variation of 1σ, were as follows: dMEO → CHOL (0.181), dMEO → BHB (-0.149), dMEO → urea (0.038), glucose → BHB (-0.55), glucose → urea (-0.194), CHOL → urea (0.175), BHB → urea (-0.049), and NEFA → urea (-0.097). Heritabilities for traits of concern obtained with SEM ranged from 0.09 to 0.2. Genetic correlations with a minimum 95% probability (P) of the posterior mean being >0 for positive means or <0 for negative means include those between dMEO and glucose (-0.583, P = 100), dMEO and BHB (0.349, P = 99), glucose and CHOL (0.325, P = 100), glucose and NEFA (-0.388, P = 100), and NEFA and BHB (0.759, P = 100). The results of this analysis revealed the existence of recursive relationships among the energy-related blood metabolites and dMEO. Understanding these connections is paramount for establishing effective genetic selection strategies, enhancing production and animal welfare.

Dairy cows face significant metabolic challenges during the different phases of the production cycle. One of their primary issues revolves around meeting the elevated energy demands, encompassing essential maintenance requirements and the energy required for milk production. In this context, monitoring of energy-related blood metabolites along lactation may be useful to detect metabolic disturbances. To date, no study is available on the investigation of putative recursive relationships among energy-related blood metabolites and milk daily energy output. The use of structural equation models presents an innovative approach to enhance our understanding of these complex relationships. This approach allows us to gain a better insight into the metabolic pathways and processes involved in energy metabolism. The findings of this study uncovered recursive relationships between energy-related blood metabolites and daily milk energy output. Grasping these interactions is crucial for developing effective breeding strategies focused on selecting more resistant and resilient cows, which have the capacity to better cope with metabolic distress along lactation.

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle.

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

Lactobacillus gasseri LG-G12 Restores Gut Microbiota and Intestinal Health in Obesity Mice on Ceftriaxone Therapy.

Gut microbiota imbalance is associated with the occurrence of metabolic diseases such as obesity. Thus, its modulation is a promising strategy to restore gut microbiota and improve intestinal health in the obese. This paper examines the role of probiotics, antimicrobials, and diet in modulating gut microbiota and improving intestinal health. Accordingly, obesity was induced in C57BL/6J mice, after which they were redistributed and fed with an obesogenic diet (intervention A) or standard AIN-93 diet (intervention B). Concomitantly, all the groups underwent a treatment phase with Lactobacillus gasseri LG-G12, ceftriaxone, or ceftriaxone followed by L. gasseri LG-G12. At the end of the experimental period, the following analysis was conducted: metataxonomic analysis, functional profiling of gut microbiota, intestinal permeability, and caecal concentration of short-chain fatty acids. High-fat diet impaired bacterial diversity/richness, which was counteracted in association with L. gasseri LG-G12 and the AIN-93 diet. Additionally, SCFA-producing bacteria were negatively correlated with high intestinal permeability parameters, which was further confirmed via functional profile prediction of the gut microbiota. A novel perspective on anti-obesity probiotics is presented by these findings based on the improvement of intestinal health irrespective of undergoing antimicrobial therapy or not.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle.

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

Improvement of Genomic Predictions in Small Breeds by Construction of Genomic Relationship Matrix Through Variable Selection.

Genomic selection has been increasingly implemented in the animal breeding industry, and it is becoming a routine method in many livestock breeding contexts. However, its use is still limited in several small-population local breeds, which are, nonetheless, an important source of genetic variability of great economic value. A major roadblock for their genomic selection is accuracy when population size is limited: to improve breeding value accuracy, variable selection models that assume heterogenous variance have been proposed over the last few years. However, while these models might outperform traditional and genomic predictions in terms of accuracy, they also carry a proportional increase of breeding value bias and dispersion. These mutual increases are especially striking when genomic selection is performed with a low number of phenotypes and high shrinkage value-which is precisely the situation that happens with small local breeds. In our study, we tested several alternative methods to improve the accuracy of genomic selection in a small population. First, we investigated the impact of using only a subset of informative markers regarding prediction accuracy, bias, and dispersion. We used different algorithms to select them, such as recursive feature eliminations, penalized regression, and XGBoost. We compared our results with the predictions of pedigree-based BLUP, single-step genomic BLUP, and weighted single-step genomic BLUP in different simulated populations obtained by combining various parameters in terms of number of QTLs and effective population size. We also investigated these approaches on a real data set belonging to the small local Rendena breed. Our results show that the accuracy of GBLUP in small-sized populations increased when performed with SNPs selected via variable selection methods both in simulated and real data sets. In addition, the use of variable selection models-especially those using XGBoost-in our real data set did not impact bias and the dispersion of estimated breeding values. We have discussed possible explanations for our results and how our study can help estimate breeding values for future genomic selection in small breeds.

Sustainable Intensification of Beef Production in the Tropics: The Role of Genetically Improving Sexual Precocity of Heifers.

Increasing productivity through continued animal genetic improvement is a crucial part of implementing sustainable livestock intensification programs. In Zebu cattle, the lack of sexual precocity is one of the main obstacles to improving beef production efficiency. Puberty-related traits are complex, but large-scale data sets from different "omics" have provided information on specific genes and biological processes with major effects on the expression of such traits, which can greatly increase animal genetic evaluation. In addition, genetic parameter estimates and genomic predictions involving sexual precocity indicator traits and productive, reproductive, and feed-efficiency related traits highlighted the feasibility and importance of direct selection for anticipating heifer reproductive life. Indeed, the case study of selection for sexual precocity in Nellore breeding programs presented here show that, in 12 years of selection for female early precocity and improved management practices, the phenotypic means of age at first calving showed a strong decreasing trend, changing from nearly 34 to less than 28 months, with a genetic trend of almost -2 days/year. In this period, the percentage of early pregnancy in the herds changed from around 10% to more than 60%, showing that the genetic improvement of heifer's sexual precocity allows optimizing the productive cycle by reducing the number of unproductive animals in the herd. It has a direct impact on sustainability by better use of resources. Genomic selection breeding programs accounting for genotype by environment interaction represent promising tools for accelerating genetic progress for sexual precocity in tropical beef cattle.

In-line near-infrared analysis of milk coupled with machine learning methods for the daily prediction of blood metabolic profile in dairy cattle.

Precision livestock farming technologies are used to monitor animal health and welfare parameters continuously and in real time in order to optimize nutrition and productivity and to detect health issues at an early stage. The possibility of predicting blood metabolites from milk samples obtained during routine milking by means of infrared spectroscopy has become increasingly attractive. We developed, for the first time, prediction equations for a set of blood metabolites using diverse machine learning methods and milk near-infrared spectra collected by the AfiLab instrument. Our dataset was obtained from 385 Holstein Friesian dairy cows. Stacking ensemble and multi-layer feedforward artificial neural network outperformed the other machine learning methods tested, with a reduction in the root mean square error of between 3 and 6% in most blood parameters. We obtained moderate correlations (r) between the observed and predicted phenotypes for γ-glutamyl transferase (r = 0.58), alkaline phosphatase (0.54), haptoglobin (0.66), globulins (0.61), total reactive oxygen metabolites (0.60) and thiol groups (0.57). The AfiLab instrument has strong potential but may not yet be ready to predict the metabolic stress of dairy cows in practice. Further research is needed to find out methods that allow an improvement in accuracy of prediction equations.

Comparison of Single-Breed and Multi-Breed Training Populations for Infrared Predictions of Novel Phenotypes in Holstein Cows.

In general, Fourier-transform infrared (FTIR) predictions are developed using a single-breed population split into a training and a validation set. However, using populations formed of different breeds is an attractive way to design cross-validation scenarios aimed at increasing prediction for difficult-to-measure traits in the dairy industry. This study aimed to evaluate the potential of FTIR prediction using training set combining specialized and dual-purpose dairy breeds to predict different phenotypes divergent in terms of biological meaning, variability, and heritability, such as body condition score (BCS), serum ß-hydroxybutyrate (BHB), and kappa casein (k-CN) in the major cattle breed, i.e., Holstein-Friesian. Data were obtained from specialized dairy breeds: Holstein (468 cows) and Brown Swiss (657 cows), and dual-purpose breeds: Simmental (157 cows), Alpine Grey (75 cows), and Rendena (104 cows), giving a total of 1461 cows from 41 multi-breed dairy herds. The FTIR prediction model was developed using a gradient boosting machine (GBM), and predictive ability for the target phenotype in Holstein cows was assessed using different cross-validation (CV) strategies: a within-breed scenario using 10-fold cross-validation, for which the Holstein population was randomly split into 10 folds, one for validation and the remaining nine for training (10-fold_HO); an across-breed scenario (BS_HO) where the Brown Swiss cows were used as the training set and the Holstein cows as the validation set; a specialized multi-breed scenario (BS+HO_10-fold), where the entire Brown Swiss and Holstein populations were combined then split into 10 folds, and a multi-breed scenario (Multi-breed), where the training set comprised specialized (Holstein and Brown Swiss) and dual-purpose (Simmental, Alpine Grey, and Rendena) dairy cows, combined with nine folds of the Holstein cows. Lastly a Multi-breed CV2 scenario was implemented, assuming the same number of records as the reference scenario and using the same proportions as the multi-breed. Within-Holstein, FTIR predictions had a predictive ability of 0.63 for BCS, 0.81 for BHB, and 0.80 for k-CN. Using a specific breed (Brown Swiss) as the training set for prediction in the Holstein population reduced the prediction accuracy by 10% for BCS, 7% for BHB, and 11% for k-CN. Notably, the combination of Holstein and Brown Swiss cows in the training set increased the predictive ability of the model by 6%, which was 0.66 for BCS, 0.85 for BHB, and 0.87 for k-CN. Using multiple specialized and dual-purpose animals in the training set outperforms the 10-fold_HO (standard) approach, with an increase in predictive ability of 8% for BCS, 7% for BHB, and 10% for k-CN. When the Multi-breed CV2 was implemented, no improvement was observed. Our findings suggest that FTIR prediction of different phenotypes in the Holstein breed can be improved by including different specialized and dual-purpose breeds in the training population. Our study also shows that predictive ability is enhanced when the size of the training population and the phenotypic variability are increased.

Prediction of meat quality traits in Nelore cattle by near-infrared reflectance spectroscopy.

The main definition for meat quality should include factors that affect consumer appreciation of the product. Physical laboratory analyses are necessary to identify factors that affect meat quality and specific equipment is used for this purpose, which is expensive and destructive, and the analyses are usually time consuming. An alternative method to performing several beef analyses is near-infrared reflectance spectroscopy (NIRS), which permits to reduce costs and to obtain faster, simpler, and nondestructive measurements. The objective of this study was to evaluate the feasibility of NIRS to predict shear force [Warner-Bratzler shear force (WBSF)], marbling, and color (*a = redness; b* = yellowness; and L* = lightness) in meat samples of uncastrated male Nelore cattle, that were approximately 2-yr-old. Samples of longissimus thoracis (n = 644) were collected and spectra were obtained prior to meat quality analysis. Multivariate calibration was performed by partial least squares regression. Several preprocessing techniques were evaluated alone and in combination: raw data, reduction of spectral range, multiplicative scatter correction, and 1st derivative. Accuracies of the calibration models were evaluated using the root mean square error of calibration (RMSEC), root mean square error of prediction (RMSEP), coefficient of determination in the calibration (R2C), and prediction (R2P) groups. Among the different preprocessing techniques, the reduction of spectral range provided the best prediction accuracy for all traits. The NIRS showed a better performance to predict WBSF (RMSEP = 1.42 kg, R2P = 0.40) and b* color (RMSEP = 1.21, R2P = 0.44), while its ability to accurately predict L* (RMSEP = 1.98, R2P = 0.16) and a* (RMSEP = 1.42, R2P = 0.17) was limited. NIRS was unsuitable to predict subjective meat quality traits such as marbling in Nelore cattle.

Circulating leptin and its muscle gene expression in Nellore cattle with divergent feed efficiency.

BACKGROUND: Leptin has a strong relation to important traits in animal production, such as carcass composition, feed intake, and reproduction. It is mainly produced by adipose cells and acts predominantly in the hypothalamus. In this study, circulating leptin and its gene expression in muscle were evaluated in two groups of young Nellore bulls with divergent feed efficiency. Individual dry matter intake (DMI) and average daily gain (ADG) of 98 Nellore bulls were evaluated in feedlot for 70 d to determinate the residual feed intake (RFI) and select 20 animals for the high feed efficient (LRFI) and 20 for the low feed efficient (HRFI) groups. Blood samples were collected on d 56 and at slaughter (80 d) to determine circulating plasma leptin. Samples of Longissimus dorsi were taken at slaughter for leptin gene expression levels. RESULTS: DMI and RFI were different between groups and LRFI animals showed less back fat and rump fat thickness, as well as less pelvic and kidney fat weight. Circulating leptin increased over time in all animals. Plasma leptin was greater in LRFI on 56 d and at slaughter (P = 0.0049). Gene expression of leptin were greater in LRFI animals (P = 0.0022) in accordance with the plasma levels. The animals of the LRFI group were leaner, ate less, and had more circulating leptin and its gene expression. CONCLUSION: These findings demonstrated that leptin plays its physiological role in young Nellore bulls, probably controlling food intake because feed efficient animals have more leptin and lower residual feed intake.