Cattle production is dependent upon fertility because it results in producing offspring to offset production costs. A number of semen attributes are believed to affect fertility and are frequently ...measured as part of routine breeding soundness exams or semen collection procedures. The objective of this study was to perform a single-step genome-wide association study (ssGWAS) for beef bull semen attributes. Beef bull fertility phenotypes including volume (VOL), concentration (CONC), number of spermatozoa (NSP), initial motility (IMot), post-thaw motility (PTMot), three-hour post-thaw motility (3HRPTMot), percentage of normal spermatozoa (%NORM), primary abnormalities (PRIM), and secondary abnormalities (SEC) were obtained from two artificial insemination (AI) centers. A total of 1819 Angus bulls with 50,624 collection records were used for ssGWAS. A five-generation pedigree was obtained from the American Angus Association and consisted of 6521 sires and 17,136 dams. Genotypes on 1163 bulls were also obtained from the American Angus Association and utilized in ssGWAS.
A multi-trait animal model was used for the estimation of single nucleotide polymorphism (SNP) effects. Significant SNP were those with a -log
P-value threshold greater than 4.0. Volume, CONC, NSP, IMot, PTMot, 3HRPTMot, %NORM, PRIM, and SEC have five, three, six, seven, two, six, six, and two genome-wide significant SNP, respectively.
Several significant SNP were determined to be near or within quantitative trait loci (QTL) associated with beef bull semen attributes. In addition, genes associated with fertility were found to contain or be near the significant SNP found in the study. The results indicate there are regions of the genome that impact fertility, proving inclusion of genomic information into genetic evaluation should be advantageous for genetic improvement of male fertility traits.
Estimated breeding values for average daily feed intake (AFI; kg/day), residual feed intake (RFI; kg/day) and average daily gain (ADG; kg/day) were generated using a mixed linear model incorporating ...genomic relationships for 698 Angus steers genotyped with the Illumina BovineSNP50 assay. Association analyses of estimated breeding values (EBVs) were performed for 41 028 single nucleotide polymorphisms (SNPs), and permutation analysis was used to empirically establish the genome‐wide significance threshold (P < 0.05) for each trait. SNPs significantly associated with each trait were used in a forward selection algorithm to identify genomic regions putatively harbouring genes with effects on each trait. A total of 53, 66 and 68 SNPs explained 54.12% (24.10%), 62.69% (29.85%) and 55.13% (26.54%) of the additive genetic variation (when accounting for the genomic relationships) in steer breeding values for AFI, RFI and ADG, respectively, within this population. Evaluation by pathway analysis revealed that many of these SNPs are in genomic regions that harbour genes with metabolic functions. The presence of genetic correlations between traits resulted in 13.2% of SNPs selected for AFI and 4.5% of SNPs selected for RFI also being selected for ADG in the analysis of breeding values. While our study identifies panels of SNPs significant for efficiency traits in our population, validation of all SNPs in independent populations will be necessary before commercialization.
ABSTRACT
Feed costs are a major economic expense in finishing and developing cattle; however, collection of feed intake data is costly. Examining relationships among measures of growth and intake, ...including breed differences, could facilitate selection for efficient cattle. Objectives of this study were to estimate genetic parameters for growth and intake traits and compare indices for feed efficiency to accelerate selection response. On-test ADFI and on-test ADG (TESTADG) and postweaning ADG (PWADG) records for 5,606 finishing steers and growing heifers were collected at the U.S. Meat Animal Research Center in Clay Center, NE. On-test ADFI and ADG data were recorded over testing periods that ranged from 62 to 148 d. Individual quadratic regressions were fitted for BW on time, and TESTADG was predicted from the resulting equations. We included PWADG in the model to improve estimates of growth and intake parameters; PWADG was derived by dividing gain from weaning weight to yearling weight by the number of days between the weights. Genetic parameters were estimated using multiple-trait REML animal models with TESTADG, ADFI, and PWADG for both sexes as dependent variables. Fixed contemporary groups were cohorts of calves simultaneously tested, and covariates included age on test, age of dam, direct and maternal heterosis, and breed composition. Genetic correlations (SE) between steer TESTADG and ADFI, PWADG and ADFI, and TESTADG and PWADG were 0.33 (0.10), 0.59 (0.06), and 0.50 (0.09), respectively, and corresponding estimates for heifers were 0.66 (0.073), 0.77 (0.05), and 0.88 (0.05), respectively. Indices combining EBV for ADFI with EBV for ADG were developed and evaluated. Greater improvement in feed efficiency can be expected using an unrestricted index versus a restricted index. Heterosis significantly affected each trait contributing to greater ADFI and TESTADG. Breed additive effects were estimated for ADFI, TESTADG, and the efficiency indices.
A study was performed to evaluate differences in thermoregulatory ability of 2 Bos taurus breeds with known differences in heat tolerance. Nine Angus (AG; 304 +/- 7 kg of BW) and 9 Romosinuano (RO; ...285 +/- 7.5 kg of BW) steers were transported to the Brody Environmental Center at the University of Missouri. Steers were housed for 18 d at thermoneutrality (TN; 21 degrees C) before initiation of heat stress (HS), which consisted of daily cyclic air temperature (26 degrees C, night; 36 degrees C, day) for 14 d. Rectal temperature and respiration rate were measured 6 times daily throughout the study. Sweat rates at shaved skin sites were recorded on specific days. Blood samples were taken once per week. Angus steers maintained rectal temperature 0.5 degrees C greater than RO at TN (P < 0.001). Likewise, respiration and sweat rates were greater (P < 0.001) in AG than RO at TN (P < 0.05). Rectal temperature increased during HS for both breeds with AG maintaining greater temperatures (P < 0.001). Both breeds increased respiration rate during HS, with AG steers exhibiting the greater rate (P < 0.001). Sweat rate increased more than 4-fold during HS (P < 0.001), followed by reduction after 7 d. Even after HS acclimation, AG exhibited the greater sweat rate (P < 0.001). Breed differences for serum leptin, creatinine, and cholesterol were found throughout the study with AG being greater than RO. Although there were no breed differences (P = 0.21) at TN, only AG steers exhibited a HS-induced increase (P < 0.05) in prolactin, creatinine, and cholesterol concentrations to suggest that an increase in rectal temperature is required for this effect. Use of rectal temperature along with endocrine markers, such as prolactin, may aid in the identification of B. taurus sensitivity to heat.
Accuracy of sire selection is limited by how well animals are characterized for their environment. The objective of this study was to evaluate the presence of genotype × environment interactions ...(G×E) for birth weight (BiW) and weaning weight (WW) for Red Angus in the United States. Adjusted weights were provided by the Red Angus Association of America. Environments were defined as 9 regions within the continental United States with similar temperature-humidity indices. Mean weights of calves were determined for each region and for each sire's progeny within each region. A reaction norm (RN) for each bull was estimated by regressing the sire means on the region means weighted for the number of progeny of each sire. The range for BiW and WW RN was -1.3 to 4.0 and -1.7 to 2.8, respectively. The heritabilities of BiW and WW RN were 0.40 and 0.39, respectively. Phenotypic and genetic correlations between BiW and WW RN were 0.19 and 0.54, respectively. The phenotypic correlation of the progeny mean to the RN was -0.20 ( <0.05) and suggests that sires with higher means are more stable in progeny performance across environments. Weights in different regions were considered separate traits and genetic correlations were estimated between all pairs of regions as another method to determine G×E. Genetic correlations < 0.80 indicate G×E at a level for concern, but existed for only 2 of 36 estimates for BiW and 12 of 36 estimates for WW. Genetic correlations between different regions ranged from 0.74 to 0.96 for BiW and 0.62 to 0.99 for WW and indicate that sires tend to rank similarly across environments for these traits.
Accuracy of sire selection is limited by how well animals are characterized for their environment. The objective of this study was to evaluate the presence of genotype × environment interactions ...(G×E) for birth weight (BiW) and weaning weight (WW) for Red Angus in the United States. Adjusted weights were provided by the Red Angus Association of America. Environments were defined as 9 regions within the continental United States with similar temperature-humidity indices. Mean weights of calves were determined for each region and for each sire's progeny within each region. A reaction norm (RN) for each bull was estimated by regressing the sire means on the region means weighted for the number of progeny of each sire. The range for BiW and WW RN was -1.3 to 4.0 and -1.7 to 2.8, respectively. The heritabilities of BiW and WW RN were 0.40 and 0.39, respectively. Phenotypic and genetic correlations between BiW and WW RN were 0.19 and 0.54, respectively. The phenotypic correlation of the progeny mean to the RN was -0.20 ( <0.05) and suggests that sires with higher means are more stable in progeny performance across environments. Weights in different regions were considered separate traits and genetic correlations were estimated between all pairs of regions as another method to determine G×E. Genetic correlations < 0.80 indicate G×E at a level for concern, but existed for only 2 of 36 estimates for BiW and 12 of 36 estimates for WW. Genetic correlations between different regions ranged from 0.74 to 0.96 for BiW and 0.62 to 0.99 for WW and indicate that sires tend to rank similarly across environments for these traits.
Brangus 3/8 Brahman (Bos indicus) × 5/8 Angus (Bos taurus); n ≈ 800 heifers from 67 sires were used to estimate heritability and conduct a genome-wide association study (GWAS) for 2 binary fertility ...traits: first service conception (FSC) and heifer pregnancy (HPG). Genotypes were from 53,692 loci on the BovineSNP50 (Infinium Bead Chips, Illumina, San Diego, CA). Yearling heifers were estrous synchronized, bred by AI, and then exposed to natural service breeding. Reproductive ultrasound and DNA-based parentage testing were used to determine if the heifer conceived by AI or natural service, and code for FSC and HPG traits. Success rates for FSC and HPG were 53.3% and 78.0% ± 0.01%, and corresponding heritability estimates were 0.18 ± 0.07 and 0.10 ± 0.06, respectively. The models used in obtaining these heritability estimates and GWAS included fixed effects of year (i.e., 2005 to 2007), birth location, calving season, age of dam, and contemporary group. In GWAS, simultaneous associations of 1 Mb SNP windows with phenotype were undertaken with Bayes C analyses using GenSel software. The 1 Mb windows contained 21.3 ± 1.1 SNP. Analyses fitted a mixture model that treated SNP effects as random, with an assumed fraction pi = 0.9995 having no effect on phenotype. The windows that accounted for 1.0% of genetic variance were considered as QTL associated with FSC or HPG. Eighteen QTL existed on 15 chromosomes for the 2 traits. On average, each QTL accounted for 2.43% ± 0.2% of the genetic variance. Chromosome 8 harbored 2 QTL for FSC and 1 for HPG; however, these regions did not overlap. Chromosomes 3, 15, 16, 19, 24, 26, 27, 29, and X included QTL only for FSC, whereas chromosomes 2, 4, 10, 13, and 20 contained QTL only for HPG. The multitude of QTL detected for FSC and HPG in this GWAS involving Brangus heifers exemplifies the complex regulation of variation in heifer fertility traits of low heritability.
Deoxyribonucleic acid-based tests were used to assign paternity to 625 calves from a multiple-sire breeding pasture. There was a large variability in calf output and a large proportion of young bulls ...that did not sire any offspring. Five of 27 herd sires produced over 50% of the calves, whereas 10 sires produced no progeny and 9 of these were yearling bulls. A comparison was made between the paternity results obtained when using a DNA marker panel with a high (0.999), cumulative parentage exclusion probability (PE) and those obtained when using a marker panel with a lower PE (0.956). A large percentage (67%) of the calves had multiple qualifying sires when using the lower resolution panel. Assignment of the most probable sire using a likelihood-based method based on genotypic information resolved this problem in approximately 80% of the cases, resulting in 75% agreement between the 2 marker panels. The correlation between weaning weight, on-farm EPD based on pedigrees inferred from the 2 marker panels was 0.94 for the 24 bulls that sired progeny. Partial progeny assignments inferred from the lower resolution panel resulted in the generation of EPD for bulls that actually sired no progeny according to the high-PE panel, although the Beef Improvement Federation accuracies of EPD for these bulls were never greater than 0.14. Simulations were performed to model the effect of loci number, minor allele frequency, and the number of offspring per bull on the accuracy of genetic evaluations based on parentage determinations derived from SNP marker panels. The SNP marker panels of 36 and 40 loci produced EPD with accuracies nearly identical to those EPD resulting from use of the true pedigree. However, in field situations where factors including variable calf output per sire, large sire cohorts, relatedness among sires, low minor allele frequencies, and missing data can occur concurrently, the use of marker panels with a larger number of SNP loci will be required to obtain accurate on-farm EPD.
The objective was to estimate genetic parameters for udder traits in Hereford cattle. American Hereford Association (AHA) members initially recorded an overall score based on all udder ...characteristics. In 2008, the Beef Improvement Federation established guidelines, which were subsequently adopted by the AHA, for evaluating udder suspension and teat size. Therefore, a female was scored for either overall score or udder suspension and teat size for a single lactation, and females may be evaluated for overall score for a parity and then for udder suspension and teat size at a later parity. In all cases, subjective scores were assigned at parturition and ranged from 1 to 9, with a score of 9 considered ideal. Records on 48,191 animals and a 3-generation pedigree with 126,814 animals were obtained from the AHA, Kansas City, MO. These records contained repeated observations for overall score (n = 73,469), suspension (n = 38,412), and teat size (n = 38,412). Because the distribution of scores for all traits peaked at 7, a linear approximation was used in the analysis. Data were modeled using a multiple-trait animal model with random effects of additive genetic and permanent environment, fixed effect of contemporary group (herd-year-season), and a linear covariate for age in days. Heritability estimates (SE) for overall score, suspension, and teat size were 0.32 (0.01), 0.32 (0.01), and 0.28 (0.01), respectively. Through genetic selection for these traits, beef producers could improve udder traits. Repeatability estimates (SE) for overall score, suspension, and teat size were 0.45 (0.005), 0.47 (0.01), and 0.44 (0.01), respectively. Producers should continue evaluating udder traits repeatedly throughout a cow's lifetime. The phenotypic correlation (SE) between suspension and teat size was 0.64 (0.004) with 57% of records for suspension and teat size having the same score for both traits. The genetic correlations (SE) between teat size and suspension, overall score and teat size, and overall score and suspension were 0.81 (0.01), 0.71 (0.03), and 0.69 (0.03), respectively, and selection for one trait should result in correlated responses in the other traits. In conclusion, traits were moderately repeatable with scores from a parity being informative for subsequent parities. Because overall score, udder suspension, and teat size were moderately heritable with strong, positive genetic correlations, genetic improvement for these traits can be achieved through selection.
The objective of this study was to estimate parameters required for genetic evaluation of Simmental carcass merit using carcass and live animal data. Carcass weight, fat thickness, longissimus muscle ...area, and marbling score were available from 5,750 steers and 1,504 heifers sired by Simmental bulls. Additionally, yearling ultrasound measurements of fat thickness, longissimus muscle area, and estimated percentage of intramuscular fat were available on Simmental bulls (n = 3,409) and heifers (n = 1,503). An extended pedigree was used to construct the relationship matrix (n = 23,968) linking bulls and heifers with ultrasound data to steers and heifers with carcass data. All data were obtained from the American Simmental Association. No animal had both ultrasound and carcass data. Using an animal model and treating corresponding ultrasound and carcass traits separately, genetic parameters were estimated using restricted maximum likelihood. Heritability estimates for carcass traits were 0.48 ± 0.06, 0.35 ± 0.05, 0.46 ± 0.05, and 0.54 ± 0.05 for carcass weight, fat thickness, longissimus muscle area, and marbling score, respectively. Heritability estimates for bull (heifer) ultrasound traits were 0.53 ± 0.07 (0.69 ± 0.09), 0.37 ± 0.06 (0.51 ± 0.09), and 0.47 ± 0.06 (0.52 ± 0.09) for fat thickness, longissimus muscle area, and intramuscular fat percentage, respectively. Heritability of weight at scan was 0.47 ± 0.05. Using a bivariate weight model including scan weight of bulls and heifers with carcass weight of slaughter animals, a genetic correlation of 0.77 ± 0.10 was obtained. Models for fat thickness, longissimus muscle area, and marbling score were each trivariate, including ultrasound measurements on yearling bulls and heifers, and corresponding carcass traits of slaughter animals. Genetic correlations of carcass fat thickness with bull and heifer ultrasound fat were 0.79 ± 0.13 and 0.83 ± 0.12, respectively. Genetic correlations of carcass longissimus muscle area with bull and heifer ultrasound longissimus muscle area were 0.80 ± 0.11 and 0.54 ± 0.12, respectively. Genetic correlations of carcass marbling score with bull and heifer ultrasound intramuscular fat percentage were 0.74 ± 0.11 and 0.69 ± 0.13, respectively. These results provide the parameter estimates necessary for genetic evaluation of Simmental carcass merit using both data from steer and heifer carcasses, and their ultrasound indicators on yearling bulls and heifers.