Abstract

Spaceflight results in a number of adaptations to skeletal muscle, including atrophy and shifts toward faster muscle fiber types. To identify changes in gene expression that may underlie these adaptations, we used both microarray expression analysis and real-time polymerase chain reaction to quantify shifts in mRNA levels in the gastrocnemius from mice flown on the 11-day, 19-h STS-108 shuttle flight and from normal gravity controls. Spaceflight data also were compared with the ground-based unloading model of hindlimb suspension, with one group of pure suspension and one of suspension followed by 3.5 h of reloading to mimic the time between landing and euthanization of the spaceflight mice. Analysis of microarray data revealed that 272 mRNAs were significantly altered by spaceflight, the majority of which displayed similar responses to hindlimb suspension, whereas reloading tended to counteract these responses. Several mRNAs altered by spaceflight were associated with muscle growth, including the phosphatidylinositol 3-kinase regulatory subunit p85α, insulin response substrate-1, the forkhead box O1 transcription factor, and MAFbx/atrogin1. Moreover, myostatin mRNA expression tended to increase, whereas mRNA levels of the myostatin inhibitor FSTL3 tended to decrease, in response to spaceflight. In addition, mRNA levels of the slow oxidative fiber-associated transcriptional coactivator peroxisome proliferator-associated receptor (PPAR)-γ coactivator-1α and the transcription factor PPAR-α were significantly decreased in spaceflight gastrocnemius. Finally, spaceflight resulted in a significant decrease in levels of the microRNA miR-206. Together these data demonstrate that spaceflight induces significant changes in mRNA expression of genes associated with muscle growth and fiber type.

spaceflight has many consequences on the musculoskeletal system. Specifically, the unloading due to microgravity results in decreases in muscle mass and force production capability as well as shifts in muscle fiber type toward a faster contracting and less endurable phenotype (10). Qualitatively similar decreases in muscle mass and shifts in fiber type occur in both rats (35) and mice (14), making them an excellent model with which to elucidate the molecular signaling mechanisms involved in the skeletal muscle alterations in response to spaceflight. Thus exposure to microgravity results in decrements in muscle strength and endurance in both humans and rodents that in humans could adversely affect the ability to perform critical functions while in space, particularly on long missions.

Ground-based models also have been used to simulate aspects of the spaceflight environment in rodents. Hindlimb suspension has been used extensively in rats and mice to approximate the musculoskeletal unloading as well as the cephalic fluid shift that occurs during spaceflight (27). Rats subjected to hindlimb suspension approximate the reduction in cross-sectional area and slow-to-fast fiber type transitions in slow-twitch muscles that are observed in the rats flown in space (reviewed in Ref. 26). However, other studies suggest that spaceflight induces gene expression changes in rat skeletal muscle that are not elicited in suspension models (28). This is in some ways unsurprising given that there are significant differences between spaceflight and hindlimb suspension. Specifically, one major difference is that organisms flown in space typically experience a few hours of reloading before death, and this reloading may also influence early events such as signal transduction and gene transcription.

Recently, much progress has been made in elucidating the intracellular signaling pathways governing muscle mass and fiber type. Flux through the phosphatidylinositol 3 (PI3)-kinase/Akt/mammalian target of rapamycin (mTOR) pathway has been demonstrated to increase protein synthesis (reviewed in Ref. 11) and to decrease protein degradation through inactivation of members of the FoxO family of transcription factors (39). This pathway is activated by insulin-like growth factor-I (IGF-I) and can also be influenced by mechanical loading (15), but definitive studies have yet to be done on the effects of spaceflight on activity and/or expression of genes involved in PI3-kinase/Akt/mTOR signaling.

Skeletal muscle mass is also negatively regulated by the myostatin pathway. Inactivation of myostatin through naturally occurring or targeted mutations results in dramatic muscle hypertrophy in mice, cattle, and humans (16). Myostatin acts in a paracrine/autocrine loop by binding to members of the activin RII receptor family, most notably ActRIIb (19), and decreases myoblast proliferation (42), differentiation (18, 33), and protein synthesis (40) while increasing expression of protein degradation genes in vitro (24). Increased myostatin expression has been reported in a number of conditions of unloading atrophy, including spaceflight and hindlimb suspension in rats (7, 17), but to date little is known about the effects of unloading on other myostatin-associated genes.

In addition, an exciting recent finding has been the identification of microRNAs and their role in regulating gene expression in a number of species and a variety of tissues. MicroRNAs are small (∼21 nt), evolutionarily conserved, noncoding RNAs with the ability to regulate gene expression on a wide scale through the binding of the microRNA to sequence-specific sites within target mRNAs (6). Several microRNAs have been identified that regulate skeletal muscle differentiation and gene expression (6) and that demonstrate muscle fiber type-specific expression (23). Currently, however, little is known regarding the expression of microRNAs during conditions of unloading atrophy such as spaceflight.

The purpose of the present study was threefold. First, to better understand the gene expression changes underlying the functional decrements in response to spaceflight, we examined the effects of spaceflight on mRNA expression in gastrocnemius muscle by utilizing both a global microarray approach and a candidate gene approach using quantitative real-time polymerase chain reaction (QRT-PCR). Second, we compared the effects of spaceflight on global mRNA expression to those of hindlimb suspension for a comparable length of time. Third, because all spaceflight studies contain the possible confound of several hours of reloading upon reentry and landing, we evaluated the effects of the 3.5 h of reloading (experienced by the spaceflight mice before death) on mRNA expression of suspended mouse muscle. We hypothesized that both spaceflight and hindlimb suspension would be associated with shifts in mRNA expression that would favor myostatin signaling and decrease PI3-kinase/Akt/mTOR signaling and thus would favor protein degradation over protein synthesis and decreased myogenic cell proliferation and differentiation over cell growth, and that reloading would reverse these effects.

MATERIALS AND METHODS

Microgravity exposure.

All animals were housed individually in animal enclosure modules [AEM; developed by National Aeronautics and Space Administration (NASA) Ames Research Center (Moffett Field, CA)] with access to food and water ad libitum. The muscles used in the spaceflight study came from the same mice as those previously described (12, 14, 30): female C57BL/6J mice were divided into AEM ground controls and spaceflight (SF) groups. Animals were maintained on a 12:12-h light-dark cycle with all animal care, treatment, and housing density within National Institutes of Health guidelines. All protocols were approved by the appropriate institutions (University of Colorado at Boulder, Amgen, Inc., NASA Ames Research Center, and NASA Kennedy Space Center). SF animals were flown on the middeck of the space shuttle Endeavour (STS-108/UF-1) for 11 days and 19 h beginning December 5, 2001. AEM animals were housed at Cape Canaveral Air Force Station's Hangar L orbital environmental simulator for 11 days and 19 h with conditions mimicking the shuttle's middeck temperature, humidity, and CO2 levels. As mentioned previously, the module containing the SF mice experienced an increase in temperature during flight that may have contributed to differences in eating and drinking behavior, as previously mentioned (12, 14, 30). AEM and SF animals were killed at 77 days of age, with SF animals being killed 3.5–4 h after landing. The triceps surae (gastrocnemius, plantaris, and soleus) of SF and AEM animals were removed and attached to corkboard with optimal cutting temperature medium (Sakura Finetek, Torrance, CA) and then frozen in liquid nitrogen-cooled isopentane for other purposes (12). The medial section of the gastrocnemius was isolated for RNA analysis. Gastrocnemius muscle (n = 4/group) was subsequently used for DNA microarray analysis.

Hindlimb suspension and reloading.

For the ground-based unloading/reloading studies, mice were divided into three groups (n = 5/group): control, suspension, and suspension plus reloading. The control group was normally housed three mice per cage. The suspension and reloading groups were suspended for 12 days to match the STS-108 flight. Tail suspension was performed using previously described methods (37). Briefly, the mice were suspended at an approximately 30° angle by taping the tail to a plastic dowel attached to a swivel apparatus. This apparatus was attached to a guide wire running the length of the cage. The mice were able to access all areas of the cage because of the swivel between the dowel and the wire. The suspension animals were killed while suspended at the end of the 12 days. Suspended mice were unsuspended for 3.5 h at the end of the 12 days of suspension and allowed to engage in normal loading for this period to match the time between landing and death for the STS-108 mice. All animals were 77 days of age at the time of death. This protocol was approved by the Clemson University Animal Care and Use Committee. The triceps surae of all animals were removed by dissection, weighed, and snap frozen in liquid nitrogen.

Preparation of total RNA.

Total RNA was prepared from tissue samples by Trizol extraction followed by purification using RNeasy columns (Qiagen, Valencia, CA). Quality and purity of the RNA preparations were assessed by spectrophotometric determination of the ratio of absorbance at 260/280 nm and by quantification of the ratios of 28S:18S ribosomal RNA using a Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA).

Synthesis of biotin-labeled cRNA targets and hybridization to Affymetrix GeneChips.

Synthesis of double-stranded cDNA, in vitro transcription of biotin-labeled cRNA targets, and fragmentation of target cRNA were performed as outlined by Affymetrix protocols (Affymetrix, Santa Clara, CA). Fragmented cRNA samples were hybridized overnight at 45°C to Affymetrix Mouse Expression 430A GeneChips. Posthybridization washing and phycoerythrin-streptavidin staining and fluorescence scanning were performed using Affymetrix instrumentation in accordance with manufacturer protocols. The resulting scans (CEL files) were deposited in the MUSC DNA microarray database (2); these can be accessed from the database using the project identifiers _1112123225.418684 (gastrocnemius) and _1120601249.362999 (suspension/reload). Microarray data were also deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus repository under the accession number GSE10533.

Analysis of DNA microarray data.

Gene hybridization intensities were normalized using Robust Multichip Average (5) implemented with the ArrayQuest web-based analysis system (3). Normalized hybridization values were imported into the analysis tool dChip 2005 (21) for subsequent analysis. mRNAs differentially expressed in response to spaceflight were identified using both statistical (P < 0.05 for unpaired Student's t-test) and fold change (fold change >1.5) thresholds. The fold change metric was chosen to ensure that the detected expression differences were sufficiently large to be validated by PCR-based approaches. False discovery rate, estimated by iterative comparison of randomized sample groupings, approximated 1.5% for these differential expression criteria. Differential expression in response to suspension, reloading, or the combination of suspension plus reloading was evaluated by pairwise comparison (Student's t-test, unpaired). Statistical power analysis using Power Atlas (29) determined that α (P) <0.05 provided the most robust outcomes in terms of proportion of expected differentially expressed genes discovered (>25%) versus proportion of true positives in the genes discovered (>75%). Gene functional categories were assigned following review of gene ontology annotations as well as functional information available through NCBI Entrez Gene. Hierarchical clustering of differentially expressed mRNAs was conducted using the centroid linkage model with distance metric 1-correlation.

