Small nucleolar RNA signatures of lung tumor-initiating cells
© Mannoor et al.; licensee BioMed Central Ltd. 2014
Received: 25 February 2014
Accepted: 17 April 2014
Published: 6 May 2014
Skip to main content
© Mannoor et al.; licensee BioMed Central Ltd. 2014
Received: 25 February 2014
Accepted: 17 April 2014
Published: 6 May 2014
Non-small cell lung cancer (NSCLC) is the number one cancer killer. Tumor-initiating cells (TICs) are responsible for tumor progression and recurrence. Emerging evidences suggest that small nucleolar RNAs (snoRNAs) play malfunctioning roles in lung tumorigenesis. This study aims to determine if snoRNAs have important function in lung TICs by: 1) profiling and comparing snoRNA expression patterns in lung ALDH1+/- cells of 28 primary NSCLC tissues to identify new signatures of TICs; 2) determining prognostic significance of the snoRNA signatures by analyzing the expression in 82 NSCLC tissues with different stages and histological types using quantitative PCR; 3) functionally investigating if the snoRNAs contribute to stemness of lung TICs using in vitro and in vivo assays.
Twenty-two snoRNAs were identified whose changes were specific to the TICs. The expression of two snoRNAs (snoRA3 and snoRA42) was inversely associated with survival of NSCLC patients (P = 0.002, p = 0.001, respectively). Functional analysis indicated that snoRA42 was upregulated in CD133+ cells isolated from NSCLC cell lines compared with the CD133- counterparts. snoRA42 knockdown reduced the proliferation and self-renewal of TICs in vitro. However, ectopic expression of snoRA42 in non-TICs enhanced the potentials of cell proliferation and self-renewal. snoRA42 expression was associated with expression of stem cell-core transcription factors in lung TICs. Blocking snoRA42 expression in TIC xenografts decreased tumorigenesis in mice.
The snoRNA signatures of lung TICs provide potential biomarkers for predicting outcome of NSCLC. snoRA42 is one of the important snoRNAs in regulating features of lung TICs, and thus contributes to lung tumorigenesis.
Non-small cell lung cancer (NSCLC) is the leading cause of cancer death for men and women worldwide. With availability of more sensitive radiological imaging studies, more NSCLC patients will be diagnosed while the disease is still at early stage [1, 2]. The standard of care for NSCLC is surgery, often followed by chemotherapy . However, approximately 84% of those diagnosed with lung cancer will die within five years . Furthermore, the current chemotherapies often have toxicity in normal host tissues, whereas tumor cells rapidly develop resistance to anticancer agents. The developments of biomarkers for identifying NSCLC patients at high risk of recurrence after surgery, and therapeutic targets for safe and effective treatment of lung cancer are clinically important.
The existence of tumor-initiating cells (TICs), also known as cancer stem cells, could explain why the current chemotherapies cannot consistently eradicate tumors, because the therapies only target the bulk of tumor cells and are unable to eliminate TICs . Furthermore, residual lung TICs may regenerate a cancer cell population, leading to tumor relapse after therapy. TICs have been identified in lung cancer by using several approaches such as CD133, a cell surface marker . We have recently characterized ALDH1+ cancer cells are TICS, as the ALDH1+ cancer cells have extensive self-renewal, proliferative, and in vivo tumorigenic potentials [5–7]. The analysis of molecular aberrations that characterize TICs would deep our understanding of lung tumor biology. Furthermore, the molecular changes could be developed as a new diagnostic system for monitoring outcome of NSCLC. In addition, the TIC-related molecular changes may enable the development of specific agents for eradicating the tumor-maintaining cells, and thus provide efficient therapeutic approaches for lung cancer.
Non-coding RNAs (ncRNAs) are functional transcripts that do not code for proteins, however, play a crucial role in regulating gene expression . ncRNAs include small nucleolar RNAs (snoRNAs), microRNAs (miRNAs), short interfering RNAs (siRNAs), piwi-associated RNAs, small Cajal body-specific RNAs (scaRNAs), snRNAs (small nuclear RNAs), and long ncRNAs . Of the small ncRNAs, miRNAs have extensively been studied in carcinogenesis [9, 10]. Dysregulation of some miRNAs is associated with features of TICs. For example, elevated miR-181 clusters were identified as vital regulators in hepatic TICs . Furthermore, up-regulation of miR-199b-5p impaired the development of TICs of medulloblastoma though repression of HES1 . In addition, the restoration of miR-34 expression suppressed the self-renewal of pancreatic TICs .
Recently, new and unexpected functions of other types of small ncRNAs have been discovered, revealing that the molecules have highly diverse roles and are actively involved in the processes of carcinogenesis than previously thought . In particular, several studies including our own data suggest that snoRNAs exhibit differential expression patterns in lung tumor and have capability to affect cell transformation, tumorigenesis, and metastasis of NSCLC [10, 14, 15]. snoRNAs range between 60–150 nucleotides in length . There are two classes of snoRNAs: box C/D snoRNAs (snoRDs) and box H/ACA snoRNAs (snoRAs) . snoRDs serve as guides for the 2′-O-ribose methylation of rRNAs or snRNAs, whereas snoRAs are guides for the isomerization of uridine residues into pseudouridine . Accumulated evidence suggests that snoRNAs can target other RNAs including snRNAs and messenger RNAs . For instance, HBII-52 regulates alternative splicing of 5-HT2CR by binding to a silencing element in exon Vb [18, 19]. Therefore, profiling the expression patterns of snoRNAs in lung TICs may help comprehend the functional repercussions of the ncRNAs in the development and progression of NSCLC, and hence provide diagnostic biomarkers and therapeutic targets for the disease.
We propose to determine if aberrant snoRNAs have important function in lung TICs by: 1) profiling and comparing snoRNA expression patterns in lung ALDH1+/- cells of primary NSCLC tissues to identify signatures of TICs; 2) determining prognostic significance of the signatures by analyzing the expression in NSCLC tissues; 3) functionally investigating if the snoRNA signatures may contribute to stemness of lung TICs.
Up-regulated and down-regulated snoRNAs in TICs vs. non-TICs
Name of snoRNA
We explored whether snoRA42 downregulation could inhibit cell growth and proliferation of lung TICs by transfecting CD133+ cells isolated from NSCLC cell lines with snoRA42-siRNA. snoRA42-siRNA produced more than 75% reduction of snoRA42 expression in the CD133+ cells compared with mock-transfected CD133+ cells or CD133+ cells transfected with the scrambled siRNA (Figure 3B). The results suggested that the snoRA42-siRNA could block snoRA42 expression in TICs. Furthermore, there was a significant decrease of numbers of CD133+ cancer cells transfected with snoRA42-siRNA compared with those with scrambled siRNA or mock transfection (Figure 3C). The inhibition of snoRA42 silencing on cell growth of TICs was confirmed by methylthiazol tetrazolium (MTT) assay. The findings suggest that snoRA42 knockdown could have a significant inhibition on growth and proliferation of lung TICs. In addition, CD133+ cancer cells transfected with snoRA42-siRNA gave rise to significantly fewer mammospheres (Figure 3D) compared with CD133 + cancer cells treated with controls, implying the downregulation of snoRA42 decreased the capacity of mammosphere formation of TICs. Therefore, reduced snoRA42 expression in TICs could inhibit the self-renewing capacity of lung TICs.
To investigate possible mechanism underlying the inhibition of in vitro tumorigenicity of TICs, we first evaluated whether there was apoptosis in CD133+ cancer cells treated with snoRA42-siRNA. As shown in Figure 4D, percentages of apoptotic cells were significantly increased upon snoRA42 knockdown compared with scrambled siRNA- and mock-treated CD133+ cancer cells. To confirm the finding, we applied enzyme-linked immunosorbent assay (ELISA) to analyze cell lysate for determining change of caspase-3. As shown in Figure 4E, cleaved caspase-3 was activated prominently in CD133+ cells treated with snoRA42-siRNA. Therefore, snoRA42 knockdown could trigger apoptosis in CD133+ cells isolated from NSCLC cell lines. The inhibition of cell proliferation and self-renewal of lung TICs by snoRA42 suppression may be achieved at least partially through inducing apoptosis.
To determine if snoRA42 suppression could affect expression of core stem cell transcription factors in CD133+ cells isolated from a panel of NSCLC cell lines, we quantified expression levels of stem cell-associated genes by qRT-PCR. siRNA-mediated knockdown of snoRA42 in CD133+ cells led to a significant decrease in transcript levels of the genes, including Oct4, Nanog, Sox2, Notch1, Smo, and ABCG2 compared to scrambled siRNA-treated CD133+ cancer cells (Figure 6). Furthermore, to examine if enforced expression of snoRA42 could have an opposite effects to that of snoRA42 knockdown, we transduced CD133- cancer cells with pCMV-snoRA42. Although the expression levels of different stem cell-associated genes in snoRA42-transduced CD133- cancer cells could not reach to that in CD133+ cancer cells, transcriptional level of the stem cell-associated genes in CD133- cancer cells with pCMV-snoRA42 significantly increased compared to CD133- cancer cells with empty vector. Altogether, snoRA42 is associated with expression of the stem cell-core transcription factors in lung TICs.
Cancer arises from a tumorigenic subpopulation of TICs. The identification of molecular signatures of lung TICs provide a key standpoint for better understanding tumorigenesis and developing prognostic biomarkers and targeted therapy . For instance, methylation of TIC-associated Wnt target genes were identified in colorectal tumors and showed the potential for predicting a poor prognosis of colon cancer patients . Yet few molecular signatures of TICs in lung cancer have been identified. By conducting transcriptomic profiling ALDH1+/-NSCLC cells of primary tumor tissues, here we successfully identify 22 novel snoRNA signatures of lung TICs.
Although only two (snoRA3 and snoRNA42) of 22 snoRNAs are evaluated in the present study, the results indicate that the aberrant expression of the TIC-snoRNAs is associated with aggressive biological behavior of NSCLC. Importantly, expression of the snoRNAs was inversely associated with overall survival of the patients. Therefore, the assessment of the TIC-snoRNA signatures might provide a potential approach for predicting outcome of NSCLC. However, as lung tumor is a heterogeneous disease, the analysis of the two TIC-snoRNA signatures cannot achieve the accuracy required to move forward for a clinical trial. Multiple biomarkers may provide more accuracy in prediction of outcome of lung cancer. From the identified 22 snoRNA signatures of lung TICs, we are optimizing a panel of biomarkers that can more precisely predict recurrence of NSCLC in our ongoing study.
snoRNAs have long been believed to modify and stabilize rRNAs for ribosome production. Several recent studies including our own observations demonstrate that dysregulation of snoRNAs is a common molecular event in carcinogenesis [14, 15]. We previously found that snoRA42 activation promoted lung cancer development and progression . The critic role of snoRA42 in TICs is first supported by that expression of snoRA42 is significantly elevated in CD133+ lung cancer cells compared with CD133- lung cancer cells. Furthermore, CD133- lung cancer cells with forced snoRA42 expression have the capacity to enhance the in vitro tumorigenicity, whereas the CD133- cells without forced snoRA42 expression do not. In addition, siRNA-mediated depletion of snoRA42 causes reduction of expression of stem cell-related genes, suppression of xenograft tumor formation, and inhibition of cellular proliferation and self-renewal of TICs in vitro. The inhibition of cell proliferation and self-renewal of lung TICs by snoRA42 dysregulation may be attained through inducing caspase 3-dependent apoptosis. Moreover, forced expression of snoRA42 in CD133- lung cancer cells produces higher cell proliferation and expression levels of stem cell-associated genes compared with CD133- lung cancer cells transduced with empty vector. Therefore, the present study extends the previous discoveries by finding that snoRA42 plays an important role in regulating stemness of lung TICs, and hence contributes to invasiveness and metastasis of NSCLC. Targeting snoRA42 in future may serve as an effective and specific therapeutic approach for lung cancer treatment. Nevertheless, investigating how the stem cell-core factors are regulated by snoRA42, and if other snoRNA signatures are also responsible to stemness features of TICs is warranted. In fact, we are functionally characterizing snoRA3 by using in vitro and in vivo approaches to determine if snoRA3 has a significant role in regulating stemness of lung TICs.
We identified novel snoRNA signatures of lung TICs. Two of the signatures were evaluated in a cohort of NSCLC specimens. The results showed that the overexpression of the snoRNAs was inversely associated with survival of NSCLC patients. In vitro and in vivo functional analyses highlighted the importance of snoRA42 for the features of lung TICs. The snoRNA signatures of lung TICs may provide potential biomarkers for predicting outcome of NSCLC and novel therapeutic targets.
Demographic and clinical characteristics of 82 NSCLC patients
No. of patients (%)
Age at diagnosis
65.3 ± 9.5
Squamous cell carcinoma
The surgical tumor specimens were processed for preparing single dissociated cells as previously described . Briefly, the surgical specimens were washed three times and left 12 hours in DMEM-F12 medium supplemented with penicillin/streptomycin and amphotericin B. Tissue dissociation was performed by adding enzymatic digestion with collagenase II (Gibco-Invitrogen, Carlsbad, CA) for two hours at 37°C. The dissociated cells were processed for flow cytometric isolation of TICs based on ALDH enzymatic activity as described below.
All NSCLC cell lines including A549, H299, H1792, H226, H522, H-1944, Calu-1, and H460 were purchased from the ATCC (Manassas, VA, USA). The cell lines were maintained as monolayer cultures in high glucose Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 50 μg/mL normocin (Invivogen, San Diego, CA) at 37°C in a humidified incubator with 5%.
We isolated ALDH1 + cells from single dissociated cells of lung tumors using an Aldefluor assay kit (Aldagen Inc., Durham, NC) as described in our previous reports [5–7]. 1 × 106 cells per well were suspended in Aldefluor assay buffer containing ALDH substrate, BAAA (1 μmol/L). Cells stained with 50 mmol/L of the specific ALDH inhibitor diethylaminobenzaldehyde (DEAB) were used as a negative control. Flow cytometric sorting was carried out using a FACSAria (Becton Dickinson, Mountain View, CA). Aldefluor fluorescence was excited at 488 nm, and fluorescence emission was detected using a fluorescein isothiocyanate (FITC) 530/30-nm band-pass filter by a FACS calibre machine (BD Biosciences, San Jose, CA).
CD133+ and CD133- cells were isolated from NSCLC cell lines as previously described . Briefly, cells were derived after trypsinization of adherent NSCLC cells. Cells were washed twice with magnetic automated cell sorting (MACS) buffer (Miltenyi Biotec Inc., San Diego, CA). Cell pellets were then resuspended in 300 μl of MACS buffer followed by addition of 100 μl of FcR blocking reagent (Miltenyi Biotec Inc.), 100 μl CD133 of microbeads (Miltenyi Biotec Inc.), and 50 μl of CD133/2-PE antibody (Miltenyi Biotec Inc.). The cells were washed with MACS buffer and MACS-assisted purification was performed according to manufacturer’s instructions. Purity of CD133+ cells and CD133- cells were checked by flowcytometry.
For sphere-forming cultures, cells were plated at a density of 2 × 104 cells per well in six-well ultra-low attachment plates (Corning Inc., Corning, NY) in serum free Bronchial Epithelial Cell Basal Medium (BEBM), supplemented with the following growth supplements: BPE, Hydrocortisone, hEGF, Epinephrine, Transferrin, Insulin, Retinoic Acid, Triiodothyronine, and GA-1000, provided as BEGM prepackaged SingleQuots (Lonza, Walkersville, MD) plus rhEGF (Sigma, St Louis, MO), bFGF (Sigma). Fresh aliquots of EGF and bFGF were added every three days. Floating spheres were subjected to mechanical dissociation, followed by replating of single cells for propagation into a second, and subsequently into a third batch of suspension culture as described in our previous reports [5–7].
snoRNA profiling was performed by using GeneChipR Array (Affymetrix, Inc, Santa Clara, CA) as described in our previous report . Microarray experiments were done with matched ALDH1+ and ALDH1- cells. One μg total RNA was labeled with Biotin FlashTag Biotin Labeling Kit (Affymetrix, Inc). The labeling reaction was hybridized on the arrays in Affymetrix Hybridization Oven 640 (Affymetrix, Inc) for 16 hours. The arrays were stained with Fluidics Station 450 using fluidics script FS450_0003 (Affymetrix, Inc), and then scanned on a microarray scanner (Axon Instruments Inc, Foster City, CA). snoRNA probe outliers were defined as per the manufacturer’s instructions (Affymetrix, Inc) and further analyzed for data summarization, normalization, and quality control by using the web-based QC Tool software (http://www.affymetrix.com). The normalized microarray data underwent further analysis as described in our previous study . We performed tree visualization by using Java Treeview 1.0 (Stanford University School of Medicine, Stanford, CA).
Expressions of snoRNAs were determined by using real-time SYBR green RT-qPCR assay as described in our previous studies [14, 15]. All PCR reactions were run on a CFX96 Real-Time qPCR detection system (Bio-Rad, Hercules, CA) using Brilliant SYBR green qPCR master mix (Agilent Technologies, Chelmsford, MA). U6 was used as an internal control gene. ΔCt was calculated by subtracting the Ct values of U6 from those of the snoRNA tested, and fold-change of each snoRNA was determined by the equation 2–ΔΔCt. qRT-PCR was also performed to decide expression levels of stem cell-associated genes, which included ACTBF, ABCB1, ABCG2, CD133, Nanog, Nestin, Notch, Oct4, Smo, and Sox2. The primer sequences are included in Additional file 1: Table S3.
RNA interference assays were performed as described in our previous study . Briefly, specific target site of designed siRNA in RNA sequence of snoRA42 (NCBI accession number NR_002974) was a 19 nucleotide stretch as follows: 5′-GTACCCATGCCATAGCAAA-3′. Corresponding non-targeting scrambled sequence was also designed. Transfection was performed by using Opti-MEM and Lipofectamine™ RNAiMAX Reagent (Invitrogen, CA, USA) as described elsewhere .
pCMV (pSilencer CMV 4.1; Ambion, Carlsbad, CA)-snoRA42 vector was constructed for ectopic snoRA42 expression as described previously . A vector comprising sequence that showed no significant homology to the human genome databases was used as a control. Transduction of cells with the vectors was performed as described in our previous study .
MTT and BrdU incorporation assays were performed as previously described [24, 14]. Cells were seeded in 96-well plates at 1 × 104 cells per well. Colony formation was determined by using 3D soft gar matrix with a 0.8% of soft agar (noble agar) base layer followed by 0.4% top agar layer containing cells at a density of 1 × 104 with different treatments in 24-well culture plates. Colonies were counted using inverted light microscopy. For transmigration assay, the ability of migration and invasion of cells was examined using BD BioCoat Matrigel Invasion Chambers (BD Biosciences) according to manufacturer’s instruction. After 24 h, the migrated cells were pelleted and counted using a hematocytometer. To further confirm cell migration, a wound healing assay was performed in 12 well plates (1 × 105 per well). When the cells were grown to 90 to 95% confluences, tranfection of CD133+ cells with siRNA-snoRA42 and scrambled siRNA, as well as mock transfection were performed. Wound lines were created manually by scratching the monolayer with a sterile 200 ml pipette tip and migration of the cells was assessed after 24 hours. Pictures were taken using Nikon inverted phase-contrast microscope (Nikon, Melville, NY). The distance between the parallel lines was measured using ImageJ software. All experiments were carried out at least three times.
Cell apoptosis was determined by two-color immunofluorescence staining followed by flowcytometric analysis. Cells were trypsinized without EDTA (Sigma) followed by two washing steps with phosphate buffered saline. 2 × 105 cells per well were then incubated with 5 μl of FITC Annexin V (BD Pharmingen) and 5 μl of propidium iodide (BD Pharmingen) in 400 μl of 1 × binding buffer for 15 min at room temperature. Data acquisition and analysis were performed on a FACScan flowcytometer (Becton-Dickinson) using BD analysis software as described previously .
Cleaved caspase 3 ELISA was performed by using PathScan Sandwich ELISA kit (Cell signaling Technology, Inc., MA) according to manufacturer’s instructions. Cytochrome c-treated jurkat cell lysate was used as cleave caspase 3 positive control, whereas untreated jurkat cell lysate was used as cleaved caspase 3 negative control (Cell signaling Technology, Inc.). Cleaved caspase index was calculated as follows: (sample O.D. - negative control O.D./positive control O.D.-negative control O.D.)/100.
The animal study protocol was performed in accordance with the guidelines of University of Maryland School of Medicine Institutional Animal Care and Use Committee. Six 6-week old female SCID/Beige mice per group were subcutaneously inoculated with 5 × 105 Calu-1-derived CD133+ cells with different treatments. The mice were observed for six weeks and then euthanized with CO2. Tumor volume was calculated by using the formula: width × width × length × 0.52 [5–7, 14].
Both univariate and multivariate Cox proportional hazard models were applied to assess the effect of various clinically and histopathologically interesting variants, and snoRNA expression on overall survival data of patients. Association of expression of the snoRNAs with overall survival rate was also analyzed using the Kaplan-Meier method. All P-values were determined by two-sided tests.
This work was supported in part by NCI R01CA161837, VA merit Award I01 CX000512, LUNGevity/Upstage Foundation Early Detection Award, and University of Maryland Cancer Epidemiology Alliance Seed Grant (F.J.).Multidisciplinary Research Career Development Program, and an exploratory research grant from Maryland Stem Cell Fund (F. J.).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.