Sex-Related Differences in Gene Expression in Human Skeletal Muscle

doi:10.1371/journal.pone.0001385

PLoS ONE. 2008; 3(1): e1385.

Published online 2008 January 2. doi: 10.1371/journal.pone.0001385

PMCID: PMC2148100

Copyright Welle et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Sex-Related Differences in Gene Expression in Human Skeletal Muscle

Stephen Welle,¹^* Rabi Tawil,² and Charles A. Thornton²

¹Department of Medicine, University of Rochester, Rochester, New York, United States of America

²Department of Neurology, University of Rochester, Rochester, New York, United States of America

Cecile Fairhead, Academic Editor

Pasteur Institute, France

* To whom correspondence should be addressed. E-mail: stephen_welle/at/urmc.rochester.edu

Conceived and designed the experiments: SW. Performed the experiments: SW RT CT. Analyzed the data: SW. Wrote the paper: SW.

Received November 27, 2007; Accepted December 9, 2007.

This article has been cited by other articles in PMC.

Abstract

There is sexual dimorphism of skeletal muscle, the most obvious feature being the larger muscle mass of men. The molecular basis for this difference has not been clearly defined. To identify genes that might contribute to the relatively greater muscularity of men, we compared skeletal muscle gene expression profiles of 15 normal men and 15 normal women by using comprehensive oligonucleotide microarrays. Although there were sex-related differences in expression of several hundred genes, very few of the differentially expressed genes have functions that are obvious candidates for explaining the larger muscle mass of men. The men tended to have higher expression of genes encoding mitochondrial proteins, ribosomal proteins, and a few translation initiation factors. The women had >2-fold greater expression than the men (P<0.0001) of two genes that encode proteins in growth factor pathways known to be important in regulating muscle mass: growth factor receptor-bound 10 (GRB10) and activin A receptor IIB (ACVR2B). GRB10 encodes a protein that inhibits insulin-like growth factor-1 (IGF-1) signaling. ACVR2B encodes a myostatin receptor. Quantitative RT-PCR confirmed higher expression of GRB10 and ACVR2B genes in these women. In an independent microarray study of 10 men and 9 women with facioscapulohumeral dystrophy, women had higher expression of GRB10 (2.7-fold, P<0.001) and ACVR2B (1.7-fold, P<0.03). If these sex-related differences in mRNA expression lead to reduced IGF-1 activity and increased myostatin activity, they could contribute to the sex difference in muscle size.

Introduction

There is sexual dimorphism of skeletal muscle in overall mass, size of individual fibers, activities of several metabolic enzymes, lipid content and oxidation, relative expression of different myosin isoforms, fatigability, and expression of a number of genes [1]–[9]. Although all types of muscle fibers are larger in men, the sex difference is especially pronounced in type 2 fibers so that there is a greater ratio of type 2 fiber mass to type 1 fiber mass in men [1], [5], [10]. The sex difference in muscle mass is presumed to be mediated by higher testosterone levels in men, because of the well known anabolic effect of testosterone [11]–[13] and because estrogens and progestins do not reduce muscle mass [14]–[17]. The relative enlargement of muscle in males develops after the pubertal increase in testosterone production. After that, men require testosterone to maintain a normal muscle mass [11]–[13], [18], [19]. Testosterone, like all steroid hormones, exerts its effects by influencing gene expression. It has not been established which genes are responsible for its anabolic effects. While some effects of testosterone on gene expression might be limited to the period of rapid muscle growth after puberty, there must be some permanent effects to maintain the larger muscle mass in men.

There have been few comparisons of broad gene expression profiles in men and women. Roth et al. reported differences between men and women in the muscle expression of ~20% of ~1,000 transcripts that yielded reliable signals on cDNA arrays [4]. Yoshioka et al. [20] used serial analysis of gene expression (SAGE) to compare muscle gene expression in male and female mice, but SAGE is better suited for studying highly-expressed genes, such as those encoding metabolic and contractile proteins, than for studying the majority of transcripts that are expressed at lower levels. We have used comprehensive oligonucleotide arrays to study the effect of aging on expression profiles of both men and women [21], [22], but have not previously reported the sex differences. The purpose of this report is to summarize the key features of the sex-related differences in gene expression in these subjects.

Results

