Changes in skeletal muscle gene expression following clenbuterol administration
© Spurlock et al; licensee BioMed Central Ltd. 2006
Received: 19 May 2006
Accepted: 20 December 2006
Published: 20 December 2006
Beta-adrenergic receptor agonists (BA) induce skeletal muscle hypertrophy, yet specific mechanisms that lead to this effect are not well understood. The objective of this research was to identify novel genes and physiological pathways that potentially facilitate BA induced skeletal muscle growth. The Affymetrix platform was utilized to identify gene expression changes in mouse skeletal muscle 24 hours and 10 days after administration of the BA clenbuterol.
Administration of clenbuterol stimulated anabolic activity, as indicated by decreased blood urea nitrogen (BUN; P < 0.01) and increased body weight gain (P < 0.05) 24 hours or 10 days, respectively, after initiation of clenbuterol treatment. A total of 22,605 probesets were evaluated with 52 probesets defined as differentially expressed based on a false discovery rate of 10%. Differential mRNA abundance of four of these genes was validated in an independent experiment by quantitative PCR. Functional characterization of differentially expressed genes revealed several categories that participate in biological processes important to skeletal muscle growth, including regulators of transcription and translation, mediators of cell-signalling pathways, and genes involved in polyamine metabolism.
Global evaluation of gene expression after administration of clenbuterol identified changes in gene expression and overrepresented functional categories of genes that may regulate BA-induced muscle hypertrophy. Changes in mRNA abundance of multiple genes associated with myogenic differentiation may indicate an important effect of BA on proliferation, differentiation, and/or recruitment of satellite cells into muscle fibers to promote muscle hypertrophy. Increased mRNA abundance of genes involved in the initiation of translation suggests that increased levels of protein synthesis often associated with BA administration may result from a general up-regulation of translational initiators. Additionally, numerous other genes and physiological pathways were identified that will be important targets for further investigations of the hypertrophic effect of BA on skeletal muscle.
Anabolic effects of β-adrenergic receptor agonists (BA) have been widely studied for potential applications in the prevention of muscle atrophy [1, 2] and improvement of the efficiency of muscle growth in production livestock [3–5]. Clenbuterol is a β2-adrenergic receptor agonist that has been shown to have a significant effect on muscle metabolism in a variety of muscle atrophy models, including hind-limb suspension atrophy [6, 7], starvation induced atrophy , and denervation induced atrophy [9, 10]. Additionally, clenbuterol is known to induce a significant repartitioning effect by increasing the growth of skeletal muscle at the expense of fat tissues in most livestock species [3, 5]. Although it is known that clenbuterol initiates these effects via activation of the β2-adrenergic receptor [3, 5], the downstream mechanisms by which activation of these receptors results in increased muscle growth or decreased muscle atrophy are not clear. To date, the expression and activity of specific genes have been investigated in a variety of models in order to implicate specific pathways with the skeletal muscle response to BA. For example, increased abundance of myofibrillar and structural proteins has been demonstrated and appears to result from increases in both transcription and translation of these genes [11–13]. Additionally, endogenous proteinases including genes of the ubiquitin-proteasome pathway and calcium-dependent proteolytic enzymes have been reported to mediate protein turnover in skeletal muscle after administration of BA . Finally, changes in skeletal muscle expression of IGF1 and IGF2 mRNA have been observed shortly after the administration of clenbuterol to rodents , suggesting the regulation of these growth factors may be important in the initial response of skeletal muscle to BA.
Although investigations of candidate genes have been informative in terms of implicating individual pathways in the skeletal muscle response to BA, they have not provided a global view of changes occurring in the tissue, and have been limited to the investigation of genes with known functions. Thus, the objective of the current research was to define changes in the global gene expression profile of skeletal muscle in response to administration of the BA clenbuterol. Two time points relative to clenbuterol administration were investigated in order to compare gene expression profiles following short- (24 hour; 24 h) and long- (10 day; 10 D) term clenbuterol administration in mice. The Affymetrix platform was chosen for the analysis of gene expression in order to investigate the most comprehensive collection of genes available.
Clenbuterol stimulated an anabolic response
Clenbuterol stimulated changes in mRNA abundance
Significantly differentially expressed genes (10% FDR threshold).
Change 24 hc
Change 10 Dd
Change 10 D/24 he
Cyclin-dependent kinase inhibitor 1A (P21)
TSC22 domain family, member 1
Growth hormone receptor
Tubulointerstitial nephritis antigen-like
Ankyrin repeat and SOCS box-containing protein 15
Myeloid/lymphoid or mixed-lineage leukaemia
Eukaryotic translation initiation factor 4A1
Chemokine ligand 12
Actin-binding Rho activating protein
Microfibrilar-associated protein 3-like
Down syndrome critical region gene 1-like 1
YY1 transcription factor
Thimet oligopeptidase 1
RIKEN cDNA 4631423 FO2 gene
Actin related protein 2/3 complex, subunit 3
Tubulin, beta 5
Fructose bisphosphatase 2
Interleukin 6 signal transducer
Enhancer of yellow 2 homolog
Serine (or cysteine) preptidase inhibitor, clade A, member 1b
Contactin associated protein-like 2
Kelch repeat and BTB (POZ) domain containing 5
Myosin, light polypeptide kinase 2, skeletal muscle
Tumor necrosis factor receptor superfamily, member 12a
Nuclear receptor subfamily 4, group A, member 1
ATPase, H+ transporting, V0 subunit E
Ortholog of human Ras association and pleckstrin homolog domains 1
Solute carrier family 7 (cationic amino acid transporter, y+ system), member 5
Homebox only domain
Protein kinase, cGMP-dependent, type 1
Ring finger protein 128
Tubulin, beta 6
S-adenosylmethionine decarboxylase 1
Growth arrest and DNA-damage-inducible 45 gamma
Ribonucleotide reductase M2B 9TP53 inducible
Bromodomain containing 4
Insulin-like growth factor 1
Dehydrogenase/reductase (SDR family) member 7
RIKEN cDNA 1190002C06 gene
TSR 1, 20S rRNA accumulation, homolog
DNA-damage-inducible transcript 4-like
Ankyrin repeat and SOCS box-containing protein 10
Lady bird-like homeobox 1 homolog (Drosophila)
Aryl hydrocarbon receptor nuclear translocator-like
Summary of results from cluster analysis.
Number of Genes
Percentage of Genes
Differentially expressed genes of interest (P < 0.01) categorized according to Gene Ontology biological process terms (level 2).
Gene Ontology Biological Process Categorya
Number of Genes
Cellular physiological process
Regulation of physiological process
Response to stimulus
Organismal physiological process
Regulation of cellular process
Regulation of development
Clenbuterol stimulated changes in families and categories of genes
Selecteda genes of interest.
Change 24 hd
Change 10 De
Change 10 D/24 hf
Ornithine decarboxylase, structural
Insulin-like growth factor binding protein 5
Solute carrier family 7 (cationic amino acid transporter, y+ system), member 2
Insulin-like growth factor 2 receptor
Proteasome (prosome, macropain) subunit, beta type 1
Ubiquitin protein ligase E3 component n-recognin 1
Ankyrin repeat and SOCS box-containing protein 14
Eukaryotic translation initiation factor 2B, subunit 4 delta
Myocyte enhancer factor 2C
Seven in absentia 2
Mitogen activated protein kinase kinase 4
Mitogen activated protein kinase kinase 6
Actinin alpha 4
Eukaryotic translation initiation factor 5
Protein phosphatase 1 (formerly 2C)-like
Eukaryotic translation initiation factor 2B, subunit 2 beta
Eukaryotic translation initiation factor 4 gamma, 3
Evaluation of differentially expressed genes of interest (P < 0.01) to identify over-represented families of genesa.
Gene Ontology Biological Process Categoryb
Homeostasis (Cell homeostasis, Cation homeostasis, Cell ion homeostasis, Ion homeostasis)
Amino acid and derivative metabolism (Amine metabolism, Amino acid metabolism)
Intracellular signalling cascade
Transcription from Pol II promoter (Positive regulation of transcription)
Translation (Translational initiation)
Validation of differential gene expression
Confirmation of differential gene expression 24 hours after clenbuterol administration by quantitative PCR.
Log of starting copy number
Although it has been known for decades that BA stimulate muscle hypertrophy , the data described herein are the first to provide an overview of global changes in gene expression associated with this biological model. These results also contribute to a growing body of literature describing gene expression changes associated with various models of skeletal muscle hypertrophy and atrophy [20–22]. Together, these data facilitate the opportunity to define genes across multiple models of muscle physiology, as well as to identify genes specific to the effects of BA. One unique aspect of our experiment is the investigation of gene expression at two time points relative to the administration of clenbuterol. Previous studies have focused on altered gene and/or protein expression after significant changes in muscle mass have occurred [12, 13, 23]. However, results of our study indicate that changes in gene expression are equally abundant at an early time point relative to the initiation of clenbuterol administration. These initial changes may represent important alterations of physiological pathways that culminate in altered protein turnover and/or recruitment of satellite cells to support muscle growth.
Our analysis showed that four of the 56 differentially expressed genes (10% FDR) met an even more stringent 5% FDR threshold: cyclin dependent kinase inhibitor 1A (p21) (Cdkn1a), TSC22 domain family, member 1 (Tsc22d1), growth hormone receptor (GHR), and tubulointerstitial nephritis antigen-like (Tinagl). Two of these genes, Cdkn1a and GHR, are of particular interest to this experimental model because of their well-established roles in muscle growth and development [24–28]. Cdkn1a inhibits cyclin dependent kinase 2 activity, contributing to the irreversible withdrawal from the cell cycle and terminal differentiation of myocytes . Although the presence of Cdkn1a is not required for normal development of mice , the absence of this gene product prevents normal regeneration of muscle through satellite cells . Therefore, the significant upregulation of Cdkn1a mRNA abundance following both short and long-term clenbuterol administration in our experiment may suggest an increased potential for terminal differentiation and recruitment of myogenic precursor cells in support of muscle hypertrophy. GHR is the transmembrane receptor for growth hormone (GH), a hormone essential for normal growth . Although GH effects are largely mediated through the stimulation of synthesis and secretion of insulin-like growth factor 1 (IGF1), it has been demonstrated that GH has important roles in the postnatal regulation of skeletal muscle growth that are independent of IGF1 . This work demonstrates GH signalling from GHR influences muscle hypertrophy by facilitating the fusion of myoblasts with myotubes . Our data reveal a significant down-regulation of GHR in skeletal muscle following one and ten days of clenbuterol administration. This change in GHR expression is consistent with previous reports of decreased GHR mRNA in a compensatory overload model , and fiber-type specific increased GHR mRNA in a hindlimb suspension model of muscle atrophy . Together, these reports clearly implicate GHR as a potential regulator of skeletal muscle growth and atrophy across multiple experimental models.
The EASE classification of potential genes of interest also identified overrepresented functional categories known to be important in muscle physiology and growth. Closer evaluation of genes within these categories reveals multiple pathways that potentially contribute to BA induced muscle hypertrophy. Each of the genes discussed below is included in Table 1 or 2, which include the level of significance, fold and direction of change for each gene. Additionally, all genes grouped in the categories of interest described below are identified in Supplementary Table 1.
Amino acid and derivative metabolism
The group of 16 genes categorized as being involved in 'amino acid and derivative metabolism' includes three genes critical to polyamine metabolism (S-adenosylmethionine decarboxylase 1 [Amd1; P = 0.00015], ornithine decarboxylase, structural 1 [Odc1; P = 0.00037], and spermine oxidase [Smox; P = 0.00745]). Polyamines, including spermine, spermidine, and putrescine, are polycationic compounds found in both prokaryotic and eukaryotic cells that are known to be crucial to growth and proliferation of mammalian cells [31–34]. The three genes differentially expressed in our experiment represent critical steps in the biosynthesis of polyamines. Polyamines have previously been associated with cardiac hypertrophy stimulated by BA [35–37]. For example, the Odc1 inhibitor DFMO successfully blocked cardiac hypertrophy normally associated with clenbuterol [35, 36], and overexpression of Odc1 in a transgenic model resulted in a significant increase in BA stimulated cardiac hypertrophy relative to non-transgenic mice . These experiments demonstrate a clear interaction between polyamines and beta adrenergic receptor stimulated cardiac hypertrophy. However, the relationship between polyamine metabolism and beta adrenergic receptors remains to be defined for skeletal muscle. Our data suggest that clenbuterol administration increases mRNA abundance of Odc1 while decreasing the abundance of Amd1 and Smox. This differential regulation is difficult to interpret since each of these genes participate in the synthesis of polyamines. However, increased Odc1 expression is consistent with cardiac models that demonstrate an interaction between Odc1 activity and BA. It is also of interest to note that mRNA abundance was altered for two genes that function in the transport of Arginine (solute carrier family 7 [cationic amino acid transporter, y+ system], members 2 and 5 [Slc7a2; P = 0.00051 and Slc7a5; P = 0.00012]), the precursor of ornithine. However, mRNA abundance for these transporters changed in opposite directions following clenbuterol administration. The identification of differential expression of multiple genes involved in polyamine metabolism, combined with previous data linking polyamines to BA induced cardiac hypertrophy define polyamine metabolism as a critical target for further investigation of BA induced skeletal muscle growth.
This gene ontology category was represented by 13 genes, including four eukaryotic translation initiation factors. Eukaryotic initiation factor 2 (Eif2) is involved in the first step of translation through the formation of a ternary complex between the initiator tRNA and GTP [38, 39]. This complex then binds to the 40S ribosomal subunit to initiate translation. Following start codon recognition, GTP is hydrolyzed and the resulting GDP-Eif2 complex is released. A new cycle of initiation of translation requires Eif2b to catalyze the exchange of Eif2-bound GDP for GTP, which is an important step in the regulation of translation initiation. Because Eif2b is present in low amounts, it is an important factor controlling the global rate of protein synthesis [38, 39]. Our data show that mRNA abundance of two subunits of Eif2b, Eif2b4 (P = 0.00325) and Eif2b2 (P = 0.00690), are upregulated in response to clenbuterol administration. Increased abundance of these subunits may contribute to increased activity of Eif2b and the global upregulation of protein synthesis in skeletal muscle previously associated with clenbuterol administration [40, 41]. An association between Eif2b and adrenergic receptors has also been described as Eif2b1 directly interacts with the beta-2 adrenergic receptor . Overexpression of Eif2b1 in 293 cells caused a small but significant increase in beta adrenergic receptor signalling activity. Although expression of Eif2b1 was not altered in our experiment, it is not known if the effect of Eif2b1 on beta receptor activity depends on association with other proteins, such as Eif2b2 and Eif2b4 that may be present in limiting quantities. Therefore, the observed upregulation of Eif2b2 and Eif2b4 in our study may represent a mechanism by which clenbuterol administration enhances signalling activity from beta adrenergic receptors. In addition to these Eif2b subunits, elongation initiation factor 5 (Eif5; P = 0.00679) and elongation initiation factor 4g3 (Eif4g3; P = 0.00971) also had increased mRNA abundance 24 hours after clenbuterol administration. Together, these results are consistent with a general upregulation of translational machinery occurring shortly after the administration of clenbuterol.
Transcription from Pol II promoter
A total of 18 genes were found in this category, including four genes with known functions in muscle growth and development: homeobox only domain (Hod; P = 0.00012), YY1 transcription factor (Yy1; P = 0.00004), myocyte enhancer factor 2C (Mef2c; P = 0.00368), and myogenin (Myog; P = 0.00703). The YY1 transcription factor has been associated with myogenic differentiation . In skeletal muscle, Yy1 inhibits transcription of the alpha actin gene, and down-regulation of Yy1 is necessary for alpha actin expression and myogenic differentiation to proceed . Our data show a significant increase in Yy1 mRNA abundance following 24 hours and 10 days of clenbuterol administration, potentially contributing to decreased differentiation of myogenic cells. However, alterations in Yy1 protein abundance have been shown to occur through proteolytic regulation, independent of changes in mRNA . Thus, future investigation of the regulation of Yy1 in response to BA stimulation will require careful consideration of changes at both the mRNA and protein level. Another transcription factor that was differentially expressed in our experiment was Hod, an unusual homeodomain protein that modulates cardiac growth and development [44, 45]. Hod functions by interacting with serum response factor (SRF) to inhibit activation of SRF-dependent transcription . Inactivation of Hod in transgenic mouse models confirm that an absence of Hod results in an imbalance between proliferation and differentiation of cardiomyocytes, culminating in impaired cardiac development [44, 45]. In contrast, transgenic mice that overexpress Hod develop severe cardiac hypertrophy that is mediated by inhibition of SRF-dependent transcriptional activity. Our data show a significant increase in Hod mRNA 24 hours after clenbuterol administration. Although the effect of Hod on skeletal muscle growth and differentiation has not been described, its effect on cardiac hypertrophy makes it an intriguing candidate as a potential mediator of the growth promoting effects of clenbuterol. Additionally, the regulation of Hod by BA may have important implications regarding cardiac hypertrophy commonly observed in response to these compounds . Myogenin and Mef2c are well characterized transcription factors known to be essential for myoblast differentiation [47, 48]. Although significant alteration of satellite cell differentiation and recruitment into muscle fibers in response to BA has not been described, the observed increase in mRNA of transcription factors that contribute to muscle cell differentiation suggest the recruitment and differentiation of pre-myogenic cells may be involved in the physiological response of skeletal muscle to clenbuterol administration.
Intracellular signalling cascade
Of the 34 genes of interest categorized as intracellular signalling molecules, four belong to the mitogen activated protein (MAP) kinase signalling pathway (growth arrest and DNA-damage-inducible 45 gamma [Gadd45g; P = 0.00015], mitogen activated protein kinase kinase 4 [Map2k4; P = 0.00492], mitogen activated protein kinase kinase 6 [Map2k6; P = 0.00552], and protein phosphatase 1 (formerly 2C)-like [Ppm11; P = 0.00689]). The MAP kinase signal transduction pathway is of particular interest to this skeletal muscle model because it is a key mediator of the cellular response to insulin-like growth factor 1 (IGF1), a known regulator of myogenic cell proliferation, differentiation, and protein turnover [49, 50]. In our experiment, mRNA abundance of IGF1 increased 24 hours after clenbuterol administration (P = 0.00018), but was not different from C after 10 days. Inconsistent results describing the effect of BA on IGF1 have been reported, but these likely reflect differences in the timing and sampling (local IGF1 production versus circulating levels) across experiments. Additionally, IGF2 receptor (P = 0.00065) and IGF binding protein 5 (P = 0.00051) had decreased mRNA abundance at both 24 h and 10 D time points. Although our experiment was not designed to clarify the relationship between IGF1 and BA stimulated muscle growth, it does provide data suggesting an interaction between these pathways. One potential explanation is that cross-talk between these pathways facilitates complex regulation of the recruitment of myogenic precursor cells or protein turnover within muscle fibers.
The intracellular signalling cascade category also includes three members of the ankyrin repeat and SOCS-box (Asb) containing gene family: Asb10 (P = 0.00022), Asb14 (P = 0.00245) and Asb15 (P = 0.00001). Members of this gene family have been shown to act as an E3 ligase to target specific proteins for degradation through the ubiquitin-proteasome degradation pathway [51–53]. Although specific proteins targeted for degradation through interactions with Asb10, Asb14 and Asb15 have not been identified, we have previously reported that Asb15 mRNA is down-regulated in response to BA in other animal models . Additionally, we have demonstrated that localized over-expression of Asb15 stimulates muscle hypertrophy in vivo, and that its over-expression causes a delay in differentiation and increase in protein synthesis in C2C12 cells . The current experiment provides additional support for the hypothesis that Asb family members participate in the regulation of intracellular signalling via beta adrenergic receptors by showing that multiple members of the Asb family are coordinately regulated at the mRNA level. Although much remains to be learned about the Asb gene family and its role in targeting proteins for ubiquitination, the down-regulation of three members of the Asb family supports the idea of changes in the turnover of specific, targeted proteins in response to BA.
The ubiquitin-proteasomal pathway is known to be a major contributor to protein degradation in skeletal muscle, and decreases in components of the ubiquitin-proteasome pathway have been associated with decreased protein degradation leading to muscle hypertrophy . Additionally, administration of clenbuterol has been reported to alter expression of ubiquitin-proteasome conjugates in mice following hindlimb suspension . Although the EASE analysis of our data did not identify ubiquitin related genes to be significantly over-represented, three genes involved in the ubiquitin-proteasomal pathway, including ubiquitin protein ligase E3 component n-recognin 1 [Ubr1; P = 0.00165], seven in abstentia 2 [Siah2; P = 0.00460], and proteasome subunit, beta type 1 [Psmb1; P = 0.00066], were differentially expressed. It should also be noted that the Asb gene family members previously discussed were not associated with gene ontology terms related to the ubiquitin pathway at the time of analysis. The mRNA abundance of the three ubiquitin related genes was upregulated following clenbuterol administration in our experiment. This result was unexpected because an increase in components of the ubiquitin-proteasomal pathway is associated with increased protein degradation [57, 58], while clenbuterol is known to decrease protein degradation . One potential explanation for this finding is that skeletal muscle undergoes significant fiber-type transition in response to clenbuterol administration, and the ubiquitin-proteasome system has been implicated in this process . Therefore, our results may reflect changes in skeletal muscle specifically associated with fiber-type transition, rather than a general increase in protein degradation.
Few examples of structural proteins typically associated with skeletal muscle hypertrophy were found in the list of differentially expressed genes, although myosin IB [Myo1b; P = 0.00578], actinin alpha 4 [Actn4; P = 0.00613], and myosin XVIIIa [Myo18a; P = 0.00377] had altered mRNA abundance associated with clenbuterol administration. Additionally, genes encoding proteins that interact with structural proteins were regulated. For example, actin related protein 2/3 complex subunit [Arpc3; P = 0.00005], and vinculin [Vcl; P = 0.00002] increased in mRNA abundance after administration of clenbuterol. Vinculin is an integrin-associated protein located at junctions between actin and the plasma membrane, while Arpc3 is a subunit of the Arp2/3 complex which functions in actin remodeling [60–62]. Vinculin is thought to facilitate actin organization by recruiting the Arp2/3 complex which binds to actin [60–62]. Both Arpc3 and Vcl had increased mRNA abundance following clenbuterol administration. Thus, it is reasonable to hypothesize that Arpc3 and Vcl are upregulated in the present study in order to support an increase in structural proteins associated with muscle hypertrophy.
Global evaluation of gene expression after administration of clenbuterol identified changes in gene expression and overrepresented functional categories that may regulate BA-induced muscle hypertrophy. Changes in the mRNA abundance of multiple genes associated with myogenic differentiation may indicate an important effect of BA on the proliferation, differentiation, and/or recruitment of satellite cells into muscle fibers to promote muscle hypertrophy. Additionally, increased mRNA abundance of genes involved in the initiation of translation suggests that increased levels of protein synthesis often associated with BA administration may result from a general up-regulation of translational initiators, rather than a sustained up-regulation of gene expression at the transcriptional level. Finally, numerous other genes and physiological pathways were identified that will be important targets for further investigations of the hypertrophic effect of BA on skeletal muscle.
Animal management and tissue collection
A total of 36 C57BL/6J male mice (Jackson Labs, Bar Harbor, ME), 3 to 5 wk of age, were utilized in the experiment. Mice were housed three per cage and maintained at 25°C with a 12:12-h light-dark cycle and ad libitum access to water and feed (Rodent Laboratory Chow, Ralston Purina, St. Louis, MO) throughout the experiment. All animals were handled in accordance with the protocol approved by the Purdue Animal Care and Use Committee. Four cages (three mice per cage) were randomly assigned to one of three experimental treatments (n = 12 mice per treatment): vehicle control (C), 24 hour clenbuterol treatment (24 h), or 10 day clenbuterol treatment (10 D). Clenbuterol (Sigma, St. Louis, MO) was prepared by dissolving 2.5 mg in 1 mL of ethyl alcohol, diluting to 0.25 mg/ml with 1:1 PEG200:phosphate buffered saline, and sterilized by filtration. The vehicle control treatment was prepared in the same manner, except that no clenbuterol was added. Following a seven day acclimation period, all mice were administered daily intraperitoneal injections of control or clenbuterol treatments for ten days. The C group received the control treatment all ten days, the 24 h group received the control treatment for nine days and clenbuterol on day 10, and the 10 D group received clenbuterol for all ten days. The clenbuterol treatment was administered at a dosage of 1 mg/kg of body weight, and an equivalent volume of the control treatment was injected. All mice were euthanized by CO2 asphyxiation 24 hours following the tenth injection. The gracillus muscle was collected from each animal, immediately frozen in liquid nitrogen, and stored at -80°C pending RNA extraction. The gracillus muscle was selected for analysis since it contains a mixed fiber type of both Type I and Type II fibers.
Body weight change and blood urea nitrogen level (BUN) were measured as indicators of a physiological response to the clenbuterol treatment. Individual body weights were recorded daily at the time of treatment administration (approximately 10 am), and body weight gain (WG) was calculated as the difference in body weight between days ten and one. Truncal blood was collected at the time of tissue collection, immediately following euthanasia. Serum was isolated and stored at -20°C pending analysis. Blood urea nitrogen was measured via methods of Kerscher and Ziegenhorn . WG and BUN were evaluated by analysis of variance using the GLM procedure of SAS  to test the main effect of clenbuterol treatment. Means of the 24 h and 10 D groups were compared to the C group using contrasts to determine if differences were significant.
RNA isolation and hybridization
A total of 27 mice (9 per treatment) were selected in order to maximize the difference in BUN and WG between the 24 h and 10 D treatments relative to C. Total RNA was extracted from the gracillus muscle of each mouse using the Qiagen RNeasy Mini kit following the manufacturer's recommended protocol, including an additional step of protein kinase digestion for muscle tissue (Qiagen, Valencia, CA). Quality and quantity of individual RNA samples were examined using a 2100 Bioanalyzer (Agilent, Palo Alto, CA). Three pools of RNA representing each treatment were made. Mice were randomly assigned to pools within treatment (three mice per pool) such that each pool contained an equal quantity of RNA from each mouse, for a total of greater than 50 μg of RNA at a concentration of at least 2.5 μg/μl.
Messenger RNA, cRNA synthesis and labelling reactions were performed independently for each replicate following the recommendations of the Gene Chip® Expression Analysis technical manual (Affymetrix, Santa Clara, CA). The MOE430A and B chips were hybridized to the fragmented cRNA, stained and washed according to the recommendations of the Gene Chip® Expression Analysis technical manual (Affymetrix, Santa Clara, CA) in the Purdue University Genomics Core facility. Image data were quantified using GeneChip Analysis Suite/Microarray Suite 5.0 (MAS 5.0).
If all replicates for a particular probeset were deemed 'absent', that probeset was removed from further consideration (n = 22,478) and the remaining probesets were analyzed (n = 22,605). Transcript levels were normalized to the chip median and log transformed. For each probeset, which represents the combined expression data from all relevant probe pairs on the chip, the generalized linear model Y ij = μ + B1T i + εij was fit. In each ANOVA, Y ij is a the log normalized transcript level for the ith treatment and the jth replicate, μ is the overall mean expression for the feature and T i represents the ith treatment (C, 24 h, and 10 D). An F test of the effect of treatment for each probeset was conducted as the ratio of the mean squares for treatment over the mean squares for error, and the P-value for the test of the null hypothesis t2 = t1 = t0 (i.e., mean expression not different among the three treatments) was calculated. We examined the model for conformation to the assumption of normality of the residuals testing the null hypothesis that the residuals for each gene were normally distributed using the Shapiro-Wilkes Test. All analyses were performed in SAS (SAS Institute, Cary NC ). An FDR  level of 10% was used for declaring findings significant, and a stringent rate of 5% was also examined. A nominal threshold of P < 0.01 was defined to identify genes of potential interest for the EASE analysis. If the test of the null hypothesis of difference across treatments was rejected, and we had no evidence for departure from normality of the residuals, we declared the gene differentially expressed across treatments and qualitatively assessed additional contrasts comparing the effect of the treatments (C versus 24 h and C versus 10 D).
We used the orthogonal polynomials method recently described in Yi and Xu 2006  to group the 22,605 probesets into clusters. Briefly, genes are clustered based on non-linear association between gene expression and time using orthogonal polynomials under the Gaussian mixture model. The whole data set is assumed to be a mixture of K clusters. The gene expression for individual probesets is the sum of a cluster mean (fixed effect), a normally distributed random effect that specifies deviation from the cluster mean, and a random measurement error. All genes are centered so that average expression is zero before clustering. In this way, clustering is based solely on the pattern of association of their expression with time and not the magnitude of gene expression. Bayesian Information Criterion (BIC) is used to choose the optimal number of the clusters.
Functional characterization of differentially expressed genes
Gene Ontology biological process categories were defined for genes of interest using DAVID  with the Entrez Gene ID as the primary identifier. Over-represented functional categories were identified using the EASE score , as determined by the Expression Analysis Systematic Explorer (EASE) software , with the Entrez Gene ID as the primary identifier. An EASE score <0.05 was used to define over-represented categories .
Validation of GeneChip results
Ten additional independently reared C57BL/6J male mice (Jackson Labs, Bar Harbor, ME), 3 to 5 wk of age, were utilized in the experiment. Mice were housed two to three per cage and maintained at 25°C with a 12:12-h light-dark cycle and ad libitum access to water and feed (Rodent Laboratory Chow, Ralston Purina, St. Louis, MO). All animals were handled in accordance with the protocol approved by the Purdue Animal Care and Use Committee. Mice were randomly assigned to one of two experimental treatments: vehicle control (C) or 24 hours (24 h) of clenbuterol treatment. Clenbuterol (Sigma, St. Louis, MO) was prepared as described previously. All mice were administered intraperitoneal injections of C or 24 h treatments at approximately 10 am. The clenbuterol treatment was administered at a dosage of 1 mg/kg of body weight, and an equivalent volume of the control treatment was injected. All mice were euthanized by CO2 asphyxiation 24 hours following the injection. The gracillus muscle was collected from each animal, immediately frozen in liquid nitrogen, and stored at -80°C pending RNA extraction by TRIzol (Invitrogen, Carlsbad, California).
Sequences of PCR primers used for quantitative PCR.
Forward Primer (5'-3')
Reverse Primer (5'-3')
This work was supported by the Purdue Agricultural Research Program. Martin Gonzalo assisted with the cluster analyses.
- Pellegrino MA, D'Antona G, Bortolotto S, Boschi F, Pastoris O, Bottinelli R, Polla B, Reggiani C: Clenbuterol antagonizes glucocorticoid-induced atrophy and fiber type transformation in mice. Exp Physiol. 2004, 89: 89-100. 10.1113/expphysiol.2003.002609.PubMedView ArticleGoogle Scholar
- Hinkle RT, Hodge KM, Cody DB, Sheldon RJ, Kobilka BK, Isfort RJ: Skeletal muscle hypertrophy and anti-atrophy effects of clenbuterol are mediated by the beta2-adrenergic receptor. Muscle Nerve. 2002, 25: 729-734. 10.1002/mus.10092.PubMedView ArticleGoogle Scholar
- Anderson DB, Veenhuizen EL, Schroeder AL, Jones DJ, Hancock DL: The use of phenethanolamines to reduce fat and increase leanness in meat animals. Proceedings of Symposium on Fat and Cholesterol Reduced; Advances in Applied Biotechnology Series. Edited by: Haberstroh C, Morris CE. 1991, New Orleans: Portfolio Publishing Company, 43-73.Google Scholar
- Mersmann HJ: Overview of the effects of b-adrenergic receptor agonists on animal growth including mechanisms of action. J Anim Sci. 1998, 76: 160-172.PubMedGoogle Scholar
- Moody DE, Hancock DL, Anderson DB: Phenethanolamine Repartitioning Agents. Farm Animal Metabolism and Nutrition. Edited by: d'Mello JPF. 2000, New York, NY: CABI Publishing, 65-96.View ArticleGoogle Scholar
- Dodd SL, Koesterer TJ: Clenbuterol attenuates muscle atrophy and dysfunction in hindlimb-suspended rats. Aviat Space Environ Med. 2002, 73: 635-639.PubMedGoogle Scholar
- Herrera NM, Zimmerman AN, Dykstra DD, Thompson LV: Clenbuterol in the prevention of muscle atrophy: a study of hindlimb-unweighted rats. Arch Phys Med Rehabil. 2001, 82: 930-934. 10.1053/apmr.2001.23756.PubMedView ArticleGoogle Scholar
- Higgins JA, Lasslett YV, Bardsley RG, Buttery PJ: The relation between dietary restriction or Clenbuterol (a selective beta 2 agonist) treatment on muscle growth and calpain proteinase (EC 184.108.40.206) and calpastatin activities in lambs. Br J of Nutr. 1988, 60: 645-652. 10.1079/BJN19880134.View ArticleGoogle Scholar
- Cockman MD, Jones MB, Prenger MC, Sheldon RJ: Magnetic resonance imaging of denervation-induced muscle atrophy: effects of clenbuterol in the rat. Muscle Nerve. 2001, 24: 1647-1658. 10.1002/mus.1200.PubMedView ArticleGoogle Scholar
- Fitton AR, Berry MS, McGregor AD: Preservation of denervated muscle form and function by clenbuterol in a rat model of peripheral nerve injury. J Hand Surg. 2001, 26: 335-346.View ArticleGoogle Scholar
- Elsasser TH, Rumsey TS, Kahl S, Czerwinski SM, Moseley WM, Ono Y, Solomon MB, Harris F, Fagan JM: Effects of Synovex-S and recombinant bovine growth hormone (Somavubove) on growth responses of steers: III. Muscle growth and protein responses. J Anim Sci. 1998, 76: 2346-2353.PubMedGoogle Scholar
- Grant AL, Skjaerlund DM, Gelferich WG, Bergen WG, Merkel RA: Skeletal muscle growth and expression of skeletal muscle α-actin mRNA and insulin-like growth factor I mRNA in pigs during feeding and withdrawal of ractopamine. J Anim Sci. 1993, 71: 3319-3326.PubMedGoogle Scholar
- Helferich WG, Jump DB, Anderson DB, Skjaerlund DM, Merkel RA, Bergen WG: Skeletal muscle alpha-actin synthesis is increased pretranslationally in pigs fed the phenethanolamine ractopamine. Endocrinology. 1990, 126: 3096-3100.PubMedView ArticleGoogle Scholar
- Yimlamai T, Dodd SL, Borst SE, Park S: Clenbuterol induces muscle-specific attenuation of atrophy through effects on the ubiquitin-proteasome pathway. J Appl Physiol. 2005, 99: 71-80. 10.1152/japplphysiol.00448.2004.PubMedView ArticleGoogle Scholar
- Lo HC, Ney DM: GH and IGF-I differentially increase protein synthesis in skeletal muscle and jejunum of parenterally fed rats. Am J Physiol. 1996, 271: E872-E878.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the flase discovery rate – a practical and powerful approach to multiple testing. J R Stat Soc. 1997, B 57: 289-300.Google Scholar
- Verhoeven KJF, Simonsen KL, McIntyre LM: Implementing flase discovery rate control: increasing your power. OIKAS. 2005, 108: 643-647. 10.1111/j.0030-1299.2005.13727.x.View ArticleGoogle Scholar
- Yi , Xu : Quantitative trait association microarray gene expression data analysis. Molecular Biology and Evolution.
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biology. 2003, 4: P4-10.1186/gb-2003-4-6-p4.View ArticleGoogle Scholar
- Carson JA, Nettleton D, Reecy JM: Differential gene expression in the rat soleus muscle during early work overload-induced hypertrophy. FASEB J. 2002, 16: 207-209.PubMedGoogle Scholar
- Komamura K, Shirotani-Ikejima H, Tatsumi R, Tsujita-Kuroda Y, Kitakaze M, Miyatake K, Sunagawa K, Miyata T: Differential gene expression in the rat skeletal and heart muscle in glucocorticoid-induced myopathy: analysis by microarray. Cardiovasc Drugs Ther. 2003, 17: 303-310. 10.1023/A:1027352703783.PubMedView ArticleGoogle Scholar
- Batt J, Bain J, Goncalves J, Michalski B, Plant P, Fahnestock M, Woodgett J: Differential gene expression profiling of short and long term denervated muscle. FASEB J. 2006, 20: 115-117.PubMedGoogle Scholar
- Smith SB, Garcia DK, Davis SK, Anderson DB: Elevation of a specific mRNA in longissimus muscle of steers fed ractopamine. J Anim Sci. 1989, 67: 3495-3502.PubMedGoogle Scholar
- Guo K, Wang J, Andres V, Smith RC, Walsh K: MyoD-induced expression of p21 inhibits cyclin-dependent kinase activity upon myocyte terminal differentiation. Mol Cell Biol. 1995, 15: 3823-3829.PubMedPubMed CentralView ArticleGoogle Scholar
- Deng C, Zhang P, Harper JW, Elledge SJ, Leder P: Mice lacking p21CIP1/WAF1 undergo normal development, but are defective in G1 checkpoint control. Cell. 1995, 82: 675-84. 10.1016/0092-8674(95)90039-X.PubMedView ArticleGoogle Scholar
- Hawke TJ, Meeson AP, Jiang N, Graham S, Hutcheson K, DiMaio JM, Garry DJ: p21 is essential for normal myogenic progenitor cell function in regenerating skeletal muscle. Am J Physiol Cell Physiol. 2003, 285: C1019-C1027.PubMedView ArticleGoogle Scholar
- Mukherjee A, Murray RD, Shalet SM: Impact of growth hormone status on body composition and the skeleton. Horm Res. 2004, 62: 35-41. 10.1159/000080497.PubMedView ArticleGoogle Scholar
- Sotiropoulos A, Ohanna M, Kedzia C, Menon RK, Kopchick JJ, Kelly PA, Pende M: Growth hormone promotes skeletal muscle cell fusion independent of insulin-like growth factor 1 up-regulation. Proc Natl Acad Sci. 2006, 103: 7315-7320. 10.1073/pnas.0510033103.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamaguchi A, Ikeda Y, Hirai T, Fujikawa T, Morita I: Local changes of IGF-I mRNA, GH receptor mRNA, and fiber size in rat plantaris muscle following compensatory overload. Jpn J Physiol. 2003, 53: 53-60. 10.2170/jjphysiol.53.53.PubMedView ArticleGoogle Scholar
- Casse AH, Desplanches D, Mayet-Sornay MH, Raccurt M, Jegou S, Morel G: Growth hormone receptor expression in atrophying muscle fibers of rats. Endocrinology. 2003, 144: 3692-3697. 10.1210/en.2003-0311.PubMedView ArticleGoogle Scholar
- Janne J, Alhonen L, Pietila M, Keinanen TA, Uimari A, Hyvonen MT, Pirinen E, Jarvinen A: Genetic manipulation of polyamine catabolism in rodents. J Biochem. 2006, 139: 155-160. 10.1093/jb/mvj035.PubMedView ArticleGoogle Scholar
- Nishimura K, Murozumi K, Shirahata A, Park MH, Kashiwagi K, Igarashi K: Independent roles of eIF5A and polyamines in cell proliferation. Biochem J. 2005, 385 (Pt 3): 779-785.PubMedPubMed CentralView ArticleGoogle Scholar
- Igarashi K, Kashiwagi K: Polyamine Modulon in Escherichia coli: Genes Involved in the Stimulation of Cell Growth by Polyamines. J Biochem. 2006, 139: 11-16. 10.1093/jb/mvj020.PubMedView ArticleGoogle Scholar
- Igarashi K, Kashiwagi K: Polyamines: mysterious modulators of cellular functions. Biochem Biophys Res Commun. 2000, 271: 559-564. 10.1006/bbrc.2000.2601.PubMedView ArticleGoogle Scholar
- Cubria JC, Reguera R, Balana-Fouce R, Ordonez C, Ordonez D: Polyamine-mediated heart hypertrophy induced by clenbuterol in the mouse. J Pharm Pharmacol. 1998, 50: 91-96.PubMedView ArticleGoogle Scholar
- Pegg AE: Effect of alpha-difluoromethylornithine on cardiac polyamine content and hypertrophy. J Mol Cell Cardiol. 1981, 13: 881-887. 10.1016/0022-2828(81)90287-X.PubMedView ArticleGoogle Scholar
- Shantz LM, Feith DJ, Pegg AE: Targeted overexpression of ornithine decarboxylase enhances beta-adrenergic agonist-induced cardiac hypertrophy. Biochem J. 2001, 358: 25-32. 10.1042/0264-6021:3580025.PubMedPubMed CentralView ArticleGoogle Scholar
- Proud CG: Role of mTOR signalling in the control of translation initiation and elongation by nutrients. Curr Top Microbiol Immunol. 2004, 279: 215-244.PubMedGoogle Scholar
- Hershey JW: Translational control in mammalian cells. Annu Rev Biochem. 1991, 60: 717-755. 10.1146/annurev.bi.60.070191.003441.PubMedView ArticleGoogle Scholar
- Bardocz S, Brown DS, Grant G, Pusztai A, Stewart JC, Palmer RM: Effect of the beta-adrenoceptor agonist clenbuterol and phytohaemagglutinin on growth, protein synthesis and polyamine metabolism of tissues of the rat. Br J Pharmacol. 1992, 106: 476-482.PubMedPubMed CentralView ArticleGoogle Scholar
- Rehfeldt C, Schadereit R, Weikard R, Reichel K: Effect of clenbuterol on growth, carcass and skeletal muscle characteristics in broiler chickens. Br Poult Sci. 1997, 38: 366-373.PubMedView ArticleGoogle Scholar
- Klein U, Ramirez MT, Kobilka BK, von Zastrow M: A novel interaction between adrenergic receptors and the alpha-subunit of eukaryotic initiation factor 2B. J Biol Chem. 1997, 272: 19099-19102. 10.1074/jbc.272.31.19099.PubMedView ArticleGoogle Scholar
- Walowitz JL, Bradley ME, Chen S, Lee T: Proteolytic regulation of the zinc finger transcription factor YY1, a repressor of muscle-restricted gene expression. J Biol Chem. 1998, 273: 6656-6661. 10.1074/jbc.273.12.6656.PubMedView ArticleGoogle Scholar
- Shin CH, Liu ZP, Passier R, Zhang CL, Wang DZ, Harris TM, Yamagishi H, Richardson JA, Childs G, Olson EN: Modulation of cardiac growth and development by HOP, an unusual homeodomain protein. Cell. 2002, 110: 725-735. 10.1016/S0092-8674(02)00933-9.PubMedView ArticleGoogle Scholar
- Chen F, Kook H, Milewski R, Gitler AD, Lu MM, Li J, Nazarian R, Schnepp R, Jen K, Biben C, Runke G, Mackay JP, Novotny J, Schwartz RJ, Harvey RP, Mullins MC, Epstein JA: Hop is an unusual homeobox gene that modulates cardiac development. Cell. 2002, 110: 713-723. 10.1016/S0092-8674(02)00932-7.PubMedView ArticleGoogle Scholar
- Kook H, Lepore JJ, Gitler AD, Lu MM, Wing-Man Yung W, Mackay J, Zhou R, Ferrari V, Gruber P, Epstein JA: Cardiac hypertrophy and histone deacetylase-dependent transcriptional repression mediated by the atypical homeodomain protein Hop. J Clin Invest. 2003, 112: 863-871. 10.1172/JCI200319137.PubMedPubMed CentralView ArticleGoogle Scholar
- Ridgeway AG, Wilton S, Skerjanc IS: Myocyte enhancer factor 2C and myogenin up-regulate each other's expression and induce the development of skeletal muscle in P19 cells. J Biol Chem. 2000, 275: 41-46. 10.1074/jbc.275.1.41.PubMedView ArticleGoogle Scholar
- Ma K, Chan JK, Zhu G, Wu Z: Myocyte enhancer factor 2 acetylation by p300 enhances its DNA binding activity, transcriptional activity, and myogenic differentiation. Mol Cell Biol. 2005, 25: 3575-3582. 10.1128/MCB.25.9.3575-3582.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Florini JR, Ewton DZ, Magri K: Hormones, growth factors, and myogenic differentiation. Annu Rev Physiol. 1991, 53: 201-216. 10.1146/annurev.ph.53.030191.001221.PubMedView ArticleGoogle Scholar
- Florini J, Magri K: Effects of growth factors on myogenic differentiation. Am J Physiol. 1989, 256: C701-C711.PubMedGoogle Scholar
- Chung AS, Guan YJ, Yuan ZL, Albina JE, Chin YE: Ankyrin repeat and SOCS box 3 (ASB3) mediates ubiquitination and degradation of tumor necrosis factor receptor II. Mol Cell Biol. 2005, 25: 4716-26. 10.1128/MCB.25.11.4716-4726.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Heuze ML, Guibal FC, Banks CA, Conaway JW, Conaway RC, Cayre YE, Benecke A, Lutz PG: ASB2 is an Elongin BC-interacting protein that can assemble with cullin 5 and rbx1 to reconstitute an E3 ubiquitin ligase complex. J Biol Chem. 2005, 280: 5468-5474. 10.1074/jbc.M413040200.PubMedView ArticleGoogle Scholar
- Wilcox A, Katsanakis KD, Bheda F, Phillay TS: Asb6, an adipocyte-specific ankyrin and SOCS box protein, interacts with APS to enable recruitment of Elongins B and C to the insulin receptor signalling complex. J Biol Chem. 2004, 139: 1794-1800.Google Scholar
- McDaneld TG, Hancock DL, Moody DE: Altered mRNA abundance of ASB15 and four other genes in skeletal muscle following administration of beta-adrenergic receptor agonists. Physiol Genomics. 2004, 16: 275-283.PubMedView ArticleGoogle Scholar
- McDaneld TG, Hannon K, Moody DE: Ankyrin repeat andSOCS box protein (ASB) 15 regulates protein synthesis in skeletalmuscle. Am J Physiol Regul Integr Comp Physiol. 2006,Google Scholar
- Glass DJ: Skeletal muscle hypertrophy and atrophy signalling pathways. Int J Biochem Cell Biol. 2005, 37: 1974-1984.PubMedView ArticleGoogle Scholar
- Attaix D, Ventadour S, Codran A, Bechet D, Taillandier D, Combaret L: The ubiquitin-proteasome system and skeletal muscle wasting. Essays Biochem. 2005, 41: 173-186.PubMedView ArticleGoogle Scholar
- Taillandier D, Combaret L, Pouch MN, Samuels SE, Bechet D, Attaix D: The role of ubiquitin-proteasome-dependent proteolysis in the remodelling of skeletal muscle. Proc Nutr Soc. 2004, 63: 357-361. 10.1079/PAR2004358.PubMedView ArticleGoogle Scholar
- Ordway GA, Neufer PD, Chin ER, DeMartino GN: Chronic contractile activity upregulates the proteasome system in rabbit skeletal muscle. J Appl Physiol. 2000, 88: 1134-1141.PubMedGoogle Scholar
- DeMali KA, Barlow CA, Burridge K: Recruitment of the Arp2/3 complex to vinculin: coupling membrane protrusion to matrix adhesion. J Cell Biol. 2002, 159: 881-891. 10.1083/jcb.200206043.PubMedPubMed CentralView ArticleGoogle Scholar
- Blanchoin L, Pollard TD, Hitchcock-DeGregori SE: Inhibition of the Arp2/3 complex-nucleated actin polymerization and branch formation by tropomyosin. Curr Biol. 2001, 11: 1300-1304. 10.1016/S0960-9822(01)00395-5.PubMedView ArticleGoogle Scholar
- Chopard A, Pons F, Marini JF: Vinculin and meta-vinculin in fast and slow rat skeletal muscle before and after hindlimb suspension. Pflugers Arch. 2002, 444: 627-633. 10.1007/s00424-002-0872-3.PubMedView ArticleGoogle Scholar
- Kerscher L, Ziegenhorn J: Urea. In: Methods of Enzymatic Analysis, 3rd edition, edited by Bergmeyer HU. Deerfield Beach, FL. Verlag Chemie. 1985, 8: 444-453.Google Scholar
- SAS: SAS User's Guide, (Release 6.03). 1998, Cary, NC:SAS Inst. IncGoogle Scholar
- DAVID: [http://david.abcc.ncifcrf.gov/]
- EASE: [http://david.abcc.ncifcrf.gov/ease/ease1.htm]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.