Quantitative PCR.

Quantitative RT-PCR (QRT-PCR) was used to confirm expression changes observed in the microarray analysis. First-strand cDNA was prepared from total RNA (1 μg) using the iScript cDNA synthesis kit (Bio-Rad, Hercules, CA) according to the manufacturer's specifications. cDNA preparations were diluted to 200 μl, and 2 μl were used in 25-μl reactions with the iQ SYBR green Supermix reagent (Bio-Rad). Reactions were amplified in an iCycler real-time PCR detection system (Bio-Rad). The resulting data were analyzed with the Gene Expression Analysis macro (Bio-Rad), which derives an unscaled expression value for the gene of interest by first calculating a relative expression value and then normalizing it by division with the geometric mean of relative expression values found for control mRNAs (Hprt1 and β-actin).

QRT-PCR was also used to quantify mRNA levels of myostatin (MSTN), its receptor, ActRIIb, and its primary binding/inhibiting protein, FSTL3, in all of the SF and AEM samples, as previously described (1). The RT reaction was carried out using 0.5 μg of RNA with the cDNA Archive kit (Advanced Biosystems) according to the manufacturer's protocol. Primer and probe sets for MSTN, ActRIIb, FSTL3, and β-actin TaqMan analysis were obtained from Applied Biosystems (Foster City, CA). These real-time PCR procedures were run in triplicate to correct for variances in loading. In addition, a standard curve ranging from 10- to 0.001-μg dilutions of mouse TA cDNA was run in duplicate for every assay to produce a standard curve for quantification. All values are expressed as the mean of the triplicate measure for the experimental (MSTN, ActRIIb, or FSTL3) divided by the statistical mean of the triplicate measure of β-actin, which was run in duplex, as a normalization control. β-Actin mRNA levels were unaffected by spaceflight (data not shown).

Finally, for analysis of miR-1, miR-133a, and miR-206, RT was performed using a microRNA-specific primer according to the manufacturer's instructions (Applied Biosystems). After RT, samples were analyzed using a microRNA-specific real-time PCR kit according to the manufacturer's instructions (Applied Biosystems). All microRNA PCR results were determined using the relative (ΔΔCt) quantification method with the small nuclear U6 RNA (snRNA U6) run in singleplex as the normalization control, since this is a more appropriate control for normalization of microRNA levels than an mRNA. Levels of the snRNA U6 were not significantly altered by spaceflight (data not shown). For both the QRT-PCR for the myostatin-associated mRNAs and for the miRs, an independent t-test was used to evaluate statistical significance, with α < 0.05 taken as significant.

RESULTS

Muscle atrophy induced by unloading.

As previously reported (14), body weight was significantly decreased in the SF mice compared with both preflight (18.3 g preflight vs. 17.4 g postflight) and AEM mice (19.4 g for AEM vs. 17.4 g for SF). Similarly, body mass in suspended mice was also lower postsuspension compared with presuspension (18.2 g presuspension vs. 17.8 g postsuspension) and unsuspended mice (19.2 g for unsuspended vs. 17.8 g for suspended). Moreover, gastrocnemius samples from mice flown on STS-108 showed significant reductions in fiber cross-sectional area and a decrease in citrate synthase activity, as reported previously (14). Analysis of triceps surae collected from the hindlimb suspended mice and suspended plus reloaded mice in the present study confirmed that the ground simulation unloading models also experienced muscle atrophy: mean weights of the triceps surae taken from the ground-based models were 96.6 ± 1.1, 85.4 ± 2.9, and 88.4 ± 1.3 mg for control, suspended, and reloaded tissues, respectively, with both suspended and reloaded weights significantly smaller than control weights (P < 0.05).

Global microarray analysis.

Analysis of DNA microarray data for STS-108 gastrocnemius revealed that 272 mRNAs were differentially expressed compared with AEM controls (fold change >1.5, P < 0.05), with 151 mRNAs upregulated and 121 mRNAs downregulated (Fig. 1A and Table 1). To evaluate how these 272 mRNAs responded to unloading and reloading stimuli, we performed DNA microarray analysis on calf muscle prepared from ground-based mouse models subjected to hindlimb suspension or hindlimb suspension plus brief reloading. Hierarchical clustering of the expression data for the 272 spaceflight-affected mRNAs revealed many similarities in the responses to spaceflight and hindlimb suspension (Fig. 1, A–C). Many of the mRNAs significantly affected by spaceflight showed similar trends in response to suspension, as indicated by the heat map (Fig. 1A), particularly in the area highlighted by the bracket at right. Examination of microarray data revealed that most mRNAs upregulated by spaceflight showed upregulation in response to hindlimb suspension (113 of the 151; Fig. 1B), whereas most mRNAs downregulated by spaceflight were also downregulated by suspension (92 of 121; Fig. 1C). Furthermore, it can be seen that reloading counteracted the effect of hindlimb suspension for many of these mRNAs such that mRNAs that were upregulated by suspension tended to be downregulated by reloading, and vice versa; only 22 upregulated mRNAs and 8 downregulated mRNAs did not follow this pattern (Fig. 1, B and C).

Fig. 1.

Fig. 1.mRNAs identified as differentially expressed in spaceflight samples and in response to ground-based unloading models. A: expression patterns of the 272 mRNAs differentially expressed in mouse gastrocnemius in response to spaceflight. Hierarchical clustering was performed for the 272 mRNAs affected by spaceflight using expression data for spaceflight (SF) and the animal enclosure module (AEM) controls as well as for hindlimb suspension (Sus), hindlimb suspension plus reloading (Sus+Rel), and control (Con). The bracket at right highlights mRNAs that generally indicate a pattern in which hindlimb suspension mimics the effects of spaceflight and for which reloading alleviates the effect of hindlimb suspension. B and C: relative expression patterns of mRNAs regulated by spaceflight and by suspension in a consistent manner. B shows expression patterns for 113 mRNAs upregulated by spaceflight and by suspension; and C shows expression patterns for 92 mRNAs downregulated by spaceflight and by suspension. D: Venn diagram schematically depicting the numbers and associations of mRNAs found to be significantly affected by suspension, reloading, suspension plus reloading, and/or spaceflight stimuli.


Table 1. mRNAs significantly affected in spaceflight gastrocnemius

