Genetic architecture of gene expression in ovine skeletal muscle

BackgroundIn livestock populations the genetic contribution to muscling is intensively monitored in the progeny of industry sires and used as a tool in selective breeding programs. The genes and pathways conferring this genetic merit are largely undefined. Genetic variation within a population has potential, amongst other mechanisms, to alter gene expression via cis- or trans-acting mechanisms in a manner that impacts the functional activities of specific pathways that contribute to muscling traits. By integrating sire-based genetic merit information for a muscling trait with progeny-based gene expression data we directly tested the hypothesis that there is genetic structure in the gene expression program in ovine skeletal muscle.ResultsThe genetic performance of six sires for a well defined muscling trait, longissimus lumborum muscle depth, was measured using extensive progeny testing and expressed as an Estimated Breeding Value by comparison with contemporary sires. Microarray gene expression data were obtained for longissimus lumborum samples taken from forty progeny of the six sires (4-8 progeny/sire). Initial unsupervised hierarchical clustering analysis revealed strong genetic architecture to the gene expression data, which also discriminated the sire-based Estimated Breeding Value for the trait. An integrated systems biology approach was then used to identify the major functional pathways contributing to the genetics of enhanced muscling by using both Estimated Breeding Value weighted gene co-expression network analysis and a differential gene co-expression network analysis. The modules of genes revealed by these analyses were enriched for a number of functional terms summarised as muscle sarcomere organisation and development, protein catabolism (proteosome), RNA processing, mitochondrial function and transcriptional regulation.ConclusionsThis study has revealed strong genetic structure in the gene expression program within ovine longissimus lumborum muscle. The balance between muscle protein synthesis, at the levels of both transcription and translation control, and protein catabolism mediated by regulated proteolysis is likely to be the primary determinant of the genetic merit for the muscling trait in this sheep population. There is also evidence that high genetic merit for muscling is associated with a fibre type shift toward fast glycolytic fibres. This study provides insight into mechanisms, presumably subject to strong artificial selection, that underpin enhanced muscling in sheep populations.

[1]  P. Buttery,et al.  Effects of trenbolone acetate and zeranol on protein metabolism in male castrate andfemale lambs , 1983, British Journal of Nutrition.

[2]  R. Scarpulla,et al.  Activation of the human mitochondrial transcription factor A gene by nuclear respiratory factors: a potential regulatory link between nuclear and mitochondrial gene expression in organelle biogenesis. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[3]  N. Fogarty,et al.  Genetic parameters for ultrasound fat depth and eye muscle measurements in live Poll Dorset sheep , 1994 .

[4]  P. Jongen,et al.  Ca2+ homeostasis in Brody's disease. A study in skeletal muscle and cultured muscle cells and the effects of dantrolene an verapamil. , 1994, The Journal of clinical investigation.

[5]  Alan J. McComas,et al.  Skeletal Muscle: Form and Function , 1996 .

[6]  Michele Harms,et al.  Skeletal Muscle — Form and Function , 1997 .

[7]  V. H. Oddy,et al.  EFFECT OF DIVERGENT SELECTION FOR YEARLING GROWTH RATE ON PROTEIN METABOLISM IN HIND-LIMB MUSCLE AND WHOLE BODY OF ANGUS CATTLE , 1998 .

[8]  B. Spiegelman,et al.  Cross-regulation of C/EBP alpha and PPAR gamma controls the transcriptional pathway of adipogenesis and insulin sensitivity. , 1999, Molecular cell.

[9]  B. Spiegelman,et al.  Cross-Regulation of C/EBPα and PPARγ Controls the Transcriptional Pathway of Adipogenesis and Insulin Sensitivity , 1999 .

[10]  A. Simard,et al.  Calcineurin and skeletal muscle growth , 2002, Nature Cell Biology.

[11]  A. Barabasi,et al.  Hierarchical Organization of Modularity in Metabolic Networks , 2002, Science.

[12]  Brad T. Sherman,et al.  DAVID: Database for Annotation, Visualization, and Integrated Discovery , 2003, Genome Biology.

[13]  Zhijin Wu,et al.  Preprocessing of oligonucleotide array data , 2004, Nature Biotechnology.

[14]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[15]  S. Horvath,et al.  A General Framework for Weighted Gene Co-Expression Network Analysis , 2005, Statistical applications in genetics and molecular biology.

[16]  Eric E Schadt,et al.  Cis-acting expression quantitative trait loci in mice. , 2005, Genome research.

[17]  C. Haley,et al.  Genetical genomics in humans and model organisms. , 2005, Trends in genetics : TIG.

[18]  Jihua Fan,et al.  Runx1 prevents wasting, myofibrillar disorganization, and autophagy of skeletal muscle. , 2005, Genes & development.

[19]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[20]  L. Kruglyak,et al.  Genetics of global gene expression , 2006, Nature Reviews Genetics.

[21]  Michael Watson,et al.  CoXpress: differential co-expression in gene expression data , 2006, BMC Bioinformatics.

[22]  A. E. Rossi,et al.  Sarcoplasmic reticulum: The dynamic calcium governor of muscle , 2006, Muscle & nerve.

