Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Research article

Association, effects and validation of polymorphisms within the NCAPG - LCORL locus located on BTA6 with feed intake, gain, meat and carcass traits in beef cattle

Amanda K Lindholm-Perry1*, Andrea K Sexten12, Larry A Kuehn1, Timothy PL Smith1, D Andy King1, Steven D Shackelford1, Tommy L Wheeler1, Calvin L Ferrell13, Thomas G Jenkins13, Warren M Snelling1 and Harvey C Freetly1

Author affiliations

1 United States Department of Agriculture, Agricultural Research Service, U.S. Meat Animal Research Center, State Spur 18D, Clay Center, NE 68933, USA

2 Kansas State University, Manhattan, KS, USA

3 F&J Livestock Consulting, LLC., Saronville, NE, USA

For all author emails, please log on.

Citation and License

BMC Genetics 2011, 12:103  doi:10.1186/1471-2156-12-103


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/12/103


Received:15 September 2011
Accepted:14 December 2011
Published:14 December 2011

© 2011 Lindholm-Perry et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

In a previously reported genome-wide association study based on a high-density bovine SNP genotyping array, 8 SNP were nominally associated (P ≤ 0.003) with average daily gain (ADG) and 3 of these were also associated (P ≤ 0.002) with average daily feed intake (ADFI) in a population of crossbred beef cattle. The SNP were clustered in a 570 kb region around 38 Mb on the draft sequence of bovine chromosome 6 (BTA6), an interval containing several positional and functional candidate genes including the bovine LAP3, NCAPG, and LCORL genes. The goal of the present study was to develop and examine additional markers in this region to optimize the ability to distinguish favorable alleles, with potential to identify functional variation.

Results

Animals from the original study were genotyped for 47 SNP within or near the gene boundaries of the three candidate genes. Sixteen markers in the NCAPG-LCORL locus displayed significant association with both ADFI and ADG even after stringent correction for multiple testing (P ≤ 005). These markers were evaluated for their effects on meat and carcass traits. The alleles associated with higher ADFI and ADG were also associated with higher hot carcass weight (HCW) and ribeye area (REA), and lower adjusted fat thickness (AFT). A reduced set of markers was genotyped on a separate, crossbred population including genetic contributions from 14 beef cattle breeds. Two of the markers located within the LCORL gene locus remained significant for ADG (P ≤ 0.04).

Conclusions

Several markers within the NCAPG-LCORL locus were significantly associated with feed intake and body weight gain phenotypes. These markers were also associated with HCW, REA and AFT suggesting that they are involved with lean growth and reduced fat deposition. Additionally, the two markers significant for ADG in the validation population of animals may be more robust for the prediction of ADG and possibly the correlated trait ADFI, across multiple breeds and populations of cattle.

Background

Feed intake and average daily gain are economically important phenotypes in beef cattle. The ability to genetically select for animals that better suit a producer's needs in terms of the amount of feed an animal consumes and the gain that can be expected would benefit the industry. However, it is critical that markers used to select for feed and growth traits do not adversely affect other traits including meat quality, carcass traits or reproduction.

Quantitative trait loci (QTL) for cattle growth traits such as birth weight, body length, carcass weight, and longissimus muscle area (LMA) have been detected on bovine chromosome 6 in several different populations of animals. Nkrumah et al. [1] discovered a QTL at 42 cM for average daily gain (ADG) in a composite population of steers. Body length at birth and birth weight QTL were identified by Gutiérrez-Gil et al. [2] at 38 and 39 Mb, respectively, in a Charolais × Holstein crossbred population. Several other QTL have been detected for birth weight in this same region [3-6]. In addition to live animal growth, length and weight QTL, others have found that QTL for carcass and meat quality traits are located in this region. Takasuga et al. [7] detected a QTL for carcass weight at 38 cM, Setogutchi et al. [8] found QTL for carcass weight and LMA at 39.1 cM and McClure et al. [6] identified a QTL for ribeye area (REA) at 41 cM. There is clear evidence that the region on BTA6 between 30 and 42 cM is involved in growth and size phenotypes across multiple breeds of cattle.

