Previous Article | Next Article ![]()
Molecular and Cellular Biology, March 2009, p. 1538-1553, Vol. 29, No. 6
0270-7306/09/$08.00+0 doi:10.1128/MCB.01375-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
and
Rob Hooft van Huijsduijnen*
Merck Serono SA, 1211 Geneva, Switzerland
Received 29 August 2008/ Returned for modification 11 December 2008/ Accepted 22 December 2008
|
|
|---|
|
|
|---|
In order to better understand which genes play a role in multiple sclerosis, a number of proteomic (13), genomic (25, 29), and genetic (35) approaches have been utilized. However, diseased lesions consist of different admixtures (e.g., contain immune and glia infiltrates) of cell types compared to controls, so such postmortem samples may present differential cell, rather than gene or protein, expression. Therefore, understanding how such genes fit into pathophysiological processes (and in which) is problematic. From population genetic surveys thus far, only a few genes have been found to be associated with multiple sclerosis (31).
Based on these considerations, we chose to examine genome-wide gene expression changes in a murine oligodendroglial precursor cell line, Oli-neu (18), as these cells underwent differentiation into myelin basic protein (MBP)-producing cells. Even for a pure cell line, this type of transition is typically associated with changes in the expression of thousands of genes, by itself providing little meaningful information. We therefore decided to look at oligodendrocyte differentiation as induced by different agents, in the hope that from this combined data set one might extract "core genes" whose modulation is closely linked to the differentiation process. The enriched set of genes was further evaluated for their functional involvement in the differentiation process.
|
|
|---|
Differentiation medium was composed of DMEM-F-12 supplemented with 50 µg of apo-transferrin/ml, 5 ng of sodium selenite/ml, 100 µM putrescine, 130 nM progesterone, 200 nM T3, 1% horse serum, 1% penicillin-streptomycin, 2 mM glutamine, 0.1% sodium bicarbonate, and 15 mM HEPES. Chemical treatments for microarray analysis included insulin at 10 µg/ml, 10 µM forskolin (Sigma F6886), 1 µM 13-cis-retinoic acid (LKT Laboratories R1779), 1 µM dexamethasone (Sigma D4902), and 1 µM PD174256 (Alexis ALX-270-323).
OPC isolation. Cortical oligodendrocyte precursors (OPCs) were purified from newborn Oncins France strain A rats. Isolated cortex was roughly dissociated with a scalpel and then incubated 10 min at 37°C in 10 ml of 1x Hanks balanced salt solution (Invitrogen catalog no. 14025)-0.01% trypsin (Invitrogen catalog no. 25300)-0.0024% DNase (Sigma D-5025). After centrifugation at 800 rpm for 5 min, the pellet was resuspended by pipetting in 2 ml of DMEM (Invitrogen catalog no. 31966)-1% penicillin-streptomycin-20% fetal bovine serum (Invitrogen catalog no. 10108), filtered through a 70-µm-pore-size cell strainer (BD-352350), and centrifuged again at 800 rpm for 10 min. One cortex per 75-cm2 poly-D-lysine flask (BD-356537) was put in culture in DMEM supplemented with 20% fetal calf serum. The medium was changed every 3 days. At day 10, flasks were shaken for 1 h at 150 rpm, and the medium (containing the microglia) was removed. Then, 10 ml of fresh medium was added, and the flasks were shaken overnight at 280 rpm.
After filtration through a 150-µm-pore-size Nitex filter, the OPCs were counted and plated at the required density on poly-D-lysine in SATO serum-free medium composed of DMEM and 60 µg of N-acetylcysteine (Sigma A-9165), 5 µg of insulin (Sigma I-0516), 10 ng of biotin (Sigma B-4639), 20 ng of forskolin (Sigma F-6886), 100 µg of bovine serum albumin (Sigma A-4161), 100 µg of apo-transferrin (Sigma T-1147), 60 ng of progesterone (Sigma P-8783), 16 µg of putrescine (Sigma P-5780), and 40 ng of NaSeO3 (Sigma S-9133)/ml.
RNA interference (RNAi) experiments. Annealed small interfering RNA (siRNA) duplexes (1 nM) were introduced into cells by using INTERFERin transfection reagent (catalog no. 409-10; Polyplus Transfection, Inc.). Transfection complexes were prepared in OptiMEM medium (Invitrogen catalog no. 51985) for 15 min at room temperature according to the manufacturer's instructions. The complexes were added to the cell suspension, and the cells were grown for 72 h.
RNA extraction and qPCR analysis. For the microarrays, total RNA extraction and DNase digestion were performed according to the manufacturer's instructions (RNeasy minikit; Qiagen). Nucleic acid concentrations of the samples were determined by using a NanoDrop spectrophotometer. Concentration and quality of the samples were confirmed by Bioanalyzer 2100 (Agilent Technologies). For all other quantitative (real-time) PCR (qPCR) assays, total RNA was prepared with TRIzol reagent (Invitrogen catalog no. 15596) according to the manufacturer's protocol. A total of 1 µg of total RNA was used as a template for reverse transcription with Superscript III (Invitrogen catalog no. 18080) and a 1:1 mix of oligo(dT)20 (Invitrogen catalog no. 18418) and random primers (Invitrogen catalog no. 48190). cDNA was then diluted 50-fold, and real-time PCR was performed in triplicates using the QuantiFast SYBR green PCR master mix (Qiagen catalog no. 204054) and the respective QuantiTect primer assays.
Microarray procedure. A one-cycle eukaryotic target labeling assay (one-cycle target labeling and control reagents; Affymetrix catalog no. P/N 900493) was used to produce fragmented biotinylated cRNA from total RNA as described in the standard Affymetrix protocol (P/N 701021 Rev.5).
In this study, 3 µg of total RNA isolated were first reverse transcribed using a T7-oligo(dT) promoter primer in the first-strand cDNA synthesis reaction. After RNase H-mediated second-strand cDNA synthesis, the double-stranded cDNA was purified and served as a template in the subsequent in vitro transcription reaction. The in vitro transcription reaction was carried out in the presence of T7 RNA polymerase and a biotinylated nucleotide analog-ribonucleotide mix for cRNA amplification and biotin labeling. The biotinylated cRNA targets were then cleaned up, an aliquot of 20 µg was fragmented, and finally 15 µg was hybridized to a GeneChip Mouse Genome 430 2.0 array (Affymetrix, P/N 900495, 900496, and 900497). Chips were washed and labeled in a fluidics station, using protocol EuKGE-WS2v_45 After scanning, CEL files were analyzed by using GCOS software to check poly(A) and hybridization controls.
Data analysis. Analysis of the microarray experiment was carried using Resolver software from Rosetta. Differentially expressed genes were analyzed comparing each treatment and time point versus the respective time-matched untreated samples. Quality control measures were performed by using a two-dimensional agglomerative Resolver hierarchical clustering analysis (HCA) procedure using z-score transformed intensity experiments with a cosine correlation similarity measure, and Ward's minimum variance heuristic criteria were applied on two probe sets selections for untreated, dexamethasone, retinoic acid, and PD174265 conditions at 10 h and 72 h.
GCOS and MAS5.0 software from Affymetrix were in the range recommended by the manufacturer. The intensity profile pipeline from Rosetta was used to create an intensity (expression) profile for each sample. Then, the intensity experiment pipeline from Rosetta was used to derive, from a triplicate of intensity profiles for a given condition, a single expression value for each probe set. Following that step the ratio experiment pipeline from Resolver was used to obtain n-fold changes between treatment versus time-matched untreated intensity experiments. A probe set with intensity P values given by the Affymetrix error model from Resolver below 10–6 was considered present. n-Fold changes in P significance values for probe sets were derived from Resolver using a two-sample t test with no false discovery rate correction. Moreover, a principal component analysis (PCA) was performed using the intensities profiles obtained for the untreated, dexamethasone, retinoic acid, and PD174265 samples at each time point (10, 24, and 72 h). Only probe sets called present in at least three arrays and with a variation coefficient greater than 0.3 were kept in the PCA step.
A two-dimensional agglomerative Resolver HCA procedure using z-score-transformed intensity experiments with a cosine correlation similarity measure and the Ward's minimum variance heuristic criteria was applied on two probe sets selections for untreated, dexamethasone, retinoic acid, and PD174265 conditions at 10 and 72 h.
The first probe set selection used in the HCA corresponds to the genes listed in Tables 1 and 2.
|
View this table: [in a new window] |
TABLE 1. "Core" differentiation genes selected for functional analysisa
|
|
View this table: [in a new window] |
TABLE 2. Additional genes selected for knockdown analysis (see the text)a
|
Western blot. Sodium dodecyl sulfate-polyacrylamide gels (Invitrogen NuPAGE Novex Bis-Tris gels,12%) were run according to the manufacturer's instructions, and proteins were transferred to nitrocellulose membrane (Invitrogen LC2001) using a semidry transfer apparatus (Hoefer SemiPhor; Amersham Pharmacia Biotech). The membrane was blocked in washing buffer (phosphate-buffered saline [PBS], 0.2% Tween 20) plus 5% dissolved nonfat milk powder. It was incubated overnight at 4°C with anti-CNP (M-300) Santa Cruz sc-30158 rabbit polyclonal antibody or anti-MBP (C-16) Santa Cruz sc-13914 goat polyclonal antibody at a 1:1,000 dilution. The blot was washed, incubated for 1 h at room temperature with a corresponding horseradish peroxidase (HRP)-conjugated antiserum (goat anti-rabbit Ig-HRP [Bio-Rad 170-6515] or donkey anti-goat Ig-HRP [Santa Cruz sc-2020] at 1:2,000 dilution), washed again, and visualized by chemiluminescence (ECL kit; Amersham Pharmacia Biotech RPN 2106).
CNP enzymatic assay.
siRNA-transfected cells in 96-well plate were lysed, after 3 days (
80% confluence), in 80 µl of triple-detergent buffer (50 mM Tris [pH 8], 150 mM NaCl, 0.1% sodium dodecyl sulfate, 1% NP-40, 0.5% sodium deoxycholate) with complete protease inhibitor cocktail (Roche catalog no. 11836170001) and then incubated for 15 min at room temperature. The protein content was measured using a BCA protein assay kit (Pierce catalog no. 23225). Then, 2.5 µg of total protein was plated in a black Costar 96-well plate, and the cyclic nucleotide phosphodiesterase (CNP) reaction was run in 100 µl of 200 mM MES free acid hydrate (Sigma M8250)-30 mM MgCl2-1.11 mM 2'-3'-cNADP (Sigma N5257)-5.55 mM glucose-6-phosphate (Sigma G7879)-0.06 U of glucose-6-phosphate-dehydrogenase (Sigma G8164) for 15 min at room temperature. The fluorescence intensity (FI) was read on a Perkin-Elmer Victor 2 spectrofluorimeter (excitation at 355 nm, emission at 460 nm, for 0.1 s) and normalized for protein content.
Cellomics readout using O4 staining. siRNA-transfected cells in 96-well plates were fixed with 4% paraformaldehyde 15 min at room temperature, washed twice with 1x PBS, and incubated with 50 µl per well of mouse anti-O4 monoclonal antibody diluted 20 times in 1x PBS at 4°C overnight. Cells were washed twice with 1x PBS and incubated with 50 µl per well of anti-mouse immunoglobulin M-Alexa Fluor 488 (Invitrogen A21042) diluted 200 times and incubated for 1 h at room temperature. Nuclei were stained with Hoechst dye diluted 10,000-fold for a further 15 min. Cells were washed twice with 1x PBS before fluorescent cell images were obtained on an ArrayScan HCS reader with a x5 or x10 objective lens, and effective image analysis was done by Cellomics Technologies (Thermo Scientific) for determination of the "neurite outgrowth" (NOG) index.
MBP enzyme-linked immunosorbent assay (ELISA). siRNA-transfected cells in 96-well plate were fixed with 1% glutaraldehyde (Fluka catalog no. 49631) in 1x PBS for 15 min room temperature, washed twice with 1x PBS, and incubated 1 h at room temperature in blocking buffer (PBS, 0.2% Tween 20, 5% nonfat milk powder). Cells were incubated with anti-MBP (C-16; Santa-Cruz sc-13914) goat polyclonal antibody at a 1:400 dilution in 5% milk-TBS (Tris-borate-EDTA)-0.1% Tween 20 overnight 4°C. After two washes in 1x PBS-0.2% Tween, cells were incubated with donkey anti-goat Ig-HRP (Santa Cruz sc-2020) antibody at a 1:400 dilution in 5% milk-TBS-0.1% Tween 20 for 1 h at room temperature. Cells were washed twice in 1x PBS-0.2% Tween and once in 1x PBS before we assessed the HRP activity with its chromogenic substrate ABTS [2,2'azinobis(3-ethylbenzthiazolinesulfonic acid); Sigma A1888] and obtained a reading at 405 nm.
Phosphatase assays. Unless otherwise stated, chemicals were purchased from Sigma-Fluka (Saint-Louis, MO). DiFMUP (6,8-difluoro-4-methyumbelliferyl phosphate) was obtained from Molecular Probes/Invitrogen (Carlsbad, CA). N-terminally fluorescein-labeled (Fl) phosphopeptide [EAIY(PO3)AAPFA] was purchased from Caliper Life Sciences (Hopkinton, MA). The 96- and 384-well plates were from Corning Inc. (Corning, NY). The glutathione S-transferase-tagged catalytic domain of the human phosphatases MKP5, PAC1, PTP1B, PTPβ, PTP-H1, PTPµ, SHP1, SHP2, TC-PTP, VHR, and Glepp1 were cloned, expressed in Escherichia coli, and purified by affinity analysis on a glutathione column as described previously (32, 41).
Mobility shift assay for MKP5 inhibitors. Assays were performed in 384-well plates. A total of 5 µl of fluorescein-labeled peptide substrate at 3 µM in reaction buffer (20 mM Bis-Tris-HCl [pH 6.2], 0.01% Igepal, 5 mM dimethyl sulfoxide [DMSO]), corresponding to a 1 µM peptide final concentration, was added to 5-µl portions of compounds in 5% DMSO. The reaction was started by the addition of 5 µl of the enzyme at a concentration of 0.66 µg/ml (corresponding to a final concentration of 0.22 µg/ml) in reaction buffer containing 3 mM DL-dithiothreitol. After 2 h of incubation, 15 µl of termination buffer (100 mM HEPES [pH 7.5], 0.01% Igepal, 0.1% coating reagent 3, 5% DMSO) was added. Aliquots from each well were sipped by a 12 sippers chip model 760137-0372R on a HTS 250 drug discovery system instrument (Caliper Life Sciences) for electrophoretic separation of the substrate and product peaks and measurement of the FI of the peaks (1). The relative peak heights of the substrate and product were measured, and ratios were calculated using HTS Well Analyzer Software 4.1 from the manufacturer. The ratio (r) was defined as the height of the product peak divided by the sum of the peak height of the product and substrate. The low controls corresponded to the absence of enzyme in columns 1 and 2 (A-H) and 24, and the high controls corresponded to the absence of compounds in columns 1 and 2 (I-P). The percent inhibition was calculated relative to the high and low controls as 100 x [1 – (r – rlow control/rhigh control – rlow control)].
DiFMUP-based assay for phosphatase selectivity. For selectivity phosphatase assays were carried out using a DiFMUP concentration corresponding to the Km of the enzyme studied (2). In 96-well plates containing 5 µl of diluted compound or solvent (100% DMSO) in each well, 55 µl of DiFMUP diluted in PTP buffer (20 mM Bis-Tris-HCl [pH 7.5], 0.1% Brij 35, 1 mM DL-dithiothreitol) was added, followed by 40 µl of recombinant enzyme diluted in PTPB buffer in order to start the reaction. After 45 min at room temperature, the FI was measured (excitation at 355 nm and emission at 460 nm for 0.2 s) on a Perkin-Elmer Fusion spectrofluorimeter (Perkin-Elmer Life Sciences). Negative controls were performed in the absence of enzyme, and positive controls were carried out in the presence of enzyme without compound. The percentage of activity was calculated according to the following formula: % activity = 100 x (FIcompound – FIlow control)/(FIhigh control – FIlow control). The 50% inhibitory doses (IC50) were determined in triplicate in two independent experiments.
Gene Expression Omnibus database accession number. The expression data described here have been deposited in the NCBI Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo) under accession number GSE14406.
|
|
|---|
![]() View larger version (35K): [in a new window] |
FIG. 1. Experimental setup of gene expression analysis. (Left) Oli-neu cells were induced by using four different chemical treatments (forskolin, retinoic acid, dexamethasone, and PD174265) compared to untreated or insulin-treated controls. Differentiation results in variable synthesis of MBP and dendrocyte outgrowth, measured using O4-stained immunofluorescence (for PD174265; shown at the bottom). (Right) Chip hybridization protocol. Each whole-genome Affymetrix chip was hybridized to one of three biological triplicate of cell cultures from six treatments over four time points (61 hybridizations; t = 0 shared).
|
![]() View larger version (49K): [in a new window] |
FIG. 2. (A) Intensity profile PCA for untreated, retinoic acid, dexamethasone, and PD174265 induced gene expression data. The numbers indicate time points (in hours). (B) HCA for 4,260 probe set expression values at 10 and 72 h for untreated, dexamethasone, retinoic acid, and PD174265 conditions. Five clusters were defined; clusters 3 (1,695 probe sets) and 5 (875 probe sets) were retained for further analysis (listed and analyzed in in the supplemental material).
|
Overall, these analyses indicated that different treatments affected different gene sets, even though they induced a similar phenotype. We counted regulated genes common for multiple treatments for various time points using Venn diagrams (Fig. 3A to E; see Materials and Methods for data analysis). Genes in overlaps were counted, while including forskolin (Fig. 3A and D) or insulin (control; Fig. 3C to E) or including only a single steroid treatment (Fig. 3C and D). Counting only genes associated with differentiation (red intersections, listed numbers) showed a clear, time-dependent increase for all diagrams. This trend was not seen for the totals of modulated genes for each treatment (Fig. 3F), thus demonstrating true convergence of the phenotype with a specific, complex transcriptional signature. We conclude that the oligodendrocyte differentiation process is characterized by an expanding set of commonly modulated genes, while individual treatments produce additional, likely unrelated transcriptional "noise." However, we cannot exclude that this "noise" represents in fact alternative, inducer-specific paths toward the differentiated state. We reasoned that the shared gene sets (Venn diagram intersections) might be enriched for key players in the differentiation process.
![]() View larger version (34K): [in a new window] |
FIG. 3. Genes specifically regulated by various differentiating agents. (A to D) Number of genes (probesets) significantly up- or downregulated (>2-fold, compared to untreated cells at the same time point) common to various combination of treatments. The numbers below each Venn diagram refer to the red-colored overlapping section and are given for different time points. (E) Overlap and criteria for the set of 48 genes selected for functional analysis (see also Table 1 and the text). (F) Total numbers of genes modulated (>2-fold) per treatment at different time points.
|
2-fold by PD174256 and >2-fold and >1.5-fold by dexamethasone and retinoic acid, respectively. The reason for looking at the early 10-h time point is that this is the "process-forming" differentiation stage induced by all treatments (Fig. 1); thus, this time point reflects true phenotypic convergence. The 48 genes in this overlap are listed in Table 1, with expression changes indicated for dexamethasone and PD174265 for the three time points. Functional involvement of genes whose activity is modulated through differentiation regardless of induction mode. Having identified a set of "signature" genes whose expression changes specifically accompany Oli-neu differentiation as induced by various agents, we next sought to understand which of these genes are functionally involved in the process, using siRNA-mediated gene silencing. In addition to the genes listed in Table 1, we also included two other sets of genes that were regulated during differentiation but did not fulfill all inclusion criteria. Such genes were included (Table 2, upper panel) because they also have chromosomal single-nucleotide polymorphisms that were statistically associated with multiple sclerosis (H. Abderrahim, unpublished data). We also included additional, regulated genes (Table 2, lower panel) from mitogen-activated protein kinase (MAPK) signaling pathways and associated dual-specificity phosphatases (DUSPs) and oligodendrocyte-specific genes as reported in the literature.
For functional interrogation of the gene set listed in Tables 1 and 2, we considered two experimental approaches. The first sought to determine which genes are essential for differentiation. This was attempted by treating Oli-neus with siRNAs, followed by treatment with the chemical inducers. Unfortunately, we failed to identify such gene candidates with confidence; since only a percentage of cells is RNAi transfected and yet all respond to the inducing chemical, our positive controls failed (I. Foucault, unpublished data). We then looked for genes whose silencing might drive spontaneous differentiation. The readout in this protocol (see below) is quite sensitive, and this approach did produce a number of potential positives.
For each of the 88 selected genes (Tables 1 and 2), two siRNAs were tested, targeting exons found in each (or in a majority of) transcripts for a given gene (all siRNAs and their sequences are given in the supplemental material). In the initial step of our screening procedure, Oli-neus were transfected in duplicate (for each siRNA) in 96-well plates, and after 72 h the lysates were tested for CNP activity. A surprisingly large subset of 47 siRNAs showed activity, and these were retested in six-well plates using CNP, MBP, and PLP qPCR analysis as readouts (data not shown). Finally, for 16 genes that were selected from this set, both siRNAs (for each gene) were tested in a 12- or 24-well format. Differentiation was measured by CNP enzymatic activity (Fig. 4A); MBP ELISA (Fig. 4B); qPCR analysis for PLP (Fig. 4C), MBP (Fig. 4D), or CNP mRNA (Fig. 4E); and in a 96-well format with Cellomics readout using O4 monoclonal antibody staining for process formation (Fig. 4F). Induction of differentiation was compared to luciferase siRNA (negative control). As a positive control, we used siRNA for ErbB2 (Neu), the target of the PD174265 inhibitor in these cells. As shown in Fig. 4, a number of siRNAs induced differentiation, as judged by mRNA and protein expression of markers. In contrast, Dscaml1 (named for Down syndrome cell adhesion molecule-like 1) produced inconsistent results (with only one siRNA being active and the other showing opposite behavior). A notable and unexpected "positive" is CNP/CNP1; while its siRNAs reduced expression of CNP protein and mRNA (Fig. 4A and E, leftmost bars); this silencing was accompanied by Oli-neu differentiation induction, as seen for all other tested markers (Fig. 4B, C, D, and F). Induction of differentiation markers was accompanied by changes in microscopic morphology (process formation), as shown for selected genes in Fig. 4G. The lower panel in this figure represents active siRNAs; the top panel shows negative siRNAs.
![]() View larger version (37K): [in a new window] |
FIG. 4. Induction of oligodendrocyte differentiation markers after siRNA-mediated silencing of genes in Oli-neu (3 days). Two different siRNAs were used for each gene (blue and pink bars). The positive control was ErbB2/Neu (green bars); the negative control (for ratio calculation) was luciferase siRNA. The data were combined from multiple experiments as indicated (e.g., two separate experiments using duplicates for panel A). The readouts were CNP enzymatic activity (A), MBP ELISA (B), PLP qPCR (C), MBP qPCR (D), CNP qPCR (E), and dendrocyte outgrowth (F) measured using O4-antibody staining (Cellomics). (G) Phase-contrast imaging of morphological changes accompanying Oli-neu differentiation as induced by siRNAs for indicated genes. Luciferase and Ptdss2 (top row) are inactive siRNAs.
|
![]() View larger version (29K): [in a new window] |
FIG. 5. RNAi-induced of differentiation of rat OPCs. The readouts were MBP and PLP mRNA induction as measured by qPCR (A) and MBP and PLP protein as measured by Western blotting (autoradiogram shown in lower panel) (B). mRNA knockdown achieved by the siRNAs (measured by qPCR) is indicated at the bottom of panel A. Statistical significance (versus the luciferase siRNA negative control): ***, P < 0.001; **, P < 0.05; *, P < 0.1.
|
) are all known to exert negative control on this cascade, while Nlgn3 (neuroligin 3) and Trpc4 (transient receptor potential cation channel 4) may enter the system via Src and Dlg4/PSD-95, the latter being a central nervous system-restricted guanylate kinase.
![]() View larger version (74K): [in a new window] |
FIG. 6. Pathway mapping of OPC differentiation inducer genes. Genes for siRNAs that induced differentiation in Oli-neu and confirmed in OPCs (see the list) were groupwise analyzed in the context of established gene-gene interactions using Ingenuity Pathway analysis software. Lines correspond to direct (physical) interactions (enzyme-substrate, yeast two-hybrid, or immunoprecipitation; dashed lines represent indirect interactions, e.g., induction).
|
![]() View larger version (28K): [in a new window] |
FIG. 7. Induction of OPC differentiation using an enzymatic MKP5 inhibitor. (A) Potent induction of MBP and CNP by siRNA knockdown on MKP5 in OPCs. (B) Selectivity among protein tyrosine phosphatases of an organic, novel MKP5 inhibitor, AS077234-4 (IC50 given in µM). (C) Treatment of rat OPCs with AS077234-4 results in dose-dependent induction of MBP and CNP differentiation marker proteins. DMSO, negative control with compound solvent only. (D) Western blot used for the graphic representation (scan) of panel C.
|
10% of non-OPCs in the preparation used in the Dugas study provide negative-feedback signaling to the OPCs as they differentiate, slowing down and segregating differentiation into distinct processes. Combined with our results, this suggests that OPCs are subject to negative control in mixed cultures and probably in situ as well. Known examples for such inhibitors of differentiation include the Lingo1 receptor and its ligands NogoA, Mag (myelin-associated glycoprotein), and Omgp (oligodendrocyte-myelin glycoprotein [recently reviewed in reference 28]).
![]() View larger version (24K): [in a new window] |
FIG. 8. Induction by PD174265 in Oli-neu of a subset of genes earlier identified (9) as "early" and "late" myelin genes in rat OPCs. "Early" genes are colored in panel A; "late" genes are colored in panel B. The suffixes "-a," "-b," and "-c" refer to different probe sets for the same gene. Induction for the Mag probe was off-scale, possibly due to off-target hybridization. The induction of all genes is highly significant (P < 0.01).
|
|
|
|---|
for Th17, and FoxP3 for Treg cells in the immune system, while embryonic stem cells can be reprogrammed using defined sets of just a few transcription factors (reviewed in references 23 and 42). The oligodendrocyte lineage is under the control of Olig1/2; Sox (SRY-box containing gene) -5, -6, and -10; Oct6; and Nkx2.2 (NK2 homeobox 2) transcription factors (26, 42). How these transcription factors interact with environmental stimuli and cascade a myriad of transcriptional activities is poorly understood. The differentiation process has been viewed as one generating stable attractors from an uncommitted metastable state to another that involves the creation of autostimulation and cross-inhibition (15). Such a view is compatible with the generation of multiple redundant, negative controls seen in our Oli-neu study. Our work focused on identifying relevant downstream effectors of OPC differentiation, starting out from comprehensive sets of genes that were regulated during Oli-neu differentiation as induced after treatment with various agents. We found that each differentiating treatment induced changes in expression in a characteristic set of genes and that only a small subset of these genes was shared between treatments. As differentiation progressed over time, however, the set of commonly modulated genes expanded.
The same data set can also be analyzed for genes that specifically associate with myelination, as opposed to process formation, which we examined earlier. We identified genes that were modulated by PD174265 (a treatment that resulted in MBP synthesis) but were not, or differently regulated by forskolin (which induced process, but not MBP formation). Analysis of the resulting gene set did not match to a single, well-characterized network, presumably because in this case we did not have a sufficiently large number of different treatments that would allow us to filter the gene sets.
Our results demonstrate that pure Oli-neu cells present an excellent model for oligodendrocyte differentiation, as witnessed by the induction of various myelin-associated markers and dendritic outgrowth. For instance, Tmem10/Opalin was recently identified as a novel oligodendrocyte differentiation marker (19), and this gene's expression increased 43-fold in PD174265-treated Oli cells. Likewise, known critical transcription factors Sox8 and -17 (reviewed in reference 42) were up- and downregulated, respectively, in our datasets, as described previously. Furthermore, we demonstrate that differentiation is under considerable negative control both in Oli-neu cells and in primary OPCs. An earlier example of oligodendrocyte negative control is Sirt2 (sirtuin 2), an oligodendroglial cell-specific gene whose repression results in CNP and MBP induction (24). It was proposed that Sirt2 functions "to prevent overdifferentiation or early aging of the [oligodendroglial] cells" (24). In our data set, Sirt2 gene expression is strongly induced by PD174265.
CNP, a well-studied constituent of uncompacted myelin, had not been identified earlier as (yet another) myelin component that inhibits oligodendrocyte differentiation. Mice that overexpress CNP lack myelin compaction and undergo precocious maturation (12, 43). Surprisingly, mice that lack CNP present normal myelin and yet undergo severe neurodegeneration due to "uncoupled oligodendrocyte functions" (21). Our work showing that CNP negatively regulates oligodendrocyte differentiation provides a basis for understanding these in vivo observations and reveals a role for CNP in appropriate timing and control of OPC differentiation.
The seven genes that we identified as autoinhibitory (listed in Fig. 6, along with MBP) are all expressed in various areas of the murine central nervous system (see the expression data in the supplemental material as retrieved from the Internet Allen Brain Atlas, which is described in reference 22)—with the exception of PTPRR, which is nevertheless known to play important neural functions (8). Among murine neural cell lines, PTPRJ, PTPRR, Dusp10, and Nlgn3 are all overexpressed in oligodendrocytes, whereas Tspan2 (tetraspanin 2) displays 28-fold increased expression in myelinating oligodendrocytes compared to OPCs (see reference 5 and online material associated with this reference).
Tspan2 has been discovered in an oligodendrocyte/OPC subtractive cDNA library (4) and was later postulated to represent "a molecular link between surface integrins and a CD9, Tspan-2 molecular web during the differentiation of oligodendrocytes" (39). Our results support models in which this complex is involved in cell-cell interactions that control oligodendrocyte differentiation.
Three of the other "repressors" that we identified are more usually associated with neurons: TRPC4 is a potential cation channel expressed in the developing (murine) brain (44), and Dscaml1 may, based on its expression pattern, be involved in formation and maintenance of neural networks (1). Nlgn3 is a neuron-specific adhesion molecule, mutations in which are associated with autism-spectrum disorders and mental retardation (7), diseases which have also been associated with abnormal myelination (10).
Three repressors (PTPRj, PTPRr, and MKP5) are protein tyrosine phosphatases, whose members are known to exert negative feedback on "forward" kinase phosphorylation cascades. For instance, PTP1B dampens insulin receptor signaling (reviewed in reference 14), while mice that lack PTPH1 display enhanced growth hormone sensitivity and systemic growth (34). PTPRj, also named density-enhanced phosphatase (DEP), has an antiproliferative effect (3), may use Erk1/2 as a substrate (27), and has been linked to cell adhesion (33). MKP5 inactivates p38 and JNK in vitro (38, 40), and PTPRr also inactivates Erk1/2. Mice that lack the PTPRr gene display impaired motor coordination (8). Thus, all PTPs that we identified here have earlier been shown to negatively regulate Erk1/2 signaling cascades. The Erk1/2 pathway has been implicated earlier in "promoting dedifferentiation of myelinating cells" (17). In a different context, Erk1/2 activating stimuli were recently shown to increase mRNA levels for eight DUSPs, and an siRNA screen showed that 12 of 16 DUSPs (including DUSP5) influence ERK2 responses (6).
Based on these recent observations, it appears plausible that the induction of PTP/DUSP mRNAs that we observed as a characteristic of early oligodendrocyte differentiation (shared by multiple treatments) may be a consequence of Erk activation and that interfering with induction of some of these PTP mRNAs prevents negative feedback to Erk signaling. Although genes that are activated during differentiation of oligodendrocytes or other cell types are often tacitly assumed to promote the differentiation process, our study indicates that the transcriptional machinery that drives differentiation also generates a number of "self-regulatory" transcripts that normally prevent cells from spontaneously changing their state.
Our finding that an MKP5 inhibitor induces OPC differentiation demonstrates that, at least for this gene, its "repressor function" is related to its enzymatic function. Earlier work on MKP5 knockout mice had shown that these animals show enhanced T-cell responses with increased production of Ifn-
and tumor necrosis factor alpha upon lymphocytic choriomeningitis virus infection and yet were, perhaps surprisingly, partially protected from experimental autoallergic encephalitis (45). Our study opens the possibility that MKP5 knockout animals are (also) protected from disease through enhanced oligodendrocyte-mediated remyelination.
More generally, our experimental approach demonstrates how functionally relevant genes can be extracted from the thousands that are modulated in a typical differential gene expression experiment, by intersecting datasets obtained for phenotypically converging conditions.
R.H.V.H. receives funding from the European Union Research Training Network MRTN-CT-2006-035830. AS077234 is available to members of the scientific community for noncommercial purposes upon request.
Published ahead of print on 12 January 2009. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
Present address: Intellikine, Inc., 10931 North Torrey Pines Road, La Jolla, CA 92037. ![]()
|
|
|---|
-tubulin. J. Neurosci. 27:2606-2616.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»