Functional CategoryGene NameGene SymbolSpaceflight
Suspension
Reloading
Suspension + Reloading
FCP valueFCP valueFCP valueFCP value
ApoptosisCell death-inducing DFFA-like effector cCidec3.10.0011.40.031−1.10.3451.20.227
Growth arrest and DNA-damage-inducible 45 βGadd45b2.40.0131.30.124−1.20.2201.10.453
Nuclear factor of κ light polypeptide gene enhancer in B cells inhibitor, αNfkbia2.30.0051.10.359−1.20.087−1.10.235
Transformation related protein 53 inducible nuclear protein 1Trp53inp12.30.0071.30.072−1.40.033−1.10.353
Pleckstrin homology-like domain, family A, member 1Phlda12.00.004−1.10.4171.10.2941.00.732
Lectin, galactose binding, soluble 3Lgals31.90.0191.10.320−1.00.7441.10.571
Angiopoietin-like 4Angptl41.70.003−1.10.410−1.10.622−1.20.118
Programmed cell death 4Pdcd41.6<0.0011.10.375−1.10.367−1.00.916
Structural maintenance of chromosomes 6Smc6−1.60.005−1.10.6661.20.5021.10.613
Phosphatase and tensin homologPten−1.80.004−1.00.8751.10.3221.10.095
RNA binding motif protein 5Rbm5−2.2<0.0011.10.6991.10.5531.20.179
ATPaseATPase, aminophospholipid transporter, class I, type 8A, 1Atp8a1−1.50.0191.20.2831.20.2011.40.029
ATPase, Na+/K+ transporting, β1 polypeptideAtp1b1−1.6<0.001−1.40.0361.50.0301.00.742
ATPase, H+ transporting, lysosomal V1 subunit AAtp6v1a−1.6<0.001−1.00.9401.30.0121.30.012
Cell cycleB-cell translocation gene 2, antiproliferativeBtg23.00.011−1.30.0072.00.0101.60.049
Cyclin-dependent kinase inhibitor 1A (P21)Cdkn1a2.7<0.0011.10.322−1.10.5901.00.909
v-maf musculoaponeurotic fibrosarcoma oncogene family, protein FMaff2.60.024−1.10.3341.20.0801.10.198
SNF1-like kinaseSnf1lk1.90.004−1.10.1981.30.0111.20.029
Proviral integration site 3Pim31.80.0041.10.444−1.00.9361.10.248
Sestrin 1Sesn11.70.007−1.10.7581.20.3071.10.281
RIKEN cDNA 1190002H23 gene1190002H23Rik1.6<0.001−1.30.0621.40.0291.10.285
Transducer of ERBB2, 2Tob21.60.0101.00.972−1.30.011−1.30.017
ELK3, member of ETS oncogene familyElk3−1.6<0.0011.00.874−1.00.9411.00.873
v-crk sarcoma virus CT10 oncogene homologCrk−1.6<0.001−1.10.5591.10.4961.00.994
RB1-inducible coiled-coil 1Rb1cc1−1.6<0.001−1.10.4301.20.0371.10.199
Avian musculoaponeurotic fibrosarcoma AS42 oncogene homologMaf−1.90.020−1.10.1171.10.363−1.10.401
Protein phosphatase 1, regulatory subunit 10Ppp1r10−1.90.001−1.20.0161.10.205−1.10.265
Cell growth regulationPhosphatidylinositol 3-kinase, regulatory subunit, polypeptide 1Pik3r12.10.0061.20.1241.10.6801.30.095
HtrA serine peptidase 1Htra12.00.004−1.10.2491.20.0941.10.593
Cysteine-rich protein 61Cyr611.8<0.0011.10.6551.80.0101.80.007
HtrA serine peptidase 3Htra31.5<0.0011.10.082−1.20.041−1.10.309
Suppressor of cytokine signaling 2Socs21.50.0321.20.460−1.40.091−1.20.591
Cysteine rich transmembrane BMP regulator 1Crim1−1.60.002−1.10.6571.30.0331.20.261
CircadianPeriod homolog 1Per13.00.0091.00.725−1.20.098−1.20.045
Cryptochrome 2 (photolyase-like)Cry21.60.0401.10.0941.10.3621.20.021
CytoskeletalNeuroepithelial cell transforming gene 1Net11.8<0.0011.10.111−1.00.5241.10.245
Tensin like C1 domain-containing phosphataseTenc11.70.0031.00.8341.01.0001.00.779
Myristoylated alanine-rich protein kinase C substrateMarcks−1.60.003−1.20.1671.10.692−1.20.120
Tubulin, α1ATuba1a−1.60.0091.10.3591.00.9411.10.272
Caldesmon 1Cald1−1.60.001−1.30.3641.40.1591.20.131
LIM domain and actin binding 1Lima1−1.60.0051.00.6111.10.2211.10.044
Cysteine-rich protein 2Crip2−1.60.007−1.10.1601.20.0161.10.012
Tropomyosin 1, αTpm1−1.70.010−1.00.9611.20.2781.20.244
PDZ and LIM domain 7Pdlim7−1.70.003−1.10.5941.10.4791.00.928
Tropomodulin 3Tmod3−1.70.004−1.00.8401.10.4741.10.559
LIM and SH3 protein 1Lasp1−1.70.004−1.10.3181.20.1471.10.424
Kinesin family member 1BKif1b−1.9<0.001−1.20.1671.10.461−1.10.243
DevelopmentNeuronatinNnat2.40.014−1.00.988−1.30.391−1.30.393
Tumor necrosis factor, α-induced protein 2Tnfaip21.70.0351.00.908−1.40.076−1.40.042
Activin A receptor, type 1BAcvr1b1.60.0021.30.003−1.30.0061.00.962
Jumonji, AT-rich interactive domain 2Jarid21.50.047−1.10.5951.00.976−1.10.414
Vascular endothelial growth factor AVegfa−1.50.004−1.30.0191.20.059−1.10.307
Kit ligandKitl−1.6<0.0011.00.9331.20.0661.20.022
Frizzled homolog 9Fzd9−1.60.003−1.50.0101.10.523−1.40.022
DNA/chromatin maintenanceH3 histone, family 3BH3f3b1.70.0101.10.427−1.00.8591.10.413
High mobility group box 2Hmgb21.60.006−1.10.5481.00.931−1.00.492
Tousled-like kinase 2 (Arabidopsis)Tlk21.60.0161.10.4561.10.6161.20.028
Centromere protein BCenpb−1.60.007−1.10.2071.10.514−1.00.883
Bromodomain adjacent to zinc finger domain, 1BBaz1b−1.70.005−1.10.2401.10.1341.00.971
Extracellular matrix/adhesionA disintegrin-like and metallopeptidase with thrombospondin type 1 motif, 15Adamts151.90.016−1.00.8811.00.8651.00.959
Basal cell adhesion moleculeBcam−1.50.018−1.10.2561.10.0571.10.215
Collagen, type III, α1Col3a1−1.50.019−1.50.0611.20.239−1.20.176
Collagen, type IV, α3 binding proteinCol4a3bp−1.50.0051.10.360−1.00.8661.10.371
Procollagen-proline, 2-oxoglutarate 4-dioxygenase, α1 polypeptideP4ha1−1.60.040−1.20.0191.10.394−1.10.207
Tissue inhibitor of metalloproteinase 4Timp4−1.70.008−1.00.9831.40.1781.40.048
Integrin β6Itgb6−1.80.004−1.20.0541.20.196−1.00.713
A disintegrin-like and metallopeptidase with thrombospondin type 1 motif, 1Adamts1−1.9<0.001−1.30.0692.40.0041.90.018
Collagen, type I, α2Col1a2−1.90.035−2.30.0661.70.201−1.40.258
Collagen, type I, α1Col1a1−2.10.021−2.20.0651.60.211−1.30.275
MetabolismAdiponectin, C1Q and collagen domain containingAdipoq2.8<0.0011.90.006−1.40.1041.30.282
Cysteine dioxygenase 1, cytosolicCdo12.20.0031.90.014−1.50.0531.30.226
Fatty acid synthaseFasn2.1<0.0013.10.009−1.30.2992.40.029
TransketolaseTkt2.00.0072.00.003−1.20.1931.70.010
Indolethylamine N-methyltransferaseInmt1.9<0.0011.40.051−1.20.3331.20.136
ATP citrate lyaseAcly1.80.0161.40.2421.10.7241.60.025
Liver glycogen phosphorylasePygl1.70.0081.30.040−1.10.2941.20.193
Phosphoenolpyruvate carboxykinase 1, cytosolicPck11.70.0021.40.005−1.20.1391.20.063
Sulfotransferase family 1A, phenol-preferring, member 1Sult1a11.70.0071.30.023−1.20.0351.10.618
Aldehyde dehydrogenase 1 family, member L1Aldh1l11.7<0.0011.20.038−1.20.0111.10.485
3-Phosphoglycerate dehydrogenasePhgdh1.6<0.0011.20.056−1.20.0631.00.914
Amine oxidase, copper containing 3Aoc31.60.0021.20.287−1.20.1921.00.932
Carbonic anhydrase 4Car41.6<0.0011.30.023−1.30.0241.00.699
Dehydrogenase/reductase (SDR family) member 3Dhrs31.50.0031.10.261−1.10.4681.00.966
Hexose-6-phosphate dehydrogenaseH6pd1.5<0.0011.10.115−1.10.1381.00.989
LeptinLep1.50.0231.10.534−1.20.093−1.10.337
Carbonic anhydrase 14Car141.5<0.0011.10.343−1.00.8021.10.409
Phosphate cytidylyltransferase 1, choline, α-isoformPcyt1a−1.60.014−1.00.6931.00.6971.00.905
CDP-diacylglycerol synthase 2Cds2−1.60.020−1.10.5031.40.0831.30.057
High-density lipoprotein binding proteinHdlbp−1.80.008−1.10.3731.10.3921.00.903
RIKEN cDNA 2610507B11 gene2610507B11Rik−2.20.007−1.20.1191.10.470−1.00.866
MitochondrionCytochrome P450, family 2, subfamily e, polypeptide 1Cyp2e12.90.0152.1<0.001−1.30.1031.60.005
Solute carrier family 25, member 33Slc25a332.20.026−1.10.2651.20.1701.10.613
Glutamate-ammonia ligaseGlul2.0<0.0011.20.3801.10.5001.40.039
CytochromeP-450, family 2, subfamily d, polypeptide 22Cyp2d221.80.0111.20.030−1.10.0901.10.305
Solute carrier family 25, member 25Slc25a251.80.017−1.30.136−1.10.403−1.40.078
Uncoupling protein 3 (mitochondrial, proton carrier)Ucp31.80.018−1.30.4271.30.3381.00.991
UDP glucuronosyltransferase 1 family, polypeptide A1Ugt1a11.7<0.0011.30.023−1.10.4771.20.056
Proline dehydrogenaseProdh1.7<0.0011.20.027−1.10.0491.10.411
Solute carrier family 35, member B3Slc35b31.60.015−1.20.0241.10.160−1.00.934
Translocator proteinTspo1.50.016−1.10.6421.00.778−1.00.816
Zinc finger protein 106Zfp106−1.60.017−1.10.367−1.10.204−1.20.142
3-Hydroxybutyrate dehydrogenase, type 1Bdh1−2.70.027−2.6<0.0011.30.138−2.00.003
mRNA processingSplicing factor, arginine/serine−rich 7Sfrs72.30.0131.30.053−1.20.1241.10.344
Transcribed locus: accession no. BB2033482.00.0081.10.421−1.10.227−1.00.488
Splicing factor 3b, subunit 1Sf3b11.80.014−1.00.7391.10.4061.10.532
Splicing factor, arginine/serine-rich 5 (SRp40, HRS)Sfrs51.70.0171.10.562−1.10.251−1.10.450
Transcribed locus: accession no. AA717264−1.60.0071.00.634−1.10.240−1.00.384
Poly (A) polymerase αPapola−1.60.009−1.10.4461.20.1341.10.627
DEAD (Asp-Glu-Ala-Asp) box polypeptide 6Ddx6−1.80.0151.10.6061.20.3381.30.165
Oxidative damageHaptoglobinHp1.90.0221.40.082−1.30.1131.10.544
ProteaseComplement factor D (adipsin)Cfd2.60.0161.90.007−1.40.0421.30.139
Signal transductionMetallothionein 2Mt29.30.0131.10.7891.10.5681.20.189
Zinc finger protein 36Zfp363.10.0051.00.994−1.00.661−1.00.666
Metallothionein 1Mt12.50.009−1.10.1171.10.514−1.00.735
SplA/ryanodine receptor domain and SOCS box containing 1Spsb12.10.012−1.10.5391.10.603−1.00.934
ADP-ribosylation factor-like 4DArl4d2.10.0051.00.855−1.00.631−1.00.570
Membrane protein, palmitoylated 3Mpp32.00.0111.30.048−1.00.8821.30.033
ERBB receptor feedback inhibitor 1Errfi12.00.001−1.20.3041.40.0251.20.113
Membrane-spanning 4-domains, subfamily A, member 6BMs4a6b1.90.0451.10.099−1.10.077−1.00.780
WD repeat domain 26Wdr261.9<0.0011.00.7911.00.7231.10.467
Serine (or cysteine) peptidase inhibitor, clade E, member 1Serpine11.9<0.0011.00.699−1.10.423−1.10.520
Ras homolog gene family, member BRhob1.80.005−1.00.8751.10.6161.00.713
Mitogen-activated protein kinase kinase kinase 6Map3k61.80.0071.00.7681.00.8191.10.584
Interleukin 6 receptor, αIl6ra1.80.0031.20.015−1.10.0951.10.261
G protein-coupled receptor 146Gpr1461.60.002−1.00.592−1.00.615−1.00.348
Heparin-binding EGF-like growth factorHbegf1.60.0161.10.5631.80.0021.90.001
Leucine-rich repeat (in FLII) interacting protein 1Lrrfip11.60.0151.10.576−1.20.336−1.10.556
Adrenergic receptor, β2Adrb21.50.022−1.10.635−1.10.641−1.20.167
IQ motif containing GTPase activating protein 1Iqgap1−1.50.003−1.00.7021.20.0261.20.065
Activin A receptor, type 1Acvr1−1.50.010−1.10.4861.00.954−1.10.240
WNK lysine-deficient protein kinase 1Wnk1−1.50.0391.00.9321.10.3811.10.422
Insulin receptor substrate-1Irs1−1.6<0.001−1.30.0431.40.0261.00.638
Frizzled homolog 7Fzd7−1.6<0.001−1.10.4011.10.403−1.00.870
Sprouty-related, EVH1 domain containing 2Spred2−1.60.0051.10.550−1.10.3211.00.997
Follistatin-like 1Fstl1−1.60.0131.00.778−1.20.034−1.20.034
Casein kinase 2, α1 polypeptideCsnk2a1−1.60.011−1.10.6411.40.2151.20.475
Regulator of G protein signaling 5Rgs5−1.70.005−1.10.2971.10.2751.00.691
Growth factor receptor-bound protein 10Grb10−1.70.016−1.10.2931.00.964−1.10.235
RAB22A, member RAS oncogene familyRab22a−1.7<0.0011.00.953−1.00.478−1.00.231
Calcium/calmodulin-dependent protein kinase IIαCamk2a−1.70.023−1.10.4101.10.2191.10.617
Tetraspanin 2Tspan2−1.70.002−1.10.5221.10.5381.00.977
Regulator of G-protein signaling 4Rgs4−1.7<0.001−1.00.9611.10.1981.10.263
Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, γYwhag−1.80.020−1.20.2051.20.195−1.00.885
A kinase (PRKA) anchor protein (yotiao) 9Akap9−1.80.0061.20.040−1.10.0881.10.181
Inhibitor of Bruton agammaglobulinemia tyrosine kinaseIbtk−1.80.002−1.20.2741.10.585−1.10.537
Phosphodiesterase 4D interacting protein (myomegalin)Pde4dip−2.10.004−1.00.9511.20.3391.20.353
Angiotensin receptor-like 1Agtrl1−2.6<0.001−1.00.914−1.10.653−1.10.618
Transcription factorV-maf musculoaponeurotic fibrosarcoma oncogene family, protein BMafb3.20.0081.10.634−1.30.004−1.20.096
SRY-box containing gene 4Sox43.20.0011.10.286−1.10.4001.00.917
CCAAT/enhancer binding protein, δCebpd3.0<0.0011.40.125−1.40.1451.00.902
Forkhead box O1Foxo12.3<0.0011.10.4991.20.2921.30.044
Nuclear factor, interleukin 3, regulatedNfil32.10.005−1.30.0861.80.0161.40.091
Kruppel-like factor 10Klf102.00.0091.00.980−1.20.367−1.20.608
MyogeninMyog1.90.021−1.00.827−1.10.595−1.10.222
Paired-like homeodomain transcription factor 2Pitx21.80.001−1.10.1651.00.776−1.10.139
Nuclear receptor interacting protein 1Nrip11.80.0321.40.0201.30.1411.70.008
CCAAT/enhancer binding protein, αCebpa1.7<0.0011.40.073−1.30.0621.10.782
Kruppel-like factor 15Klf151.70.0021.10.170−1.20.064−1.10.244
Transcription factor CP2-like 1Tcfcp2l11.7<0.0011.40.020−1.10.3871.30.031
Forkhead box P1Foxp11.60.021−1.00.5921.00.985−1.00.333
CCAAT/enhancer binding protein, βCebpb1.60.0011.00.9661.00.7731.00.755
Kruppel-like factor 9Klf91.60.0201.10.323−1.10.4491.00.530
Prolactin regulatory element bindingPreb1.5<0.001−1.00.8471.10.5121.10.379
Jun oncogeneJun1.50.006−1.10.2952.2<0.0012.0<0.001
HairlessHr1.50.0011.10.3321.10.5531.10.287
Muscle, skeletal, receptor tyrosine kinaseMusk1.50.0301.00.5091.20.1181.20.039
Kruppel-like factor 3 (basic)Klf3−1.50.0021.10.1931.00.5141.20.106
Nuclear receptor subfamily 1, group D, member 2Nr1d2−1.50.018−1.20.3671.20.414−1.00.781
MAD homolog 5Smad5−1.50.0251.10.3221.10.2271.20.052
Nuclear factor I/XNfix−1.60.0401.00.8751.00.9371.00.819
Nuclear factor of activated T cells, cytoplasmic, calcineurin-dependent 3Nfatc3−1.60.023−1.10.4381.10.4111.00.876
Ankyrin repeat and SOCS box-containing protein 15Asb15−1.60.0121.10.6161.00.7731.10.390
Inhibitor of DNA binding 3Id3−1.60.006−1.20.0411.20.114−1.10.332
Kruppel-like factor 7 (ubiquitous)Klf7−1.70.003−1.10.3291.20.1041.10.372
SRY-box containing gene 7Sox7−1.7<0.0011.00.6321.20.0741.20.021
Paired related homeobox 1Prrx1−1.7<0.001−1.30.0821.30.0691.00.800
Kruppel-like factor 4 (gut)Klf4−1.80.013−1.40.0111.30.127−1.10.337
Mesenchyme homeobox 2Meox2−1.8<0.001−1.10.2441.20.0141.10.297
Peroxisome proliferator activated receptor αPpara−1.80.007−1.30.0991.30.029−1.00.853
E26 avian leukemia oncogene 1, 5′ domainEts1−1.90.043−1.30.2281.20.277−1.10.671
Zinc finger protein 191Zfp191−1.90.029−1.10.3671.00.995−1.10.363
Inhibitor of DNA binding 1Id1−2.0<0.001−1.80.0011.50.013−1.20.232
Nuclear factor I/BNfib−2.10.0031.00.7511.10.6431.10.382
SRY-box containing gene 18Sox18−2.2<0.001−1.30.0151.10.136−1.10.272
Peroxisome proliferative activated receptor, γ, coactivator 1αPgc1a−2.40.002−1.00.9252.70.0322.70.033
Transporter activityLipocalin 2Lcn23.10.0061.00.907−1.10.408−1.10.617
Double C2, βDoc2b2.2<0.0011.10.362−1.10.6281.10.451
Solute carrier family 10 (sodium/bile acid cotransporter family), member 6Slc10a61.80.0031.00.7501.00.9091.00.502
Solute carrier family 7 (cationic amino acid transporter, y+ system), member 8Slc7a81.80.0061.30.051−1.10.2801.10.340
Retinol binding protein 4, plasmaRbp41.70.0011.40.023−1.20.0931.10.225
Solute carrier family 1 (neutral amino acid transporter), member 5Slc1a51.60.0021.40.011−1.20.2141.20.177
Solute carrier family 38, member 2Slc38a2−1.50.008−1.10.5732.8<0.0012.6<0.001
Solute carrier family 6 (neurotransmitter transporter, creatine), member 8Slc6a8−1.50.012−1.3<0.0011.10.092−1.10.106
Transferrin receptorTfrc−1.60.029−1.40.0331.10.509−1.20.106
Target of myb1-like 2 (chicken)Tom1l2−1.70.029−1.00.7421.00.903−1.00.858
UbiqutinRing finger protein 144BRnf144b1.80.0061.00.562−1.00.9281.00.574
F-box protein 32Fbxo321.60.001−1.00.934−1.40.152−1.40.116
Ubiquitin-conjugating enzyme E2 variant 2Ube2v21.60.0031.10.209−1.10.3831.00.443
Male-specific lethal 2-like 1Msl2l11.60.0021.00.436−1.10.178−1.10.344
3′-Phosphoadenosine 5′-phosphosulfate synthase 2Papss21.50.0091.10.174−1.00.5721.10.375
Neural precursor cell expressed, developmentally downregulated gene 4-likeNedd4l1.5<0.0011.10.6001.00.9931.10.449
Cullin 2Cul21.50.0391.10.109−1.10.3911.00.531
Tripartite motif-containing 25Trim251.50.003−1.00.957−1.00.782−1.00.693
Topoisomerase I binding, arginine/serine-richTopors1.5<0.0011.00.795−1.10.128−1.10.104
CUE domain containing 1Cuedc1−1.60.019−1.10.7111.10.780−1.00.918
Ubiquitin protein ligase E3AUbe3a−1.70.001−1.10.3891.20.1211.10.074
Neuralized-like homologNeurl−2.10.002−1.00.763−1.10.528−1.10.372
Unfolded protein bindingFK506 binding protein 5Fkbp51.8<0.0011.30.082−1.20.2271.10.588
FUS interacting protein (serine-arginine rich) 1Fusip11.70.0181.10.257−1.10.4311.10.618
Sarcolemma associated proteinSlmap−1.60.011−1.30.1431.10.643−1.20.237
Serine (or cysteine) peptidase inhibitor, clade H, member 1Serpinh1−1.70.017−1.40.0461.40.0421.00.883
Other/UnknownDNA-damage-inducible transcript 4Ddit46.60.0371.20.599−1.30.382−1.10.650
SarcolipinSln2.80.005−1.20.4271.10.680−1.10.725
Stanniocalcin 2Stc22.8<0.0011.10.270−1.20.041−1.10.126
Thyroid hormone responsive SPOT14 homolog (Rattus)Thrsp2.5<0.0012.00.009−1.30.1561.50.086
Accession BC0040652.4<0.0011.10.2591.10.4961.10.100
TCDD-inducible poly(ADP-ribose) polymeraseTiparp2.30.0041.00.835−1.20.019−1.20.340
Cold inducible RNA binding proteinCirbp2.20.0041.30.049−1.10.4701.20.035
Son cell proliferation proteinSon2.10.0051.10.073−1.00.6731.10.156
RNA binding motif protein 39Rbm392.10.0061.20.057−1.10.2961.10.608
Plasma membrane associated protein, S3-12S3-122.10.0101.00.809−1.00.513−1.00.660
ResistinRetn2.00.0041.60.007−1.30.1171.20.285
Interferon, α-inducible protein 27Ifi272.00.0261.40.038−1.40.1101.00.991
DNA segment, Chr 19, Wayne State University 162, expressedD19Wsu162e2.00.0031.10.102−1.00.7731.10.283
Phenylethanolamine-N-methyltransferasePnmt2.00.0061.20.132−1.10.4201.10.340
Amyotrophic lateral sclerosis 2 (juvenile) chromosome region, candidate 13Als2cr131.8<0.0011.10.348−1.20.057−1.10.178
Protein-tyrosine sulfotransferase 2Tpst21.80.0031.00.838−1.10.150−1.10.431
HIV-1 Rev binding protein-likeHrbl1.8<0.0011.30.027−1.30.020−1.00.914
Melanocortin 2 receptor accessory proteinMrap1.80.0401.10.277−1.00.6331.10.522
Small nucleolar RNA host gene (nonprotein coding) 1Snhg11.80.0061.10.5801.10.2381.20.124
Matrix-remodelling associated 8Mxra81.70.0011.00.2501.00.9191.10.352
Hypothetical gene LOC433762LOC4337621.70.0051.10.577−1.00.7131.00.816
Accession BG1434611.70.0141.20.184−1.20.222−1.00.927
Immediate early response 5Ier51.60.004−1.10.1771.90.0101.80.014
RIKEN cDNA 2410002F23 gene2410002F23Rik1.6<0.0011.10.452−1.10.077−1.10.335
Complement component 3C31.60.0121.10.357−1.00.9541.10.308
Arrestin domain containing 4Arrdc41.60.0041.10.154−1.10.1251.00.853
Similar to zinc finger and BTB domain containing 161.60.0141.00.9941.30.0821.30.061
Ribosomal protein S25Rps251.60.0171.10.019−1.10.0651.00.460
Prominin 1Prom11.60.0011.20.072−1.20.1501.00.749
SERTA domain containing 2Sertad21.60.0081.10.2601.00.8851.10.171
Coiled-coil domain containing 80Ccdc801.60.0111.00.9701.10.5501.10.493
CD1d1 antigenCd1d11.5<0.0011.10.063−1.10.3651.10.305
Translocated promoter regionTpr−1.50.0121.00.8461.10.1671.20.214
MOB1, Mps One Binder kinase activator-like 1BMobkl1b−1.50.0131.10.329−1.10.002−1.00.490
DNA segment, Chr 3, University of California at Los Angeles 1D3Ucla1−1.50.011−1.20.1131.20.120−1.00.865
Zinc finger protein 354BZfp354b−1.50.0071.10.3781.00.9171.10.201
Cd200 antigenCd200−1.60.0011.00.8851.00.9831.00.785
RIKEN cDNA 6720475J19 gene6720475J19Rik−1.60.018−1.20.2051.10.264−1.00.675
Transmembrane 4 superfamily member 1Tm4 sf1−1.60.006−1.00.812−1.00.746−1.00.510
ThrombomodulinThbd−1.6<0.001−1.00.9081.40.0031.30.003
N-deacetylase/N-sulfotransferase 1Ndst1−1.60.0221.00.784−1.10.603−1.00.761
Accession AK084266−1.60.024−1.10.1051.10.300−1.00.379
Transmembrane protein 140Tmem140−1.60.008−1.30.2361.50.0831.20.277
Carboxypeptidase DCpd−1.60.0161.00.887−1.00.9041.00.992
Endothelin converting enzyme 1Ece1−1.60.005−1.00.8661.10.4281.00.782
Accession BG976607−1.70.016−1.10.218−1.00.471−1.20.069
Myeloid/lymphoid or mixed-lineage leukemia translocated to, 3Mllt3−1.7<0.001−1.10.1721.20.0691.10.271
Pleckstrin homology, Sec7 and coiled-coil domains 1Pscd1−1.70.042−1.10.3331.00.641−1.00.558
Zinc finger protein 422, related sequence 1Zfp422-rs1−1.7<0.0011.00.924−1.10.343−1.10.636
Coiled-coil domain containing 50Ccdc50−1.70.020−1.10.5031.00.723−1.00.687
Eukaryotic translation initiation factor 4, γ1Eif4 g1−1.70.025−1.00.6801.10.4601.00.658
Musashi homolog 2Msi2−1.70.0031.10.2521.20.0101.30.005
DNA segment, Chr 4, Wayne State University 53, expressedD4Wsu53e−1.70.002−1.00.702−1.00.600−1.10.481
RIKEN cDNA 2010106G01 gene2010106G01Rik−1.80.0081.10.2551.10.4321.20.079
FolliculinFlcn−1.9<0.001−1.00.3771.10.0841.10.198
Dual specificity phosphatase 18Dusp18−1.9<0.001−1.20.1651.00.811−1.20.180
Myeloid-associated differentiation markerMyadm−2.1<0.001−1.20.0131.30.0201.00.579
Serum deprivation responseSdpr−2.1<0.001−1.30.1381.20.355−1.10.181
SH3 domain protein D19Sh3d19−2.20.001−1.00.8581.10.3901.10.381
QuakingQk−2.2<0.001−1.30.1521.30.0991.10.461
Solute carrier family 9, member 3 regulator 2Slc9a3r2−2.5<0.001−1.00.8261.10.0361.10.059

