RNAseq analysis of fast skeletal muscle in restriction-fed transgenic coho salmon (Oncorhynchus kisutch): an experimental model uncoupling the growth hormone and nutritional signals regulating growth

BMC Genomics

Ian A Johnston²  Robert H Devlin¹  Daniel Garcia de la serrana² 

[1] Department of Fisheries and Oceans, Centre for Aquaculture and Environmental Research, 4160 Marine Drive, West Vancouver V7V 1N6, BC, Canada;Scottish Oceans Institute, School of Biology, University of St Andrews, St Andrews, KY16 8LB, Scotland, UK

关键词: Fish nutrition;    Skeletal muscle transcriptome;    Growth hormone transgenics;    teleost fish;    Muscle growth;   

Others  :  1222457
DOI  :  10.1186/s12864-015-1782-z

 received in 2015-03-31, accepted in 2015-07-15,  发布年份 2015

【 摘 要 】

Background

Coho salmon (Oncorhynchus kisutch) transgenic for growth hormone (Gh) express Gh in multiple tissues which results in increased appetite and continuous high growth with satiation feeding. Restricting Gh-transgenics to the same lower ration (TR) as wild-type fish (WT) results in similar growth, but with the recruitment of fewer, larger diameter, muscle skeletal fibres to reach a given body size. In order to better understand the genetic mechanisms behind these different patterns of muscle growth and to investigate how the decoupling of Gh and nutritional signals affects gene regulation we used RNA-seq to compare the fast skeletal muscle transcriptome in TR and WT coho salmon.

Results

Illumina sequencing of individually barcoded libraries from 6 WT and 6 TR coho salmon yielded 704,550,985 paired end reads which were used to construct 323,115 contigs containing 19,093 unique genes of which >10,000 contained >90 % of the coding sequence. Transcripts coding for 31 genes required for myoblast fusion were identified with 22 significantly downregulated in TR relative to WT fish, including 10 (vaspa, cdh15, graf1, crk, crkl, dock1, trio, plekho1a, cdc42a and dock5) associated with signaling through the cell surface protein cadherin. Nineteen out of 44 (43 %) translation initiation factors and 14 of 47 (30 %) protein chaperones were upregulated in TR relative to WT fish.

Conclusions

TR coho salmon showed increased growth hormone transcripts and gene expression associated with protein synthesis and folding than WT fish even though net rates of protein accretion were similar. The uncoupling of Gh and amino acid signals likely results in additional costs of transcription associated with protein turnover in TR fish. The predicted reduction in the ionic costs of homeostasis in TR fish associated with increased fibre size were shown to involve multiple pathways regulating myotube fusion, particularly cadherin signaling.

【 授权许可】

   
2015 de la serrana et al.

【 图 表 】

Fig. 1.

Fig. 2.

Fig. 3.

Fig. 4.

Fig. 5.

Fig. 6.

Fig. 7.

Fig. 8.

【 参考文献 】

