Genome-wide copy number variation in the bovine genome detected using low coverage sequence of popular beef breeds.

Copy number variations (CNVs) are large insertions, deletions or duplications in the genome that vary between members of a species and are known to affect a wide variety of phenotypic traits. In this study, we identified CNVs in a population of bulls using low coverage next-generation sequence data. First, in order to determine a suitable strategy for CNV detection in our data, we compared the performance of three distinct CNV detection algorithms on benchmark CNV datasets and concluded that using the multiple sample read depth approach was the best method for identifying CNVs in our sequences. Using this technique, we identified a total of 1341 copy number variable regions (CNVRs) from genome sequences of 154 purebred sires used in Cycle VII of the USMARC Germplasm Evaluation Project. These bulls represented the seven most popular beef breeds in the United States: Hereford, Charolais, Angus, Red Angus, Simmental, Gelbvieh and Limousin. The CNVRs covered 6.7% of the bovine genome and spanned 2465 protein-coding genes and many known quantitative trait loci (QTL). Genes harbored in the CNVRs were further analyzed to determine their function as well as to find any breed-specific differences that may shed light on breed differences in adaptation, health and production.

Linkage disequilibrium among commonly genotyped SNP variants detected from bull sequence.

Genomic prediction utilizing causal variants could increase selection accuracy above that achieved with SNPs genotyped by currently available arrays used for genomic selection. A number of variants detected from sequencing influential sires are likely to be causal, but noticeable improvements in prediction accuracy using imputed sequence variant genotypes have not been reported. Improvement in accuracy of predicted breeding values may be limited by the accuracy of imputed sequence variants. Using genotypes of SNPs on a high-density array and non-synonymous SNPs detected in sequence from influential sires of a multibreed population, results of this examination suggest that linkage disequilibrium between non-synonymous and array SNPs may be insufficient for accurate imputation from the array to sequence. In contrast to 75% of array SNPs being strongly correlated to another SNP on the array, less than 25% of the non-synonymous SNPs were strongly correlated to an array SNP. When correlations between non-synonymous and array SNPs were strong, distances between the SNPs were greater than separation that might be expected based on linkage disequilibrium decay. Consistently near-perfect whole-genome linkage disequilibrium between the full array and each non-synonymous SNP within the sequenced bulls suggests that whole-genome approaches to infer sequence variants might be more accurate than imputation based on local haplotypes. Opportunity for strong linkage disequilibrium between sequence and array SNPs may be limited by discrepancies in allele frequency distributions, so investigating alternate genotyping approaches and panels providing greater chances of frequency-matched SNPs strongly correlated to sequence variants is also warranted. Genotypes used for this study are available from https://www.animalgenome.org/repository/pub/;USDA2017.0519/.

Ruminal expression of the NQO1, RGS5, and ACAT1 genes may be indicators of feed efficiency in beef steers.

Ruminal genes differentially expressed in crossbred beef steers from USMARC with variation in gain and feed intake were identified in a previous study. Several of the genes identified with expression patterns differing between animals with high gain-low feed intake and low gain-high feed intake were evaluated in a separate, unrelated population of Angus × Hereford beef steers from the University of Wyoming that was classified to differ in residual feed intake (RFI). Of the 17 genes tested, two were differentially expressed by RFI class in the Angus × Hereford animals. These genes included NAD(P)H dehydrogenase, quinone 1 (NQO1; P = 0.0009) and regulator of G-protein signaling 5 (RGS5; P = 0.01). A third gene, acetyl-CoA acetyltransferase 1 (ACAT1; P = 0.06), displayed a trend toward association with RFI. These data suggest that some of the genes identified in a previous rumen transcriptome discovery study may have utility for identifying or selecting for animals with superior feed efficiency phenotypes across cattle breeds and populations.

Genetic markers that influence feed efficiency phenotypes also affect cattle temperament as measured by flight speed.

