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.

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.

Publisher Correction: Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Whole-genome sequencing provides new insights into genetic mechanisms of tropical adaptation in Nellore (Bos primigenius indicus).

Most of the knowledge about genetic variants at the sequence level in cattle is for Bos primigenius taurus populations. Here, we presented a complete genomic characterization of 52 Nellore (Bos primigenius indicus) bulls, revealing specific zebu DNA variants with putative impact in tropical adaptation and productive traits. Single nucleotide polymorphisms (SNPs) and insertion/deletion (INDELs) mutations were identified using the newest bovine reference genome ARS_UCD1.2, and variant functional consequences were predicted using the Ensembl VEP software. A total of 35,753,707 SNPs and 4,492,636 INDELs were detected and annotated to their functional effects. We identified 400 genes that comprised both, a SNP and an INDEL, of high functional impact on proteins (i.e. variants that cause protein truncation, loss of function or triggering nonsense-mediated decay). Among these, we highlight the following genes: BoLA, associated with cattle immune response to infections and reproduction aspects; HSPA8, DNAJC27, and DNAJC28, involved with thermoregulatory protective mechanisms in mammals; and many olfactory signaling pathway related genes that are important genetic factors in the evolution of mammalian species. All these functional aspects are directly related to cattle adaptability to tropical environments.

Genome-wide prediction for complex traits under the presence of dominance effects in simulated populations using GBLUP and machine learning methods.

The aim of this study was to compare the predictive performance of the Genomic Best Linear Unbiased Predictor (GBLUP) and machine learning methods (Random Forest, RF; Support Vector Machine, SVM; Artificial Neural Network, ANN) in simulated populations presenting different levels of dominance effects. Simulated genome comprised 50k SNP and 300 QTL, both biallelic and randomly distributed across 29 autosomes. A total of six traits were simulated considering different values for the narrow and broad-sense heritability. In the purely additive scenario with low heritability (h2 = 0.10), the predictive ability obtained using GBLUP was slightly higher than the other methods whereas ANN provided the highest accuracies for scenarios with moderate heritability (h2 = 0.30). The accuracies of dominance deviations predictions varied from 0.180 to 0.350 in GBLUP extended for dominance effects (GBLUP-D), from 0.06 to 0.185 in RF and they were null using the ANN and SVM methods. Although RF has presented higher accuracies for total genetic effect predictions, the mean-squared error values in such a model were worse than those observed for GBLUP-D in scenarios with large additive and dominance variances. When applied to prescreen important regions, the RF approach detected QTL with high additive and/or dominance effects. Among machine learning methods, only the RF was capable to cover implicitly dominance effects without increasing the number of covariates in the model, resulting in higher accuracies for the total genetic and phenotypic values as the dominance ratio increases. Nevertheless, whether the interest is to infer directly on dominance effects, GBLUP-D could be a more suitable method.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle.

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

Genome-wide association studies for growth traits in buffaloes using the single step genomic BLUP.

Growth traits are important for the profitability of buffalo breeding systems, since in general, these animals are raised both for meat and milk. In this study, the single-step genomic BLUP method was employed to prospect the genomic regions' associated with weight at standard ages of 100, 210, 365, and 550 days in a buffalo population, aiming to identify genes with stronger expression for those characteristics. We found 6, 1, 2, and 5 SNPs significantly associated (p value < 10-5) with weight at 100, 210, 365, and 550 days of age, respectively, where those SNPs respectively explained 0.164, 0.040, 0.044, and 0.213% of the additive variance of each trait. SNP AX-85099682 (BBU24) was significant for weight at 100, 210, and 365 days, indicating the existence of a possible QTL affecting the initial growth rate of buffaloes. All told, eight genes (CBLB, TRNAG-UCC, GADD45B, LOC112583811, MGAT4C, KCNMA1, SLC5A2, and TGFB1I1) were identified as candidates for the growth traits of buffaloes. However, molecular and gene expression studies are necessary to validate these genes for subsequent use in programs for genetic improvement of the species.

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.