Co-Expression Analysis of Fetal Weight-Related Genes in Ovine Skeletal Muscle during Mid and Late Fetal Development Stages

Background: Muscle development and lipid metabolism play important roles during fetal development stages. The commercial Texel sheep are more muscular than the indigenous Ujumqin sheep. Results: We performed serial transcriptomics assays and systems biology analyses to investigate the dynamics of gene expression changes associated with fetal longissimus muscles during different fetal stages in two sheep breeds. Totally, we identified 1472 differentially expressed genes during various fetal stages using time-series expression analysis. A systems biology approach, weighted gene co-expression network analysis (WGCNA), was used to detect modules of correlated genes among these 1472 genes. Dramatically different gene modules were identified in four merged datasets, corresponding to the mid fetal stage in Texel and Ujumqin sheep, the late fetal stage in Texel and Ujumqin sheep, respectively. We further detected gene modules significantly correlated with fetal weight, and constructed networks and pathways using genes with high significances. In these gene modules, we identified genes like TADA3, LMNB1, TGF-β3, EEF1A2, FGFR1, MYOZ1, and FBP2 correlated with fetal weight. Conclusion: Our study revealed the complex network characteristics involved in muscle development and lipid metabolism during fetal development stages. Diverse patterns of the network connections observed between breeds and fetal stages could involve some hub genes, which play central roles in fetal development, correlating with fetal weight. Our findings could provide potential valuable biomarkers for selection of body weight-related traits in sheep and other livestock.

[1]  Y. Kataoka,et al.  Differential expression of nuclear lamin, the major component of nuclear lamina, during neurogenesis in two germinal regions of adult rat brain , 2007, The European journal of neuroscience.

[2]  Aiqing He,et al.  Identification of inflammatory gene modules based on variations of human endothelial cell responses to oxidized lipids , 2006, Proceedings of the National Academy of Sciences.

[3]  Antonio Reverter,et al.  A gene coexpression network for bovine skeletal muscle inferred from microarray data. , 2006, Physiological genomics.

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

[5]  B. Dave,et al.  Alteration/deficiency in activation-3 (Ada3) plays a critical role in maintaining genomic stability , 2012, Cell cycle.

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

[7]  A. Reverter,et al.  Gene expression studies of developing bovine longissimus muscle from two different beef cattle breeds , 2007, BMC Developmental Biology.

[8]  Leif Andersson,et al.  Domestic-animal genomics: deciphering the genetics of complex traits , 2004, Nature Reviews Genetics.

[9]  Goberdhan P Dimri,et al.  An Essential Role of Human Ada3 in p53 Acetylation* , 2007, Journal of Biological Chemistry.

[10]  G. Breur,et al.  FGFR1 inhibits skeletal muscle atrophy associated with hindlimb suspension , 2007, BMC musculoskeletal disorders.

[11]  J. Wen,et al.  Identification of differentially expressed genes and pathways for intramuscular fat deposition in pectoralis major tissues of fast-and slow-growing chickens , 2012, BMC Genomics.

[12]  M. Bergo,et al.  A carboxyl-terminal interaction of lamin B1 is dependent on the CAAX endoprotease Rce1 and carboxymethylation , 2003, The Journal of cell biology.

[13]  S. McWilliam,et al.  Skeletal muscle specific genes networks in cattle , 2010, Functional & Integrative Genomics.

[14]  H. Newbery,et al.  Haploinsufficiency for Translation Elongation Factor eEF1A2 in Aged Mouse Muscle and Neurons Is Compatible with Normal Function , 2012, PloS one.

[15]  Se-Jin Lee,et al.  Regulation of myostatin activity and muscle growth , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[16]  R. Derynck,et al.  Differentiation plasticity regulated by TGF-β family proteins in development and disease , 2007, Nature Cell Biology.

[17]  E. Androphy,et al.  hAda3 regulates p14ARF-induced p53 acetylation and senescence , 2007, Oncogene.

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

[19]  R. Derynck,et al.  Differentiation plasticity regulated by TGF-|[beta]| family proteins in development and disease , 2007 .

