B. P. Sollero

Genomic prediction for tick resistance in Braford and Hereford cattle.

One of the main animal health problems in tropical and subtropical cattle production is the bovine tick, which causes decreased performance, hide devaluation, increased production costs with acaricide treatments, and transmission of infectious diseases. This study investigated the utility of genomic prediction as a tool to select Braford (BO) and Hereford (HH) cattle resistant to ticks. The accuracy and bias of different methods for direct and blended genomic prediction was assessed using 10,673 tick counts obtained from 3,435 BO and 928 HH cattle belonging to the Delta G Connection breeding program. A subset of 2,803 BO and 652 HH samples were genotyped and 41,045 markers remained after quality control. Log transformed records were adjusted by a pedigree repeatability model to estimate variance components, genetic parameters, and breeding values (EBV) and subsequently used to obtain deregressed EBV. Estimated heritability and repeatability for tick counts were 0.19 ± 0.03 and 0.29 ± 0.01, respectively. Data were split into 5 subsets using k-means and random clustering for cross-validation of genomic predictions. Depending on the method, direct genomic value (DGV) prediction accuracies ranged from 0.35 with Bayes least absolute shrinkage and selection operator (LASSO) to 0.39 with BayesB for k-means clustering and between 0.42 with BayesLASSO and 0.45 with BayesC for random clustering. All genomic methods were superior to pedigree BLUP (PBLUP) accuracies of 0.26 for k-means and 0.29 for random groups, with highest accuracy gains obtained with BayesB (39%) for k-means and BayesC (55%) for random groups. Blending of historical phenotypic and pedigree information by different methods further increased DGV accuracies by values between 0.03 and 0.05 for direct prediction methods. However, highest accuracy was observed with single-step genomic BLUP with values of 0.48 for -means and 0.56, which represent, respectively, 84 and 93% improvement over PBLUP. Observed random clustering cross-validation breed-specific accuracies ranged between 0.29 and 0.36 for HH and between 0.55 and 0.61 for BO, depending on the blending method. These moderately high values for BO demonstrate that genomic predictions could be used as a practical tool to improve genetic resistance to ticks and in the development of resistant lines of this breed. For HH, accuracies are still in the low to moderate side and this breed training population needs to be increased before genomic selection could be reliably applied to improve tick resistance.

Analyses of reaction norms reveal new chromosome regions associated with tick resistance in cattle

Despite single nucleotide polymorphism (SNP) availability and frequent cost reduction has allowed genome-wide association studies even in complex traits as tick resistance, the use of this information source in SNP by environment interaction context is unknown for many economically important traits in cattle. We aimed at identifying putative genomic regions explaining differences in tick resistance in Hereford and Braford cattle under SNP by environment point of view as well as to identify candidate genes derived from outliers/significant markers. The environment was defined as contemporary group means of tick counts, since they seemed to be the most appropriate entities to describe the environmental gradient in beef cattle. A total of 4363 animals having tick counts (n=10 673) originated from 197 sires and 3966 dams were used. Genotypes were acquired on 3591 of these cattle. From top 1% SNPs (410) having the greatest effects in each environment, 75 were consistently relevant in all environments, which indicated SNP by environment interaction. The outliers/significant SNPs were mapped on chromosomes 1, 2, 5, 6, 7, 9, 11, 13, 14, 15, 16, 18, 21, 23, 24, 26 and 28, and potential candidate genes were detected across environments. The presence of SNP by environment interaction for tick resistance indicates that genetic expression of resistance depends upon tick burden. Markers with major portion of genetic variance explained across environments appeared to be close to genes with different direct or indirect functions related to immune system, inflammatory process and mechanisms of tissue destruction/repair, such as energy metabolism and cell differentiation.

Genomewide association study for production and meat quality traits in Canchim beef cattle

