|
|
||||||||
Research Institute of Biotechnology, Histostem Co., Kangdong-Gu, Seoul, Korea
Key Words. Mesenchymal stem cells • Human umbilical cord blood • DNA microarray
Correspondence: Hoeon Kim, Ph.D., Research Institute of Biotechnology, Histostem Co. 518-4 Taijul Bldg., Doonchun-dong, Kangdong-gu, Seoul 134-060, Korea. Telephone: 82-2-488-8154; Fax: 82-2-470-6342; e-mail: hoeonkim{at}seoulcord.co.kr
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Self-renewal is one of the most fundamental properties of stem cells, and the cellular and molecular mechanisms underlying it have been a subject of extensive studies. Elaborated searches for key transcriptional players led to discoveries of Oct3/4 [3], leukemia inhibitory factoractivated Stat3 [46], and Nanog [7, 8] in ESCs, and a recent finding of Bmi-1 in both hematopoietic stem cells (HSCs) [9] and neural stem cells (NSCs) [10]. Also in efforts to envisage a molecular entity of stem cells in the genomic scale, the DNA microarray-based analyses were recently employed in human or mouse ESCs, HSCs, and NSCs, leading to identification of commonly expressed genes, called stemness genes or a stem cell molecular signature [1114].
As in any other stem cells, the self-renewal of MSCs is likely to be operated by a defined set of molecular factors, but no molecular factor relevant to this function has been identified to date. The gene expression profile of MSCs has been previously investigated through serial analysis of gene expression (SAGE) [1517] and restriction fragment differential display [18]. Although the studies provided us with a plausible framework to define the MSCs in the genetic level, the presence of abundant housekeeping genes prevented the correct assessment of MSC-specific genetic messages.
In the study reported here, with the specific aim to generate an MSC-specific transcriptome, we performed a DNA microarray-based differential gene expression analysis between a fraction of human umbilical cord blood (UCB)derived mononuclear cells (MNCs) and its MSC subpopulation. UCB-derived cells were proven to be more advantageous in cell procurement, storage, and transplantation than their bone marrow (BM) counterpart and therefore better suited in tissue engineering and development of cell-based therapeutics. A number of reports from different laboratories [1923], including ours [2426], indicate that UCB-derived MSCs are highly similar to the cells of BM origin with respect to cell characteristics and multilineage differentiation potential. Therefore, this study may lead us to reveal the molecular signature that is specific to human MSCs but independent of their origins, and it will assist further studies on molecular mechanisms controlling various core stem cell properties.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Microarrays and Target Preparation
The Amersham CodeLink system (Amersham Biosciences, Chandler, AZ, http://www.amersham.com) using the UniSet Human 20K Bioarray, containing approximately 20,289 probes on a single glass slide, was used to generate the gene expression profile of the cells. Microarrays were processed in parallel using the CodeLink Shaker Kit and CodeLink Parallel Processing Kit. For each microarray, double-stranded cDNA and subsequent cRNA were synthesized from 1 µg of total RNA using the CodeLink Expression Assay Kit, according to the manufacturers instructions. Briefly, first-strand cDNA was generated using SuperScript II reverse transcriptase and a T7-oligo(dT) primer. Subsequently, second-strand cDNA was produced using Escherichia coli DNA polymerase 1 and RNase H. The resultant double-stranded cDNA was purified on a QIAquick column (Qiagen), and cRNA was generated via an in vitro transcription reaction using T7 RNA polymerase and biotin-11-uridine-5'-triphosphate, tetralithium salt (Perkin-Elmer, Boston, http://www.perkinelmer.com), then purified on an RNeasy column and quantified by UV spectrophotometry. A total of 10 µg of cRNA was subjected to fragmentation by heating at 94°C for 20 minutes in the presence of Mg ions.
Hybridization, Processing, and Scanning
The fragmented cRNA in 260 ml of hybridization solution was added to each slide and incubated overnight at 37°C in a shaking incubator (Vision Scientific Co., Kyunggi-do, Korea, http://www.visionsci.co.kr) at 300 rpm. After hybridization, the slides were washed in 0.75x TNT buffer (1x TNT: 0.1 mol/L Tris-HCl, pH 7.6, 0.15 mol/L NaCl, and 0.05% Tween 20) at 46°C for 1 hour, followed by incubation with Cy5-streptavidin at room temperature for 30 minutes in the dark. Slides were then washed in 1x TNT twice for 5 minutes each, followed by a rinse in 0.05% Tween 20 in water. The slides were then dried by centrifugation and kept in the dark until scanning. Images were captured on an Axon GenePix Scanner (Molecular Devices Co., Union City, CA, http://www.axon.com).
Microarray Data Processing and Hierarchical Clustering
Scanned data images were processed using CodeLink Expression Analysis Software. The mean intensity is taken for each spot and background corrected by subtracting the surrounding median local background intensity. The intensities were global linearly normalized according to the standard normalization procedure of the software. The normalized intensity of each gene probe was separately averaged over the MNC and MSC populations, and the probes were ranked by the MSC-to-MNC ratio of the average intensity. The full list of the gene probes whose average intensity ratios are greater than 1.5 is available, along with their respective normalized intensity in each sample, at http://callisto.snu.ac.kr/hoeonkim/microarray. For hierarchical clustering, the raw intensity data were exported to GeneSpring software version 6.0 (Silicongenetics, Redwood, CA, http://www.silicongenetics.com). Gene expression data were normalized to the 50th percentile expression level. Rigorous filtration of flagged probes resulted in a total of 11,662 gene probes, which were then subjected to hierarchical clustering using the standard correlation as a similarity measure.
RT-PCR Confirmation of Microarray Data
To confirm the gene expression profile determined by microarrays, a number of select genes were subjected to RT-PCR analysis, using total RNAs derived from the four cell samples used for a microarray experiment, as well as an additional pair of fresh MNC and MSC samples. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) mRNA amplified from the same samples served as an internal control. The sense and antisense primers used for each gene were as follows:
CTGF: sense5'-TTCACTTGCCACAAGCTGTC-3'and antisense5'-CTGCACCAGCATGAAGACAT-3'
COL1A1: sense5'-CCATGTGAAATTGTCTCCCA-3'and antisense5'-GGGGCAAGACAGTGATTGAA-3'
NNMT: sense5'-GGTCAAAGGAATTGCTTTAATTG-3'and antisense5'-CTACTACATGATTGGTGAGCAGA-3'
TAGLN: sense6'-AGAGGGCTAGCCCTGAGC-3'and antisense5'-CTTCTTCAATGGGCTTTTGC-3'
SERPINH2: sense5'-CAACTATAAAACTAGGTGCTG-3'and antisense5'-CCCTATCTTTGAAGAACTAAA-3'
COL6A3: sense5'-AAGCCTGTGTATCGTGGAG-3'and antisense5'-ATTAAGGCATTGGTCCCAAC-3'
OSF-2: sense5'-AGTTCTGGCTAACTTTGGAATCC-3'and antisense5'-TAAAAAATATGCATTGCAAGAAGC-3'
COL1A2: 5'-CAGGTCCTTGGAAACCTTGA-3'and antisense5'-TTTACAAGAGGAAACTGTAAGA-3'
SERPINE1: sense5'-AGATCAGCACCACAGACGC-3'and antisense5'-TTTTGTGTGTGTCTTCACCCA-3'
COL3A1: sense5'-GACTTCCAAGACCTCCTCTTT-3'and antisense5'-CCACAAGGATTACAAGGCTTG-3'
The thermocycler conditions used for amplification were initial denaturation at 95°C for 5 minutes, followed by 35 cycles at 95°C for 30 seconds, 55°C for 30 seconds, 72°C for 1 minute, and finally 72°C for 7 minutes. The amplified PCR products were resolved in 1% agarose gel, stained with ethidium bromide, visualized and photographed with Chemi Doc XRS (Bio-Rad Laboratories, Hercules, CA, http://www.bio-rad.com).
Integration of Bioarray Data with SAGE Data
The SAGE tag frequency table was downloaded from the Web-site (http://bit.fmrp.usp.br/msc_tags) provided by M.A. Zagos laboratory [16], inspected, and integrated with our microarray data. To facilitate an integration process between the two data-sets, microarray gene probes and SAGE tags were mapped to the NCBI Human Unigene clusters. If more than one tag or probe corresponded to the same Unigene cluster, then the higher ranked one was selected. This integration process led to selection of a total of 7,898 uniquely represented Unigene clusters. The clusters were then sorted by the average normalized intensity in MSCs. A full list of the integrated data can be found at http://callisto.snu.ac.kr/hoeonkim/microarray.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
1, I
2, III
1, IV
1, IV
2, VI
2, and VI
3), CTHRC1, CRTL1, lumican, and fibulin 4, as well as ECM biogenesis factors, including two types of serpins (SERPINE1 and SERPINH2) and three lysyl oxidases (LOXL1, LOX and its variant), whereas the latter consists of three different types of IGFBPs (IGFBP-6, IGFBP-7, and CTGF), PRSS11, OSF-2, CYR61, Wnt5B, PLAB, FAP, GAS6, and follistatin. It also contains four genes encoding membrane proteinsThy-1, KDELR3, NDUFA4, and BNIP3in which the former two proteins are destined to plasma and endoplasmic reticulum, respectively, whereas the latter two belong to mitochondrial membrane proteins. There are also five genes for cytoskeleton-associated proteins: transgelin,
-B-crystallin, tropomyosin I, EPLIN, and RIL; four genes for soluble proteins: NNMT, PKCBP, P311, and PYCR1; and finally two nuclear genes: necdin and a LIM domain protein FHL2. It is noteworthy that among those genes, collagens types I, III, IV, and VI [27], transgelin [28], Thy-1 [29], and FAP [30] were known as characteristics of MSCs and previously identified in human MSCs. Of the five novel genes, MGC3047 and MGC17528 were recently predicted to encode an immunoglobulin superfamily protein limitrin and an S100 calcium-binding protein, A16, respectively. And our blast analysis of MGC3278 and FLJ12442, encoding hypothetical proteins containing 563 and 520 amino acids, respectively, showed that the former contains a DUF719 domain of unknown function but conserved in several eukaryotic proteins while the latter belong to a 5' nucleotidase protein family. However, the identity of AGENCOURT_6683145 could not be resolved.
Confirmation of Gene Expression by RT-PCR
To verify the gene expression profile determined by our microarray analysis, the expression levels of the top 10 genes in the subset were analyzed by RT-PCR, using total RNAs obtained from the four cell samples. The result showed that all tested genes were expressed highly in MSCs, but either weakly or negligibly in MNCs (Fig. 4A
). This differential expression pattern was in a good agreement with that from the microarray analysis, confirming the high fidelity in microarray data and analytical methods. Moreover, when an additional pair of fresh MNC and MSC samples was analyzed by RT-PCR, they exhibited a differential expression pattern (Fig. 4B
) that was consistent with those of the former samples (Fig. 4A
). This finding implies that the gene expression profile determined in this study can be extrapolated to most, if not all, UCB from healthy individuals.
|
-actin; two heat shock proteins: HSPA8 and HSPB1; and other cytosolic or extracellular proteins. Among known products, CTGF, TAGLN, COL1A1, and IGFBP7 belong to MSC-specific molecules, as mentioned earlier. Most of the rest are housekeeping genes whose expression patterns are more or less constant in all proliferating cells. A most abundant gene is EEF1A1, which is responsible for the enzymatic delivery of aminoacyl tRNAs to the ribosome.
|
|
| DISCUSSION |
|---|
|
|
|---|
For this study, we used UCB-derived cells because UCB, when compared with other sources of the stem cells, has profound benefits in cell procurement and storage. Furthermore, cells in the neonatal blood are less mature than adult cells so that they do not trigger an immense immune reaction in unrelated donor transplantation. Accordingly, they are more readily applicable to stem cellbased therapy and transplantation than are cells of other origins.
Before undertaking an analysis of the differential expression profile between the cells, we carefully examined the reproducibility in microarray experiments, as well as interdonor variation in the transcriptional profile. As it turned out, both microarray systems and all experimental steps from target preparation to data analysis were highly reproducible. Moreover, no significant donor-to-donor variation was found in either MNC or MSC transcriptome. Slightly greater variability in the MSC transcriptome could reflect divergent culture conditions that the cells might experience during cell isolation and expansion. Taken together, these preparatory results provided a basis for recognition and interpretation of differentially expressed genes between UCB-derived MSCs and MNCs, which was the specific goal of this study.
To identify genes preferentially expressed by MSCs versus MNCs, we ranked genes by a ratio of the intensity score of the two cell populations. The RT-PCR analysis confirmed the fidelity of these data by showing that the top 10 genes were all expressed in accordance with their differential expression patterns in microarray data. When the genes were chosen with a ratio above 50, a total of 47 different genes were pooled out. A majority of the genes were found to belong to ECM components or cytokines, indicating that both unique ECM environment and specific cytokine signaling are important determinants of the functional states of the cells. Apparently, seven different collagens mainly constitute the ECM of MSCs, while a number of the serpin and lysyl oxidase family proteins are actively involved in its biogenesis. Connective tissue growth factor, also called IGFBP8, is a most preferentially expressed gene by MSCs. It belongs to an IGFBP family and is known to bind insulin-like growth factors with relatively low affinity. Its high expression in MSCs was previously shown in SAGE experiments [1517].
Relatively few in number are genes encoding membrane or nuclear proteins. Among membrane proteins, Thy-1, previously known as a T cell and an MSC-related cell surface marker, is differentially expressed by MSCs, suggesting that this antigen can be used in efficient immunoselection of the cells from the neonatal blood, using magnetic bead technology or fluorescence-activated cell sorting. As for nuclear genes, necdin and an LIM protein FHL2 are uniquely identified. The former was previously known as a neuronal growth suppressor [32] and its disruption confers lethality to mouse [33], whereas the latter was recently related to epithelial ovarian cancer [34]. Neither of them has been previously described in MSCs, and therefore more studies are needed to determine their functional significance in stem cell properties such as self-renewal and differentiation.
Finally, we found that DNA microarray- and SAGE-generated transcriptomes were weakly correlated. When we examined the 50 most highly expressed genes in terms of microarray intensity, we found that all of them were high in frequency in the SAGE experiment. Since the two datasets were obtained not only by different analytical methods but also from cells of different sources (BM versus UCB), it might indicate that transcriptomes of BM- and UCB-derived MSCs are somewhat similar to each other, which is in good agreement with growing evidence for their equivalent cell characteristics and multilineage potential. This finding also indicates that the two high-throughput techniques are both relevant in absolute quantification of gene expression, especially for highly expressed genes, and they can be used in parallel in accurate assessment of a given transcriptome.
In conclusion, a genome-wide differential expression analysis of human MSCs was performed using UCB-derived cells as a model, and their specific gene expression profile was elucidated, for the first time to our knowledge, by this study. The data will provide the basis for further studies on the molecular mechanisms controlling various core stem cell properties of human MSCs.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. J. Roberts, D. Howard, L. D. Buttery, and K. M. Shakesheff Clinical applications of musculoskeletal tissue engineering Br. Med. Bull., June 1, 2008; 86(1): 7 - 22. [Abstract] [Full Text] [PDF] |
||||
![]() |
M.-S. Tsai, S.-M. Hwang, K.-D. Chen, Y.-S. Lee, L.-W. Hsu, Y.-J. Chang, C.-N. Wang, H.-H. Peng, Y.-L. Chang, A.-S. Chao, et al. Functional Network Analysis of the Transcriptomes of Mesenchymal Stem Cells Derived from Amniotic Fluid, Amniotic Membrane, Cord Blood, and Bone Marrow Stem Cells, October 1, 2007; 25(10): 2511 - 2523. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Baharvand, A. Fathi, D. van Hoof, and G. H. Salekdeh Concise Review: Trends in Stem Cell Proteomics Stem Cells, August 1, 2007; 25(8): 1888 - 1903. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ohnishi, T. Yasuda, S. Kitamura, and N. Nagaya Effect of Hypoxia on Gene Expression of Bone Marrow-Derived Mesenchymal Stem Cells and Mononuclear Cells Stem Cells, May 1, 2007; 25(5): 1166 - 1177. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Jeong, K.-M. Ko, S. Bae, C.-J. Jeon, G. Young Koh, and H. Kim Genome-Wide Differential Gene Expression Profiling of Human Bone Marrow Stromal Cells Stem Cells, April 1, 2007; 25(4): 994 - 1002. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Song, N. E. Webb, Y. Song, and R. S. Tuan Identification and Functional Analysis of Candidate Genes Regulating Mesenchymal Stem Cell Self-Renewal and Multipotency Stem Cells, July 1, 2006; 24(7): 1707 - 1718. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Aslan, Y. Zilberman, L. Kandel, M. Liebergall, R. J. Oskouian, D. Gazit, and Z. Gazit Osteogenic Differentiation of Noncultured Immunoisolated Bone Marrow-Derived CD105+ Cells Stem Cells, July 1, 2006; 24(7): 1728 - 1737. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Yu, J. Fuhr, E. Boye, S. Gyorffy, S. Soker, A. Atala, J. B. Mulliken, and J. Bischoff Mesenchymal Stem Cells and Adipogenesis in Hemangioma Involution Stem Cells, June 1, 2006; 24(6): 1605 - 1612. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| STEM CELLS | THE ONCOLOGIST | CME | ALPHAMED PRESS JOURNALS |