[20]  Joshua M. Stuart,et al.  A Gene-Coexpression Network for Global Discovery of Conserved Genetic Modules , 2003, Science.

[21]  M. Gerstein,et al.  Genomic analysis of regulatory network dynamics reveals large topological changes , 2004, Nature.

[22]  C. Berri,et al.  Muscle fibre ontogenesis in farm animal species. , 2002, Reproduction, nutrition, development.

[23]  Ting Wang,et al.  hADA3 is required for p53 activity , 2001, The EMBO journal.

[24]  Qianjun Zhao,et al.  [Variation of MSTN gene UTR in eleven sheep breeds]. , 2009, Yi chuan = Hereditas.

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

[26]  D. Birnbaum,et al.  Endogenous retroviral sequence is fused to FGFR1 kinase in the 8p12 stem-cell myeloproliferative disorder with t(8;19)(p12;q13.3). , 2003, Blood.

[27]  Peter Langfelder,et al.  Weighted gene co-expression network analysis of the peripheral blood from Amyotrophic Lateral Sclerosis patients , 2009, BMC Genomics.

[28]  Advances in livestock genomics: opening the barn door. , 2005 .

[29]  Peter Langfelder,et al.  Eigengene networks for studying the relationships between co-expression modules , 2007, BMC Systems Biology.

[30]  C. Guda,et al.  Mammalian Alteration/Deficiency in Activation 3 (Ada3) Is Essential for Embryonic Development and Cell Cycle Progression* , 2012, The Journal of Biological Chemistry.

[31]  T. P. Smith,et al.  Mutations in myostatin (GDF8) in double-muscled Belgian Blue and Piedmontese cattle. , 1997, Genome research.

[32]  V. Kaartinen,et al.  Abnormal lung development and cleft palate in mice lacking TGF–β3 indicates defects of epithelial–mesenchymal interaction , 1995, Nature Genetics.

[33]  Joaquín Dopazo,et al.  Serial Expression Analysis: a web tool for the analysis of serial gene expression data , 2010, Nucleic Acids Res..

[34]  Haja N Kadarmideen,et al.  Building gene co-expression networks using transcriptomics data for systems biology investigations: Comparison of methods using microarray data , 2012, Bioinformation.

[35]  B. Buendia,et al.  Caspase-dependent proteolysis of integral and peripheral proteins of nuclear membranes and nuclear pore complex proteins during apoptosis. , 1999, Journal of cell science.

[36]  J. Hocquette.  Endocrine and metabolic regulation of muscle growth and body composition in cattle. , 2010, Animal : an international journal of animal bioscience.

[37]  H. Su,et al.  Histological and transcriptome-wide level characteristics of fetal myofiber hyperplasia during the second half of gestation in Texel and Ujumqin sheep , 2011, BMC Genomics.

[38]  E. Koonin,et al.  Conservation and coevolution in the scale-free human gene coexpression network. , 2004, Molecular biology and evolution.

[39]  S. Friend,et al.  A network view of disease and compound screening , 2009, Nature Reviews Drug Discovery.

[40]  Yan Wang,et al.  VisANT 3.5: multi-scale network visualization, analysis and inference based on the gene ontology , 2009, Nucleic Acids Res..

[41]  J. Hickford,et al.  Association between a g.+6723G-A SNP in the myostatin gene (MSTN) and carcass traits in New Zealand Texel sheep. , 2009, Journal of animal science.

[42]  David Bryant,et al.  DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists , 2007, Nucleic Acids Res..

[43]  J. McEwan,et al.  Early stages of myogenesis in a large mammal: Formation of successive generations of myotubes in sheep tibialis cranialis muscle , 1992, Journal of Muscle Research & Cell Motility.

[44]  S. Horvath,et al.  Statistical Applications in Genetics and Molecular Biology , 2011 .

[45]  S. Kaufmann,et al.  Differential expression of nuclear envelope lamins A and C in human lung cancer cell lines. , 1991, Cancer research.

[46]  T. Vuocolo,et al.  Genetic architecture of gene expression in ovine skeletal muscle , 2011, BMC Genomics.