• [1]Bjornsson BT, Johanson V, Benedet S, Einarsdottir IE, Hildahl J, Agustsson T, Jonsson E. Growth hormone endocrinology of salmonids: regulatory mechanism and mode of action. Fish Physiol Biochem. 2002; 27:227-242.
• [2]Mori T, Devlin RH. Transgene and host growth hormone gene expression in pituitary and nonpituitary tissues of normal and growth hormone transgenic salmon. Mol Cel Endocrinol. 1999; 149:129-139.
• [3]Raven PA, Sakhrani MU, Beckman BR, Cooper K, Pinter J, Leder EH. Endocrine effects of growth hormone overexpression in transgenic coho salmon. Gen Comp Endocrinol. 2008; 159:26-37.
• [4]Devlin BH, Timothy TT, Edward MD, Choy-Leong H. Transmission and phenotypic effect of an antifreeze/GH gene construct in coho salmon (Oncorhynchus kisutch). Aquaculture. 1995; 137:161-169.
• [5]Devlin BH, D’Andrade M, Uh M, Biagi CA. Population effects of growth hormone transgenic coho salmon depend on food availability and genotype by environment interactions. Proc Natl Acad Sci U S A. 2004; 101:9303-9308.
• [6]Devlin RH, Sakhrani D, Tymchuk WE, Rise ML, Goh B. Domestication and growth hormone transgenesis cause similar changes in gene expression in coho salmon (Oncorhynchus kisutch). Proc Natl Acad Sci U S A. 2009; 106:3047-3052.
• [7]Bower NI, Johnston IA. Transcriptional regulation of the IGF signaling pathway by amino acids and insulin-like growth factors during myogenesis in Atlantic salmon. PLoS ONE. 2010; 42: Article ID e11100
• [8]Garcia de la serrana D, Johnston IA. Expression of heat shock protein 90 (Hsp90) paralogues is regulated by amino acids in skeletal muscle of Atlantic salmon. PLoS ONE. 2013; 8:e74295.
• [9]Rowlerson A, Mascarello F, Radaelli G, Veggetti A. Differentiation and growth of muscle in the fish Sparus aurata (L): II. Hyperplastic and hypertrophic growth of lateral muscle from hatching to adult. J Mus Res Cell Mot. 1995; 3:223-236.
• [10]Johnston IA, Cole NJ, Abercromby M, Vieira VLA. Embryonic temperature modulates muscle growth characteristics in larval and juvenile herring. J Exp Biol. 1998; 201:623-646.
• [11]Weatherley AH, Gill HS, Lobo AF. Recruitment and maximum diameter of axial muscle fibres in teleost and their relationship to somatic growth and ultimate size. J Fish Biol. 1988; 33:851-859.
• [12]Hill JA, Kiessling A, Devlin RH. Coho salmon (Oncorhynchus kisutch) transgenic for growth hormone gene construct exhibit increased rates of hyperplasia and detectable levels of differential gene expression. Can J Aquat Sci. 2000; 57:939-950.
• [13]Johnston IA, Garcia de la serrana D, Devlin RH. Muscle fibre size optimisation provides flexibility to energy budgeting in caloric-restricted Coho salmon transgenic for growth hormone. J Exp Biol. 2014; 217:3392-95.
• [14]Johnston IA, Kristjánsson BK, Paxton CG, Vieira VLA, Macqueen DJ, Bell MA. Universal scaling rules predict evolutionary patterns of myogenesis is species with indeterminate growth. Proc R Sco B. 2012; 279:2255-61.
• [15]Jimenez AG, Dillaman RM, Kinsey ST. Large fibre size in skeletal muscle is metabolically advantageous. Nature Comm. 2013; 4:2150.
• [16]Glasauer SMK, Neuhauss SCF. Whole-genome duplication in teleost fishes and its evolutionary consequences. Mol Gen Genom. 2014; 289:1045-1060.
• [17]Garcia de la serrana D, Mareco EA, Johnston IA. Systematic variation in the pattern of gene paralog retention between the teleost superorders Ostariophysi and Acanthopterygii. Gen Biol Evol. 2014; 6:981-7.
• [18]Macqueen DJ, Johnston IA. A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification. Proc R Soc B Biol Scie. 2014; 281:20132881.
• [19]Macqueen DJ, Johnston IA. Evolution of the multifaceted eukaryotic akirin gene family. BMC Evol Biol. 2009; 9:34. BioMed Central Full Text
• [20]Macqueen DJ, Garcia de la serrana D, Johnston IA. Evolution of ancient functions in the vertebrate insulin-like growth factor system uncovered by study of duplicated salmonid fish genomes. Mol Biol Evol. 2013; 30:1060-1076.
• [21]Lindholm ME, Huss M, Solnestam BW, Kjellqvist S, Lunderberg J, Sundberg C. The human skeletal muscle transcriptome: sex differences, alternative splicing, and tissue homogeneity assessed with RNA sequencing. FASEB J. 2014; 28:4571-81.
• [22]Choo AY, Yoon SO, Kim SG, Roux PP, Blenis J. Rapamycin differentially inhibits S6Ks and E4-BP1 to mediate cell-type-specific repression of mRNA translation. Proc Natl Acad Sci. 2008; 105:17414-17419.
• [23]Gringas AC, Raught B, Gygi SP, Niedzwiecka A, Miron M, Burley SK, et al. Hierarchical phosphorylation of the translation inhibitor 4E-BP1. Gen Dev. 2001;15:2852–64.
• [24]Huang K, Fingar DC. Growing knowledge of the mTOR signaling network. Sem Cell Dev Biol. 2014; 36:79-90.
• [25]Sancak Y, Bar-Peled L, Zoncu R, Markhard AL, Nada S, Sabatini DM. Ragulator-Rag compex target mTORC1 to the lysosomal surface and is necessary for its activation by amino acids. Cell. 2010; 141:290-303.
• [26]Roccio M, Bos JL, Zwarkruis FJT. Regulation of the small GTPase Rheb by amino acids. Oncogene. 2006; 25:657-664.
• [27]Du J, Bernick E. Molecular chaperones Hsp90a1 and Unc45b are required for sarcomere assembly in skeletal muscles of zebrafish embryos. FASEB. 2011; 28:10.
• [28]Hindi SM, Marjan M, Kumar A. Signaling Mechanisms in Mammalian Myoblast Fusion. Sci Signal. 2013; 6:272.
• [29]Moore CA, Parkin CA, Bidet Y, Ingham PW. A role for the Myoblast city homologues Dock1 and Dock5 and the adaptors Crk and Crk-like in zebrafish myoblast fusion. Develop. 2007; 134:3145-3153.
• [30]Charrasse S, Comunale F, Grumbach Y, Poulat F, Blangy A, Gauthier-Rouviére C. RhoA GTPase regulates m-cahderin activity and myoblast fusion. Mol Biol Cell. 2006; 17:749-759.
• [31]Millay DP, O’Rourke JR, Sutherland LB, Bezprozvannya S, Shelton JM, Bassel-Duby R, et al. Myomaker is a membrane activator of myoblast fusion and muscle formation. Nature. 2013;499:301–5.
• [32]Quach NL, Biressi S, Reichardt LF, Keller C, Rando TA. Focal adhesion kinase signalling regulates the expression of caveolin 3 and beta integrin, genes essential for normal myoblast fusion. Mol Biol Cell. 2009; 14:3422-3455.
• [33]Devlin RH, Biagi CA, Yesaki TY. Growth, viability and genetic characteristics of GH transgenic coho salmon strains. Aquaculture. 2004; 236:607-632.
• [34]Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genom Biol. 2002;3:research0034–research0034.11.
• [35]Hass BJ, Papanicolau A, Yassour M, Grabherr M, Blood PD, Bowden, et al. De novo transcript sequence reconstruction from RNA-seq: reference generation and analysis with Trinity. Nat Protoc. 2013;8:1494–512.
• [36]Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucl Acids. 2008;36:3420–35.
• [37]Moriya I, Itoh M, Okuda S, Yoshizama A, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucl Acids. 2007; 35:182-185.
• [38]www.ensembl.org (last entrance March 2014).
• [39]Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinformatics. 2011; 12:323. BioMed Central Full Text
• [40]Anders S, Huber W. Differential expression analysis of sequence count data. Genom Biol. 2010; 11:106. BioMed Central Full Text
• [41]Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increase coverage and integration. Nucl Acids. 2013;41:808–15.
• [42]Penn O, Privman E, Ashkenazy H, Landan G, Graur D, Pupko T. GUIDANCE: a web server for assessing alignment confidence scores. Nucl Acids Res. 2010; 1:W23-W28.
• [43]Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013; 30:2725-2729.
• [44]Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F et al.. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucl Acids Res. 2008; 36:465-9.
• [45]Garcia de la serrana D, Estevez A, Andree K, Johnston IA. Fast skeletal muscle transcriptome of the Gilthead sea bream (Sparus aurata) determined by next generation sequencing. BMC Genomics. 2012; 13:181. BioMed Central Full Text
• [46]Mareco EA, Garcia de la serrana D, Johnston IA, Dal-Pai-Silva M. Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus). BMC Genomics. 2015.
• [47]Caraux G, Pinloche S, Permutmatrix. A graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005; 21:1280-1281.