- Research article
- Open Access
Single nucleotide polymorphisms for feed efficiency and performance in crossbred beef cattle
BMC Genetics volume 15, Article number: 14 (2014)
This study was conducted to: (1) identify new SNPs for residual feed intake (RFI) and performance traits within candidate genes identified in a genome wide association study (GWAS); (2) estimate the proportion of variation in RFI explained by the detected SNPs; (3) estimate the effects of detected SNPs on carcass traits to avoid undesirable correlated effects on these economically important traits when selecting for feed efficiency; and (4) map the genes to biological mechanisms and pathways. A total number of 339 SNPs corresponding to 180 genes were tested for association with phenotypes using a single locus regression (SLRM) and genotypic model on 726 and 990 crossbred animals for feed efficiency and carcass traits, respectively.
Strong evidence of associations for RFI were located on chromosomes 8, 15, 16, 18, 19, 21, and 28. The strongest association with RFI (P = 0.0017) was found with a newly discovered SNP located on BTA 8 within the ELP3 gene. SNPs rs41820824 and rs41821600 on BTA 16 within the gene HMCN1 were strongly associated with RFI (P = 0.0064 and P = 0.0033, respectively). A SNP located on BTA 18 within the ZNF423 gene provided strong evidence for association with RFI (P = 0.0028). Genomic estimated breeding values (GEBV) from 98 significant SNPs were moderately correlated (0.47) to the estimated breeding values (EBVs) from a mixed animal model. The significant (P < 0.05) SNPs (98) explained 26% of the genetic variance for RFI. In silico functional analysis for the genes suggested 35 and 39 biological processes and pathways, respectively for feed efficiency traits.
This study identified several positional and functional candidate genes involved in important biological mechanisms associated with feed efficiency and performance. Significant SNPs should be validated in other populations to establish their potential utilization in genetic improvement programs.
As feed costs are a major factor influencing the profitability of beef cattle production, there are many endeavors to reduce these costs. Improving feed efficiency can be achieved by novel feeding strategies and genetic improvement technologies. Although residual feed intake (RFI) has emerged as one of the important feed efficiency traits for beef cattle , there are limitations with RFI for direct selection to improve feed efficiency industry-wide. These limitations are the expense and difficulty of recording an animal’s daily feed intake. Genomic approaches offer opportunities to select cattle that are more efficient, as once the relationships between genetic markers and feed efficiency are determined, this prediction can be applied to animals that are genotyped, but are not phenotyped with costly feed intake measurements .
Since 2000, advances in high-throughput genotyping and sequencing techniques have resulted in high density SNP chips, such as the Illumina BovineSNP50 BeadChip  being available. The use of the Bovine SNP50 in dairy cattle has increased the accuracy for predicting the genetic value of animals . In beef cattle, the use of such developments will benefit most traits such as feed efficiency and carcass traits which are difficult to measure or require the animals to be slaughtered for recording their phenotypes . Several genome wide association studies (GWAS) indicated that many genes affect feed efficiency traits and that the majority of these effects are small [6–11]. These studies reported many SNPs conferring genetic variation in feed efficiency. Nonetheless, although many SNPs were studied, the genetic architecture of feed efficiency was not completely explained.
Results from fine-mapping by Abo-Ismail et al. () suggested a list of candidate genes for further investigation to identify the causal mutations for feed efficiency within these genes . Discovery of the causal mutations within these genes could help explain the genetic architecture of feed efficiency. Furthermore, this approach could provide a panel of the most informative SNPs that could be used to predict feed efficiency accurately and affordably for producers. Therefore, the objectives of this study were to: (1) identify new SNPs for RFI and performance traits within candidate genes identified in previous GWAS studies; (2) estimate the proportion of variation in feed efficiency traits explained by the detected SNPs; (3) estimate the effect of detected SNPs on carcass traits to avoid undesirable correlated effects when selecting for feed efficiency; and (4) map the corresponding genes to a biological process and pathway to understand the biological meaning behind the detected associations. In this way it was hoped to identify causal mutations or to identify markers in strong linkage disequilibrium with such mutations.
Animals and phenotypic data
The study was approved from The University of Guelph Animal Care Committee based on the recommendations outlined in the Canadian Council on Animal Care (1993) guidelines.
Feed efficiency traits
Average daily dry matter intake (DMI), average daily gain (ADG), midpoint metabolic weight (MMWT), RFI and feed conversion ratio (FCR) phenotypes were measured on 726 crossbred beef cattle, heifers (38), steers (387), and bulls (301) at the University of Guelph’s Elora Beef Research Center (EBRC). Average breed compositions were formed by Angus (45.9%), Simmental (20.7%), Piedmontese (5%), Gelbvieh (4.2%), Charolais (2%) and Limousin (1.4%). Animals primarily originated from one of two University of Guelph herds (EBRC and NLARS), the Agriculture and Agri-Food Canada Kapuskasing Research Centre (KAP) or were purchased from producers in Ontario, Canada. Calves were weaned at approximately 200 days of age, and were involved in various post-weaning trials at the EBRC with different nutrition treatments. The body weights of the animals were recorded a number of times over the trials with most trials recording weights at least every four weeks.
The ADG for individual animals was calculated as a linear regression coefficient of their live weights on the actual days of measurement using the nlme package from R software . The MMBW was calculated as the midpoint body weight (kg) to the power 0.75. The DMI was calculated for each animal as total DMI divided by number of days for the test period. The RFI was calculated from the difference between the average of the animal’s actual daily DMI and its expected daily DMI . Expected DMI was determined through the regression coefficients estimated from the data through a multiple phenotypic regression model as follows:
Where, yijk is the total DMI for animal k during the feeding period, μ is the overall mean, β1 is the regression coefficient of the linear regression on ADG as determined through a linear regression of weights on days on trial as described, β2 is the regression coefficient of the linear regression on MMWT, sexi is the effect of ith sex, TTYj is the effect of jth treatment × trial × year (42 levels) and eijk is the residual random effect associated with the animal k and is the resulting RFI used in further analyses.
Carcass and meat quality traits
The association analysis of carcass and meat quality traits was carried out on 693-990 (depending on the trait) crossbred animals, including heifers (n = 33), steers (n = 705), and bulls (n = 252). In total 698 of these animals have RFI measures. All cattle were slaughtered at the University of Guelph Meat Science Laboratory Abattoir. Hot carcass weight (HCW) was measured just before the carcass was placed in the cooler. Meat Laboratory staff assessed the longissimus muscle interface (i.e. muscle surface) between the 12th and 13th ribs to obtain the following carcass measurements: subcutaneous fat depths between the 1st and 2nd, 2nd and 3rd, and 3rd and 4th quadrants of longissimus muscle (recorded as F1, F2 and F3, respectively), the grade fat (GRF), the minimum measurement of subcutaneous fat depth within the 4th quadrants of longissimus muscle and longissimus muscle area, measured using an electronic planimeter (MOP-3; Carl Zeiss Inc., Thornwood, NY) after acetate tracing (Bergen et al.). Canadian Beef Grading Agency formulae (http://www.beefgradingagency.ca/) were used to determine lean yield (LY), an estimate of the percentage of the carcass that is red meat. Marbling was assessed to determine the average amount, size and distribution of fat particles or deposits within longissimus muscle andwas scored as ≤3.0 = devoid; 3.1 to 4.0 = traces; 4.1 to 5.9 = slight; 6.0 to 7.0 = small to moderate; and ≥7.0 = slightly abundant to abundant. Rib dissection traits were also measured using a 4-6 rib section depending on the trial and year (physical separation of ribs 8-12 or 6-12, respectively). This procedure determines the amount of lean meat and bone, and a quantitative and qualitative assessment of fat depots (body, subcutaneous and intermuscular) within the rib to evaluate carcass composition. A complete description of carcass measurements was discussed by .
SNP discovery, DNA isolation and genotyping
Messenger RNA from seven tissue types (adipose, muscle, hypothalamus, duodenum, liver, lung and kidney) was extracted using TRIzol (Invitrogen). The tissue samples were collected from beef cattle at the Lacombe Research Centre in Alberta (Canada). RNA from 7 to 14 animals was pooled for each tissue before sequencing. Sequencing libraries were constructed from each RNA pool according to a standard protocol (mRNA Sequencing Sample Preparation Guide, Illumina, USA). Sequencing was performed on the Illumina Genome Analyzer II following the manufacturer’s recommendations. The resulting reads (more than 140 M) were mapped to transcript sequences from the reference bovine genome assembly (Btau4.0)  using maq 0.6.6 . More than 1.2 million SNPs were detected by comparing the aligned reads to the reference transcripts. From this list a subset of 300 SNPs from 215 candidate genes was selected based on SNP functional consequences assigned by NGS-SNP . An additional 158 coding SNPs were chosen from publicly available SNPs within the same candidate genes (Additional file 1). These genes were selected based on their proximity (on average distance 116,963 base pair) to significant SNPs identified in a previous study .
Tissue or blood samples were prepared and sent to Laboratory Services, University of Guelph, for genomic DNA extraction. Then, prepared DNA samples were sent to GeneSeek, Inc. for genotyping using a commercial platform for high-throughput SNP genotyping. In total, 1,032 animals, as assessed by the numerator relationship matrix using CFC, born subsequent to the animals used in the GWAS population  were genotyped for 458 SNP. The 300 SNPs identified through this work that were verified through genotyping have been submitted to dbSNP under the handle name “UALG”.
Quality control (QC) was done using the GenABEL package  in R software. Animals (n = 14) and individual SNPs (n = 5) with a low call-rate (<90%) were excluded from the analysis. Mean Identical By State (IBS) was 0.783 ± 0.0327. Animals (n = 1) with high estimation of IBS (≥0.95) were excluded. SNPs (n = 114) with a minor allele frequency (MAF) (< 1%) were excluded from the analysis of feed efficiency traits. Mean autosomal heterozygosity (HET) was 0.27 ± 0.036; animals (n = 6) with high HET (≥0.446) were excluded. Three hundred thirty nine SNPs and 727 animals passed all QC criteria where these SNPs were mapped to 180 corresponding genes (83,58, 24, 9, 4, and 2 genes including only 1, 2, 3, 4, 5 and 6 SNPs, respectively). The distribution of genotyped SNPs (339) across 29 chromosomes of the bovine genome is presented in Figure 1.
Single locus regression model (SLRM)
Genotypic data were coded as 0, 1, 2 corresponding to the number of minor alleles using GenABEL. In this model, phenotypes were regressed on the number of copies of a minor allele (0, 1, or 2) for estimating the allele’s substitution effect using ASReml 3 software . For feed efficiency traits, the univariate animal model was fitted as follows:
in which Yijkl is the trait measured in the lth animal of the jth herd-year of birth and the kth treatment-trial-year group; μ is the overall mean for the trait; Sexi is the fixed effect of the ith sex of lth animal; HYj is the fixed effect of the jth (17 level) herd-year of birth group; TTYk is the fixed effect of the kth (42 level) treatment trial-year of the test group; β1 is the regression coefficient of the linear regression on the number of copies of a minor allele; β2 is the regression coefficient of the linear regression on age at the end of the test period (AET) of the lth animal; βl is the regression coefficient of the linear regressions on proportion of AN, CH, LM, SM, PI, and GV breeds in the lth animal; β9 is the regression coefficient of the linear regression on the percentage of heterozygosity of the lth animal; al is the random additive genetic (polygenic) effect of the lth animal; and eijklm is the residual random effect associated with the lth animal. The TTY level that had less than three animals was excluded from the analysis. Phenotypes that were not within the mean ± 3 standard deviations for the respective trait were excluded from the analysis.
For carcass traits, the previous model (2) was modified to include the effect of the treatment trial-year-sex group instead of TTY and to include the fixed effect of the herd-year slaughter season instead of HY. Also, the effect of age at the end of the test period (day) was substituted by the age at slaughter (day).
The significance of associations was determined by an overall value of P < 0.05. To allow for multiple hypothesis-testing, chromosome wise false discovery rate (FDR) was used . A threshold of 5 and 20% FDR were used for strong and suggestive associations, respectively.
This model was fitted only for feed efficiency traits to consider genetic effects other than the additive effect. The model included the same effects in the SLRM, except that the allele substitution effect was replaced with the genotype effect. This model was not fitted for carcass traits to reduce the volume of results as the trait of primary interest was feed efficiency for this study.
Estimation of genetic variance explained by identified SNPs
The proportion of phenotypic variance in RFI explained by the full set of SNPs (339) that passed QC was estimated using the BayesC algorithm implemented in GenSel 3.13 software . Also, the proportion of the genetic variance of RFI explained by the set of significant SNPs for at least one of the feed efficiency traits using SLRM and/or the genotypic model was estimated. Missing genotypes were inferred using fastPHASE . Estimated breeding value (EBV) was determined with the SLRM without the regression on SNPs by ASReml. BayesC was then used to run the analysis with the two sets of SNPs (the full set [339 SNPs] and significant SNPs from the two models [98 SNPs]). Posterior residual and genetic variances were estimated after 41,000 iterations including 1,000 burn-in cycles. The proportion of genetic variance explained by the set of SNPs was estimated as the posterior genetic variance divided by phenotypic variance (posterior residual plus posterior genetic variance). In addition, the correlations between genomic breeding values predicted by estimated solutions and EBVs were estimated.
The significant (P < 0.05) SNPs (98) from the SLRM and genotype models for at least one feed efficiency trait from the association analysis were mapped to 74 genes. The list of the genes was submitted to DAVID 6.7 Beta software  for an in silico functional analysis. In DAVID, Gene ontology (GO) was used to identify functionally related genes. The genes were also mapped to biological pathways using web software in the Kyoto Encyclopedia of Genes and Genomes (KEGG) .
Results and discussion
Our goal in the current study was to identify informative or causal mutations for feed efficiency traits for use in Marker Assisted Selection (MAS). This would accelerate genetic improvement in beef cattle by improving the accuracy of selection and shortening intervals between generations . Genetic Improvement of feed efficiency could subsequently minimize methane production  while optimizing beef production. In this study, a crossbred population was used to evaluate the relationship between potential genes identified from fine mapping and RFI. The descriptive statistics of feed efficiency, performance and carcass traits are given in Table 1. Using the single trait animal model in ASReml, the estimates of heritabilities are given in Table 2. Heritability estimates for feed efficiency traits are in the range reported in the literature. The estimated heritability for RFI (0.19) is within the reported range from 0.16 to 0.45 [28, 29], whereas heritability (0.35) for ADG is in agreement with . Estimated heritability for DMI (0.42) is within the reported range, from 0.31 to 0.44 [30, 31]. FCR (0.25) is also within the reported range (0.17 to 0.37) [28, 30] as is MMWT (0.48) from 0.36 to 0.69 [28, 29]. The genetic variation and moderate heritabilities in feed efficiency traits indicate effective selection would be possible, where the trait is measured. However, the detected genetic variation also indicates that MAS could be effective where the genetic markers are closely linked to or is the causative mutation and that have repeatable effects across independent populations.
In the current study, SNP effects were estimated using an allele substitution effect model or the genotypic model. To avoid population stratification effects from influencing the estimated SNP effect, the phenotypes were adjusted for breed proportion, and the polygenic effect was fitted using the animal model to account for possible family effects . The population used in the GWAS using the Illumina BovineSNP50 was different animals to those used in the current study. However, the two populations are not independent as the animals in the current study were born subsequent to the animals used in the GWAS population from the same primary herds. The average relatedness among individuals between the two populations was estimated to be low at 0.0005 on average using the numerator relationship matrix calculated using CFC . There was zero pedigree-based inbreeding among the animals used in the current study. There was no separate dataset for feed efficiency traits on the same 339 SNP chip that could be used for validation. Thus, these significant associations require validation in other independent populations.
Results indicated 15 SNPs were significantly (at less 5% FDR) associated with at least one feed efficiency trait phenotype using the allele substitution effect model (Table 3). These findings reveal several candidate genes that provide highly significant evidence of association with RFI (Table 3). These promising candidate genes are located on Bos taurus autosomes (BTA) 8, 15, 16, 18, 19, 21, and 28. The strongest evidence of association with RFI and DMI was in SNP (8: 10674426) in the three prime untranslated region (3′ UTR) of gene elongation protein 3 homolog (ELP3). Gene ELP3 modulates transcription by working as a catalytic histone acetyltransferase subunit of the RNA polymerase II elongator complex involved in transcriptional elongation [34, 35]. In Drosophila, reduction in ELP3 expression during the development of the nervous system increases activity and decreases sleep  and the growth of adult flies (or could be lethal for the pupa) .
In the current study on BTA 16, the splice site intronic mutation (rs41820824) and the missense mutation (rs41821600) within gene hemicentin 1(HMCN1) were associated with RFI, where the substitution with the minor allele was associated with increased RFI and decreased F1. In addition, the minor allele of SNP rs41824268, within gene HMCN1, was associated with decreasing HCW, whereas in SNP rs41820800, it was associated with decreasing F2. Gene HMCN1 is known to be involved in age-related, macular degeneration , and polymorphisms within gene HMCN1 were associated with diabetes in man .
Synonymous coding SNP (18: 17150858), within the gene encoding zinc finger protein 423 (ZNF423), was associated with RFI, DMI and MMWT and located near a reported QTL (ID: 4449) for DMI . In addition, the minor allele of SNP (18: 17152044) was associated with decreasing GRF, F3, and F2, and increased LY, and was located near a reported QTL (ID: 11062) for LMA and body weight (ID: 11061) . Gene ZNF423 is a transcription factor involved in metal ion-binding. Down regulation of ZNF423 expression increases cell growth and retards differentiation as a consequence of its important role with the Vitamin A metabolite, retinoic acid .
On BTA 6, the SNP (6: 37288379) at 3′ UTR, within gene protein phosphatase, Mg2+/Mn2+ dependent, 1 K (PPM1K), was associated with increased MMWT and HCW and decreased RFI, FCR, marbling, and IFR and was located near a reported QTL (ID: 10761) for fat thickness at the 12th rib and a QTL (ID: 10758) for marbling score (EBV)  and a QTL (ID: 1753) for milk fat percentage . Gene PPM1K is involved in the phosphorus metabolic process or in amino acid dephosphorylation. In addition, PPM1K plays a key role in cellular survival and development by regulating mitochondrial permeability transition pore function . However, different genes are involved in mitochondrial adenosine triphosphate (ATP) synthesis efficiency and associated with differences in RFI [44–46], therefore, the effect of gene PPM1K on mitochondrial ATP synthesis is not clear .
For ADG, the most significant (at less than 5% FDR; P = 0.0009) SNP (rs41574929) was located on BTA 6, at 5′ UTR, within gene family with sequence similarity 190, member (A FAM190A; ID: 616908) (Table 3). The SNP rs41574929 was also associated ssignificantly at less than 5% FDR with HCW (P = 0.006). This result is in agreement with the function described for FAM190A where it is a necessary regulator for normal mitosis . A deletion mutation in FAM190A causes a cell division defect .
Allele substitution effect estimates of SNPs influencing (P ≤ 0.05) growth and efficiency traits, but which did not pass chromosome wise false discovery rate (FDR) threshold q = 0.2 were listed in Additional file 2. Also, all the SNPs associated at P-value < 0.05 using the genotypic model for growth and feed efficiency traits were listed in Additional file 3.
The association analysis using SLRM indicated that 59 SNPs were strongly (Table 4) or suggestively (Additional file 4) associated at 5% or 20% FDR test, respectively, for at least one carcass trait phenotype. Results indicated that the majority of the strong or suggestive associations were for intermuscular fat % (IFR) (14 indications). Thirteen SNPs were strongly associated with marbling, whereas 12 SNPs were associated with longissimus muscle area (LMA), and as follows HCW (9), F3 (8), GRF (8), F2 (6), LY (6), body cavity fat within the rib section (BFR) (5), % lean meat within the rib section (LR) (4), % subcutaneous fat within the rib section (SQFR) (4), and F1 (2).
Significant effects (at less 5% FDR) were found in 27 genes where gene ERCC5 (ID: 509602) had the highest proportion of the effects, revealing 8 of the significant associations with carcass traits (Table 4). The newly discovered SNP on BTA 12 (76889667 bp), within gene ERCC5 (ID: 509602), provided evidence of association with 5 of the studied carcass traits where the substitution of the minor allele was associated with increases of LYR and decreases in marbling, F3, GRF, IFR and SQFR (Table 4). Another newly discovered SNP in gene ERCC5 (76885563 bp) was strongly associated with three carcass traits where the substitution of the minor allele was associated with increases in F3, and GRF and decreases in LR (Table 4). Gene ERCC5 is involved in response to abiotic stimulus and negative regulation of programmed cell death and nucleotide excision repair pathway. In mice selected for high muscle mass, ERCC5 was located in QTL for lean mass .
SNP on BTA 27 (39712547 bp), within gene solute carrier family 20 (phosphate transporter), member 2 (SLC20A2; ID: 518905), was associated with one carcass trait where the substitution of the minor allele was associated with an increase in BFR (Table 4). The SLC20A2 is involved in ion and cation transport. In human, mutations within SLC20A2 are associated with idiopathic basal ganglia calcification .
SNP rs43702346 on BTA 6, within gene polycystic kidney disease 2 (PKD2; ID: 530393), was significantly associated with two carcass traits where substitution with the minor allele was associated with a decrease in HCW and IFR (Table 4). The PKD2 gene is involved in negative regulation of G1/S transition of mitotic cell cycle process. Gene PKD2 is near an identified QTL for bone percentage, fat percentage, meat percentage, meat-to-bone ratio, moisture content and subcutaneous fat . In human, polymorphisms within PKD2 may take part in the development of gout .
The in silicofunctional analysis
In the current study, the 74 genes containing significant (P < 0.05) SNPs were submitted to DAVID for enrichment analysis. In total 39 genes out of the 74 genes were enriched in 35 biological process terms (Table 5). Ion transport and cation transport mechanisms contained the highest number of genes associated with feed efficiency traits. In addition, some genes affecting feed efficiency traits in the current study were involved in proteolysis, protein complex biogenesis, and protein amino acid glycosylation. The ion transport mechanism in conjunction with protein turnover and metabolism account for 37% of the variation in RFI .
In ruminants protein synthesis accounts for 23% of total energy use in the whole body and protein turnover accounts for 42% of total gastrointestinal tract energy use . In the current study, some genes were involved in phosphorus metabolic processes, phosphorylation, and amino acid phosphorylation. Protein metabolism can be controlled by changing the phosphorylation status . Genes involved in phosphorus metabolic processes and phosphorylation mechanisms regulate the metabolism of energy . In the current study, regulation of transcription mechanisms contributed to variation in feed efficiency traits. The connection between a functional mutation in a specific transcription factor can increase or decrease expression of genes involved in glucose, amino acid, lipid, and cholesterol metabolism . Other studies have demonstrated that genes that up-regulate in response to nutritional restriction are involved in transcription control .
The in silico functional study of genes having significant SNPs revealed potential pathways likely to contribute to variation in feed efficiency traits (Table 6). Mitogen-activated protein kinases (MAPK) signaling pathway included three of the identified genes (RASA1, CACNA1G and STK3). In a study of the differences in global gene expression between high and low RFI animals, the majority of up-regulated genes in low RFI animals were stimulated by MAPKs , where the MAPKs were involved in signal transduction pathways to activate different cellular processes, such as cell division, differentiation, and cell death as a response to hormones and stress . The TYR gene is involved in Riboflavin metabolism, melanogenesis, tyrosine metabolism, and catecholamine biosynthesis, and the minor allele of SNP rs42402428, within gene TYR (ID: 280951) was associated with decreasing FCR. Polymorphisms in gene TYR have been associated with changing the coat colour of Braunvieh cattle . Gene GALNT13, affecting ADG, MMWT, DMI, F2, GRF, HCW, LMA, LY, and F3, is involved in mucin type O-Glycan biosynthesis. Gene ATP6V1E2 (ID: 540113), which affects DMI and MMWT, plays an important role in various pathways and biological mechanisms. Gene ATP6V1E2 is near an identified QTL for mycobacterium avium spp. Paratuberculosis resistance in Holstein cattle . Gene GTF2F2 (ID: 509259) affected RFI and is involved in basal transcription factors pathways, which regulate glucose, amino acids and protein, lipid metabolism and many other important metabolic processes. Changes in the function of GTF2F2 would be associated with feed efficiency or metabolic diseases . The minor allele of a newly discovered SNP (6: 37386084), within gene ABCG2 (ID: 536203), was associated with decreasing IFR and marbling. The in silico functional analysis showed that gene ABCG2 is involved in ATP-binding cassette (ABC) transporters and bile secretion pathways. The results of gene ABCG2 in the current study agree with reported gene ABCG2 as QTL for increasing milk yield and decreasing milk fat and protein [62–64]. The analysis also indicated that insulin-like growth factor 1 receptor gene IGF1R (ID: 281848) affecting ADG and marbling is involved in seven different pathways. Nonetheless, there was no association between production traits and the genotypes of IGF-IR/TaqI polymorphism [65–67]. This might be because a small number of animals was used to test the association in those analyses. Functional analysis allows a better understanding of the underlying mechanisms contributing to the genetic variation in feed efficiency, and it sheds light on potential pathways to target in future investigations.
Genetic variation in RFI explained by candidate genes
The accuracy of a DNA panel to predict a trait like feed efficiency depends on the amount of genetic variation explained. The 98 SNP set associated (P < 0.05) with at least one feed efficiency trait included SNPs that did not pass the FDR threshold, although they significantly contributed towards building the prediction equation in GWAS. The 98 SNP set explained 26% of the genetic variance in RFI whereas the proportion explained by the set of 339 SNPs was 29.6%. The correlation between EBVs of RFI using ASReml and GEBV were 0.52 and 0.66 from the 98 and 399 SNP sets, respectively. Based on the proportion of the genetic variance explained by the 98 SNPs (26%), the corresponding Beef Improvement Federation (BIF) accuracy is 0.127. Nonetheless, the estimated genetic variance by the 98 SNPs might be overestimated as the additive polygenic animal effect was not included in the model. To improve the accuracy of the SNP panel developed from a crossbred population, a large number of phenotypes is required (~2000 animals) . This might partially explain the relatively low estimated accuracy in the current study. In addition, large numbers of identified genes (83 out of 180) from fine mapping RFI were genotyped for only one SNP, and that decreases the probability of detecting the functional mutations. Nonetheless, combining validated SNPs from further fine mapping and the identified 98 SNPs may help develop a DNA test panel for commercial use.
This study reported SNPs that are significantly associated with RFI, performance, and carcass traits. We postulated that the identified significant SNPs, genes, biological mechanisms and pathways could be the direct cause of the variations in feed efficiency traits and carcass traits. The ability of the significant SNP to predict the genetic merit of feed efficiency and carcass traits should be measured in another population.
Single nucleotide polymorphism
Quantitative trait loci
False discovery rate
Bos taurus autosome
Residual feed intake
Average daily gain
Average daily dry matter intake, MMWT, Mid-point metabolic weight
Feed efficiency conversion ratio (kg gain kg-1 DM)
Subcutaneous fat depth between the 1st and 2nd quarter of the longissimus
Subcutaneous fat depth between 2nd and 3rd quarter of the longissimus
Subcutaneous fat depth between the 3rd and 4th quarter of the longissimus
Hot carcass weight (kg)
Longissimus dorsi muscle area (cm2)
Lean meat within the rib section (%)
Lean yield grade (%)
Grade fat (mm)
Intermuscular fat (%)
Body cavity fat within the rib section (%)
Proportion of subcutaneous fat from the rib section (%)
Elora Beef Research Center
the Agriculture and Agri-Food Canada Kapuskasing Research Centre
New Liskeard Agriculture Research Station.
Wulfhorst JD, Ahola JK, Kane SL, Keenan LD, Hill RA: Factors affecting beef cattle producer perspectives on feed efficiency. J Anim Sci. 2010, 88 (11): 3749-3758. 10.2527/jas.2010-2907.
Meuwissen TH, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001, 157 (4): 1819-1829.
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connell J, Moore SS, Smith TP, Sonstegard TS, et al: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4 (4): e5350-10.1371/journal.pone.0005350.
Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME: Invited review: genomic selection in dairy cattle: progress and challenges. J Dairy Sci. 2009, 92 (2): 433-443. 10.3168/jds.2008-1646.
Dekkers JC, Hospital F: The use of molecular genetics in the improvement of agricultural populations. Nat Rev Genet. 2002, 3 (1): 22-32. 10.1038/nrg701.
Nkrumah JD, Sherman EL, Li C, Marques E, Crews DH, Bartusiak R, Murdoch B, Wang Z, Basarab JA, Moore SS: Primary genome scan to identify putative quantitative trait loci for feedlot growth rate, feed intake, and feed efficiency of beef cattle. J Anim Sci. 2007, 85 (12): 3170-3181. 10.2527/jas.2007-0234.
Barendse W, Reverter A, Bunch RJ, Harrison BE, Barris W, Thomas MB: A validated whole-genome association study of efficient food conversion in cattle. Genetics. 2007, 176 (3): 1893-1905. 10.1534/genetics.107.072637.
Sherman EL, Nkrumah JD, Moore SS: Whole genome single nucleotide polymorphism associations with feed intake and feed efficiency in beef cattle. J Anim Sci. 2010, 88 (1): 16-22. 10.2527/jas.2008-1759.
Bolormaa S, Hayes BJ, Savin K, Hawken R, Barendse W, Arthur PF, Herd RM, Goddard ME: Genome-wide association studies for feedlot and growth traits in cattle. J Anim Sci. 2011, 89 (6): 1684-1697. 10.2527/jas.2010-3079.
Mujibi FD, Nkrumah JD, Durunna ON, Grant JR, Mah J, Wang Z, Basarab J, Plastow G, Crews DH, Moore SS: Associations of marker panel scores with feed intake and efficiency traits in beef cattle using preselected single nucleotide polymorphisms. J Anim Sci. 2011, 89 (11): 3362-3371. 10.2527/jas.2010-3362.
Snelling WM, Allan MF, Keele JW, Kuehn LA, Thallman RM, Bennett GL, Ferrell CL, Jenkins TG, Freetly HC, Nielsen MK, et al: Partial-genome evaluation of postweaning feed intake and efficiency of crossbred beef cattle. J Anim Sci. 2011, 89 (6): 1731-1741. 10.2527/jas.2010-3526.
Abo-Ismail MK, Squires EJ, Swanson KC, Lu D, Wang Z, Mah J, Plastow G, Moore SS, Miller SP: Fine mapping QTL and candidate genes discovery for residual feed intake on Chromosomes 5, 15, 16, and 19 in beef cattle. Proceedings of the 9th World Congress on Genetics Applied to Livestock Production: 1-6 August 2010; Leipzig. Edited by: German Society for Animal Science. 2010, 0783-
Pinheiro J, Bates D, DebRoy S, Sarkar D, Team. tRDC: nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–102. 2011
Arthur PF, Renand G, Krauss D: Genetic and phenotypic relationships among different measures of growth and feed efficiency in young Charolais bulls. Livest Prod Sci. 2001, 68 (2): 131-139.
Bergen R, Miller SP, Mandell IB, Robertson WM: Use of live ultrasound, weight and linear measurements to predict carcass composition of young beef bulls. Can J Anim Sci. 2005, 85 (1): 23-35. 10.4141/A04-011.
Sequencing TBG, Consortium A, Elsik CG, Tellam RL, Worley KC: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324 (5926): 522-528.
Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18 (11): 1851-1858. 10.1101/gr.078212.108.
Grant JR, Arantes AS, Liao X, Stothard P: In-depth annotation of SNPs arising from resequencing projects using NGS-SNP. Bioinformatics. 2011, 27 (16): 2300-2301. 10.1093/bioinformatics/btr372.
Aulchenko YS, Ripke S, Isaacs A, van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics (Oxford, England). 2007, 23 (10): 1294-1296. 10.1093/bioinformatics/btm108.
Gilmour AR, Gogel BJ, Cullis BR, Welham SJ, Thompson R: ASReml User Guide Release 3.0. 2009, Hemel Hempstead, HP1 1ES, UK: VSN International Ltd
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995, 57 (1): 289-300.
Fernando RL, Garrick D: GenSel—User manual for a portfolio of genomic. Selection related analyses. 2009, [http://taurus.ansci.iastate.edu]
Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78 (4): 629-644. 10.1086/502802.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38 (Database issue): D355-D360.
Davis GP, DeNise SK: The impact of genetic markers on selection. J Anim Sci. 1998, 76 (9): 2331-2339.
Nkrumah JD, Okine EK, Mathison GW, Schmid K, Li C, Basarab JA, Price MA, Wang Z, Moore SS: Relationships of feedlot feed efficiency, performance, and feeding behavior with metabolic rate, methane production, and energy partitioning in beef cattle. J Anim Sci. 2006, 84 (1): 145-153.
Herd RM, Bishop SC: Genetic variation in residual feed intake and its association with other production traits in British Hereford cattle. Livest Prod Sci. 2000, 63 (2): 111-119. 10.1016/S0301-6226(99)00122-0.
Crowley JJ, McGee M, Kenny DA, Crews DH, Evans RD, Berry DP: Phenotypic and genetic parameters for different measures of feed efficiency in different breeds of Irish performance-tested beef bulls. J Anim Sci. 2010, 88 (3): 885-894. 10.2527/jas.2009-1852.
Schenkel FS, Miller SP, Wilton JW: Genetic parameters and breed differences for feed efficiency, growth, and body composition traits of young beef bulls. Can J Anim Sci. 2004, 84 (2): 177-185. 10.4141/A03-085.
Arthur PF, Archer JA, Johnston DJ, Herd RM, Richardson EC, Parnell PF: Genetic and phenotypic variance and covariance components for feed intake, feed efficiency, and other postweaning traits in Angus cattle. J Anim Sci. 2001, 79 (11): 2805-2811.
Goddard ME, Hayes BJ: Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009, 10 (6): 381-391. 10.1038/nrg2575.
Sargolzaei M, Iwaisaki H, Colleau JJ: CFC: A tool for monitoring genetic diversity. The 8th World Congress on Genetics Applied to Livestock Production: Aug. 13-18, 2006; Belo Horizonte, Brazil. Edited by: Instituto Prociência. 2006, 27-28.
Wittschieben BO, Otero G, de Bizemont T, Fellows J, Erdjument-Bromage H, Ohba R, Li Y, Allis CD, Tempst P, Svejstrup JQ: A novel histone acetyltransferase is an integral subunit of elongating RNA polymerase II holoenzyme. Mol Cell. 1999, 4 (1): 123-128. 10.1016/S1097-2765(00)80194-X.
Karam CS, Kellner WA, Takenaka N, Clemmons AW, Corces VG: 14-3-3 mediates histone cross-talk during transcription elongation in Drosophila. PLoS Genet. 2010, 6 (6): e1000975-10.1371/journal.pgen.1000975.
Singh N, Lorbeck MT, Zervos A, Zimmerman J, Elefant F: The histone acetyltransferase Elp3 plays in active role in the control of synaptic bouton expansion and sleep in Drosophila. J Neurochem. 2010, 115 (2): 493-504. 10.1111/j.1471-4159.2010.06892.x.
Walker J, Kwon SY, Badenhorst P, East P, McNeill H, Svejstrup JQ: Role of elongator subunit Elp3 in Drosophila melanogaster larval development and immunity. Genetics. 2011, 187 (4): 1067-1075. 10.1534/genetics.110.123893.
Schmutz SM, Berryere TG, Ciobanu DC, Mileham AJ, Schmidtz BH, Fredholm M: A form of albinism in cattle is caused by a tyrosinase frameshift mutation. Mamm Genome. 2004, 15 (1): 62-67. 10.1007/s00335-002-2249-5.
Kim S, Abboud HE, Pahl MV, Tayek J, Snyder S, Tamkin J, Alcorn H, Ipp E, Nast CC, Elston RC, et al: Examination of association with candidate genes for diabetic nephropathy in a Mexican American population. Clin J Am Soc Nephrol. 2010, 5 (6): 1072-1078. 10.2215/CJN.06550909.
McClure MC, Morsci NS, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, et al: A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle. Anim Genet. 2010, 41 (6): 597-607. 10.1111/j.1365-2052.2010.02063.x.
Huang S, Laoukili J, Epping MT, Koster J, Holzel M, Westerman BA, Nijkamp W, Hata A, Asgharzadeh S, Seeger RC, et al: ZNF423 is critically required for retinoic acid-induced differentiation and is a marker of neuroblastoma outcome. Cancer Cell. 2009, 15 (4): 328-340. 10.1016/j.ccr.2009.02.023.
Ashwell MS, Heyen DW, Sonstegard TS, Van Tassell CP, Da Y, VanRaden PM, Ron M, Weller JI, Lewin HA: Detection of quantitative trait loci affecting milk production, health, and reproductive traits in Holstein cattle. J Dairy Sci. 2004, 87 (2): 468-475. 10.3168/jds.S0022-0302(04)73186-0.
Lu G, Ren S, Korge P, Choi J, Dong Y, Weiss J, Koehler C, Chen J-n, Wang Y: A novel mitochondrial matrix serine/threonine protein phosphatase regulates the mitochondria permeability transition pore and is essential for cellular survival and development. Genes Dev. 2007, 21 (7): 784-796. 10.1101/gad.1499107.
Kelly AK, Waters SM, McGee M, Fonseca RG, Carberry C, Kenny DA: mRNA expression of genes regulating oxidative phosphorylation in the muscle of beef cattle divergently ranked on residual feed intake. Physiol Genomics. 2011, 43 (1): 12-23. 10.1152/physiolgenomics.00213.2009.
Kelly DP, Scarpulla RC: Transcriptional regulatory circuits controlling mitochondrial biogenesis and function. Genes Dev. 2004, 18 (4): 357-368. 10.1101/gad.1177604.
Kolath WH, Kerley MS, Golden JW, Keisler DH: The relationship between mitochondrial function and residual feed intake in Angus steers. J Anim Sci. 2006, 84 (4): 861-865.
Patel K, Scrimieri F, Ghosh S, Zhong J, Kim M-S, Ren YR, Morgan RA, Iacobuzio-Donahue CA, Pandey A, Kern SE: FAM190A deficiency creates a cell division defect. Am J Pathol. 2013, 183 (1): 296-303. 10.1016/j.ajpath.2013.03.020.
Kärst S, Cheng R, Schmitt AO, Yang H, de Villena FPM, Palmer AA, Brockmann GA: Genetic determinants for intramuscular fat content and water-holding capacity in mice selected for high muscle mass. Mamm Genome. 2011, 22 (9–10): 530-543.
Zhang Y, Guo X, Wu A: Association between a Novel Mutation in < italic > SLC20A2</italic > and Familial Idiopathic Basal Ganglia Calcification. PLoS One. 2013, 8 (2): e57060-10.1371/journal.pone.0057060.
Gutiérrez-Gil B, Williams JL, Homer D, Burton D, Haley CS, Wiener P: Search for quantitative trait loci affecting growth and carcass traits in a cross population of beef and dairy cattle. J Anim Sci. 2009, 87 (1): 24-36.
Lee YH, Song GG: Pathway analysis of genome-wide association studies on uric acid concentrations. Hum Immunol. 2012, 73 (8): 805-810. 10.1016/j.humimm.2012.05.004.
Richardson EC, Herd RM: Biological basis for variation in residual feed intake in beef cattle. 2. Synthesis of results following divergent selection. Aust J Exp Agric. 2004, 44 (5): 431-440. 10.1071/EA02221.
Caton JS, Bauer ML, Hidari H: Metabolic components of energy expenditure in growing beef cattle - review. Asian-Australasian J Animal Sci. 2000, 13 (5): 702-710.
Baldwin RL, Sainz RD: Energy partitioning and modeling in animal nutrition. Annu Rev Nutr. 1995, 15: 191-211. 10.1146/annurev.nu.15.070195.001203.
Berndt T, Kumar R: Novel mechanisms in the regulation of phosphorus homeostasis. Physiology (Bethesda, Md). 2009, 24: 17-25. 10.1152/physiol.00034.2008.
Brown DR, DeNise SK, McDaniel RG: Mitochondrial respiratory metabolism and performance of cattle. J Anim Sci. 1988, 66 (6): 1347-1354.
Desvergne B, Michalik L, Wahli W: Transcriptional Regulation of Metabolism. Physiological Reviews. 2006, 86 (2): 465-514. 10.1152/physrev.00025.2005.
Byrne KA, Wang YH, Lehnert SA, Harper GS, McWilliam SM, Bruce HL, Reverter A: Gene expression profiling of muscle tissue in Brahman steers during nutritional restriction. J Anim Sci. 2005, 83 (1): 1-12.
Chen Y, Gondro C, Quinn K, Herd RM, Parnell PF, Vanselow B: Global gene expression profiling reveals genes expressed differentially in cattle with high and low residual feed intake. Anim Genet. 2011, 42 (5): 475-490. 10.1111/j.1365-2052.2011.02182.x.
Pearson G, Robinson F, Beers Gibson T, Xu BE, Karandikar M, Berman K, Cobb MH: Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions. Endocr Rev. 2001, 22 (2): 153-183.
Pant SD, Schenkel FS, Verschoor CP, You Q, Kelton DF, Moore SS, Karrow NA: A principal component regression based genome wide analysis approach reveals the presence of a novel QTL on BTA7 for MAP resistance in holstein cattle. Genomics. 2010, 95 (3): 176-182. 10.1016/j.ygeno.2010.01.001.
Cohen-Zinder M, Seroussi E, Larkin DM, Loor JJ, Wind AE-v, Lee J-H, Drackley JK, Band MR, Hernandez AG, Shani M, et al: Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle. Genome Res. 2005, 15 (7): 936-944. 10.1101/gr.3806705.
Olsen HG, Nilsen H, Hayes B, Berg PR, Svendsen M, Lien S, Meuwissen T: Genetic support for a quantitative trait nucleotide in the ABCG2 gene affecting milk composition of dairy cattle. BMC Genet. 2007, 8: 32-
Seroussi E: The concordance test emerges as a powerful tool for identifying quantitative trait nucleotides: lessons from BTA6 milk yield QTL. Anim Genet. 2009, 40 (2): 230-234. 10.1111/j.1365-2052.2008.01809.x.
Akis I, Oztabak K, Gonulalp I, Mengi A, Un C: IGF-1 and IGF-1r gene polymorphisms in East Anatolian Red and South Anatolian Red cattle breeds. Genetika. 2010, 46 (4): 497-501.
Curi RA, Krauskopf MM, Hadlich JC, Fortes MRS, Vankan DM, Silva JAIIV, Oliveira HN, Mota MDS, XUhwsbspssa, pid SX, nrm = iso: Candidate SNPs for carcass and meat traits in Nelore animals and in their crosses with Bos taurus. Pesq Agrop Brasileira. 2012, 47: 294-301. 10.1590/S0100-204X2012000200019.
Zhang R, Li X: Association between IGF-IR, m-calpain and UCP-3 gene polymorphisms and growth traits in Nanyang cattle. Mol Biol Rep. 2011, 38 (3): 2179-2184. 10.1007/s11033-010-0346-1.
Goddard ME, Hayes BJ, Meuwissen T: Genomic selection in farm animal species - lessons learnt and future perspectives. Proceedings of the 9th World Congress on Genetics Applied to Livestock Production: 1-6 August 2010; Leipzig. Edited by: German Society for Animal Science. 2010, 701-
The authors wish to thank: The staff at the Elora Beef Research Centre and the University of Guelph Meat Science laboratory and Tim Caldwell for their assistance with data collection and analysis. Financial support was provided by: Ontario Cattlemen’s Association, Agriculture Adaptation Council (FIP), Agriculture and Agri-food Canada (Growing Forward), Canadian Cattlemen’s Beef Cattle Research Council, Ontario Ministry of Agriculture Food and Rural Affairs, Alberta Beef Producers and Alberta Livestock and Meat Agency.
The authors declare that they have no competing interests.
MKA contributed in designing the study, preparing the phenotypes and genotypes, performing the statistical and enrichment analysis, and drafting the manuscript. RV performed the data editing and statistical analysis. GV and IBM provided help in collection of data, analysis and manuscript editing. EJS and KCS participated in designing the study, collection of data and manuscript editing. PS, SM, GP participated in designing the study, preparing the genotypes, and editing the manuscript. SPM helped in design, data collection, analysis, and draft the manuscript. All authors’ read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Abo-Ismail, M.K., Vander Voort, G., Squires, J.J. et al. Single nucleotide polymorphisms for feed efficiency and performance in crossbred beef cattle. BMC Genet 15, 14 (2014). https://doi.org/10.1186/1471-2156-15-14
- Candidate genes
- Single nucleotide polymorphism
- Feed efficiency
- Carcass traits