Pairwise comparison data are shown for spaceflight (vs. animal enclosure module), suspension (vs. control), reloading (vs. suspension), and suspension plus reloading (vs. control). Gene names in bold type show significant differences for 1 or more of the ground-based stimuli. FC, fold change; P value, unpaired Student's t-test.

Further examination of the expression data revealed that 88 of the mRNAs significantly affected by spaceflight were also significantly affected by one or more of the ground-based stimuli (Table 1). Specifically, 46 of the spaceflight-affected mRNAs were significantly affected by suspension, 42 by reloading, and 37 by the combination of suspension plus reloading, with overlaps apparent among these groups (Fig. 1D and Table 1). The majority of these mRNAs (72 of 88) showed a response to suspension that matched the direction of change observed in the SF samples. Furthermore, most of these mRNAs (74 of 88) showed expression patterns that were consistent with sensitivity to both unloading and reloading stimuli, with reloading causing a change in expression that opposed the change induced by hindlimb suspension.

Validation of differential expression

Expression changes observed by microarray experimentation were confirmed for selected representative genes by QRT-PCR analysis. Generally, QRT-PCR results closely matched the microarray findings in terms of magnitude and direction of change (Table 2). Specifically, QRT-PCR supported the effect of spaceflight for all nine mRNAs. Furthermore, all mRNAs that were identified by microarray analysis as significantly affected by suspension or reloading also showed significant changes by QRT-PCR (Table 2).

