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

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][4][5][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.

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 F 1 , born in 1999, 2000, and 2001. Female F 1 and the 2001-born F 1 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 F 1 2 . The F 1 2 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 F 1 2 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) F 1 2 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 11 th 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). 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 MgCl 2 ; 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.

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 (r 2 ) was defined for the 47 SNP on chromosome 6 using Haploview 4.0 software [25] http://www.broad.mit.edu/mpg/haploview/index.php. 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.

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).

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][27][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).
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).

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 r 2 = 0.89 than the LD between these two markers in the SNP discovery population (r 2 = 0.96) (Figure 2). 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  3 Position based on the Btau 4.0 Bovine Genome Assembly. 4 Number of genotypic calls used for statistical analysis. 5 ADFI and ADG effects are displayed as kg/d. Effects are presented for the minor allele. 6 Correction for multiple testing performed by multiplying the P-values by the number of markers tested (n = 47). NS = not significant, P-value > 0.05. 7 Markers evaluated in Setoguchi et al. [8] and Eberlein et al. [11]. 8 Marker evaluated in Charolais and German Holstein population in Eberlein et al. [11]. 9 Markers in bold font, Hapmap28589-BTC-041045 and 81405_282, were significantly associated with ADG in the validation population of animals (P = 0.04 and 0.03, respectively).  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.  [27], and Bovine Genome Sequencing and Analysis Consortium et al. [28]. 3 Effects of phenotypes are displayed in the following units: AFT = cm, HCW = kg, REA = cm 2 , Marbling = MSU, and SSF = kg. The effects presented correspond to the minor allele. 4 Markers evaluated in Setoguchi et al. [8] and Eberlein [11]. 5 Marker evaluated in Charolais and German Holstein population in Eberlein et al. [11].
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,  [27], and Bovine Genome Sequencing and Analysis Consortium et al. [28]. 2 Markers with check marks are from the BovineSNP50 Illumina BeadChip. These markers were significantly associated with ADFI and/or ADG in the discovery population of crossbred steers [9]. 3 Position based on the Btau 4.0 Bovine Genome Assembly. 4 NP corresponding to the effects presented. 5 ADFI and ADG effects are displayed as kg/d. Effects are shown that correspond to the minor allele. 6 Markers in bold font, Hapmap28589-BTC-041045 and 81405_282, were significantly associated with ADG in the validation population of animals (P = 0.04 and 0.03, respectively).
ADG and the carcass/meat quality traits HCW, AFT and REA. All of these SNP are in LD (r 2 ≤ 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 r 2 = 0.96 in our SNP discovery population and r 2 = 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 12 th 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 Figure 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 r 2 is shown between markers with shading. The darker shading represents higher LD, black diamonds are r 2 = 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]. 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. Figure 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 r 2 is shown between markers by shading. The darker shading represents higher LD, black diamonds are r 2 = 1. Markers in bold signify those that are located in a haplotype block.