|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
CHARACTERIZATION SERIES |
a Gerontology Research Center; Stem Cell Biology Unit, Laboratory of Neurosciences, National Institute on Aging, National Institutes of Health, Department of Health and Human Services, Baltimore, Maryland, USA;
b Laboratory of Molecular Neurobiology, Medical Biochemistry and Biophysics, Retzius Laboratory, Karolinska Institutet, Stockholm, Sweden;
c Buck Institute for Age Research, Novato, California, USA;
d SuperArray Bioscience Corporation, Frederick, Maryland, USA;
e Department of Neuroscience, Johns Hopkins University, Baltimore, Maryland, USA;
f Invitrogen Corporation, Carlsbad, California, USA
Key Words. Human embryonic stem cells • Differentiation • Oligodendrocytes • Astrocytes • Dopamienrgic neurons
Correspondence: Mahendra S. Rao, M.D., Ph.D., Invitrogen Corporation, 1600 Faraday Road, Carlsbad, California 92008, USA. Telephone: 240-344-1781; Fax: 410-558-8249; e-mail: mahendra.rao{at}invitrogen.com
Received August 15, 2005;
accepted for publication December 5, 2005.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Multiple classes of cells have been considered for cell therapy, including neural stem cells (NSCs), glial restricted precursor cells, embryonic stem cells (ESCs), mesenchymal stem cells, and transdifferentiated cells [7, 8]. In each case, questions concerning the character of the transplanted population, signals directing differentiation, and specificity of differentiation have been raised. Addressing these issues has been difficult due to the lack of consensus over surrogate markers of efficacy and the availability and cost associated with assessing a reasonable number of markers. Multiple analytical techniques to assess gene expression in defined cell types have been developed, including microarray analysis, EST (expressed sequence tag) enumeration, SAGE (serial analysis of gene expression), and MPSS (massively parallel signature sequencing). Focused microarray analysis offers an advantage by allowing one to develop a customized array containing genes specific to particular cell populations or signaling pathways. Recently, focused microarrays have been designed to assess the state of ESC differentiation and characteristics of the undifferentiated ESC state [9, 10]. Application of these arrays has allowed researchers to evaluate the state of cells in a relatively inexpensive, rapid, and reliable way. We reasoned that a focused array encompassing a substantial, but not overwhelming, number of genetic markers for oligodendrocytes and dopaminergic neurons would allow rapid assessment of the state of a cell population prior to their use in cell replacement therapy and would be as good as other, more expensive large-scale array methods.
We have developed a focused array of approximately 280 genes that includes markers for dopaminergic neurons, glia, neural progenitors, pluripotent cells, and signaling molecules thought to be important in regulating neural differentiation. Additionally, we included cytokines, chemokines, and their receptors to provide a global view of the potential signaling pathways that may regulate the process of differentiation. Using human adult brain substantia nigra, human ESC (hESC)derived NSCs, and pluripotent hESCs, we show that the array is capable of distinguishing these populations by their gene expression profiles. Moreover, we have used this array to monitor NTera2 differentiation into dopaminergic neurons. We show that such an array can monitor the process of differentiation and provide important clues to factors that regulate the differentiation process.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Cultures of hESCs and NTera2 Cells
hESC lines BG01 and BG03 were obtained from BresaGen (Thebarton, SA, Australia) and grown using published protocols described elsewhere [11]. Briefly, hESCs (20,000 cells/cm2) were cultured on inactivated mouse embryonic fibroblast (MEF) feeder cells in Dulbeccos modified Eagles medium/Hams F-12 (DMEM/F-12; 1:1) supplemented with 20% knockout serum replacement, 2 mM nonessential amino acids, 2 mM L-glutamine, 50 µg/ml Penn-Strep (all from Invitrogen, Carlsbad, CA, http://www.invitrogen.com), 0.1 mM ß-mercaptoethanol (Specialty Media, Phillipsburg, NJ, http://www.specialtymedia.com), and 4 ng/ml of basic fibroblast growth factor (bFGF; Sigma, St. Louis, http://www.sigmaaldrich.com). After 45 days, cells were then dissociated by using Cell Dissociation Buffer (Invitrogen) and passaged on mitotically inactivated MEF at a density of 20,000 cells per cm2. The culture medium was changed every other day. Under these conditions, these cells remained Oct3/4- and SSEA4-positive undifferentiated stem cell status [11].
The NTera2 cell line was obtained from American Type Culture Collection (Manassas, VA, http://www.atcc.org) and cultured according to the manufacturers protocol with minor modifications. Briefly, cells were plated at a density of 10,00015,000 cells per cm2 at 37°C in 5% CO2 in DMEM supplemented with 10% fetal bovine serum (FBS) and 50 µg/ml penicillin/streptomycin (all from Invitrogen). The media was changed every 23 days until cells become confluent. The cells were either passaged using a Costar cell scraper (Krackler Scientific Inc., Albany, NY, http://www.krackeler.com) or used for RNA extraction.
Derivation and Cultures of NSC
NSCs used were derived from the BG03 hESC line. Briefly, after 1 week of culture on mouse feeder cell layer, hESCs were fed every other day with derivation medium containing DMEM/ F-12 medium (1:1; Gibco, Grand Island, NY, http://www.invitrogen.com) supplemented with 2 mM L-glutamine, 50 U/ml penicillin, 50 µg/ml streptomycin, 1X N2 (Gibco), and 4 ng/ml of basic fibroblast growth factor (Sigma) for 7 days. The mouse feeder layer was then removed physically using a pipette, allowing neuroepithelial precursor to attach to the culture dish and develop rosettes after 3 days in derivation medium [12]. Derived NSCs were propagated further on 20 µg/ml of poly-ornithine and 1 µg/ml of laminin-coated dishes in neurobasal medium (Gibco) supplemented with 2 mM L glutamine, 50 U/ml penicillin, 50 µg/ml streptomycin, 1 X B27 (Gibco), 20 ng/ml bFGF (Sigma), and 10 ng/ml of leukemia inhibitory factor (Chemicon, Temecula, CA, http://www.chemicon.com). The culture medium was changed every other day, and cells were then passaged by mechanical triturating or by trypsin when confluent. The NSCs used in this paper had undergone four passages and were uniformly Oct3/4- and SSEA4-negative and Nestin- and Sox2-positive.
Preparation of PA6-Induced PSA-NCAMSorted Differentiated NTera2 Cells
Neural differentiation of NTera2 cells was induced using the mouse stromal cell line PA6 as described by Schwartz et al. with minor modifications [13]. Briefly, PA6 cells were grown and maintained in minimum essential medium-
supplemented with 10% FBS and 50 U/ml penicillin/streptomycin (all from Invitrogen) until confluent. NTera2 cells were then seeded onto this confluent monolayer of PA6 feeder cells at a density of 2,000 cells per cm2. The day at which NTera2 cells were seeded onto PA6 cells was taken to be day 1. After 12 days in culture, these NTera2-PA6 co-cultures were used to isolate polysialic acid (PSA)-NCAMpositive cells by flow cytometry sorting as previously described [13]. These PSA-NCAMpositive cells were then used for total RNA extraction and subsequent analysis.
RNA Preparation, Reverse TranscriptionPolymerase Chain Reaction Amplification, and Focused Oligo-Microarray Analysis
Total RNAs from indicated cells were isolated using RNA STAT-60 (Tel-Test Inc., Friendswood, TX, http://www.isotexdiagnostics.com) by following the manufacturers instructions. Total RNA derived from adult human substantia nigra was purchased from Clontech (Palo Alto, CA, http://www.clontech.com). The cDNAs were synthesized using a Superscript II reverse transcriptase kit with 1 µg total RNA and 500 ng oligo(dT)1218 (both from Invitrogen). The polymerase chain reaction (PCR) was performed in a 20 µl reaction solution containing 2 µl 10 x PCR buffer, 150 nmol MgCl2, 10 nmol dNTP, 20 pmol primer, 1 µl 50x diluted cDNA, and 1 U RedTag DNA polymerase (Sigma). The PCRs were run as follows: 35 cycles of 94°C for 30 seconds, 55°C for 30 seconds, and 72°C for 30 seconds, and final extension for 10 minutes at 72°C. Primer sequences are shown in Table 2
.
|
For data analysis, the positive and negative spots were independently identified and verified by at least two people. Only matched positive and negative results from two independent experiments were used for analysis. For quantification, intensity of signal was first measured using ImageQuant 5.2 software (Molecular Dynamics, Sunnyvale, CA, http://www.mdy.com) with a local background subtraction method. These subtracted intensities were then divided by the average of intensities from glyceraldehyde-3-phosphate dehydrogenase (GAPD) (three spots in each array) to obtain a relative intensity for each spot. These relative intensities were used to calculate correlation value (R2) related to input amounts of cRNAs or RNAs or to compare samples.
| RESULTS |
|---|
|
|
|---|
The gene specific 60-mer 3' biased oligonucleotides for these 281 genes were designed and arrayed in a 25 x 12 format. These oligo probes were grouped and printed based on their function (Fig. 1A
). Housekeeping genes, such as GAPD, ß-actin, RPS27A, and B2M, served as positive controls. GAPD was spotted in triplicate in the upper right corner, and hybridization intensities were used as a reference to normalize data. RPS27A and B2M, known to be expressed high and low, respectively, were spotted on the top and bottom of the left side and served as additional positive hybridization controls. Biotinylated artificial sequence 2 complementary (BAS2C) sequences were spotted at an increasing gradient on the right corner to serve as a positive detection control. Negative controls, such as a blank, plasmid PUC18 DNA, or artificial sequences not expected to be present in cDNA, were also included as detection and background hybridization controls. This array format allows us to easily monitor quality of array hybridization.
|
Validation Testing
Next, we examined the ability of focused array to selectively detect genes highly expressed in dopaminergic neurons, glial, and embryonic stem, and neural progenitors by using cell-specific samples. We chose to use adult human substantia nigra as a source of material that would include dopaminergic neurons and glial cells, and the hESC line BG01V and NSCs derived from the hESC line BG03 as samples that would contain pluripotent and progenitor cells.
The patterns of gene expression profiles easily distinguished these three cell types (Fig. 2A
). In the human substantia nigra sample, the array detected high gene expression of dopaminergic and glial markers with no or low expression of pluripotent markers. In undifferentiated BG01V and NSCs, the array showed high expression levels of pluripotent markers and little expression of dopaminergic and glial markers. Correlation analysis among these samples with normalized intensities indicated that the adult substantia nigra sample poorly correlated with BG01V (R2 = 0.50) and NSC (R2 = 0.59) samples, whereas BG01V correlated well with NSC (R2 = 0.81) samples. We detected nine dopaminergic neuron markers and 18 glial markers in the adult substantia nigra samples; however, few of these markers were detected in the BG01V and NSC samples. For example, tyrosine hydroxylase (TH), dopa decarboxylase (DDC), and dopamine transport (DAT), enzymes required for either synthesis or transport of the neurotransmitter dopamine in dopaminergic neurons, were highly detected in the adult substantia nigra (Fig. 2A
, arrows) but not observed in BG01V and NSC samples. Likewise, glial fibrillary acidic protein (GFAP), S100-ß, myelin basic protein, and oligodendrocyte transcription factor 1 and 2 (Olig1 and Olig2), all known glial markers for either astrocytes or oligodendrocytes, were also detected in the adult substantia nigra but not observed in BG01V and NSC samples. Pluripotent markers, including Nanog and Oct3/4, were detected in BG01V samples and not detected in the adult substantia nigra samples. Sox2 and nestin, markers typically expressed in NSCs or progenitors, were highly expressed in NSC samples but low or absent in the adult substantia nigra. These tissue-specific genes detected in the array are summarized in Figure 2B
. Additionally, some of these tissue-specific genes were used to perform RT-PCR analysis to confirm array data (Fig. 2C
). Thus, this array can distinguish dopaminergic neurons and glial cells from ESC and NSC populations.
|
Although numerous markers and signaling molecules were detected in a tissue-specific fashion, we noted that some genes, which are known to be expressed in these tissues, were not detected by this array. For example, LMX1B, a transcriptional factor promoting dopaminergic neuron development, and Dpp5a, an hESC-specific pluripotent marker, were both detected by RT-PCR (Fig. 2C
) but not detected by the array. This is likely caused by low sensitivity of the array compared with RT-PCR and highlights the importance of direct testing and redesigning of these probes further to improve array quality.
Titration Testing
Given the ability of the array to distinguish dopaminergic neurons, glia, and hESCs, we next performed titration experiments to determine whether the array could detect tissue-specific genes in a concentration-dependent manner. These experiments included titration of both the hybridization and labeling processes. For titration experiments in testing hybridization, we first generated biotionylated cRNA targets from human substantia nigra and/or hESC line BG01V. Then, we held the total cRNA (5 µg/array) constant and altered the ratio of cRNA from human substantia nigra to hESC BG01V as follows: 1:0, 0.8:0.2, 0.5: 0.5, 0.2:0.8, and 0:1. The hybridizations were performed, and intensity data were normalized to the housekeeping gene GAPD. The results showed a positive linear relationship of detected hybridization signals of the above tissue-specific genes to inputs of cRNAs with R2 values from 0.6 to 0.99 (RSQHyb in Fig. 3C
), indicating that the hybridization condition was optimal. Next, we tested whether the amplification labeling system was linear to total RNA input (2 µg/labeling) by using the same ratio of human substantia nigra to human BG01V samples as those in hybridization experiments. Visually, intensity changes of specific genes are correlated to their input amounts of RNA, such as DAT1, TH, GFAP, and Oct3/4 (Fig. 3A
). As shown in Figure 3B
, there was a high correlation between detected signals of cell typespecific genes and RNA inputs. The R2 values and slopes obtained from both titration experiments in hybridization alone and hybridization plus labeling process are summarized in Figure 3C
. The results showed that the labeling and hybridization system was optimal to detect these cell-specific markers in a dose-dependent manner. Considering these data, we suggest that these dose-response tissue-specific genes may serve to develop a reporting list of candidate differentially expressed genes for detection of dopaminergic neurons and glial cells in mixed populations.
|
The array results indicated that the sorted PSA-NCAM+ NTera2 cells expressed markers for dopaminergic neurons. The hybridization image profile of the sorted PSA-NCAM+ NTera2 cells (Fig. 4A
) showed more hybridized spots in the dopaminergic marker group (first five rows of array) than did those of undifferentiated NTera2 and FA2B5+ populations. Dopaminergic markers, such as TH and DDC, showed higher expression in the PSA-NCAM+ NTera2 samples than in either the undifferentiated NTera2 or the FA2B5+ samples (Fig. 4B
). Engrailed 1, a transcriptional factor, and Smoh, a receptor for Shh, were also expressed in the PSA-NCAM+ NTera2 sample. RT-PCR and immunostaining results further verified the array results (Fig. 4C, 4D
, first row).
|
| DISCUSSION |
|---|
|
|
|---|
Several important parameters were assessed to ensure that the produced array be reliable and reproducible. We used sophisticated selection criteria to identify oligonucleotides of the same length which showed minimal cross-hybridization and had roughly the same annealing temperatures, allowing for an increased sensitivity and reduced background. We tested the uniformity of loading and arraying using labeled oligonucleotides and measuring overall spot intensity. Additionally, hybridization with single probes showed that each probe identified its cognate partner and not another gene on the array (data not shown). Loading controls at three different concentrations measured saturation and provided a rough measure of the relative level of gene expression. For additional control purposes, we included random oligomers as negative controls and included oligonucleotides on the array, such that spiking with probes would allow one to assess the quality of labeling and hybridization. These controls serve to provide immediate feedback on the hybridization experiment and on whether one can compare between two independent experiments.
It is useful to emphasize the advantage of being able to format the array such that visual information can be provided even before a detailed analysis is performed. We have previously shown that results are identical when probes are randomly placed or organized into useful subsets [9]. In these experiments, we separated potential markers for pluripotent and progenitor cells, oligodendrocytes, astrocytes, and dopaminergic neurons and placed cytokines and their receptors in separate groups (see Results). Visually, it was clear simply from the pattern of gene expression which cell type was being analyzed. This provides a quick feedback prior to a detailed quantitative analysis and the relative levels of expression (compared with the controls spotted at different concentrations) and allows one to more carefully adjust the scanning parameters for more quantitative studies.
We used human substanatia niagra samples to assess the sensitivity and reliability of our arrays. Human substanatia niagra was chosen because it contains a subset of dopaminergic neurons in a mixed population and thus represents a common situation in the laboratory. Examining gene expression identified markers, such as TH, DDC, and DAT (see Results), in which the magnitude of change was sufficient to be readily assessed. These results were further confirmed by RT-PCR using standardized primer sets we developed for all genes present on the array. Testing a variety of samples showed that this focused array could reliably distinguish dopaminergic neurons and glial cells from ESCs and NSCs. The lack of any mismatch of cell typespecific markers by the array indicated the absence of problems in probe design, printing contamination, or cross-hybridization, and these were key criteria in establishing the validity of the focused array. Moreover, our titration experiments, in which we gradually mixed the substantia nigra sample with the hESC sample, showed a linear relationship of hybridization intensity with input RNA, further supporting the reliability of this array. The titration experiments were very helpful in identifying genes that distinguish one cell type from another on the focused array.
Although we very carefully chose a list of genes for dopaminergic neurons and glial cells based on the published literature, not all of genes were differentially expressed at levels that could be detected by the array. This is likely due to their expression levels and sensitivity of the array hybridization process as their expression could be confirmed by RT-PCR (data not shown). For example, LMX1B, a transcription factor required for dopaminergic neuronal development, was detected by RT-PCR in human substantia nigra; however, the array could not detect LMX1B in the human substania niagra sample. Although one could redesign the probes and further optimize hybridization, it would be impractical to regenerate new lists with increased sensitivity given that this array was sufficiently sensitive to distinguish populations and monitor dopaminergic differentiation. The responses of reporting genes in titration experiments were concentration-dependent in both the hybridization and labeling processes with high linear correlation efficiency (0.71). The slopes of their response to inputs reflected the sensitivities of detection in this array. However, we generated a more restricted list of genes that shows the largest range of difference; we suggest that these be used to develop a more sensitive quantitative PCR assay if such an assay is required.
In addition, we included an additional 118 genes that are known to play essential roles in neural development which could assist in dissecting molecular events occurring during neuronal development. These include signaling molecules for Wnt-Fzd, TGF-ß, Notch, fibroblast growth factor (FGF), and BMP (bone morphogenetic protein) pathways. By comparison of gene expression profiling in signaling groups between undifferentiated NTera2 cells and PA6-induced sorted cells, we have noted that some known and required pathways for induction and formation of dopaminergic neurons were activated in the NCAM+ NTera2 sorted population. We observed enhanced expression of Smoh, Nr4a2, En1, GFRA2, and Fzds in the sorted differentiating cells that have previously been reported as signaling pathways in dopaminergic differentiation. Shh-Smoh activation and FGF8 signaling are known to be key players in midbrain pattering and genesis of dopaminergic neurons [15, 16]. Nr4a2 is a transcription factor and is expressed in both dopaminergic precursors and neurons in ventral midbrain, and deletion of Nr4a2 results in a loss of dopaminergic neurons in ventral midbrain [1719]. Engrailed genes (En1 and En2) were shown to be involved in dopaminergic neuron survival and maintenance [20]. Activation of Wnt-Fzds pathways has multiple functions, including promoting proliferation of NSCs and dopaminergic precursors and differentiation from dopaminergic precursors to their mature neurons depending on the members of Wnts involved [2123]. Although further dissection of signaling pathways involved in promoting dopaminergic formation of NTera2 cells is required, our results suggest that similar pathways are activated in PA6-induced dopaminergic differentiation and in midbrain dopaminergic neuron formation during development.
Many have debated the utility of focused arrays versus a global array, which contains all genes present on a focused array and provides substantial additional information. In the past, we have argued that most large-scale arrays are not complete and often do not contain the genes that are of immediate relevance. Although still true, newer whole genome arrays from providers such a NimbleGen Systems, Inc. (Madison, WI, http://www.nimblegen.com), Illumina (San Diego, http://www.illumina.com), and Affymetrix (Santa Clara, CA, http://www.affymetrix.com) are closer to ensuring the presence of all genes of interest. Similarly, as technology has advanced, global arrays have become less expensive than previous arrays. A large array, although 10 times more expensive ($500$600 [USD]), incurs costs in reagents and personnel time similar to a focused array and provides more information. Nevertheless, we feel that a focused array offers several advantages. A focused array is still significantly cheaper to develop and run, and it can be run easily by any laboratory. Also, the data analysis and reporting issues are easily manageable without elaborate bioinformatics support. It is potentially useful for some screening when a number of samples have to be treated for a short time. Posting and comparing data across laboratories is also relatively simple given that criteria of sensitivity, cutoffs, and hybridization efficiency are much easier to address for a few genes than for larger numbers of genes. The sensitivity of assessment using selected genes is often higher than in global arrays, in which statistical methods and large numbers of replicates have to be run to extract meaningful information [24, 25]. Adding new genes (e.g., markers for other types of neurons, such as GABA-type), rearranging formats, and adding additional controls is trivial in a focused array, which allows the maintenance of continuity with old data sets.
In summary, our goal was to develop a list of genes that could be used effectively in a focused array format for routine assessment of the process of differentiation. We desired an inexpensive method that was robust and reliable and could be used on a routine basis to monitor differentiation. Our results showed that a focused array fulfills these criteria and permits the monitoring of at least 300 genes (an order of magnitude larger than before) as cells differentiate. Careful selection of the genes has permitted one to distinguish stages of differentiation and degree of contamination of undifferentiated cells and identify the crucial signaling pathways that direct the process of differentiation. As data accumulate with different populations and different methods of differentiation, one will perhaps be able to identify the key regulators and biomarkers that may allow further reduction of the number of genes needed to monitor specific populations of neural derivatives.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| DISCLOSURES |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. M.-S. Tay, W.-L. Tam, Y.-S. Ang, P. M. Gaughwin, H. Yang, W. Wang, R. Liu, J. George, H.-H. Ng, R. J. Perera, et al. MicroRNA-134 Modulates the Differentiation of Mouse Embryonic Stem Cells, Where It Causes Post-Transcriptional Attenuation of Nanog and LRH1 Stem Cells, January 1, 2008; 26(1): 17 - 29. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Rao One Successful Series Begets Another Stem Cells, October 1, 2006; 24(10): 2160 - 2161. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| STEM CELLS | THE ONCOLOGIST | CME | ALPHAMED PRESS JOURNALS |