Table 2. Validation of microarray mRNA expression measurements by QRT-PCR

GeneForward Primer/Reverse Primer (5′-3′)Spaceflight
Suspension
Reloading
ArrayQRT-PCRArrayQRT-PCRArrayQRT-PCR
AdipoqAACGACTCTACATTTACTGGCT/CCTGGTGTATGGGCTATG2.8*3.71.9*1.9−1.4−1.2
FasnCCACACCTGCCACCAGAG/CGCGTGATGGTCTGTTTGTA2.1*1.93.1*5.3 −1.3−1.2
Id1CGCACTGAGGGACCAGATG/GCTGCGGTAGTGTCTTTCCCA−2.0*−2.8−1.8*−1.31.5*1.3
Pgc1-αGGCCAGTGCTACCTGAGAGA/GGATGACCGAAGTGCTTGT−2.4*−1.8−1.0−1.32.7*3.2
Cyr61AATAAGATTAGTTGGGCCAAGT/CGGGAAGGAATGAGAGG1.8*1.51.11.61.8*2.0
Col1a2CACCTGGTCCTGTTGGAAGT/ACCAGGGAAGCCAGTCATAC−1.9*−3.3−2.3−2.01.71.5
Foxo1GCGCATAGCACCAAGTCTTC/AGCGTGACACAGGGCAT2.4*2.31.11.21.11.7
Pitx2CGCCGACTCCTCCGTA/ATACTGGCAAGCACTCAGGTT1.8*2.1−1.11.11.0−1.2
Zfp36AAGCAAATAGCCAAAGCCA/AACGGAGTAGACGACCCAATA3.1*3.3−1.01.2−1.01.1
QRT-PCR control primers
ActbTGACAGGATGCAGAAGGAGA/CGCTCAGGAGGAGCAATG
Hprt1TCCTCCTCAGACCGCTTTT/CCTGGTTCATCATCGCTAATC

Fold change values determined by microarray analysis (Array) or quantitative RT-PCR (QRT-PCR) are shown for spaceflight (vs. animal enclosure module), suspension (vs. control), and reloading (vs. suspension).

*P < 0.05, significant difference by microarray analysis (fold change >1.5, unpaired t-test).

P < 0.05, significant difference by QRT-PCR analysis (unpaired t-test).

Spaceflight and PI3-kinase/Akt/mTOR signaling genes.

Spaceflight was associated with significant shifts in expression of mRNAs involved with the PI3-kinase/Akt/mTOR pathway. Levels of mRNA of the PI3-kinase regulatory subunit polypeptide 1, pi3kr1/p85α (2.13-fold), the forkhead box O1 (FoxO1) transcription factor (2.31-fold), the muscle-specific ubiquitin ligase F-box only protein 32 (MAFbx/atrogin1; 1.61-fold), and the ubiquitin-conjugating enzyme E2 variant 2 (1.57-fold) genes were significantly increased in SF muscle compared with control. Conversely, insulin receptor substrate-1 (IRS-1) mRNA levels were significantly decreased 1.55-fold in SF gastrocnemius.

Spaceflight and expression of myostatin-associated genes.

Expression of mRNA for the myostatin-associated gene activin Ib receptor (Acvr1b; also known as ALK4), was significantly increased 1.59-fold with spaceflight (Table 1). The myostatin mRNA is not represented on the 430A GeneChip. However, the primary myostatin receptor ActRIIb and the myostatin binding/inhibitory protein FSTL3 were both present on the GeneChip, although neither were detected as differentially expressed. Nevertheless, given the established importance of these three genes in regulating muscle mass, we examined their mRNA expression using QRT-PCR as an alternate means of analysis. As shown by QRT-PCR, myostatin mRNA levels tended to increase in the SF samples, whereas mRNA levels of FSTL3 tended to decrease in the SF samples, although neither reached statistical significance (Fig. 2). Conversely, mRNA levels of the myostatin receptor ActRIIb were significantly decreased in SF gastrocnemius muscle compared with AEM control (Fig. 2).

Fig. 2.

