The objective of this study was to identify genomic regions that are associated with meat quality traits in the Nellore breed. Nellore steers were finished in feedlots and slaughtered at a commercial slaughterhouse. This analysis included 1,822 phenotypic records of tenderness and 1,873 marbling records. After quality control, 1,630 animals genotyped for tenderness, 1,633 animals genotyped for marbling, and 369,722 SNPs remained. The results are reported as the proportion of variance explained by windows of 150 adjacent SNPs. Only windows with largest effects were considered. The genomic regions were located on chromosomes 5, 15, 16 and 25 for marbling and on chromosomes 5, 7, 10, 14 and 21 for tenderness. These windows explained 3,89% and 3,80% of the additive genetic variance for marbling and tenderness, respectively. The genes associated with the traits are related to growth, muscle development and lipid metabolism. The study of these genes in Nellore cattle is the first step in the identification of causal mutations that will contribute to the genetic evaluation of the breed.

The objective of this study was to investigate the association of single nucleotide polymorphisms (SNPs) with birth weight, weight gain from birth to weaning and from weaning to yearling, yearling height and cow weight in Nelore cattle. Data from 5064 animals participating in the DeltaGen and PAINT breeding programs were used. The animals were genotyped with a panel of 777 962 SNPs (Illumina BovineHD BeadChip) and 412 993 SNPs remained after quality control analysis of the genomic data. A genome-wide association study was performed using a single-step methodology. The analyses were processed with the BLUPF90 family of programs. When applied to a genome-wide association studies, the single-step GBLUP methodology is an iterative process that estimates weights for the SNPs. The weights of SNPs were included in all analyses by iteratively applying the single-step GBLUP methodology and repeated twice so that the effect of the SNP and the effect of the animal were recalculated in order to increase the weight of SNPs with large effects and to reduce the weight of those with small effects. The genome-wide association results are reported based on the proportion of variance explained by windows of 50 adjacent SNPs. Considering the two iterations, only windows with an additive genetic variance >1.5% were presented in the results. Associations were observed with birth weight on BTA 14, with weight gain from birth to weaning on BTA 5 and 29, with weight gain from weaning to yearling on BTA 11, and with yearling height on BTA 8, showing the genes TMEM68 (transmembrane protein 8B) associated with birth weight and yearling height, XKR4 (XK, Kell blood group complex subunit-related family, member 4) associated with birth weight, NPR2 (natriuretic peptide receptor B) associated with yearling height, and REG3G (regenerating islet-derived 3-gamma) associated with weight gain from weaning to yearling. These genes play an important role in feed intake, weight gain and the regulation of skeletal growth.

BackgroundTraditional single nucleotide polymorphism (SNP) genome-wide association analysis (GWAA) can be inefficient because single SNPs provide limited genetic information about genomic regions. On the other hand, using haplotypes in the statistical analysis may increase the extent of linkage disequilibrium (LD) between haplotypes and causal variants and may also potentially capture epistastic interactions between variants within a haplotyped locus, providing an increase in the power and robustness of the association studies. We performed GWAA (413,355 SNP markers) using haplotypes based on variable-sized sliding windows and compared the results to a single-SNP GWAA using Warner-Bratzler shear force measured in the longissimus thorasis muscle of 3161 Nelore bulls to ascertain the optimal window size for identifying the genomic regions that influence meat tenderness.ResultsThe GWAA using single SNPs identified eight variants influencing meat tenderness on BTA 3, 4, 9, 10 and 11. However, thirty-three putative meat tenderness QTL were detected on BTA 1, 3, 4, 5, 8, 9, 10, 11, 15, 17, 18, 24, 25, 26 and 29 using variable-sized sliding haplotype windows. Analyses using sliding window haplotypes of 3, 5, 7, 9 and 11 SNPs identified 57, 61, 42, 39, and 21% of all thirty-three putative QTL regions, respectively; however, the analyses using the 3 and 5 SNP haplotypes, cumulatively detected 88% of the putative QTL. The genes associated with variation in meat tenderness participate in myogenesis, neurogenesis, lipid and fatty acid metabolism and skeletal muscle structure or composition processes.ConclusionsGWAA using haplotypes based on variable-sized sliding windows allowed the detection of more QTL than traditional single-SNP GWAA. Analyses using smaller haplotypes (3 and 5 SNPs) detected a higher proportion of the putative QTL.Electronic supplementary materialThe online version of this article (10.1186/s12863-019-0713-4) contains supplementary material, which is available to authorized users.

No abstract