We obtained muscle samples from the vastus lateralis from normal adult subjects, including 15 men and 15 women 20–75 years old. To limit variability of activity and diet prior to the biopsies, subjects were admitted to the University of Rochester General Clinical Research Center for 3 days, where they were provided a standard weight-maintaining diet and were instructed not to perform any activity more strenuous than walking. Table 1 shows their mean body composition (by dual-energy X-ray absorptiometry) and isometric knee extension strength. The men had more lean tissue mass overall (mean 41%) and in the legs (mean 48%) than the women. There did not appear to be a sex difference in physical fitness—isometric knee extension strength and maximal oxygen consumption were similar in the men and women when expressed per kg lean body mass or lean tissue mass of the legs.

Table 1

Mean±SEM body composition, muscle performance, and MYH gene expression patterns in 15 men and 15 women with normal muscles who donated tissue for gene expression study.

Table 1 also shows the expression of genes encoding types 1, 2a, and 2× myosin heavy chains (MYH). Expression of these genes was determined by RT-PCR because these highly expressed mRNAs saturated the microarray probes. Based on the fact that women have a higher ratio of type 1 fiber mass to type 2 fiber mass, the higher ratio of type 1 MYH to types 2a and 2× MYH was expected. (To avoid confusion, we did not use the official gene symbols in Table 1; types 1, 2a, and 2× fibers express the MYH7, MYH2, and MYH1 genes, respectively.)

Affymetrix U133A and U133B high-density oligonucleotide arrays, containing ~45,000 probe sets targeting all known human mRNAs and expressed sequence tags at the time of UniGene build 133, were used to examine gene expression. After elimination of probe sets that did not generate signals significantly above cross-hybridization (“mismatch”) control oligos, as well as those with signals too strong to be accurately quantified, there were ~17,000 probe sets included in the analysis. The sex-related differences for these probe sets are available in a supporting information file (Excel File S1). The CEL files and GCRMA-normalized data for every probe set are available from the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) via accession number GSE9676.

The Significance Analysis of Microarrays (SAM) program [23], which evaluates the distribution of the t statistic with numerous random permutations of group assignments (male vs. female in this study), indicated that there were more than 3,000 sex-related differences when limiting the false discovery rate to 5%. This method did not account for age as a source of variance, so we used analysis of variance (ANOVA) as the primary method to rank genes according to the likelihood that there is a sex difference in expression after accounting for overall age effects (which have been reported previously). According to ANOVA there were more than 2,400 probe sets with a sex difference at P<0.01, more than 800 at P<0.001, and more than 300 at P<0.0001 (Figure 1). These nominal P values do not account for the fact that 17,000 comparisons were made simultaneously, but by comparison with the number of differences suggested by SAM it appears that the proportion of false discoveries at nominal P<0.01 should be very low.

Figure 1