We have also detected a region on chromosome 6 that is correlated with feed intake and gain located within 37.96 - 38.53 Mb in a crossbred population of cattle [9]. There are seven genes on the Btau 4.0 genome assembly in this region, including leucine aminopeptidase 3 (LAP3), mediator complex subunit 28 (MED28), LOC523874, C6H4Orf30, non-SMC condensin I complex, subunit G (NCAPG), and LOC540095 (bovine ortholog of ligand dependent nuclear receptor corepressor-like, LCORL). LAP3 is an aminopeptidase that catalyzes the removal of amino acids from intracellular proteins and peptides. LAP3 is most highly expressed in the bovine pineal gland, kidney, skin, intestine, mammary, and adipose tissues (UGID: 1959690, UniGene Bt.56962). MED28 is a subunit of the Mediator complex that interacts with RNA polymerase II and can interact with co-activators and co-repressors [10]. MED28 is most highly expressed in the bovine pineal gland, lyphoreticular tissue, uterus, abomassum, and expressed in lower quantities in the liver, reticulum, intestine and kidney (UGID:1230902, UniGene: Bt.29898). LOC523874 is homologous to the human family with sequence similarity 184, member B (FAM184B) for which there is little information. Recent SNP association studies suggest that NCAPG has a function in cattle growth ([8,11] possibly through a role in cell proliferation [12]. NCAPG is expressed in bovine pancreas, skin, intestine, liver, kidney, and ovary (UGID:3093465, UniGene Bt.100379). LCORL has recently been associated with human skeletal size [13] and height growth in infancy [14]. LCORL is also expressed in the liver, intestine in the bovine fetus and calf, and in the extraembryonic tissue of the bovine fetus (UGID:1484658, UniGene Bt.38533). Based on this information, we chose to evaluate LAP3, NCAPG and LCORL as potential positional candidate genes for cattle feed intake and growth.

The purpose of this study was to further refine the region on chromosome 6 that was previously identified by a SNP association as significant for feed intake and average daily gain in our population of steers [9]. While genetic markers predictive for feed intake and growth may be useful for the identification and selection of animals that are more efficient, attention must be given to the impact of these SNP on other production traits. Accordingly, we evaluated the markers that were highly significant for feed intake and gain for effects on carcass and meat quality traits. Here, we present a number of markers that may be useful for the genetic selection of animals of multiple breeds for feed and gain traits and illustrate that these markers appear to play a role in lean growth and fat deposition.

Methods

Animals

All animal procedures were reviewed and approved by the U.S. Meat Animal Research Center (USMARC) Animal Care and Use Committee. Procedures for handling cattle complied with those specified in the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching [15].

Discovery population

The breeding scheme used for the SNP discovery population of animals with feed efficiency phenotypes is described in Snelling et al. [9]. Briefly, purebred Angus (AN), Hereford (HH), Simmental (SM), Limousin (LM), Charolais (CH), Gelbvieh (GV) and Red Angus (AR) sires were mated by artificial insemination to composite MARC III (1/4 AN, 1/4 HH, 1/4 Pinzgauer, 1/4 Red Poll), AN and HH cows to produce progeny designated as F1, born in 1999, 2000, and 2001. Female F1 and the 2001-born F1 males (AN, HH dams only) were kept for breeding, and mated in multiple-sire pastures to produce 2-, 3-and 4-breed cross progeny designated F12. The F12 calves were born in March through May of 2003 to 2007, from 3-year-old and older dams. Male calves were castrated within 24 h after birth. Calves were weaned in September at approximately 165 d of age [9]. Animals were screened and excluded for medical or health issues that may have affected either feed intake or gain phenotypes.

SNP identification population

Animals (n = 24) used for resequencing of candidate genes LAP3, NCAPG, and LCORL were selected from the discovery population of animals. These animals represented extreme feed intake and gain phenotypes.

Validation population

Validation animals (n = 406) were derived from progeny in the USMARC continuous Germplasm Evaluation Program [16]. These progeny were part of a breeding program to develop purebred cattle of each of 14 breeds. In addition to the seven mentioned previously, these include Braunvieh, Shorthorn, Brahman, Maine Anjou, Santa Gertrudis, Salers, and Chi-Angus. Sires of each of these 14 breeds were mated to female descendants of cattle in the discovery population. This population did not include full or partial sibs from the discovery population. Approximately 60% of the validation animals were half-bloods sired by these seven new breeds and 40% were half- or three-quarter-bloods from the original seven breeds. Progeny were produced in three breeding seasons (Fall 2007, Spring 2008, and Fall 2008). Animals (n = 29) were removed from the analysis due to chronic diagnoses of pneumonia, foot-rot, and bloat).

Feed efficiency phenotypes