Flight speed is a predictive indicator of cattle temperament and is associated with feed efficiency phenotypes. Genetic markers associated with both traits may assist with selection of calmer animals with improved economic value. A preliminary genome-wide association study determined chromosomal regions on BTA9, and 17 were associated with flight speed. The genes quaking (QKI), glutamate receptor, ionotropic, AMPA 2 (GRIA2) and glycine receptor ß (GLRB) were identified in these regions as potential functional candidates. Beef steers (n = 1057) were genotyped with SNPs located within and flanking these genes. One SNP located near QKI and one near GRIA2 were nominally associated with flight speed (P ≤ 0.05) although neither was significant after Bonferroni correction. Several studies have shown a correlation between flight speed and feed intake or gain; therefore, we also analyzed SNPs on BTA6:38-39 Mb known to be associated with average daily gain (ADG) and average daily feed intake (ADFI) for association with flight speed. Several SNPs on BTA6 were associated with flight speed (P ≤ 0.005), and three were significant after Bonferroni correction. These results suggest that the genes tested are unlikely to contribute to flight speed variation for our cattle population, but SNPs on BTA6 associated with ADG and ADFI may influence temperament. Use of these markers to select for economically important feed efficiency phenotypes may produce cattle with more desirable temperaments.

A region on BTA14 that includes the positional candidate genes LYPLA1, XKR4 and TMEM68 is associated with feed intake and growth phenotypes in cattle(1).

Feed cost for beef cattle is the largest expense incurred by cattle producers. The development of genetic markers to enhance selection of more efficient animals that require less feed while still achieving acceptable levels of production has the potential to substantially reduce production costs. A genome-wide marker association approach based on the Illumina BovineSNP50 BeadChip™ was used to identify genomic regions affecting average daily feed intake (ADFI), average daily gain (ADG) and residual feed intake traits in a population of 1159 crossbred steers. This approach identified a region on BTA14 from 22.02 to 23.92 Mb containing several single-nucleotide polymorphisms (SNPs) that have significant association with at least one of the traits. Two genes in this region, lysophospholipase 1 (LYPLA1) and transmembrane protein 68 (TMEM68), appeared to be logical positional and functional candidate genes. LYPLA1 deacylates ghrelin, a hormone involved in the regulation of appetite in the rat stomach, while TMEM68 is expressed in bovine rumen, abomasum, intestine and adipose tissue in cattle, and likely affects lipid biosynthetic processes. SNPs lying in or near these two genes were identified by sequencing a subset of animals with extreme phenotypes. A total of 55 SNPs were genotyped and tested for association with the same population of steers. After correction for multiple testing, five markers within 22.79-22.84 Mb, located downstream of TMEM68, and between TMEM68 and the neighbouring gene XKR4, were significant for both ADFI and ADG. Genetic markers predictive of feed intake and weight gain phenotypes in this population of cattle may be useful for the identification and selection of animals that consume less feed, although further evaluation of these markers for effects on other production traits and validation in additional populations will be required.

Genetic markers on BTA14 predictive for residual feed intake in beef steers and their effects on carcass and meat quality traits.

With the high cost of feed for animal production, genetic selection for animals that metabolize feed more efficiently could result in substantial cost savings for cattle producers. The purpose of this study was to identify DNA markers predictive for differences among cattle for traits associated with feed efficiency. Crossbred steers were fed a high-corn diet for 140 days and average daily feed intake (ADFI), average daily gain (ADG), and residual feed intake (RFI) phenotypes were obtained. A region on chromosome 14 was previously associated with RFI in this population of animals. To develop markers with the highest utility for predicting an animal's genetic potential for RFI, we genotyped additional markers within this chromosomal region. These polymorphisms were genotyped on the same animals (n = 1066) and tested for association with ADFI, ADG and RFI. Six markers within this region were associated with RFI (P ≤ 0.05). After conservative correction for multiple testing, one marker at 25.09 Mb remained significant (P = 0.02) and is responsible for 3.6% of the RFI phenotypic variation in this population of animals. Several of these markers were also significant for ADG, although none were significant after correction. Marker alleles with positive effects on ADG corresponded to lower RFI, suggesting an effect increasing growth without increasing feed intake. All markers were also assessed for their effects on meat quality and carcass traits. All of the markers associated with RFI were associated with adjusted fat thickness (AFT, P ≤ 0.009) and three were also associated with hot carcass weight (HCW, P ≤ 0.003). Marker alleles associated with lower RFI were also associated with reduced AFT, and if they were associated for HCW, the effect was an increase in weight. These markers may be useful as prediction tools for animals that utilize feed more efficiently; however, validation with additional populations of cattle is required.

Porcine NAMPT gene: search for polymorphism, mapping and association studies.

