Molecular Cancer Gene Expression Profiling of Cancer Stem Cell in Human Lung Adenocarcinoma A549 Cells

Background: The 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.


Background
Cancer 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][4][5][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].
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].

The distinct gene regulations in SP cells
We 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.

Validation of gene regulations
To 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).

Discussion
Based 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]. Figure 1 Sorting of SP and non-SP cells by FACSVantage SE. 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].

Sorting of SP and non-SP cells by FACSVantage SE
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 upregulated NR0B1, though up-regulated NR0B1 gene was required for the transformed phenotype of certain sarcoma [45].

Conclusion
It 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.

Fluorescence activated cell sorting (FACS) for SP and non-SP cells
The

RNA isolation and oligonucleotide microarray analysis
We prepared total RNA from approximately 2 × 10 6 SP or non-SP cells with Qiagen RNeasy Mini kit (Qiagen, CA, Gene clustering of up-regulated genes in SP and non-SP cells Figure 2 Gene clustering of up-regulated genes in SP and non-SP cells. After normalizing each chip to the 50 th percentile of the measurements taken that chip, gene-probes scored less than 0.1 either in SP or non-SP were excluded from data analysis. Only matched up-regulated genes in SP compared to non-SP cells are selected in each step of chip data analysis. The 12 genes and 46 genes were considered as up-regulated in SP and non-SP cells, respectively. 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; Fprimer 5'GTGGAAGCTGACCAGGTTGT3' and R-primer 5'AAGCCGTGTTCTTTCTGCTG3'.

Microarray data analysis
We 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 50 th 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 upregulated 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 interests
The author(s) declare that they have no competing interests.

Authors' contributions
DS 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 prep- The up-regulation of AKR1C1 and TM4SF1 in smokers (derived from GEO public database) Figure 4 The up-regulation of AKR1C1 and TM4SF1 in smokers (derived from GEO public database). AKR1C1 (A) and TM4SF1 (B) were increased in current smoker compared to former or never smokers. TM4SF1 was increased in metastatic form compared to primary form of colon cancer (C).
Relative fold changes of AKR1C1/C2 gene in SP and non-SP cells Figure 3 Relative fold changes of AKR1C1/C2 gene in SP and non-SP cells. The fold changes of AKR1C1/C2 between SP and non-SP cells were compared using GeneChip data and quantitative real time-reverse transcriptase PCR. The data (n = 3) were presented as mean ± SD. aration of the manuscript. All authors read and approved the final manuscript.