Fig. 2.Quantitative RT-PCR analysis of myostatin, ActRIIb, and FSTL3 in ground control (AEM) and SF gastrocnemius muscle. Bars represent means ± SE for n = 4 mice per group. All values were normalized to a β-actin duplexed internal normalizing control. A: myostatin mRNA levels. B: ActRIIb mRNA levels. C: FSTL3 mRNA levels. Myostatin increased and ActRIIb and FSTL3 mRNA levels decreased in SF gastrocnemius, but only the changes in ActRIIb mRNA levels were statistically significant. *P < 0.05 compared with AEM control.


Spaceflight and other growth-associated genes.

Spaceflight also affected mRNA levels for other growth regulatory genes. For example, mRNA levels of suppressor of cytokine signaling 2 (Socs2), p21, and Btg2 were significantly increased in gastrocnemius in response to spaceflight (1.51-, 2.68- and 2.95-fold, respectively). Other pathways that have been linked to muscle growth and adaptation include the TNF-α/NF-κB signaling pathway and the calcineurin/nuclear factor of activated T cells (NFAT) pathway. The TNF-α downstream target TNF-α-induced protein 2 mRNA was significantly increased in SF gastrocnemius (1.68-fold), but the NF-κB inhibitor nuclear factor κB light chain gene enhancer in B cells inhibitor α (Nfkbia/IκBα) mRNA was also significantly increased (2.28-fold). Moreover, mRNA levels of the NFAT cytoplasmic, calcineurin-dependent 3 (Nfatc3) transcription factor were significantly decreased (−1.57-fold) in SF gastrocnemius. Finally, mRNAs for three members of the CAAT/enhancer binding protein (C/EBP) family of transcription factors, C/EBP-α (1.72-fold), C/EBP-β (1.58-fold), and C/EBP-δ (3.03-fold) were significantly increased in SF gastrocnemius.

Spaceflight and fiber type genes.

The gastrocnemius muscle of the STS-108 mice examined in this study exhibited a decrease in citrate synthase activity consistent with a shift toward a less oxidative phenotype (14). Among the 272 mRNAs differentially expressed in response to spaceflight, two encode proteins associated with regulating oxidative phenotype through activation of oxidative gene expression. Peroxisome proliferator activated receptor (PPAR)-α and the PPAR-γ coactivator 1α (PGC1-α) were both significantly decreased (1.83-fold and −2.38-fold, respectively).

mRNAs affected by spaceflight, suspension, and reloading.

As mentioned above, the majority of mRNAs affected by spaceflight showed a qualitatively similar response to suspension, and reloading tended to diminish or even reverse the effects of suspension. This was true for pi3kr1/p85α, p21, Socs2, C/EBP-α and -γ, and PPAR-α, which were affected by spaceflight and suspension in a matching fashion and for which reloading counteracted the suspension effect (Table 1). Moreover, 13 other mRNAs affected by spaceflight were significantly by suspension and by reloading (Fig. 1D). Notably, all 13 of these mRNAs showed counteracting responses to suspension versus reloading. Furthermore, all except one showed a significant response to suspension that matched the direction of the significant response to spaceflight, the exception being Btg2, which significantly increased with both spaceflight and reloading but significantly decreased with suspension. Among the 13 mRNAs was the key growth and metabolism regulator IRS-1, which significantly decreased in response to spaceflight and suspension and significantly increased in response to reloading. Conversely, Acvr1b/ALK4 mRNA levels were significantly increased by spaceflight and suspension but significantly decreased in response to reloading (Table 1). Finally, levels of PGC-1α mRNA were significantly decreased by spaceflight and, although not significantly affected by suspension, were significantly increased by reloading relative to suspension (Table 1).

MicroRNA affected by spaceflight.

As shown in Fig. 3, levels of miR-206 relative to the endogenous control snRNA U6 were significantly decreased by nearly 50% in gastrocnemius muscle from SF mice compared with AEM controls (Fig. 3C). In contrast, neither miR-1 nor miR-133a were significantly affected by spaceflight, although both showed a trend toward a decrease, with the decrease in miR-133a nearly reaching significance (Fig. 3, A and B). However, the ratio of miR-1 to miR-133a was significantly increased in SF gastrocnemius (Fig. 3D).

Fig. 3.

Fig. 3.TaqMan analysis of microRNA (miRNA) expression in ground control (AEM) and SF gastrocnemius muscle for miR-1 (A), miR-133a (B), miR-206 (C), and the ratio of miR-1 to miR-133a (D). Total RNA was subjected to reverse transcription using gene-specific primers miR-1, miR-133a, miR-206, and small nuclear RNA (snRNA) U6 and then quantified using quantitative real-time PCR. Fold change was determined using the ΔΔCt relative quantitation method with snRNA U6 as the endogenous control gene. *P < 0.05 compared with AEM control.


DISCUSSION

The data previously published for mice from STS-108 showed that murine skeletal muscle exhibited several adaptations in response to spaceflight (14). All fiber types analyzed showed a reduction in cross-sectional area. Moreover, the gastrocnemius showed a decrease in citrate synthase activity consistent with a shift toward a less oxidative phenotype. In this study, analysis of mRNA expression by both gene array and QRT-PCR has shown that spaceflight affects mRNA levels of genes involved in numerous cellular processes, including metabolism, cell cycle, apoptosis, and cytoskeletal and mitochondrial function. Moreover, spaceflight impacts expression of several signal transduction pathways affecting both muscle growth, including those related to growth signaling via PI3-kinase/Akt/mTOR, calcineurin/NFAT, and myostatin and those related to muscle fiber type determination, such as PPAR-α and PGC1-α.

Spaceflight and genes associated with PI3-kinase/Akt/mTOR signaling.

Numerous studies have now demonstrated that a primary pathway regulating muscle protein metabolism, and therefore muscle growth, is the PI3-kinase/Akt/mTOR pathway. This pathway represents a primary node for regulating both protein synthesis and protein degradation in response to various stimuli (11, 20, 32). In the present study, we observed that spaceflight altered mRNA expression of several genes involved in this pathway: the PI3-kinase regulatory subunit p85α, which negatively regulates PI3-kinase signaling by sequestration of the insulin receptor substrate adaptor-1 (13), was increased, IRS-1 was decreased, and mRNA levels of the FoxO1 transcription factor were increased, as were those of one of its primary targets, MAFbx/atrogin1. Together these data support the hypothesis that spaceflight results in a shift in gene expression that favors reduced flux through the PI3-kinase/Akt/mTOR pathway.

Spaceflight and myostatin signaling.

Both the gene array and QRT-PCR data suggested that spaceflight is associated with a shift in gene expression favoring increased myostatin signaling. Levels of myostatin mRNA tended to increase in spaceflight gastrocnemius (Fig. 2), consistent with previous findings from both spaceflight and hindlimb suspension muscle (7, 17). Similarly, mRNA levels of the myostatin binding/inhibiting protein gene FSTL3 also tended to decrease in spaceflight gastrocnemius. Expression of the myostatin-ActRIIb partnering receptor Acvr1b /ALK4 was significantly increased by both spaceflight and hindlimb suspension (and decreased by reloading), although ActRIIb mRNA levels significantly decreased in spaceflight. With the exception of the decrease in ActRIIb levels, these shifts, if translated into changes in the associated proteins, would be consistent with a shift toward increased myostatin signaling in spaceflight muscle. Consistent with this interpretation, levels of p21, a cell cycle inhibitor that is a target of myostatin signaling (42), were significantly upregulated in response to spaceflight, consistent with a previous gene array study on spaceflight rat muscle (41).

In addition, there is considerable feedback between the PI3-kinase/Akt/mTOR/FoxO pathway and the myostatin pathway. Myostatin activates protein degradation by inhibiting Akt phosphorylation of the FoxO transcription factors, resulting in increased expression of ubiquitin ligases (24). Moreover, myostatin expression is itself positively regulated by the FoxO transcription factors (1). Thus the changes in mRNA expression associated with reduced Akt activity described above, as well as the increase in FoxO1 mRNA expression, may also have contributed to the increase in myostatin transcription observed in this study.

Spaceflight and other growth pathways.

Spaceflight significantly affected mRNA levels of several other genes associated with growth regulation. mRNA encoding Socs2, a suppressor of cytokine and growth hormone signaling, was increased, as were mRNAs encoding the cell cycle inhibitors p21 and Btg2. These are consistent with a shift in gene expression favoring reduced cell proliferation, consistent with results from previous gene array studies on spaceflight and hindlimb suspension (41, 44). In addition, expression of mRNAs encoding TNF-α-induced protein 2 and Nfatc3 were altered in ways consistent with shifts in activity of the NF-κB and calcineurin/NFAT pathways, respectively, that would favor muscle atrophy. Finally, there were increases in mRNAs encoding three members of the C/EBP family of transcription factors: C/EBP-α, C/EBP-β, and C/EBP-δ. Recent evidence has suggested that members of the C/EBP family are regulated by glucocorticoids, particularly C/EBP-δ (25), and their binding activity is increased during periods of muscle catabolism (31). Given the present results, the potential role of these transcription factors in regulating muscle atrophy in response to spaceflight and/or other modes of unloading needs to be further addressed.

Spaceflight and metabolism.

Muscle adaptation to disuse is often characterized by fiber type alterations in both contractile and metabolic proteins and protein isoforms, and previous gene array studies on rat soleus have demonstrated metabolic gene expression shifts with spaceflight and/or hindlimb suspension (8, 38, 43). In the present study, expression of both PPAR-α and PGC1-α mRNA was decreased in spaceflight gastrocnemius and was significantly increased in suspended plus reloaded muscle. Dapp et al. (8) also showed a significant decrease in PPAR-α expression with 7 days of hindlimb suspension in the soleus, and expression of both these transcription factors is increased by endurance exercise (36). In addition, PPAR-α agonists induce a fiber-selective transcriptional response in rats (9), whereas transgenic expression of PGC1-α in type II fibers has been shown to stimulate changes in muscle morphology, gene expression, and transition toward a type I fiber phenotype in mice (22). Thus the decrease in expression of these two genes may contribute to the shift toward a less oxidative phenotype observed in response to spaceflight (14).

Spaceflight and microRNA expression.