NAMPT encodes an enzyme catalysing the rate-limiting step in NAD biosynthesis. The extracellular form of the enzyme is known as adipokine visfatin. We detected SNP AM999341:g.669T>C (referred to as 669T>C) in intron 9 and SNP FN392209:g.358A>G (referred to as 358A>G) in the promoter of the gene. RH mapping linked the gene to microsatellite SW944. Linkage analysis placed the gene on the current USDA ­ USMARC linkage map at position 92 cM on SSC9. Association analyses were performed in a wild boar × Meishan F2 family (W × M), with 45 traits recorded (growth and fattening, fat deposition, muscling, meat quality, stress resistance and other traits), and in a commercial Landrace × Chinese-European (LCE) synthetic population with records for 15 traits (growth, fat deposition, muscling, intramuscular fat, meat colour and backfat fatty acid content). In the W × M, SNP 669T>C was associated with muscling, fat deposition, growth and fattening, meat quality and other traits and in the LCE with muscling, meat quality and backfat fatty acid composition. In the W × M, SNP 358A>G was associated with muscling, fat deposition, growth and other traits. After correction for multiple testing, the NAMPT haplotypes were associated in the W × M with, in descending order, muscling (q = 0.0056), growth (q = 0.0056), fat deposition (q = 0.0109), fat-to-meat ratio (q = 0.0135), cooling losses (q = 0.0568) and longissimus pHU (q = 0.0695). The SNPs are hypothesized to be in linkage disequilibrium with a causative mutation affecting energy metabolism as a whole rather than fat metabolism alone.

Genetic variance and covariance and breed differences for feed intake and average daily gain to improve feed efficiency in growing cattle.

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.

Breed effects and genetic parameter estimates for calving difficulty and birth weight in a multibreed population.

Birth weight (BWT) and calving difficulty (CD) were recorded on 4,579 first-parity females from the Germplasm Evaluation Program at the U.S. Meat Animal Research Center (USMARC). Both traits were analyzed using a bivariate animal model with direct and maternal effects. Calving difficulty was transformed from the USMARC scores to corresponding -scores from the standard normal distribution based on the incidence rate of the USMARC scores. Breed fraction covariates were included to estimate breed differences. Heritability estimates (SE) for BWT direct, CD direct, BWT maternal, and CD maternal were 0.34 (0.10), 0.29 (0.10), 0.15 (0.08), and 0.13 (0.08), respectively. Calving difficulty direct breed effects deviated from Angus ranged from -0.13 to 0.77 and maternal breed effects deviated from Angus ranged from -0.27 to 0.36. Hereford-, Angus-, Gelbvieh-, and Brangus-sired calves would be the least likely to require assistance at birth, whereas Chiangus-, Charolais-, and Limousin-sired calves would be the most likely to require assistance at birth. Maternal breed effects for CD were least for Simmental and Charolais and greatest for Red Angus and Chiangus. Results showed that the diverse biological types of cattle have different effects on both BWT and CD. Furthermore, results provide a mechanism whereby beef cattle producers can compare EBV for CD direct and maternal arising from disjoined and breed-specific genetic evaluations.

Genomewide association study of liver abscess in beef cattle.

Fourteen percent of U.S. cattle slaughtered in 2011 had liver abscesses, resulting in reduced carcass weight, quality, and value. Liver abscesses can result from a common bacterial cause, , which inhabits rumen lesions caused by acidosis and subsequently escapes into the blood stream, is filtered by the liver, and causes abscesses in the liver. Our aim was to identify SNP associated with liver abscesses in beef cattle. We used lung samples as a DNA source because they have low economic value, they have abundant DNA, and we had unrestricted access to sample them. We collected 2,304 lung samples from a beef processing plant: 1,152 from animals with liver abscess and 1,152 from animals without liver abscess. Lung tissue from pairs of animals, 1 with abscesses and another without, were collected from near one another on the viscera table to ensure that pairs of phenotypically extreme animals came from the same lot. Within each phenotype (abscess or no abscess), cattle were pooled by slaughter sequence into 12 pools of 96 cattle for each phenotype for a total of 24 pools. The pools were constructed by equal volume of frozen lung tissue from each animal. The DNA needed to allelotype each pool was then extracted from pooled lung tissue and the BovineHD Bead Array (777,962 SNP) was run on all 24 pools. Total intensity (TI), an indicator of copy number variants, was the sum of intensities from red and green dyes. Pooling allele frequency (PAF) was red dye intensity divided TI. Total intensity and PAF were weighted by the inverse of their respective genomic covariance matrices computed over all SNP across the genome. A false discovery rate ≤ 5% was achieved for 15 SNP for PAF and 20 SNP for TI. Genes within 50 kbp from significant SNP were in diverse pathways including maintenance of pH homeostasis in the gastrointestinal tract, maintain immune defenses in the liver, migration of leukocytes from the blood into infected tissues, transport of glutamine into the kidney in response to acidosis to facilitate production of bicarbonate to increase pH, aggregate platelets to liver injury to facilitate liver repair, and facilitate axon guidance. Evidence from the 35 detected SNP associations combined with evidence of polygenic variation indicate that there is adequate genetic variation in incidence rate of liver abscesses, which could be exploited to select sires for reduced susceptibility to subacute acidosis and associated liver abscess.

