Genome-wide association analysis of feed intake and residual feed intake in Nellore cattle
BMC Genetics volume 15, Article number: 21 (2014)
Feed intake plays an important economic role in beef cattle, and is related with feed efficiency, weight gain and carcass traits. However, the phenotypes collected for dry matter intake and feed efficiency are scarce when compared with other measures such as weight gain and carcass traits. The use of genomic information can improve the power of inference of studies on these measures, identifying genomic regions that affect these phenotypes. This work performed the genome-wide association study (GWAS) for dry matter intake (DMI) and residual feed intake (RFI) of 720 Nellore cattle (Bos taurus indicus).
In general, no genomic region extremely associated with both phenotypic traits was observed, as expected for the variables that have their regulation controlled by many genes. Three SNPs surpassed the threshold for the Bonferroni multiple test for DMI and two SNPs for RFI. These markers are located on chromosomes 4, 8, 14 and 21 in regions near genes regulating appetite and ion transport and close to important QTL as previously reported to RFI and DMI, thus corroborating the literature that points these two processes as important in the physiological regulation of intake and feed efficiency.
This study showed the first GWAS of DMI to identify genomic regions associated with feed intake and efficiency in Nellore cattle. Some genes and QTLs previously described for DMI and RFI, in other subspecies (Bos taurus taurus), that influences these phenotypes are confirmed in this study.
Feed intake plays an important economic role in cattle growth and may represent the greatest costs in beef cattle, both in beef cattle finishing systems and calves production . Feed intake is evaluated by dry matter intake (DMI) in cattle, which is associated with weight gain, carcass traits and feed efficiency . It is, therefore, a relevant variable for the entire meat production system and, possibly, a trait to be included in genetic breeding programs.
Although genetic breeding foci primarily on growth and reproductive traits, other variables must be taken into account because increased weight gain rates may lead to a concomitant increase in the adult size generating higher maintenance costs of animals . Thus, the objectives of genetic breeding must be adequately delineated to attain balance between weight gain rates and other essential characteristics, such as reproduction , meat quality and feeding efficiency .
In recent decades, many efforts have been made to better balance the relationship of weight gain with feed intake in beef cattle; however, some traits are less effective to minimize the negative response of increased adult animal size, such as gross feeding efficiency [3, 5]. In this respect, the residual feed intake (RFI) was proposed in the 1960s  and has gained more notability for being independent from growth and body size, and designed towards intake reduction [3, 5–8]. This independence is attributed to the fact that RFI is calculated as the difference between observed and estimated intake by a regression equation of DMI over the average daily gain (ADG) and the mid body weight (MBW) .
Both DMI as RFI in cattle have enough variability and heritability to respond to genetic selection [2, 7, 8]; however, there is still no consensus on how these traits should be considered in the indices of selection. In this context, two major limitations comprise the difficulties inherent to obtain phenotypes and how to use this information in the breeding process . In Nellore cattle (Bos taurus indicus), these difficulties are even greater, given the few animals with phenotype available for DMI and RFI, and the scarce knowledge about genetic parameters of these traits in Nellore. However, the phenotypic variability of RFI in Nellore show a standard deviation ranging between 0.31 and 0.69 kg/DM per day [10, 11], is similar to that observed in Bos taurus taurus (often referred as taurine) animals and crossbreed [2, 3, 6, 7].
The use of genomic information can be a strategy to improve the selection of phenotypes such as RFI and DMI, if the marker effects are estimated accurately. Genome-wide association study (GWAS) allowed to identify subsets of markers that explain an important portion of the variation of these characteristics [9, 12, 13]. The use of the information obtained from these markers along the chromosomes (BTA) can improve the accuracy of young animals candidates for genetic selection, and thus improve the genetic gain by reducing the generation interval.
Several studies have reported the viability of using the information from single nucleotide polymorphism (SNP) to identify regions of the genome that affect phenotypes of interest, aiming at improving breeding schemes for weight gain, reproduction and carcass traits in beef cattle [14–16]. Additionally, studies on molecular markers in cattle were enhanced with the recent release of the reference bovine genome  and with the improvement of beadchip technologies that perform fast and automated analyses of hundreds of thousands of SNPs and with the decreasing cost per SNP analyzed. The development of high-density commercial panels of SNPs opened a range of opportunities for GWAS . Furthermore, the imputation of genotypes has proven to be an effective tool in enhancing the power of GWAS by increasing the number of genotyped animals and can be a valuable strategy for reducing even more the genotyping cost .
However, the vast majority of GWAS has been performed in animals of the taurine subspecies. Also, the first beadchip of thousands of SNPs were developed based on this subspecies, which causes several SNPs, described as being polymorphic in taurines, to be non-informative in zebu cattle (Bos taurus indicus), especially Nellore [19, 20]. Only in recent years, GWAS was carried out more often in zebu from the development of optimized beadchips also taking into account this subspecies (e.g. ). The objectives of this study were: 1) to identify SNPs associated with DMI and RFI in Nellore cattle, using medium (Illumina® BovineSNP50 v2 BeadChip), high density (Illumina® BovineHD BeadChip) and a combined of medium to high-density panels by imputation; and 2) to explore the regions surrounding the identified markers seeking possible genes with known function near these SNPs.
Phenotype, quality assurance and imputation
For the three datasets of the two phenotypic variables, no evidences were found to deviate from normality and homoscedasticity of model residuals in the Shapiro-Wilk and Breusch-Pagan tests (P > 0.05), respectively. The mean, additive genetic variances and residual variances were 8.76 ± 1.96, 0.29 and 0.42, respectively to DMI and 0.00 ± 0.89, 0.20 and 0.33, respectively to RFI. However, four samples were considered outliers in the DMI distribution and removed from the HD dataset and HDimp dataset, whereas two samples were considered outliers for RFI. Regarding the criteria for exclusion of samples and SNPs, the results of the quality control of 50 k and HD are shown in Table 1.
After the quality control criteria, the final datasets were: 672 samples in 50 k with 28,231 SNPs and 365 samples in HD with 436,588 SNPs.
The number of SNPs imputed with over 95% accuracy was 290,620 and the number of remaining samples was 672 for the dataset HDimp. The median imputation accuracy was 97.2% with an average of 94.6%. Only those genotypes imputed with over 95% accuracy were used in the GWAS.
Genome-wide association study
The deflation/inflation factor (λ) calculated for all association analyses was lower than 1.1 (Additional files 1, 2, 3, 4, 5 and 6), which was considered acceptable and can be used on genomic control correction approach (GC). The P-values of SNPs along the chromosomes are shown in the form of Manhattan Plots for DMI and RFI in Figures 1 and 2, respectively, with the threshold represented as the Bonferroni significance line.
The SNPs significant in the Bonferroni adjustment in each DMI and RFI analysis are described in Table 2.
The exploration of the region around the five SNPs that were associated with DMI and RFI are shown in Table 3, which shows all the genes with known functions located around (100 kb) these SNPs. Several other genes are at shorter distances; however, these distances were not included because their function is still unknown according to Ensembl genes 72 UMD v 3.1 . Additionally, the distance of the QTLs mapped for both phenotypic variables closer to these SNPs is shown in the same table.
The thresholds for the SNP exclusion are still questionable. The accuracy load of each criterion depends heavily on the dataset and information that it is expected from each type of GWAS. The criterion that eliminates most SNPs and, possibly, the most controversial is MAF. In general, there is a tendency to use 5% for GWAS [9, 14, 16, 22–24]; however, in this study, we adopted 2% of MAF, even with approximately 25% of SNPs being eliminated in this criterion. Other studies also adopted the same level [20, 25, 26].
The GWAS across the 29 autosomal chromosomes showed no genomic region associated with DMI and RFI of feedlot Nellore young bulls and steers, as expected for the variables that have their regulation controlled by many genes. This fact was reported in other studies that also assessed DMI and RFI in different cattle subspecies [9, 14] and RFI in pigs .
The benefits of using imputation were observed in this study because, in the dataset HDimp, we found regions significantly associated with phenotypes, which was not noted in the dataset HD even though both had a very close Bonferroni threshold. Although the number of SNPs of the dataset of imputed animals is much lower than that of the dataset HD (290,620 vs 436,588), this tool enabled a much more robust association study due to the considerable increase in the number of samples (672 vs 365).
The two variables showed high genetic and phenotypic correlation among each other [2, 8, 28], but not always the same regions showed strong association with these characteristics. This can be partly explained by the difference between the physiological mechanisms that regulate RFI are not exactly the same that regulate DMI. On the other hand, regions with important effect on the two traits suggest the existence of pleiotropic effects on these variables [29, 30]. However, some regions are well evidenced in both analysis of the same trait and, in some cases, we can observe genomic regions that relate to both, such as in BTA4, BTA8 and BTA14.
Three SNPs surpassed the threshold for the Bonferroni multiple test for DMI and two SNPs for RFI. Several markers have been associated with these two variables in the literature [9, 12, 14, 22–24, 27, 29, 30]; however, the methodologies used for this purpose are diverse and populations assessed are extremely distinct, which may imply that associations made in a particular breed may not be applied in others . These SNPs can explain part of the phenotypic variance, insomuch that few markers explain more than 30% of the variation in RFI [9, 27, 30]. However, this calculation takes into account allele frequencies, the allele substitution effect and phenotypic variance of the trait. This prediction can be overrated depending on these factors, mainly when it assumes independence between the markers considered in this calculation. The allele substitution effect of the SNPs varied between the panels, and this effect in DMI was higher for markers in the 50 k panel, for RFI in the HDimp panel.
Regarding the location of SNPs related to DMI, the SNP rs109784719 (BTA14) is at 27.4 kb of the single gene (STMN2); however, it is in the region of QTL #4365 described previously for DMI. Other studies found SNP associated with DMI in the BTA14 in beef cattle [9, 24] and in chromosome 14 in pigs , both in a region surrounding the genes such as PLAG1, RDHE2 and CHCHD7 that notably influence the stature of various species [20, 24, 33]. The SNP rs29024524, in the BTA8, is in gene CCDC171 and surrounding other genes; moreover, it lies next to the QTL #4425, which is a genomic region that seems to affect RFI and DMI.
Other two QTLs (#4353 and #5274) involved with RFI surround this SNP and SNP rs41660853, associated with RFI. This SNP linked to RFI is located near gene CLCN3, related to the ionic transport processes already reported as influential in RFI , corroborating a previous study that described the metabolic differences of RFI . Three other association studies indicated the importance of this region in RFI and DMI [14, 23, 29].
The ionic transport system can account for more than 10% of all ruminant energy expenditure ; therefore, for animals with lower energy expenditure to maintain this system, they can, at the same time, direct the energy consumed to other processes, which can influence feeding efficiency. The SNP rs135777172 also lies in a region of QTL (#4462) already described for RFI at final part of BTA21, where there is an abundance of genes encoding various types of proteins .
Finally, SNP rs134003539 is in a locus described as important, for both DMI and RFI [12, 14, 23, 29, 36, 37], at 4 Mb of the QTL (#10584). Other two genes widely studied in both phenotypic variables are also located in this locus. The neuropeptide Y (NPY, Gene ID 504216) positioned at 1.4 Mb of distance and the leptin (LEP, Gene ID 280836) at 19 Mb from the SNP and inserted in this QTL. The neuropeptides and hormones that control appetite, energy expenditure and metabolism of fat and glucose have a relationship to each other, which seems to be mediated by leptin . These compounds can increase appetite such as the neuropeptide Y, ghrelin, AGRP, MCH, orexines and noradrenalines, or they can even have a reducing effect on appetite such as leptin, POMC, CART, CRH, α-MSH and serotonin . In addition, the rs134003539 is in gene ZNF804B that is a form of protein “zinc finger”, characterized by coordination and stabilization of one or more zinc ions in several ionic exchange processes .
In general, the GWAS presents itself as an interesting tool to identify genomic regions that can influence these phenotypes. Another potential advantage is the creation of less dense panels designed for a breed or subspecies of interest containing SNPs directed to a certain characteristic. This could reduce costs of a possible genetic selection with the use of genomic data, improving the accuracy of estimates of genetic value in animals.
The identification of regions associated with DMI and RFI may elucidate loci that influence these variables and highlight possible important physiological mechanisms. In this study, some promising regions were identified, with important effect on DMI and RFI. However, the strict significance level adopted and the number of genotyped animals may have contributed to reducing the power of the present study, given that important regions may not have been identified, which suggests the need for further studies aimed at overcoming such restrictions.
The current study showed the results of GWAS in high and medium-density panels to identify genomic regions associated with feed intake and efficiency. It is believed to be the first study of this kind conducted in Nellore cattle to these traits. The results suggest that RFI and DMI are influenced by loci previously described and these results obtained in zebu cattle are supported by the literature.
The study was conducted from a compilation of 11 experiments of feed efficiency and intake conducted in Brazil from 2007 to 2011. One experiment was performed in the South , eight in the Southeast [10, 19] and two in the Midwest region in Brazil. In these experiments, we obtained phenotypes of 720 young bulls and steers at 550 ± 115 days old, with body weight of 380 ± 51 kg. The sample size of each trial ranged from 50 to 120 animals. Three different facilities were used to measure the phenotypes. We used automated systems of collective stalls (Calan gates and GrowSafe™) and an individual pens system. The experiment lasted at least 70 days, and the dietary intake was measured daily with subsequent adjustment for dry matter content.
During the experiments, the MBW of each animal was obtained by periodic weighing without fasting and ADG was calculated as the slope of the linear regression of weights compared to the testing days. To calculate the RFI, we considered the regression equation residue of DMI on ADG and MBW ; however, because of the different experiments, we formed contemporary groups (CG), which means that animals of the same experiment were considered as a CG. Within three of the CG, there were young bulls and castrated steers; therefore, the sexual condition (SC) was also included as a covariate in the statistical model. The RFI was calculated in the PROC REG procedure from the statistical package SAS 9.3 and the general equation was:
RFI and DMI were tested for data normality (Shapiro-Wilk) and homoscedasticity of model residues (Breusch-Pagan). We also performed the control of data outliers (DMI, ADG, MBW and RFI), in which the records outside the mean range ± 3 standard deviations, classified as possible measurement error, were excluded from the analysis.
DNA extraction, sample assurance and SNP quality control
Blood was collected from all animals evaluated by puncturing the jugular. We used tubes containing K2 EDTA as anticoagulant. The samples were stored at 4°C for late DNA extraction by NaCl precipitation . After extraction, the quality of the DNA samples was assessed by determining the ratio A260/280 in biophotometer. The samples were only accepted when values remained between 1.8 and 2.0 concomitantly, the DNA was, then, quantified and samples were diluted to a minimum concentration of 50 ng/μL and maximum of 150 ng/μl for subsequent genotyping.
Genotyping was performed in two types of DNA beadchip: the Illumina® BovineHD BeadChip (777,962 SNPs) , and the Illumina® BovineSNP50 v2 BeadChip (54,609 SNPs) , both in the standard test Infinium Assay II for the HiScanSQ® platform (Illumina Inc., San Diego, CA). Genotype calls were determined using the validated standard cluster file provided by the manufacturer, with GenCall Scores greater than 0.70. In total, 720 animals were genotyped, including 336 young bulls and steers in BovineSNP50 (50 k dataset) and 384 young bulls in BovineHD (HD dataset). As most of the BovineSNP50 markers are embedded in the BovineHD panel, the HD samples were also included in the 50 k dataset in order to increase sample size. We assessed the occurrence of duplicate samples by calculating the proportion of alleles identical by state (IBS) between all possible pairs of individuals. For IBS calculation, genotypes were considered for 5,000 and 20,000 markers taken randomly for BovineSNP50 and BovineHD, respectively, and any pair of samples with IBS ≥ 95% were excluded.
For GWAS, SNPs were subjected to a quality control in which only autosomal SNPs with known genomic coordinate were considered. Samples with Call rate (IDCR) lower than 90% were removed from the study. Markers were removed if they presented minor allele frequency (MAF) ≤ 0.02, Call rate (SNPCR) ≤ 0.95 and P-value for Fisher’s exact test for Hardy-Weinberg Equilibrium (HWE) ≤ 1 x 10-5. This quality control was performed in R v2.15.2 software using scripts developed for this purpose  and the GenABEL v1.7-6 package .
In order to increase sample size, an imputation analysis was performed to combine the available information of animals genotyped with BovineSNP50 and BovineHD. The HD panel information was used to verify the imputation efficiency from the 50 k Panel to the HD. The quality control of SNPs was performed again, which excluded SNPs located in non-autosomal chromosomes, SNPCR ≤ 0.97, HWE ≤ 1 × 10-6, IDCR ≤ 0.90.
Afterwards, 290 animals randomly sampled (genotyped in high-density panel) were considered as population reference for the imputation analysis, while the remaining animals of the panel were included in a validation set. Except for the markers in common between the two technologies, the genotyped animals in the validation BovineHD had their genotypes masked and imputed, which allowed to simulate the imputation of 50 k to HD. These analyses were performed using the FImpute 2.2 software , using HD genotypes in combination with pedigree information. The imputation efficiency was assessed using the proportion of genotypes imputed correctly. Subsequently, imputation was performed similarly to all animals in the 50 k, in which all genotyped animals in HD that passed the quality control criteria (362) were considered as reference population. The final dataset (HDimp dataset) was composed of 672 animals and 290,620 SNPs that were allocated with accuracy greater than 95%.
The association analysis was based on a variance-components method, namely Grammar-Gamma . This method is a computationally efficient unbiased approximation of the gold standard likelihood ratio test (LRT) [45, 46], which corrects the association analysis for confounding due to genetic substructure and relatedness.
The first step in the association analysis was to use the variance-covariance matrix as genomic relationship matrix, to correct for relatedness and substructure and Grammar-Gamma factor calculation. To ensure the reliability of the estimates, we calculated the inflation/deflation factor (λ) for the correction of the GC .
The second step consisted of associating the phenotype with the genotypes without the variance-covariance matrix and the estimates of SNPs effects adjusted by the Grammar-Gamma factor . The model used was the polygenic  and the association test varied depending on the phenotype used. The polygenic model for RFI included age and SC as covariate and, for DMI was included SC and MBW as covariates. The results were presented as Manhattan Plots in which the -log10(P- values) were plotted, corrected for λ and the Bonferroni correction (α = 0.05/number of SNPs) was considered as significance threshold. All the procedures described in this item were carried out in R v2.15.2, using the GenABEL v1.7-6 package .
Region surrounding significant SNPs
The SNPs that surpassed the threshold of the Bonferroni adjustment were described and their allelic substitution effects were reported. The proportion of the phenotypic variance explained by the SNPs was also estimated as:
Where p and q are the allele frequencies of the analyzed SNP, β2 is the square of the allele substitution effect and S2 is the total variance of phenotype.
Additionally, they were investigated for their genomic location (genes surrounding and possible QTLs already mapped for DMI and RFI). The exploration of the region searched for genes with known functions located at the maximum 100 kb distant from the SNPs in Ensembl genes 72 using the assembly UMD v3.1 . The search for QTLs was examined in cattle QTLdb database .
Anderson RV, Rasby RJ, Klopfenstein TJ, Clark RT: An evaluation of production and economic efficiency of two beef systems from calving to slaughter. J Anim Sci. 2005, 83: 694-704.
Nkrumah JD, Basarab JA, Wang Z, Li C, Price MA, Okine EK, Crews DH, Moore SS: Genetic and phenotypic relationships of feed intake and measures of efficiency with growth and carcass merit of in beef cattle. J Anim Sci. 2007, 85: 2711-2720. 10.2527/jas.2006-767.
Archer JA, Richardson EC, Herd RM, Arthur PF: Potential for selection to improve efficiency of feed use in beef cattle: a review. Aust J Agric Res. 1999, 50: 147-161. 10.1071/A98075.
Phocas F, Bloch C, Chapelle P, Bécherel F, Renand G, Ménissier F: Developing a breeding objective for a French purebred beef cattle selection programme. Livest Prod Sci. 1998, 57: 49-65. 10.1016/S0301-6226(98)00157-2.
Kennedy BW, van der Werf JHJ, Meuwissen THE: Genetic and statistical properties of residual feed intake. J Anim Sci. 1993, 71: 3239-3250.
Koch RM, Swiger LA, Chambers D, Gregory KE: Efficiency of feed use in beef cattle. J Anim Sci. 1963, 22: 486-494.
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: 111-119. 10.1016/S0301-6226(99)00122-0.
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: 131-139. 10.1016/S0301-6226(00)00243-8.
Rolf MM, Taylor JF, Schnabel RD, McKay SD, McClure MC, Northcult SL, Kerley MS, Weaber RL: Genome-wide association analysis for feed efficiency in Angus cattle. Anim Genet. 2011, 43: 367-374.
Gomes RC, Sainz RD, Silva SL, César MC, Bonin MN, Leme PR: Feedlot performance, feed efficiency reranking, carcass traits, body composition, energy requirements, meat quality and calpain system activity in Nellore steers with low and high residual feed intake. Livest Sci. 2012, 150: 265-273. 10.1016/j.livsci.2012.09.012.
Sobrinho TL, Branco RH, Bonilha SFM, Castilhos AM, Figueiredo LA, Razook AG, Mercadante MEZ: Residual feed intake and relationships with performance of Nellore cattle selected for post weaning weight. Rev Bras Zoot. 2011, 40: 929-937. 10.1590/S1516-35982011000400030.
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: 1893-905. 10.1534/genetics.107.072637.
Moore SS, Mujibi FD, Sherman EL: Molecular basis for residual feed intake in beef cattle. J Anim Sci. 2009, 87 (E Suppl): 41-47.
Bolormaa S, Hayes BJ, Savin K, Hawken R, Barendse W, Arthur PF, Johnston DJ, Herd RM, Goddard ME: Genomewide association studies for feedlot and growth traits in cattle. J Anim Sci. 2011, 89: 1684-1697. 10.2527/jas.2010-3079.
Cole JB, Van Raden PM, O’Connell JR, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Wiggans GR: Distribution and location of genetic effects for dairy traits. J Dairy Sci. 2009, 92: 2931-2946. 10.3168/jds.2008-1762.
Snelling WM, Allan MF, Keele JW, Kuehn LA, McDaneld T, Smith TPL, Sonstegard TS, Thallman RM, Bennett GL: Genome-wide association study of growth in crossbred beef cattle. J Anim Sci. 2010, 88: 837-848. 10.2527/jas.2009-2257.
Bovine Genome Sequencing and Analysis Consortium: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324: 522-528.
VanRaden PM, O’Connell JR, Wiggans GR, Weigel KA: Genomic evaluations with many more genotypes. Genet Sel Evol. 2011, 43: 10-10.1186/1297-9686-43-10.
Gomes RC, Silva SL, Carvalho ME, Rezende FM, Pinto LFB, Santana MHA, Stella TR, Meirelles FV, Rossi Junior P, Leme PR, Ferraz JBS: Protein synthesis and degradation gene SNPs related to feed intake, feed efficiency, growth and ultrasound carcass traits in Nellore cattle. Gen Mol Res. 2013, 12: 2923-2936. 10.4238/2013.August.12.8.
Utsunomiya YT, Carmo AS, Carvalheiro R, Neves HHR, Matos MC, Zavarez LB, O’Brien AMP, Sölkner J, McEwan JC, Cole JB, Van Tassel CP, Schenkel FS, Silva MVGB, Porto Neto LR, Sonstegard TS, Garcia JF: Genome-wide association study for birth weight in Nellore cattle points to previously described orthologous genes affecting human and bovine height. BMC Genetics. 2013, 14: 52-64.
Kinsella RJ, Kähäri A, Haider S, Zamora J, Proctor G, Spudich G, Almeida-King J, Staines D, Derwent P, Kerhornou A, Kersey P, Flicek P: Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database (Oxford). 2011, 2011: Bar030-
Snelling WM, Allan MF, Keele JW, Keuhn 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. 10.2527/jas.2010-3526.
Mujibi FDN, 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: 3362-3371. 10.2527/jas.2010-3362.
Pryce JE, Arias J, Bowman PJ, Davis SR, Macdonald KA, Waghorn GC, Wales WJ, Williams YJ, Spelman RJ, Hayes BJ: Accuracy of genomic predictions of residual feed 14 intake and 250 day bodyweight in 15 growing heifers using 625,000 SNP markers. J Dairy Sci. 2012, 95: 2108-2119. 10.3168/jds.2011-4628.
Wang WYS, Barratt BJ, Clayton DG, Todd JA: Genome wide association studies: theoretical and practical concerns. Nat Rev Genet. 2005, 6: 109-118. 10.1038/nrg1522.
Andersson CA, Petterson FM, Clarke GM, Cardon LR, Morris AP, Zondervan KT: Data quality control in genetic case–control association studies. Nat Protoc. 2010, 5 (9): 1564-1573. 10.1038/nprot.2010.116.
Onteru SK, Gorbach DM, Young JM, Garrick DJ, Dekkers JC, Rothschild MF: Whole genome association studies of residual feed intake and related traits in the pig. PloS One. 2013, 8 (6): e61756-10.1371/journal.pone.0061756.
Santana MHA, Rossi P, Almeida R, Cucco DC: Feed efficiency and its correlations with carcass traits measured by ultrasound in Nellore bulls. Livest Sci. 2012, 145: 252-257. 10.1016/j.livsci.2012.02.012.
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: 16-22. 10.2527/jas.2008-1759.
Sherman EL, Nkrumah JD, Murdoch BM, Moore SS: Identification of polymorphisms influencing feed intake and efficiency in beef cattle. Anim Genet. 2008, 39: 225-231. 10.1111/j.1365-2052.2008.01704.x.
Hayes BJ, Bowman PJ, Chamberlain AC, Verbyla K, Goddard ME: Accuracy of genomic breeding values in multi-breed dairy cattle populations. Genet Sel Evol. 2009, 41: 51-10.1186/1297-9686-41-51.
Do DN, Strathe AB, Ostersen T, Jensen J, Mark T, Kadarmideen HN: Genome-wide association study reveals genetic architecture of eating behavior in pigs and its implications for humans obesity by comparative mapping. PLoS ONE. 2013, 8 (8): e71509-10.1371/journal.pone.0071509.
Fortes MR, Kemper K, Sasazaki S, Reverter A, Pryce JE, Barendse W, Bunch R, McCulloch R, Harrison B, Bolormaa S, Zhang YD, Hawken RJ, Goddard ME, Lehnert SA: Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle. Anim Genet. 2013, 44 (6): 636-647. 10.1111/age.12075.
Richardson EC, Herd RM, Archer JA, Arthur PF: Metabolic differences in Angus steers divergently selected for residual feed intake. Aust J Exp Agric. 2004, 44: 441-452. 10.1071/EA02219.
McBride BW, Kelly JM: Energy cost of absorption and metabolism in the ruminant gastrointestinal tract and liver: a review. J Anim Sci. 1990, 68: 2997-3006.
Banos G, Woolliams JA, Woodward BW, Forbes AB, Coffey MP: Impact of single nucleotide polymorphisms in leptin, leptin receptor, growth hormone receptor, and diacylglycerol acyltransferase (DGAT1) gene loci on milk production, feed, and body energy traits of UK dairy cows. J Dairy Sci. 2008, 91: 3190-3200. 10.3168/jds.2007-0930.
Nkrumah JD, Li C, Basarab JB, Guercio S, Meng Y, Murdoch B, Hansen C, Moore SS: Association of a single nucleotide polymorphism in the bovine leptin gene with feed intake, feed efficiency, growth, feeding behaviour, carcass quality and body composition. Can J Anim Sci. 2004, 84: 211-219. 10.4141/A03-033.
Campfield LA, Smith FJ, Burn P: The OB protein (leptin) pathway - a link between adipose tissue mass and central neural networks. Horm Metab Res. 1996, 28: 619-32. 10.1055/s-2007-979867.
Klug A, Rhodes D: 'Zinc fingers’: a novel protein motif for nucleic acid recognition. Trends Biochem Sci. 1987, 12: 464-467.
Olerup O, Zetterquist H: HLA-DR typing by PCR amplification with sequence-specific primers (PCRSSP) in 2 hours: an alternative to serological DR typing in clinical practice including donor-recipient matching in cadaveric transplantation. Tis Antig. 1992, 39: 225-235. 10.1111/j.1399-0039.1992.tb01940.x.
Matukumalli LK, Schroeder S, DeNise SK, Sonstegard TS, Lawley CT, Georges M, Coppieters W, Gietzen K, Medrano JF, Rincon G, Lince D, Eggen A, Glaser L, Cam G, Van Tassel CP: Analyzing LD blocks and CNV segments in cattle: novel genomic features identified using the BovineHD BeadChip. 2011, San Diego, CA: Illumina Inc, Pub. No. 370–002
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connel J, Moore SS, Smith TPL, Sonstegard TS, Van Tassel CP: Development and characterization of a high density SNP genotyping assay for cattle. PLoS ONE. 2009, 4: e5350-10.1371/journal.pone.0005350.
Aulchenko YS, Ripke S, Isaacs A, Van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23 (10): 1294-1296. 10.1093/bioinformatics/btm108.
Sargolzaei M, Chesnais J, Schenkel F: Efficient combined family and population imputation in large data sets. Open Ind Session. 2012, http://www.cdn.ca/Articles/GEBOCT2012/Mehdi_FImputeV2.1.pdf,
Svishcheva GR, Axenovich TI, Belonogova NM, van Duijn CM, Aulchenko YS: Rapid variance components-based method for whole-genome association analysis. Nat Genet. 2012, 44: 1166-1170. 10.1038/ng.2410.
Aulchenko YS, de Koning DJ, Haley C: Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis. Genet. 2007, 177: 577-585. 10.1534/genetics.107.075614.
Devlin B, Roeder K, Wasserman L: Genomic control for association studies.: a semiparametric test to detect excess-haplotype sharing. Biostat. 1999, 55: 997-1004.
Thompson EA, Shaw RG: Pedigree analysis for quantitative traits: variance components without matrix inversion. Biomet. 1990, 46: 399-413. 10.2307/2531445.
Hu ZL, Park CA, Wu XL, Reecy JM: Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era. Nucleic Acids Res. 2013, 41: 871-879. 10.1093/nar/gks1150.
The study was funded in part by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, 2011/16643-2, #2012/02039-9 and 2012/03551-5). The contributions of Núcleo de Criadores de Nelore do Norte do Paraná, Luciano Borges (Rancho da Matinha) and Eduardo Penteado Cardoso (Faz. Mundo Novo) are gratefully acknowledged.
JBSF coordinated and designed the study. MHAS, YTU, RCG, JFG, HF, SLL, PRL, TGL and JBSF contributed in study design and obtaining phenotypes. YTU, MHAS and HHRN helped in data analysis and imputation. MHAS, PAA, RCG and HF did the extraction and standardization of DNA samples from the animals. LLC and RAB performed the genotyping of animals. MHAS performed data analyzes, phenotype collection and preparation of the manuscript. JBSF, YTU, MHAS, RCG, GAOJ, JFG FVM HF and interpreted results and edited the manuscript. All authors approved the final version of the manuscript.
Electronic supplementary material
Additional file 1: Quantile-quantile plot for the test statistics used in the association analysis for DMI (50K).(TIFF 945 KB)
Additional file 2: Quantile-quantile plot for the test statistics used in the association analysis for DMI (HD).(TIFF 941 KB)
Additional file 3: Quantile-quantile plot for the test statistics used in the association analysis for DMI (HDimp).(TIFF 947 KB)
Additional file 4: Quantile-quantile plot for the test statistics used in the association analysis for RFI (50K).(TIFF 943 KB)
Additional file 5: Quantile-quantile plot for the test statistics used in the association analysis for RFI (HD).(TIFF 940 KB)
Additional file 6: Quantile-quantile plot for the test statistics used in the association analysis for RFI (HDimp).(TIFF 943 KB)
About this article
Cite this article
Santana, M.H., Utsunomiya, Y.T., Neves, H.H. et al. Genome-wide association analysis of feed intake and residual feed intake in Nellore cattle. BMC Genet 15, 21 (2014). https://doi.org/10.1186/1471-2156-15-21