Abstract
The molecular changes underlying neural progenitor differentiation are essentially unknown. We applied cDNA microarrays with 13,627 clones to measure dynamic gene expression changes during the in vitro differentiation of neural progenitor cells that were isolated from the subventricular zone of postnatal day 7 mice and grown in vitro as neurospheres. In two experimental series in which we withdrew epidermal growth factor and added the neurotrophins Neurotrophin-4 or BDNF, four time points were investigated: undifferentiated cells grown as neurospheres, and cells 24, 48, and 96 hr after differentiation. Expression changes of selected genes were confirmed by semiquantitative RT-PCR. Ten different groups of gene expression dynamics obtained by cluster analysis are described. To correlate selected gene expression changes to the localization of respective proteins, we performed immunostainings of cultured neurospheres and of brain sections from adult mice. Our results provide new insights into the genetic program of neural progenitor differentiation and give strong hints to as yet unknown cellular communications within the adult subventricular zone stem cell niche.
- development
- differentiation
- hybridization
- ventricle
- stem cell
- microarray
- subventricular zone
- rostral migratory stream
- neurosphere
Introduction
Adult neural progenitor/stem cells in rodents and primates are well characterized in two brain regions: the subventricular zone (SVZ) of the lateral ventricles and the subgranular cell layer of the hippocampal dentate gyrus (McKay, 1997; Eriksson et al., 1998; Kornack and Rakic, 1999; Gage, 2000; Temple 2001; Alvarez-Buylla et al., 2002; Sanai et al., 2004). Stem cells in the adult SVZ are self-renewing and give rise to neurons that incorporate as interneurons into the olfactory bulb. Within the rodent SVZ, four major cell types have been described: ependymal cells lining the lateral ventricle, astrocytes (type B cells), and rapidly dividing transit-amplifying type C cells that generate neuroblasts (type A cells) (Doetsch et al., 1997). SVZ neuroblasts migrate as chains in the rostral migratory stream (RMS) to the olfactory bulb through tunnels formed by type B astrocytes (Lois et al., 1996; Peretto et al., 1997; Wichterle et al., 1997). The identity of the SVZ stem cell has been controversial. Originally proposed to correspond to a rare population of relatively quiescent cells (Morshead et al., 1994), ciliated ependymal cells were later suggested to represent the SVZ stem cell (Johansson et al., 1999). Findings of other groups, however, indicate that SVZ astrocytes (type B cells) are neural stem cells (Chiasson et al., 1999; Doetsch et al., 1999a,b; Laywell et al., 2000; Capela and Temple, 2002).
Cultured cells from the adult subventricular zone form floating multicellular spheres, so-called neurospheres, in the presence of epidermal growth factor (EGF) or basic fibroblast growth factor (bFGF) under appropriate conditions (Reynolds and Weiss, 1992; Richards et al., 1992; Morshead et al., 1994). These cells self-renew and are multipotent as they give rise to neuronal, astrocytic, and oligondendrocytic cells (Reynolds and Weiss, 1992). Doetsch et al. (2002) suggested that most of the EGF-responsive neurosphere-forming cells correspond to the rapidly dividing Distal-less homeobox 2-positive type C cells. Morshead and van der Kooy (2003) proposed that precursors of these cells give rise to neurospheres.
In principle, the in vitro differentiation of neurosphere-forming cells includes five steps: (1) cessation of proliferation, (2) attachment of the neurosphere, (3) detachment of cells from the neurosphere cell cluster, (4) migration of these cells away from the sphere, and (5) their actual differentiation into different cell types. These steps are comparable with adult neurogenesis of rodents in vivo; however, few of the underlying molecular changes are known. To learn more about the early and probably decisive mechanisms of this process, we have studied dynamic gene expression changes during the first 4 d of in vitro neurosphere differentiation with cDNA microarrays. We identified genes with relevant expression changes during this process and classified the results according to functional categories. Expression changes of selected genes were confirmed by semiquantitative RT-PCR. Because neurosphere cultures are not homogeneous (Suslov et al., 2002) and give rise to different cell types in vitro, we performed immunostainings to assign some of the most interesting results to specific cell types. We correlate our findings with the in vivo situation in the SVZ by staining tissue sections of adult mouse brains.
Our study revealed genes known or suggested to play a role in neurogenesis, but in addition revealed many new interesting genes with potential relevance for the maintenance, migration, and differentiation of neural progenitor cells.
Materials and Methods
For the nomenclature of mouse genes and proteins, we followed the guidelines of the international committee on standardized genetic nomenclature for mice (http://www.informatics.jax.org/mgihome/nomen/).
Cell culture. Neural progenitor cells were prepared from postnatal day 7 (P7) CD1 mice (Charles River, Wilmington, MA). Cells obtained from 10 animals were used for each series [BDNF and Neurotrophin-4 (NT4), see below]. Briefly, animals were killed, and tissue including the subventricular zone was dissected out from frontal sections. The tissue was digested with Accumax (PAA Laboratories, Coelbe, Germany) for 10 min at 37°C, fragmented with a cell strainer (40 μm; BD Falcon, Heidelberg, Germany), and plated out in suspension cell flasks (Nunc, Wiesbaden, Germany) at 50,000-100,000 cells per milliliter in DMEM/F12 (1:1) with Glutamax, B27, HEPES (10 mm), insulin (20 μg/ml), penicillin (100 U/ml), streptomycin (100 U/ml), fungizone (2.5 μg/ml), and human recombinant EGF (20 ng/ml) (all from Invitrogen, Karlsruhe, Germany). Cells were incubated at 37°C in 5% CO2. After 8-11 d in culture, neurospheres were dissociated by incubation with Accumax, as above, followed by careful trituration, and replated at 50,000-100,000 cells per milliliter. For differentiation, tertiary spheres were used after 27 d in culture. Spheres were plated out in cell culture dishes coated with poly-l-lysine (Sigma, Munich, Germany) in the same medium without EGF, but with human recombinant BDNF or NT4 (20 ng/ml; Promega, Mannheim, Germany). For immunofluorescence analysis, spheres were grown on poly-l-lysine-coated coverslips.
For differentiation in Matrigel, spheres were resuspended in ice-cold Matrigel (200 μl; BD Biosciences, Heidelberg, Germany) diluted with 75 μl of L15-medium (Sigma) and spread on coverslips. After incubation at 37°C for 30 min, Neurobasal medium supplemented with B27, HEPES (10 mm), penicillin (100 U/ml), streptomycin (100 U/ml), fungizone (2.5μg/ml), glutamine (2 mm), and human recombinant EGF (20 ng/ml) (all from Invitrogen) was added, and cells were incubated for 3 d at 37°C in 5% CO2.
Pictures of cultured cells were taken with an Axiovert 100 microscope (Zeiss, Oberkochen, Germany) and a digital CCD camera DSC-S75 (Sony, Tokyo, Japan).
Generation of cDNA microarrays. Murine cDNA clones (n = 10,080; arrayTAG clone collection) were purchased from LION Bioscience (Heidelberg, Germany), and 3510 murine cDNA clones representing brain-expressed transcripts were obtained from the resource center of the German Human Genome Project (RZPD Berlin). Clone inserts were amplified under standard conditions. Primers Lion 3′ (TCG AGC GGC CGC CCG GGC AGG T) and Lion 5′ (AGC GTG GTC GCG GCC GAG GT) were used for clones from LION bioscience, and primers M13for (GTA AAA CGA CGG CCA G) and M13rev (CAG GAA ACA GCT ATG AC) were used for RZPD clones. Primers were obtained from MWG Biotech (Ebersberg, Germany). In addition, 34 plant cDNA sequences and 3 murine cDNA sequences were amplified with specific primers. All PCR products were evaluated by agarose gel electrophoresis. Purified PCR products were resuspended in 8 μl of 3× SSC and printed on Corning GAPS II slides by using a robotic spotting device (SDDC-2 MicroArrayer, ESI, Toronto, Canada/ChipWriter Pro, Bio-Rad, Munich, Germany) with SMP3 pins from TeleChem International (Sunnyvale, CA).
RNA isolation and microarray hybridization. Total RNA was isolated with Trizol reagent (Invitrogen). Microarray labeling and hybridization reactions were performed using the 3DNA Array 50 Expression Array Detection kit (Genisphere, Hatfield, PA) according to the manufacturer's instructions. Briefly, 20 μg of total RNA of undifferentiated and differentiated cells (24, 48, or 96 hr) was reverse transcribed using primers with a specific capture sequence. The two cDNA samples were pooled and hybridized to the microarrays in a humidified chamber at 42°C for 16 hr. Dye-swap repeats were performed for each time point. Thus, six arrays were used for each series (NT4 and BDNF). Slides were washed at room temperature in 2× SSC/0.2% SDS for 10 min, 2× SSC for 10 min, and 0.2× SSC for 10 min. Bound cDNA was visualized by hybridization with 3DNA (dendrimers specific for respective capture sequences conjugated to either Cy3 or Cy5) in a humidified chamber at 42°C for 3 hr. Subsequently, slides were washed as before.
Image acquisition and data analysis. Fluorescence intensities of Cy3 and Cy5 were measured separately at 532 and 635 nm with a laser scanner (Affymetrix 428 Array Scanner, Affymetrix). The resulting 16-bit data files were imported into Microarray Suite image analysis software (version 2.0), which runs as an extension of IPLab Spectrum software (Scanalytics, Fairfax, VA). Raw spot intensities of Cy3 and Cy5 were locally background subtracted by the MicroArray Suite software. Empty spots and spots carrying plant gene sequences were excluded from further analysis. Each dye swap experiment was normalized by applying variance stabilization (Huber et al., 2002) using the vsn package of bioconductor (http://www.bioconductor.org). Means of normalized log-products and log-ratios of each dye swap experiment pair were used for further analysis. Normalization procedures were performed using R (http://cran.R-project.org).
Identification of differentially expressed genes and cluster analysis. First, the dataset was reduced to 454 clones that showed an estimated fold change of at least 2.0 at one differentiation time-point in both experiments (BDNF and NT4).
The threshold of estimated 2.0-fold change was determined by a statistical analysis aimed at minimizing the percentage of false positives. To do so, a variance estimation using a pooled estimate of the variance over all genes of three self-to-self comparisons with RNA from undifferentiated cells was performed. A similar approach has been described by Sabatti et al. (2002). To use a robust variance measurement, we determined the median of absolute deviation (MAD) as variance estimator. The MAD in all three independent experiments of self-to-self comparison was very similar: 0.297 ± 0.021. On the basis of this analysis, we can assume a rate of 2% false positives when applying a universal threshold of 2.17-fold change. A rate of 5% false positives can be assumed when applying a threshold of 1.8-fold. We therefore chose to consider all genes above a 2.0-fold change as relevant. Thus, we can assume a false positive rate between 2 and 5%, which we consider tolerable.
Of the clones shown in Table 1, we sequence-verified those that were present only once on the array. If a gene was represented by more clones on the array, which showed comparable expression changes, we did not sequence-verify. In this case, one representative clone is listed in Table 1.
Gene expression changes in neural progenitor cells during differentiation
Of the 454 clones, 441 clones showing a consistent upregulation or downregulation at the three differentiation time points were included in a cluster analysis. We clustered the 441 datasets, each representing a single time course, by applying the k-means algorithm and using a refined Euclidean distance measure. Ten cluster centers were randomly defined. The distance was defined as the weighted sum of k-means assignment and a similarity of shapes between cluster centers (gradient). The distance measure (D) was defined as follows: D(x,y) = a d1(x,y) + (1 - a) d2(x,y), where a = 0.5, x,y = two gene profiles to be compared, d1 represents the Euclidean distance between x and y, and d2 represents the gradient of x and y. These calculations and the respective visualization were performed using MATLAB (version 6.0.0.88, release 12, MathWorks).
Reverse transcription and semiquantitative PCR. Two micrograms of total RNA from each of the four time points were reverse-transcribed into cDNA by incubation with 5 μl of Deoxyribonucleoside triphosphate solution (dNTPs) (10 mm of dATP, dCTP, dTTP, dGTP each), 2.5 μl of random primer (500 ng/μl), nuclease-free water, 10 μl of Moloney murine leukemia virus (M-MLV) 5× reaction buffer, 2 μl of RNasin, and 800 U of M-MLV RT (all from Promega) for 1 hr at 37°C. PCR reactions were performed using equal amounts of cDNA as template, 2.5 μl of dNTPs (2 mm each), 2.5 μl of 10× reaction buffer, 1.5 μl of each primer (10 pmol/μl), water, and 2 μl of Red Taq (Promega). Amplifications consisted of 21-26 cycles (94°C for 45 sec, 54-57°C for 60 sec, and 72°C for 90 sec). PCR products were visualized by agarose electrophoresis. Gene-specific primers and annealing temperatures are provided as supplementary material.
Immunofluorescence. Cells were grown on coverslips, fixed by sequential incubation in ice-cold methanol (5 min) and acetone (30 sec), and stored at -20°C. Cryostat sections (10 μm) of brains from 3-month-old CD1 mice were placed on glass slides, fixed in -20°C cold acetone for 10 min, and air dried. If not used immediately, slides were stored at -80°C until use. Cells used for galactosyl cerebroside (GALC) staining were fixed in 4% paraformaldehyde/PBS.
Sections and coverslips were treated with Triton X-100 (0.1% in PBS) for 5 min (except for GALC staining), washed with PBS three times (5 min each), and blocked with 5% donkey serum (Dianova, Hamburg, Germany) in PBS for 30 min at room temperature. Primary antibodies were applied for 16 hr at 4°C or for 45 min at room temperature. PBS without primary antibody was used as a negative control. The following primary antibodies were used: tubulin-β-III (1:150; Upstate, Charlottesville, VA), glial fibrillary acidic protein (GFAP) (polyclonal, 1:100; Santa Cruz Biotechnology, Santa Cruz, CA), GFAP (1:100, MAB3402 monoclonal; Chemicon, Temecula, CA), nestin (1:100; Transduction Laboratories, Lexington, KY), Ki-67 (1:100; Dako Cytomation, Hamburg, Germany), protein tyrosine phosphatase, nonreceptor type substrate 1 (PTPNS1) (1:200; BD Biosciences), CD47 antigen (CD47) (1:200; BD Biosciences), FABP7 (1:100; kindly provided by T. Mueller, C. Birchmeier, MDC, Berlin, Germany), pleiotrophin (1:100; R & D Systems, Minneapolis, MN), HUC/D (1:100; Molecular Probes, Eugene, OR), CD24 antigen (CD24) (1:100; Santa Cruz Biotechnology), polysialylated neural cell adhesion molecule (PSA-NCAM) (1:100; AbCys, Paris, France), galactosylceramidase (GALC) (1:100; Roche). After three washes in PBS (15 min each), secondary antibodies coupled to Cy3 or fluorescein isothiocyanate (Dianova) were applied for 30 min at room temperature. Specimens were then washed three times in PBS (15 min each), rinsed briefly in water, dehydrated in ethanol, and air dried. A mounting medium containing 4′,6′-diamidino-2-phenylindole (DAPI) nuclear stain was used. Immunofluorescence was documented with a Zeiss Axioskop epifluorescence microscope (Zeiss) or with a Leica DM IRE2 (Leica, Wetzlar, Germany).
Results
In vitro differentiation of neurosphere cells
Soon after the transfer of neurospheres into poly-l-lysine-coated dishes with a medium that contained no EGF, but did contain BDNF or NT4, they attached to the plate. By 24 hr a few cells had migrated out of the spheres. This number increased with time, and migrating cells underwent morphological changes (Fig. 1). Cells in close proximity to the sphere appeared as large flat glia-like cells, and they were immunopositive for the astrocyte marker GFAP (Fig. 2). In contrast, some cells that migrated farther away from the sphere acquired a neuronal morphology with small cell bodies and long, slim processes. These stained positive for the neuron-specific marker tubulin-β-III (Fig. 2). Most cells within the spheres were positive for GFAP; however, after 96 hr the compact sphere structure had frequently loosened up, and some of these cells had changed their morphology as well. Few cells inside the spheres produced tubulin-β-III after induction of differentiation, whereas most differentiating cells were GFAP positive. Very few cells were positive for the oligodendrocyte marker GALC (data not shown).
Photographs depicting morphological changes in the course of neural progenitor cell differentiation in vitro. Undifferentiated neurospheres before and cells at 24, 48, and 96 hr after induction of differentiation are shown. Scale bar, 50 μm.
Immunofluorescence microscopy of undifferentiated neurospheres and of cells at 24, 48, and 96 hr after induction of differentiation. Cells were immunolabeled with antibodies against tubulin-β-III, Ki-67, Nestin, PTPNS1 (SHPS-1), CD47, FABP7, Pleiotrophin (all red), and GFAP (green). Nuclei are stained blue. Note the different staining of undifferentiated neurosphere clusters and differentiating cells. Scale bar, 100 μm.
To test for cell proliferation, we immunostained with an antibody against Ki-67, a marker for dividing cells (Gerdes et al., 1983). Before induction of differentiation, some cells within neurospheres were positive for Ki-67; however, 24 hr after induction of differentiation, no more dividing cells were found (Fig. 2). This is in line with our microarray results (see Gene expression changes related to cell proliferation). We also immunostained the cells for nestin, a marker present in neuroepithelial precursor cells. Before and after induction of differentiation, neurosphere clusters were nestin-positive; however, cells that had migrated away from the spheres no longer expressed nestin (Fig. 2).
Gene expression changes: general findings
We performed cDNA microarray co-hybridizations with cDNA samples derived from cells of two experimental series (BDNF and NT4). Neurosphere cells were kept in culture for two passages and then pooled and redistributed before differentiation with either BDNF or NT4. Cells were harvested at time points 0 (undifferentiated cells), 24, 48, and 96 hr after induction of differentiation. For each series, three co-hybridizations, 0 versus 24 hr, 0 versus 48 hr, and 0 versus 96 hr were performed and repeated with fluorescent dyes swapped to compensate for dye-specific effects. The six resulting single datasets per series were subjected to further analysis that involved a normalization step (see Materials and Methods), resulting in three datasets per series, one for each co-hybridization (24, 48, and 96 hr). To identify relevant differentially expressed genes, all clones with more than twofold estimated differences in expression in at least one of the three datasets were considered for further evaluation. This threshold is based on a statistical analysis as described in Materials and Methods.
Seven hundred twenty-two clones of the BDNF experiment and 624 clones of the NT4 experiment fulfilled this criterion. Eighty-three percent of the clones from the BDNF series showed at least a 1.6-fold estimated change in the NT4 experiment, and 84% of the clones from the NT4 series showed at least a 1.6-fold estimated change in the BDNF experiment. Four hundred fifty-four clones exhibited relevant (more than twofold) expression changes under both conditions. Of these, 97.1% were consistently upregulated (271 clones) or downregulated (170 clones) at the three time points (negative or positive log ratio of signal intensity) compared with undifferentiated cells in both series. Thirteen clones (2.9%) showed an upregulation at one time point and subsequent downregulation or vice versa (e.g., phosphoglucose isomerase 1) (Table 1, metabolism, NT4 series). The overlap of results from these two independent experiments suggests that the neurotrophins BDNF and NT4, both of which bind to the TrkB receptor, lead to mostly overlapping gene expression changes. All further analyses were limited to the 454 clones of the intersection list, thus concentrating only on the commonalities associated with the TrkB receptor activation in the two experimental series. They were grouped into 10 different categories (Table 1) (see below).
The complete dataset (MIAMI standards) including remaining clones with relevant expression changes, which are not listed in these groups, can be found on our homepage (http://www.molgen.mpg.de/~dna_microarrays/neural_differentiation/neural.html).
To extract the dynamic behavior of the relevant gene expression changes in the course of neural progenitor differentiation, cluster analysis was performed using the 441 consistently upregulated or downregulated clones from the intersection list (see Materials and Methods). We grouped the 441 clones into 10 clusters that show a similar onset and course of transcriptional changes (Fig. 3A): clusters 1-5 containing upregulated and clusters 6-10 containing downregulated genes, both in the BDNF and NT4 series. With the exception of clusters 7 and 8, there is a large fraction of clones in each cluster that occur in both series (Fig. 3B, Table 1). Most genes fall into clusters that display only slight expression changes (clusters 2, 4, 5, 6, 7, 10). Clusters 1 and 3 contain strongly upregulated genes, and clusters 8 and 9 contain strongly downregulated genes. Concerning the dynamics, genes that change late in expression are found especially in cluster 2. Genes that change very early (between time point 0 and 24 hr) are found in two of the four clusters with strongest changes (clusters 1 and 9); in addition, they are found in cluster 6. Downregulated genes of the functional category “chromatin-associated components and nuclear factors,” in particular, fall into clusters 6 and 10, which contain slightly decreasing transcript levels. Cluster 9 includes numerous cell cycle-related genes. Cluster 2 (late upregulated transcripts) contains many genes that localize to the cell membrane or are found in cellular extensions. Apart from these findings, single clusters do not preferentially contain genes of a certain category. The 10 genes with highest expression changes at one time point were Igfbp2 (16.28/25.79) (highest fold-change in BDNF experiment/highest fold-change in NT4 experiment), 6330403K07Rik (10.91/14.01), PDZ-binding kinase (10.49/7.24), Cdc2/p34 (10.07/9.49), Igfbp4 (9.87/8.85), Fxyd1/Plm (7.85/7.77), Cyclin D1 (6.36/7.69), Complement component 3 (7.61/6.89), Cdc20 (7.10/6.96), and an EST adjacent to cadherin 13 (6.69/5.87).
A, Clusters of clones with similar expression changes over time. The dynamic behavior of the 441 clones that are consistently upregulated or downregulated in both the BDNF and the NT4 series is shown separately. For each series, 10 clusters were found. Clusters 1-5 contain upregulated clones; clusters 6-10 contain downregulated clones. The x-axis depicts the three time points of differentiation. The y-axis shows relative fold changes; that is, expression changes referred to the undifferentiated state (0 hr). These relative numbers are estimated logarithmic fold changes and therefore not identical to the estimated fold changes listed in Table 1. The logarithmic values allow a better graphic presentation. Note that y-axis values at 24 hr, which are different from the value 0, indicate that the expression of a clone changed already between 0 and 24 hr. B, Numbers of clones contained in each cluster of the NT4 and BDNF series.
Cell type-specific gene expression changes
We investigated very early differentiation stages of neurosphere cells to determine the basic molecular changes that underlie the conversion of the progenitor cell toward a cell that enters a differentiation process. In vitro, neurosphere cells differentiate into three major types of cells: neurons, astrocytes, and oligodendrocytes (Reynolds and Weiss, 1992). In addition to a functional categorization, we also screened our data for cell type-specific genes and found expression changes that already point to a cellular specification. Several genes typical for neurons and oligodendrocytes were found, all of them upregulated, and many only at a late time point of differentiation (96 hr). With the exception of the myelin basic protein gene in the NT4 series, none of them falls into cluster 1. Two genes characteristic for astrocytes were found: S100B is upregulated, deiodinase iodothyronine type II is downregulated.
Concerning oligodendrocyte-specific genes, we found an increased expression of many myelin components. Of note is the temporal sequence of myelin component upregulation: transcript levels of myelin basic protein increase first (already at 24 hr), whereas myelin-associated glycoprotein and myelin-associated oligodendrocyte basic protein transcripts increase more than twofold only by 96 hr.
Of neuron-specific genes, increased expression was found for Nsg2 (neuron-specific gene family member 2), Chromogranin B, the RIKEN cDNA 4930565N16 gene homologous to rat Tomosyn, Syntaphilin, and Reticulon 1. Additional genes confined to neurons in the nervous system were Neurexin II β, RhoN, Neuronal tropomodulin, 4.1N, tubulin-β-III, and Darpp-32 (see below).
Gene expression changes related to cell adhesion
Although the molecular basis for the formation of neurospheres in vitro is not known, it could be mediated by direct cell-cell adhesion or indirectly via extracellular matrix (ECM) proteins that act as adhesive substance between them. Searching our dataset for genes involved in cell adhesion and defined cell junction structures, we found Cadherin13 (H-cadherin, T-cadherin), an EST adjacent to cadherin 13, and the gene encoding Chl1 (close homolog of the neural adhesion molecule L1) to be downregulated in the course of differentiation, indicating their expression in undifferentiated neurospheres.
The most striking decrease, approximately sixfold at 96 hr, was seen in the case of the EST adjacent to cadherin 13. In addition, a cadherin 13 clone showed a strong transcript decrease. Cadherin 13 is a glycosylphosphatidyl inositol (GPI)-linked membrane protein that can mediate homophilic cell-cell adhesion (for review, see Takeuchi and Ohtsuki, 2001). Moreover, anti-adhesive properties and a function in cell cycle regulation have been reported (Takeuchi et al., 2000; Huang et al., 2003; Ivanov et al., 2004). Chl1 belongs to the L1 family of neural recognition molecules and promotes neurite outgrowth and neuronal survival and enhances integrin-mediated cell migration (Hillenbrand et al., 1999; Buhusi et al., 2003).
Some genes of the cell adhesion category were upregulated, suggesting a role in cellular contacts between differentiating cells: Cd24, Emp3, Neurexin II, and Plakophilin-4. The first two fall into cluster 2, which contains very late gene expression changes. None of them belongs to cluster 1. CD24 is a membrane-anchored glycoprotein expressed in the SVZ by neuroblasts (type A cells), in the SVZ ependymal cells, and in the brain parenchyma by differentiated neurons (Calaora et al., 1996; Doetsch et al., 2002). It binds to L1 on neighboring cells and has been shown to influence neurite outgrowth. The Neurexin II cDNA clone on our array shows highest similarity to a Neurexin II β variant (for review, see Missler and Sudhof, 1998). β Neurexins have been implicated in heterophilic cell-cell adhesion and assembly of presynaptic terminals (Dean et al., 2003). The gene encoding Plakophilin-4 (p0071 protein/Nprap), an intracellular plaque protein at adherens junctions and desmosomes in epithelial cells, shows gradually increasing transcript levels (cluster 5). The Emp3 (Hnmp-1) product belongs to the PMP22 family of tetra-spanning membrane proteins and is thought to be involved in cell-cell interactions and cell proliferation. The Sh3 domain protein 4 (vinexin) binds vinculin and localizes to sites of cell-matrix interactions (focal adhesions) and at cell-cell adhesions. It will be of interest to investigate the upregulated genes on the protein level to assign them to specific cell types.
Gene expression changes related to extracellular matrix components
ECM molecules serve as a scaffold for the aggregation of cells within tissues. In addition, they constitute a cellular microenvironment that influences cell behavior, including cell migration and differentiation. These processes are controlled in part by biologically active molecules such as growth factors, which are stored in the ECM.
Our study revealed three transcripts encoding ECM proteins that are downregulated in differentiated cells compared with undifferentiated, floating cells: Tenascin C, Collagen XVIII, and Laminin receptor 1.
Tenascin C is present in the SVZ, RMS, and olfactory bulb. Specifically, it corresponds to GFAP-positive processes that form the tunnel through which SVZ neuroblasts migrate (Jankovski and Sotelo, 1996). In the embryonic brain, it is expressed by radial glial cells (Malatesta et al., 2003). Tenascin mediates adhesion of cells, either increasing or decreasing it (Jones and Jones, 2000). Collagen XVIII is a basement membrane component. One splice variant is expressed in brain (Muragaki et al., 1995; Saarela et al., 1998; Sertie et al., 2000). The downregulation of the gene encoding Laminin receptor 1 is in line with its expression in undifferentiated, mitotically active progenitor and tumor cells in contrast to differentiated forms of these cells (for review, see Menard et al., 1997).
ECM genes that were upregulated in our study encode the Adipocyte enhancer binding protein 1 (Aortic carboxypeptidase-like protein), Extracellular matrix protein 2 (Sparcl1, Mast9, Hevin), Collagen IX α 3, an extracellular metalloprotease (Adamts10), and several inhibitors of proteinases that participate in ECM degradation: TIMP2, TIMP3, Serpina3k, and Inter-α trypsin inhibitor heavy chain 3.
The increased expression of these inhibitors suggests a stabilization of certain ECM components of differentiating cells. In adult rats TIMP3 is localized in the SVZ, the rostral migratory stream, and the olfactory bulb. It has been discussed as a mediator of neuronal migration by regulating axon guidance and process outgrowth (Jaworski and Fager, 2000).
Gene expression changes of complement components
Complement regulatory proteins and receptors are expressed by neurons and glial cells of the CNS (Nataf et al., 1999; Thomas et al., 2000; Yu et al., 2002; Jauneau et al., 2003). We found increased expression during neurosphere differentiation of four genes belonging to the complement cascade. These genes encode C3 and C4 precursors, the C1qb component, and Complement factor I precursor, which inactivates C3b and C4b. C3 and C4 belong to cluster 1, which contains strongly upregulated genes.
Gene expression changes related to the actin cytoskeleton
We found many differentially expressed genes, mostly upregulated, that are involved in morphological changes expected to take place during the switch from a round neurosphere cell that is held in a floating cell ball to an adherent cell that extends cellular processes. The building and breakdown of actin filaments (actin turnover) at the cellular periphery is an indispensable event for the establishment of cellular protrusions, including the extension of dendrites and axons. Furthermore, it is important for the extension of leading edges as one mechanism of neuronal migration. Many factors control actin turnover at these structures: upstream acting key mediators including Rho family GTPases, plexins, which regulate axon guidance through these GTPases, and actin-binding and actin-remodeling proteins (Meyer and Feldman, 2002; Small et al., 2002).
A member of the Rho family GTPases, Rho7 (RhoN), and genes coding for proteins that influence GTPase activity, e.g., PlexinB3, Arhgap5, MgcRacgap, Active bcr-related gene, and Cdc42 effector protein 4 gene increased in expression, whereas Arhgap18 and N-chimerin decreased. Several genes the products of which directly or indirectly bind actin, some of which influence actin polymerization, were upregulated: Non-muscle α actinin 1, Cysteine-rich protein 1 (Crp1), Coronin actin binding protein 1C, Gelsolin (ADF/brevin), Neuronal tropomodulin (tropomodulin 2), and an EST homologous to the human actin binding LIM protein 1 and similar to Caenorhabditis elegans UNC-115. Prothymosin β 10 transcripts decreased.
An increase in transcription was also observed for three genes encoding proteins implicated in the linkage of actin to the cell membrane: Ezrin (Villin2), Neuronal protein 4.1, and Adducin 3 (Gamma adducin). Ezrin localizes to fine, peripheral astrocytic processes but not to neurons (Derouiche and Frotscher, 2001). A function in neuronal development, however, was proposed on the basis of its localization at neuronal growth cone filopodia and perikarya in vitro (Goslin et al., 1989; Birgbauer et al., 1991; Paglini et al., 1998) and on its neuronal localization in the chick embryo (Everett and Nichol, 1990; Takahashi et al., 1999). Adducin is important for the regulation of membrane ruffling and cell motility. It localizes to dendritic spines of neurons and is associated with regions of morphogenetic cell movement in the chick embryo (Matsuoka et al., 1998; Akai and Storey, 2002).
Furthermore, four genes encoding calcium-binding proteins of the S100 family increase in expression: S100A1, S100A4, S100A6, and S100B. Calcium-binding proteins of this family exert multiple intracellular functions (regulation of enzyme activities, modification of the actin cytoskeleton, calcium homeostasis, cell growth, and differentiation). Several of them also have extra-cellular functions, including regulatory effects on neural cells (for review, see Donato, 2003). All but one of the identified S100 genes (S100A6/Calcyclin) promote neurite extension as an extracellular effect. Furthermore, S100A4 regulates metastasis and is intracellularly present at the leading edge of lamellipodia where it interacts with Myosin heavy chain IIA (Kim and Helfman, 2003). S100B is strongly expressed by astrocytes, and the secreted form also influences neuronal survival.
Gene expression changes related to the tubulin cytoskeleton
Besides the actin cytoskeleton, components of the microtubule filament system changed in expression during neurosphere differentiation. All of them were upregulated and, with one exception, changed late in our experiments (none above twofold before 48 hr). They include the tubulin-β-III, β-4, α-3/7 genes and the Sperm-associated antigen 6 (Spag6). Furthermore, they encode microtubule motor proteins, the cytoplasmic dynein 2B light chain, the axonemal dynein heavy polypeptide 7, and the kinesin heavy chain member 5B.
SPAG6 is found in motile structures (motor cilia and flagella) that contain an axoneme (“9 + 2” structure) built from nine outer microtubule pairs and one inner microtubule pair (called the central apparatus) together with dyneins and is required for their structural integrity and function (Sapiro et al., 2002). Within the axoneme, SPAG6 is localized to the central apparatus. Of CNS cells, ependymal cells lining the ventricles carry multiple motor cilia. Type B SVZ astrocytes often have a single cilium with a “9 + 0” microtubule structure. Both types of cilia extend into the lateral ventricle space (Doetsch et al., 1997, 1999). The latter structure is found in sensory cilia and in the unique motile monocilia of embryonic node cells (Brody et al., 2000).
Microtubule motor proteins fulfill a specialized function in motor cilia and flagella. In addition, these proteins mediate the transport of cargos along microtubules, including synaptic vesicle transport within axons.
Gene expression changes related to cell proliferation
Our microarray results indicate that already during the first 24 hr after induction of differentiation, neurosphere cells cease proliferating and exit mitosis. All genes related to DNA synthesis and cell cycle progression identified in this study are contained in clusters of downregulated genes. Many fall into cluster 9, which describes the strongest decrease. Among these are Cyclins B1, D1, D2, the Cell division control protein 2 (Cdc2/p34), and the Cell division cycle protein 20 (Cdc20). These are important factors for cell cycle progression. We confirmed the transcript level decrease of Cdc2/p34 by semiquantitative RT-PCR (Fig. 4). The PDZ-binding kinase gene, a serine/threonine kinase, was also strongly downregulated. The corresponding human protein was identified as an associated protein of the tumor suppressor hDlg and is an in vitro substrate of CDC2/Cyclin B, which is phosphorylated in a cell cycle-dependent manner (Gaudet et al., 2000). Furthermore, transcripts of genes encoding proteins involved in DNA replication decreased. These include genes for structural components of the replicating enzyme (DNA polymerase δ 1, DNA polymerase δ subunit p12) as well as the DNA replication licensing factors Mcm3 and Mcm6. In addition, we found decreased expression of Chromatin assembly factor 1 subunit C (Caf-1 subunit C), Kinesin family member 23 (Kif23), Pituitary tumor-transforming protein 1 (Pttg/Securin), and Fusion 1 protein (Fus1). CAF-1 is thought to mediate chromatin assembly in DNA replication by recruiting newly synthesized histones H3 and H4 to replicating DNA (Verreault et al., 1996). KIF23 is homologous to human kinesins involved in mitosis. Securin is a key regulator of sister chromatid separation during anaphase, and Fus1 is a candidate tumor suppressor gene possibly involved in G1 arrest (Kondo et al., 2001).
Semiquantitative RT-PCR analysis of selected genes. Relative gene expression levels before and after induction of differentiation are shown. The respective transcripts were upregulated or downregulated as expected from the microarray results. Expression of Txrnd1 and Mapk3/Erk-1 was used as a reference.
Gene expression changes of signal transduction molecules
Most of the large group of differentially regulated genes involved in signal transduction pathways fall into clusters that describe slight expression changes.
We found three downregulated and two upregulated genes encoding components of tyrosine kinase-mediated pathways. These are Igfbp2, Igfbp4, and Cnil (downregulated), and Igfbp5 and Ptpns1 (Shps-1) (upregulated). The expression change of Igfbp2 was confirmed by semiquantitative RT-PCR (Fig. 4). The Igfbp genes are exceptional for the signal transduction category because they fall into clusters 1 and 9 of strongest and earliest expression change. Of all 454 clones with relevant (more than twofold) changes in both series, those representing the Igfbp2 gene show the highest fold changes in our study: up to 25.79-fold decrease in transcript level. A dramatic decrease is especially seen between 24 and 48 hr. In addition to Igfbp2, Igfbp4 is strongly downregulated. Both IGF-binding proteins inhibit IGF actions. Igfbp5, on the other hand, increases in expression and is known to potentiate IGF-I effects. Thus, differentiating cells are apparently converted into an IGF-responsive state. IGF-I stimulates neurogenesis, oligodendrogenesis, synaptogenesis, and cell migration, and it inhibits apoptosis (D'Ercole et al., 2002; Hsieh et al., 2004). In addition, IGFBP5 can exert IGF-independent effects and has also been found in the nucleus (for review, see Firth and Baxter, 2002; Mohan and Baylink, 2002).
Ptpns1 (Shps-1/Sirpα/Bit/Mfr/p84 neural adhesion molecule) encodes a transmembrane protein, which acts as an important regulator of IGF-I signaling (Clemmons and Maile, 2003). Our microarray study revealed an expression increase of Ptpns1 (Shps-1) already at 24 hr and even more at 48 and 96 hr. This result was confirmed by semiquantitative RT-PCR (Fig. 4).
The PTPNS1 (SHPS-1) ligand gene Cd47 [Integrin-associated protein (IAP)], which also encodes a transmembrane protein, slightly decreased in expression (Fig. 4). Immunostainings of neurospheres revealed that both the CD47 and PTPNS1 (SHPS-1) protein are synthesized by undifferentiated neurosphere cells that form tight cellular clusters (Fig. 2). After differentiation, PTPNS1 (SHPS-1) is still present in the sphere, but is also synthesized by cells that have dissociated from the sphere and surround the compact cell cluster. Cells that have migrated farther away do not stain with the anti-PTPNS1 (SHPS-1) antibody. In addition to the signal of compact neurosphere cells, we detected a faint CD47 staining of single cells that had migrated away from the sphere. In our in vitro system, cells migrate as single cells. In vivo, however, rodent neuroblast cells migrate as packed chains with tight cell-cell contacts. Using an in vitro chain migration assay of neurosphere cells by plating them in Matrigel (Wichterle et al., 1997), we demonstrated that cells migrating in chains synthesized CD47 as detected by immunofluorescence staining (Fig. 5a). Migrating cells also showed PTPNS1 (SHPS-1) immunoreactivity (Fig. 5b,c), albeit less intensely so than CD47.
Immunolocalization of CD47 (a, red signals) and PTPNS1 (b, c, red signals) in differentiating neurosphere cells that were embedded in Matrigel. Blue signals indicate DAPI-stained nuclei. Cellular chains are indicated by arrows. Scale bars, 50μm.
Adult sagittal brain sections immunoreacted with PTPNS1 (SHPS-1), and CD47 antibodies revealed a band of intense staining immediately underlying the less intensely stained SVZ (Fig. 6a-d, arrowheads). Weaker signals were detected in what appeared to be the cytoplasm of GFAP-positive astrocytes (Fig. 6a,b, arrows), and of HUC/D-positive neuroblasts (Fig. 6c,d, arrows). In the RMS, GFAP staining partially overlapped with that of PTPNS1 and CD47 (Fig. 6e,f, arrows), which was also the case for PSA-NCAM with CD47 and PTPNS1, respectively (Fig. 6g,h, arrows)
Localization of PTPNS1 (a, c, e, g) and CD47 (b, d, f, h) on sagittal adult brain sections as determined by immunofluorescence microscopy. Regions of the lateral ventricle wall (a-d) and the RMS (e-h) are shown. Coimmunostainings with GFAP (a, b, e, f, green signals), HUC/D (c, d, red signals), and PSA-NCAM (g, h, red signals) were performed. c, d, g, h, PTPNS1 and CD47 signals are shown in green. a, b, e, f, PTPNS1 and CD47 signals are shown in red. Note the band-like PTPNS1 and CD47 staining pattern (a-d, arrowheads). Examples of cells positive for GFAP and PTPNS1 (a) or GFAP and CD47 (b) are marked by arrows. c, d, Arrows depict examples of cells positive for PTPNS1 and HUC/D (c) or CD47 and HUC/D (d). e-h, Arrows indicate an overlapping staining pattern of PTPNS1 and GFAP (e), CD47 and GFAP (f), PTPNS1 and PSA-NCAM (c) as well as CD47 and PSA-NCAM (h). Scale bars: a, e, 100 μm; b-d, f, h, 50 μm.
Four genes coding for factors involved in G-protein-linked pathways increase in expression: Edg3, Neurotensin receptor 2, Rgs2, and Rgs16. Edg3 (Endothelial differentiation gene 3/S1p3), a widely expressed G-protein-coupled receptor for the sphingolipid Sphingosine-1-phosphate (S1P), was upregulated in the course of neurosphere differentiation. EDG3 is coupled to ERK/MAPK (extracellular signal-related kinase/mitogen-activated protein kinase), RHO (ras homolog), and RAC (ras-related C3 botulinum toxin substrate), Ca mobilization, transcription factors, and serum response factor (An et al., 2000). Of special interest for our study is the finding of EDG3-mediated induction of membrane ruffling, activation of RAC, and chemotaxis toward S1P (Kon et al., 1999; Wang et al., 1999; Okamoto et al., 2000). The ligand of Neurotensin receptor 2 has various functions in the CNS (Tyler-McMahon et al., 2000). On the cellular level, Neurotensin can lead to tyrosine phosphorylation of the focal adhesion kinase, a protein that localizes to integrin clusters and influences cell migration, adhesion, survival, and proliferation (Leyton et al., 2002).
Differentially expressed downstream components of the inositol phospholipid pathway were also found: Pten-induced kinase, PLC δ, and Cdp-diacylglycerol synthase increased and Pkc β decreased.
Two genes, cGMP-stimulated phosphodiesterase 2A and Darpp-32, both involved in cyclic monophosphate signaling, decreased.
We found three differentially expressed genes that belong to calcium-dependent signal transduction cascades: Calmodulin was upregulated, whereas the Calcium/Calmodulin-dependent protein kinase II β and the catalytic subunit of the Calmodulin-regulated protein phosphatase-2B (Calcineurin) decreased. Camk2b is striking because it belongs to cluster 9 of early and very strongly downregulated genes. We confirmed its transcript level change by semiquantitative RT-PCR (Fig. 4).
Two genes implicated in the thyroid hormone system changed their expression: Deiodinase iodothyronine type II (D2) was downregulated, and Thyroid hormone receptor interactor 6 (Trip6) increased. Deiodinase iodothyronine type II (D2) catalyzes the conversion of T4 into the active form of thyroid hormone, T3. In the brain, ∼80% of T3 is generated locally from T4 by the action of D2. This enzyme is expressed by glial cells in normothyroid postnatal rats (Guadano-Ferraz et al., 1999). The effects of T3 in the brain are diverse and include neuronal differentiation, migration, neurite outgrowth, synapse formation, and oligodendrocyte differentiation and myelination. Astrocytes treated with T3 secrete growth factors, including EGF, that in turn act on neurons (Bernal and Nunez, 1995; Gomes et al., 2001a,b; Lima et al., 2001; Forrest et al., 2002; Konig and Moura Neto, 2002; Martinez and Gomes, 2002).
Genes involved in the NF-κB-dependent signaling pathway were also noticeable: Tnfrsf21 and NF-κB target gene Ier3 transcripts decrease and levels of the Tnf receptor-associated factor 1 increase. In line with its designation, the immediate early gene Ier3 is contained in cluster 6, which represents slight but early expression changes (between time point 0 and 24 hr). Deltex1, a component of the Notch signaling pathway, increased slightly. The Notch signaling pathway is fundamental in cell-fate determination, including the development of neural progenitor cells into neurons and glial cells. Cytoplasmic Deltex proteins are positive regulators of this pathway that act independently of the Notch cascade, which involves RBP-J (recombining binding protein suppressor of hairless) and Hes proteins. Previous reports implicate the Notch pathway in astroglial versus neuronal cell-fate determination (Tanigaki et al., 2001); however, a recent publication demonstrates its function in the maintenance of neural stem cells (Hitoshi et al., 2002). Mammalian Deltex1 has been shown to act as a negative modulator of Notch1 signals (Sestan et al., 1999; Izon et al., 2002). It interacts with the intracellular domain of Notch1 and inhibits the action of neurogenic basic helix-loop-helix (bHLH) transcription factors Mash1 and E2A (Matsuno et al., 1998; Kishi et al., 2001; Yamamoto et al., 2001). Overexpressed Deltex1 inhibits neural progenitor cell differentiation and opposes the effect of Notch1 on neurite outgrowth (Sestan et al., 1999; Yamamoto et al., 2001). A Notch/Deltex1 signaling pathway that involves F3/contactin promotes oligodendrocyte maturation and myelination (Hu et al., 2003). Murine Deltex1 is expressed in the embryonic neuroepithelium. In later developmental stages its expression is restricted to postmitotic differentiating neurons and absent from the proliferative ventricular zone (Mitsiadis et al., 2001).
Finally, the Latent transforming growth factor-β binding protein (Ltbp) 4S, implicated in the Tgfβ pathway, was upregulated.
Gene expression changes related to neural progenitors
We identified three growth factor encoding genes, all of which were downregulated during differentiation: Hb-Egf (Heparin-binding EGF-like growth factor/Diphtheria toxin receptor), Pleiotrophin (Heparin-binding growth-associated molecule/Hb-Gam), and Nell2 (Neural epidermal growth factor-like-like 2). Thus, they are expressed by neural progenitor cells. The transcriptional changes of these genes were confirmed by semiquantitative RT-PCR (Fig. 4). The Pleiotrophin receptor gene Syndecan-3 also decreased in expression. In addition, we found a downregulation of the gene encoding brain fatty acid-binding protein (FABP7/BFABP/BLBP), a brain-specific member of the FABP family of cytosolic proteins that bind fatty acids, eicosanoids, and retinoids. FABP7 is characteristic for radial glial cells.
Immunostainings of neurosphere cells revealed that Pleiotrophin was restricted to compact cells in the center of the neurosphere and disappeared as cells dissociated from the sphere and differentiated (Fig. 2). In the rostral migratory stream, cytoplasmic Pleiotrophin immunoreactivity was associated with small vessels (Fig. 7a, thin arrow) as has been described previously (Yeh et al., 1998), but not with PSA-NCAM-positive neuroblasts (Fig. 7a). In the subventricular zone lining the walls of the lateral ventricles, we found evidence for cellular colocalization of Pleiotrophin with the neuronal marker HUC/D (Fig. 7c, arrowhead) and CD24 (Fig. 7d), which is present in ependymal cells and type A neuroblasts (Calaora et al., 1996). Pleiotrophin immunoreactivity was also co-detected in the same cellular SVZ layer as the glial marker GFAP (Fig. 7b), but because of the tight association of the glial processes with neuroblasts and ependymal cells in this layer, only electron microscopic analysis can unambiguously ascertain whether Pleiotrophin is expressed by glial cells. Nuclear Pleiotrophin staining in the superior-posterior ventricle wall (Fig. 7b-d) was mostly punctate (Fig. 7d, arrow), whereas more homogeneous nuclear staining was characteristically seen in other regions of the SVZ. Pleiotrophin has previously been shown to bind to nucleolin, a major nucleolar protein, which functions as a shuttle protein between the nucleus and cytoplasm (Take et al., 1994). In addition, we consistently found distinct Pleiotrophin-positive spots in the anterior wall abutting the striatum, the most anterior portion of the superior ventricle wall, as well as in its rostral extension (Fig. 7a,e, arrowheads). These spots appeared cytosolic and were difficult to assign to specific cell types.
Localization of Pleiotrophin (a-e) and FABP7 (f) on sagittal adult brain sections as determined by immunofluorescence microscopy. A nuclear staining of the Pleiotrophin antibody can be seen in cells of the lateral ventricle wall (b, red signal; c, d, green signals). b, Green signals indicate GFAP. c, Red signals indicate HUC/D. d, Red signals indicate CD24. b, Arrow points at a GFAP-positive cell. c, Arrow points at an HUC/D-negative cell. a, e, Arrowheads point at a spot-like Pleiotrophin pattern. a, Green signals indicate Pleiotrophin; red signals indicate PSA-NCAM. e, Green signals indicate GFAP; red signals indicate Pleiotrophin. In the RMS, the Pleiotrophin antibody stains small vessels (a, single small arrow pointing at green signals). a, Thick arrows indicate PSA-NCAM-positive cells in the RMS. f, Colocalization of FABP7 (red) and GFAP (green) in the SVZ and RMS. a-f, Blue signals indicate DAPI-stained nuclei. LV, Lateral ventricle; RMS, rostral migratory stream. Scale bars, 50 μm.
Pleiotrophin exerts many different cellular effects. In the context of our report, some are of special interest (Li et al., 1990; for review, see Deuel et al., 2002; Muramatsu 2002). It stimulates the migration of embryonic neurons, is present in radial glial processes, and localizes along neuronal migration routes (Matsumoto et al., 1994; Rauvala et al., 1994; Wewetzer et al., 1995; Nolo et al., 1996; Silos-Santiago et al., 1996; Kinnunen et al., 1998, 1999; Maeda and Noda, 1998). The migratory effect is mediated through the binding of Pleiotrophin to the transmembrane protein-tyrosine phosphatase PTPζ on migrating neurons (Maeda and Noda, 1998). Furthermore, it promotes neurite outgrowth, likely mediated through another receptor, Syndecan-3 (Raulo et al., 1994). Syndecan-3 is a transmembrane proteoglycan that binds to extracellular ligands and regulates signaling and cytoskeletal organization during developmental processes. It is localized to axons in the developing brain and to oligodendrocyte progenitors (Kinnunen et al., 1998; Hsueh and Sheng, 1999; Winkler et al., 2002) (for review, see Rapraeger, 2001).
Antibodies against FABP7 stain the neurosphere as well as the cell body and cellular processes of differentiating cells close to the cell cluster (Fig. 2). In the brain, FABP7 was found in the SVZ and along the RMS. The FABP signals coincided entirely with the GFAP staining, which strongly suggests that SVZ astrocytes synthesize FABP7 (Fig. 7f). Because the monoclonal GFAP antibody that we used is known to specifically recognize type B cells (neurogenic astrocytes) but not type C cells (transit-amplifying cells) (Doetsch et al., 2002), we conclude that FABP7 is expressed by type B cells.
FABPs are believed to play a role in the cellular solubilization, transport, and metabolism of their ligands, and the spatiotemporal production of mouse Fabp7 mRNA and protein correlates with differentiation of neuronal and glial cells (Feng et al., 1994; Kurtz et al., 1994). At embryonic day 10 (E10), it is expressed in the CNS columnar neuroepithelium, and from E12 to E14 becomes enriched in radial glial cell bodies at the ventricular surface, which extend FABP7-positive radial processes through the ventricular wall to the pial surface. This expression parallels the appearance of a broad layer of MAP2-positive neuronal cells in the intermediate zone of the ventricular wall. At later stages, FABP7 is no longer present in the ventricular zone; it is restricted to cells attached to the pial surface. This change in expression reflects a transition from radial glial cells to astrocytes. In the adult CNS, FABP7 has been detected in Bergmann glial cells of the cerebellum, in radial glial cells spanning the granule cell layer of the hippocampus, in the glia limitans, and in some astrocytes throughout the CNS. The application of anti-FABP7 antibodies to primary cell cultures led to an inhibition of glial and neuronal differentiation but did not show any effect on cell proliferation or adhesion (Feng et al., 1994).
Hb-Egf had transcriptional changes characteristic of cluster 6, with the strongest downregulation occurring at 48 hr. HB-EGF exists as a membrane-bound form (proHB-EGF) that is growth inhibitory. Through ectodomain shedding, a secreted form (sHB-EGF) is generated. The secreted form acts as a mitogen for multipotent neural progenitor cells in vitro as well as in vivo, particularly within the SVZ, and enhances neuronal survival (Kornblum et al., 1999; Caric et al., 2001; Jin et al., 2002, 2003). In situ hybridization studies on the expression of Hb-Egf in postnatal brain are contradictory. Nakagawa et al. (1998) found high levels within proliferative regions including the subventricular zone; Kornblum et al. (1999) did not. Caric et al. (2001) reported Hb-Egf expression in several migration pathways within the CNS.
Nell2 belongs to cluster 10, showing steadily increasing downregulation with time. It encodes a glycoprotein found in the rough endoplasmic reticulum and is secreted into the extracellular space (Kuroda et al., 1999). Immunocytochemical studies revealed its presence in brain, but not in peripheral tissues, and expression by neurons, but not by glial cells. NELL2 protein is abundant in the hippocampus and cerebral cortex, and moderate levels are found in the olfactory bulb (Oyasu et al., 2000). In the embryonic rat, Nell2 mRNA is detected in ventricular and differentiating zone of the spinal cord, medulla, and pons (Kim et al., 2002). In chick embryos, transcripts are found in the spinal cord ventricular zone and differentiating areas, whereas in the CNS Nell2 mRNA is restricted to postmitotic neurons as they start differentiation (Nelson et al., 2002).
Gene expression changes of chromatin-associated components and transcription factors
The change of an undifferentiated progenitor cell to a specialized cell type requires switches in gene expression that alter the pattern of inactive and active genes of these two cellular states. This switch in gene expression involves changes of chromatin states that differ in the accessibility of genomic DNA regions for transcription regulatory factors. In addition, the expression of individual transcription factors that positively or negatively control certain target genes is involved in this process. Strikingly, all downregulated genes in this category fall into clusters 6 and 10, with the exception of Fosb in the NT4 experiment. Both of these clusters represent the slightest gene expression changes of our study. Within cluster 6, genes decrease the most from time point 0 to 24 hr; genes of cluster 10 show a continuous downregulation over time. In parallel, most of the upregulated genes in this category are found in cluster 4, which describes slight alterations with the highest change between time point 0 and 24 hr. Exceptions are Id4 and Meiosis expressed gene 1; both fall into cluster 1, which contains the strongest upregulated genes.
Among the differentially expressed genes encoding chromatin and chromatin-associated proteins Histone 2Az (H2A.Z), Histone 2Ax (H2A.X), the high mobility group gene Hmgb2 (formerly Hmg2), and Thymopoietin (Lap2) were downregulated, whereas transcript levels of Histone 2B increased. The downregulation of Histone 2Az was confirmed by semiquantitative RT-PCR (Fig. 4).
Among the transcription factors with relevant expression changes were five that decreased in expression. Three of these, Transcription factor 12 (Htf4/Heb/Reb), Stra13 (bHlhb2/Clast5/Sharp2/Dec1), and Olig2, belong to the helix-loop-helix class. Furthermore, we found a downregulation of Etv5 (Erm) that belongs to the PEA3 group of transcription factors, Fosb, and Yb-1.
Increasing mRNA levels of eight transcription factors were identified. These include four genes encoding zinc finger proteins related to transcription: the RIKEN cDNA 3110024A21 gene, ENSMUSG00000013419, and genes encoding ZDHHC9, and DAN (neuroblastoma suppressor of tumorigenicity 1, N03). Furthermore, the bHLH-leucine zipper transcription factor Srebp-1, Id4, a negative regulator of HLH transcription factors, an EST similar to the human p66 α encoding gene, and the forkhead transcriptional activator Foxj1 (Hfh-4/Fkhl13/Hnf-3) were upregulated. The transcript level increase of Id4 was confirmed by semiquantitative RT-PCR. Finally, we found a distinct increase in the Meiosis expressed gene 1 expression. Its product associates specifically with meiotic chromosomes (Steiner et al., 1999); the relevance with respect to neural differentiation is unclear.
Our finding of an altered expression of H2A.X has to be looked at with caution, because the coding sequence of this gene differs only slightly from the 10 H2A1 and the single H2A2 variants, in contrast to H2A.Z, which differs considerably. Therefore, a cross-hybridization of H2A1, H2A2, and H2A.X transcripts can occur.
Hmgb2 is widely expressed during embryogenesis but exhibits tissue-restricted expression in the adult (Ronfani et al., 2001). Of note is its localization to the SVZ and intermediate zone at E17 but reported absence in the adult brain. Blockage of Hmgb2 expression with antisense oligonucleotides inhibits cell division (Yamazaki et al., 1995).
Htf4 belongs to the E-protein family of bHLH transcription factors that are widely expressed (Zhang et al., 1995) and decreases during neural differentiation (Klein et al., 1993; Neuman et al., 1993). Stra13 belongs to a distinct bHLH class similar to the Hes class. Its product acts as a transcriptional repressor at E-box sequence elements and promotes neuronal but inhibits mesodermal and endodermal differentiation (Boudjelal et al., 1997; Sun and Taneja, 2000; St. Pierre et al.,2002). Olig2 is a transcription factor that acts as a key mediator of oligodendrocyte specification. It is important for the self-renewal of cultured neurosphere cells and their differentiation into neurons and oligodendrocytes (Hack et al., 2004). Oligodendrocyte precursors, differentiated oligodendrocytes, and also neural progenitor cells express Olig2 (Lu et al., 2000; Takebayashi et al., 2000; Zhou et al., 2000; Tekki-Kessaris et al., 2001; Karsten et al., 2003).
Inhibitor of DNA binding proteins act through the formation of heterodimers with E-protein bHLH factors, thereby sequestering them from bHLH factors that require E-proteins such as members of the Mash, Neurogenin, NeuroD, and Olig family (Norton, 2000). In this way, they inhibit the differentiation of progenitor cells into neurons and oligodendrocytes. Of interest is also the elevated transcript level of an EST similar to the human p66 α encoding gene. The latter belongs to a new family of potent transcriptional repressors interacting with the methyl-CpG-binding proteins MBD2 and MBD3 (Brackertz et al., 2002).
In the context of our findings concerning genes for motor cilia and flagella components (see Gene expression changes related to the tubulin cytoskeleton), the slight upregulation of the forkhead transcriptional activator Foxj1 is of special interest. Foxj1 is expressed in tissues containing cells with motor cilia and flagella such as testis, oviduct, lung, and brain (Clevidence et al., 1993; Hackett et al., 1995; Lim et al., 1997; Murphy et al., 1997). It is present in motor cilia and flagella as well as in the unique motile monocilium of the embryonic node, but not in sensory cilia (Brody et al., 2000). The failure of directed left-right determination in Foxj1 (Hfh-4) null mice, which is normally regulated by embryonic node cells, indicates a Foxj1 (Hfh-4)-based defect of their motile monocilium (Brody et al., 2000).
Erm is expressed in various tissues at different developmental stages (Chotteau-Lelievre et al., 1997). Erm expression becomes restricted to cells lining the ventricles rather than differentiating cells of the neural wall as neural development proceeds.
Gene expression changes related to metabolism
Some of the genes with differential expression encode proteins involved in cellular metabolism. The downregulated genes were found in clusters 6 and 10 describing very slight expression changes. With the exception of Nat-1, none of the upregulated genes belongs to cluster 1, which contains strongly and very early increasing transcripts. Most of the genes of this category are related to glucose and energy metabolism, for example, Lactate dehydrogenase, Glucose-6-phosphate dehydrogenase, Phosphoglycerate dehydrogenase, and Brain glycogen phosphorylase. The latter is present in astrocytes and ependymal cells, but rarely found in neurons of the rat brain (Pfeiffer-Guglielmi et al., 2003).
Other differentially regulated genes identified are involved in amino acid transport (N system amino acids transporter Nat-1), carbon dioxide transport and acid-base balance (Carbonic anhydrase 2 and 14), and nucleic acid metabolism (Thymidine kinase 1, Dutpase). Thymidine kinase 1 is known to be present in proliferating cells with a peak in S phase. It is low in resting cells, which is in line with the decreased proliferation in our differentiation culture (see Gene expression changes related to cell proliferation).
An interesting gene that was downregulated is Phosphoglucose isomerase 1 (Gpi1). GPI1 plays a central role in both glycolysis and gluconeogenesis; however, it is also known as Autocrine motility factor or Neuroleukin and seems to fulfill several functions. Purified Neuroleukin from Bacillus stearothermophilus promotes neurite outgrowth of rat EGF-responsive neural embryonic progenitor cells (Sun et al., 1999).
Several genes coding for enzymes involved in steroid metabolism were also identified, all of them with increasing expression levels: UDP-glucuronosyltransferase 1, Oxysterol binding protein 2, Prostaglandin D2 synthase, the Leukotriene B4 ω-hydroxylase (CYP4F14), and Retinal short-chain dehydrogenase/reductase 2.
Gene expression changes related to ion channels, transporters, and lipoproteins
In addition, some genes encoding ion channels, transporters, and lipoproteins were differentially regulated, most of them upregulated. An exception was the Fxyd1 gene, which belongs to cluster 1 (strongly and very early increasing genes). It encodes a regulator of the Na,K-ATPase (Crambert and Geering, 2003). The GABA-A receptor α-5 subunit slightly increased, but only by 96 hr. Increasing transcripts of the genes Na, K transporting ATPase subunit α, the Abca3 transporter, and the Monocarboxylate transporter 1 were also detected. Of the genes related to lipoproteins, Lipid transfer protein II (Pltp) expression increases. PLTP is an important regulator of plasma high-density lipoprotein (HDL) levels and HDL particle distribution and is secreted by neurons, microglia, and astrocytes (Vuletic et al., 2003).
A slight and late increase is also noted for the gene encoding sorLA/LR11, a protein with features typical for low density lipoprotein receptor family and the Vps10p domain receptor family, which is present in neurons (Jacobsen et al., 1996; Motoi et al., 1999). On the other hand, transcript levels of Scavenger receptor class B1, encoding a lysosomal membrane protein that can bind HDL, decrease. Finally, gene expression of the Apolipoprotein J increases.
Discussion
The transition of multipotent neural progenitor cells into differentiated neuronal and glial cells involves complex gene expression changes of which only few are known to date.
To learn more about very early and probably decisive transcriptional changes, we compared the gene expression of undifferentiated neurosphere cells from the SVZ of P7 mice with cells at 24, 48, and 96 hr after induction of differentiation using BDNF or NT4. This time course allows us to ascertain the onset of transcriptional changes and follow subsequent alterations. Previously, neurosphere cells have only been compared with one cellular state after induction of differentiation using large scale expression methods. Geschwind et al. (2001) investigated undifferentiated neurosphere cells from P0 mouse cortex and cells differentiated for 24 hr using a combination of representational difference analysis and microarrays. Karsten et al. (2003) compared progenitor cells from prenatal day 1 and cells that were differentiated for 24 hr. Easterday et al. (2003) compared P1 cortical neuropheres differentiated for 24 hr with three different undifferentiated neurosphere populations from mice. Cells derived from the striatum of rat embryos were used in two other differentiation studies by Zhou et al. (2001) and Wen et al. (2002).
We have provided detailed information on many genes in Results. Our discussion therefore covers selected new and important aspects. Almost all of the differentially expressed genes showed consistent upregulation or downregulation throughout the time course. From an upregulation or downregulation at 24 or 48 hr, one can therefore usually assume a further increase or decrease at later time points for the vast majority of genes; however, we also identified genes that change only 96 hr after differentiation was initiated. Thus, by comparing undifferentiated cells with cells at only one early time point of differentiation, e.g., 24 hr, not all relevant transcriptional changes can be identified.
We describe 10 major patterns (clusters) of gene expression changes in the course of neural progenitor differentiation. The highest expression changes were found for downregulated genes. Of these, cluster 9 contains strong and very early expression changes, most of which involved genes with a role in cell proliferation; however, most clusters do not correspond to single functional gene categories.
A comparison of in vitro differentiation steps with biological processes in vivo is shown in Table 2.
Comparison of in vitro and in vivo findings
Molecular profile of neural progenitor cells
Because differentiated neuronal and glial cells were not investigated independently, we cannot correlate expression changes to cell types, except for the known cell type-specific genes; however, we can assign the downregulated genes to the undifferentiated neurosphere state, although spheres are also heterogeneous. This approach has been used successfully previously to identify a molecular profile of neural progenitor cells by comparing only two states (undifferentiated vs differentiated) (Geschwind et al., 2001; Easterday et al., 2003; Karsten et al., 2003). Several genes with higher expression in progenitor cells that were reported in one of these studies were also identified by us. They include the cell cycle-related genes Cyclin B1, D1, D2, Cdc2a, and PDZ-binding kinase. Other genes that decreased in expression at 24 hr in our study and the above-mentioned studies are Fabp7, Gpi, Hmgb2, Lrrfip1, Mcm6, Olig2, Pttg1, Tnc, Racgap1, Histone H2A.Z, and Cadherin 13 (T-cadherin).
Almost no overlap exists between genes that we identified and those found by Zhou et al. (2001) and Wen et al. (2002). These groups compared undifferentiated neurosphere cells from the striatum of rat embryos with differentiated cells. Only a DNA polymerase subunit was found by us and Zhou et al. (2001) and the Siva-1 (CD27) gene (not listed in Table 1) was found by us and by Wen et al. (2002).
In addition, many of the genes that we identified have been found previously to be specific to or enriched in neural stem-progenitor cells derived from embryonic (Ivanova et al., 2002; D'Amour and Gage, 2003; Fortunel et al., 2003) or adult (Ramalho-Santos et al., 2002) tissue by comparing their microarray profiles with those of other stem cell types (Table 1).
Many of the genes downregulated during differentiation are localized to the neural stem cell compartments. Thus, they are characteristic for neural progenitor cells in vivo and their expression is maintained under in vitro conditions. We detected the CD47, PTPNS1 (SHPS-1), Pleiotrophin, and FABP7 proteins in the SVZ and RMS of adult mouse brain sections by immunostaining. Among other cell types, all of them were found to be expressed by GFAP-positive cells in this region. In addition, some of our identified genes have been reported previously to be expressed within neural germinal zones. These are Tenascin C (Jankovski and Sotelo, 1996), Hb-Egf (Nakagawa et al., 1998), Nell2 (Kim et al., 2002), Transcription factor 12 (Uittenbogaard and Chiaramello, 2002), Cyclin D1, Cyclin D2 (Geschwind et al., 2001), and Lrrfip1, Hmg2b, Racgap1, Olig2, Fabp7, and Pttg1 (Karsten et al., 2003).
Similarities to radial glial cells
From our microarray data and immunostaining results we noticed that the undifferentiated neurosphere cells share certain markers with radial glial cells, namely FABP7, Tenascin C, Nestin, and Pleiotrophin. During embryonic brain development, radial glial cells appear in the ventricular zone and extend a long radial process to the pia (for review, see Parnavelas and Nadarajah, 2001; Campbell and Gotz, 2002; Tamamaki, 2002). They disappear from this region during the first 2 postnatal weeks (Tramontin et al., 2003). In contrast to songbirds, radial glial cells are absent from the intact adult mammalian SVZ (Alvarez-Buylla et al., 2002). Mammalian radial glial cells transform into mature astrocytes but also give rise to neurons in vitro and in vivo (Malatesta et al., 2000; Miyata et al., 2001; Noctor et al., 2001; Tamamaki et al., 2001; Heins et al., 2002). In the adult mammalian SVZ, astrocytes serve as stem cells. These cells and radial glial cells show an overlapping but not identical expression of molecular markers. Both express Nestin and Vimentin (Doetsch et al., 1997). In contrast to primates, radial glial cells of rodents do not express GFAP, with the exception of Bergmann radial glial cells of the cerebellum (Levitt and Rakic, 1980; Choi, 1981; Rakic, 1995; Sancho-Tello et al., 1995; Hatten, 1999), whereas SVZ adult stem cells (type B cells) are GFAP positive (Doetsch et al., 1997). FABP7, Tenascin C, Pleiotrophin, RC1, RC2, and GLAST are additional proteins synthesized by radial glial cells (Malatesta et al., 2003). Tenascin C is also present in the adult SVZ and at elongated astrocytes in the RMS (Jankovski and Sotelo, 1996). This localization strongly indicates that adult neural progenitor cells synthesize Tenascin C.
It has recently been suggested that there is a continuum of a neural stem cell compartment that links radial glial cells to astrocytes, which serve as neural stem cells in the adult mammalian brain (Alvarez-Buylla et al., 2001; Doetsch, 2003; Goldman, 2003; Gotz, 2003; Tramontin et al., 2003). The expression of GFAP, Nestin, FABP7, Tenascin C, and Pleiotrophin by our neurosphere cells places them between radial glial cells and SVZ astrocytes and could be explained by a transitional state from radial glial cells to SVZ astrocytes, because they were isolated at P7. Of interest in this context is our finding of FABP7 expression by SVZ and RMS astrocytes and of Pleiotrophin expression by different cell types in the SVZ, apparently also by SVZ astrocytes, of adult mice (Fig. 7), indicating that not only radial glial cells and possibly cells in a transitional state but also neural progenitor cells in the adult brain share these markers. Although one has to interpret the presence of molecular markers with caution because they might specify a cellular state rather than lineage, our results are consistent with the proposed lineage relationship of radial glial cells and adult SVZ stem cells.
Cellular interactions within the SVZ and RMS
It is essentially unknown how the various cells in the SVZ stem cell niche and the RMS communicate through cell-cell contacts or locally secreted factors. We found that undifferentiated neurosphere cells express Hb-Egf, and the secreted form of the respective protein stimulates neural progenitor cell proliferation and enhances neuronal survival (Kornblum et al., 1999). Hb-Egf has also been detected in several migration pathways (Caric et al., 2001). HB-EGF could therefore act as an autocrine factor on the neural progenitor cells themselves and at the same time support the survival of migrating neuroblasts. We also found that the Pleiotrophin gene is expressed by undifferentiated neurosphere cells. Moreover, Pleiotrophin appears to be expressed by different cell types, including astrocytes and neuronal cells, in the SVZ. Several features of this protein in the embryonic brain correspond well to the situation in the adult SVZ-RMS compartment. Pleiotrophin is known to stimulate the migration of embryonic neurons in vitro, and it has been localized to radial glial processes along which embryonic neurons migrate (Maeda and Noda, 1998). We conclude that Pleiotrophin is secreted by neural progenitor cells, which as outlined above, share similarities with radial glial cells.
Finally, we found that the transmembrane proteins PTPNS1 (SHPS-1) and CD47, which can bind each other through their extracellular domains, are both produced by undifferentiated neurosphere cells. Moreover, an enrichment of these proteins at the interface of SVZ cells and the underlying extracellular space was noticed. In addition, we localized CD47 and PTPNS1 to GFAP-positive astrocytes and neuronal cells in the SVZ and RMS. In view of the very tight cellular contacts in this brain region, a clear assignment to the membrane of specific cell types, however, must await analysis by electron microscopy. SHPS-1 and CD47 are a cell-cell communication and signal transduction system with multiple functions in immune cells (Han et al., 2000; Oldenborg et al., 2000; Seiffert et al., 2001). In the brain, SHPS-1 and CD47 are associated with synapses, and SHPS-1 (P84) has been shown to promote neurite outgrowth (Chuang and Lagenaur, 1990; Mi et al., 2000). Two aspects of the CD47-SHPS-1 system are especially interesting. First, they are involved in IGF-I receptor signaling, and we found striking expression changes of genes encoding IGF-I modulating proteins (IGFBPs), which indicate that very early in differentiation neurosphere cells are converted into an IGF-responsive state. Second, the CD47-SHPS-1 system regulates cell aggregation and migration. This function has not been described for neural cells so far, but has been described for immune and melanoma cells (Babic et al., 2000; Liu et al., 2002; Motegi et al., 2003). The reported localization in tightly associated neurosphere cells, cells migrating as chains in vitro, and in the SVZ and RMS would be compatible with a role in cell-cell/cell-matrix contact and migration. Further investigations on the CD47-SHPS-1-system are underway in a separate study.
In summary, we have provided a dataset of dynamic gene expression changes during the first 96 hr of postnatal neural progenitor differentiation. In addition to the investigation of transcript levels, we could correlate selected proteins to the cells that produce them. Our analysis has identified many new candidate molecules and confirmed several others that likely play important roles during the maintenance, differentiation, and migration of precursors in the SVZ. For several, the functional profile in other cellular systems strongly suggests a plausible role also in the SVZ; for others we have shown an association with cell types in the SVZ. Thus, our data should serve as a basis for future analyses of the molecular and cellular dissection of SVZ function in adult neurogenesis.
Footnotes
This work was supported by Deutsche Forschungsgemeinschaft Grant SFB577. We thank Andrew Fanous for technical assistance and Isabelle Caille for support in establishing neurosphere cultures.
Correspondence should be addressed to Ulrike A. Nuber, Max Planck Institute for Molecular Genetics, Human Molecular Genetics, Ihnestrasse 73, Berlin 14195, Germany. E-mail: nuber{at}molgen.mpg.de.
Copyright © 2004 Society for Neuroscience 0270-6474/04/245982-21$15.00/0
References
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