Economic values of reproductive, growth, feed efficiency and carcass traits in Nellore cattle.

A bioeconomic model was developed to calculate the economic value (ev) of reproductive and growth performance, feed efficiency and carcass traits of a seedstock Nellore herd. Data from a full-cycle cattle operation (1,436 dams) located in the Brazilian Cerrado were assessed. The ev was calculated by the difference in profit before and after one-unit improvement in the trait, with others remaining unchanged. The ev was standardized by the phenotypic standard deviation of each trait. Preweaning average daily gain (ADG) was the most economically important trait evaluated (R$ 58.04/animal/year), followed by age at first calving (R$ 44.35), postweaning ADG (R$ 31.43), weight at 450 days (R$ 25.36), accumulated productivity (R$ 21.43), ribeye area (R$ 21.35), calving interval (R$ 19.97), feed efficiency (R$ 15.24), carcass dressing per cent (R$ 8.27), weight at 120 days (R$ 6.22), weight at 365 days (R$ 6.06), weight at weaning (210 days, R$ 5.82), stayability (R$ 5.70) and the probability of early calving (R$ 0.32). The effects of all traits on profits are evidence that their selection may result in the economic and genetic progress of the herd if there is genetic variability.

Genomic predictions combining SNP markers and copy number variations in Nellore cattle.

BACKGROUND: Due to the advancement in high throughput technology, single nucleotide polymorphism (SNP) is routinely being incorporated along with phenotypic information into genetic evaluation. However, this approach often cannot achieve high accuracy for some complex traits. It is possible that SNP markers are not sufficient to predict these traits due to the missing heritability caused by other genetic variations such as microsatellite and copy number variation (CNV), which have been shown to affect disease and complex traits in humans and other species. RESULTS: In this study, CNVs were included in a SNP based genomic selection framework. A Nellore cattle dataset consisting of 2230 animals genotyped on BovineHD SNP array was used, and 9 weight and carcass traits were analyzed. A total of six models were implemented and compared based on their prediction accuracy. For comparison, three models including only SNPs were implemented: 1) BayesA model, 2) Bayesian mixture model (BayesB), and 3) a GBLUP model without polygenic effects. The other three models incorporating both SNP and CNV included 4) a Bayesian model similar to BayesA (BayesA+CNV), 5) a Bayesian mixture model (BayesB+CNV), and 6) GBLUP with CNVs modeled as a covariable (GBLUP+CNV). Prediction accuracies were assessed based on Pearson's correlation between de-regressed EBVs (dEBVs) and direct genomic values (DGVs) in the validation dataset. For BayesA, BayesB and GBLUP, accuracy ranged from 0.12 to 0.62 across the nine traits. A minimal increase in prediction accuracy for some traits was noticed when including CNVs in the model (BayesA+CNV, BayesB+CNV, GBLUP+CNV). CONCLUSIONS: This study presents the first genomic prediction study integrating CNVs and SNPs in livestock. Combining CNV and SNP marker information proved to be beneficial for genomic prediction of some traits in Nellore cattle.

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.

In-depth exploration of the high and normal pH beef proteome: First insights emphasizing the dynamic protein changes in Longissimus thoracis muscle from pasture-finished Nellore bulls over different postmortem times.

This study aimed to evaluate for the first time the temporal dynamic changes in early postmortem proteome of normal and high ultimate pH (pHu) beef samples from the same cattle using a shotgun proteomics approach. Ten selected carcasses classified as normal (pHu < 5.8; n = 5) or high (pHu ≥ 6.2; n = 5) pHu beef from pasture-finished Nellore (Bos taurus indicus) bulls were sampled from Longissimus thoracis muscle at 30 min, 9 h and 44 h postmortem for proteome comparison. The temporal proteomics profiling quantified 863 proteins, from which 251 were differentially abundant (DAPs) between high and normal pHu at 30 min (n = 33), 9 h (n = 181) and 44 h (n = 37). Among the myriad interconnected pathways regulating pH decline during postmortem metabolism, this study revealed the pivotal role of energy metabolism, cellular response to stress, oxidoreductase activity and muscle system process pathways throughout the early postmortem. Twenty-three proteins overlap among postmortem times and may be suggested as candidate biomarkers to the dark-cutting condition development. The study further evidenced for the first time the central role of ribosomal proteins and histones in the first minutes after animal bleeding. Moreover, this study revealed the disparity in the mechanisms underpinning the development of dark-cutting beef condition among postmortem times, emphasizing multiple dynamic changes in the muscle proteome. Therefore, this study revealed important insights regarding the temporal dynamic changes that occur in early postmortem of high and normal muscle pHu beef, proposing specific pathways to determine the biological mechanisms behind dark-cutting determination.

The Quality of Salted Sun-Dried Meat from Young Nellore Bulls Fed Diets with Lauric Acid.