[47]  N. da Costa,et al.  Development of a porcine skeletal muscle cDNA microarray: analysis of differential transcript expression in phenotypically distinct muscles , 2003, BMC Genomics.

[48]  C. Cowled,et al.  Genetic architecture of gene expression in the chicken , 2013, BMC Genomics.

[49]  P. von Rohr,et al.  From genetical genomics to systems genetics: potential applications in quantitative genomics and animal breeding , 2006, Mammalian Genome.

[50]  Yang Li,et al.  Differential transcriptional analysis between red and white skeletal muscle of Chinese Meishan pigs , 2010, International journal of biological sciences.

[51]  S. Horvath,et al.  Weighted gene coexpression network analysis strategies applied to mouse weight , 2007, Mammalian Genome.

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

[53]  Shu-lin Yang,et al.  Characterization of different expression patterns of calsarcin-1 and calsarcin-2 in porcine muscle. , 2006, Gene.

[54]  J. McEwan,et al.  Myosin heavy chain composition of single fibres and their origins and distribution in developing fascicles of sheep tibialis cranialis muscles , 1992, Journal of Muscle Research & Cell Motility.

[55]  S. Horvath,et al.  Evidence for anti-Burkitt tumour globulins in Burkitt tumour patients and healthy individuals. , 1967, British Journal of Cancer.

[56]  H. Newbery,et al.  Translation Elongation Factor eEF1A2 Is Essential for Post-weaning Survival in Mice* , 2007, Journal of Biological Chemistry.

[57]  Michel Georges,et al.  A deletion in the bovine myostatin gene causes the double–muscled phenotype in cattle , 1997, Nature Genetics.

[58]  A. Barabasi,et al.  Network biology: understanding the cell's functional organization , 2004, Nature Reviews Genetics.

[59]  J. McEwan,et al.  Meat quality traits were unaffected by a quantitative trait locus affecting leg composition traits in Texel sheep. , 2005, Journal of animal science.

[60]  Ana Conesa,et al.  Gene expression maSigPro : a method to identify significantly differential expression profiles in time-course microarray experiments , 2006 .

[61]  S. Scholpp,et al.  Zebrafish fgfr1 is a member of the fgf8 synexpression group and is required for fgf8 signalling at the midbrain-hindbrain boundary , 2004, Development, Genes and Evolution.

[62]  E. Olson,et al.  Calsarcin-2 deficiency increases exercise capacity in mice through calcineurin/NFAT activation. , 2008, The Journal of clinical investigation.

[63]  T. Ideker,et al.  Integrative approaches for finding modular structure in biological networks , 2013, Nature Reviews Genetics.

[64]  E. Olson,et al.  Calsarcins, a novel family of sarcomeric calcineurin-binding proteins. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[65]  Biochemical and transcriptomic analyses of two bovine skeletal muscles in Charolais bulls divergently selected for muscle growth. , 2005, Meat science.

[66]  Antonio Reverter,et al.  Inferring the Transcriptional Landscape of Bovine Skeletal Muscle by Integrating Co-Expression Networks , 2009, PloS one.

[67]  Florian Caiment,et al.  A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep , 2006, Nature Genetics.

[68]  M. Gustin,et al.  Functional meta-analysis of double connectivity in gene coexpression networks in mammals. , 2008, Physiological genomics.

[69]  S. Gordon,et al.  Differential global gene expression in red and white skeletal muscle. , 2001, American Journal of Physiology - Cell Physiology.

[70]  A. Rizzino,et al.  Expression of transforming growth factor-β3 by embryonal carcinoma cells, parietal endoderm-like cells and early mouse embryos , 1990, In Vitro Cellular & Developmental Biology.

[71]  Marc S. Williams.  Myostatin mutation associated with gross muscle hypertrophy in a child. , 2004, The New England journal of medicine.

[72]  C. Bustamante,et al.  A Mutation in the Myostatin Gene Increases Muscle Mass and Enhances Racing Performance in Heterozygote Dogs , 2007, PLoS genetics.

[73]  Antonio Reverter,et al.  A Differential Wiring Analysis of Expression Data Correctly Identifies the Gene Containing the Causal Mutation , 2009, PLoS Comput. Biol..