Accuracy of genomic breeding values and predictive ability for postweaning liveweight and age at first calving in a Nellore cattle population with missing sire information.

The multiple sire system (MSS) is a common mating scheme in extensive beef production systems. However, MSS does not allow paternity identification and lead to inaccurate genetic predictions. The objective of this study was to investigate the implementation of single-step genomic BLUP (ssGBLUP) in different scenarios of uncertain paternity in the evaluation for 450-day adjusted liveweight (W450) and age at first calving (AFC) in a Nellore cattle population. To estimate the variance components using BLUP and ssGBLUP, the relationship matrix (A) with different proportions of animals with missing sires (MS) (scenarios 0, 25, 50, 75, and 100% of MS) was created. The genotyped animals with MS were randomly chosen, and ten replicates were performed for each scenario and trait. Five groups of animals were evaluated in each scenario: PHE, all animals with phenotypic records in the population; SIR, proven sires; GEN, genotyped animals; YNG, young animals without phenotypes and progeny; and YNGEN, young genotyped animals. The additive genetic variance decreased for both traits as the proportion of MS increased in the population when using the regular REML. When using the ssGBLUP, accuracies ranged from 0.13 to 0.47 for W450 and from 0.10 to 0.25 for AFC. For both traits, the prediction ability of the direct genomic value (DGV) decreased as the percentage of MS increased. These results emphasize that indirect prediction via DGV of young animals is more accurate when the SNP effects are derived from ssGBLUP with a reference population with known sires. The ssGBLUP could be applied in situations of uncertain paternity, especially when selecting young animals. This methodology is shown to be accurate, mainly in scenarios with a high percentage of MS.

Genomic selection for meat quality traits in Nelore cattle.

The objective of this study was to present heritability estimates and accuracy of genomic prediction using different methods for meat quality traits in Nelore cattle. Approximately 5000 animals with phenotypes and genotypes of 412,000 SNPs, were divided into two groups: (1) training population: animals born from 2008 to 2013 and (2) validation population: animals born in 2014. A single-trait animal model was used to estimate heritability and to adjust the phenotype. The methods of GBLUP, Improved Bayesian Lasso and Bayes Cπ were performed to estimate the SNP effects. Accuracy of genomic prediction was calculated using Pearson's correlations between direct genomic values and adjusted phenotypes, divided by the square root of heritability of each trait (0.03-0.19). The accuracies varied from 0.23 to 0.73, with the lowest accuracies estimated for traits associated with fat content and the greatest accuracies observed for traits of meat color and tenderness. There were small differences in genomic prediction accuracy between methods.

Association study between copy number variation and beef fatty acid profile of Nellore cattle.

The aim of this study was to analyze the association between the copy number variation regions (CNVRs) and fatty acid profile phenotypes for saturated (SFA), monosaturated (MUFA), polyunsaturated (PUFA), ω6 and ω3 fatty acids, PUFA/SFA and ω6/ω3 ratios, as well as for their sums, in Nellore cattle (Bos primigenius indicus). A total of 963 males were finished in feedlot and slaughtered with approximately 2 years of age. Animals were genotyped with the BovineHD BeadChip (Illumina Inc., San Diego, CA, USA). The copy number variation (CNV) detection was performed using the PennCNV algorithm. Log R ratio (LRR) and allele B frequency (BAF) were used to estimate the CNVs. The association analyses were done using the CNVRuler software and applying a logistic regression model. The phenotype was adjusted using a linear model considering the fixed effects of contemporary group and the animal age at slaughter. The fatty acid profile was analyzed on samples of longissimus thoracis muscle using gas chromatography with a 100-m capillary column. For the association analysis, the adjusted phenotypic values were considered for the traits, while the data was adjusted for the effects of the farm and year of birth, management groups at birth, weaning, and superannuation. A total of 186 CNVRs were significant for SFA (43), MUFA (42), PUFA (66), and omega fatty acid (35) groups, totaling 278 known genes. On the basis of the results, several genes were associated with several fatty acids of different saturations. Olfactory receptor genes were associated with C12:0, C14:0, and C18:0 fatty acids. The SAMD8 and BSCL2 genes, both related to lipid metabolic process, were associated with C12:0. The RAPGEF6 gene was found to be associated with C18:2 cis-9 cis-12 n-6, and its function is related to regulation of GTPase activity. Among the results, we highlighted the olfactory receptor activity (GO:0004984), G-protein-coupled receptor activity (GO:0004930), potassium:proton antiporter activity (GO:0015386), sodium:proton antiporter activity (GO:0015385), and odorant-binding (GO:0005549) molecular functions. A large number of genes associated with fatty acid profile within the CNVRs were identified in this study. These findings must contribute to better elucidate the genetic mechanism underlying the fatty acid profile of intramuscular fat in Nellore cattle.

