The embryonic genes Dkk3, Hoxd8, Hoxd9 and Tbx1 identify muscle types in a diet-independent and fiber-type unrelated way
BMC Genomics volume 11, Article number: 176 (2010)
The mouse skeletal muscle is composed of four distinct fiber types that differ in contractile function, number of mitochondria and metabolism. Every muscle type has a specific composition and distribution of the four fiber types. To find novel genes involved in specifying muscle types, we used microarray analysis to compare the gastrocnemius with the quadriceps from mice fed a low fat diet (LFD) or high fat diet (HFD) for 8 weeks. Additional qPCR analysis were performed in the gastrocnemius, quadriceps and soleus muscle from mice fed an LFD or HFD for 20 weeks.
In mice fed the 8-week LFD 162 genes were differentially expressed in the gastrocnemius vs. the quadriceps. Genes with the strongest differences in expression were markers for oxidative fiber types (e.g. Tnni1) and genes which are known to be involved in embryogenesis (Dkk3, Hoxd8,Hoxd9 and Tbx1). Also Dkk2, Hoxa5, Hoxa10, Hoxc9, Hoxc10, Hoxc6 and Tbx15 were detectably, but not differentially expressed in adult muscle tissue. Expression of differentially expressed genes was not influenced by an 8-week or 20-week HFD. Comparing gastrocnemius, quadriceps and soleus, expression of Hoxd8 and Hoxd9 was not related with expression of markers for the four different fiber types. We found that the expression of both Hoxd8 and Hoxd9 was much higher in the gastrocnemius than in the quadriceps or soleus, whereas the expression of Dkk3 was high in quadriceps, but low in both gastrocnemius and soleus. Finally, expression of Tbx1 was high in quadriceps, intermediate in soleus and low in gastrocnemius.
We found that genes from the Dkk family, Hox family and Tbx family are detectably expressed in adult mouse muscle. Interestingly, expression of Dkk3, Hoxd8, Hoxd9 and Tbx1 was highly different between gastrocnemius, quadriceps and soleus. In fact, every muscle type showed a unique combination of expression of these four genes which was not influenced by diet. Altogether, we conclude that genes important for embryogenesis identify mouse muscle types in a diet-independent and fiber type-unrelated manner.
The mouse skeletal muscle is composed of four distinct fiber types, i.e., the type I, IIa, IIx and IIb fibers that differ in respect to contractile function, the number of mitochondria and metabolism. The slow type I fibers can sustain prolonged low power work, contain more mitochondria and exhibit higher rates of fat oxidative metabolism. The fast type IIx and IIb fibers are adapted to brief and intense contractions, contain fewer mitochondria and generate energy mainly through glycolysis. The type IIa fibers exhibit an intermediate contractile function and are oxidoglycolytic [1–4]. Every muscle type has a specific composition and distribution of the four different fiber types. The soleus contains the highest number of type I fibers, whereas both the quadriceps and gastrocnemius are muscle groups that contain predominantly type II fibers [5, 6]. However, varying percentages of type I fibers ranging from 0-45% for the quadriceps and 1-8% for the gastrocnemius are reported .
Numerous genes have been identified playing a role in the generation of more oxidative muscle types [6, 8–13]. The most extensively studied gene is calcineurin (CnA). In the C2C12 cell line, CnA up-regulates the mRNA levels of genes that are markers for the slow fiber types. Additionally, studies with CnA transgenic mice have reported an increase of oxidative fibers in the skeletal muscle [8, 9]. PGC1α is another well-established factor that induces remodeling of skeletal muscle fiber type composition. Lin et al.  showed that when PGC1α is expressed at physiological levels in muscle groups normally rich in type II fibers a fiber type conversion occurs from type II fibers to type I fibers . In line with this, Mortensen et al.  showed that PGC1α overexpression in primary rat skeletal muscle cells induces a switch towards a more oxidative fiber type. Finally, muscle-specific overexpression of PPARδ resulted in the increase of the number of oxidative fibers and in higher expression levels of several markers of oxidative metabolism [6, 12].
The first aim of the present study was to find novel genes that may play a role in specifying muscle types. Therefore we compared gene expression profiles of the gastrocnemius with the quadriceps at the level of the whole-transcriptome. Functional implications were assessed by the analyses of predefined gene sets based on Gene Ontology, biochemical, metabolic and signaling pathways. Recently, we studied the effects of a short-term high fat diet (HFD) on skeletal muscle gene and protein expression. Since both gene and protein levels of markers of the more oxidative fiber types were increased in the quadriceps of HFD mice we hypothesized that an HFD can induce a conversion towards a more oxidative fiber type via a transcriptional mechanism . The second aim of the present study therefore was to explore the diet-sensitivity of genes that are involved in the determination of muscle types. We compared expression patterns of genes that are differentially expressed between muscle types in mice fed a low fat diet (LFD) with mice fed an HFD.
Microarray analysis shows that the majority of genes involved in the determination of muscle types are not influenced by diet
To obtain more insight in genes involved in specifying muscle types we compared the gastrocnemius transcriptome with the quadriceps transcriptome of 8-week LFD mice. When using the criteria of a fold change > 1.3 and an FDR < 0.05, we found 91 genes with an increased expression and 71 genes with a decreased expression in the gastrocnemius as compared to the quadriceps. The ten genes with the strongest increased and decreased expression levels are shown in Table 1 and Fig. 1. The ten genes with the strongest increased expression included markers of oxidative fiber types (Tnnc1, Tnni1, Tnnt1 and Myh7), but also genes from the Hox family (Hoxd8, Hoxd9 and Hoxd10). In fact, Hoxd8, Hoxd9 and Hoxd10 were the genes with the strongest increased expression in the gastrocnemius. The ten genes with the strongest decreased expression contained a variety of genes such as Scd1 and Cidec, which are involved in lipogenesis, and Pck1, which plays a role in gluconeogenesis, but also the Dkk3 and Tbx1 genes, which are important during embryogenesis. A complete list of all differentially expressed genes is available in supplement 1. qPCR analyses were performed for 13 genes validating the microarray results (supplement 2).
To relate changes in gene expression to biological changes we applied ErmineJ and GSEA. By using Ermine J we identified 68 GO classes that were overrepresented in the gastrocnemius (Table 2). These overrepresented GO classes were mainly descriptors for contraction, morphogenesis, remodeling and development. Other overrepresented GO classes included descriptors for calcium homeostasis, immune function and metabolic processes. A parallel GSEA was used to identify up- or down-regulated processes. GSEA revealed that ten gene sets were down-regulated including three gene sets related to the cell cycle, two gene sets describing amino acid metabolism and two gene sets corresponding to immune function. Only the gene sets "striated muscle contraction" and "ribosome" were up-regulated in the gastrocnemius as compared with the quadriceps. Among the core-enriched genes in the gene set "striated muscle contraction" were genes that are markers for oxidative fiber types, whereas all core-enriched genes in the gene set "ribosome" were ribosomal proteins which are involved in translation (Table 3).
The diet-sensitivity of differentially expressed genes, overrepresented GO classes and regulated gene sets was studied by comparing the gastrocnemius with the quadriceps under HFD conditions (supplement 3). In the HFD mice we found 215 differentially expressed genes in the gastrocnemius as compared to the quadriceps. A total of 128 genes was overlapping with the differentially expressed genes in the LFD mice. Interestingly, these overlapping genes were all regulated in the same direction with a comparable fold change. A total of 87 genes was only differentially expressed in HFD mice, whereas 34 genes were differentially expressed in LFD mice suggesting that these genes are diet-dependent. Under HFD conditions, ErmineJ revealed that 193 GO classes were overrepresented in the gastrocnemius as compared to the quadriceps. A total of 67 GO classes showed an overlap with the overrepresented GO classes under LFD conditions. The other overrepresented GO classes corresponded also to contraction, morphogenesis, remodeling, development, calcium homeostasis, immune function and metabolic processes. GSEA showed that 11 gene sets were down-regulated including five gene sets describing cell cycle and four gene sets related to amino acid metabolism. Only one gene set "striated muscle contraction" was up-regulated in the gastrocnemius vs. the quadriceps. A total of six gene sets (five down-regulated and one up-regulated) were overlapping between the LFD and HFD condition.
To summarize, we showed that a considerable number of genes is differentially expressed in the gastrocnemius as compared with the quadriceps. The genes with the strongest regulations were markers for oxidative fiber types (Tnnc1, Tnni1, Tnnt1 and Myh7) and genes which are known to be involved in embryogenesis (Dkk3, Hoxd8, Hoxd9, Hoxd10 and Tbx1). The differentially expressed genes corresponded to a variety of process including morphogenesis and contraction, but also amino acid metabolism, cell cycle and immune function. Finally, we found more differentially expressed genes and more overrepresented GO classes under HFD conditions than under LFD conditions. However, this did not result in the identification of other regulated processes indicating that diet has only a very small effect on the difference between gastrocnemius and quadriceps.
Numerous genes involved in embryogenesis are detectably expressed in adult mouse muscle tissue
One of the most interesting findings of the comparison between gastrocnemius and quadriceps was the fact that genes of the Dkk family (Dkk3), Hox family (Hoxd8, Hoxd9 and Hoxd10) and Tbx family (Tbx1) were detectably expressed in adult mouse muscle tissue. To find out if also other genes from these three gene families are detectably expressed in adult mouse muscle tissue we analyzed the microarray signal intensities of all members, which were present on the microarray, in the gastrocnemius and quadriceps. An average microarray signal intensity > 20 was used as cut off. Fig. 2A shows that four members of the Dkk gene family were present on the microarray. In addition to Dkk3, we found that also Dkk2 was detectably expressed in both gastrocnemius and quadriceps (22.5 vs. 23.4 in gastrocnemius vs. quadriceps). From the Hox family 34 members were present on the microarray. Eight genes of the Hox family, including Hoxd8, Hoxd9 and Hoxd10, were expressed in adult muscle tissue. Hoxa5, Hoxa10, Hoxc9, Hoxc10, Hoxd8 and Hoxd9 were detectably expressed in both the gastrocnemius and quadriceps. Furthermore, Hoxc6 was detectably expressed in the quadriceps and not in the gastrocnemius (29.1 versus 16.9), whereas Hoxd10 was detectably expressed in the gastrocnemius and not in the quadriceps (36.2 versus 4.9). The gene with the highest expression level was identified as the Hoxc10 gene (Fig. 2B). Fig. 2C shows that from the Tbx family 12 members were present on the microarray. Additionally to Tbx1, we found that Tbx15 was highly expressed in the gastrocnemius as well as in the quadriceps (654.8 vs. 836.4, respectively). Altogether, these results show that some genes that are known to play an important role in embryogenesis are also detectably expressed in adult mouse muscle tissue.
The Tbx1 gene is a muscle-type specific marker
Dkk3, Hoxd8, Hoxd9 and Tbx1 were relatively high and differentially expressed in the gastrocnemius as compared to the quadriceps of mice fed an 8-week LFD. To extrapolate our result to a third muscle type we analyzed the expression level of these genes in quadriceps, gastrocnemius and soleus muscle of 20-week LFD mice and 20-week HFD mice. None of the genes was significantly influenced by diet which is in line with the earlier described observations (Fig. 3A-D). Hoxd8 and Hoxd9 were highly expressed in the gastrocnemius, but low gene expression levels were observed in the quadriceps and soleus (Fig. 3A-B). Dkk3 was significantly higher expressed in the quadriceps than in the gastrocnemius and soleus. Although not significantly, Dkk3 gene expression was lower in the soleus than in the gastrocnemius (Fig. 3C). The highest expression levels of the Tbx1 gene was found in the quadriceps. Intermediate expression levels were measured in the soleus and the lowest expression level was found in the gastrocnemius. The Tbx1 expression levels between the three muscle types were all significantly different (Fig. 3D).
Gene expression patterns of Hoxd8 and Hoxd9 do not correlate with expression patterns of markers for the different fiber types
We found that markers for the oxidative fiber types as well as members of the Hoxd gene cluster were higher expressed in the gastrocnemius than in the quadriceps. To find out if these Hoxd genes play a role in specifying more oxidative muscle types we compared genes expression levels of Hoxd8, Hoxd9 and markers for different fiber types (type I: Myh7; type IIa: Myh2; type IIx: Myh1 and type IIb: Myh4) in quadriceps, gastrocnemius and soleus muscle of 20-week LFD mice and 20-week HFD mice. Expression levels of markers for the different fiber types were not affected by diet. Whereas Hoxd8 and Hoxd9 showed increased expression in the gastrocnemius as compared to quadriceps and soleus, expression levels of Myh7 was significantly higher in the soleus than in the quadriceps and gastrocnemius, respectively. Also the Myh2 expression level was higher in the soleus than in the quadriceps and gastrocnemius. For both Myh7 and Myh2 no differences between the quadriceps and the gastrocnemius were observed. Myh1 gene expression levels were comparable between the three muscle types. Finally, the expression level of the Myh4 gene was significantly lower in the soleus than in the quadriceps and gastrocnemius, respectively. Myh4 gene expression levels were similar in the quadriceps and gastrocnemius (Fig. 4). Altogether, the expression patterns of Hoxd8 and Hoxd9 did not correspond with the expression patterns of one of the markers for the different fiber types.
In the present study we searched for novel genes that are involved in specifying muscle types. Additionally, we studied the diet-sensitivity of differentially expressed genes. By comparing the gastrocnemius with the quadriceps we identified 162 differentially expressed genes corresponding to a variety of biological processes. Especially processes involved in cell cycle, contraction, development, differentiation, morphogenesis and remodeling were differentially regulated in the gastrocnemius vs. the quadriceps. Both the gastrocnemius and the quadriceps are muscle groups that consist of predominantly type II fibers. In the present study we found that Myh7, Tnnc1, Tnni1 and Tnnt1, which are all markers for the more oxidative fiber types, are strongly up-regulated in the gastrocnemius as compared to the quadriceps. Western blotting showed that protein levels of myosin heavy chain, slow fiber type protein (marker for the oxidative type I fibers) were similar between the gastrocnemius and the quadriceps, whereas the protein level of Myh2 (marker for the oxidoglycolytic type IIA fibers) was higher in the gastrocnemius than in the quadriceps (data not shown). Therefore, we propose that in this study the gastrocnemius had a more oxidative character than the quadriceps.
Recently, we showed that a 4-week HFD intervention results in increased gene expression and protein levels of markers for the oxidative fiber types in the quadriceps of mice. These findings suggested that the skeletal muscle can adapt to an increased lipid load by inducing a switch to a more oxidative phenotype . Although the adult skeletal muscle is known to have the capacity to adapt to functional demands including exercise, hormones, innervation and mechanical loading [15, 16], a link between HFD and a switch of fiber type was not yet described. However, in the present study we showed that an 8-week HFD has none or only minor effects on expression levels of genes that play a role in specifying different muscle types including the markers for oxidative fiber types. Therefore we speculate that the up-regulation of markers for the oxidative fiber types is an early response of the skeletal muscle to the HFD which is not maintained in the long run.
Among the strongest regulated genes we found five genes (Dkk3, Hoxd8, Hoxd9, Hoxd10 and Tbx1) which are especially known for their critical role during embryogenesis. The skeletal muscle is a dynamic tissue containing different types of stem cells such as satellite cells, side-population (SP) cells and muscle-derived stem cells (MDSC) [17, 18]. To exclude the possibility that expression of these five embryogenesis-related genes originated from other cell types than myocytes we analyzed the expression of genes that specify satellite cells, SP cells and MDSC (satellite cells, pax7, myf5, cd34 and c-met; SP cells, sca1, sdc4 and pax3; MDSC, cd34 and Blc2 [17, 18]). We found that only three genes (cd34, sca1 and sdc4) were detectably expressed (microarray signal intensity > 20) in the skeletal muscle of these mice. However, none of these detectably expressed markers differed more than 1.3-fold between the gastrocnemius and the quadriceps. Altogether, we cannot exclude the presence of cells other than myocytes, however, it is most unlikely that these cells have contributed to the different expression levels of Dkk3, Hoxd8, Hoxd9, Hoxd10 and Tbx1.
Two genes, Dkk3 and Tbx1, were expressed at a lower level in the gastrocnemius than in the quadriceps. Extrapolating our data to the soleus we found that Dkk3 expression in the soleus was lower than in the quadriceps or gastrocnemius. Interestingly, Tbx1 expression levels in the soleus were lower than in the quadriceps, but higher than in the gastrocnemius. During embryogenesis Dkk3 and Tbx 1 are expressed in a variety of organs including the heart and limb buds [19, 20]. Whereas the function of Dkk3 is still poorly understood , Tbx1 is well-studied especially in relation to heart development . Recently it was shown that Tbx1 is also involved in development of the limb buds by regulating the number of myocytes . Since further research revealed that Tbx1 is not necessary for skeletal muscle specification, differentiation, patterning or activation of the embryonic myogenic program, this is the only described function of Tbx1 in relation to embryonic muscle development . Both Dkk3 and Tbx1 are reported to be detectably expressed in the adult skeletal muscle [19, 24], but until now nothing is known about their possible function. We speculate that the relatively high expression levels of Dkk3 and Tbx1 in the quadriceps is involved in maintaining identity of this muscle type.
Three genes of the Hoxd gene cluster (Hoxd8, Hoxd9 and Hoxd10) were up-regulated in the gastrocnemius as compared with the quadriceps. The Hox gene family consists of at least 39 Hox genes organized in four gene clusters (A, B, C and D) which are subdivided into parallel groups numbered 1 to 13. The major function of Hox genes is the regulation of the formation of anterior-posterior patterning during embryonic development . Different groups have studied the expression levels of the Hox gene family in adult tissues of human origin. They found that each tissue displays a unique combination of detectable Hox gene expression levels which is altered when normal cells changes to malignant cancer cells [26–28]. Thus, in addition to their critical role in embryonic development, Hox genes play an important role in adult cells by controlling critical processes like cellular identity and differentiation . Houghton et al. . have studied expression patterns of the Hox gene family in skeletal muscle tissue of adult mice. In line with our observations they reported that the Hoxa10, Hoxc6, Hoxc9 and Hoxc10 genes were detectably expressed. Furthermore, they could not detect Hoxd11 which also is consistent with our findings. Other members of the Hoxd gene cluster were not studied . Thus to our knowledge, we are the first to show that Hoxd8, Hoxd9 and Hoxd10 are detectably expressed in adult murine muscle tissue, with highly increased expression levels in the gastrocnemius as compared with the quadriceps.
Together with the higher Hoxd8 and Hoxd9 gene expression levels in the gastrocnemius we found increased gene expression levels of markers for the more oxidative fiber types. Therefore we hypothesized that Hoxd8 and Hoxd9 might be involved in the determination of more oxidative muscle types. However, when we extrapolated our data to the soleus muscle we could not find a corresponding expression pattern between Hoxd8, Hoxd9 and Myh7 (marker for type I fibers) and Tnni1 (marker for type I fibers; data not shown). In fact, Myh7 and Tnni1 were highly expressed in the soleus, whereas Hoxd8 and Hoxd9 expression in the soleus was as low as in the quadriceps. Also the expression patterns of Myh2, Myh1 and Myh4, which are markers for type IIA, IIX and IIB fibers, respectively, did not show any resemblance with the expression pattern of Hoxd8 and Hoxd9. Although we did not find any resemblances between expression patterns of Hoxd8, Hoxd9 and any of the fiber type-specific markers we could distinguish the gastrocnemius, quadriceps and soleus from each other. Altogether, we suggest that the high expression levels of the Hoxd8, Hoxd9 and possibly Hoxd10 genes are a specific characteristic of the molecular profile of the gastrocnemius and are not related to an oxidative phenotype.
We found that especially genes that are markers for oxidative fiber types and genes known to be important for embryogenesis were differentially expressed in the gastrocnemius vs. the quadriceps. The expression levels of these genes were not influenced by diet. Furthermore, expression of Hoxd8 and Hoxd9 was low and expression of Dkk3 and Tbx1 was high in the quadriceps. In the gastrocnemius we observed high expression of Hoxd8 and Hoxd9 and low expression of Dkk3 and Tbx1. Finally, in the soleus we found low expression of Hoxd8, Hoxd9 and Dkk3, but intermediate expression of Tbx1. Comparing gastrocnemius, quadriceps and soleus, the gene expression pattern of Hoxd8 and Hoxd9 did not correspond to the gene expression pattern of any of the fiber type-specific markers. Therefore, we conclude that the highly expressed embryonic genes Dkk3, Hoxd8, Hoxd9 and Tbx1 are involved in identifying muscle types in a diet-independent and fiber type-unrelated manner.
Mice and diet
Male C57BL/6J mice were obtained from Harlan (Horst, The Netherlands). At 9 weeks of age mice were switched to a run-in diet consisting of an LFD (10 kcal% fat) for 3 weeks. Following the run-in period mice were randomly assigned to the LFD or the HFD (45 kcal% fat). Both diets contained fat in the form of palm oil (based on D12450B and D12451; Research Diet Services, Wijk bij Duurstede, the Netherlands) as described previously .
After 8 weeks of diet intervention, mice (n = 10 per diet) were fasted for 6 hours, anesthetized with a mixture of isofluorane (1.5%), nitrous oxide (70%) and oxygen (30%) and killed by cervical dislocation. Gastrocnemius and quadriceps muscles were dissected, snap-frozen in liquid nitrogen and stored at -80°C until further analysis.
After 20 weeks of diet intervention, mice (n = 5 per diet) were fasted for 6 hour, anaesthetized and killed as described. Quadriceps, gastrocnemius and soleus muscles were dissected, snap-frozen in liquid nitrogen and stored at -80°C until further analysis.
Both animal studies were approved by the Local Committee for Care and Use of Laboratory Animals at Wageningen University.
Mouse total RNA was isolated from gastrocnemius, quadriceps and soleus muscles by mechanically homogenization in Trizol reagent (Invitrogen, Breda, the Netherlands). RNA was purified with the RNeasy Mini Kit (Qiagen, Venlo, the Netherlands) and contaminating genomic DNA was removed with the RNase-free DNase set (Qiagen). RNA quantity and purity was measured with the ND-1000 spectrophotometer (Isogen Life Science B.V., IJsselstein, the Netherlands). An A260/A280 ratio between 2.04 and 2.16 was found for all RNA samples. RNA integrity was checked on an Agilent 2100 BioAnalyzer (Agilent Technologies, Amsterdam, The Netherlands) using nanochips according to the manufacturer's instructions. All RNA samples showed an RIN value between 7.5 and 9.0.
Affymetrix microarray analysis
RNA isolated from the gastrocnemius and quadriceps of mice fed an 8-week LFD or HFD was used for performing Affymetrix GeneChip® Mouse Genome 430 2.0 arrays as described previously . Array images were processed using packages from the Bioconductor project  and probe sets were redefined according to Dai et al. . In this method probes are annotated using up-to-date databases and assigned to unique gene identifiers, in this case Entrez ID's, instead of the 'classic' GeneChip probe sets. This results in a less ambiguous and more accurate annotation. Arrays were normalized using quantile normalization, and expression estimates were calculated using GC-RMA, implementing the empirical Bayes estimate for non-specific binding . Differentially expressed probe sets were identified using linear models, applying moderated t-statistics that implement empirical Bayes regularization of standard errors . Comparisons were made between the gastrocnemius and quadriceps muscles using the fold change threshold method. In this method, most false-positives are caused by genes with low signal intensity and/or absent call (~50% of all genes). Thus, filtering out this constant noise before performing data analysis will improve the power to discriminate true changes from noise in the fold-change threshold method . We manually set the signal threshold at 20, across all arrays, filtering out 60% of the genes. Probe sets that satisfied the criteria of a fold change > 1.3 and an FDR < 0.05, correcting for multiple-testing , were considered as differentially expressed.
To relate changes in gene expression to functional changes two complementary methods were used. The first method is based on overrepresentation of Gene Ontology (GO) terms, which uses a gene score resampling (GSR) method . Full resampling was run with 200,000 iterations. Only classes with a false discovery rate (FDR) < 0.001 and with minimal 8 and maximal 125 genes were taken into account. The second method, gene set enrichment analysis (GSEA), is focused on predefined gene sets, that is, groups of genes that share biological function, chromosomal location or regulation . The 'functional catalogue' constructed by Subramanian et al. was modified to contain only 505 well-defined murine, biochemical, metabolic and signal pathways compiled from the following publicly available, curated databases: BioCarta (BioCarta. BioCarta, 2005, p. http://www.biocarta.com.), GenMAPP , Kyoto Encyclopedia of Genes and Genomes (KEGG) , Sigma-Aldrich pathways (Sigma-Aldrich. Sigma-Aldrich Metabolic and Cell Signaling pathways, 2005.) and Signal Transduction Knowledge Environment (STKE; Signal Transduction Knowledge Environment. In: http://stke.sciencemag.org/. 2005.). The analysis was run using 1000 permutations per gene set. Gene sets with a FDR < 0.05 were considered as significantly regulated. The advantage is that both methods are unbiased and a score is computed based on all genes in a GO term or gene set.
The diet-sensitivity of differentially expressed genes, overrepresented GO classes and regulated gene sets was studied by comparing the gastrocnemius with the quadriceps under HFD conditions. Also in this comparison we used the criteria of a fold change > 1.3 and an FDR < 0.05. A heatmap of log-transformed microarray signal intensities was generated by using GeneMaths XT software (Sint-Martens-Latem, Belgium). Array data have been submitted to the Gene Expression Omnibus, accession number GSE18127.
Quantitative real-time PCR for RNA
To validate microarray-detected changes as well as extrapolating microarray-detected changes to the gastrocnemius, quadriceps and soleus of 20-week LFD mice and HFD mice quantitative real-time PCR (qPCR) was performed using individual cDNA samples. RNA (1 μg in 20 μl) was reverse transcribed using the iScript cDNA synthesis kit (Biorad, Veenendaal, the Netherlands) containing RNase H+ iScript reverse transcriptase, a premixed RNase inhibitor to prevent indiscriminate degradation of RNA template, and a unique blend of oligo(dT) and random primers. cDNA was synthesized using a 3-step program (5 minutes (min) at 25°C, 30 min at 42°C and 5 min at 85°C). Primer sequences were retrieved from the online PrimerBank database , or otherwise designed using Beacon Designer5 (Biorad). Primers were tested for specificity by BLAST analysis. The qPCR reactions were performed in a volume of 25 μl containing 12.5 ng cDNA, 1× IQ SYBR Green Supermix (100 mM KCl, 40 mM Tris-HCI, 6 mM MgCl2, 0.4 mM of each dNTP, 50 units/ml iTaq DNA polymerase, SYBR Green I and 20 nM fluoresein; Biorad) and 400 nM of gene-specific forward and reverse primers (Biolegio, Nijmegen, the Netherlands). cDNA was amplified using a 2-step program (40 cycles of 10 seconds (sec) at 95°C and 45 sec at 60°C) with a MyiQ system (Biorad). Specificity of amplification was verified by melt curve analysis and evaluation of efficiency of PCR amplification. PCR reactions were performed in duplicate and gene expression levels were determined using iQ5 software (Biorad) using a Δ-Cq relative quantification model with PCR efficiency correction and multiple reference gene normalization [43, 44]. As the Canx, Hprt1 and Arbp genes were identified as most stably expressed reference genes, we used the geometric mean of these three reference genes as normalization factor. Supplement 4 shows all relevant primer information.
All data are expressed as means ± SE. Microarray data were analyzed as described above. All other statistical analyses were performed using Prism software (GraphPad Software, San Diego, CA, USA). Two-way ANOVA was used to test gene expression levels of Dkk3, Hoxd8, Hoxd9, Myh1, Myh2, Myh4, Myh7 and Tbx1 between muscle types and diets. When significant differences were found, a Bonferroni post hoc test was used to determine the exact location of the differences. A p value ≤ 0.05 was considered as statistically significant.
Liang H, Ward WF: PGC-1alpha: a key regulator of energy metabolism. Advances in physiology education. 2006, 30 (4): 145-151. 10.1152/advan.00052.2006.
Scott W, Stevens J, Binder-Macleod SA: Human skeletal muscle fiber type classifications. Physical therapy. 2001, 81 (11): 1810-1816.
Berchtold MW, Brinkmeier H, Muntener M: Calcium ion in skeletal muscle: its crucial role for muscle function, plasticity, and disease. Physiological reviews. 2000, 80 (3): 1215-1265.
Booth FW, Thomason DB: Molecular and cellular adaptation of muscle in response to exercise: perspectives of various models. Physiological reviews. 1991, 71 (2): 541-585.
Chanseaume E, Malpuech-Brugere C, Patrac V, Bielicki G, Rousset P, Couturier K, Salles J, Renou JP, Boirie Y, Morio B: Diets high in sugar, fat, and energy induce muscle type-specific adaptations in mitochondrial functions in rats. The Journal of nutrition. 2006, 136 (8): 2194-2200.
Wang YX, Zhang CL, Yu RT, Cho HK, Nelson MC, Bayuga-Ocampo CR, Ham J, Kang H, Evans RM: Regulation of muscle fiber type and running endurance by PPARdelta. PLoS biology. 2004, 2 (10): e294-10.1371/journal.pbio.0020294.
Kho AT, Kang PB, Kohane IS, Kunkel LM: Transcriptome-scale similarities between mouse and human skeletal muscles with normal and myopathic phenotypes. BMC musculoskeletal disorders. 2006, 7: 23-10.1186/1471-2474-7-23.
Naya FJ, Mercer B, Shelton J, Richardson JA, Williams RS, Olson EN: Stimulation of slow skeletal muscle fiber gene expression by calcineurin in vivo. The Journal of biological chemistry. 2000, 275 (7): 4545-4548. 10.1074/jbc.275.7.4545.
Talmadge RJ, Otis JS, Rittler MR, Garcia ND, Spencer SR, Lees SJ, Naya FJ: Calcineurin activation influences muscle phenotype in a muscle-specific fashion. BMC cell biology. 2004, 5: 28-10.1186/1471-2121-5-28.
Lin J, Wu H, Tarr PT, Zhang CY, Wu Z, Boss O, Michael LF, Puigserver P, Isotani E, Olson EN: Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature. 2002, 418 (6899): 797-801. 10.1038/nature00904.
Mortensen OH, Frandsen L, Schjerling P, Nishimura E, Grunnet N: PGC-1alpha and PGC-1beta have both similar and distinct effects on myofiber switching toward an oxidative phenotype. Am J Physiol Endocrinol Metab. 2006, 291 (4): E807-816. 10.1152/ajpendo.00591.2005.
Luquet S, Lopez-Soriano J, Holst D, Fredenrich A, Melki J, Rassoulzadegan M, Grimaldi PA: Peroxisome proliferator-activated receptor delta controls muscle development and oxidative capability. Faseb J. 2003, 17 (15): 2299-2301.
Hakimi P, Yang J, Casadesus G, Massillon D, Tolentino-Silva F, Nye CK, Cabrera ME, Hagen DR, Utter CB, Baghdy Y: Overexpression of the Cytosolic Form of Phosphoenolpyruvate Carboxykinase (GTP) in Skeletal Muscle Repatterns Energy Metabolism in the Mouse. Journal of Biological Chemistry. 2007, 282 (45): 32844-32855. 10.1074/jbc.M706127200.
de Wilde J, Mohren R, Berg van den S, Boekschoten M, Dijk KW, de Groot P, Muller M, Mariman E, Smit E: Short-term high fat-feeding results in morphological and metabolic adaptations in the skeletal muscle of C57BL/6J mice. Physiol Genomics. 2008, 32 (3): 360-369.
Pette D, Staron RS: Myosin isoforms, muscle fiber types, and transitions. Microscopy research and technique. 2000, 50 (6): 500-509. 10.1002/1097-0029(20000915)50:6<500::AID-JEMT7>3.0.CO;2-7.
Pette D, Staron RS: Transitions of muscle fiber phenotypic profiles. Histochemistry and cell biology. 2001, 115 (5): 359-372.
Otto A, Collins-Hooper H, Patel K: The origin, molecular regulation and therapeutic potential of myogenic stem cell populations. Journal of anatomy. 2009, 215 (5): 477-497. 10.1111/j.1469-7580.2009.01138.x.
Figeac N, Daczewska M, Marcelle C, Jagla K: Muscle stem cells and model systems for their investigation. Development dynamics. 2007, 236 (12): 3332-3342. 10.1002/dvdy.21345.
Monaghan AP, Kioschis P, Wu W, Zuniga A, Bock D, Poustka A, Delius H, Niehrs C: Dickkopf genes are co-ordinately expressed in mesodermal lineages. Mech Dev. 1999, 87 (1-2): 45-56. 10.1016/S0925-4773(99)00138-0.
Stennard FA, Harvey RP: T-box transcription factors and their roles in regulatory hierarchies in the developing heart. Development (Cambridge, England). 2005, 132 (22): 4897-4910.
Niehrs C: Function and biological roles of the Dickkopf family of Wnt modulators. Oncogene. 2006, 25 (57): 7469-7481. 10.1038/sj.onc.1210054.
Dastjerdi A, Robson L, Walker R, Hadley J, Zhang Z, Rodriguez-Niedenfuhr M, Ataliotis P, Baldini A, Scambler P, Francis-West P: Tbx1 regulation of myogenic differentiation in the limb and cranial mesoderm. Dev Dyn. 2007, 236 (2): 353-363. 10.1002/dvdy.21010.
Grifone R, Jarry T, Dandonneau M, Grenier J, Duprez D, Kelly RG: Properties of branchiomeric and somite-derived muscle development in Tbx1 mutant embryos. Dev Dyn. 2008, 237 (10): 3071-3078. 10.1002/dvdy.21718.
Chieffo C, Garvey N, Gong W, Roe B, Zhang G, Silver L, Emanuel BS, Budarf ML: Isolation and characterization of a gene from the DiGeorge chromosomal region homologous to the mouse Tbx1 gene. Genomics. 1997, 43 (3): 267-277. 10.1006/geno.1997.4829.
Svingen T, Tonissen KF: Hox transcription factors and their elusive mammalian gene targets. Heredity. 2006, 97 (2): 88-96. 10.1038/sj.hdy.6800847.
Calvo R, West J, Franklin W, Erickson P, Bemis L, Li E, Helfrich B, Bunn P, Roche J, Brambilla E: Altered HOX and WNT7A expression in human lung cancer. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (23): 12776-12781. 10.1073/pnas.97.23.12776.
Cantile M, Pettinato G, Procino A, Feliciello I, Cindolo L, Cillo C: In vivo expression of the whole HOX gene network in human breast cancer. Eur J Cancer. 2003, 39 (2): 257-264. 10.1016/S0959-8049(02)00599-3.
Cillo C: HOX genes in human cancers. Invasion & metastasis. 1994, 14 (1-6): 38-49.
Cillo C, Cantile M, Faiella A, Boncinelli E: Homeobox genes in normal and malignant cells. Journal of cellular physiology. 2001, 188 (2): 161-169. 10.1002/jcp.1115.
Houghton L, Rosenthal N: Regulation of a muscle-specific transgene by persistent expression of Hox genes in postnatal murine limb muscle. Dev Dyn. 1999, 216 (4-5): 385-397. 10.1002/(SICI)1097-0177(199912)216:4/5<385::AID-DVDY7>3.0.CO;2-G.
de Wit NJ, Bosch-Vermeulen H, de Groot PJ, Hooiveld GJ, Bromhaar MM, Jansen J, Muller M, Meer van der R: The role of the small intestine in the development of dietary fat-induced obesity and insulin resistance in C57BL/6J mice. BMC medical genomics. 2008, 1: 14-10.1186/1755-8794-1-14.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome biology. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33 (20): e175-10.1093/nar/gni179.
Wu Z, Irizarry R, Gentleman R, Martinez-Murillo F, Spencer F: A model-based background adjustment for oligonucleotide expression arrays. Journal of the American Statistical Association. 2004, 99 (468): 909-917. 10.1198/016214504000000683.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3
Xinmin L, Kim J, Zhou J, Gu W, Quigg R: Use of signal thresholds to determine significant changes in microarray data analyses. Genetics and Molecular Biology. 2005, 28: 191-200. 10.1590/S1415-47572005000200002.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.
Lee HK, Braynen W, Keshav K, Pavlidis P: ErmineJ: tool for functional analysis of gene expression data sets. BMC Bioinformatics. 2005, 6: 269-10.1186/1471-2105-6-269.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR: GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways. Nat Genet. 2002, 31 (1): 19-20. 10.1038/ng0502-19.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.
Wang X, Seed B: A PCR primer bank for quantitative gene expression analysis. 2003, 31: e154-
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome biology. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome biology. 2007, 8 (2): R19-10.1186/gb-2007-8-2-r19.
We thank Freek Bouwman, Mechteld Grootte Bromhaar, Shohreh Keshtkar Ghiasabadi and Bert Weijers for excellent technical assistance. This study was funded by the Top Institute Food and Nutrition, with financial support by the Dutch government.
JdW, ES and EM participated in the design and supervision of the study. JdW and MH performed the RNA isolations and qPCRs. PDG and MB assessed the quality control of the microarrays and provided support for microarray analysis. JDW performed microarray analysis and drafted the manuscript. EM and ES provided valuable feedback on the initial draft. All authors read and approved the final manuscript
Electronic supplementary material
Additional file 1: List of differentially expressed genes. Additional file 1 gives the differentially expressed genes in the gastrocnemius as compared to the quadriceps from 10 LFD mice. To find these differentially expressed genes we used the criteria of an average signal intensity > 20 over all arrays, a fold change > 1.3 and an FDR < 0.05. (XLS 55 KB)
Additional file 2: Overview of qPCR analyses for 13 genes validating the microarray results. Additional file 2 gives an overview of the microarray signal intensity, relative expression (qPCR), normalized expression (qPCR) and corresponding fold changes of 13 genes. (XLS 22 KB)
Additional file 3: Differentially expressed genes, overrepresented GO classes (ErmineJ) and regulated gene sets (GSEA) between gastrocnemius and quadriceps in either LFD or HFD conditions. Additional file 3 gives an overview of the differentially expressed genes, overrepresented GO classes (ErmineJ) and regulated gene sets (GSEA) in the gastrocnemius as compared to the quadriceps in HFD and LFD conditions. (XLS 114 KB)
About this article
Cite this article
de Wilde, J., Hulshof, M.F., Boekschoten, M.V. et al. The embryonic genes Dkk3, Hoxd8, Hoxd9 and Tbx1 identify muscle types in a diet-independent and fiber-type unrelated way. BMC Genomics 11, 176 (2010). https://doi.org/10.1186/1471-2164-11-176