Bos taurus haplotypes segregating in Nellore (Bos indicus) cattle.

Brazil is the largest exporter of beef in the world, and most of that beef derives from Nellore cattle. Although considered a zebu breed (Bos indicus), the history of Nellore cattle in Brazil is marked by the importation of bulls from India, the use of a Creole taurine (Bos taurus) maternal lineage to quickly expand the herds and backcrossing to Nellore bulls to recover zebu ancestry. As a consequence, the current Brazilian Nellore population carries an average taurine ancestry of approximately 1%. Although that percentage seems small, some taurine variants deviate substantially from that average, with the better-known cases being the PLAG1-Q haplotype involved with body size variation and the Guarani (PG ) polled variant producing hornless animals. Here, we report taurine haplotypes in 9074 Nellore animals genotyped for 539 657 imputed SNP markers. Apart from PLAG1-Q and PG , our analysis further revealed common taurine haplotypes (>3%) spanning genes related to immunity, growth, reproduction and hair and skin phenotypes. Using data from 22 economically important traits, we showed that many of the major QTL previously reported in the breed are at least partially driven by taurine haplotypes. As B. taurus and B. indicus haplotypes are highly divergent, presenting widely different sets of functional variants, our results provide promising targets for future scrutiny in Nellore cattle.

Single-Nucleotide Polymorphism Variations Associated With Specific Genes Putatively Identified Enhanced Genetic Predisposition for 305-Day Milk Yield in the Girolando Crossbreed.

Milk production phenotypes are the main focus of genetic selection in dairy herds, and although there are many genes identified as related to the biology of these traits in pure breeds, little is known about crossbreed animals. This study aimed to identify potential genes associated with the 305-day milk yield in 337 crossbreed Gir × Holstein (Girolando) animals. Milk production records were genotyped for 45,613 single-nucleotide polymorphisms (SNPs). This dataset was used for a genome-wide association study (GWAS) using the 305-day milk yield adjusted for the fixed effects of herd and year and linear and quadratic effects of age at calving (in days) and calving factor averaged per animal. Genes within the significant SNPs were retrieved from the Bos taurus ARS-UCD1.2 assembly (bosTau9) for gene ontology analysis. In summary, the GWAS identified 52 SNPs associated [p ≤ 10-4, false discovery rate (FDR) = 8.77%] with milk production, including NUB1 and SLC24A2, which were previously described as related to milk production traits in cattle. The results suggest that SNPs associated mainly with NUB1 and SLC24A2 could be useful to understand milk production in Girolando and used as predictive markers for selecting genetic predisposition for milk yield in Girolando.

Crossbreeding effects on growth and efficiency in beef cow-calf systems: evaluation of Angus, Caracu, Hereford and Nelore breed direct, maternal and heterosis effects.

The objective of this study was to determine breed additive and heterosis effects on growth curves, total milk yield (TMY), calf weaning weight (WW), predicted energy intake (EI), and cow efficiency (CE) of purebred and crossbred beef cows raised in Southern Brazil. The data were from 175 purebred and crossbred cows representing eight genetic groups: Angus (A), Hereford (H), Nelore (N), A × H (AH), H × A (HA), A × N (AN), N × A (NA), and Caracu (C) × A (CA). Growth of the cows was modeled using the nonlinear Brody function and machine milking was used to assess TMY. WW was linearly adjusted to 210 d. EI was predicted with an equation in which the independent variables were estimates of parameters of the Brody function and TMY. The ratio of WW to EI estimated CE. Taurine-indicine heterosis effects were significant for all traits, and greater than those for taurine breed crosses. In general, crossbred cows were heavier at maturity, matured earlier, produced more milk, weaned heavier calves, and were predicted to consume more energy. Thus, they were more efficient than purebred cows, despite their greater predicted feed intake. Among the purebreds, A cows matured most rapidly, weighed the least at maturity, produced the most milk, weaned the heaviest calves, were predicted to consume the least energy; and were therefore most efficient among the breeds that were evaluated. These results are useful as inputs to bioeconomic models that can be used to predict productive and economic outcomes from crossbreeding and to facilitate recommendations for beef producers of southern Brazil and other similar subtropical climatic regions.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits.