This study aims to evaluate the quality of salted sun-dried meat from young bulls (Nellore cattle) fed with a diet containing 0.0, 0.5, 1.0 and 1.5% of lauric acid in the total dry matter (DM). Thirty-two Nellore bulls with initial body weight of 368 ± 32 kg were used. A linear decrease (p < 0.05) in pH and protein content of the salted sun-dried meat was observed with the inclusion of lauric acid. The moisture, ash, lipid, collagen content, water-holding capacity, cooking loss, color indexes (L*, a*, b*, C*), and shear force were not affected. Lipid oxidation at 7 days of storage increased linearly in the salted sun-dried meat. Most of the fatty acid composition of the salted sun-dried meat from the semimembranosus muscle of young bulls was not influenced (p > 0.05) by the lauric acid inclusion in the bulls' diet. However, there was a linear increase (p < 0.05) in the SFA lauric acid (C12:0), PUFAn-3 EPA (C20:5n − 3) and DHA (C22:6n − 3), and a quadratic increase in the PUFAn-6 arachidonic (C20:4n − 6) due to lauric acid addition from palm kernel oil in the diet. There was a liner increase (p < 0.05) in the total ∑PUFA, ∑n − 6, ∑n − 3 contents of salted sun-dried meat from the semimembranosus muscle of young bulls and the h:H health index of the level of lauric acid inclusion in bull's diet. In contrast, the thrombogenicity health index (TI) and ∑n − 6:∑n − 3 ratio content in salted sun-dried meat from the semimembranosus muscle of young bulls presented a linear decrease (p < 0.05) due to lauric acid addition in the bulls' diet. Lauric acid (C12:0) inclusion up to 1.5% in the diet of young Nellore bull improved the fatty acid composition of the salted sun-dried meat, increasing EPA, DHA, n − 6 and n − 3, TI, and h:H indexes, which are associated with a better lipid quality of meat products, and further improves tenderness at the highest concentration.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness.

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR < 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

An assessment of genomic connectedness measures in Nellore cattle.

An important criterion to consider in genetic evaluations is the extent of genetic connectedness across management units (MU), especially if they differ in their genetic mean. Reliable comparisons of genetic values across MU depend on the degree of connectedness: the higher the connectedness, the more reliable the comparison. Traditionally, genetic connectedness was calculated through pedigree-based methods; however, in the era of genomic selection, this can be better estimated utilizing new approaches based on genomics. Most procedures consider only additive genetic effects, which may not accurately reflect the underlying gene action of the evaluated trait, and little is known about the impact of non-additive gene action on connectedness measures. The objective of this study was to investigate the extent of genomic connectedness measures, for the first time, in Brazilian field data by applying additive and non-additive relationship matrices using a fatty acid profile data set from seven farms located in the three regions of Brazil, which are part of the three breeding programs. Myristic acid (C14:0) was used due to its importance for human health and reported presence of non-additive gene action. The pedigree included 427,740 animals and 925 of them were genotyped using the Bovine high-density genotyping chip. Six relationship matrices were constructed, parametrically and non-parametrically capturing additive and non-additive genetic effects from both pedigree and genomic data. We assessed genome-based connectedness across MU using the prediction error variance of difference (PEVD) and the coefficient of determination (CD). PEVD values ranged from 0.540 to 1.707, and CD from 0.146 to 0.456. Genomic information consistently enhanced the measures of connectedness compared to the numerator relationship matrix by at least 63%. Combining additive and non-additive genomic kernel relationship matrices or a non-parametric relationship matrix increased the capture of connectedness. Overall, the Gaussian kernel yielded the largest measure of connectedness. Our findings showed that connectedness metrics can be extended to incorporate genomic information and non-additive genetic variation using field data. We propose that different genomic relationship matrices can be designed to capture additive and non-additive genetic effects, increase the measures of connectedness, and to more accurately estimate the true state of connectedness in herds.

Use of nonlinear mixed models for describing testicular volume growth curve in Nellore bulls.

This study aimed to describe longitudinal testicular volume (TV) data of Nellore bulls by using nonlinear mixed models. Dataset consisted of 2,294 TV measurements from 505 bulls with ages ranging from 563 to 4,307 days. Nine nonlinear models were evaluated: Brody, Gompertz, Hill, Logistic I and II, Meloun, Michaelis-Menten, Mitscherlich and von Bertalanffy. Goodness of fit was evaluated by Akaike's information criterion (AIC), Bayesian information criterion (BIC), adjusted R2, percentage of convergence, error sum of squares (ESS), mean absolute deviation (MAD), and average prediction error (APE). These criteria were used to select the best model, then the absolute growth rate (AGR) for TV was estimated by the first derivate of the adjusted model related to time (∂Y/∂t). The values of adjusted R2, ESS, MAD and APE were similar among models. Percentage of convergence was higher for the Logistic I (76.8%), Logistic II (75.5%) and Mitscherlich (78.6%) models, but Logistic I and II showed the lowest values of AIC and BIC, indicating a better fit, so the Logistic I model was chosen for subsequent analyses. The TV growth occurred at a high rate until the inflection point, which was estimated at approximately 22 months of age; it stabilized and reached a plateau at approximately 2,500 days of age. This may suggest that TV is more related to sexual maturity than precocity. Additionally, the Logistic I model was used to estimate the growth curve and the AGR of the testicular length and width. As a result, testicular length increased at a higher rate than testicular width until approximately 1,600 days of age, indicating that the testes become longer with increasing age.

