Multi-Trait GWAS and New Candidate Genes Annotation for Growth Curve Parameters in Brahman Cattle

Understanding the genetic architecture of beef cattle growth cannot be limited simply to the genome-wide association study (GWAS) for body weight at any specific ages, but should be extended to a more general purpose by considering the whole growth trajectory over time using a growth curve approach. For such an approach, the parameters that are used to describe growth curves were treated as phenotypes under a GWAS model. Data from 1,255 Brahman cattle that were weighed at birth, 6, 12, 15, 18, and 24 months of age were analyzed. Parameter estimates, such as mature weight (A) and maturity rate (K) from nonlinear models are utilized as substitutes for the original body weights for the GWAS analysis. We chose the best nonlinear model to describe the weight-age data, and the estimated parameters were used as phenotypes in a multi-trait GWAS. Our aims were to identify and characterize associated SNP markers to indicate SNP-derived candidate genes and annotate their function as related to growth processes in beef cattle. The Brody model presented the best goodness of fit, and the heritability values for the parameter estimates for mature weight (A) and maturity rate (K) were 0.23 and 0.32, respectively, proving that these traits can be a feasible alternative when the objective is to change the shape of growth curves within genetic improvement programs. The genetic correlation between A and K was -0.84, indicating that animals with lower mature body weights reached that weight at younger ages. One hundred and sixty seven (167) and two hundred and sixty two (262) significant SNPs were associated with A and K, respectively. The annotated genes closest to the most significant SNPs for A had direct biological functions related to muscle development (RAB28), myogenic induction (BTG1), fetal growth (IL2), and body weights (APEX2); K genes were functionally associated with body weight, body height, average daily gain (TMEM18), and skeletal muscle development (SMN1). Candidate genes emerging from this GWAS may inform the search for causative mutations that could underpin genomic breeding for improved growth rates.

[1]  E. J. Pollak,et al.  QTLs associated with dry matter intake, metabolic mid-test weight, growth and feed efficiency have little overlap across 4 beef cattle studies , 2014, BMC Genomics.

[2]  S. Lagarrigue,et al.  Re-Sequencing Data for Refining Candidate Genes and Polymorphisms in QTL Regions Affecting Adiposity in Chicken , 2014, PloS one.

[3]  L. Kiemeney,et al.  A Comparison of Multivariate Genome-Wide Association Methods , 2014, PloS one.

[4]  F. Schenkel,et al.  Genome-Wide Association for Growth Traits in Canchim Beef Cattle , 2014, PloS one.

[5]  D. Riley,et al.  Mapping genomic markers to closest feature using the R package Map2NCBI , 2014 .

[6]  M. Goddard,et al.  A Multi-Trait, Meta-analysis for Detecting Pleiotropic Polymorphisms for Stature, Fatness and Reproduction in Beef Cattle , 2014, PLoS genetics.

[7]  José Fernando Garcia,et al.  Genome-wide association analysis of feed intake and residual feed intake in Nellore cattle , 2014, BMC Genetics.

[8]  B. Penninx,et al.  Resequencing Three Candidate Genes for Major Depressive Disorder in a Dutch Cohort , 2013, PloS one.

[9]  S. Rodriguez-Zas,et al.  Bivariate Genome-Wide Association Analysis of the Growth and Intake Components of Feed Efficiency , 2013, PloS one.

[10]  M. Stephens,et al.  Integrated Enrichment Analysis of Variants and Pathways in Genome-Wide Association Studies Indicates Central Role for IL-2 Signaling Genes in Type 1 Diabetes, and Cytokine Signaling Genes in Crohn's Disease , 2013, PLoS genetics.

[11]  S. Moore,et al.  Genome-wide association analyses for carcass quality in crossbred beef cattle , 2013, BMC Genetics.

[12]  R. Lôbo,et al.  Genetic Analysis of Growth Traits in Polled Nellore Cattle Raised on Pasture in Tropical Region Using Bayesian Approaches , 2013, PloS one.

[13]  Jie Tian,et al.  A Longitudinal Study of Hand Motor Recovery after Sub-Acute Stroke: A Study Combined fMRI with Diffusion Tensor Imaging , 2013, PloS one.

[14]  D. Plenz,et al.  powerlaw: A Python Package for Analysis of Heavy-Tailed Distributions , 2013, PloS one.

[15]  S. Chien,et al.  Involvement of Rab28 in NF-κB Nuclear Transport in Endothelial Cells , 2013, PloS one.

[16]  Yuandan Zhang,et al.  The IGF1 pathway genes and their association with age of puberty in cattle. , 2013, Animal genetics.

[17]  A. Reverter,et al.  Breeding and Genetics Symposium: networks and pathways to guide genomic selection. , 2013, Journal of animal science.

[18]  M. McGowan,et al.  Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits , 2013 .

[19]  M. McGowan,et al.  Male traits and herd reproductive capability in tropical beef cattle. 1. Experimental design and animal measures , 2013 .

[20]  Hong Chen,et al.  Novel SNPs in the bovine Transmembrane protein 18 gene, their linkage and their associations with growth traits in Nanyang cattle , 2012, Genes & Genomics.

[21]  Xiaolin Wu,et al.  Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era , 2012, Nucleic Acids Res..

[22]  R. Snell,et al.  Genetic variation in PLAG1 associates with early life body weight and peripubertal weight and growth in Bos taurus. , 2012, Animal genetics.

[23]  K. Meyer,et al.  “SNP Snappy”: A Strategy for Fast Genome-Wide Association Studies Fitting a Full Mixed Model , 2012, Genetics.

[24]  Kiranmoy Das,et al.  Genome-Wide Association Studies for Bivariate Sparse Longitudinal Data , 2011, Human Heredity.