BACKGROUND: The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (FROH), genomic relationship matrix (FGRM), and pedigree-based coefficient (FPED). RESULTS: The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of FPED and FROH, and their correlations (- 0.05 to 0.26) were low. Estimates of correlation between FGRM-FPED was zero, while the correlation (- 0.01 to - 0.07) between FGRM-FROH decreased as a function of ROH length, except for FROH > 8Mb (- 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands. CONCLUSIONS: Low FPED-FROH correlation estimates indicate that FPED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and FROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Accuracy of genomic predictions in Gyr (Bos indicus) dairy cattle.

Genomic selection may accelerate genetic progress in breeding programs of indicine breeds when compared with traditional selection methods. We present results of genomic predictions in Gyr (Bos indicus) dairy cattle of Brazil for milk yield (MY), fat yield (FY), protein yield (PY), and age at first calving using information from bulls and cows. Four different single nucleotide polymorphism (SNP) chips were studied. Additionally, the effect of the use of imputed data on genomic prediction accuracy was studied. A total of 474 bulls and 1,688 cows were genotyped with the Illumina BovineHD (HD; San Diego, CA) and BovineSNP50 (50K) chip, respectively. Genotypes of cows were imputed to HD using FImpute v2.2. After quality check of data, 496,606 markers remained. The HD markers present on the GeneSeek SGGP-20Ki (15,727; Lincoln, NE), 50K (22,152), and GeneSeek GGP-75Ki (65,018) were subset and used to assess the effect of lower SNP density on accuracy of prediction. Deregressed breeding values were used as pseudophenotypes for model training. Data were split into reference and validation to mimic a forward prediction scheme. The reference population consisted of animals whose birth year was ≤2004 and consisted of either only bulls (TR1) or a combination of bulls and dams (TR2), whereas the validation set consisted of younger bulls (born after 2004). Genomic BLUP was used to estimate genomic breeding values (GEBV) and reliability of GEBV (R2PEV) was based on the prediction error variance approach. Reliability of GEBV ranged from ∼0.46 (FY and PY) to 0.56 (MY) with TR1 and from 0.51 (PY) to 0.65 (MY) with TR2. When averaged across all traits, R2PEV were substantially higher (R2PEV of TR1 = 0.50 and TR2 = 0.57) compared with reliabilities of parent averages (0.35) computed from pedigree data and based on diagonals of the coefficient matrix (prediction error variance approach). Reliability was similar for all the 4 marker panels using either TR1 or TR2, except that imputed HD cow data set led to an inflation of reliability. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information. A reduced panel of ∼15K markers resulted in reliabilities similar to using HD markers. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information.

Variation in pestivirus growth in testicle primary cell culture is more dependent on the individual cell donor than cattle breed.

The causes of bovine respiratory disease complex (BRDC) are multifactorial and include infection with both viral and bacterial pathogens. Host factors are also involved as different breeds of cattle appear to have different susceptibilities to BRDC. Infection with bovine pestiviruses, including bovine viral diarrhea virus 1 (BVDV1), BVDV2 and 'HoBi'-like viruses, is linked to the development of BRDC. The aim of the present study was to compare the growth of different bovine pestiviruses in primary testicle cell cultures obtained from taurine, indicine and mixed taurine and indicine cattle breeds. Primary cells strains, derived from testicular tissue, were generated from three animals from each breed. Bovine pestivirus strains used were from BVDV-1a, BVDV-1b, BVDV-2a and 'HoBi'-like virus. Growth was compared by determining virus titers after one passage in primary cells. All tests were run in triplicate. Virus titers were determined by endpoint dilution and RT-qPCR. Statistical analysis was performed using one way analysis of variance (ANOVA) followed by the Tukey's Multiple Comparison Test (P˂0.05). Significant differences in virus growth did not correlate with cattle breed. However, significant differences were observed between cells derived from different individuals regardless of breed. Variation in the replication of virus in primary cell strains may reflect a genetic predisposition that favors virus replication.