Assessing different metrics of pedigree and genomic inbreeding and inbreeding effect on growth, fertility, and feed efficiency traits in a closed-herd Nellore cattle population.

BACKGROUND: The selection of individuals based on their predicted breeding values and mating of related individuals can increase the proportion of identical-by-descent alleles. In this context, the objectives of this study were to estimate inbreeding coefficients based on alternative metrics and data sources such as pedigree (FPED), hybrid genomic relationship matrix H (FH), and ROH of different length (FROH); and calculate Pearson correlations between the different metrics in a closed Nellore cattle population selected for body weight adjusted to 378 days of age (W378). In addition to total FROH (all classes) coefficients were also estimated based on the size class of the ROH segments: FROH1 (1-2 Mb), FROH2 (2-4 Mb), FROH3 (4-8 Mb), FROH4 (8-16 Mb), and FROH5 (> 16 Mb), and for each chromosome (FROH_CHR). Furthermore, we assessed the effect of each inbreeding metric on birth weight (BW), body weights adjusted to 210 (W210) and W378, scrotal circumference (SC), and residual feed intake (RFI). We also evaluated the chromosome-specific effects of inbreeding on growth traits. RESULTS: The correlation between FPED and FROH was 0.60 while between FH and FROH and FH and FPED were 0.69 and 0.61, respectively. The annual rate of inbreeding was 0.16% for FPED, 0.02% for FH, and 0.16% for FROH. A 1% increase in FROH5 resulted in a reduction of up to -1.327 ± 0.495 kg in W210 and W378. Four inbreeding coefficients (FPED, FH, FROH2, and FROH5) had a significant effect on W378, with reductions of up to -3.810 ± 1.753 kg per 1% increase in FROH2. There was an unfavorable effect of FPED on RFI (0.01 ± 0.0002 kg dry matter/day) and of FROH on SC (-0.056 ± 0.022 cm). The FROH_CHR coefficients calculated for BTA3, BTA5, and BTA8 significantly affected the growth traits. CONCLUSIONS: Inbreeding depression was observed for all traits evaluated. However, these effects were greater for the criterion used for selection of the animals (i.e., W378). The increase in the genomic inbreeding was associated with a higher inbreeding depression on the traits evaluated when compared to pedigree-based inbreeding. Genomic information should be used as a tool during mating to optimize control of inbreeding and, consequently, minimize inbreeding depression in Nellore cattle.

Genome-wide association study between copy number variation and feeding behavior, feed efficiency, and growth traits in Nellore cattle.

BACKGROUND: Feeding costs represent the largest expenditures in beef production. Therefore, the animal efficiency in converting feed in high-quality protein for human consumption plays a major role in the environmental impact of the beef industry and in the beef producers' profitability. In this context, breeding animals for improved feed efficiency through genomic selection has been considered as a strategic practice in modern breeding programs around the world. Copy number variation (CNV) is a less-studied source of genetic variation that can contribute to phenotypic variability in complex traits. In this context, this study aimed to: (1) identify CNV and CNV regions (CNVRs) in the genome of Nellore cattle (Bos taurus indicus); (2) assess potential associations between the identified CNVR and weaning weight (W210), body weight measured at the time of selection (WSel), average daily gain (ADG), dry matter intake (DMI), residual feed intake (RFI), time spent at the feed bunk (TF), and frequency of visits to the feed bunk (FF); and, (3) perform functional enrichment analyses of the significant CNVR identified for each of the traits evaluated. RESULTS: A total of 3,161 CNVs and 561 CNVRs ranging from 4,973 bp to 3,215,394 bp were identified. The CNVRs covered up to 99,221,894 bp (3.99%) of the Nellore autosomal genome. Seventeen CNVR were significantly associated with dry matter intake and feeding frequency (number of daily visits to the feed bunk). The functional annotation of the associated CNVRs revealed important candidate genes related to metabolism that may be associated with the phenotypic expression of the evaluated traits. Furthermore, Gene Ontology (GO) analyses revealed 19 enrichment processes associated with FF. CONCLUSIONS: A total of 3,161 CNVs and 561 CNVRs were identified and characterized in a Nellore cattle population. Various CNVRs were significantly associated with DMI and FF, indicating that CNVs play an important role in key biological pathways and in the phenotypic expression of feeding behavior and growth traits in Nellore cattle.

Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits.