Evaluation of single nucleotide polymorphisms in chromosomal regions impacting pregnancy status in cattle.

Reproductive success is an important component of commercial beef cattle production, and identification of DNA markers with predictive merit for reproductive success would facilitate accurate prediction of mean daughter pregnancy rate, enabling effective selection of bulls to improve female fertility. A previous study identified SNP associated with beef cattle reproductive efficiency based on a genomewide association analysis approach using genotyping multiple-animal pools of DNA to increase the number of animals that could be genotyped with available resources. For the current study, we expand on this previous study by individually genotyping cattle from the pooling study for 89 SNP that were associated with female pregnancy rate. The aims of the study were to confirm the results of the pooling study and, more specifically, identify modes of gene action and DNA variations such as haplotypes that would not be possible with pooled genotyping. Eighty-nine SNP selected from the pooling study were evaluated using the Sequenom MassARRAY system to individually genotype animals from populations evaluated in the pooling study, including both and breeds. From this research, regions on chromosomes 5 (26.3-48.1 Mb; UMD3.1 assembly) and 9 (37,436,575 bp; UMD3.1 assembly), first identified in the previous pooling study, were shown through individual genotyping to harbor genetic variation ( < 0.05 genomewide significance) affecting reproductive efficiency in interspecific crosses ( and ) of cattle. Each of these markers exhibited additive (vs. dominant) gene action. Additionally, a haplotype block harboring an allele of origin with negative effects on reproduction was identified on chromosome 5 in interspecific composite breeds of × composites.

Differences in transcript abundance of genes on BTA15 located within a region associated with gain in beef steers.

Using results from a previous GWAS, we chose to evaluate seven genes located within a 229Kb region on BTA15 for variation in RNA transcript abundance in a library of tissue samples that included adipose, liver, rumen papillae, spleen, muscle, and small intestine epithelial layers from the duodenum, ileum and jejunum collected from steers (n = 14) with positive and negative residual GN near mean dry matter intake (DMI). The genes evaluated were two olfactory receptor-like genes (LOC525033 and LOC618173), RRM1, STIM1, RHOG, PGAP2, and NUP98. The rumen papillae transcript abundance of RHOG was positively correlated with residual GN (P = 0.02) and ruminal STIM1 exhibited a trend towards an association with residual GN (P = 0.08). The transcript abundance of one olfactory receptor (LOC618173) in the ileum was also positively associated with residual GN (P = 0.02) and PGAP2 and LOC525033 in the ileum displayed trends for association with GN (P ≤ 0.1). To further evaluate the differential expression detected in the ileum and rumen of these animals, the transcript abundance of STIM1 and RHOG in the rumen and of PGAP2 and the olfactory receptors in the ileum were assessed in an additional group of 32 animals with divergent average daily gain (ADG) and average daily feed intake (ADFI) collected over two groups. The olfactory receptor, LOC525033, was not expressed in the ileum for the majority of these animals. Only RHOG showed a slight, but non-significant trend towards greater expression in animals with greater gain. We have detected differences in the transcript abundance of genes within this region in the rumen and ileum of animals selected for greater and less residual gain; however, we were unable to validate the expression of these genes in the larger group of cattle possibly due to the differences in phenotype or contemporary group.

A survey of polymorphisms detected from sequences of popular beef breeds.