Discovery population calves were familiarized to the test ration as described in Snelling et al. [9]. Briefly, beginning at 276 ± 15 d of age, individual feed intake measurements were acquired with Calan Broadbent Feeding Systems (American-Calan-Broadbent, Northwood, NH). Animals were provided feed at 0800 h in the morning and were given ad libitum access. The mixed ration on a dry matter basis was 82.0% dry rolled corn, 10.6% ground alfalfa hay, 5.66% Soybean meal, 1.25% limestone, 0.40% urea, 0.06% salt, 0.007% trace minerals, 0.008% Vitamins A, D, and E, and 0.015% Rumensin 80. Feed refusals were obtained each week.

Feed efficiency phenotypes for F12 steers (discovery population) are described in detail in Snelling et al. [9]. Briefly, dry matter intake (DMI) was equal to cumulative dry matter intake for the 140-d feeding period. Average daily feed intake (ADFI) is the average daily dry matter intake (DMI/140). Individual animal quadratic regressions were fitted for body weight on time, and gain was calculated as the difference of BW predicted at 140 d and the intercept [17].

The steers used for the validation study were fed for either 84 or 148 d. Feed intake for those fed for 84 d was measured using an Insentec system (Marknesse, The Netherlands). Diets were similar to those used in the discovery population. Dry matter intake was equal to the total DMI consumed for either the 84- or 148-d period and averaged to a daily basis. A linear regression was used to predict BW in steers fed for 84 d and a quadratic equation was used to predict BW in steers fed for 148 d. Gain for both groups was calculated for the respective feeding periods as described in the discovery population above.

Meat quality and carcass phenotypes (discovery population)

F12 steers were serially harvested at a commercial processing plant as described by King et al. [18]. USDA yield and quality grade data were obtained by trained USMARC personnel after 36 h in the cooler at 0°C. Wholesale ribs were obtained as described by King et al. [18]. A posterior section of the ribeye was frozen at 14 d postmortem and a 2.54 cm thick steak was cut from the 11th rib region. For slice shear force (SSF), steaks were thawed at 5°C for 24 h then cooked on a conveyorized electric belt grill to a final internal temperature of 71°C as described by Wheeler et al. [19]. SSF was determined as described by Shackelford et al. [20].

Identification of SNP markers in LAP3, NCAPG and LCORL (SNP identification population)

DNA was isolated from whole blood or buffy coats collected from steers selected from the discovery population of animals. Animals (n = 24) with extreme ADG or ADFI phenotypes were chosen for SNP discovery. Twelve animals with low ADG or ADFI values and twelve animals with high ADG or ADFI values were used. All seven breeds were represented in the animals with high and low phenotypes. None of the breeds were over-represented in the high and low groups. Angus and Hereford breeds each represented ~30% of the breed composition in both groups. Simmental represented 12.5% in the low group and 16% in the high. Composition of the other breeds (Limousin, Gelbvieh, Red Angus, Charolais, MARCIII) ranged from 3-8%.