BACKGROUND: Feed efficiency (FE) related traits play a key role in the economy and sustainability of beef cattle production systems. The accurate knowledge of the physiologic background for FE-related traits can help the development of more efficient selection strategies for them. Hence, multi-trait weighted GWAS (MTwGWAS) and meta-analyze were used to find genomic regions associated with average daily gain (ADG), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), and residual feed intake (RFI). The FE-related traits and genomic information belong to two breeding programs that perform the FE test at different ages: post-weaning (1,024 animals IZ population) and post-yearling (918 animals for the QLT population). RESULTS: The meta-analyze MTwGWAS identified 14 genomic regions (-log10(p -value) > 5) regions mapped on BTA 1, 2, 3, 4, 7, 8, 11, 14, 15, 18, 21, and 29. These regions explained a large proportion of the total genetic variance for FE-related traits across-population ranging from 20% (FCR) to 36% (DMI) in the IZ population and from 22% (RFI) to 28% (ADG) in the QLT population. Relevant candidate genes within these regions (LIPE, LPL, IGF1R, IGF1, IGFBP5, IGF2, INS, INSR, LEPR, LEPROT, POMC, NPY, AGRP, TGFB1, GHSR, JAK1, LYN, MOS, PLAG1, CHCD7, LCAT, and PLA2G15) highlighted that the physiological mechanisms related to neuropeptides and the metabolic signals controlling the body's energy balance are responsible for leading to greater feed efficiency. Integrated meta-analysis results and functional pathway enrichment analysis highlighted the major effect of biological functions linked to energy, lipid metabolism, and hormone signaling that mediates the effects of peptide signals in the hypothalamus and whole-body energy homeostasis affecting the genetic control of FE-related traits in Nellore cattle. CONCLUSIONS: Genes and pathways associated with common signals for feed efficiency-related traits provide better knowledge about regions with biological relevance in physiological mechanisms associated with differences in energy metabolism and hypothalamus signaling. These pleiotropic regions would support the selection for feed efficiency-related traits, incorporating and pondering causal variations assigning prior weights in genomic selection approaches.

Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle.

BACKGROUND: This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection. RESULTS: Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (FST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4. CONCLUSIONS: The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.

Benchmarking machine learning and parametric methods for genomic prediction of feed efficiency-related traits in Nellore cattle.

Genomic selection (GS) offers a promising opportunity for selecting more efficient animals to use consumed energy for maintenance and growth functions, impacting profitability and environmental sustainability. Here, we compared the prediction accuracy of multi-layer neural network (MLNN) and support vector regression (SVR) against single-trait (STGBLUP), multi-trait genomic best linear unbiased prediction (MTGBLUP), and Bayesian regression (BayesA, BayesB, BayesC, BRR, and BLasso) for feed efficiency (FE) traits. FE-related traits were measured in 1156 Nellore cattle from an experimental breeding program genotyped for ~ 300 K markers after quality control. Prediction accuracy (Acc) was evaluated using a forward validation splitting the dataset based on birth year, considering the phenotypes adjusted for the fixed effects and covariates as pseudo-phenotypes. The MLNN and SVR approaches were trained by randomly splitting the training population into fivefold to select the best hyperparameters. The results show that the machine learning methods (MLNN and SVR) and MTGBLUP outperformed STGBLUP and the Bayesian regression approaches, increasing the Acc by approximately 8.9%, 14.6%, and 13.7% using MLNN, SVR, and MTGBLUP, respectively. Acc for SVR and MTGBLUP were slightly different, ranging from 0.62 to 0.69 and 0.62 to 0.68, respectively, with empirically unbiased for both models (0.97 and 1.09). Our results indicated that SVR and MTGBLUBP approaches were more accurate in predicting FE-related traits than Bayesian regression and STGBLUP and seemed competitive for GS of complex phenotypes with various degrees of inheritance.

Genotype by Prenatal Environment Interaction for Postnatal Growth of Nelore Beef Cattle Raised under Tropical Grazing Conditions.

The prenatal environment is recognized as crucial for the postnatal performance in cattle. In tropical regions, pregnant beef cows commonly experience nutritional restriction during the second half of the gestation period. Thus, the present study was designed to analyze the genotype by prenatal environment interaction (G × Epn) and to identify genomic regions associated with the level and response in growth and reproduction-related traits of beef cattle to changes in the prenatal environment. A reaction norm model was applied to data from two Nelore herds using the solutions of contemporary groups for birth weight as a descriptor variable of the gestational environment quality. A better gestational environment favored weights until weaning, scrotal circumference at yearling, and days to first calving of the offspring. The G × Epn was strong enough to result in heterogeneity of variance components and genetic parameters in addition to reranking of estimated breeding values and SNPs effects. Several genomic regions associated with the level of performance and specific responses of the animals to variations in the gestational environment were revealed, which harbor QTLs and can be exploited for selection purposes. Therefore, genetic evaluation models considering G × Epn and special management and nutrition care for pregnant cows are recommended.

Tannin-based product in feedlot diet as a strategy to reduce enteric methane emissions of Nellore cattle finished under tropical conditions.

A total of 120 Nellore bulls, [initial body weight (BW) = 307 ± 11.6 kg and 12 mo of age] were allocated into 12 collective pens (10 bulls per pen) in a commercial feedlot to evaluate the effects of a specific blend of tannin and saponins on enteric methane (CH4) emissions. The study was a completely randomized design, in which pens were considered the experimental units (N = 6 pens per treatment) and were randomly allocated into one of two treatments: 1) Control (CON), a basal diet with monensin supplementation (25 mg/kg dry matter [DM]; Rumensin, Elanco Animal Health, Greenfield, IN, USA), or 2) Control + a specific blend of tannin and saponins (TAN; 7 g/kg DM; composed of quebracho and chestnut tannin extracts along with carriers from cereals rich in saponins; SilvaFeed BX, Silvateam, San Michele Mondovi, CN, Italy). After the adaptation period (20 d), the experiment was divided into two phases: growing phase (21 to 53 d; total of 33 d) and fattening phase (54 to 139 d; total of 86 d). Enteric methane emissions were estimated using the sulfur hexafluoride (SF6) tracer gas technique. Interactions between treatment and period (growing vs. fattening) were detected for daily CH4 emissions, in which animals fed TAN reduced CH4 emissions by 17.3% during the fattening period compared to bulls fed CON (P = 0.05). In addition, bulls fed TAN had lower CH4 emissions expressed by dry matter intake (DMI) during the fattening period compared to bulls fed CON (P = 0.06). The findings presented herein indicate that a specific blend of tannin and saponins can be used as a strategy to reduce enteric CH4 emissions and its intensity of Nellore bulls finished in feedlot systems under tropical conditions.