The objective of this study was to determine whether visual scores used as selection criteria in Nellore breeding programs are effective indicators of carcass traits measured after slaughter. Additionally, this study evaluated the effect of different structures of the relationship matrix ( and ) on the estimation of genetic parameters and on the prediction accuracy of breeding values. There were 13,524 animals for visual scores of conformation (CS), finishing precocity (FP), and muscling (MS) and 1,753, 1,747, and 1,564 for LM area (LMA), backfat thickness (BF), and HCW, respectively. Of these, 1,566 animals were genotyped using a high-density panel containing 777,962 SNP. Six analyses were performed using multitrait animal models, each including the 3 visual scores and 1 carcass trait. For the visual scores, the model included direct additive genetic and residual random effects and the fixed effects of contemporary group (defined by year of birth, management group at yearling, and farm) and the linear effect of age of animal at yearling. The same model was used for the carcass traits, replacing the effect of age of animal at yearling with the linear effect of age of animal at slaughter. The variance and covariance components were estimated by the REML method in analyses using the numerator relationship matrix () or combining the genomic and the numerator relationship matrices (). The heritability estimates for the visual scores obtained with the 2 methods were similar and of moderate magnitude (0.23-0.34), indicating that these traits should response to direct selection. The heritabilities for LMA, BF, and HCW were 0.13, 0.07, and 0.17, respectively, using matrix and 0.29, 0.16, and 0.23, respectively, using matrix . The genetic correlations between the visual scores and carcass traits were positive, and higher correlations were generally obtained when matrix was used. Considering the difficulties and cost of measuring carcass traits postmortem, visual scores of CS, FP, and MS could be used as selection criteria to improve HCW, BF, and LMA. The use of genomic information permitted the detection of greater additive genetic variability for LMA and BF. For HCW, the high magnitude of the genetic correlations with visual scores was probably sufficient to recover genetic variability. The methods provided similar breeding value accuracies, especially for the visual scores.

Carcass traits measured after slaughter are economically relevant traits in beef cattle. In general, the slaughter house payment system is based on HCW. Ribeye area (REA) is associated with the amount of the meat in the carcass, and a minimum of backfat thickness (BFT) is necessary to protect the carcass during cooling. The aim of this study was to identify potential genomic regions harboring candidate genes affecting those traits in Nellore cattle. The data set used in the present study consisted of 1,756 Nellore males with phenotype records. A subset of 1,604 animals had both genotypic and phenotypic information. Genotypes were generated based on a panel with 777,962 SNPs from the Illumina Bovine HD chip. The SNP effects were calculated based on the genomic breeding values obtained by using the single-step GBLUP approach and a genomic matrix re-weighting procedure. The proportion of the variance explained by moving windows of 100 consecutive SNPs was used to assess potential genomic regions harboring genes with major effects on each trait. The top 10 non-overlapping SNP-windows explained 8.72%, 11.38%, and 9.31% of the genetic variance for REA, BFT, and HCW, respectively. These windows are located on chromosomes 5, 7, 8, 10, 12, 20, and 29 for REA; chromosomes 6, 8, 10, 13, 16, 17, 18, and 24 for BFT; and chromosomes 4, 6, 7, 8, 14, 16, 17, and 21 for HCW. For REA, there were identified genes ( and ) involved in the cell cycle biological process which affects many aspects of animal growth and development. The and genes, both from AA transporter family, was also associated with REA. The AA transporters are essential for cell growth and proliferation, acting as carriers of tissue nutrient supplies. Various genes identified for BFT (, , , , , and ) have been associated with lipid metabolism in different mammal species. One of the most promising genes identified for HCW was the . There is evidence, in the literature, that this gene is located in putative QTL affecting carcass weight in beef cattle. Our results showed several genomic regions containing plausible candidate genes that may be associated with carcass traits in Nellore cattle. Besides contributing to a better understanding of the genetic control of carcass traits, the identified genes can also be helpful for further functional genomic studies.

Multitrait meta-analyses are a strategy to produce more accurate genome-wide association studies, especially for complex phenotypes. We carried out a meta-analysis study for traits related to sexual precocity in tropical beef cattle (Nellore and Brahman) aiming to identify important genomic regions affecting these traits. The traits included in the analyses were age at first calving (AFC), early pregnancy (EP), age at first corpus luteum (AGECL), first postpartum anoestrus interval (PPAI), and scrotal circumference (SC). The traits AFC, EP, and SCN were measured in Nellore cattle, while AGECL, PPAI, and SCB were measured in Brahman cattle. Meta-analysis resulted in 108 significant single-nucleotide polymorphisms (SNPs), at an empirical threshold P-value of 1.39 × 10-5 (false discovery rate [FDR] < 0.05). Within 0.5 Mb of the significant SNP, candidate genes were annotated and analyzed for functional enrichment. Most of the closest genes to the SNP with higher significance in each chromosome have been associated with important roles in reproductive function. They are TSC22D2, KLF7, ARHGAP29, 7SK, MAP3K5, TLE3, WDR5, TAF3, TMEM68, PPP1R15B, NR2F2, GALR1, SUFU, and KCNU1. We did not observe any significant SNP in BTA5, BTA12, BTA17, BTA18, BTA19, BTA20, BTA22, BTA23, BTA25, and BTA28. Although the majority of significant SNPs are in BTA14, it was identified significant associations in multiple chromosomes (19 out of 29 autosomes), which is consistent with the postulation that reproductive traits are complex polygenic phenotypes. Five proposed association regions harbor the majority of the significant SNP (76%) and were distributed over four chromosomes (P < 1.39 × 10-5, FDR < 0.05): BTA2 (5.55%) from 95 to 96 Mb, BTA4 (5.55%) from 94.1 to 94.8 Mb, BTA14 (59.26%) from 24 to 25 Mb and 29 to 30 Mb, and BTA21 (5.55%) from 6.7 Mb to 11.4 Mb. These regions harbored key genes related to reproductive function. Moreover, these genes were enriched for functional groups associated with immune response, maternal-fetal tolerance, pregnancy maintenance, embryo development, fertility, and response to stress. Further studies including other breeds and precocity traits could confirm the importance of these regions and identify new candidate regions for sexual precocity in beef 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.