[25]  G. Chrousos,et al.  Association of TMEM18 variants with BMI and waist circumference in children and correlation of mRNA expression in the PFC with body weight in rats , 2011, European Journal of Human Genetics.

[26]  Tom Druet,et al.  Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature , 2011, Nature Genetics.

[27]  Runze Li,et al.  A dynamic model for genome-wide association studies , 2011, Human Genetics.

[28]  F. Schick,et al.  Novel Obesity Risk Loci Do Not Determine Distribution of Body Fat Depots: A Whole‐body MRI/MRS study , 2010, Obesity.

[29]  Robert Fredriksson,et al.  The obesity gene, TMEM18, is of ancient origin, found in majority of neuronal cells in all major brain regions and associated with obesity in severely obese children , 2010, BMC Medical Genetics.

[30]  M. Calus,et al.  QTLMAS 2009: simulated dataset , 2010, BMC proceedings.

[31]  R. Pong-Wong,et al.  A two-step approach combining the Gompertz growth model with genomic selection for longitudinal data , 2010, BMC proceedings.

[32]  G. L. Bennett,et al.  Genome-wide association study of growth in crossbred beef cattle. , 2010, Journal of animal science.

[33]  Cedric Gondro,et al.  Immunological Profiles of Bos taurus and Bos indicus Cattle Infested with the Cattle Tick, Rhipicephalus (Boophilus) microplus , 2009, Clinical and Vaccine Immunology.

[34]  S. Wells,et al.  Sfrp1 and Sfrp2 are required for normal male sexual development in mice. , 2009, Developmental biology.

[35]  L. G. Albuquerque,et al.  Comparison of different nonlinear functions to describe Nelore cattle growth. , 2009, Journal of animal science.

[36]  Ermias Kebreab,et al.  Mathematical Modelling in Animal Nutrition , 2008 .

[37]  L. G. Albuquerque,et al.  Analysis of beef cattle longitudinal data applying a nonlinear model. , 2007, Journal of animal science.

[38]  Graydon B. Gonsalvez,et al.  A Drosophila melanogaster model of spinal muscular atrophy reveals a function for SMN in striated muscle , 2007, The Journal of cell biology.

[39]  R. Wu,et al.  Functional mapping — how to map and study the genetic architecture of dynamic complex traits , 2006, Nature Reviews Genetics.

[40]  F. Casas,et al.  Coactivation of nuclear receptors and myogenic factors induces the major BTG1 influence on muscle differentiation , 2005, Oncogene.

[41]  Mark Daly,et al.  Haploview: analysis and visualization of LD and haplotype maps , 2005, Bioinform..

[42]  M. Ohno,et al.  Growth retardation and dyslymphopoiesis accompanied by G2/M arrest in APEX2-null mice. , 2004, Blood.

[43]  S. Moore,et al.  Identification and fine mapping of quantitative trait loci for growth traits on bovine chromosomes 2, 6, 14, 19, 21, and 23 within one commercial line of Bos taurus. , 2004, Journal of animal science.

[44]  D. Gianola,et al.  Bayesian inference about parameters of a longitudinal trajectory when selection operates on a correlated trait. , 2003, Journal of animal science.

[45]  G. Casella,et al.  Functional mapping of quantitative trait loci underlying the character process: a theoretical framework. , 2002, Genetics.

[46]  S. Gabriel,et al.  The Structure of Haplotype Blocks in the Human Genome , 2002, Science.

[47]  J. Magaud,et al.  Identification of functional domains involved in BTG1 cell localization , 2001, Oncogene.

[48]  J. Melki,et al.  Deletion of Murine SMN Exon 7 Directed to Skeletal Muscle Leads to Severe Muscular Dystrophy , 2001, The Journal of cell biology.

[49]  M. Herrero,et al.  Modeling extended lactations of dairy cows. , 2000, Journal of dairy science.

[50]  R. Maas,et al.  Msx2 deficiency in mice causes pleiotropic defects in bone growth and ectodermal organ formation , 2000, Nature Genetics.

[51]  C. Chase,et al.  Evaluation of growth curves of Brahman cattle of various frame sizes. , 1996, Journal of animal science.

[52]  J. Takahashi,et al.  Postweaning growth of brahman and santa gertrudis steers under feedlots in the subtropics , 1988 .

[53]  R S DeNise,et al.  Genetic and environmental aspects of the growth curve parameters in beef cows. , 1985, Journal of animal science.

[54]  H. Fitzhugh,et al.  Analysis of growth curves and strategies for altering their shape. , 1976, Journal of animal science.

[55]  J. E. Brown,et al.  A Comparison of Nonlinear Models for Describing Weight-Age Relationships in Cattle , 1976 .

[56]  J. Nelder.  The Fitting of a Generalization of the Logistic Curve , 1961 .

[57]  F. J. Richards.  A Flexible Growth Function for Empirical Use , 1959 .

[58]  L. Bertalanffy.  Quantitative Laws in Metabolism and Growth , 1957 .

[59]  Samuel Brody,et al.  Bioenergetics and growth. With special reference to the efficiency complex in domestic animals. , 1946 .

[60]  Peter Dalgaard,et al.  R Development Core Team (2010): R: A language and environment for statistical computing , 2010 .

[61]  C. Haley,et al.  Search for quantitative trait loci affecting growth and carcass traits in a cross population of beef and dairy cattle. , 2009, Journal of animal science.

[62]  Ncbi.  National Center for Biotechnology Information , 2008 .

[63]  M. Goddard,et al.  Genomic selection. , 2007, Journal of animal breeding and genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie.

[64]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[65]  C. Winsor,et al.  The Gompertz Curve as a Growth Curve. , 1932, Proceedings of the National Academy of Sciences of the United States of America.