The commercial value of the bovine carcass is determined by a set of traits, such as weight, yield, back fat thickness, and marbling; therefore, the genetic improvement of growth, meat, and carcass quality traits is an important tool to add value to the supply chain. Genomewide association studies (GWAS) enable the identification of loci that control phenotypic expression of quantitative traits (QTL). Therefore, the objective of this work was to perform a GWAS to identify genomic regions and genes associated with growth, carcass traits, and meat quality in Canchim beef cattle. These traits were yearling weight (YW), rib eye area (REA), back fat thickness (BFT), and marbling (MARB). To increase sample size and marker density, genotype imputation was performed, and only markers imputed with greater than 95% accuracy were used. Genomewide association study was performed using a Bayesian approach, by the Bayes B statistical method, incorporating genotypes and phenotypes from 614 animals from both the Canchim breed and the MA genetic group (offspring of Charolais bulls and one-half Canchim + one-half Zebu cows). This investigation identified 1 and 4 genomic regions explaining 0.23 and 7.35% of the genetic variance for REA and YW, respectively. These regions harbor a total of 19 genes, 7 of which were classified for biological functions by functional analysis. Significant associations were not observed for BFT and MARB. The identification of QTL that had been previously described in the literature reinforces associations found in this study.

Genomic prediction using different estimation methodology, blending and cross-validation techniques for growth traits and visual scores in Hereford and Braford cattle

The objective of the present study was to evaluate the accuracy and bias of direct and blended genomic predictions using different methods and cross-validation techniques for growth traits (weight and weight gains) and visual scores (conformation, precocity, muscling, and size) obtained at weaning and at yearling in Hereford and Braford breeds. Phenotypic data contained 126,290 animals belonging to the Delta G Connection genetic improvement program, and a set of 3,545 animals genotyped with the 50K chip and 131 sires with the 777K. After quality control, 41,045 markers remained for all animals. An animal model was used to estimate (co)variance components and to predict breeding values, which were later used to calculate the deregressed estimated breeding values (DEBV). Animals with genotype and phenotype for the traits studied were divided into 4 or 5 groups by random and k-means clustering cross-validation strategies. The values of accuracy of the direct genomic values (DGV) were moderate to high magnitude for at weaning and at yearling traits, ranging from 0.19 to 0.45 for the k-means and 0.23 to 0.78 for random clustering among all traits. The greatest gain in relation to the pedigree BLUP (PBLUP) was 9.5% with the BayesB method with both the k-means and the random clustering. Blended genomic value accuracies ranged from 0.19 to 0.56 for k-means and from 0.21 to 0.82 for random clustering. The analyses using the historical pedigree and phenotypes contributed additional information to calculate the GEBV, and in general, the largest gains were for the single-step (ssGBLUP) method in bivariate analyses with a mean increase of 43.00% among all traits measured at weaning and of 46.27% for those evaluated at yearling. The accuracy values for the marker effects estimation methods were lower for k-means clustering, indicating that the training set relationship to the selection candidates is a major factor affecting accuracy of genomic predictions. The gains in accuracy obtained with genomic blending methods, mainly ssGBLUP in bivariate analyses, indicate that genomic predictions should be used as a tool to improve genetic gains in relation to the traditional PBLUP selection.

Threshold and linear models for genetic evaluation of visual scores in Hereford and Braford cattle

Data from 127 539 Hereford and Braford cattle were used to compare estimates of genetic parameters for navel, conformation, precocity, muscling and size visual scores at yearling, using linear and threshold animal models. In a second step, these models were cross-validated using a multinomial logistic regression in order to quantify the association between phenotype and genetic merit for each trait. For navel score, higher heritability was obtained with the threshold model (0.42 ± 0.02) in relation to the linear model (0.22 ± 0.02). However, similar heritability was estimated in both models for conformation, precocity, muscling and size, with values of 0.18 ± 0.01, 0.19 ± 0.01, 0.19 ± 0.01 and 0.26 ± 0.01, respectively, using linear model, and of 0.19 ± 0.01, 0.19 ± 0.01, 0.20 ± 0.01, and 0.29 ± 0.01, respectively, using threshold model. For navel score, Spearman correlations between sires’ breeding values predicted using linear and threshold models ranged from 0.60 (1% of the best sires are selected) to 0.96 (all sires are selected). For conformation, precocity, muscling and size scores, low changes in sires’ rank are expected using these models (Spearman correlations >0.86), regardless of the proportion of sires selected. Except for navel with the linear model, the direction of the associations between phenotype and genetic merit were in accordance with its expectation, as there were increases in the phenotype per unit of change in the breeding value. Thus, the threshold model would be recommended to perform genetic evaluation of navel score in this population. However, linear and threshold models showed similar predictive ability for conformation, precocity, muscling and size scores.