Volcano plot illustrating statistical significance of sex-related differences (shown as negative log of the ANOVA P level, so that the most significant differences are on top) in relation to magnitude of mean differences (shown as log of the female/male (more ...)

As expected, the largest sex differences in signal intensity were observed with Y-chromosome genes and XIST, which mediates inactivation of one copy of the X chromosome in female cells [24] (Figure 1).

Table 2 lists the sex differences that stood out in the present study both in magnitude (≥2-fold) and consistency (P<0.0001 by ANOVA). Presentation of only these genes does not necessarily mean that these are the most important genes. This short list is presented to illustrate the fact that genes with various functions are differentially expressed, and that most of these genes have not been studied by researchers interested in the regulation of muscle mass. An expanded list (less stringent P and fold-difference criteria) would show the same thing.

Table 2

Genes with ≥2-fold differential expression in muscles of men and women at P<0.0001, excluding XIST and Y-chromosome genes.

Two of the genes expressed at higher levels in women, GRB10 and ACVR2B, encode proteins that are in signaling pathways of growth factors known to regulate muscle mass. The sex difference in GRB10 expression was especially convincing, demonstrating a 2.3- to 2.5-fold difference with three probe sets as highlighted in Figure 1 (ACVR2B has only one probe set). The increased insulin receptor gene expression in women could be related to the higher GRB10 expression (see Discussion).

Because of the potential significance of GRB10 and ACVR2B genes as regulators of muscle mass, we used quantitative RT-PCR as an independent method to examine their expression. The HPRT1 transcript was selected as the reference gene on an empirical basis—HPRT1 expression in these men (1863±111 arbitrary units, mean±SEM) and women (1881±138 arbitrary units) was similar according to the microarray data (this choice does not necessarily mean that HPRT1 expression always is constant in muscle). According to this method, GRB10 expression was 2.6-fold greater in the women (P<0.0001), and ACVR2B expression was 1.6-fold greater (P<0.01). GRB10/HPRT1 ratios by microarray in individual subjects correlated closely with these ratios determined by RT-PCR (r

0.91, P<0.0001). ACVR2B/HPRT1 ratios by microarray and PCR also correlated well (r

0.69, P<0.01).

We also examined our database of expression profiles of 10 men and 9 women with facioscapulohumeral dystrophy [25] for a sex difference in expression of GRB10 and ACVR2B. The women had higher expression of GRB10 (2.7-fold, P<0.001) and ACVR2B (1.7-fold, P<0.03).

Two bioinformatics resources, Gene Set Enrichment Analysis (GSEA) [26] and Expression Analysis Systematic Explorer (EASE) [27], were used to relate differentially expressed genes to functional categories. Rather than focusing on the largest differences among individual genes, GSEA integrates similar trends in expression levels among multiple genes related to one another (functionally, by co-regulation in previous experiments, by chromosomal location, or by similarity of transcription factor binding domains). According to the GSEA approach, genes expressed at higher levels in women were enriched (false discovery rate <5%) in a Chromosome 6q24 gene set. We have no explanation for this cluster. Genes expressed at higher levels in men were enriched in a “testis-expressed” set, with numerous Y-chromosome genes, and the “UVB_NHEK2_UP” gene set. The top member of the latter set in the current study was a Y-chromosome gene, and it is likely that the similarity of genes sets would not have been statistically significant without inclusion of this gene.

EASE sorts genes according to Gene Ontology categories. This program compares a list of “significant” genes to the entire list of genes included on the array to determine if certain categories are represented more than expected by chance. Use of very strict criteria to define differential expression limits the utility of this approach, so we selected genes with sex differences at nominal P<0.01 with no fold-difference criterion. With genes expressed at higher levels in women, there were no categories enriched at P<0.05 (corrected for multiple tests by the Bonferroni method). With genes expressed at higher levels in men, there were several enriched categories (Figure 2). Genes encoding proteins involved in protein synthesis accounted for enrichment of the “macromolecule biosynthesis” category. Numerous ribosomal protein genes were the main reason for enrichment of the “RNA binding” and “structural constituent of ribosome” categories, and also were a major part of the “macromolecule biosynthesis” category. Several translation initiation factor genes also were in the biosynthesis category. The “H⁺ transporter” list was mostly genes encoding proteins of the electron transport chain and mitochondrial ATP synthase. These genes also contributed to enrichment of the “oxidative phosphorylation” and “mitochondrion” categories. Genes encoding mitochondrial ribosomal proteins contributed to enrichment of both the “biosynthesis” and “mitochondrion” categories.

Figure 2

Summary of functional categories of genes expressed at higher levels in men than in women (nominal P<0.01).

Discussion

It is well accepted that IGF-1 has an anabolic effect on muscle [28]. Although IGF1 gene expression in muscle is sensitive to testosterone depletion or administration in men [18], [29], no sex difference in plasma IGF-1 levels or muscle IGF1 gene expression has been observed [10]. The microarrays used in the present study also did not detect any sex difference in IGF1 gene expression in muscle. However, IGF-1 signaling might be reduced in women because of their elevated expression of GRB10, which encodes a protein that interferes with IGF-1 signal transduction [30], [31]. GRB10 also interferes with insulin signaling [31], [32], in part by mediating degradation of the insulin receptor [33], but the 3-fold increase in insulin receptor gene expression in women (Table 2) could serve to minimize any negative effect on insulin sensitivity that increased GRB10 expression might have. In contrast, there is not a sex difference in expression of the IGF-1 receptor according to the present microarray study and a previous study in which IGFR gene expression was determined by RT-PCR analysis [10]. Thus, the greater GRB10 expression in women might restrain IGF-1 signaling and contribute to the sex difference in muscularity. The plausibility of this hypothesis is supported by recent studies demonstrating that GRB10 knockout in mice was associated with increased muscularity [31], [34].

Myostatin has a major role in determining muscle size. Constitutive myostatin knockout in mice leads to a hypermuscular phenotype caused by an increased number of fibers and fiber enlargement [35]. Myostatin appears to function similarly in humans [36]. Post-developmental inhibition of myostatin activity causes muscle hypertrophy in mice without increasing the number of muscle fibers [37]–[39]. Although there is no evidence for differential expression of myostatin mRNA in men and women either from the present study or previous analyses by RT-PCR [40], [41], this does not necessarily mean that the myostatin pathway has the same level of activity in men and women. For example, myostatin protein levels in muscle are higher in female than male mice even though mRNA levels are similar in males and females [42]. Moreover, myostatin activity could differ in men and women if there is a difference in any of the proteins that process myostatin, bind to myostatin, or mediate myostatin signal transduction. Here, we report that a myostatin receptor gene, ACVR2B, is expressed at a higher level in women. In mice, knockout of either activin A receptor IIA (ACVR2A) or ACVR2B causes muscle hypertrophy [38]. Hemizygous knockout of ACVR2B accentuates muscle growth in ACVR2A-knockout mice, although it has no effect in mice with normal ACVR2A expression [38]. The microarrays indicated that the ACVR2A gene is expressed in human muscle, with no significant difference between men and women. Thus, understanding the functional significance of higher ACVR2B expression in women will require elucidation of the relative importance of ACVR2A and ACVR2B in myostatin signal transduction in humans, and whether ACVR2B expression is a limiting factor. It was reported recently that women with different ACVR2B haplotypes had different levels of quadriceps strength, although this was not observed in men, and there was no difference in leg muscle mass as assessed by DEXA [43]. It is not known whether the polymorphisms in the ACVR2B genes in the different haplotypes affect ACVR2B expression in muscle.

Because testosterone is assumed to be responsible for the greater muscle mass in men, an obvious issue is whether it regulates expression of GRB10 or ACVR2B genes. Consensus androgen responsive elements [44], [45] are not near (within 5 kb) the transcription start sites of these genes. It must be emphasized that many, if not most, functional androgen responsive elements differ in several bases from the consensus motif [46]. Hence it is not currently possible to prove or rule out by computer algorithms whether testosterone regulates any particular gene. Even if there are no androgen responsive DNA elements regulating these genes directly, testosterone could influence their expression indirectly by its effects on expression of other genes. Studies of men undergoing pharmacologic inhibition of testosterone production would be the best approach to determine whether these sex differences in gene expression are mediated by testosterone.

One of the genes listed in Table 2, FOXO3, encodes a forkhead box transcription factor that can influence muscle catabolism. In its active non-phosphorylated state, this transcription factor stimulates expression of the ubiquitin ligase atrogin-1 (also known as MAFbx or F-box protein 32) [47]. Atrogin-1 is upregulated under conditions associated with muscle atrophy [48]. Atrogin-mediated ubiquitination of proteins increases the rate of proteasomal proteolysis. Although FOXO3 expression was ~2-fold higher in women than in men, there was no evidence from the microarray data for increased atrogin-1 expression in women. Moreover, there is no evidence for increased fractional muscle protein turnover in women [49]–[51], suggesting that increased muscle proteolysis is not responsible for the smaller muscle mass of women.

Most investigators using microarrays focus on the largest and most consistent differences in gene expression. However, more subtle differences in gene expression might have significant consequences, particularly if several genes involved in the same pathway are involved. We used the GSEA and EASE methods to search for sex-related differences in pathways and functional categories of genes that might not be obvious based on the largest fold differences or lowest P levels. Given the greater ratio of type 1 muscle fiber mass to type 2 mass in women, one might expect a tendency for greater expression in women of genes encoding mitochondrial proteins. However, the opposite trend was observed. There also was a tendency for the men to have higher expression of ribosomal proteins (both cytosolic and mitochondrial) and translation initiation factors. Higher expression of these genes might support an increased rate of protein synthesis in men. As mentioned above, men and women have similar fractional muscle protein breakdown rates. Therefore, men must have a higher rate of protein synthesis per muscle fiber or else they would not have a greater protein mass per fiber. (Even though synthesis per fiber is greater in men, the fractional rate is similar in men and women because protein mass per muscle fiber is greater in men.)

The GSEA program did not find enrichment in any of the Gene Ontology categories that were found by the EASE program to be more highly expressed in men. This discrepancy can be explained by the difference in the approaches of the two methods. EASE does not impose a penalty if differences opposite to the overall trend are present, whereas GSEA does because it takes into consideration all differences in both directions regardless of the statistical significance for individual genes. For example, there were 9 translation initiation factor genes expressed at a higher level in men (at nominal P<0.01), but there were 4 expressed at a higher level in women. For some other categories shown in Figure 2, the ratio of significant effects in men versus women were 22/4 (mitochondrial ribosomal proteins), 17/3 (ribosomal proteins), and 7/2 (lysosomal H⁺ transporters).

The functional significance of sex-related differences in mRNA levels depends on whether they cause differences in protein levels. Unfortunately, we do not have any more muscle tissue from these subjects to examine protein levels. We encourage others who are interested in the sexual dimorphism of muscle to examine mRNA and protein expression of GRB10, ACVR2B, or other genes that were differentially expressed in men and women in the present study.

Methods

Subjects

Muscle was obtained by needle biopsy, with local anesthesia, from 30 subjects with normal muscles [7 young men (21–27 yr old), 8 older men (65–75 yr old), 7 young women (20–29 yr old), 8 older women (65–71 yr old)]. Quadriceps strength, maximal oxygen consumption, and body composition measures were done as described elsewhere [21], [22]. All subjects were fully informed of risks and gave written consent for participation. The study was approved by the University of Rochester Research Subjects Review Board.

Expression profiles

Affymetrix U133A and U133B GeneChip arrays were used to examine gene expression. Details of RNA extraction and microarray methods have been described elsewhere [21], [22]. The gene expression database was examined for differential expression in men and women after re-calculating expression scores with the GCRMA method [52] as implemented by ArrayAssist Lite software (Stratagene, Version 3.3). This method was chosen after comparing it the standard Affymetrix method, a modified version of the Affymetrix comparison analysis algorithm that was used for our analyses of age-related differences [53], and RMA [54], because it had the best performance in accounting for non-specific hybridization (expression scores in women close to zero for Y chromosome genes; in men close to zero for XIST gene) and had the highest correlation with quantitative RT-PCR for mean age-related differences across 8 different transcripts (r

0.83). Probe sets with high (P>0.10) “detection P” values (Affymetrix algorithm) for more than half of the arrays from the group with higher mean expression were excluded from further analysis, which reduces the number of false positives [53], [55]. Probe sets that were saturated by very abundant transcripts also were excluded from the analysis. ANOVA for the probe sets that passed these filters was done with a Visual Basic script to calculate F ratios, which were exported to a Microsoft Excel spreadsheet for determination of nominal P values. Excel also was used for computation of mean differences between men and women. The SAM program (Version 2.21), which runs within Excel, was provided by Stanford University [23]. The GSEA program (Version 2.0) was provided by the Broad Institute [26]. The EASE program (Version 2.0) was provided by the National Institutes of Health [27].

Quantitative RT-PCR

Expression of the MYH genes was determined with a competitive PCR method [56]. This was done prospectively, so RNA was available from all subjects for these assays. GRB10 and ACVR2B transcripts were examined by quantitative real time PCR. These assays were done after the microarray analyses were completed, at which time RNA samples were available from 8 men and 7 women. HPRT1 mRNA served as the reference transcript. RNA was treated with DNase I, then reverse transcribed by MMLV reverse transcriptase with oligo-dT as the primer. Real-time PCR quantification of cDNAs was done in triplicate reactions with the Applied Biosystems (ABI) Prism 7900HT Sequence Detection System. Primer and probe sets for TaqMan assays were purchased from ABI (for GRB10, assay number Hs00959287_m1; for ACVR2B, Hs00609603_m1; for HPRT1, catalog number 4333768F ). The unit of measure in these assays is the difference in cycle threshold (ΔCt) between the cDNA of interest and the reference cDNA, with an increase in ΔCt of one cycle reflecting a decrease of ~2-fold in the GRB10/HPRT1 or ACVR2B/HPRT1 ratio.

Supporting Information

Excel File S1

Sex differences for all transcripts detected in muscle. Affymetrix probe set IDs, female/male expression ratios, P values (from ANOVA), and gene names and IDs for all transcripts detected in skeletal muscle.

(3.07 MB XLS)

Click here for additional data file.^{(2.9M, xls)}

Acknowledgments

We thank Nancy Needler, Don Henderson, Andrew Cardillo, Kirti Bhatt, Bharati Shah, and the staff of the University of Rochester General Clinical Research Center for their assistance.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: The project was supported by grants from the National Institutes of Health (AG17891, AG18254, AR/NS48143, RR00044). The funding agency had no other role in the conduct of this study or preparation of this paper.

References

1. Simoneau JA, Bouchard C. Human variation in skeletal muscle fiber-type proportion and enzyme activities. Am J Physiol. 1989;257:E567–E572. [PubMed]

2. Coggan AR, Spina RJ, King DS, Rogers MA, Brown M, et al. Histochemical and enzymatic comparison of the gastrocnemius muscle of young and elderly men and women. J Gerontol A Biol Sci Med Sci. 1992;47:B71–B76.

3. Forsberg AM, Nillson E, Wernerman J, Bergstrom J, Hultman E. Muscle composition in relation to age and sex. Clin Sci. 1991;81:249–256. [PubMed]

4. Roth SM, Ferrell RE, Peters DG, Metter EJ, Hurley BF, et al. Influence of age, sex, and strength training on human muscle gene expression determined by microarray. Physiol Genomics. 2002;10:181–190. [PMC free article] [PubMed]

5. Staron RS, Hagerman FC, Hikida RS, Murray TF, Hostler DP, et al. Fiber type composition of the vastus lateralis muscle of young men and women. J Histochem Cytochem. 2000;48:623–629. [PubMed]

6. Steffensen CH, Roepstorff C, Madsen M, Kiens B. Myocellular triacylglycerol breakdown in females but not in males during exercise. Am J Physiol Endocrinol Metab. 2002;282:E634–E642. [PubMed]

7. Roepstorff C, Steffensen CH, Madsen M, Stallknecht B, Kanstrup IL, et al. Gender differences in substrate utilization during submaximal exercise in endurance-trained subjects. Am J Physiol Endocrinol Metab. 2002;282:E435–E447. [PubMed]

8. Hunter SK, Enoka RM. Sex differences in the fatigability of arm muscles depends on absolute force during isometric contractions. J Appl Physiol. 2001;91:2686–2694. [PubMed]

9. Kent-Braun JA, Ng AV, Doyle JW, Towse TF. Human skeletal muscle responses vary with age and gender during fatigue due to incremental isometric exercise. J Appl Physiol. 2002;93:1813–1823. [PubMed]

10. Bamman MM, Hill VJ, Adams GR, Haddad F, Wetzstein CJ, et al. Gender differences in resistance-training-induced myofiber hypertrophy among older adults. J Gerontol A Biol Sci Med Sci. 2003;58:108–116. [PubMed]

11. Bhasin S, Storer TW, Berman N, Yarasheski KE, Clevenger B, et al. Testosterone replacement increases fat-free mass and muscle size in hypogonadal men. J Clin Endocrinol Metab. 1997;82:407–413. [PubMed]

12. Brodsky IG, Balagopal P, Nair KS. Effects of testosterone replacement on muscle mass and muscle protein synthesis in hypogonadal men--a Clinincal Research Center study. J Clin Endocrinol Metab. 1996;81:3469–3475. [PubMed]

13. Snyder PJ, Peachey H, Berlin JA, Hannoush P, Haddad G, et al. Effects of testosterone replacement in hypogonadal men. J Clin Endocrinol Metab. 2000;85:2670–2677. [PubMed]

14. Jensen J, Christiansen C, Rodbro P. Oestrogen-progestogen replacement therapy changes body composition in early post-menopausal women. Maturitas. 1986;8:209–216. [PubMed]

15. Hassager C, Christiansen C. Estrogen/gestagen therapy changes soft tissue body composition in postmenopausal women. Metabolism. 1989;38:662–665. [PubMed]

16. Sipila S, Taaffe DR, Cheng S, Puolakka J, Toivanen J, et al. Effects of hormone replacement therapy and high-impact physical exercise on skeletal muscle in post-menopausal women: a randomized placebo-controlled study. Clin Sci (Lond) 2001;101:147–157. [PubMed]

17. Taaffe DR, Newman AB, Haggerty CL, Colbert LH, de Rekeneire N, et al. Estrogen replacement, muscle composition, and physical function: The Health ABC Study. Med Sci Sports Exerc. 2005;37:1741–1747. [PubMed]

18. Mauras N, Hayes V, Welch S, Rini A, Helgeson K, et al. Testosterone deficiency in young men: marked alterations in whole body protein kinetics, strength and adiposity. J Clin Endocrinol Metab. 1998;83:1886–1892. [PubMed]

19. Smith MR, Finkelstein JS, McGovern FJ, Zietman AL, Fallon MA, et al. Changes in body composition during androgen deprivation therapy for prostate cancer. J Clin Endocrinol Metab. 2002;87:599–603. [PubMed]

20. Yoshioka M, Boivin A, Bolduc C, St Amand J. Gender difference of androgen actions on skeletal muscle transcriptome. J Mol Endocrinol. 2007;39:119–133. [PubMed]

21. Welle S, Brooks AI, Delehanty JM, Needler N, Thornton CA. Gene expression profile of aging in human muscle. Physiol Genomics. 2003;14:149–159. [PubMed]

22. Welle S, Brooks AI, Delehanty JM, Needler N, Bhatt K, et al. Skeletal muscle gene expression profiles in 20-29 year old and 65-71 year old women. Exp Gerontol. 2004;39:369–377. [PubMed]

23. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001;98:5116–5121. [PubMed]

24. Brown CJ, Greally JM. A stain upon the silence: genes escaping X inactivation. Trends Genet. 2003;19:432–438. [PubMed]

25. Osborne RJ, Welle S, Venance SL, Thornton CA, Tawil R. Expression profile of FSHD supports a link between retinal vasculopathy and muscular dystrophy. Neurology. 2007;68:569–577. [PubMed]

26. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–15550. [PubMed]

27. Hosack DA, Dennis G, Jr, Sherman BT, Lane HC, Lempicki RA. Identifying biological themes within lists of genes with EASE. Genome Biol. 2003;4:R70. [PMC free article] [PubMed]

28. Glass DJ. Skeletal muscle hypertrophy and atrophy signaling pathways. Int J Biochem Cell Biol. 2005;37:1974–1984. [PubMed]

29. Urban RJ, Bodenburg YH, Gilkison C, Foxworth J, Coggan AR, et al. Testosterone administration to elderly men increases skeletal muscle strength and protein synthesis. Am J Physiol. 1995;269:E820–E826. [PubMed]

30. Dufresne AM, Smith RJ. The adapter protein GRB10 is an endogenous negative regulator of insulin-like growth factor signaling. Endocrinology. 2005;146:4399–4409. [PubMed]

31. Smith FM, Holt LJ, Garfield AS, Charalambous M, Koumanov F, et al. Mice with a disruption of the imprinted Grb10 gene exhibit altered body composition, glucose homeostasis, and insulin signaling during postnatal life. Mol Cell Biol. 2007;27:5871–5886. [PMC free article] [PubMed]

32. Langlais P, Dong LQ, Ramos FJ, Hu D, Li Y, et al. Negative regulation of insulin-stimulated mitogen-activated protein kinase signaling by Grb10. Mol Endocrinol. 2004;18:350–358. [PubMed]

33. Ramos FJ, Langlais PR, Hu D, Dong LQ, Liu F. Grb10 mediates insulin-stimulated degradation of the insulin receptor: a mechanism of negative regulation. Am J Physiol Endocrinol Metab. 2006;290:E1262–E1266. [PubMed]

34. Wang L, Balas B, Christ-Roberts CY, Kim RY, Ramos FJ, et al. Peripheral disruption of the grb10 gene enhances insulin signaling and sensitivity in vivo. Mol Cell Biol. 2007;27:6497–6505. [PMC free article] [PubMed]

35. McPherron AC, Lawler AM, Lee S-J. Regulation of skeletal muscle mass in mice by a new TGF-β superfamily member. Nature. 1997;387:83–90. [PubMed]

36. Schuelke M, Wagner KR, Stolz LE, Hubner C, Riebel T, et al. Myostatin mutation associated with gross muscle hypertrophy in a child. N Engl J Med. 2004;350:2682–2688. [PubMed]

37. Whittemore LA, Song K, Li X, Aghajanian J, Davies M, et al. Inhibition of myostatin in adult mice increases skeletal muscle mass and strength. Biochem Biophys Res Commun. 2003;300:965–971. [PubMed]

38. Lee SJ, Reed LA, Davies MV, Girgenrath S, Goad ME, et al. Regulation of muscle growth by multiple ligands signaling through activin type II receptors. Proc Natl Acad Sci U S A. 2005;102:18117–18122. [PubMed]

39. Welle S, Bhatt K, Pinkert CA, Tawil R, Thornton CA. Muscle growth after postdevelopmental myostatin gene knockout. Am J Physiol Endocrinol Metab. 2007;292:E985–E991. [PubMed]

40. Kim JS, Cross JM, Bamman MM. Impact of resistance loading on myostatin expression and cell cycle regulation in young and older men and women. Am J Physiol Endocrinol Metab. 2005;288:E1110–E1119. [PubMed]

41. Roth SM, Martel GF, Ferrell RE, Metter EJ, Hurley BF, et al. Myostatin gene expression is reduced in humans with heavy-resistance strength training: a brief communication. Exp Biol Med. 2003;228:706–709.

42. McMahon CD, Popovic L, Jeanplong F, Oldham JM, Kirk SP, et al. Sexual dimorphism is associated with decreased expression of processed myostatin in males. Am J Physiol Endocrinol Metab. 2003;284:E377–E381. [PubMed]

43. Walsh S, Metter EJ, Ferrucci L, Roth SM. Activin-type II receptor B (ACVR2B) and follistatin haplotype associations with muscle mass and strength in humans. J Appl Physiol. 2007;102:2142–2148. [PMC free article] [PubMed]

44. Nelson CC, Hendy SC, Shukin RJ, Cheng H, Bruchovsky N, et al. Determinants of DNA sequence specificity of the androgen, progesterone, and glucocorticoid receptors: evidence for differential steroid receptor response elements. Mol Endocrinol. 1999;13:2090–2107. [PubMed]

45. Verrijdt G, Haelens A, Claessens F. Selective DNA recognition by the androgen receptor as a mechanism for hormone-specific regulation of gene expression. Mol Genet Metab. 2003;78:175–185. [PubMed]

46. Reid KJ, Hendy SC, Saito J, Sorensen P, Nelson CC. Two classes of androgen receptor elements mediate cooperativity through allosteric interactions. J Biol Chem. 2001;276:2943–2952. [PubMed]

47. Sandri M, Sandri C, Gilbert A, Skurk C, Calabria E, et al. Foxo transcription factors induce the atrophy-related ubiquitin ligase atrogin-1 and cause skeletal muscle atrophy. Cell. 2004;117:399–412. [PubMed]

48. Lecker SH, Jagoe RT, Gilbert A, Gomes M, Baracos V, et al. Multiple types of skeletal muscle atrophy involve a common program of changes in gene expression. FASEB J. 2004;18:39–51. [PubMed]

49. Welle S, Thornton C, Statt M, McHenry B. Postprandial myofibrillar and whole-body protein synthesis in young and old human subjects. Am J Physiol. 1994;267:E599–E604. [PubMed]

50. Welle S, Thornton C, Statt M. Myofibrillar protein synthesis in young and old human subjects after three months of resistance training. Am J Physiol. 1995;268:E422–E427. [PubMed]

51. Jahn LA, Barrett EJ, Genco ML, Wei L, Spraggins TA, et al. Tissue composition affects measures of postabsorptive human skeletal muscle metabolism: comparison across genders. J Clin Endocrinol Metab. 1999;84:1007–1010. [PubMed]

52. Wu Z, Irizarry RA. Preprocessing of oligonucleotide array data. Nat Biotechnol. 2004;22:656–658. [PubMed]

53. Welle S, Brooks AI, Thornton CA. Computational method for reducing variance with Affymetrix microarrays. BMC Bioinformatics. 2002;3:23. [PMC free article] [PubMed]

54. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–264. [PubMed]

55. McClintick JN, Edenberg HJ. Effects of filtering by Present call on analysis of microarray experiments. BMC Bioinformatics. 2006;7:49. [PMC free article] [PubMed]

56. Welle S, Bhatt K, Thornton CA. Stimulation of myofibrillar synthesis by exercise is mediated by more efficient translation of mRNA. J Appl Physiol. 1999;86:1220–1225. [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science