The genome sequence was obtained from 270 sires used in the Germplasm Evaluation (GPE) project. These bulls included 154 purebred AI sires from GPE Cycle VII breeds (Hereford, Angus, Simmental, Limousin, Charolais, Gelbvieh, and Red Angus), 83 F crosses of those breeds, and 33 AI sires from 8 other breeds. The exome capture sequence targeting coding regions of the genome was obtained from 176 of these bulls. Sequence reads were mapped to the UMD 3.1 bovine genome assembly; a mean of 2.5-fold (x) coverage per bull was obtained from the genomic sequence, and the targeted exons were covered at a mean of 20.0x. Over 28.8 million biallelic sequence variants were detected where each allele was present in at least 3 different bulls. These included 22.0 million previously reported variants and 94.1% of the 774,660 autosomal and BTA X SNP on the BovineHD BeadChip assay (HD). More than 92% of the variants detected in targeted exons were also detected from the low-coverage genome sequence. Less than 1% of the variants detected from the combined genome and exome sequence occurred in annotated protein-coding sequences and 5' and 3' untranslated regions (UTR) surrounding the 19,994 annotated protein coding regions. Variation was detected in the coding sequence or UTR of 96.8% of the genes: loss-of-function variants were predicted for 3,298 genes, 14,973 contained nonsynonymous variants, 11,276 had variation in UTR, and 17,721 genes contained synonymous variants. Minor allele frequencies (MAF) were <0.05 for 47.8% of the coding sequence and UTR variants, and MAF distributions were skewed toward low MAF. In contrast, 11.1% of the HD SNP detected in these bulls had MAF < 0.05, and the distribution was skewed toward higher MAF. Genes involved in immune system processes and immune response were overrepresented among those genes containing high MAF loss-of-function and nonsynonymous polymorphisms. Detected variants were submitted to the National Center for Biotechnology Information genetic variation database (dbSNP) under the handle MARC, batch GPE_Bull_GenEx.

Systematic error in genetic evaluation of miles city line 1 hereford cattle resulting from preadjustment for age of dam.

Differences in preweaning growth of calves nursing 2- and 3-yr-old dams compared with contemporaries nursing older dams are accentuated in the Miles City Line 1 Hereford herd relative to age-of-dam (AOD) effects implied by preadjustment of 205-d weight in national cattle evaluation. Mixed-model analyses of 205-d weight that fit random individual direct effects and maternal genetic and permanent environmental effects on 4,998 calves were conducted to 1) determine the magnitude of residual AOD effects after preadjustment (PA) using industry-standard procedures and 2) compare changes in genetic predictions resulting from either PA or simultaneous adjustment (SA) for AOD. Expressed as differences from the 5- to 10-yr-old age effect, simultaneously estimated AOD effects were 45 +/- 1, 19 +/- 1, 6 +/- 1, and 19 +/- 3 kg for 2, 3, 4, and 11+ AOD classes, respectively. Comparable estimates of residual AOD effects after PA were 20 +/- 1, 6 +/- 1, 1 +/- 1, and 14 +/- 3 kg. Rank correlations of direct (BVd) and maternal (BVm) breeding values (BV) for 205-d weight from the analysis using PA with BV predicted using SA for AOD were .98 and .77, respectively. Estimated genetic trends were also affected by the method of accounting for AOD effects. One hundred fifty replicate simulations of 205-d weights with pedigree, fixed effect, and variance-covariance structures corresponding to the experimental population were used to establish correlations (r) of predicted BV with underlying true values. The r of predicted BVd with true values were reduced less than .02 by PA compared to SA in accounting for AOD. However, r of predicted BVm with true values were reduced more than .13 by PA compared to SA in accounting for AOD. These data indicate potential for systematic error in genetic evaluations that apply standard adjustments for AOD to 205-d weight.

Within-herd genetic analyses of stayability of beef females.

Nonlinear mixed-model procedures for analysis of binary data were used to estimate heritability (h2), predict individual genetic merit, and determine genetic and environmental trends for four measures of stayability of beef females. Traits considered were probabilities of a female having 2 [S(2/1)], 5 [S(5/1)], 8 [S(8/1)] and 11 [S(11/1)] calves, given that she calved once. Colorado State University Beef Improvement Center (BIC) and Beckton Stock Farm (BSF) provided data for the analyses. Heritability was estimated using animal model marginal maximum likelihood (AM MML), sire model marginal maximum likelihood (SM MML), and animal model Method R (AM MR). Individual genetic merit was predicted using single-trait animal models with each h2 estimate. Birth year was treated as fixed in all analyses. Only AM MML yielded h2 estimates for all traits in both herds. The AM MML h2 estimates for S(2/1), S(5/1), S(8/1), and S(11/1) were .09, .11, .07, and .20, respectively, for BSF data and .02, .14, .09, and .07, respectively, for BIC data. Differing h2 estimates did not substantially influence rank of individual predictions. Genetic trends in stayability were positive in both herds, although birth year solutions indicated variable or negative environmental trends. Genetic improvement of stayability may be accelerated by incorporating predictions of genetic merit for stayability in selection criteria. S(5/1) may be the most useful trait for consideration in national cattle evaluations.

