Genetic Variants and Their Putative Effects on microRNA-Seed Sites: Characterization of the 3' Untranslated Region of Genes Associated with Temperament.

The 3' untranslated region has an important role in gene regulation through microRNAs, and it has been estimated that microRNAs regulate up to 50% of coding genes in mammals. With the aim of allelic variant identification of 3' untranslated region microRNA seed sites, the 3' untranslated region was searched for seed sites of four temperament-associated genes (CACNG4, EXOC4, NRXN3, and SLC9A4). The microRNA seed sites were predicted in the four genes, and the CACNG4 gene had the greatest number with 12 predictions. To search for variants affecting the predicted microRNA seed sites, the four 3' untranslated regions were re-sequenced in a Brahman cattle population. Eleven single nucleotide polymorphisms were identified in the CACNG4, and eleven in the SLC9A4. Rs522648682:T>G of the CACNG4 gene was located at the predicted seed site for bta-miR-191. Rs522648682:T>G evidenced an association with both exit velocity (p = 0.0054) and temperament score (p = 0.0097). The genotype TT had a lower mean exit velocity (2.93 ± 0.4 m/s) compared with the TG and GG genotypes (3.91 ± 0.46 m/s and 3.67 ± 0.46 m/s, respectively). The allele associated with the temperamental phenotype antagonizes the seed site, disrupting the bta-miR-191 recognition. The G allele of CACNG4-rs522648682 has the potential to influence bovine temperament through a mechanism associated with unspecific recognition of bta-miR-191.

Comparison of breeding value prediction for two traits in a Nellore-Angus crossbred population using different Bayesian modeling methodologies.

The objectives of this study were to 1) compare four models for breeding value prediction using genomic or pedigree information and 2) evaluate the impact of fixed effects that account for family structure. Comparisons were made in a Nellore-Angus population comprising F2, F3 and half-siblings to embryo transfer F2 calves with records for overall temperament at weaning (TEMP; n = 769) and Warner-Bratzler shear force (WBSF; n = 387). After quality control, there were 34,913 whole genome SNP markers remaining. Bayesian methods employed were BayesB (π̃ = 0.995 or 0.997 for WBSF or TEMP, respectively) and BayesC (π = 0 and π̃), where π̃ is the ideal proportion of markers not included. Direct genomic values (DGV) from single trait Bayesian analyses were compared to conventional pedigree-based animal model breeding values. Numerically, BayesC procedures (using π̃) had the highest accuracy of all models for WBSF and TEMP (ρ̂gg = 0.843 and 0.923, respectively), but BayesB had the least bias (regression of performance on prediction closest to 1, ß̂y,x = 2.886 and 1.755, respectively). Accounting for family structure decreased accuracy and increased bias in prediction of DGV indicating a detrimental impact when used in these prediction methods that simultaneously fit many markers.

Alternative parameterizations of relatedness in whole genome association analysis of pre-weaning traits of Nelore-Angus calves.

Gestation length, birth weight, and weaning weight of F2 Nelore-Angus calves (n = 737) with designed extensive full-sibling and half-sibling relatedness were evaluated for association with 34,957 SNP markers. In analyses of birth weight, random relatedness was modeled three ways: 1) none, 2) random animal, pedigree-based relationship matrix, or 3) random animal, genomic relationship matrix. Detected birth weight-SNP associations were 1,200, 735, and 31 for those parameterizations respectively; each additional model refinement removed associations that apparently were a result of the built-in stratification by relatedness. Subsequent analyses of gestation length and weaning weight modeled genomic relatedness; there were 40 and 26 trait-marker associations detected for those traits, respectively. Birth weight associations were on BTA14 except for a single marker on BTA5. Gestation length associations included 37 SNP on BTA21, 2 on BTA27 and one on BTA3. Weaning weight associations were on BTA14 except for a single marker on BTA10. Twenty-one SNP markers on BTA14 were detected in both birth and weaning weight analyses.