[23]  Rainer Breitling,et al.  A verification protocol for the probe sequences of Affymetrix genome arrays reveals high probe accuracy for studies in mouse, human and rat , 2007, BMC Bioinformatics.

[24]  P. Greenwood,et al.  Influences on the loin and cellular characteristics of the m. longissimus lumborum of Australian Poll Dorset-sired lambs , 2006 .

[25]  Eric P. Hoffman,et al.  An interactive power analysis tool for microarray hypothesis testing and generation , 2006, Bioinform..

[26]  C. Lien,et al.  A threshold of GATA4 and GATA6 expression is required for cardiovascular development. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Crispin J. Miller,et al.  Cell Culture , 2010, Cell.

[28]  B. Black,et al.  Determinants of Myogenic Specificity within MyoD Are Required for Noncanonical E Box Binding , 2007, Molecular and Cellular Biology.

[29]  T. Vuocolo,et al.  Identification of a gene network contributing to hypertrophy in callipyge skeletal muscle. , 2007, Physiological genomics.

[30]  Michael T. McManus,et al.  Essential role for Dicer during skeletal muscle development. , 2007, Developmental biology.

[31]  T. Vuocolo,et al.  Analysis of gene expression during the onset of muscle hypertrophy in callipyge lambs. , 2007, Animal genetics.

[32]  R. Alberts,et al.  Sequence Polymorphisms Cause Many False cis eQTLs , 2007, PloS one.

[33]  Chung F. Wong,et al.  MicroRNA-26a Targets the Histone Methyltransferase Enhancer of Zeste homolog 2 during Myogenesis* , 2008, Journal of Biological Chemistry.

[34]  Bin Zhang,et al.  Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R , 2008, Bioinform..

[35]  H. Kadarmideen Genetical systems biology in livestock: application to gonadotrophin releasing hormone and reproduction. , 2008, IET systems biology.

[36]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[37]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[38]  G. Owens,et al.  Kruppel-like factor 4, Elk-1, and histone deacetylases cooperatively suppress smooth muscle cell differentiation markers in response to oxidized phospholipids. , 2008, American journal of physiology. Cell physiology.

[39]  T. Gulick,et al.  Nuclear Respiratory Factor 1 Controls Myocyte Enhancer Factor 2A Transcription to Provide a Mechanism for Coordinate Expression of Respiratory Chain Subunits* , 2008, Journal of Biological Chemistry.

[40]  Edward R B McCabe,et al.  Weighted gene co-expression network analysis identifies biomarkers in glycerol kinase deficient mice. , 2009, Molecular genetics and metabolism.

[41]  Danielle G. Lemay,et al.  Unlocking the bovine genome , 2009, BMC Genomics.

[42]  J. Mallinson,et al.  Chapter 2. Calcineurin signaling and the slow oxidative skeletal muscle fiber type. , 2009, International review of cell and molecular biology.

[43]  S. Dietrich,et al.  An evolutionarily conserved Myostatin proximal promoter/enhancer confers basal levels of transcription and spatial specificity in vivo , 2009, Development Genes and Evolution.

[44]  Charles A Tilford,et al.  Gene set enrichment analysis. , 2009, Methods in molecular biology.

[45]  K. Worley,et al.  The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution , 2009, Science.

[46]  Guoping Fan,et al.  Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonic stem cells , 2009, BMC Genomics.

[47]  D. Freyssenet,et al.  A New Role for Sterol Regulatory Element Binding Protein 1 Transcription Factors in the Regulation of Muscle Mass and Muscle Cell Differentiation , 2009, Molecular and Cellular Biology.

[48]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[49]  Y. E. Chen,et al.  Smooth and Cardiac Muscle-selective Knock-out of Krüppel-like Factor 4 Causes Postnatal Death and Growth Retardation* , 2010, The Journal of Biological Chemistry.

[50]  E. Krüger,et al.  Proteasomal degradation is transcriptionally controlled by TCF11 via an ERAD-dependent feedback loop. , 2010, Molecular cell.

[51]  C. Farber,et al.  Identification of a gene module associated with BMD through the integration of network analysis and genome‐wide association data , 2010, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[52]  A. G. de la Fuente From 'differential expression' to 'differential networking' - identification of dysfunctional regulatory networks in diseases. , 2010, Trends in genetics : TIG.

[53]  Zhou Du,et al.  agriGO: a GO analysis toolkit for the agricultural community , 2010, Nucleic Acids Res..

[54]  T. Vuocolo,et al.  A gene network switch enhances the oxidative capacity of ovine skeletal muscle during late fetal development , 2010, BMC Genomics.

[55]  K. Wimmers,et al.  Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits , 2010, BMC Genomics.

[56]  Antonio Reverter,et al.  Genome-wide patterns of promoter sharing and co-expression in bovine skeletal muscle , 2011, BMC Genomics.

[57]  Jie Chen,et al.  MicroRNAs in skeletal myogenesis , 2011, Cell cycle.

[58]  Haja N Kadarmideen,et al.  Systems biology of ovine intestinal parasite resistance: disease gene modules and biomarkers. , 2011, Molecular bioSystems.

[59]  H. Asahara,et al.  Cellular and Molecular Life Sciences REVIEW The myogenic transcriptional network , 2022 .