One of the most exciting breakthroughs in modern biology has been the identification of microRNAs and their role in regulating cell function. MicroRNAs regulate RNA stability and/or translation in a wide variety of species and cells, including muscle (reviewed in Ref. 6). More recent data have demonstrated that several microRNAs are expressed in a muscle type-specific and/or differentiation-specific manner (6) and that these microRNAs appear to play critical roles in muscle gene expression, proliferation, and differentiation (6). In the present study, we found that expression of miR-206 was significantly decreased in spaceflight gastrocnemius compared with ground-based control, whereas levels of miR-1 and miR-133a were not significantly altered (Fig. 3, AC). Interestingly, however, the ratio of miR-1 to miR-133a was significantly increased with spaceflight (Fig. 3D). The significance of this is that miR-1 and miR-133a act in opposition to one another with respect to their effect on skeletal muscle, with miR-1 inhibiting proliferation and miR-133a promoting proliferation via their effects on HDAC4 and serum response factor, respectively (43). This increase in the ratio of miR-1 to miR-133a suggests a microRNA expression profile favoring inhibition of cell proliferation, a finding that is consistent with the shift in expression of key regulatory genes such as p21 and Socs2 that appears to favor decreased proliferation, as well. To our knowledge these are the first data to examine the levels of microRNAs in muscle in response to spaceflight. To date, only one study has investigated the in vivo effects of altering muscle activity/load on skeletal muscle microRNA expression (23). Seven days of functional overload (FO) of the plantaris resulted in a ∼50% decrease in miR-1 and miR-133a combined with an ∼18-fold increase in the miR-206 precursor but no change in mature miR-206 expression (23).

At present, the functional consequences of decreased miR-206 expression with microgravity exposure are unknown. miR-206 expression is known to be higher in slow, oxidative muscle than in fast, glycolytic muscle (23), and thus the decrease in spaceflight muscle might reflect the slow-to-fast shift that characterizes this state. Furthermore, miR-206 expression is low in proliferating myoblasts and increased during myoblast differentiation (6), and overexpression of miR-206 promotes cell cycle withdrawal and differentiation (6), and thus the decrease in miR-206 observed in the present study may also represent a shift in differentiation capabilities of the satellite cells and/or muscle fibers. FSTL1 is a target gene negatively regulated by miR-206 (34); however, in the present study FSTL1 mRNA levels were decreased in response to spaceflight, which is at odds with the results observed for miR-206 expression. Thus at present it is unclear what role, if any, miR-206 plays in regulating muscle gene expression during spaceflight.

Comparison of spaceflight and ground-based unloading models.

A majority of the 272 mRNAs affected by spaceflight (171 mRNAs) in the present study showed a similar trend toward increased or decreased expression with hindlimb suspension (Fig. 1B). Furthermore, 88 of these mRNAs were significantly affected by one or more of the ground-based stimuli. Generally, there was an overriding reciprocal relationship between the effect of suspension versus the effect of reloading: only 14 of the 88 mRNAs affected by one or more of the ground-based stimuli did not follow this pattern. Although the sets of mRNAs affected significantly by the ground-based stimuli are only partially overlapping, our results suggest that significance for any of the stimuli appears to act as a strong indicator of sensitivity to unloading/reloading factors. A previous study of rat muscle (28) reported that spaceflight resulted in at least an 8-fold change to 257 total genes, whereas hindlimb suspension resulted in at least an 8-fold change to just 74 genes (see Table 3 of Ref. 28), percentages fairly similar to those observed in the present study, suggesting that spaceflight appears to affect more genes than ground-based unloading in both mice and rats.

Thus the present data suggest that spaceflight elicits qualitatively similar but quantitatively different effects on muscle mRNA levels compared with hindlimb suspension. It is also possible that some or all of the mRNAs that were significantly affected by spaceflight but not by hindlimb suspension were affected by a component of spaceflight other than the unloading due to microgravity, such as radiation or hypergravity upon reentry, among others, which may have potentiated the effects of microgravity-induced unloading on spaceflight muscle. In addition, as mentioned in materials and methods, slightly different tissue preparations were used; for the spaceflight studies, the medial portion of the gastrocnemius was used, whereas in the suspension studies, the entire calf was used. It is therefore possible that at least some of the differences in mRNA expression between spaceflight and suspension are a consequence of differences in how one portion of the gastrocnemius responds compared with the combined muscles of the calf. However, although the three calf muscles differ appreciably in their fiber type composition in normal animals, we previously demonstrated (14) that all three calf muscles undergo qualitatively similar shifts in fiber size in response to spaceflight. Thus it is unlikely that changes in mRNA expression were masked by differential shifts in signaling pathways regulating fiber size between the different calf muscles. Nevertheless, differences in responses between spaceflight and suspension may reflect a dilution of the mRNAs affected in the medial gastrocnemius used in the spaceflight studies by the additional calf musculature used in the suspension studies.

Spaceflight and reloading.

An important consideration for analysis of any spaceflight data is the reloading period that occurs after landing and before any potential examination. In the present study, we found that 42 of the 272 mRNAs affected in spaceflight samples were significantly affected by the reloading stimulus, and most of these showed an antagonistic response to reloading relative to suspension. This is consistent with a previous study on acute (12 h) suspension and 4-h reloading on rat soleus muscle (4), which also showed that reloading tended to counteract the effects of suspension. In the present study, notable genes whose mRNA was significantly affected by reloading relative to spaceflight and/or suspension were IRS-1, Acvr1b/ALK4, and PGC1-α.

As described above, the general effect of reloading was to counteract the effect of suspension. Among 13 genes independently affected by suspension and by reloading, only 1 exhibited a significant difference for combined suspension plus reloading relative to the control state, and only 3 showed a net change of >1.06-fold compared with control. Therefore, reloading tended to diminish and thus mask the effect of suspension. Despite this, mRNA expression changes in spaceflight samples were qualitatively most similar to the effects of suspension, suggesting that the predominance of gene expression changes occurring in spaceflight samples are likely due to the 12 days of unloading and not the 3.5 h of reloading. Nevertheless, our finding that 3.5 h of reloading were sufficient to cause some transcriptional changes, combined with the finding that reloading may mask some portion of the unloading response, indicates that care must be taken in the interpretation of mRNA expression data taken from reloaded spaceflight samples.

In summary, spaceflight induces changes in expression of mRNAs involved in a number of cellular processes and in particular appears to alter the mRNA expression of many genes involved with regulating muscle growth and fiber type. Surprisingly, few of these were significantly altered in both spaceflight and hindlimb unloading, which suggests key differences in how muscles respond on the transcriptional level to these stimuli. Moreover, reloading is associated with changes in several mRNAs, some of which respond differentially compared with spaceflight and/or hindlimb suspension. At present it is not clear to what extent these mRNA changes are translated into differences in protein levels for these genes. Future studies are needed to follow up on exploring the effects of spaceflight and suspension/reloading on protein levels of a selected group of candidate genes from the present study.

GRANTS

This work was funded by NASA Grant NCC8-242 (to L. S. Stodieck; BioServe Space Technologies), Amgen, Inc., South Carolina Space Grant Consortium NASA Grant (SCSGCNG) NCC5-575 (to T. A. Bateman), National Space Biomedical Research Institute NASA Grant NCC 9-58 (to T. A. Bateman), a NASA Graduate Student Research Program fellowship from Kennedy Space Center (to E. R. Bandstra), National Center for Research Resources Grants P20 RR016434 (to W. S. Argraves), SCSGCNG NCC5-575 (to W. S. Argraves), and South Carolina INBRE Bioinformatics Core P20 RR16461 (to W. S. Argraves).

FOOTNOTES

  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We acknowledge Mark Rupert (BioServe/U of Colorado) for excellent payload management, support from Dr. Beverly Girten [National Aeronautics and Space Administration (NASA) Ames Research Center], and Ramona Bober (Kennedy Space Center) and resource support from the Medical University of South Carolina Proteogenomics Facility, including exceptional technical assistance from Victor Fresco.

