Multiple Country and Breed Genomic Prediction of Tick Resistance in Beef Cattle.

Ticks cause substantial production losses for beef and dairy cattle. Cattle resistance to ticks is one of the most important factors affecting tick control, but largely neglected due to the challenge of phenotyping. In this study, we evaluate the pooling of tick resistance phenotyped reference populations from multi-country beef cattle breeds to assess the possibility of improving host resistance through multi-trait genomic selection. Data consisted of tick counts or scores assessing the number of female ticks at least 4.5 mm length and derived from seven populations, with breed, country, number of records and genotyped/phenotyped animals being respectively: Angus (AN), Brazil, 2,263, 921/1,156, Hereford (HH), Brazil, 6,615, 1,910/2,802, Brangus (BN), Brazil, 2,441, 851/851, Braford (BO), Brazil, 9,523, 3,062/4,095, Tropical Composite (TC), Australia, 229, 229/229, Brahman (BR), Australia, 675, 675/675, and Nguni (NG), South Africa, 490, 490/490. All populations were genotyped using medium density Illumina SNP BeadChips and imputed to a common high-density panel of 332,468 markers. The mean linkage disequilibrium (LD) between adjacent SNPs varied from 0.24 to 0.37 across populations and so was sufficient to allow genomic breeding values (GEBV) prediction. Correlations of LD phase between breeds were higher between composites and their founder breeds (0.81 to 0.95) and lower between NG and the other breeds (0.27 and 0.35). There was wide range of estimated heritability (0.05 and 0.42) and genetic correlation (-0.01 and 0.87) for tick resistance across the studied populations, with the largest genetic correlation observed between BN and BO. Predictive ability was improved under the old-young validation for three of the seven populations using a multi-trait approach compared to a single trait within-population prediction, while whole and partial data GEBV correlations increased in all cases, with relative improvements ranging from 3% for BO to 64% for TC. Moreover, the multi-trait analysis was useful to correct typical over-dispersion of the GEBV. Results from this study indicate that a joint genomic evaluation of AN, HH, BN, BO and BR can be readily implemented to improve tick resistance of these populations using selection on GEBV. For NG and TC additional phenotyping will be required to obtain accurate GEBV.

Tag-SNP selection using Bayesian genomewide association study for growth traits in Hereford and Braford cattle.

The aim of this study was to perform a Bayesian genomewide association study (GWAS) to identify genomic regions associated with growth traits in Hereford and Braford cattle, and to select Tag-SNPs to represent these regions in low-density panels useful for genomic predictions. In addition, we propose candidate genes through functional enrichment analysis associated with growth traits using Medical Subject Headings (MeSH). Phenotypic data from 126,290 animals and genotypes for 131 sires and 3,545 animals were used. The Tag-SNPs were selected with BayesB (π = 0.995) method to compose low-density panels. The number of Tag-single nucleotide polymorphism (SNP) ranged between 79 and 103 SNP for the growth traits at weaning and between 78 and 100 SNP for the yearling growth traits. The average proportion of variance explained by Tag-SNP with BayesA was 0.29, 0.23, 0.32 and 0.19 for birthweight (BW), weaning weight (WW205), yearling weight (YW550) and postweaning gain (PWG345), respectively. For Tag-SNP with BayesA method accuracy values ranged from 0.13 to 0.30 for k-means and from 0.30 to 0.65 for random clustering of animals to compose reference and validation groups. Although genomic prediction accuracies were higher with the full marker panel, predictions with low-density panels retained on average 76% of the accuracy obtained with BayesB with full markers for growth traits. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene products related to the growth traits. The proposed Tag-SNP panels may be useful for future fine mapping studies and for lower-cost commercial genomic prediction applications.

Economic selection index in small rural dairy farms.

This study aimed to calculate economic values (EVs) and economic selection indices for milk production systems in small rural properties. The traits 305-d milk yield in kg (MY), fat (FP) and protein (PP) percentage, daily fat (FY) and protein (PY) yield, cow live weight in kg (LW), calving interval (CI), and logarithm of daily somatic cell count (SCC) in milk were considered the goals and selection criteria. The production systems were identified from 29 commercial properties based on the inventory of revenues and costs and of zootechnical field data. Later, bioeconomic models were developed to calculate the productive performance, revenues, and costs concerning milk production to estimate EVs, which were calculated as the difference in annual profit with dairy production resulting from a change in one unit of the trait while keeping the others constant and dividing the value by the number of cows. After the EVs were known, ten economic selection indices were estimated for each system so they could be compared by modifying the selection criteria and calculating the relative importance of each selection criteria, the accuracy of the economic selection index, and response expected to the selection in USD, among other parameters. One of the systems detected was called less intensive (LS) and was characterized by having ten cows in lactation that produced 13·5 l/d and consumed 1·8 kg of concentrate/d. The second system detected was called more intensive (IS) and had 22 cows in lactation that produced 17·5 l/d and consumed 3·4 kg of concentrate/d. Monthly profits per cows in lactation of USD 2·60 and USD 68·77 were recorded for LS and IS, respectively. The EVs of the traits MY, FP, and PP were all positive, while for the other traits they were all negative in all situations. The best economic selection indices were those featuring selection criteria MY, LW, and CI, while the trait LW had the greatest importance in both systems. These results indicate that animal frame must be controlled in order to maximize the system's profit.

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.