Gene expression profiling of cancer stem cell in human lung adenocarcinoma A549 cells1 Department of Veterinary Pharmacology and Toxicology College of Veterinary Medicine, Konkuk University, 1 Hwayang-dong, Kwangjin-gu, Seoul 143-701, Republic of Korea 2 Department of Public Health, College of Veterinary Medicine, Konkuk University, 1 Hwayang-dong, Kwangjin-gu, Seoul 143-701, Republic of Korea 3 Department of Infectious Diseases, College of Veterinary Medicine, Konkuk University, 1 Hwayang-dong, Kwangjin-gu, Seoul 143-701, Republic of Korea 4 Cell and Gene Therapy Research Institute, Pochon CHA University, CHA General Hospital, Seoul 135-081, Republic of Korea
Molecular Cancer 2007, 6:75doi:10.1186/1476-4598-6-75 The electronic version of this article is the complete one and can be found online at: http://www.molecular-cancer.com/content/6/1/75
©
2007 Seo et al; licensee BioMed Central Ltd. AbstractBackgroundThe studies on cancer-stem-cells (CSCs) have attracted so much attention in recent years as possible therapeutic implications. This study was carried out to investigate the gene expression profile of CSCs in human lung adenocarcinoma A549 cells. ResultsWe isolated CSCs from A549 cell line of which side population (SP) phenotype revealed several stem cell properties. After staining the cell line with Hoechst 33342 dye, the SP and non-side population (non-SP) cells were sorted using flow cytometric analysis. The mRNA expression profiles were measured using an Affymetrix GeneChip® oligonucleotide array. Among the sixty one differentially expressed genes, the twelve genes inclusive three poor prognostic genes; Aldo-keto reductase family 1, member C1/C2 (AKR1C1/C2), Transmembrane 4 L six family member 1 nuclear receptor (TM4SF1), and Nuclear receptor subfamily 0, group B, member 1 (NR0B1) were significantly up-regulated in SP compared to non-SP cells. ConclusionThis is the first report indicating the differences of gene expression pattern between SP and non-SP cells in A549 cells. We suggest that the up-regulations of the genes AKR1C1/C2, TM4SF1 and NR0B1 in SP of human adenocarcinoma A549 cells could be a target of poor prognosis in anti-cancer therapy. BackgroundCancer stem cell hypothesis is the tumoral cells which have stem cell features such as self-renewal, high migration capacity, drug resistance, and aberrant differentiation which constitute the heterogeneous population of tumor [1,2]. Tissue-specific stem cells are defined by their ability to self-renew and to produce the well differentiated and functional cells within an organ. Differentiated cells are generally short-lived; in skin and blood for example, they are produced from a small pool of long-lived stem cells that last throughout the life [3-6]. Therefore, stem cells are necessary for tissue development, replacement, and repair [7]. On the other hands, the longevity of stem cells make them susceptible to accumulating genetic damage and thereby representing the growth root for cancer recurrence following treatment [8]. It was reported that some of the tumor stem cells can survive chemotherapy and support re-growth of the tumor mass [9]. Cancer stem cells (CSCs) were first identified in 1990s in hematological malignancies, mainly acute myelogenous leukemia (AML) and also in other subtypes like AML M0, M1, M2, M4, and M5, chronic myeloid leukemia (CML), acute lymphoblastic leukemia (ALL), and multiple myeloma [10,11]. CSCs are also known in solid tumors like breast, brain, lung, prostrate, testis, ovary, stomach, colon, skin, liver, and pancreas [12-17]. A character of stem cells, termed "side population (SP)", has been identified using Hoechst 33342 dye. The flow cytometric analysis makes sorting possible either to SP or non-SP cells. The SP cells have been isolated from various types of adult tissue where they demonstrate stem cell activity [18-23]. The findings of these previous studies suggest that the SP phenotype represents a common feature of stem cells. We performed our work on human lung adenocarcinoma A549 cells (of which SP phenotype revealed several stem cell properties [24]) to identify the genes, which make the CSCs of poor prognostic phenotype and evaluate the gene expression intensities of SP and non-SP cells using oligonucleotide micro-array. The reasons why the A549 cell line was selected, because it has a relatively high proportion of SP cells compared to other cell lines [25] and is more chemo-resistant particularly to platinum drugs [26]. ResultsThe distinct gene regulations in SP cellsWe sorted A549 cell line to SP and non-SP cells (Fig. 1) and compared the gene expression intensities of both cells. Official symbols and gene names were used in accordance with the symbol and name lists approved by HUGO (Human genome organization) Gene Nomenclature Committee (Table 1) [27]. Following data analysis, 12 genes were considered as up-regulated in SP cells (TM4SF1 has 2 probe ID) (fold changes are shown in Table 2), whereas, 49 genes were down-regulated (Fig. 2). Since we focused on distinct gene regulations, the student's t'-test was not employed to prevent loss of up-regulated genes in all of three chip data, though it had large chip variations.
Table 1. The approved gene symbols and names in reference to HUGO Gene Nomenclature
Table 2. Gene list up-regulated in SP cells compared to non-SP cells Validation of gene regulationsTo confirm the fold changes of AKR1C1 in chip data, quantitative real time – reverse transcriptase PCR was employed. The relative fold changes in SP compared to non-SP cells were 3.11 ± 0.92 and 2.88 ± 0.17 in microarray and qrt- rtPCR, respectively (Fig. 3).
DiscussionBased on the cancer stem cell hypothesis, we assumed that the up-regulation of certain genes that are related to poor prognosis (high migration capacity or drug resistance) in SP of cancer cells could be a target for therapeutic index. In the present study, we found some genes that are related to drug resistance, such as AKR1C1/C2 and NR0B1, or cancer metastasis, such as TM4SF1, were up-regulated in SP cells of human lung adenocarcinoma A549 cell line. Furthermore, the up-regulated gene, ABCG2, has been noticed to be as an indicator for sorting SP cells by Hoechst 33342 staining [24]. It was reported that ABCG2 pumping out the drugs was associated with multi-drug resistance in many cancers [28,29] and/or effects higher levels of DNA repair and hence lowered the ability to apoptosis [30]. AKR1C belongs to a superfamily of monomeric, cytosolic NADP(H)-dependent oxidoreductases that catalyzes the metabolic reduction and either activate or inactivate several xenobiotics [31,32]. In humans, at least four isoforms of AKR1C (AKR1C1~4) have been identified [33]. AKR1C1/AKR1C3 was known to inactivate progesterone, which could alter the progesterone/estrogen ratio in certain cancers [34,35]. Additionally, AKR1C1/C2 inhibitors were reported as potential anti-neoplastic agents [36,37]. The high expression of AKR1C was considered as a poor prognostic factor in patients with non-small-cell lung cancer [38], and it was enriched in hepatocellular carcinoma than normal hepatic cells [39]. A previous study evaluated the relationship between the AKR1C and drug resistance and revealed that the over expression of AKR1C1/C2 led to drug resistance in non-small-cell lung cancer cells [26]. A public database (Gene Expression Omnibus (GEO), NCBI) [40] has shown significant up-regulation of AKR1C1 in smokers (Figure 4A). Similar trends were also reported in lung cancer patients [41,42]. This means that smoking can alter the regulation of certain genes related to poor prognosis such as AKR1C1. Moreover, AKR1C1 was suggested as a marker of stem-like cells in thyroid cancer cell lines, though it was not proved by the authors [41].
From the data of GEO, the TM4SF1 gene, which is believed to be involved in cancer invasion and metastasis [44] was also up-regulated in smokers (Figure 4B), and metastatic form of colon cancer patients compared to primary form of colon cancer patients (Figure 4C). TM4SF1 was also suggested as a possible marker of stem-like cells in thyroid cancer cells [41]. However, we could not find out the GEO data that match to tumor genesis of the up-regulated NR0B1, though up-regulated NR0B1 gene was required for the transformed phenotype of certain sarcoma [45]. ConclusionIt is still unclear how the cancer stem cells express poor prognostic phenotype in cancer, but we found for the first time that several genes related to chemo-resistance and metastasis are up-regulated in SP cells compared to non-SP cells. Therefore, enrichment of AKR1C1/C2, TM4SF1 and NR0B1 in SP of A549 cells might be a target of poor prognosis in cancer therapy. MethodsCell cultureA549 cells, a representative human lung adenocarcinoma cell line, were obtained from American Type Culture Collection (ATCC, VA, USA). The cells were cultured on F-12K nutrient mixture, Kaighn's modification (1×) liquid (Gibco, CA, USA) supplemented with 10% fetal bovine serum (FBS) (Genetron Life Technology, INC., CA, USA) and 1% of penicillin/streptomycin (P/S) (Gibco, CA, USA) on standard plastic tissue culture dishes (SPL Life Sciences, Seoul, Republic of Korea) and incubated in an atmosphere of 95% air/5% CO2 at 37°C. Fluorescence activated cell sorting (FACS) for SP and non-SP cellsThe cells were detached with trypsin (Gibco, CA, USA), washed with phosphate buffer solution (PBS)/2% FBS, and resuspended at 1 × 106 cells per ml in pre-warmed Hanks' balanced salt solution (HBSS; Gibco, CA, USA)/10% FBS. Hoechst 33342 dye (Sigma, St. Louis, MO, USA) was then added to this portion (final concentration: 5 μg/ml), and incubation continued for 90 min at 37°C. After washing with PBS/2% FBS, the cells re-suspended in ice-cold HBSS/10% FBS were labeled with 2 μg/ml propidium iodide (Sigma, St. Louis, MO, USA) to distinguish live from dead cells prior to analysis. SP analysis and sorting were done using a FACSVantage SE (BD Biosciences, CA, USA). The Hoechst 33342 dye excited at 350 nm using UV laser and the DM610 SP was used to distinguish the red from the blue fluorescence signals. The EF675-LP (Hoechst Red) and BP450/20 nm (Hoechst Blue) filters were then installed in front of the PMT detector. RNA isolation and oligonucleotide microarray analysisWe prepared total RNA from approximately 2 × 106 SP or non-SP cells with Qiagen RNeasy Mini kit (Qiagen, CA, USA) according to the instruction of the manufacturer. The RNAs were subjected to GeneChip® expression array in full commercial service with two-cycle target labeling (SeouLin Bioscience, Seoul, Republic of Korea). Briefly, cDNA were synthesized from total RNA using T7-Oligo (dT) primers. Using that cDNA, biotinylated cRNA was then synthesized. Fifteen μg of the labeled cRNA was hybridized to a Human Genome U133 Plus 2.0 Array (Affymetrix, CA, USA). Array image was scanned and analyzed using Genechip operating software (GCOS) (Affymetrix, CA, USA). Quantitative real-time rt-PCR (qrt-rtPCR)To validate the fold changes in the expression intensity of SP and non-SP cells in microarray data, we performed qrt-rtPCR using SYBR Premix Ex Taq Perfect Real Time kit (TaKaRa Bio, Ohtsu, Japan) in a SmartCycler (Cepheid, CA, USA) according to the manufacturer's instruction. The cycle threshold value, which was determined using second derivative, was used to calculate the normalized expression of the indicated genes using Smartcycler Software (Cepheid, CA, USA). The following primer pairs were used: GAPDH (as an internal control); F-primer 5'CGACCACTTTGTCAAGCTCA3' and R-primer 5'AGGGGAGATTCAGTGTGGTG3', AKR1C1;AKR1C2; F-primer 5'GTGGAAGCTGACCAGGTTGT3' and R-primer 5'AAGCCGTGTTCTTTCTGCTG3'. Microarray data analysisWe used GeneSpring GX 7.3.1 software (Agilent Technologies, CA, USA) to normalize and analyze the microarray data. Following normalization of each chip to the 50th percentile of the measurements taken that chip, probes (genes) scored less than 0.1 both in SP and non-SP cells were excluded from the data analysis. Only matched up-regulated genes in SP and non-SP cells were selected from each step of chip (gene expression microarray) data analysis. More than 2 fold changes in all chip data were considered as up-regulated. Standard curve method was used to analyze qrt-rtPCR for confirmation of fold changes in chip data. Competing interestsThe author(s) declare that they have no competing interests. Authors' contributionsDS and HS conceived the study. DS, JS, HC and HY carried out the sample preparation, gene expression analyses and quantitative real time RT-PCR. DK contributed to stem cell preparation. KS, IC, JK, AMAE and HS participated in the design, reviewed all data, and contributed in the preparation of the manuscript. All authors read and approved the final manuscript. AcknowledgementsThis work was supported by MOEHRD Basic Research Promotion Fund (KRF-2005-003-E00265) and the BK21 Program of Ministry of Education, Republic of Korea. References
Have something to say? Post a comment on this article! |





on Google Scholar







author email
corresponding author email
Figure 1.
Figure 2.
Figure 3.
Figure 4.