1 Commissariat à lEnergie Atomique, Service de Génomique Fonctionnelle, and 2 Genethon, Centre National de la Recherche Scientifique UMR 8115, Evry, France
ABSTRACT
With the use of Hoechst staining techniques, we have previously shown that the C2C12 myogenic cell line contains a side population (SP) that is largely increased in the presence of fibroblast growth factor 6 (FGF6). Here, we compared transcriptional profiles from SP and main population (MP) cells from either C2C12 or FGF6-expressing C2C12. Expression profiles of SPs show that these cells are less differentiated than MPs and display some similarities to stem cells. Moreover, principal component analysis made it possible to distinguish specific contributions of either FGF6 or differentiation effects on gene expression profiles. This demonstrated that FGF6-expanded SPs were similar to parental C2C12-derived SPs. Conversely, FGF6-treated MPs differed from parental MPs and were more related to SP cells. These results show that FGF6 pushed committed myogenic cells toward a more immature phenotype resulting in the accumulation of cells with a SP phenotype. We propose that FGF6 conditioning could provide a way to expand the pool of immature cells by myoblast dedifferentiation.
muscle; fibroblast growth factor 6; stem cell; microarray; principal component analysis
STEM CELLS have been identified in many adult tissues that undergo extensive cell replacement due to physiological turnover or injury. Until recently, these tissue-specific stem cells (SCs) have been considered monopotential, meaning they can give rise to cells contributing to homeostasis of the parental tissue, with little or no transdifferentiation occurring naturally. This accepted opinion has been reconsidered in the light of recent studies showing that hematopoietic SCs could participate in angiogenesis and muscle and hepatic regeneration as well as in the formation of central nervous system cell types (1, 13, 14, 25, 33, 44, 45). Conversely, it has been reported that blood cells could derive from SCs originating from different tissues (2, 26). Therefore, the bulk of information from the literature has suggested that SC populations display an "apparent" plasticity or versatility giving them the capability to differentiate into a range of distinct mature cells probably depending on the dominant effect of the local environment (for reviews, see Refs. 35 and 40). However, this hypothesis is challenged by recent findings showing that SCs used in these paradigms consisted of a mixture of SC from various origins (27, 32). In addition, it has been shown that these cells could also fuse with committed cells, thus indicating that transdifferentiation was not necessarily required to explain the plasticity of the SCs studied (49).
The rarity of these cells and the absence of specific markers have made the search for pure SC populations a challenge during the past years. Most of the studies have focused on the characterization of SCs using various well-known hematopoietic cell surface markers. The main limit of this approach is that purified cells consist of a mixture of cells from different lineages, which may share common surface markers. A way to purify the most immature cells from an adult tissue without prejudging their immunophenotypic characteristics is to use their properties to exclude vital dyes due to the fact that they highly express multidrug/ATP-binding cassette transporters (7, 52). As a result, SCs are only slightly stained in the presence of Hoechst 33342 and were collected using fluorescence-activated cell sorting (FACS) as a side population (SP) consisting of about 0.10.5% of total cells. Conversely, the more mature cells, forming the main population (MP), displayed low Hoechst excluding properties and thus were brightly stained. This technique was originally used to sort hematopoietic SCs with high repopulating activity when they were injected into lethally irradiated animals (17, 18, 46).
SP cells have been sorted from many tissues including skeletal muscle (20). However, immunophenotyping of the cells accumulated in these tissue-specific SP fractions showed a number of variations as well as a certain heterogeneity (3). Because skeletal muscle contains several types of cells (i.e., muscle, connective tissue, blood vessels, nerves, etc.), isolation of the most immature myogenic cells cannot be achieved by this method alone because the muscle-derived SP fraction would potentially consist of a mixture of putative SCs with varied differentiation properties. As an example, it has been shown that SP cells isolated from muscle biopsies may be either positive or negative for CD45 expression. All muscle-derived hematopoietic progenitor and in vivo bone marrow reconstitution activity is derived from CD45-positive cells (3). In addition, muscle-derived CD45-positive cells purified from uninjured muscle are uniformly nonmyogenic in vitro and do not form muscle efficiently in vivo compared with CD45-negative cells (3). More recently, it has been observed that CD45-positive cells purified from regenerating muscle readily gave rise to myoblasts (39). These results underline the high complexity of the biology of the muscular SC and the difficulty to isolate cell populations using solely the expression of known expressed membrane markers.
A way to gain insight the molecular characterization of true muscle SP cells is to use a pure source of myogenic cells such as the C2C12 cell line. We established that the C2C12 cell line contained a subset of SP cells that consisted in quiescent cells sharing common features with regular SCs (5). Like muscle-derived SP cells and/or cell populations obtained by serial preplating or from disaggregated primary muscle, C2C12-derived SP cells exhibited reduced adhesion to plastic relative to fully committed myoblasts and required a longer time to recover myotube-forming capacities, features that are currently accepted as being those of early muscle precursor cells. C2C12-derived SP cells are CD34 and lymphocyte antigen 6 complex, locus A (Sca1)-positive and express ATP-binding cassette, subfamily B (MDR/TAP), member 1a (mdr1a). These cells are also myoblast determination protein (MyoD) low/ and resemble "dormant" cells because they are in the G0/G1 phase. In a recent study, Kondo et al. (29) also demonstrated that SP cells are present in a cell line, the malignant C6 glioma cell line, suggesting that the presence of a SP population in the myogenic C2C12 cell line is not a unique feature.
We used this model to investigate the possibility of expanding this fraction by using appropriate tissue culture conditions. Specifically, we (23) showed in a previous study that fibroblast growth factor 6 (FGF6) increases the C2C12-derived SP fraction by around 50 times. To ensure that such an increase was not due to insufficient dye staining, SP analysis was performed using a range of Hoechst concentrations. This work also showed that FGF6 selectively upregulates the mdr1a gene (but not the mdr1b gene) and suggested a role for FGF6 in the maintenance of a reserve pool of progenitor cells in skeletal muscle.
Here, our goal was to initiate the characterization of the genetic program in SP and MP cells sorted from either C2C12 or FGF6-expressing C2C12 (C2CF6). This would first give us a standard baseline for primary muscle SCs and also make it possible to determine whether SP cells obtained in the presence of FGF6 are equivalent to canonical SPs.
Recent publications (24, 41) show that microarray technology is an important tool for the characterization of SCs. Previously, a large number of cells was required to perform such analysis, essentially excluding the use of the microarray technology for cells that are too infrequent in vivo, as it is the case for SCs. Here, we used the in vitro transcription RNA amplification procedure first described by Van Gelder et al. (48), which allowed us to perform global expression profiling of genes in SP and MP populations from C2C12 cells transduced or not by FGF6 using <100,000 cells in each sorted fraction.
With the use of microarray analysis, we characterized the transcriptome of both SP and MP cells purified from C2C12 (S2 and M2 cells) or C2CF6 (S6 and M6 cells). The differential analysis of expression patterns allowed us to dissociate the effects of FGF6 from that of differentiation and to define the transcriptional signature of these SP cells.
MATERIALS AND METHODS
Cell Cultures and RNA Preparation
C2C12, a subclone of the C2 mouse myoblast cell line (6, 50), was obtained from the American Type Culture Collection. The C2C12 cell line retrovirally transduced by a murine leukemia virus-derived vector recombinant for a bis-cistronic construction consisting of a murine FGF6 and green fluorescent protein (GFP) cDNA separated by internal ribosome entry site (IRES) as previously described (23) was named C2CF6. Cell cultures, Hoechst staining, and FACS analysis were performed as previously described (23). About 30,00050,000 SP cells were obtained from each sorting starting with a total of 10 x 106 cells. Total RNA was extracted from the FACS-sorted cells using the Absolutely RNA Microprep Kit as described in the manufacturers protocol (Stratagene). One microgram of total RNA from each sample was amplified using the MessageAmp Amplified RNA (aRNA) Kit (Ambion) according to the manufacturers protocol. The integrity of total RNA and aRNA was verified using a Bioanalyzer (Agilent).
Production of DNA Microarrays
DNA microarrays were generated using a collection of 7,200 probes. A total of 1,700 mouse cDNA clones (as bacterial colonies) was obtained from Research Genetics. Approximately 1,400 genes were represented in this set of clones (according to Unigene cluster comparison). All these clones are from Integrated Molecular Analysis of Genome and Their Expression libraries and have been sequence verified. The cDNA clone inserts were amplified using two primers complementary to flanking sequences in the cloning vector [M13(21) and M13 REV]. A total of 2,200 PCR products was obtained using specific oligonucleotides based on a selection of genes chosen as responsive to keywords like stem cell, apoptosis, growth factors, transcription factors, etc. One thousand five hundred mouse cDNA clones were produced from a subtracted library (mdx muscle/normal muscle) and 1,800 rat cDNA clones from another subtracted library (atrophied muscle/normal muscle). The 7,200 PCR products were analyzed by electrophoresis on an agarose gel for quality control and quantitation and then ethanol precipitated. Control probes were added at least in duplicate, including positive controls (GAPDH, actin, and tubulin) and negative controls (genes of plant, plasmid DNA, and TE/DMSO). Finally, a total of 7,690 probes was spotted on the slide. Microarrays were prepared by printing PCR-amplified probes, arrayed in 384-well microtiter plates, suspended in a spotting buffer composed of 50% DMSO and TE on AminoSilane slides (CMT GAPS II, Corning) with a robot Microgrid II (BioRobotics) at a density of 1,000 cells/cm2 under a controlled atmosphere.
Target Preparation and Array Hybridization
For each target preparation, 1 µg of aRNA was reverse transcribed using SuperscriptTM II reverse transcriptase (Life Technologies) in the presence of random hexamers and an amino-modified nucleotide (amino allyl dUTP). Amino-modified cDNAs were purified through a Microcon Centricon 30 microconcentrator (Amicon) and ethanol precipitated. In a second step, monofunctional forms of Cy3 and Cy5 dyes (Amersham) were coupled with the purified amino-modified cDNAs. Unincorporated fluorescent molecules and salts were removed through the Microcon Centricon 30. Labeled cDNA was mixed with 10 µg poly(A) RNA (Boehringer), 10 µg tRNA (Life Technologies), and 10 µg mouse Cot1 DNA (Life Technologies). Each slide was rehydrated over boiling water and quickly heated on a hot plate, and the probes were ultraviolet light cross-linked in a Stratalinker at 254 nm/250 mJ. The slides were prehybridized in 5x SSC-0.1% SDS-1% BSA for 30 min at 50°C. Purified labeled targets were resuspended in 40 µl of hybridization solution (70% formamide-3.6x Dehnardt-0.7% SDS-8.6x saline-sodium phosphate-EDTA) and heated at 95°C for denaturation. The labeled targets were applied to prehybridized slides and covered with a 60 x 22-mm polyethylene coverslip (Sigma). The hybridization was performed at 42°C overnight. Washings were performed at room temperature for 15 min in 0.1x SSC-0.1% SDS and for 2 x 15 min in 0.1x SCC.
Data Processing
After hybridization and washing, fluorescent signals were acquired by scanning each slide using the ScanArray 5000 scanner (Packard). A separate image was captured for each of the two fluorophores used. The resolution of the scan was 10 µm/pixel on a gray scale. The .tiff images resulting from the scan were imported into image-analysis GenePix Pro 4.0 software (Axon) to quantitate the signal for each spot. Spots with obvious blemishes were flagged. Genes with null intensity values and low mean intensity were also annotated with a specific flag. All the array elements (flagged or not flagged) were included in the statistical analysis, but flagged elements were excluded from the list of validated differentially expressed genes.
Statistical Analysis of Experimental Data
ANOVA.
Standard ANOVA was performed on the logarithm of raw gene expression data without prior background subtraction. Four main factors were included in the model as well as four two-order interaction factors, as proposed in Ref. 28. The model equation is given by Eq. 1:
![]() | (1) |
![]() | (2) |
Once these corrections had been performed, the corresponding effects were suppressed from the ANOVA model. The possibly remaining factors that could not be corrected were then accounted for in the model residues. The resulting model was then of the following form [Eq. 3 (43)]:
![]() | (3) |
Selection of differentially expressed genes.
Differentially expressed genes were then identified using a classical parametric hypothesis test of mean comparison. A differential score g was computed for each gene as the difference between the mean intensity values in the two biological conditions that are compared (Eq. 4):
![]() | (4) |
Let X be a normal random variable with mean 0 and standard deviation and a be a real value so that the probability that |X| is greater than a, P(|X| > a) =
. Then, gene g is considered as differentially expressed if
![]() |
Global analysis with principal component analysis.
Principal component analysis (PCA) is a descriptive multivariate data analysis method. It enables exploration of the structure of large data sets in which data points are described by a large number of variables that prevents visualizing it. PCA takes advantage of the correlation that exists between the variables. It provides a projection of the data set in a new space of reduced size that allows visualization in a two-dimensional space while minimizing the loss of information. The axes of this new space are called principal components. They are defined as linear combinations of the original variables. Two graphical representations result from a PCA: the "data set" representation and the "variables" one.
In the first type of graph, the data points are plotted in the new space. A restriction to the first two principal components makes it possible to visualize the structure of the data set as well the proximity between data points. Thus data points that appear to be close in this two-dimensional space may be similar with respect to the original variables that were chosen to describe them, given that the projected data set is close to the original one. Outliers are data points whose values for the original variables are far from the mean values exhibited by the whole data set.
The second type of graph shows the correlations between the original variables and principal components. It points out the variables that influence at most the structure of the data set. Both representations provide a complementary view of the data set. They may be superimposed on a single graph.
Most computations for the data analysis were carried out using R software (22) except for global analyses, which were performed using either Matlab software (Mathworks; Natick, MA) or SAS software (SAS Institute; Cary, NC). The microarray dataset used for this analysis is described and available from GEO (GSE1436 and GSM241086GSM24142).
Real-time quantitative RT-PCR analysis.
RNA samples prepared from SP and MP cells purified from C2C12 (S2 and M2 cells) or C2CF6 (S6 and M6 cells) were analyzed by RT-PCR using the One-Step RT-PCR Kit (Qiagen) with 10 ng RNA. The incorporation of SYBR green dye into PCR products was monitored in real time with an ABI PRISM 7700 sequence detection system (PE Applied Biosystems). The efficiency of the amplification was determined for each pair of primers by comparison with a standard curve generated with serially diluted cDNA. Target genes were quantified relative to a reference gene (18S) using the mathematical model described by Pfaffl (37). All PCRs were performed in triplicate. Primers were as follows: 1-actin (Acta1), forward 5'-GTGAGATTGTGCGCGACATC-3' and reverse 5'-GGCAACGGAAACGCTCATT-3'; IGF binding protein 2 (Igfbp2), forward 5'-CCTCAAGTCAGGCATGAAGGA-3' and reverse 5'-GCAGGGAGTAGAGATGTTCCA-3'; ß-actin (Actb), forward 5'-GAAATCGTGCGTGACATCAAAG-3' and reverse 5'-TGTAGTTTCATGGATGCCACAG-3'; glycoprotein 38 (GP38), forward 5'-AGAGAACACGAGAGTACAACCA-3' and reverse 5'-TGCGTTTCATCCCCTGCATT-3'; secretory leukocyte protease inhibitor (Slpi), forward 5'-CTGTTCCCATTCGCAAACCAG-3' and reverse 5'-CCACATATACCCTCACAGCACTT-3'; activated leukocyte cell adhesion molecule (Alcam), forward 5'-ACGAAGAAAAGTGTGCAGTATGA-3' and reverse 5'-ACTAGGGTAGGTGCTTCAAACA-3'; p21, forward 5'-GAAAACGGAGG-CAGACCAGC-3' and reverse 5'-CACAGCAGAAGAGGGCGGG-3'; CD34, forward 5'-CAGCAGTAAGACCACACCAGC-3' and reverse 5'-GGGGAAGTCTGTGGTTGTGAA-3'; inhibitor of DNA binding 2 (Idb2), forward 5'-AGAACCAGGCGTCCAGGAC-3' and reverse 5-CTGCAAGGACAGGATGCTGA-3'; Sca1, forward 5'-ACCTATGCTGGTGGTCTGCC-3' and reverse 5'-GACCAGAGCCTCTGGGTTGA-3'; transformed mouse 3T3 cell double minute 2 (MDM2), forward 5'-GTCTACCGAGGGTGCTGCAA-3' and reverse 5'-TCCCCAGGTAGCTCATCTGTG-3'; ATP-binding cassette, subfamily B, member 4 (Abcb4), forward 5'-CTCAACACACGTCTAACAGATGA-3' and reverse 5'-GATGAACCCCACTATGAATCCTG-3'; and ATP-binding cassette, subfamily A, member 8 (ABCA8), forward 5'-ACCCCAACAACTCAGAGGATAA-3' and reverse 5'-CCACTACATCACTGAAACGCAT-3'.
RESULTS
Isolation of C2C12 and C2CF6 SP and MP Cell Populations by Hoechst 33342 DNA Staining
Hoechst 33342 staining was used to isolate SP and MP cell populations from the C2C12 murine muscle cell line expressing FGF6 or not as previously described (24) (Fig. 1). We were able to discriminate SP (S2) and MP (M2) fractions from C2C12, giving us around 0.5% of SP cells (0.01% with verapamil) (Fig. 1A). With the use of the C2CF6 cell line, we observed the same discrimination of SP (S6) and MP (M6) cell fractions with a significant increase of the SP cell fraction in the presence of FGF6 to up to 20% in average (n = 6) compared with the untransduced cells (Fig. 1B).
|
The results showed important array, dye, biological condition, array x dye, and dye x gene effects for all biological samples that were compared. Because it is performed independently on each array, the lowess normalization method theoretically reduces the array x dye x gene effect, which does not appear explicitly in the ANOVA model used. In practice, it reduces both the dye x gene and array x dye effects (data not shown). It is used to diminish the well-known red-green imbalance in low intensities.
The second normalization process, based on the subtraction of significant main effects, can be considered as robust, because the mean values used for the correction are estimated over a large set of intensity values.
Accuracy and Robustness of the Experimental Data
Figure 2A shows the accuracy of hybridization performed with the same cDNA target labeled with either Cy3 or Cy5. This plot represents the differential scores obtained for all the probes involved in the experiment, sorted in the ascending order. The region of acceptance of H0 (no differential expression) is between the upper and lower 1, 5, 10, 15, or 20% type I error rate thresholds, according to the false positive rate that we are ready to accept. At the level of the 20% error rate, 99% of the differential scores fall within these thresholds, meaning that there are 2 of 7,200 (0.03%) false positives. The same conclusion can be drawn from the experiment using hybridization of the same sample after RNA amplification (Fig. 2B) with 10 of 7,200 (0.07%) false positives at the level of the 20% error rate. The threshold value to determine genes differentially expressed is inversely proportional to the square root of the number of dye swaps involved in the experiment. Thus, as the number of replicates increases, the decision threshold decreases, resulting in a larger number of genes that are found to be significantly differentially expressed.
|
|
Thus genes with up- and downregulated expression were identified for all hybridization conditions, and the distribution of the number of genes does not indicate a significant global imbalance of transcriptionally modulated genes between each comparison (data not shown).
The functional analysis of these genes, presented in Table 1, indicates that the majority of them are involved in biological processes such as RNA synthesis (9%), metabolism (27%), cell shape (16%), cell signaling (17%), cell division (14%), and cell defense (4%). Interestingly, the genes involved in the shape of the cell (matrix, cytoskeletal, and muscular markers) were preferentially downregulated in the M6, S2, and S6 cell fractions compared with the M2 cell fraction. This suggests a more undifferentiated phenotype of these cells. In addition, no significant variation of the expression of cell shape-related genes was observed between the S2, S6, and M6 fractions, suggesting a similar undifferentiated phenotype of these cells in relation to these molecular markers.
|
In the first case, the complete set of validated genes was used (5,654 genes). In the second case, the PCA table was restricted to genes that had been found differentially expressed in at least one of the six experiments (297 genes). In both cases, the resulting cluster was centered on a core of genes with little (or no) differential expression, whereas the most modulated genes behaved as "outliers" and were located outside, as shown in Fig. 4. Superimposed on this data cluster were the six directions associated with the six experiments (or variables). A simultaneous analysis of the data set and variables pointed out the genes that best characterized either each experiment or a set of experiments. Here, outliers define three distinct directions. The colinear axis (Fig. 4) to the direction M2/M6 and S2/S6 distinguishes genes that characterized the FGF6 effect independently of their phenotypic link (SP or MP). At one extremity are the most upregulated genes, and at the other one are the most downregulated ones. The second axis is derived from the S2/M2 and S6/M2 directions defined by genes with modulated expression between the M2 and S2/S6 cell fractions. It corresponds to a "differentiation" effect from differentiated M2 cells to less-differentiated S2/S6 cells.
|
The S2/S6 direction almost follows the FGF6 axis as defined above (Fig. 4). Thus this indicates that there are no significant differences between these two cell fraction transcriptomes in relation to the 7,690 molecular markers present in our study, showing that these cells have a similar SP phenotype.
To confirm the descriptive results obtained by PCA, we performed ANOVA and a differential analysis on the two subsets of hybridizations. The genes up- or downregulated in FGF6-treated cells compared with untreated cells were labeled as genes preferentially modulated by FGF6 or implicated in the FGF6 effect; the genes upregulated or downregulated in SP compared with MP cells were labeled as genes preferentially implicated in the differentiation effect. Thus 511 genes were highlighted with a P value <0.20 according to a FGF6 effect and 322 genes according to a differentiation effect. The results presented in Tables 2 and 3 show the 250 genes, with annotated function, with expression regulated in response to the presence of FGF6. Tables 4 and 5 present the 164 annotated genes preferentially implicated in the differentiation effect.
|
|
|
|
The distribution of these genes, accordingly to seven main cellular functions (RNA synthesis, protein synthesis, metabolism, cell shape, cell signaling, cell division, and cell defense) is presented in Fig. 5A and shows the modulation of four principal functions: metabolism (24%), cell shape (15%), cell signaling (22%), and cell division (16%).
|
Microarray experimental data validation by real-time RT-PCR.
We performed real-time RT-PCR for 14 genes to confirm both the FGF6 and differentiation effects (Table 6). To validate the FGF6 effect, we confirmed the downregulation of Acta1, Idb2, Igfbp2, and Actb in M6 versus M2 cell fractions and the upregulation of ABCA8, Slpi, Abcg4, GP38, Sca1, Alcam, MDM2, CD34, and the control gene FGF6 in M6 versus M2 cell fractions. Moreover, for Myd116, we observed weak upregulation in M6 versus M2 cell fractions and strong upregulation in S6 versus M2 cell fractions (not shown) and observed an upregulation of Sca1 expression in M6 versus M2 cell fractions not detected by the microarray.
|
Real-time RT-PCR was not suitable to confirm the expression modulations observed in S6 versus M2 cell fractions, which combines both the FGF6 and differentiation effects.
DISCUSSION
In this study, we compared gene expression profiles from cells with SP and MP phenotypes derived from the C2C12 muscle cell line in the presence or absence of FGF6.
Many factors have been shown to act upon the process of satellite cell activation and recruitment, such as IGF-1, leukemia inhibitory factor, IL-6, and FGF6. In particular, it has been shown that knockout of FGF6 in normal mice leads to severely defective muscle regeneration and, in the mdx mouse, to a considerable aggravation of the dystrophic phenotype with diminished expression of MyoD at sites of necrosis-regeneration and a concomitant accumulation of interstitial collagen (16). FGF6 is specifically produced by muscle fibers and both promotes cell proliferation and counteracts cell apoptosis.
On the other hand, FGF6 also gives rise to multiple effects depending on its concentration. At 5 ng/ml, it promotes differentiation of myogenic precursors (38). At 25 ng/ml, FGF6 acts as a muscle proliferation agent and delays differentiation of muscle precursors (38). At about 100 ng/ml, we have shown that this growth factor reduces cell growth while increasing the proportion of cells with the SP phenotype (23). These unique features designate FGF6 as a key factor in muscle regeneration. Therefore, we decided to perform global transcriptional profiling to examine FGF6 for its capacity to either increase or mobilize the pool of immature muscle cells. In the presence of high concentrations of FGF6 (up to 75100 ng/ml), the proportion of SP cells (S6) was tremendously increased. As we (23) have previously described, C2CF6 differs from the original C2C12 cell line: C2CF6 cells are more likely to divide and have increased resistance to apoptosis. Moreover, they never reach a confluent state and do not fuse to form myotubes. They also display severe downregulation of expression of some important myogenic regulatory factors such as MyoD and myogenin.
With the use of microarray technology coupled with robust statistical analysis of the data, we identified genes involved in the FGF6 cell response. Among them, well-known actors of the FGF family or serum-responsive genes such as methylenetetrahydrofolate dehydrogenase (mthFd2) (36), progressive ankylosis (ank) (19), RAB5C (Rac3) (21), and immediate-early response 3 (Ier3) (8) have been identified. Interestingly, we demonstrated modulation of FGF receptors 1 (upregulation) and 4 (downregulation) and downregulation of Igfbp2. These results corroborate the proliferative phenotype of the C2C12 cells and their dedifferentiation in the presence of FGF6, which has been previously reported (38), and show that Igfbp2 may be implicated in the differentiation process in C2C12 cells in response to FGF6. Thus the product of Igfbp2 is a well-known inhibitor of IGF-1. Some reports of the role of the IGF growth factor suggest its implication in the regulation of skeletal muscle adaptation in response to different phenotypic events such as satellite cell proliferation and differentiation (42, 47), inhibition of apoptosis (30, 34), and improvement of the overall functional properties of skeletal muscles in dystrophic mice (4, 31). This suggests a potential role of FGF6 in the IGF pathway through the modulation of the expression of Igfbp2, and we hypothesize that this regulation may be preferentially involved in subpopulations of muscular SP cells. In addition, we observed a drift of the molecular program of the C2C12 cells that corroborates important modifications of C2C12 behavior in the presence of FGF6. These modifications concern the extracellular matrix and cytoskeleton structure with a downregulation of the transcription of major actors in their organization, such as procollagen, titin, and actin. The drift of the molecular program in C2C12 cells ends in a dedifferentiated phenotype that contributes to the increase of the SP cell population. This is confirmed by the upregulation of two well-known SC/SP cell markers, Sca1 and CD34, and two molecular markers of dedifferentiation, MDM2 and Alcam. Indeed, we found that MDM2 is one of the most overexpressed genes in the presence of FGF. The protein encoded by this gene is known to inhibit MyoD function and muscular differentiation (15) and to induce G1/S arrest (for a review, see Ref. 11). On the other hand, ALCAM (or CD166) is involved in homophilic adhesion as well as binding to CD6. Recently, the human homolog of ALCAM, hematopoietic cell antigen (HCA), was isolated by Cortes et al. (10), and this protein was detected in the most primitive subset of hematopoietic SCs as well as in myeloid progenitors in bone marrow. Here, we support this conjecture by demonstrating the upregulation of ALCAM in early progenitor muscular cells. These findings imply that ALCAM might be a key adhesion molecule involved in the characterization of SCs.
Finally, it appears that C2C12-derived SP cells probably correspond to a more immature cell population compared with the MP cells, thus confirming our previous reports indicating that the C2C12 cell line displays spontaneous cell heterogeneity, although it was initially clonal (5). We show that FGF6 increases the SP fraction in the C2C12 cell line and that M6 cells were different from parental MP (M2) cells, whereas S6 cells were similar to the original SP (S2). Thus a key question remains: whether the resulting S6 population derived from the accumulation of dedifferentiated MP cells in the presence of FGF6 or from direct mobilization of the resident SP subset. In previous reports, we showed that the SP derived from C2C12 cells was predominantly in the G0/G1 phase (5) and that the presence of FGF6 did not modify this state (23). These results suggested that the increase in the SP fraction in the C2CF6 cells was due to dedifferentiation and not to SP cell proliferation. Dye efflux in SP cells is mediated by membrane-associated P-glycoproteins of the ATP-binding cassette transporter family that act as energy-dependent pumps. Previous reports had shown that these proteins were associated with an increase in the population of cells with an immature phenotype (5, 9, 41, 52). Here, we identified two members of the ATP-binding cassette transporter family upregulated in presence of FGF6. It is predictable that many ATP-binding cassette transporters might show comparable effects provided their biological activities resulted in fashioning a similar metabolic environment by pumping out substrates and factors involved in sustaining cell differentiation. Here, we add support to this notion by showing that the increase of C2CF6-derived SP cells was due to the dedifferentiation of MP into M6 and S6 cells. The molecular mechanism involved in this process is still unknown, and we propose that the pumps implicated in the dye efflux in SP cells could be involved in this phenomenon. In conclusion, we hypothesize, as described in Fig. 6, that FGF6 triggers dedifferentiation of the C2C12 MP into a dedifferentiated population MP' (M6), which is phenotypically different from the original cell population and has the capacity to "give rise" to SP cells (S6) similar to native C2C12 SP cells (S2). This finding suggests that the SC or immature phenotype of the SP subset of cells could be acquired by dedifferentiation of a main population of cells in adult tissues under environmental influence.
|
GRANTS
The research carried out in the laboratory of the authors was supported by grants from the Association Française contre les Myopathies, the Commissariat à lEnergie Atomique, and the Centre National de la Recherche Scientifique.
ACKNOWLEDGMENTS
The authors are grateful to Dr. Claude Dechesne (UMR6543-Centre National de la Recherche Scientifique, Nice, France) for providing the mouse and rat subtracted cDNA libraries. We also thank Terry Partridge, Vincent Frouin, and Carlo Lucchesi for many helpful discussions and Sophie Lemoine for technical support.
FOOTNOTES
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
Address for reprint requests and other correspondence: G. Piétu, CECS/I-STEM, Batiment Généthon, 1 rue de lInternationale, 91000 Evry, France (e-mail: gpietu{at}istem.genethon.fr).
10.1152/physiolgenomics.00141.2004
REFERENCES
|
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Visit Other APS Journals Online |