Genomic prediction ability for beef fatty acid profile in Nelore cattle using different pseudo-phenotypes.

The aim of the present study was to compare the predictive ability of SNP-BLUP model using different pseudo-phenotypes such as phenotype adjusted for fixed effects, estimated breeding value, and genomic estimated breeding value, using simulated and real data for beef FA profile of Nelore cattle finished in feedlot. A pedigree with phenotypes and genotypes of 10,000 animals were simulated, considering 50% of multiple sires in the pedigree. Regarding to phenotypes, two traits were simulated, one with high heritability (0.58), another with low heritability (0.13). Ten replicates were performed for each trait and results were averaged among replicates. A historical population was created from generation zero to 2020, with a constant size of 2000 animals (from generation zero to 1000) to produce different levels of linkage disequilibrium (LD). Therefore, there was a gradual reduction in the number of animals (from 2000 to 600), producing a "bottleneck effect" and consequently, genetic drift and LD starting in the generation 1001 to 2020. A total of 335,000 markers (with MAF greater or equal to 0.02) and 1000 QTL were randomly selected from the last generation of the historical population to generate genotypic data for the test population. The phenotypes were computed as the sum of the QTL effects and an error term sampled from a normal distribution with zero mean and variance equal to 0.88. For simulated data, 4000 animals of the generations 7, 8, and 9 (with genotype and phenotype) were used as training population, and 1000 animals of the last generation (10) were used as validation population. A total of 937 Nelore bulls with phenotype for fatty acid profiles (Sum of saturated, monounsaturated, omega 3, omega 6, ratio of polyunsaturated and saturated and polyunsaturated fatty acid profile) were genotyped using the Illumina BovineHD BeadChip (Illumina, San Diego, CA) with 777,962 SNP. To compare the accuracy and bias of direct genomic value (DGV) for different pseudo-phenotypes, the correlation between true breeding value (TBV) or DGV with pseudo-phenotypes and linear regression coefficient of the pseudo-phenotypes on TBV for simulated data or DGV for real data, respectively. For simulated data, the correlations between DGV and TBV for high heritability traits were higher than obtained with low heritability traits. For simulated and real data, the prediction ability was higher for GEBV than for Yc and EBV. For simulated data, the regression coefficient estimates (b(Yc,DGV)), were on average lower than 1 for high and low heritability traits, being inflated. The results were more biased for Yc and EBV than for GEBV. For real data, the GEBV displayed less biased results compared to Yc and EBV for SFA, MUFA, n-3, n-6, and PUFA/SFA. Despite the less biased results for PUFA using the EBV as pseudo-phenotype, the b(Yi,DGV estimates obtained for the different pseudo-phenotypes (Yc, EBV and GEBV) were very close. Genomic information can assist in improving beef fatty acid profile in Zebu cattle, since the use of genomic information yielded genomic values for fatty acid profile with accuracies ranging from low to moderate. Considering both simulated and real data, the ssGBLUP model is an appropriate alternative to obtain more reliable and less biased GEBVs as pseudo-phenotype in situations of missing pedigree, due to high proportion of multiple sires, being more adequate than EBV and Yc to predict direct genomic value for beef fatty acid profile.

Application of single step genomic BLUP under different uncertain paternity scenarios using simulated data.

The objective of this study was to investigate the application of BLUP and single step genomic BLUP (ssGBLUP) models in different scenarios of paternity uncertainty with different strategies of scaling the G matrix to match the A22 matrix, using simulated data for beef cattle. Genotypes, pedigree, and phenotypes for age at first calving (AFC) and weight at 550 days (W550) were simulated using heritabilities based on real data (0.12 for AFC and 0.34 for W550). Paternity uncertainty scenarios using 0, 25, 50, 75, and 100% of multiple sires (MS) were studied. The simulated genome had a total length of 2,333 cM, containing 735,293 biallelic markers and 7,000 QTLs randomly distributed over the 29 BTA. It was assumed that QTLs explained 100% of the genetic variance. For QTL, the amount of alleles per loci randomly ranged from two to four. The BLUP model that considers phenotypic and pedigree data, and the ssGBLUP model that combines phenotypic, pedigree and genomic information were used for genetic evaluations. Four ways of scaling the mean of the genomic matrix (G) to match to the mean of the pedigree relationship matrix among genotyped animals (A22) were tested. Accuracy, bias, and inflation were investigated for five groups of animals: ALL = all animals; BULL = only bulls; GEN = genotyped animals; FEM = females; and YOUNG = young males. With the BLUP model, the accuracies of genetic evaluations decreased for both traits as the proportion of unknown sires in the population increased. The EBV accuracy reduction was higher for GEN and YOUNG groups. By analyzing the scenarios for YOUNG (from 0 to 100% of MS), the decrease was 87.8 and 86% for AFC and W550, respectively. When applying the ssGBLUP model, the accuracies of genetic evaluation also decreased as the MS in the pedigree for both traits increased. However, the accuracy reduction was less than those observed for BLUP model. Using the same comparison (scenario 0 to 100% of MS), the accuracies reductions were 38 and 44.6% for AFC and W550, respectively. There were no differences between the strategies for scaling the G matrix for ALL, BULL, and FEM groups under the different scenarios with missing pedigree. These results pointed out that the uninformative part of the A22 matrix and genotyped animals with paternity uncertainty did not influence the scaling of G matrix. On the basis of the results, it is important to have a G matrix in the same scale of the A22 matrix, especially for the evaluation of young animals in situations with missing pedigree information. In these situations, the ssGBLUP model is an appropriate alternative to obtain a more reliable and less biased estimate of breeding values, especially for young animals with few or no phenotypic records. For accurate and unbiased genomic predictions with ssGBLUP, it is necessary to assure that the G matrix is compatible with the A22 matrix, even in situations with paternity uncertainty.

Genetic correlation estimates between beef fatty acid profile with meat and carcass traits in Nellore cattle finished in feedlot.

The objective of this study was to estimate the genetic-quantitative relationships between the beef fatty acid profile with the carcass and meat traits of Nellore cattle. A total of 1826 bulls finished in feedlot conditions and slaughtered at 24 months of age on average were used. The following carcass and meat traits were analysed: subcutaneous fat thickness (BF), shear force (SF) and total intramuscular fat (IMF). The fatty acid (FA) profile of the Longissimus thoracis samples was determined. Twenty-five FAs (18 individuals and seven groups of FAs) were selected due to their importance for human health. The animals were genotyped with the BovineHD BeadChip and, after quality control for single nucleotide polymorphisms (SNPs), only 470,007 SNPs from 1556 samples remained. The model included the random genetic additive direct effect, the fixed effect of the contemporary group and the animal's slaughter age as a covariable. The (co)variances and genetic parameters were estimated using the REML method, considering an animal model (single-step GBLUP). A total of 25 multi-trait analyses, with four traits, were performed considering SF, BF and IMF plus each individual FA. The heritability estimates for individual saturated fatty acids (SFA) varied from 0.06 to 0.65, for monounsaturated fatty acids (MUFA) it varied from 0.02 to 0.14 and for polyunsaturated fatty acids (PUFA) it ranged from 0.05 to 0.68. The heritability estimates for Omega 3, Omega 6, SFA, MUFA and PUFA sum were low to moderate, varying from 0.09 to 0.20. The carcass and meat traits, SF (0.06) and IMF (0.07), had low heritability estimates, while BF (0.17) was moderate. The genetic correlation estimates between SFA sum, MUFA sum and PUFA sum with BF were 0.04, 0.64 and -0.41, respectively. The genetic correlation estimates between SFA sum, MUFA sum and PUFA sum with SF were 0.29, -0.06 and -0.04, respectively. The genetic correlation estimates between SFA sum, MUFA sum and PUFA sum with IMF were 0.24, 0.90 and -0.67, respectively. The selection to improve meat tenderness in Nellore cattle should not change the fatty acid composition in beef, so it is possible to improve this attribute without affecting the nutritional beef quality in zebu breeds. However, selection for increased deposition of subcutaneous fat thickness and especially the percentage of intramuscular fat should lead to changes in the fat composition, highlighting a genetic antagonism between meat nutritional value and acceptability by the consumer.

Genomic prediction for beef fatty acid profile in Nellore cattle.

The objective of this study was to compare SNP-BLUP, BayesCπ, BayesC and Bayesian Lasso methodologies to predict the direct genomic value for saturated, monounsaturated, and polyunsaturated fatty acid profile, omega 3 and 6 in the Longissimus thoracis muscle of Nellore cattle finished in feedlot. A total of 963 Nellore bulls with phenotype for fatty acid profiles, were genotyped using the Illumina BovineHD BeadChip (Illumina, San Diego, CA) with 777,962 SNP. The predictive ability was evaluated using cross validation. To compare the methodologies, the correlation between DGV and pseudo-phenotypes was calculated. The accuracy varied from -0.40 to 0.62. Our results indicate that none of the methods excelled in terms of accuracy, however, the SNP-BLUP method allows obtaining less biased genomic evaluations, thereby; this method is more feasible when taking into account the analyses' operating cost. Despite the lowest bias observed for EBV, the adjusted phenotype is the preferred pseudophenotype considering the genomic prediction accuracies regarding the context of the present study.