DNA methylation profile at a satellite region is associated with aberrant placentation in cloned calves.

INTRODUCTION: Cloning via somatic cell nuclear transfer (SCNT) has been associated with a variety of pathologies, primarily in the placenta, and these alterations may be associated with aberrant epigenetic reprogramming of the donor cell genome. We tested the hypothesis that DNA methylation patterns are not appropriately established after nuclear transfer and that those altered patterns are associated with specific aberrant phenotypes. METHODS: We compared global and specific placental DNA methylation patterns between aberrant and healthy SCNT-produced calves. Foetal cotyledon samples of ten SCNT pregnancies were collected. Global DNA methylation and hydroxymethylation levels were measured using an ELISA-based assay and specific DNA methylation of satellite I, and α-satellite repeat elements were measured using bisulfite PCR. RESULTS: Our analysis revealed that the SCNT-produced calves, which showed aberrant phenotypes, exhibited a reduced methylation pattern of the satellite I region compared to that of healthy calves. In contrast, global methylation and hydroxymethylation analyses showed higher levels for both cytosine modifications in SCNT-produced female calves with aberrant phenotypes. The satellite I region showed most of the sequences to be hypermethylated in live cloned calves compared with those in deceased calves. DISCUSSION: Our results suggest that this satellite I region could be used as an epigenetic biomarker for predicting offspring viability. Studies evaluating DNA methylation patterns of this satellite region in the donor cell genome or embryo biopsies could shed light on how to improve the efficiency of SCNT cloning.

Copy number variations and genome-wide associations reveal putative genes and metabolic pathways involved with the feed conversion ratio in beef cattle.

The use of genome-wide association results combined with other genomic approaches may uncover genes and metabolic pathways related to complex traits. In this study, the phenotypic and genotypic data of 1475 Nellore (Bos indicus) cattle and 941,033 single nucleotide polymorphisms (SNPs) were used for genome-wide association study (GWAS) and copy number variations (CNVs) analysis in order to identify candidate genes and putative pathways involved with the feed conversion ratio (FCR). The GWAS was based on the Bayes B approach analyzing genomic windows with multiple regression models to estimate the proportion of genetic variance explained by each window. The CNVs were detected with PennCNV software using the log R ratio and B allele frequency data. CNV regions (CNVRs) were identified with CNVRuler and a linear regression was used to associate CNVRs and the FCR. Functional annotation of associated genomic regions was performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) and the metabolic pathways were obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG). We showed five genomic windows distributed over chromosomes 4, 6, 7, 8, and 24 that explain 12 % of the total genetic variance for FCR, and detected 12 CNVRs (chromosomes 1, 5, 7, 10, and 12) significantly associated [false discovery rate (FDR) < 0.05] with the FCR. Significant genomic regions (GWAS and CNV) harbor candidate genes involved in pathways related to energetic, lipid, and protein metabolism. The metabolic pathways found in this study are related to processes directly connected to feed efficiency in beef cattle. It was observed that, even though different genomic regions and genes were found between the two approaches (GWAS and CNV), the metabolic processes covered were related to each other. Therefore, a combination of the approaches complement each other and lead to a better understanding of the FCR.

Body condition score of Nellore beef cows: a heritable measure to improve the selection of reproductive and maternal traits.

Despite the economic importance of beef cattle production in Brazil, female reproductive performance, which is strongly associated with production efficiency, is not included in the selection index of most breeding programmes due to low heritability and difficulty in measure. The body condition score (BCS) could be used as an indicator of these traits. However, so far little is known about the feasibility of using BCS as a selection tool for reproductive performance in beef cattle. In this study, we investigated the sources of variation in the BCS of Nellore beef cows, quantified its association with reproductive and maternal traits and estimated its heritability. BCS was analysed using a logistic model that included the following effects: contemporary group at weaning, cow weight and hip height, calving order, reconception together with the weight and scores of conformation and early finishing assigned to calves at weaning. In the genetic analysis, variance components of BCS were estimated through Bayesian inference by fitting an animal model that also included the aforementioned effects. The results showed that BCS was significantly associated with all of the reproductive and maternal variables analysed. The estimated posterior mean of heritability of BCS was 0.24 (highest posterior density interval at 95%: 0.093 to 0.385), indicating an involvement of additive gene action in its determination. The present findings show that BCS can be used as a selection criterion for Nellore females.