Primer pairs for amplification of regions of LAP3, NCAPG and LCORL from genomic DNA were designed using Primer 3 [21] (code available at http://frodo.wi.mit.edu webcite). Cattle sequences used as templates for primer design were obtained from the chromosome 6 Btau 4.0 assembly. Oligonucleotide primers were synthesized by IDT (Integrated DNA Technologies, Coralville, IA). PCR was performed in a DNA engine Dyad® peltier thermal cycler (Bio-Rad, Hercules, CA). PCR reactions included 0.25 U Hot Star Taq polymerase (Qiagen, Valencia, CA); 1× supplied buffer; 1.5 mM MgCl2; 80 μM dNTPs; 0.33 μM each primer; and 25 ng genomic DNA in 12 μl reactions. Three microliters of the PCR reaction were electrophoresed in 2% agarose gels to determine quality of amplification, and the remainder was prepared for sequencing after treatment with 0.1 U exonuclease I (USB, Cleveland, OH). Sequencing reactions were precipitated with 70% isopropanol and sequenced with an ABI 3730 capillary sequencer (Applied Biosystems, Foster City, CA). Bases were called with Phred and assembled into contigs with Phrap. Polymorphisms were identified using Polyphred and assessed using Consed http://www.phrap.org webcite.

Genotyping analysis

Multiplex assays for the Sequenom MASSARRAY® instrument were designed with the MASSARRAY® Assay Design 3.0 software (Sequenom, San Diego, CA). Assays were designed for 47 SNPs in and around the LAP3, NCAPG, and LOC540095 (LCORL) gene loci. Each amplification primer contained a 10-base tag to yield amplification products with different masses. Amplicon lengths were between 90 and 120 bp. Reaction conditions were performed as recommended by Sequenom. A total of 1,137 animals with dry matter intake and average daily gain phenotypes and their sires were genotyped and used in the analyses after exclusion.

Statistical analysis

The data were analyzed using an animal model (MTDFREML) [22] that included the fixed effects of year and barn. Covariates of age and heterosis were also included. Calf and dam breed compositions were modeled with covariates for proportions AN, HH, SM, LI, CH, GV, AR, and MARC III. Covariates for expected calf heterosis were computed from parental breed composition. The number of copies of the minor allele was fitted separately in the association model. Polygenic and breed effects were included to reduce the effects of family structure on breed- and family-specific alleles [23,24]. Variance components for polygenic effects and error were estimated using MTDFREML [22]. Nominal significance values for each SNP were computed. A Bonferroni correction for multiple testing was applied to the P-values by multiplying the P-value by the number of markers tested (n = 47). Analyses of significant SNP in the validation population followed the same model definition.

Linkage disequilibrium (LD) analysis

Linkage disequilibrium (r2) was defined for the 47 SNP on chromosome 6 using Haploview 4.0 software [25]http://www.broad.mit.edu/mpg/haploview/index.php webcite. Blocks of LD were based on pairwise LD values. Haploview settings were as follows: the exclusion of animals with > 50% missing genotypes, ignoring pairwise comparisons of markers > 500 kb apart, the percentage of genotypes ≥50% and minimum minor allele frequency of 0.001.

Results

Candidate gene SNP identification

The genes LAP3, NCAPG and LCORL were chosen as positional and functional candidate genes in this region. The coding regions and some intronic regions of these genes were sequenced and polymorphisms are shown in Table 1. A polymorphism in the NCAPG gene identified previously by Setoguchi et al. [8] (p.I442M) was also detected in our population of animals (marker 83714_224 in this study). In addition to this polymorphism, we identified a second non-synonymous alteration in NCAPG (p.L877M). A 3 base pair deletion resulting in an in-frame deletion of a glutamic acid (p.E47del) and a non-synonymous alteration (p.D92E) were detected in LAP3 (Table 1).

Table 1. SNP identified in LAP3, LCORL, and NCAPG genes in crossbred steers.

Discovery population analyses

A total of 21 out of the 39 SNP identified within the LAP3 (n = 9), NCAPG (n = 8) and LCORL (n = 4) genes were multiplexed into Sequenom genotyping assays. In addition to these, another 26 publicly available SNP located up- and down-stream of these genes that were not used on the BovineSNP50 BeadChip [26-28], were also placed into the genotyping assays. The population of steers genotyped previously [9] with the BovineSNP50 BeadChip (n = 1,066) and their sires (n = 71) were genotyped with these markers (Table 2).

Table 2. SNP marker associations and estimated effects for ADFI and ADG in the discovery population of crossbred steers.

Marker associations for ADFI and ADG are presented in Table 2. A Bonferroni correction for multiple testing was applied to the P-values obtained from the genotyping analysis (Table 2). P-values were multiplied by the number of markers tested (n = 47) to obtain a corrected P. After correction, 16 markers between 38.16 and 38.37 Mb were significantly associated with ADFI (P < 0.05) and 20 markers between 38.16 and 38.41 Mb were associated with ADG (P < 0.05).

Two markers within the NCAPG gene locus were associated (after multiple correction) with ADFI (P ≤ 0.02) and seven were associated with ADG (P = 0.003; Table 2); one (83714_224) was the same marker previously identified by Setoguchi et al. [8] (p.I442M). However, several other markers were identified with lower P-values for ADFI and ADG (Table 2).

To determine whether these markers may adversely affect carcass and meat traits, they were analyzed for association with hot carcass weight (HCW), adjusted fat thickness (AFT), marbling score (MARB), ribeye area (REA), and SSF. The phenotypic data is presented in Table 3 and includes: number of animals, mean with standard deviation, and range of phenotypes. The effects and association between meat quality and carcass traits for the markers on BTA6 are presented in Table 4. The alleles for positive effects or increases in ADFI and ADG were associated with positive effects or higher HCW and REA. These same alleles had negative effects or lower AFT. Marbling and SSF were not significant for any of the markers tested (Table 4).

Table 3. Description of the USMARC resource population carcass and meat quality phenotypic means.

Table 4. Association and effects of SNP significant for ADG on BTA6 with meat quality and carcass traits in the discovery population of crossbred steers.

Validation population analysis

Twenty of the markers genotyped in the discovery population were genotyped in a validation population of crossbred animals (n = 406; Table 5). This set of markers included two markers in the LAP3 gene locus, 6 in NCAPG and 9 in the LCORL gene loci. Of the 20 markers, 13 and 17 were associated with ADFI and ADG, respectively, after Bonferroni correction in the discovery population. In addition to these 20 markers, six markers from the Illumina BovineSNP50 BeadChip that were nominally significant ADG (4 were also nominally significant for ADFI) in the discovery population were also genotyped in the validation population of animals. Association analysis produced two markers (Hapmap28589-BTC-041045, P = 0.04 and 81405_282, P = 0.03) that were associated with ADG in these animals (Table 5). Both markers are located within the LCORL gene locus. Marker 83714_224 (or p.I422M) [8], in NCAPG was not significant in this population (P = 0.09). Genotypes were analyzed with Haploview for linkage disequilibrium (LD) in this population of animals (Figure 1). The LD between marker 83714_224 in NCAPG and 81405_282 in LCORL was lower at r2 = 0.89 than the LD between these two markers in the SNP discovery population (r2 = 0.96) (Figure 2).

Table 5. SNP marker associations and estimated effects for ADG in the validation population of crossbred steers (n = 406).

thumbnailFigure 1. Haploview display of LD of the 26 genetic markers on chromosome 6 that were genotyped on the validation population of crossbred steers (n = 406). Pairwise LD in r2 is shown between markers with shading. The darker shading represents higher LD, black diamonds are r2 = 1. Markers in bold signify those that are located in a haplotype block. Markers with *_ prior to the marker name are markers from the Illumina BovineSNP50 BeadChip that were associated with ADG in the discovery population of animals [9].

thumbnailFigure 2. Haploview display of LD of the genetic markers on chromosome 6 genotyped on the SNP discovery population of crossbred steers (n = 1,066). Pairwise LD in r2 is shown between markers by shading. The darker shading represents higher LD, black diamonds are r2 = 1. Markers in bold signify those that are located in a haplotype block.

The two markers significantly associated with ADG in the validation population, Hapmap28589-BTC-041045 and 81405_282, are each responsible for 1.6% of the additive phenotypic variation of ADG (data not shown); however, they are in complete LD (Figure 1). These same markers each account for 3.6% of the phenotypic variation of ADG in the discovery population. These markers also represent 2% of the ADFI phenotypic variation in the discovery population of animals. These markers are also in complete LD in the discovery population of animals (Figure 2).

Discussion

A previous single nucleotide polymorphism association analysis identified a 570 Kb region on BTA6 from 37.96 to 38.53 Mb containing eight SNP with significant association (P ≤ 0.003) with ADG and three SNP with significant association with ADFI (P ≤ 0.002) [9]. The three markers significant for ADFI were the same markers significant for ADG. The most significant marker in this region from Snelling et al. [9] for ADG was ARS-BFGL-NGS-45457 (P, nominal = 4.8 × 10-6) and the most significant marker for ADFI was Hapmap31285-BTC-041097 (P, nominal = 8 × 10-4). Of the markers presented in this study, 14 are more significant for ADFI (P, nominal = 5 × 10-4) and 11 were more significant for ADG (P, nominal = 4 × 10-6) than the markers on the BovineSNP50 BeadChip.

We chose to evaluate LAP3, NCAPG and LCORL as potential positional and functional candidate genes on bovine chromosome 6 for steer feed intake and gain. These genes were within the region that was identified as significant for ADFI and ADG from a previous association analysis [9]. None of the markers in LAP3 or in close proximity to LAP3 were significant after correction for multiple testing suggesting that these SNP are not in high LD with the causative mutation and that this region is not likely to be contributing to variation in feed intake and gain in this population of cattle. Our data do, however, suggest that the more likely region responsible for variation in cattle feed intake and gain is further downstream in the NCAPG-LCORL locus.

NCAPG is a subunit of the condensin 1 complex that is involved in the condensation of chromosomes and also interacts with a DNA methyltransferase linking methylation and chromatin condensation; however, the molecular and cellular events that might specifically tie NCAPG to growth and appetite have not been identified.

LCORL encodes a transcription factor that is thought to function during spermatogenesis. Aside from the association of SNPs in LCORL with human height phenotypes, there is little published information regarding this gene and its potential functions relating to growth. While it is possible to speculate that as a transcription factor, SNPs that affect the transcription or translation of LCORL may result in increased or decreased regulation of gene(s) involved in growth and/or appetite pathways, these functions have yet to be elucidated.

The NCAPG-LCORL locus has been interrogated in other populations of cattle for growth, meat quality and carcass traits. The most significant marker tested within the NCAPG-LCORL region in Japanese Black and Japanese Brown cattle for carcass weight, longissimus muscle area and subcutaneous fat thickness was the I422M marker (or marker 83714_224 in this study) in NCAPG [8]. This same marker was significant for fetal growth in a crossbred population of Charolais and German Holstein [11]. In our crossbred steer population consisting of seven continental breeds of beef cattle, this marker was significant for similar traits including, HCW, REA, AFT, and also for ADG and ADFI; however, markers located further downstream and within the LCORL gene locus were even more significant. The only markers tested in all three of these studies were in the NCAPG gene. Only one marker within the LCORL gene locus (HapMap28154-BTC-041134) was tested in our population and in the study by Eberlein et al. [11], making a direct comparison of the LCORL region difficult. In our population of cattle, we detected 14 SNP in the NCAPG-LCORL locus that were significant for ADFI, ADG and the carcass/meat quality traits HCW, AFT and REA. All of these SNP are in LD (r2 ≤ 0.73), thus it is unclear which of these SNP and genes may be responsible for the phenotypic variation within this population of animals.

In an attempt to better define which SNPs are responsible for the phenotypic variation, we performed a validation of these markers in another crossbred population of animals consisting of 14 breeds of cattle. This produced two SNP associated with ADG located in the intronic regions of the LCORL gene and suggests that these SNP may be robust among differing breeds of cattle for the prediction of weight gain. The LD in the validation population of animals in this region differs from the LD in the discovery population of animals. An example is that the LD for the NCAPG SNP 83714_224 (I422M) marker and the most significant SNP for ADG in our study (HapMap28589-BTC041045) was r2 = 0.96 in our SNP discovery population and r2 = 0.91 in the validation population.

The SNPs tested in our study and in previous work by Setoguchi et al. [8] are significant for traits that suggest they are involved in lean growth and fat deposition. In support of this, Gutiérrez-Gil et al. [2] also detected a QTL for subcutaneous fat in this region at BTA6:38.5, and McClure et al. [6] found a QTL for fat thickness at the 12th rib at BTA6:36. These QTL plus the numerous QTL for growth and carcass weight seem to be harboring a gene or genes that are responsible for lean growth and fat deposition. Moreover, the NCAPG I442M SNP was associated with plasma arginine, symmetric dimethylarginine (SDMA), and linoleylcarnitine (C18:2) levels in a recent metabolomic study [29]. The allele that was associated with positive weight gain during 182-273 d (442M) was also associated with higher plasma levels of arginine, SDMA, and linoleylcarnitine [29]. The role of arginine in growth is through its activation of the mTOR pathway in skeletal muscle which increases protein synthesis [29,30], and through its function as a precursor of nitrogen oxide which effects glucose uptake and insulin responsiveness in skeletal muscle [29,31,32]. These data suggest there is some type of physiological role for growth via arginine for the gene or genes located within this chromosomal region. As the gene functions for NCAPG and LCORL do not currently seem to suggest a direct interaction with the metabolites identified in Weikard et al. [29], and the SNP in these genes are in high LD with each other in our discovery and validation populations of animals, additional functional genomic studies will be required to elucidate the mechanisms of NCAPG and/or LCORL in growth.

Conclusion

In summary, we have identified several markers that are highly significant for steer feed intake and gain phentoypes. Moreover, several of the markers evaluated in this study are more significant for ADFI and ADG than the markers on the BovineSNP50 BeadChip. Genetic selection with these markers is likely to have an impact on HCW, REA and AFT with little influence on marbling or tenderness. It is not known how these markers will affect reproductive traits. These markers need to be evaluated on additional populations of animals to ensure that they maintain utility across breeds. We were able to confirm that two of these markers were significant for ADG in a validation population of crossbred cattle that included 14 breeds. These two markers may be particularly effective for the prediction of weight gain across various populations and breeds of cattle.

Authors' contributions

ALP conceived of the study, carried out molecular genetic studies, and drafted the manuscript. AKS participated in the project design and carried out molecular genetic studies. LAK performed statistical analysis of genotypic data. WMS provided statistical advice. TPLS assisted with SNP selection. DAK, SDS and TLW managed meat quality and carcass phenotypic data collection. CLF, TGJ, HCF and LAK collected and provided feed efficiency phenotypic data. All authors read and approved the manuscript.

Acknowledgements

The authors would like to thank Linda Flathman for excellent technical assistance and Donna Griess for assistance with manuscript formatting and editing. Mention of trade name proprietary product or specified equipment does not constitute a guarantee or warranty by the USDA and does not imply approval to the exclusion of other products that may be suitable. USDA is an equal opportunity provider and employer.

References

  1. 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:3170-3181. PubMed Abstract | Publisher Full Text OpenURL

  2. 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:24-36. PubMed Abstract | Publisher Full Text OpenURL

  3. Casas E, Shackelford SD, Keele JW, Stone RT, Kappes SM, Koohmaraie M: Quantitative trait loci affecting growth and carcass composition of cattle segregating alternate forms of myostatin.

    J Anim Sci 2000, 78:560-569. PubMed Abstract | Publisher Full Text OpenURL

  4. Kneeland J, Li C, Basarab J, Snelling WM, Benkel B, Murdoch B, Hansen C, Moore SS: Identification and fine mapping of quantitative trait loci for growth traits on bovine chromosomes 2, 6, 14, 19, 21, 23 within one commercial line of Bos taurus.

    J Anim Sci 2004, 82:3405-3414. PubMed Abstract | Publisher Full Text OpenURL

  5. Maltecca C, Weigel KA, Khatib H, Cowan M, Bagnato A: Whole-genome scan for quantitative trait loci associated with birth weight, gestation length and passive immune transfer in a Holstein × Jersey crossbred population.

    Anim Genet 2008, 40:27-34. PubMed Abstract | Publisher Full Text OpenURL

  6. McClure MC, Morsci NS, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, Taylor JF: A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle.

    Anim Genet 2010, 41:597-607. PubMed Abstract | Publisher Full Text OpenURL

  7. Takasuga A, Watanabe T, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujikawa A, Chiba K, Kobayashi N, Tatsuda K, Oe T, Furukawa-Kuroiwa M, Nishimura-Abe A, Fujita T, Inoue K, Mizoshita K, Ogino A, Sugimoto Y: Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping.

    Mamm Genome 2007, 18:125-136. PubMed Abstract | Publisher Full Text OpenURL

  8. Setoguchi K, Furuta M, Hirano T, Nagao T, Watanabe T, Sugimoto Y, Takasuga A: Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate.

    BMC Genet 2009, 10:43. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Snelling WM, Allan MF, Keele JW, Kuehn LA, Thallman RM, Bennett GL, Ferrell CL, Jenkins TG, Freetly HC, Nielsen MK, Rolfe KM: Partial-genome evaluation of postweaning feed intake and efficiency of crossbred beef cattle.

    J Anim Sci 2011, 89:1731-1741. PubMed Abstract | Publisher Full Text OpenURL

  10. Malik S, Roeder RG: The metazoan Mediator co-activator complex as an integrative hub for transcriptional regulation.

    Nat Rev Genet 2010, 11:761-772. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Eberlein A, Takasuga A, Setoguchi K, Pfuhl R, Flisikowski K, Fries R, Klopp N, Furbass R, Weikard R, Kuhn C: Genetic factors modulating fetal growth in cattle indicates a potential role of the non-SMC condensin I complex subunit G (NCAPG) gene.

    Genetics 2009, 183:951-964. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Seipold S, Priller FC, Goldsmith P, Harris WA, Baier H, Abdelilah-Seyfried S: Non-SMC condensing I complex proteins control chromosomal segregation and survival of proliferating cells in the zebrafish neural retina.

    BMC Dev Biol 2009, 9:40. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Soranzo N, Rivadeneira F, Chinappen-Horsley U, Malkina I, Richards JB, Hammond N, Stolk L, Nica A, Inouye M, Hofman A, Stephens J, Wheeler E, Arp P, Gwilliam R, Jhamai PM, Potter S, Chaney A, Ghori MJ, Ravindrarajah R, Ermakov S, Estrada K, Pols HA, Williams FM, McArdle WL, van Meurs JB, Loos RJ, Dermitzakis ET, Ahmadi KR, Hart DJ, Ouwehand WH, Wareham NJ, Barroso I, Sandhu MS, Strachan DP, Livshits G, Spector TD, Uitterlinden AG, Deloukas P: Meta-analysis of genome-wide scans for human adult stature identifies novel Loci and associations with measures of skeletal frame size.

    PLoS Genet 2009, 5:e1000445. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Sovio U, Bennett AJ, Millwood IY, Molitor J, O'Reilly PF, Timpson NJ, Kaakinen M, Laitinen J, Haukka J, Pillas D, Tzoulaki I, Molitor J, Hoggart C, Coin LJ, Whittaker J, Pouta A, Hartikainen AL, Freimer NB, Widen E, Peltonen L, Elliott P, McCarthy MI, Jarvelin MR: Genetic determinants of height growth assessed longitudinally from infancy to adulthood in the northern Finland birth cohort 1966.

    PLoS Genet 2009, 5:e1000409. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. FASS: Guide for Care and Use of Agricultural Animals in Agricultural Research and Teaching. In Fed Animal Science Society. 1st edition. Savoy, IL; 1999. OpenURL

  16. Kuehn LA, Rohrer GA, Nonneman DJ, Thallman RM, Leymaster KA: Detection of single nucleotide polymorphisms associated with ultrasonic backfat depth in a segregating Meishan × White Composite population.

    J Anim Sci 2007, 85:1111-1119. PubMed Abstract | Publisher Full Text OpenURL

  17. Rolfe KM, Snelling WM, Nielsen MK, Freetly HC, Ferrell CL, Jenkins TG: Genetic and phenotypic parameter estimates for feed intake and other traits in growing beef cattle and opportunities for selection.

    J Anim Sci 2011, 89:3452-3459. PubMed Abstract | Publisher Full Text OpenURL

  18. King DA, Shackelford SD, Kuehn LA, Kemp CM, Rodriguez AB, Thallman RM, Wheeler TL: Contribution of genetic influences to animal-to-animal variation in myoglobin content and beef lean color stability.

    J Anim Sci 2010, 88:1160-1167. PubMed Abstract | Publisher Full Text OpenURL

  19. Wheeler TL, Cundiff LV, Shackelford SD, Koohmaraie M: Characterization of biological types of cattle (Cycle VIII): Carcass, yield, and longissimus palatability traits.

    J Anim Sci 2010, 88:3070-3083. PubMed Abstract | Publisher Full Text OpenURL

  20. Shackelford SD, Wheeler TL, Koohmaraie M: Evaluation of slice shear force as an objective method of assessing beef longissimus tenderness.

    J Anim Sci 1999, 77:2693-2699. PubMed Abstract | Publisher Full Text OpenURL

  21. Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by Krawetz S, Misener S. Totowa. New Jersey: Human Press; 2000:365-386. OpenURL

  22. Boldman KG, Kriese LA, Van Vleck LD, Van Tassel CP, Kachman SD: A manual for the use of MTDFREML. A set of programs to obtain estimates of variance and covariances (Draft). USDA, ARS, Lincoln, NE; 1995.

  23. Kuehn LA, Keele JW, Thallman RM: Estimation of current breed differences in multibreed genetic evaluations using quantitative and molecular approaches. Beef Improvement Federation: 9th Genetic Prediction Workshop.

    Prediction of Genetic Merit of Animals for Selection 2008, 49-60. OpenURL

  24. Goddard ME, Hayes BJ: Mapping genes for complex traits in domestic animals and their use in breeding programs.

    Nature Rev Genet 2009, 92:381-391. OpenURL

  25. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps.

    Bioinformatics 2005, 21:263-265. PubMed Abstract | Publisher Full Text OpenURL

  26. Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, Lawley CT, Haudenschild CD, Moore SS, Warren WC, Sonstegard TS: SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries.

    Nat Methods 2008, 5:247-252. PubMed Abstract | Publisher Full Text OpenURL

  27. Bovine HapMap Consortium: Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds.

    Science 2009, 324:528-532. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Bovine Genome Sequencing and Analysis Consortium: The genome sequence of taurine cattle: a window to ruminant biology and evolution.

    Science 2009, 324:522-528. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Weikard R, Altmaier E, Suhre K, Weinberger KM, Hammon HM, Albrecht E, Setoguchi K, Takasuga A, Kühn C: Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on Bos taurus growth and lipid deposition.

    Physiol Genomics 2010, 42A:79-88. PubMed Abstract | Publisher Full Text OpenURL

  30. Yao K, Yin YL, Chu WY, Li ZQ, Deng D, Li TJ, Huan RL, Zhang JS, Tan B, Wang W, Wu G: Dietary arginine supplementation increases mTOR signaling activity in skeletal muscle of neonatal pigs.

    J Nutr 2008, 138:867-872. PubMed Abstract | Publisher Full Text OpenURL

  31. Lira VA, Soltow QA, Long JHD, Betters JL, Sellman JE, Criswell DS: Nitric oxide increases GLUT4 expression and regulated AMPK signaling in skeletal muscle.

    Am J Physiol Endocrinol Metab 2007, 293:E1062-E10628. PubMed Abstract | Publisher Full Text OpenURL

  32. Roy D, Perreault M, Marette A: Insulin stimulation of glucose uptake in skeletal muscles and adipose tissues in vivo is NO dependent.

    Am J Physiol Endocrinol Metab 1998, 274:E692-E699. OpenURL