Improving the accuracy of beef cattle methane inventories in Latin America and Caribbean countries.

On-farm methane (CH4) emissions need to be estimated accurately so that the mitigation effect of recommended practices can be accounted for. In the present study prediction equations for enteric CH4 have been developed in lieu of expensive animal measurement approaches. Our objectives were to: (1) compile a dataset from individual beef cattle data for the Latin America and Caribbean (LAC) region; (2) determine main predictors of CH4 emission variables; (3) develop and cross-validate prediction models according to dietary forage content (DFC); and (4) compare the predictive ability of these newly-developed models with extant equations reported in literature, including those currently used for CH4 inventories in LAC countries. After outlier's screening, 1100 beef cattle observations from 55 studies were kept in the final dataset (∼ 50 % of the original dataset). Mixed-effects models were fitted with a random effect of study. The whole dataset was split according to DFC into a subset for all-forage (DFC = 100 %), high-forage (94 % ≥ DFC ≥ 54 %), and low-forage (50 % ≥ DFC) diets. Feed intake and average daily gain (ADG) were the main predictors of CH4 emission (g d-1), whereas this was feeding level [dry matter intake (DMI) as % of body weight] for CH4 yield (g kg-1 DMI). The newly-developed models were more accurate than IPCC Tier 2 equations for all subsets. Simple and multiple regression models including ADG were accurate and a feasible option to predict CH4 emission when data on feed intake are not available. Methane yield was not well predicted by any extant equation in contrast to the newly-developed models. The present study delivered new models that may be alternatives for the IPCC Tier 2 equations to improve CH4 prediction for beef cattle in inventories of LAC countries based either on more or less readily available data.

Genetic parameters for stayability in Murrah buffaloes.

In order to contribute to the breeding programmes of Asian water buffalo, the aim of this study was to analyse the influence of genetic effects in the stayability of Murrah dairy buffaloes. The stayability trait (ST) was defined as the female's ability to stay in the herd for one (ST1), two (ST2), three (ST3), four (ST4), five (ST5) or six years (ST6) after the first calving. The same trait was also considered as continuous and was designated stayability in days up to one (STD1), two (STD2), three (STD3), four (STD4), five (STD5) or six years (STD6) after the first calving. Data from 1016 females reared in nine herds located in the State of São Paulo, Brazil, were analysed. Statistical models included the additive genetic effect of the animal and the fixed effects of the buffalo breeding herd, birth year and birth season. Additive effects for ST were estimated by approximate restricted maximum likelihood using a threshold model, while for STD, the additive effects were estimated by restricted maximum likelihood. Heritability estimates were lower for ST, except for ST1, (0.11+/-0.07, 0.17+/-0.06, 0.23+/-0.06, 0.16+/-0.08, 0.14+/-0.09 and 0.16+/-0.10 for ST1, ST2, ST3, ST4, ST5 and ST6, respectively) when compared with STD (0.05+/-0.06, 0.18+/-0.08, 0.40+/-0.10, 0.49+/-0.11, 0.41+/-0.11 and 0.30+/-0.13, for STD1, STD2, STD3, STD4, STD5 and STD6, respectively). Considering the values of heritability and owing to the serial nature of STD to a specific age, selection for STD3 should have a favourable influence on STD to other ages.

Association of ADIPOQ, OLR1 and PPARGC1A gene polymorphisms with growth and carcass traits in Nelore cattle.

In beef cattle farming, growth and carcass traits are important for genetic breeding programs. Molecular markers can be used to assist selection and increase genetic gain. The ADIPOQ, OLR1 and PPARGC1A genes are involved in lipid synthesis and fat accumulation in adipose tissue. The objective of this study was to identify polymorphisms in these genes and to assess the association with growth and carcass traits in Nelore cattle. A total of 639 animals were genotyped by PCR-RFLP for rs208549452, rs109019599 and rs109163366 in ADIPOQ, OLR1 and PPARGC1A gene, respectively. We analyzed the association of SNPs identified with birth weight, weaning weight, female yearling weight, female hip height, male yearling weight, male hip height, loin eye area, rump fat thickness, and backfat thickness. The OLR1 marker was associated with rump fat thickness and weaning weight (P < 0.05) and the PPARGC1 marker was associated with female yearling weight.