REFERENCES

  • 1 Allen DL, Unterman TG. Regulation of myostatin expression and myoblast differentiation by FoxO and SMAD transcription factors. Am J Physiol Cell Physiol 292: C188–C199, 2007.
    Link | ISI | Google Scholar
  • 2 Argraves GL, Barth JL, Argraves WS. The MUSC DNA Microarray Database. Bioinformatics 19: 2473–2474, 2003.
    Crossref | ISI | Google Scholar
  • 3 Argraves GL, Jani S, Barth JL, Argraves WS. ArrayQuest: a web resource for the analysis of DNA microarray data. BMC Bioinformatics 6: 287, 2005.
    Crossref | ISI | Google Scholar
  • 4 Bey L, Akunuri N, Zhao P, Hoffman EP, Hamilton DG, Hamilton MT. Patterns of global gene expression in rat skeletal muscle during unloading and low-intensity ambulatory activity. Physiol Genomics 13: 157–167, 2003.
    Link | ISI | Google Scholar
  • 5 Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185–193, 2003.
    Crossref | PubMed | ISI | Google Scholar
  • 6 Callis TE, Deng Z, Chen JF, Wang DZ. Muscling through the microRNA world. Exp Biol Med (Maywood) 233: 131–138, 2008.
    Crossref | PubMed | ISI | Google Scholar
  • 7 Carlson CJ, Booth FW, Gordon SE. Skeletal muscle myostatin mRNA expression is fiber-type specific and increases during hindlimb unloading. Am J Physiol Regul Integr Comp Physiol 277: R601–R606, 1999.
    Link | ISI | Google Scholar
  • 8 Dapp C, Schmutz S, Hoppeler H, Fluck M. Transcriptional reprogramming and ultrastructure during atrophy and recovery of mouse soleus muscle. Physiol Genomics 20: 97–107, 2004.
    Link | ISI | Google Scholar
  • 9 De Souza AT, Cornwell PD, Dai X, Caguyong MJ, Ulrich RG. Agonists of the peroxisome proliferator-activated receptor alpha induce a fiber-type-selective transcriptional response in rat skeletal muscle. Toxicol Sci 92: 578–586, 2006.
    Crossref | ISI | Google Scholar
  • 10 Fitts RH, Desplanches D, Romatowski JG, Widrick JJ. Spaceflight effects on single skeletal muscle fiber function in the rhesus monkey. Am J Physiol Regul Integr Comp Physiol 279: R1546–R1557, 2000.
    Link | ISI | Google Scholar
  • 11 Glass DJ. Signalling pathways that mediate skeletal muscle hypertrophy and atrophy. Nat Cell Biol 5: 87–90, 2003.
    Crossref | PubMed | ISI | Google Scholar
  • 12 Gridley DS, Nelson GA, Peters LL, Kostenuik PJ, Bateman TA, Morony S, Stodieck LS, Lacey DL, Simske SJ, Pecaut MJ. Genetic models in applied physiology: selected contribution: effects of spaceflight on immunity in the C57BL/6 mouse. II. Activation, cytokines, erythrocytes, and platelets. J Appl Physiol 94: 2095–2103, 2003.
    Link | ISI | Google Scholar
  • 13 Hara K, Yonezawa K, Sakaue H, Ando A, Kotani K, Kitamura T, Kitamura Y, Ueda H, Stephens L, Jackson TR, Hawkins PT, Dhand R, Clark AE, Holman GD, Waterfield MD, Kasuga M. 1-Phosphatidylinositol 3-kinase activity is required for insulin-stimulated glucose transport but not for RAS activation in CHO cells. Proc Natl Acad Sci USA 91: 7415–7419, 1994.
    Crossref | PubMed | ISI | Google Scholar
  • 14 Harrison BC, Allen DL, Girten B, Stodieck LS, Kostenuik PJ, Bateman TA, Morony S, Lacey D, Leinwand LA. Skeletal muscle adaptations to microgravity exposure in the mouse. J Appl Physiol 95: 2462–2470, 2003.
    Link | ISI | Google Scholar
  • 15 Hornberger TA, Sukhija KB, Chien S. Regulation of mTOR by mechanically induced signaling events in skeletal muscle. Cell Cycle 5: 1391–1396, 2006.
    Crossref | PubMed | ISI | Google Scholar
  • 16 Joulia-Ekaza D, Cabello G. The myostatin gene: physiology and pharmacological relevance. Curr Opin Pharmacol 7: 310–315, 2007.
    Crossref | PubMed | ISI | Google Scholar
  • 17 Lalani R, Bhasin S, Byhower F, Tarnuzzer R, Grant M, Shen R, Asa S, Ezzat S, Gonzalez-Cadavid NF. Myostatin and insulin-like growth factor-I and -II expression in the muscle of rats exposed to the microgravity environment of the NeuroLab space shuttle flight. J Endocrinol 167: 417–428, 2000.
    Crossref | PubMed | ISI | Google Scholar
  • 18 Langley B, Thomas M, Bishop A, Sharma M, Gilmour S, Kambadur R. Myostatin inhibits myoblast differentiation by down-regulating MyoD expression. J Biol Chem 277: 49831–49840, 2002.
    Crossref | PubMed | ISI | Google Scholar
  • 19 Lee SJ, McPherron AC. Regulation of myostatin activity and muscle growth. Proc Natl Acad Sci USA 98: 9306–9311, 2001.
    Crossref | PubMed | ISI | Google Scholar
  • 20 Leger B, Cartoni R, Praz M, Lamon S, Deriaz O, Crettenand A, Gobelet C, Rohmer P, Konzelmann M, Luthi F, Russell AP. Akt signalling through GSK-3beta, mTOR and Foxo1 is involved in human skeletal muscle hypertrophy and atrophy. J Physiol 576: 923–933, 2006.
    Crossref | ISI | Google Scholar
  • 21 Li C, Hung Wong W. Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol 2: RESEARCH0032, 2001.
    Crossref | PubMed | Google Scholar
  • 22 Lin J, Wu H, Tarr PT, Zhang CY, Wu Z, Boss O, Michael LF, Puigserver P, Isotani E, Olson EN, Lowell BB, Bassel-Duby R, Spiegelman BM. Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature 418: 797–801, 2002.
    Crossref | PubMed | ISI | Google Scholar
  • 23 McCarthy JJ, Esser KA. MicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophy. J Appl Physiol 102: 306–313, 2007.
    Link | ISI | Google Scholar
  • 24 McFarlane C, Plummer E, Thomas M, Hennebry A, Ashby M, Ling N, Smith H, Sharma M, Kambadur R. Myostatin induces cachexia by activating the ubiquitin proteolytic system through an NF-kappaB-independent, FoxO1-dependent mechanism. J Cell Physiol 209: 501–514, 2006.
    Crossref | PubMed | ISI | Google Scholar
  • 25 Menconi M, Fareed M, O'Neal P, Poylin V, Wei W, Hasselgren PO. Role of glucocorticoids in the molecular regulation of muscle wasting. Crit Care Med 35: S602–608, 2007.
    Crossref | PubMed | ISI | Google Scholar
  • 26 Morey-Holton ER, Hill EL, Souza KA. Animals and spaceflight: from survival to understanding. J Musculoskelet Neuronal Interact 7: 17–25, 2007.
    Google Scholar
  • 27 Musacchia XJ, Fagette S. Weightlessness simulations for cardiovascular and muscle systems: validity of rat models. J Gravit Physiol 4: 49–59, 1997.
    PubMed | Google Scholar
  • 28 Nikawa T, Ishidoh K, Hirasaka K, Ishihara I, Ikemoto M, Kano M, Kominami E, Nonaka I, Ogawa T, Adams GR, Baldwin KM, Yasui N, Kishi K, Takeda S. Skeletal muscle gene expression in space-flown rats. FASEB J 18: 522–524, 2004.
    Crossref | PubMed | ISI | Google Scholar
  • 29 Page GP, Edwards JW, Gadbury GL, Yelisetti P, Wang J, Trivedi P, Allison DB. The PowerAtlas: a power and sample size atlas for microarray experimental design and research. BMC Bioinformatics 7: 84, 2006.
    Crossref | ISI | Google Scholar
  • 30 Pecaut MJ, Nelson GA, Peters LL, Kostenuik PJ, Bateman TA, Morony S, Stodieck LS, Lacey DL, Simske SJ, Gridley DS. Genetic models in applied physiology: selected contribution: effects of spaceflight on immunity in the C57BL/6 mouse. I. Immune population distributions. J Appl Physiol 94: 2085–2094, 2003.
    Link | ISI | Google Scholar
  • 31 Penner G, Gang G, Sun X, Wray C, Hasselgren PO. C/EBP DNA-binding activity is upregulated by a glucocorticoid-dependent mechanism in septic muscle. Am J Physiol Regul Integr Comp Physiol 282: R439–R444, 2002.
    Link | ISI | Google Scholar
  • 32 Reynolds TH, Bodine SC, Lawrence JC Jr. Control of Ser2448 phosphorylation in the mammalian target of rapamycin by insulin and skeletal muscle load. J Biol Chem 277: 17657–17662, 2002.
    Crossref | PubMed | ISI | Google Scholar
  • 33 Rios R, Carneiro I, Arce VM, Devesa J. Myostatin is an inhibitor of myogenic differentiation. Am J Physiol Cell Physiol 282: C993–C999, 2002.
    Link | ISI | Google Scholar
  • 34 Rosenberg MI, Georges SA, Asawachaicharn A, Analau E, Tapscott SJ. MyoD inhibits Fstl1 and Utrn expression by inducing transcription of miR-206. J Cell Biol 175: 77–85, 2006.
    Crossref | PubMed | ISI | Google Scholar
  • 35 Roy RR, Baldwin KM, Edgerton VR. Response of the neuromuscular unit to spaceflight: what has been learned from the rat model. Exerc Sport Sci Rev 24: 399–425, 1996.
    Crossref | Google Scholar
  • 36 Russell AP, Feilchenfeldt J, Schreiber S, Praz M, Crettenand A, Gobelet C, Meier CA, Bell DR, Kralli A, Giacobino JP, Deriaz O. Endurance training in humans leads to fiber type-specific increases in levels of peroxisome proliferator-activated receptor-gamma coactivator-1 and peroxisome proliferator-activated receptor-alpha in skeletal muscle. Diabetes 52: 2874–2881, 2003.
    Crossref | PubMed | ISI | Google Scholar
  • 37 Simske SJ, Broz JJ, Fleet ML, Schmeister TA, Gayles EC, Luttges MW. Contribution of dietary and loading changes to the effects of suspension on mouse femora. J Exp Zool 269: 277–285, 1994.
    Crossref | Google Scholar
  • 38 Stein T, Schluter M, Galante A, Soteropoulos P, Tolias P, Grindeland R, Moran M, Wang T, Polansky M, Wade C. Energy metabolism pathways in rat muscle under conditions of simulated microgravity. J Nutr Biochem 13: 471, 2002.
    Crossref | PubMed | ISI | Google Scholar
  • 39 Stitt TN, Drujan D, Clarke BA, Panaro F, Timofeyva Y, Kline WO, Gonzalez M, Yancopoulos GD, Glass DJ. The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. Mol Cell 14: 395–403, 2004.
    Crossref | PubMed | ISI | Google Scholar
  • 40 Taylor WE, Bhasin S, Artaza J, Byhower F, Azam M, Willard DH Jr, Kull FC Jr, Gonzalez-Cadavid N. Myostatin inhibits cell proliferation and protein synthesis in C2C12 muscle cells. Am J Physiol Endocrinol Metab 280: E221–E228, 2001.
    Link | ISI | Google Scholar
  • 41 Taylor WE, Bhasin S, Lalani R, Datta A, Gonzalez-Cadavid NF. Alteration of gene expression profiles in skeletal muscle of rats exposed to microgravity during a spaceflight. J Gravit Physiol 9: 61–70, 2002.
    PubMed | Google Scholar
  • 42 Thomas M, Langley B, Berry C, Sharma M, Kirk S, Bass J, Kambadur R. Myostatin, a negative regulator of muscle growth, functions by inhibiting myoblast proliferation. J Biol Chem 275: 40235–40243, 2000.
    Crossref | PubMed | ISI | Google Scholar
  • 43 Van Rooij E, Marshall WS, Olson EN. Toward microRNA-based therapeutics for heart disease: the sense in antisense. Circ Res 103: 919–928, 2008.
    Crossref | PubMed | ISI | Google Scholar
  • 44 Wittwer M, Fluck M, Hoppeler H, Muller S, Desplanches D, Billeter R. Prolonged unloading of rat soleus muscle causes distinct adaptations of the gene profile. FASEB J 16: 884–886, 2002.
    Crossref | PubMed | ISI | Google Scholar

AUTHOR NOTES

  • Address for reprint requests and other correspondence: J. L. Barth, Dept. of Cell Biology and Anatomy, Medical Univ. of South Carolina, Charleston, SC 29425 (e-mail: )