Additive genetic groups for animals evaluated in more than one breed association national cattle evaluation.

Additive genetic groups were included in the 1993 Red Angus Association of America national cattle evaluation for phantom parents of individuals who were registered with the American Angus Association (AAA). Genetic groups were formed for each component in two multiple-trait evaluations in which all animal effects were fit. Additive direct effects were included for birth weight, weaning weight (WW), and milk (MILK). In a second analysis the additive direct effect of 160-d postweaning gain was analyzed with WW and MILK. Of the 387,665 animals, 50,838 had at least one phantom parent assigned to one of five genetic groups fit as fixed effects for each additive component. Of these 50,838 animals, 1,324 were identified as registered with the AAA. An average of 906 animals per component had an AAA EPD available. Animals with a known AAA EPD were designated into one of three groups of equal numbers based on AAA EPD for each component (1 = low, 2 = medium, 3 = high). Animals in the fourth genetic group were those registered with the AAA but with no EPD available for the component. The fifth genetic group included all other animals with phantom parents. Grouping on AAA EPD allowed for EPD on animals out of parent(s) registered with the AAA to be more closely aligned to the AAA EPD because they were regressed from the group solution instead of zero. Grouping based on EPD from another NCE should be considered in the production of multibreed EPD.

Comparison of selection by independent culling levels for below-average birth weight and high yearling weight with mass selection for high yearling weight in line 1 Hereford cattle.

Mass selection by independent culling levels (YB subline) for below-average birth weight (BWT) and high yearling weight (YWT) was compared with single-trait mass selection (YW subline) for high YWT in the inbred population of Line 1 Hereford cattle at Miles City, Montana. There were 4.2 generations of selection in YB and YW. Heritability estimates for the base population derived from multiple-trait REML were .28 and .31 for direct effects and .16 and .06 for maternal effects on BWT and YWT, respectively. Mid-parent cumulative selection differentials for BWT of YB and YW diverged (-2.9 vs 8.2 kg, respectively), as did the associated genetic trends for direct effects (-.014 kg/yr vs .105 kg/yr, respectively). Mid-parent cumulative selection differential for YWT of YB (102.1 kg) was 64% of that attained in YW (160.7 kg). Likewise, response in YWT of YB (.91 kg/yr) was 61% of response attained in YW (1.5 kg/yr). For BWT and YWT, maternal genetic trends were similar across selection lines. Assistance at parturition of first-parity 2-yr-old heifers was consistently less frequent in YB than in YW.

Assessment of positional candidate genes myf5 and igf1 for growth on bovine chromosome 5 in commercial lines of Bos taurus.

Quantitative trait loci for growth traits in beef cattle have been previously reported and fine-mapped in three chromosomal regions of 0 to 30 cM, 55 to 70 cM, and 70 to 80 cM of bovine chromosome 5. In this study, we further examined the association between gene-specific single nucleotide polymorphisms (SNP) of two positional candidate genes, bovine myogenic factor 5 (myf5) and insulin-like growth factor-1 (igf1), in the QTL regions and the birth weight (BWT), preweaning average daily gain (PWADG), and average daily gain on feed (ADGF) in commercial lines of Bos taurus. The QTL regions for the growth traits identified using a haplotype association analysis, which included the gene-specific SNP markers for both genes in this study, were in agreement with previous studies. The gene-specific SNP marker association analysis indicated that the SNP in myf5 had a significant additive effect on PWADG in the M1 line of Beefbooster Inc. (P < 0.10), and a significant additive effect (P < 0.05) and a significant dominance effect (P < 0.10) on ADGF in the M3 line of Beefbooster Inc. When the data from the two commercial lines were pooled, the SNP in myf5 showed a significant association with PWADG (P < 0.10) and with ADGF (P < 0.05). The association between the SNP and BWT, however, did not reach a significance level in the M1 line, the M3 line, or across the lines. For igf1, no significant association between the SNP and the growth traits was detected in either the M1 line or the M3 line, whereas there was only a significant dominance effect (P < 0.10) on BWT detected for the SNP in igfl when the data from the two commercial lines were pooled. These results suggest that myf5 is a strong candidate gene that influences PWADG and ADGF in beef cattle. The SNP of igf1 may not be a causative or close to the causative mutation that affects the three growth traits in the populations of beef cattle examined in this study. Other SNP of igf1 and myf5 or other genes in their respective chromosomal regions, however, should also be studied.