PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of oncotargetLink to Publisher's site
 
Oncotarget. 2017 February 7; 8(6): 9425–9441.
Published online 2016 December 23. doi:  10.18632/oncotarget.14116
PMCID: PMC5354742

Gene expression profiling of tumor-initiating stem cells from mouse Krebs-2 carcinoma using a novel marker of poorly differentiated cells

Abstract

Using the ability of poorly differentiated cells to natively internalize fragments of extracellular double-stranded DNA as a marker, we isolated a tumorigenic subpopulation present in Krebs-2 ascites that demonstrated the features of tumor-inducing cancer stem cells. Having combined TAMRA-labeled DNA probe and the power of RNA-seq technology, we identified a set of 168 genes specifically expressed in TAMRA-positive cells (tumor-initiating stem cells), these genes remaining silent in TAMRA-negative cancer cells. TAMRA+ cells displayed gene expression signatures characteristic of both stem cells and cancer cells. The observed expression differences between TAMRA+ and TAMRA− cells were validated by Real Time PCR. The results obtained corroborated the biological data that TAMRA+ murine Krebs-2 tumor cells are tumor-initiating stem cells. The approach developed can be applied to profile any poorly differentiated cell types that are capable of immanent internalization of double-stranded DNA.

Keywords: tumor-initiating stem cells, DNA internalization, RNAseq, Real Time PCR, TAMRA

INTRODUCTION

Tumor-initiating stem cells (TISCs) were discovered in the late 1990ies [1, 2]. They have a set of specific properties that are shared between stem and cancer cells. These cells are capable of self-renewal; they maintain the pluripotency state and determine the tumorigenic potential of a graft. One of the major features of TISCs is their ability to form metastases, accounted for by the change in the expression of the genetic network controlled by the epithelial-mesenchymal transition [3].

Understanding the molecular basis of TISC functioning and development of approaches to their eradication are related to two practical tasks: identification of these cells in a pool of cancer cells and establishing a molecular portrait of gene expression in these cells. There are several ways of identifying TISCs in the pool of cancer cells [49]. One of them is based on the analysis of cell surface markers specific for this type of cells [6]. To date, no universal cell surface marker has been identified that could be used to identify TISCs in different cancers. Another way of detecting TISCs is based on their ability to efflux lipophilic dyes, such as rhodamine or Hoechst. It is these two methods that are used for isolating TISCs in the numbers required for conducting either microarray or RNAseq analyses and for retrieving the molecular portrait of this type of cells. Using these two approaches, molecular portraits of gene expression of several types of TISCs were created [915].

In our studies [1620], we characterized a new marker of poorly differentiated cells of different origin, including TISCs, – the ability to internalize fragments of extracellular double-stranded DNA (such as plasmid DNA or FITC/TAMRA-labeled PCR products) by a natural internalization mechanism.

It was shown on the model of Krebs-2 ascites tumor that TAMRA+ cells trigger growth of a new tumor with similar histological characteristics. Elimination of these cells results in elimination of the engraftment potential of the cells and in curing of the mice from developed Krebs-2 ascites [16, 20].

In this study, we used the above marker to identify a set of genes overexpressed in TAMRA+ Krebs-2 cells. Based on this, a molecular portrait of gene expression in TAMRA+ Krebs-2 tumor cells was obtained, which firmly positions these cells as TISCs.

RESULTS

Internalization of TAMRA-labeled DNA fragments by Krebs-2 cells, RNAseq

It has been previously demonstrated that Krebs-2 ascites cells internalize TAMRA-labeled DNA probe by a natural internalization mechanism [16, 20]. 0.1 to 3% of Krebs-2 cells are naturally capable of capturing extracellular double-stranded DNA fragments (Figure 1A, 1B). Using this property, the populations of cells containing TAMRA-labeled DNA and TAMRA– cells were isolated from TAMRA-DNA treated Krebs-2 cells by repeated sorting (BD FACSAria, Becton Dickinson), with TAMRA+ target cells becoming 50% to 80% enriched and with nearly 100% pure TAMRA – cells (Figure (Figure1C1C).

Figure 1
Analysis of internalization of a TAMRA-labeled DNA into the Krebs-2 cells (Axioskop 2 Plus, Zeiss)

To obtain and analyze RNA transcripts, the RNAseq was applied. As a result, 2 million reads were obtained.

Trimming, analysis and mapping of reads

Figure Figure22 shows representative snapshots and illustrates the general approach to visual analysis of mapping data that will be further referred to as “expression profiles”. Since it was practically impossible to isolate the pure population of TAMRA + cells in the amount sufficient for RNA isolation and subsequent RNA sequencing (10- to 12-fold enrichment of this population could only be achieved; additional sorting rounds drastically decreased the cell survival and, as a consequence, RNA availability), while the pure TAMRA – population was available almost unlimitedly, we selected only those genes which were exclusively expressed in TAMRA+ cells.

Figure 2
(A) The general view of the expression profiles for TAMRA+ (blue profile) and TAMRA– (green profile) cells overlaid to melanoregulin gene (Mreg, uc007bkg.2, chromosome 1, – strand) marking. The gene consists of 5 exons visualized as yellow ...

Following primary visual inspection, we obtained a set of approximately 300 genes. For these genes, we retrieved the appropriate mRNA sequences, which were further used as target sequences for re-mapping of the total reads, both from TAMRA+ and TAMRA– datasets. In total, 3 mapping rounds were performed with sequentially more stringent specificity thresholds, which allowed us to filter out a significant number of genes that did not pass the test. This step was necessary due to a small number (2 millions) of reads obtained, which in turn was the consequence of previously described difficulties. Such a small number of mapped reads is supposed to produce a number of “homeless” reads, which can occasionally be mapped to the wrong place, and thus give a false positive or false negative signal.

As a result, we obtained a set of 168 genes that were expressed in TAMRA+ cells and not expressed (or expressed at an extremely low level) in TAMRA– cells. Further, these genes were subjected to qPCR analysis for additional verification.

Systemization of genes and their breakdown into functional groups

Supplementary Table S1 contains a list of 168 genes with their GenBank accession numbers, the frequency of reads in TAMRA– and TAMRA+ cells, and qPCR expression data for select genes (TAMRA+ vs TAMRA– cells). In addition, for some genes the information on the type of protein and subcellular localization in accordance with the classification of the Ingenuity Pathway Analysis (IPA) is provided.

To determine whether a gene is involved into molecular pathways determining the basic biological features of a cell, two different approaches were used.

The first approach was based on the data obtained in the study of Dolgova et al. [16], which indicated that TAMRA+ cells possessed the features of TISCs. This in turn indicated that gene ontology (GO) terms related to “stemness” and “cancer” should be overrepresented among the genes specific for TAMRA+ cells. Whether this was indeed the case was tested by matching their properties characterized in the original papers against the above GO categories.

“Stemness” genes

A stem cell is characterized by two features: the ability to divide asymmetrically and the ability to develop into any types of cells of the body (the pluripotency feature), transmitting this property to one of the daughter cells throughout many acts of cell division.

Asymmetric division of the stem cells is ensured by the HH, NOTCH and WNT pathways [2127]. The pluripotent status of the stem cells is primarily maintained via retinol signaling system [28]. Thus, to test the “stemness” the genes specifically expressed in TAMRA+ cells, they were considered in terms of their participation in pluripotency maintenance and asymmetric division.

Asymmetric division (group one) – Asb4, Nfatc2, Ppap2b, Tcf7l2, Uty, Wnt5a, Dnaja4, Rab37, Thpo, Sox9.

Pluripotency maintenance (the molecular factors of the regulation system of retinol) (group two) – Zfp39, Abca1, Abca9, Abca13, Aldh1a1, Aldh1l1, Bmper, Crabp2, Cyp2d26, Cyp26a1, Eif2s3y, Igf1, Igf2, Pik3ip1, Rragd, Rbp7, Rrh.

Other factors of stem cells (group three) – Clec11a, Lhx4, Gtf2a1l, Tbata, Mreg, Slc43a3, Ifi27l2b.

Altogether, 34 genes were found in the group of genes determining the “stemness” properties of the TAMRA+ cells. The ability of cells to undergo asymmetric cell division is characterized by the activity of the genes of the WNT-dependent signaling cascade (10 genes). Maintenance of the pluripotent state throughout numerous divisions is determined by the signaling cascade of the retinol system (17 genes). There are genes in each group that are known as key regulators of these cascades. For group one, these include the transcription factors Sox9, Nfatc2, Tcf7l2, Ppap2b, and secreted cytokines Wnt5a, Thpo. For group two, these are growth factors Igf1 and Igf2. Besides, TAMRA+ cells are characterized by the expression of genes belonging to other molecular systems of stem cells (group three). These are the growth factor Clec11a, the transcription factor Lhx4, and the proteins Gtf2a1l, Tbata, Mreg, Slc43a3, Ifi27l2b.

“Cancer” group genes

TISCs possess the following unique properties allowing them not to depend on morphogenetic laws determining the orderliness of the cell systems of a body. When dividing, TISCs form two types of cells. A cell is formed, completely identical with the original totipotent cells, plus a committed cell, capable of multiple, though not infinite, division. A TISC is able to induce development of a new tumor characterized by the same histological features as the original tumor [2931]. Further, TISCs are capable of metastasizing [3234]. The ability of TISCs to metastasize is determined by the property of these cells to undergo epithelial-mesenchymal transition [3539]. Remodeling of the cytoskeleton occurs, and factors altering the structure of the extracellular matrix in the cellular environment become expressed. A cell is thus able to cross the basal membrane and the endothelial wall, to move with the blood flow, to extravasate, and to form a new tumor nest. In this process, a TISC consistently activates the respective gene systems, acquires the property of escaping the attack of the immune cell system [40] and activates antiapoptotic genes [4144]. In addition, TISCs possess other unique qualities, such as the ability to eliminate and neutralize toxins [4549], to interact with the cells of the surrounding stroma and to stimulate them to transition to anaerobic respiration (glycolysis), which is known as the Warburg effect [5052]. Regulation of growth and differentiation of a tumor cell is also related to metabolism of Ca2+ and is determined by calcium-dependent protein kinase (PKC). Normally, PKC functions as a modulator and balances the growth and differentiation processes. Tumor cells are invariably characterized by PKC overactivation, which in turn induces proliferation, stimulates formation of phosphotyrosine, and enhances uncontrolled multiplication of cells. PKC is known to be activated by 3-phosphoinositol-dependent protein kinase (PDK) [5355]. In addition, apoptotic blebs of cancer cells were described to fuse together around the main body to form so-called blebbishields that have the properties of a TISCs [56]. Each of the above properties is characterized by concerted activity of a certain set of genes. Just as above, the selected gene set was analyzed for possible involvement into cancer based on the literature data.

Oncogenes, oncosuppressors – Eef1a2, Per2, Rasgrp3, Sash1, Tal1.

Genes characteristic of cancer cells – Oit3, Tnxb, Zfr2, VpreB3, Trim40.

Suppression of immunity – Arg2, Vsig4.

Multiple drug resistance, detoxification – Abca1, Abca13, Aldh1l1, Gstm3.

Anti-apoptotic genes – Cd5l, Tnfrsf13c.

Cytoskeleton reorganization – Myo1b, Pdlim4, Ttll2, Tubb1, Trpv4, Upk1b, Cobl, Cp, Dock10.

BM reorganization, adhesion, mobility in tissues – Pde4d, Prg4, Nrcam, Perp, Comp, Igsf3, Pvrl1, S100a14, Adamts2, Acpp, Fcna, Fblim1, Dsc2, Col6a2, Col3a1, Hepacam, Itga9, Mmp2, Ltbp1.

Dispersal via blood flow, exit from blood vessels – Selp, Lyve1, Cldn1.

Angiogenesis – Pdzd3 (Pdzk1).

Survival in tissues – Wnk2, Gas6.

Calcium metabolism – Pdk4.

In total, 54 genes fell into the “Cancer” category, based on the literature data.

Genes related to the cAMP-dependent signaling pathway

Analysis of the genes differentially expressed in TAMRA+ cells has allowed us to identify a group of genes that one way or another are associated with the cAMP-dependent signaling pathway. This is one of the most ancient (if not the most ancient) signalling systems, which is broadly involved in a variety of cellular processes, such as cell survival and apoptosis evasion, immediate early response to mitogens, cell differentiation and dedifferentiation, metabolism of lipids and carbohydrates. The reconstructed network (Figure (Figure3)3) includes more than 20% (36 out of 168) of the identified cancer stem cell-specific genes that, inter alia, provide a sufficient positive feedback to keep the system active (Supplementary Table S2): Nt5e, Acpp, Pde4d, Pde7b, Igf1, Igf2, Pf4, Gpha2, Prok2, Nts, Wnt5a, Nppa, Adrb3, Fgfr1, Gpr97, Gpr128, Cacna1d, Cp, Pdk4, Riiad1, Trpv4, Abca1, Cd200, Gata6, Amy1, Nfatc2. Kcnq2, Itm2a, Pon1, Tal1, Per2, Slc2a4, Tdo2, Mmp2, Alox15, Mst1.

Figure 3
The cAMP-dependent signaling network

All of the genes are conditionally divided into 3 groups: 1. upstream genes, that regulate cAMP metabolism and activity of the signaling core; 2. genes of the signaling system itself that include cAMP-metabolic enzymes and PKA complex; 3. downstream genes that are affected by the cAMP-signaling system. Genes with a positive feedback are shown in red. Most of the upstream genes affect the adenylate cyclase activity, but there are a number of them that affect the PKA and the related system itself. It is interesting to note that Wnt5a utilizes both ways: it triggers cAMP elevation at the plasma membrane and is implicated in increasing the catalytic subunits of PKA. Downstream genes, in turn, form 2 groups: 1) functionally activated ones that include transcription factors (like Nfatc2 that is functionally activated by PKA-dependent glycogen synthase kinase-3β inactivation) or any cellular function effectors (like Kcnq2 that is functionally activated by PKA-dependent phosphorylation) and 2) transcriptionally activated, which also include transcription factors (like Per2) and molecules that participate in some cellular process (like Pon1). Cellular processes, regulated by the downstream, include K+ (Kcnq2) or Ca2+ (Trpv4) transport, lipid (Abca1) and carbohydrate (Amy1) metabolism, as well as cell differentiation and commitment (Itm2a) and many others.

The remaining 55 genes – Alas2, Amac1, Ankrd22, Ap1m2, Blnk, Catsperg2, Ccdc8, Ccr3, Cd55, Cplx1, Crabp2, Ddx3y, Dusp13, Fam107a, Fmnl2, Gdf6, Ggt7, Gng11, Gpr97, Gtpbp10, Hpn, Iglon5, Il10, Il17rb, Inhbe, Jsrp1, Lass4/Cers4, Lpcat2b, Maged2, Marco, Mycbpap, Nppa, Pdlim3, Ppfia4, Psd2, Pydc4, Pyroxd2, Rab15, Sec14l4, Sec31b, Sel1l3, Serpinb1a, Serpinb2, Slc40a1, Slc44a3, Slco4a1, Srms, Tdo2, Timp3, Tmem59l, Tmem82, Tnn, Upk3b, Vsig8, Xlr3a, have also been annotated according to their functions described in the original publications, without binding them to the molecular systems determining the properties of TISCs. In this group, genes can be identified that participate in metabolic pathways as activators or repressors of molecular cascades of these pathways. These include growth factors (Gdf6, Inhbe), a kinase (Srms), a transcription factor (Lass4/Cers4), and cytokines (Il10, Il17rb). A number of these genes are insufficiently characterized in terms of their function and are of great interest as they may represent potential targets in anti-TISC therapy (at least in the context of Krebs-2 cancer cells).

It is to be noted that more than 40 genes characteristic for TAMRA+ TISC Krebs-2 (distributed across all the groups) are in the list of the core set of genes implicated in epithelial-mesenchymal transition of the cancer cells [3, 38] (Supplementary Table S2). The list includes two secreted trigger molecules – Wnt5a and Igf that are key to this process.

Real Time PCR verification of differential gene expression data

To validate the results obtained in the RNAseq experiments, we performed qPCR on cDNA synthesized from polyA+ mRNA of TAMRA+ and TAMRA– cells. Expression of the main genes representative of the categories of interest was characterized. The results of this analysis are shown in Figure Figure44 and are represented as fold increase in expression in TAMRA+ cells vs TAMRA– cells.

Figure 4
Real Time PCR validation of gene expression data of select genes identified in RNAseq

The analysis performed confirmed the results of the RNAseq and allowed a group of genes to be identified that are overexpressed in Krebs-2 cancer cells. In this group, two pairs of genes stand out: the secreted growth factor Igf1 and the transcription factor Nfatc2 activated by it, and cytokine Wnt5 and its downstream target transcription factor Tcf7L2 (Figure (Figure55).

Figure 5
(A) Distribution of all gene expression of TAMRA+ Krebs-2 cells in qPCR. (B) List of 22 genes whose expression in TAMRA+ cells relative TAMRA– cells was maximal in qPCR.

WNT5 is known to be a trigger molecule of the WNT5-dependent signaling pathway, while the transcription factor TCF712 activated as a result of triggering the WNT signaling cascade launches transcription of the genes of a genetic network determining the stemness properties of the TAMRA+ Krebs-2 cells [5760]. In its turn, IGF2 is a trigger molecule of the MAPK signaling cascade, where the signaling converges on the transcription factor NFATC2 that induces transcription of the genes from a genetic network determining the cancer properties of cells [6163]. At the same time, the list of over-expressing genes does not include any of the intermediate factors of the indicated signaling pathways. We believe TISCs may control the maintenance of their stem and cancer properties in an autocrine fashion, by permanently up-regulating the expression of these molecules. Secreted factors may ensure spatial protection from different extracellular regulatory molecules or other factors, in particular, regulatory exosomes. Increased local concentration of such factors may thus form a “protective shield” around a cell and maintains permanent extracellular trigger signaling [6466]. High expression of specific transcription factors ensures mandatory activation of a strictly determined genetic network that shapes the phenotype of a TISC. It is noteworthy that NFATC2 activates the genes encoding secreted growth factors and cytokines (Wnt5, Igf1, 2, Il-10) and the genes of multifunctional transcription factors (Gata6, Tal1, Tcf7L2). Simultaneously, Nfatc2 undergoes autoactivation (Supplementary Table S3). Such a pattern of interactions between the factors analyzed suggests that both pairs of proteins (IGF1-NFATC2, WNT5-TCF7L2), IL10, GATA6 and TAL1 are interconnected in the shared activation network and so they trigger and control several signaling cascades.

In the total list of the overexpressed genes, two groups are most noticeable, constituting more than half of the list and determining the stemness features of the cells – non-symmetric division and pluripotency. Increased expression of such genes may reflect the ontological necessity of the existence of this type of cells, primarily control and preserve the stemness properties, which provide a cell with a variety of advantages compared to a committed somatic cell.

Comparative functional analysis of the genes specifically overexpressed in TAMRA+ Krebs-2 TISCs across database resources

Next, the gene set found was compared against current pathway/functional databases to evaluate their possible contribution to the regulation of Krebs-2 TISCs. To assess the functional associations, analysis was conducted using Ingenuity Pathway Analysis (IPA) and KEGG platforms that cluster the genes into functional groups, and using OPOSSUM system to analyze the potential targets of transcription factors.

When fed with 168 genes overexpressed in TAMRA+ Krebs-2 TISCs, IPA grouped them into 4 main genetic networks (Figure (Figure6).6). As expected, in agreement with the manual analysis discussed above, these genes are primarily involved in cell functions such as organism development, cell morphology, and contribute to various molecular systems determining the biological properties of a cancer cell.

Figure 6
Hypothetical interaction networks of 168 genes expressed in TAMRA+ Krebs-2 TISCs as assessed by IPA

The four identified gene networks are interconnected by functional molecules (kinase PDK4, phosphatase ACPP, and the transcription factor TCF7L2) and so they may constitute a common system of functional associations in TAMRA+ Krebs-2 TISCs (Figure (Figure6B6B).

In the total list of genes processed by IPA, groups of genes were identified that mediate, trigger, or block certain gene programs. The following categories were found: transcription factors, kinases, phosphatases, growth factors, and cytokines (Table (Table11).

Table 1
Groups of genes involved in triggering or blocking of certain gene programs

It turned out that the gene networks built by IPA encompass general cellular transcription regulators (Nfatc2, Tal1, Gata6, Tcf7l2+) as well as growth factors and kinases, that are known as key inducers of many important cellular signaling systems (Fgfr1, Igf1,2, Mst1, Wnt5).

OPOSSUM analysis produced a list of gene targets of Nfatc2, Sox9, Tal1/Gata1, Tal1/Tcf3, and Tcf7l2+ transcription factors. This analysis suggests 69 genes out of 168 are under control of these factors, including own genes and the genes of other transcription factors (Supplementary Table S3). Noteworthy is activation by regulating the genes of the growth factors and cytokines, determining the specific functions of a cell and/or participating in maintenance of these functions (Igf1,2, Mst1, Il10, Wnt5a, Gdf6), which, in their turn, activate the said transcription factors, forming a loop of the function activation.

KEGG analysis indicates that several critical signaling pathways may be active in the TAMRA+ Krebs-2 TISCs. These signaling systems determine both stem and cancer phenotypes of these cells. Two gene pairs, Igf2-Nfatc2 and Wnt5-Tcf7l2+, which are the main players in MAPK- and WNT-dependent signaling systems, were again revealed. The first gene pair is known to be activated in cancer cells, while the second one determines the stem properties of a cell. In addition, NFAT2C is the transcription factor of such signaling cell systems as cAMP and VEGF. One can also mention yet another interesting pair of up-regulated genes in TAMRA+ Krebs-2 TISCs that are involved in mTOR-dependent signaling pathway – these are Igf1,2 and Eef1a2. That activity of this gene cascade is increased in Krebs-2 TISCs is supported by the fact the largest fold difference in expression of Igf1,2 genes in TAMRA+ vs TAMRA– cells. As noted previously, in TAMRA+ Krebs-2 TISCs, activation of genes that encode other growth factors/cytokines and transcription factors and that do not appear to be involved in a common pathway, is detected by KEGG analysis. These include Mst1, Il10, Gata6, and Sox9 genes, which are the functional components of the cGTP-, cAMP-, and FOXO-dependent signaling cascades.

DISCUSSION

Since their discovery, TISCs have been actively investigated in many laboratories of the world. In addition to the importance of academic knowledge of the specific features of this type of cells, scientists are interested in the search for molecular factors distinguishing these cells from normal stem cells and from committed cells. It is assumed that, knowing certain discriminating molecules that are on top of signaling cascades determining the properties of TISCs, one can target and block these factors, thus destroying the signaling pathways and eliminating the cells from the pathological tissue, improving the therapeutic prognosis [6769]. Currently investigations are being conducted of the molecular pattern of RNA expression of genes of different types of cells, including TISCs, in the hope to determine either the molecular target common for all the TISCs types or the molecular factor specific for only this type of TISC. To make a TISC preparation, either their ability to efflux lipophilic dyes, such as rhodamine or Hoechst, or fluorescent antibodies for a characteristic surface marker are used.

The TISC status is known to be determined by the activity of the genetic platforms activated by functionally tied signaling cascades forming stem and cancer properties of these cells. The main molecules of any characterized signaling pathways are represented by a functional pair – an extracellular trigger ligand, to which primarily secreted growth factors and cytokines (and their receptors) belong, and specific transcription factors activated by them, which trigger a gene platform, which determines the specific properties or the function of a cell. Other intermediate and finite factors mostly determine signal transduction, or connection with the other molecular cascades or the specific functional response of the cell.

Analysis of the genes overexpressed in the TAMRA+ cells has allowed us to propose the idea of autocrine control of the main qualities of TISCs due to increased expression of the specific trigger secreted factors launching signaling cascades, which determine the stem and cancer properties of TISCs and manifesting themselves in activation of the specific genetic platforms by ontologically related transcription factors, also characterized by increased expression.

We assumed this property to be characteristic of other types of TISCs. We analyzed the set of genes obtained in microarrays and by RNAseq analysis of 4 transcriptome libraries of TISCs (Table (Table2,2, Supplementary Table S2).

Table 2
Comparative analysis of the functional groups of genes in transcriptome libraries obtained for cancer stem cells of several types of cancer

A gene library [9] was obtained for the cancer stem cells of the human liver carcinoma. For this type of cancer stem cells, the following ontological paired trigger molecules and TF are detected: extracellular activator/stabilizer IGF/IGFBP5 as a marker of the extracellular secreted growth factor [70] and the transcription factor FOSL1 [71], activating the specific gene platform and, in particular, triggering the synthesis of MMP9, which leads to reorganization of the extracellular matrix. IGF/IGFBP5 is a component of FOXO, mTOR and MAPK/ERK signaling cascades, FOSL1 is in the system of the mTOR and MAPK/ERK signaling pathways [72]. Both factors are responsible for activation of the gene network characterizing the cancer properties of a cell. Activity of an ontological pair determines cell invasion by remodeling the structures of the extracellular matrix, which is one of the characteristics of TISCs. In this type of cells, activated genes are associated with such biological functions as multiple drug resistance, fat metabolism, angiogenesis, proliferation, transport, response to inflammation, production of cytokines, a cell cycle, and signal transduction. Observed are suppression of the expression of genes involved in development of organs, response to hypoxia, fixation of nucleotides, translation and elongation, and hemotaxis.

The SP OCC transcriptome (the side population cell from an ovarian cancer cell line) was described in [12]. This list includes the transcripts of 5 main secreted factors, which participate in biological cellular processes determining the cancer properties of cells. These are LTBP1, LIF, NRG1, HB-EGF, IL6. Transcription of these genes suggests the active state of signaling pathways as PI3K/AKT [73], CSF1 [74], AKT/mTOR [75], TGFβ1 [76]. In addition, indicating the activity of the VEGF, EGF, WNT signaling cascades is expression of the genes of the respective receptors NRP2 [77], EGFR [78], FZD5, 7, 8 [7982]. It is the FZD receptors that are required for maintaining the pluripotent state of the cells. Besides, the list of genes contains transcription factors connected with secreted molecules by functional ties. The NFRB, TLE, NFAT5, RA\UNX, NKX3.1, NR4a2 transcription factors are identified, which are active molecules of the FOXO-, ERK-, MAPK-, WNT-dependent signaling pathways. The PLAGL2 transcription factor is present, an expressed suppressor of differentiation of stem cells and of stem cancer glial cells. At the same time, PLAGL2 determines their ability of self-renewal expressed through modulation of the WNT/catenin-dependent signaling pathway [83]. It was also shown that the set of genes characterizing the activity of the NOTCH cascade is at a high level of expression. The activity of this signaling pathway also determines the ability of self-renewal as one of the characteristics of stemness of this type of cells. Thus, trigger secreted molecules and transcription factors under their control, which determine the stem and cancer properties of cells, are the present in the SP OCC transcriptome, too.

In total, genes were clusterized characteristic of such biological processes as the cell cycle, apoptosis, proliferation, transport, signal transduction, transcription and translation, metabolism and modification of protein in the type of cells analyzed. Such a set of activities characterizes high metabolic activity of cells, which is a characteristic feature of cancer cells.

Birnie and colleagues evaluated the stem and cancer cell expression profile for the human prostate cancer [84]. Despite the small number of the characterized genes, secreted trigger molecules and ontological transcription factors, determining activation of the specific gene platforms, are also identified in the transcriptome obtained. Secreted factors are represented by the WNT5A, IL6, CSF2, IFNK u IFNGR/1, and PAPPA genes [8588]. There are two transcription factors STAT1 and NFKB1 and the modulator of many transcription factors BTG1 in this list [8991]. The following ontological functional ties can be seen between the factors. Secreted IL6, IFNK and IFNGR/1 cytokines activate two main signaling pathways JAK/STAT and MAPK and are connected with the classical ontological partner STAT1. Phosphorylation of STAT1 is accompanied by transcription of the platform of the IF/IL6-induced genes, which characterizes development of an immune response. In addition, activation of JAK/STAT and MAPK results in activation of genes related to metastasizing, apoptosis, and proliferation, indicating participation of an ontological group in functioning of the cancer cell [86, 88, 89, 92, 93]. PAPPA and WNT5 are connected with the NFKB1 transcription factor by ontological ties. Paracrine activation of the NFKB-dependent signaling pathway may be accompanied both by an immune response (WNT5 [60]) and changes in signaling cascades determining cell survival and maintenance of the stability of its genome (PAPPA [87]). WNT5, as a trigger factor, directly activates the WNT-dependent signaling pathway determining one of the main properties of the stem cells, non-symmetric division. CSF2–GM-CSF is one of the factors required for proliferation, differentiation and survival of macrophages and granulocytes. In addition, its participation in the early embryo development in the sustenance of gestation has also been demonstrated [9496], which may characterize the stem properties of the cancer cells of the human prostate cancer. CSF2 activates AKT1-, ERK1/2-, mTOR- but not JAK/STAT-dependent signaling cascades, which, as mentioned previously, are active in many cancer stem cells. Thus, the transcriptome analysis of the cancer stem cells of the human prostate cancer is in agreement with the idea that cancer stem cells are expressed by a small set of trigger molecules and by the transcription factors connected with them by functional ties, which maintain and protect their stem and cancer properties from the impact of the surrounding factors.

The combined analysis in the KEGG system indicates the activity of the clusters of genes determining mobility of the caner stem cells of the human prostate cancer. The group of genes of focal adhesion and the genes coding the proteins destroying matrix-receptor interaction have been identified.

Using an orthotopic mouse model, the pattern of gene expression profiling of the TISCs isolated from human osteosarcoma cell lines has been analyzed [10]. The following secreted trigger molecules IGF1, IHH, and NFIX, NFYB, SPI1, TFDB1 transcription factors have been identified in the expression profile obtained. All these factors are connected in one regulatory cluster, in which both the stemness and “cancer” properties of cells are determined. IGF1-IGFR1 is required for expression of IHH, a member of the family of HH mammal genes. The IGF-dependent signaling pathway modulates the expression of IHH, resulting in the control over cell proliferation, the survival, differentiation, and production of matrix in the growth plate in the course of post-natal development [97]. At the same time, IHH is connected with the WNT-dependent signaling cascade and activates it, reducing the level of the WNT antagonists (LRP, SFRP -) [98, 99]. Abnormal activation of the IHH-dependent signaling pathway due to increased expression of IGF leads to development defects and transformation of stem cells in TISCs, resulting in cancer development. The IGF-dependent signaling pathway is a component of the FOXO-, mTOR- and MAPK-dependent signaling cascades, the impaired activity of which also results in the carcinogenic transformation in the stem cells. Thus, abnormal expression of two trigger molecules IGF/IHH may be one of the causes of cancer development and of maintenance of the cancer properties of TISCs [99, 100].

Ontological interactions among the transcription factors analyzed have been established. SPI1 (PU.1) is the key regulator of hematopoiesis. Reduction of the expression of this transcription factor results in development of acute myeloid leukemia [101]. It is functionally associated with IHH. Joint expression of SPI1 and IHH results in the appearance of self-renewing leukemic cancer stem cells erythroid progenitors are infected with viruses with subsequent occurrence of leukemia [102]. E2F1 induces expression of NFYB, which together with E2F1 regulates the expression of many target genes. Excessive expression of these genes is found in osteosarcoma and is associated with non-sensitivity to chemotherapy. The NFYB/E2F1 complex becomes connected with FOXO1/3 and is the factor of the FOXO- dependent signaling cascade, which is activated in the cells of different cancers [103]. Another transcription factor, NFIX, also physically interacts with FOXA1 and regulates functioning and survival of hematopoietic stem cells and is also in the chain of signal transduction of the FOXO-dependent signaling cascade [104, 105]. The TFDP1 transcription factor is required to control the performance of a cell cycle and expression of target genes. GATA1 and TAL1 regulate the expression of TFDP1. Amplification of this gene results in deregulation of the cell cycle and in cancer development [106, 107]. Thus, the ontological pair SPI1 and IHH controls the primary property of a TISCs – self-renewal. The complex of the NFYB/E2F1- FOXO1/3 – NFIX transcription factors and the TFDP1 transcription factor regulate such properties of a TISC as aberrant proliferation, drug resistance, and survival. The combined comparison in the KEGG system indicates that the stem cell genes of the human osteosarcoma are clustered in the groups characteristic of two biological processes, development of skeletal bones and mobility (migration). No relevant differences in the expression of the known markers of Oct/Nanog embryo stem cells have been found.

Analysis of 5 transcriptomes of cancer stem cells made indicates that 1) the genetic platforms of different types of cancer stem cells differ; 2) there is no equally proportional simultaneous expression of all the genes of a certain signaling pathway; 3) TISCs autocrinely controls constancy of its stem and cancer properties, maintaining a permanently high level of expression of secreted trigger molecules, which maintain the permanent extracellular trigger signal. The high expression of the specific transcription factors of the ontological partners of trigger molecules ensures mandatory activation of a strictly defined gene platform, which forms the specific properties of a TISCs.

An inference from the analysis made is the idea that, in order to destroy the basic properties of cancer stem cells, it is necessary to eliminate the secreted fraction of the trigger molecules and of their ontological transcription partners immediately and simultaneously and to stop expression of the genes of all the participants of the manifested factors which determine “stemness’ and “cancer” features of this type of TISCs.

MATERIALS AND METHODS

Isolation of Krebs-2 ascites cells internalizing TAMRA-labeled DNA

Alu-TAMRA DNA is a 500 bp fragment of the human Alu repeat labeled with dUTP-5′-TAMRA by PCR [16]. Krebs-2 cells were incubated with Alu-TAMRA DNA (0.5 μg per 1 million cells in 200 μl RPMI-1640 medium) for 1 hour at room temperature. Separation of TAMRA+ and TAMRA– subpopulations was performed on a BD FACSAria cell sorter equipped with 100 mkm nozzle at the event rate varying between 2000 to 3500 cells. To detect TAMRA fluorescence, a 488 nm laser and a 585/42 channel (also known as a PE channel) were used.

Several rounds of flow cytometry-based cell sorting were used to obtain enough amounts of TAMRA+ cells. Due to the low content of TAMRA+ cells in the total Krebs-2 population (typically below 5%), the isolation of a reasonable number of pure TAMRA+ cells was nearly impossible. The best result we managed to achieve was 80% of TAMRA+ cells after a sorting pass. Subjecting the cells to multiple rounds of re-sorting resulted in drastic reduction in cell viability. So our population subsequently named as “TAMRA+ cells” consist of 60–70% of TAMRA+ cells and 30–40% of TAMRA– cells. In contrast to TAMRA+ cells, the pure TAMRA– subpopulation was virtually unlimited.

After sorting, microscopy analysis of the percentage of TAMRA+ cells was performed using a laser scanning microscope LSM 780 NLO (Carl Zeiss, USA) in the collective use center for microscopy of biological objects, Siberian Branch of the Russian Academy of Sciences.

RNA isolation and RNAseq

Total RNA from TAMRA+ and TAMRA– cells was isolated using a Qiagen RNeasy kit. Libraries were constructed using Illumina TruSeq RNA sample preparation kit, v.2. The libraries were loaded on Illumina Hiseq2000 using Truseq v.3 reagents. After sequencing, the reads were sorted according to their indices with the CASAVA 1.8.2. software package (laboratory of evolutionary genomics, department of bioengineering and bioinformatics, Lomonosov Moscow State University, Moscow). Two independent DNA sequencing repeats produced 2 million paired-end reads 70 bp each, that were treated as single-end reads for mapping purposes.

Reads analysis and software

The obtained reads were trimmed out of adapter sequences using a Trimmomatic tool from the Usadel Lab (http://www.usadellab.org/cms/?page=trimmomatic) and further mapped to murine chromosomes (available at the UCSC Genome Bioinformatics: http://hgdownload.soe.ucsc.edu/goldenPath/mm10/chromosomes/). As a mapping tool, the Softberry's ReadsMap (Softberry.Inc) software [108] was used. The unique feature of this tool is that it maps a certain read to the appropriate contig (a chromosome in our case) by selecting the best alignment from the pool of those found. Combined with another unique feature that allowed us to determine the splicing sites with high accuracy and thereby to drastically minimize the number of read “tails” mismapped to introns (http://www.softberry.com/berry.phtml?topic=readsmap-i&group=help&subgroup=pipelines and http://www.softberry.com/berry.phtml?topic=rnaseq), the ReadsMap tool produces precise mapping profiles that reflect the real number of reads aligned to the given position in the genome and in fact can be considered as relevant to the gene expression level. The trimmed reads are available at the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena/data/view/PRJEB15164).

Read counts were smoothed and log2-scaled to facilitate visual analysis of genes having extremely high expression differences (10 times and more). As a visualization tool for comparing the profiles and testing their correspondence to the annotated genes (KnownGenes available at UCSC Genome Bioinformatics: http://hgdownload.soe.ucsc.edu/goldenPath/mm10/database/), the Softberry's Sequence Explorer v.2.2.18 (http://www.softberry.com/) was used.

qPCR primer selection

Krebs-2 ascites tumor cells were sorted into TAMRA+ and TAMRA– subpopulations and RNA was isolated. CDNA was synthesized for each group using a GoScript Reverse Transcription System (Promega, USA). A pair of primers was selected for each gene, which corresponded to the area (exon) of the highest read density. The primers were tested for specificity using the total Krebs-2 cDNA as a template. Real-time qPCR was run on an Applied Biosystems instrument using the SYBR® Green PCR Master Mix (Applied Biosystems, USA) reagent in three replicas. Murine β-actin gene was used as a reference.

SUPPLEMENTARY MATERIALS TABLES

Acknowledgments

We express our gratitude to our colleagues from the laboratory of evolutionary genomics, the faculty of bioengineering and bioinformatics, Moscow State University, Moscow, for performing RNAseq.

The authors express their gratitude to Tatyana Padve for translating the paper and Andrey Gorchakov for smart critical comments.

Footnotes

CONFLICTS OF INTEREST

The authors declare that they have no competing interests.

GRANT SUPPORT

The work was funded by the LLC “BA-pharma”, LLC “Panagen”, RFBR grant N 15-04-03386/16 and the State scientific project N 0324-2015-0003.

REFERENCES

1. Lapidot T, Sirard C, Vormoor J, Murdoch B, Hoang T, Caceres-Cortes J, Minden M, Paterson B, Caligiuri MA, Dick JE. A cell initiating human acute myeloid leukaemia after transplantation into SCID mice. Nature. 1994;367:645–648. [PubMed]
2. Bonnet D, Dick JE. Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell. Nat Med. 1997;3:730–737. [PubMed]
3. Gröger CJ, Grubinger M, Waldhör T, Vierlinger K, Mikulits W. Meta-analysis of gene expression signatures defining the epithelial to mesenchymal transition during cancer progression. PLoS One. 2012;7:e51136. [PMC free article] [PubMed]
4. Frisan T, Levitsky V, Masucci M. Limiting dilution assay. Methods Mol Biol. 2001;174:213–216. [PubMed]
5. Rajasekhar VK. Analytical methods for cancer stem cells. Methods Mol Biol. 2007;407:83–95. [PubMed]
6. O’Brien CA, Kreso A, Jamieson CH. Cancer stem cells and self-renewal. Clin Cancer Res. 2010;16:3113–3120. [PubMed]
7. Hu Y, Smyth GK. ELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays. J Immunol Methods. 2009;347:70–78. [PubMed]
8. Cao L, Zhou Y, Zhai B, Liao J, Xu W, Zhang R, Li J, Zhang Y, Chen L, Qian H, Wu M, Yin Z. Sphere-forming cell subpopulations with cancer stem cell properties in human hepatoma cell lines. BMC Gastroenterol. 2011;11:71. [PMC free article] [PubMed]
9. Ho DW, Yang ZF, Yi K, Lam CT, Ng MN, Yu WC, Lau J, Wan T, Wang X, Yan Z, Liu H, Zhang Y, Fan ST. Gene expression profiling of liver cancer stem cells by RNA-sequencing. PLoS One. 2012;7:e37159. [PMC free article] [PubMed]
10. Rainusso N, Man TK, Lau CC, Hicks J, Shen JJ, Yu A, Wang LL, Rosen JM. Identification and gene expression profiling of tumor-initiating cells isolated from human osteosarcoma cell lines in an orthotopic mouse model. Cancer Biol Ther. 2011;12:278–287. [PMC free article] [PubMed]
11. Wang H, Huang M, Zhang DY, Zhang F. Global profiling of signaling networks: study of breast cancer stem cells and potential regulation. Oncologist. 2011;16:966–979. [PMC free article] [PubMed]
12. Vathipadiekal V, Saxena D, Mok SC, Hauschka PV, Ozbun L, Birrer MJ. Identification of a potential ovarian cancer stem cell gene expression profile from advanced stage papillary serous ovarian cancer. PLoS One. 2012;7:e29079. [PMC free article] [PubMed]
13. Karimi-Busheri F, Zadorozhny V, Carrier E, Fakhrai H. Molecular integrity and global gene expression of breast and lung cancer stem cells under long-term storage and recovery. Cell Tissue Bank. 2013;14:175–186. [PMC free article] [PubMed]
14. Pang LY, Gatenby EL, Kamida A, Whitelaw BA, Hupp TR, Argyle DJ. Global gene expression analysis of canine osteosarcoma stem cells reveals a novel role for COX-2 in tumour initiation. PLoS One. 2014;9:e83144. [PMC free article] [PubMed]
15. Zakaria N, Yusoff NM, Zakaria Z, Lim MN, Baharuddin PJ, Fakiruddin KS, Yahaya B. Human non-small cell lung cancer expresses putative cancer stem cell markers and exhibits the transcriptomic profile of multipotent cells. BMC Cancer. 2015;15:84. [PMC free article] [PubMed]
16. Dolgova EV, Alyamkina EA, Efremov YR, Nikolin VP, Popova NA, Tyrinova TV, Kozel AV, Minkevich AM, Andrushkevich OM, Zavyalov EL, Romaschenko AV, Bayborodin SI, Taranov OS, et al. Identification of cancer stem cells and a strategy for their elimination. Cancer Biol Ther. 2014;15:1378–1394. [PMC free article] [PubMed]
17. Dolgova EV, Efremov YR, Orishchenko KE, Andrushkevich OM, Alyamkina EA, Proskurina AS, Bayborodin SI, Nikolin VP, Popova NA, Chernykh ER, Ostanin AA, Taranov OS, Omigov VV, et al. Delivery and processing of exogenous double-stranded DNA in mouse CD34+ hematopoietic progenitor cells and their cell cycle changes upon combined treatment with cyclophosphamide and double-stranded DNA. Gene. 2013;528:74–83. [PubMed]
18. Dolgova EV, Potter EA, Proskurina AS, Minkevich AM, Chernych ER, Ostanin AA, Efremov YR, Bayborodin SI, Nikolin VP, Popova NA, Kolchanov NA, Bogachev SS. Properties of internalization factors contributing to the uptake of extracellular DNA into tumor-initiating stem cells of mouse Krebs-2 cell line. Stem Cell Res Ther. 2016;7:76. [PMC free article] [PubMed]
19. Dolgova EV, Proskurina AS, Nikolin VP, Popova NA, Alyamkina EA, Orishchenko KE, Rogachev VA, Efremov YR, Dubatolova TD, Prokopenko AV, Chernykh ER, Ostanin AA, Taranov OS, et al. “Delayed death” phenomenon: A synergistic action of cyclophosphamide and exogenous DNA. Gene. 2012;495:134–145. [PubMed]
20. Potter EA, Dolgova EV, Proskurina AS, Minkevich AM, Efremov YR, Taranov OS, Omigov VV, Nikolin VP, Popova NA, Bayborodin SI, Ostanin AA, Chernykh ER, Kolchanov NA, et al. A strategy to eradicate well-developed Krebs-2 ascites in mice. Oncotarget. 2016;7:11580–11594. doi: 10.18632/oncotarget.7311. [PMC free article] [PubMed] [Cross Ref]
21. Taipale J, Beachy PA. The Hedgehog and Wnt signalling pathways in cancer. Nature. 2001;411:349–354. [PubMed]
22. Merchant AA, Matsui W. Targeting Hedgehog--a cancer stem cell pathway. Clin Cancer Res. 2010;16:3130–3140. [PMC free article] [PubMed]
23. Wilson CW, Chuang PT. Mechanism and evolution of cytosolic Hedgehog signal transduction. Development. 2010;137:2079–2094. [PubMed]
24. Ryan KE, Chiang C. Hedgehog secretion and signal transduction in vertebrates. J Biol Chem. 2012;287:17905–17913. [PMC free article] [PubMed]
25. Pez F, Lopez A, Kim M, Wands JR, C Caron de Fromentel, Merle P. Wnt signaling and hepatocarcinogenesis: molecular targets for the development of innovative anticancer drugs. J Hepatol. 2013;59:1107–1117. [PubMed]
26. Morell CM, Strazzabosco M. Notch signaling and new therapeutic options in liver disease. J Hepatol. 2014;60:885–890. [PubMed]
27. Purro SA, Galli S, Salinas PC. Dysfunction of Wnt signaling and synaptic disassembly in neurodegenerative diseases. J Mol Cell Biol. 2014;6:75–80. [PMC free article] [PubMed]
28. Khillan JS. Vitamin A/retinol and maintenance of pluripotency of stem cells. Nutrients. 2014;6:1209–1222. [PMC free article] [PubMed]
29. Alison MR, Lin WR, Lim SM, Nicholson LJ. Cancer stem cells: in the line of fire. Cancer Treat Rev. 2012;38:589–598. [PubMed]
30. Baccelli I, Trumpp A. The evolving concept of cancer and metastasis stem cells. J Cell Biol. 2012;198:281–293. [PMC free article] [PubMed]
31. Peitzsch C, Kurth I, Kunz-Schughart L, Baumann M, Dubrovska A. Discovery of the cancer stem cell related determinants of radioresistance. Radiother Oncol. 2013;108:378–387. [PubMed]
32. Pavelic SK, Sedic M, Bosnjak H, Spaventi S, Pavelic K. Metastasis: new perspectives on an old problem. Mol Cancer. 2011;10:22. [PMC free article] [PubMed]
33. Valastyan S, Weinberg RA. Tumor metastasis: molecular insights and evolving paradigms. Cell. 2011;147:275–292. [PMC free article] [PubMed]
34. Wan L, Pantel K, Kang Y. Tumor metastasis: moving new biological insights into the clinic. Nat Med. 2013;19:1450–1464. [PubMed]
35. Chiarugi P. Escaping from, moving towards, following a path, squeezing through: lots of opportunities for moving cells. Cell Commun Signal. 2010;8:25. [PMC free article] [PubMed]
36. Cieślik M, Hoang SA, Baranova N, Chodaparambil S, Kumar M, Allison DF, Xu X, Wamsley JJ, Gray L, Jones DR, Mayo MW, Bekiranov S. Epigenetic coordination of signaling pathways during the epithelial-mesenchymal transition. Epigenetics Chromatin. 2013;6:28. [PMC free article] [PubMed]
37. Tsai JH, Yang J. Epithelial-mesenchymal plasticity in carcinoma metastasis. Genes Dev. 2013;27:2192–2206. [PubMed]
38. Venkov C, Plieth D, Ni T, Karmaker A, Bian A, George AL, Jr, Neilson EG. Transcriptional networks in epithelial-mesenchymal transition. PLoS One. 2011;6:e25354. [PMC free article] [PubMed]
39. Glushankova NA. Changes in the regulation of intercellular adhesion in tumor transformation: a review. [Article in Russian]. Biokhimiya. 2008;73:925–934.
40. Berezhnaya NM. The role of immunity cells in the tumor micro-environment. I. Cytokine cells – participants of inflammation. [Article in Russian]. Oncologiya. 2009;11:86–91.
41. Maeda S, Yoshida H, Mitsuno Y, Hirata Y, Ogura K, Shiratori Y, Omata M. Analysis of apoptotic and antiapoptotic signalling pathways induced by Helicobacter pylori. Mol Pathol. 2002;55:286–293. [PMC free article] [PubMed]
42. Al-Harbi S, Hill BT, Mazumder S, Singh K, Devecchio J, Choudhary G, Rybicki LA, Kalaycio M, Maciejewski JP, Houghton JA, Almasan A. An antiapoptotic BCL-2 family expression index predicts the response of chronic lymphocytic leukemia to ABT-737. Blood. 2011;118:3579–3590. [PubMed]
43. Gross KL, Oakley RH, Scoltock AB, Jewell CM, Cidlowski JA. Glucocorticoid receptor alpha isoform-selective regulation of antiapoptotic genes in osteosarcoma cells: a new mechanism for glucocorticoid resistance. Mol Endocrinol. 2011;25:1087–1099. [PubMed]
44. Halder UC, Bhowmick R, T Roy Mukherjee, Nayak MK, Chawla-Sarkar M. Phosphorylation drives an apoptotic protein to activate antiapoptotic genes: paradigm of influenza A matrix 1 protein function. J Biol Chem. 2013;288:14554–14568. [PMC free article] [PubMed]
45. Holohan C, S Van Schaeybroeck, Longley DB, Johnston PG. Cancer drug resistance: an evolving paradigm. Nat Rev Cancer. 2013;13:714–726. [PubMed]
46. Januchowski R, Wojtowicz K, Zabel M. The role of aldehyde dehydrogenase (ALDH) in cancer drug resistance. Biomed Pharmacother. 2013;67:669–680. [PubMed]
47. Housman G, Byler S, Heerboth S, Lapinska K, Longacre M, Snyder N, Sarkar S. Drug resistance in cancer: an overview. Cancers (Basel) 2014;6:1769–1792. [PMC free article] [PubMed]
48. Zhang YK, Wang YJ, Gupta P, Chen ZS. Multidrug Resistance Proteins (MRPs) and Cancer Therapy. AAPS J. 2015;17:802–812. [PMC free article] [PubMed]
49. Stavrovsskaya AA, Stromskaya TP. Transport proteins of the ABC family in multidrug resistance of human tumor cells. [Article in Russian] Biokhimiya. 2008;73:735–750.
50. Hsu PP, Sabatini DM. Cancer cell metabolism: Warburg and beyond. Cell. 2008;134:703–707. [PubMed]
51. Pavlides S, Whitaker-Menezes D, Castello-Cros R, Flomenberg N, Witkiewicz AK, Frank PG, Casimiro MC, Wang C, Fortina P, Addya S, Pestell RG, Martinez-Outschoorn UE, Sotgia F, et al. The reverse Warburg effect: aerobic glycolysis in cancer associated fibroblasts and the tumor stroma. Cell Cycle. 2009;8:3984–4001. [PubMed]
52. Bonuccelli G, Whitaker-Menezes D, Castello-Cros R, Pavlides S, Pestell RG, Fatatis A, Witkiewicz AK, Vander Heiden MG, Migneco G, Chiavarina B, Frank PG, Capozza F, Flomenberg N, et al. The reverse Warburg effect: glycolysis inhibitors prevent the tumor promoting effects of caveolin-1 deficient cancer associated fibroblasts. Cell Cycle. 2010;9:1960–1971. [PubMed]
53. Alessi DR, Deak M, Casamayor A, Caudwell FB, Morrice N, Norman DG, Gaffney P, Reese CB, MacDougall CN, Harbison D, Ashworth A, Bownes M. 3-Phosphoinositide-dependent protein kinase-1 (PDK1): structural and functional homology with the Drosophila DSTPK61 kinase. Curr Biol. 1997;7:776–89. [PubMed]
54. Alessi DR, Kozlowski MT, Weng QP, Morrice N, Avruch J. 3-Phosphoinositide-dependent protein kinase 1 (PDK1) phosphorylates and activates the p70 S6 kinase in vivo and in vitro. Curr Biol. 1998;8:69–81. [PubMed]
55. Le Good JA, Ziegler WH, Parekh DB, Alessi DR, Cohen P, Parker PJ. Protein kinase C isotypes controlled by phosphoinositide 3-kinase through the protein kinase PDK1. Science. 1998;281:2042–2045. [PubMed]
56. Jinesh GG, Choi W, Shah JB, Lee EK, Willis DL, Kamat AM. Blebbishields, the emergency program for cancer stem cells: sphere formation and tumorigenesis after apoptosis. Cell Death Differ. 2013;20:382–395. [PMC free article] [PubMed]
57. Fu H, Kesari S, Cai J. Tcf7l2 is tightly controlled during myelin formation. Cell Mol Neurobiol. 2012;32:345–352. [PMC free article] [PubMed]
58. Saito-Diaz K, Chen TW, Wang X, Thorne CA, Wallace HA, Page-McCaw A, Lee E. The way Wnt works: components and mechanism. Growth Factors. 2013;31:1–31. [PMC free article] [PubMed]
59. Zhao C, Deng Y, Liu L, Yu K, Zhang L, Wang H, He X, Wang J, Lu C, Wu LN, Weng Q, Mao M, Li J, et al. Dual regulatory switch through interactions of Tcf7l2/Tcf4 with stage-specific partners propels oligodendroglial maturation. Nat Commun. 2016;7:10883. [PMC free article] [PubMed]
60. Zhao Y, Wang CL, Li RM, Hui TQ, Su YY, Yuan Q, Zhou XD, Ye L. Wnt5a promotes inflammatory responses via nuclear factor κB (NF-κB) and mitogen-activated protein kinase (MAPK) pathways in human dental pulp cells. J Biol Chem. 2014;289:21028–21039. [PMC free article] [PubMed]
61. Church DN, Phillips BR, Stuckey DJ, Barnes DJ, Buffa FM, Manek S, Clarke K, Harris AL, Carter EJ, Hassan AB. Igf2 ligand dependency of Pten(+/−) developmental and tumour phenotypes in the mouse. Oncogene. 2012;31:3635–3646. [PMC free article] [PubMed]
62. Dai J, Buhusi M, Demyanenko GP, Brennaman LH, Hruska M, Dalva MB, Maness PF. Neuron glia-related cell adhesion molecule (NrCAM) promotes topographic retinocollicular mapping. PLoS One. 2013;8:e73000. [PMC free article] [PubMed]
63. Guillaud-Bataille M, Ragazzon B, de Reyniès A, Chevalier C, Francillard I, Barreau O, Steunou V, Guillemot J, Tissier F, Rizk-Rabin M, René-Corail F, Al Ghuzlan A, Assié G, et al. IGF2 promotes growth of adrenocortical carcinoma cells, but its overexpression does not modify phenotypic and molecular features of adrenocortical carcinoma. PLoS One. 2014;9:e103744. [PMC free article] [PubMed]
64. Yang C, Robbins PD. The roles of tumor-derived exosomes in cancer pathogenesis. Clin Dev Immunol. 2011;2011:842849. [PMC free article] [PubMed]
65. Camussi G, Deregibus MC, Cantaluppi V. Role of stem-cell-derived microvesicles in the paracrine action of stem cells. Biochem Soc Trans. 2013;41:283–287. [PubMed]
66. Abd Elmageed ZY, Yang Y, Thomas R, Ranjan M, Mondal D, Moroz K, Fang Z, Rezk BM, Moparty K, Sikka SC, Sartor O, Abdel-Mageed AB. Neoplastic reprogramming of patient-derived adipose stem cells by prostate cancer cell-associated exosomes. Stem Cells. 2014;32:983–997. [PMC free article] [PubMed]
67. Yang ZJ, Wechsler-Reya RJ. Hit‚ em where they live: targeting the cancer stem cell niche. Cancer Cell. 2007;11:3–5. [PubMed]
68. McDermott SP, Wicha MS. Targeting breast cancer stem cells. Mol Oncol. 2010;4:404–419. [PMC free article] [PubMed]
69. Insan MB, Jaitak V. New approaches to target cancer stem cells: current scenario. Mini Rev Med Chem. 2014;14:20–34. [PubMed]
70. Güllü G, Karabulut S, Akkiprik M. Functional roles and clinical values of insulin-like growth factor-binding protein-5 in different types of cancers. Chin J Cancer. 2012;31:266–280. [PMC free article] [PubMed]
71. Kent LN, Rumi MA, Kubota K, Lee DS, Soares MJ. FOSL1 is integral to establishing the maternal-fetal interface. Mol Cell Biol. 2011;31:4801–4813. [PMC free article] [PubMed]
72. Benavides-Serrato A, Anderson L, Holmes B, Cloninger C, Artinian N, Bashir T, Gera J. mTORC2 modulates feedback regulation of p38 MAPK activity via DUSP10/MKP5 to confer differential responses to PP242 in glioblastoma. Genes Cancer. 2014;5:393–406. doi: 10.18632/genesandcancer.41. [PMC free article] [PubMed] [Cross Ref]
73. Liang X, Ding Y, Zhang Y, Chai YH, He J, Chiu SM, Gao F, Tse HF, Lian Q. Activation of NRG1-ERBB4 signaling potentiates mesenchymal stem cell-mediated myocardial repairs following myocardial infarction. Cell Death Dis. 2015;6:e1765. [PMC free article] [PubMed]
74. Zhou ZN, Sharma VP, Beaty BT, Roh-Johnson M, Peterson EA, Van Rooijen N, Kenny PA, Wiley HS, Condeelis JS, Segall JE. Autocrine HBEGF expression promotes breast cancer intravasation, metastasis and macrophage-independent invasion in vivo. Oncogene. 2014;33:3784–3793. [PMC free article] [PubMed]
75. Li X, Yang Q, Yu H, Wu L, Zhao Y, Zhang C, Yue X, Liu Z, Wu H, Haffty BG, Feng Z, Hu W. LIF promotes tumorigenesis and metastasis of breast cancer through the AKT-mTOR pathway. Oncotarget. 2014;5:788–801. doi: 10.18632/oncotarget.1772. [PMC free article] [PubMed] [Cross Ref]
76. Ge G, Greenspan DS. BMP1 controls TGFbeta1 activation via cleavage of latent TGFbeta-binding protein. J Cell Biol. 2006;175:111–20. [PMC free article] [PubMed]
77. Goel HL, Chang C, Pursell B, Leav I, Lyle S, Xi HS, Hsieh CC, Adisetiyo H, Roy-Burman P, Coleman IM, Nelson PS, Vessella RL, Davis RJ, et al. VEGF/neuropilin-2 regulation of Bmi-1 and consequent repression of IGF-IR define a novel mechanism of aggressive prostate cancer. Cancer Discov. 2012;2:906–921. [PMC free article] [PubMed]
78. Baumdick M, Brüggemann Y, Schmick M, Xouri G, Sabet O, Davis L, Chin JW, Bastiaens PI. EGF-dependent re-routing of vesicular recycling switches spontaneous phosphorylation suppression to EGFR signaling. Elife. 2015;4:e12223. pii. [PMC free article] [PubMed]
79. Carmon KS, Loose DS. Wnt7a interaction with Fzd5 and detection of signaling activation using a split eGFP. Biochem Biophys Res Commun. 2008;368:285–291. [PMC free article] [PubMed]
80. Slater PG, Ramirez VT, Gonzalez-Billault C, Varela-Nallar L, Inestrosa NC. Frizzled-5 receptor is involved in neuronal polarity and morphogenesis of hippocampal neurons. PLoS One. 2013;8:e78892. [PMC free article] [PubMed]
81. Asad M, Wong MK, Tan TZ, Choolani M, Low J, Mori S, Virshup D, Thiery JP, Huang RY. FZD7 drives in vitro aggressiveness in Stem-A subtype of ovarian cancer via regulation of non-canonical Wnt/PCP pathway. Cell Death Dis. 2014;5:e1346. [PMC free article] [PubMed]
82. Fernandez A, Huggins IJ, Perna L, Brafman D, Lu D, Yao S, Gaasterland T, Carson DA, Willert K. The WNT receptor FZD7 is required for maintenance of the pluripotent state in human embryonic stem cells. Proc Natl Acad Sci USA. 2014;111:1409–1414. [PubMed]
83. Zheng H, Ying H, Wiedemeyer R, Yan H, Quayle SN, Ivanova EV, Paik JH, Zhang H, Xiao Y, Perry SR, Hu J, Vinjamoori A, Gan B, et al. PLAGL2 regulates Wnt signaling to impede differentiation in neural stem cells and gliomas. Cancer Cell. 2010;17:497–509. [PMC free article] [PubMed]
84. Birnie R, Bryce SD, Roome C, Dussupt V, Droop A, Lang SH, Berry PA, Hyde CF, Lewis JL, Stower MJ, Maitland NJ, Collins AT. Gene expression profiling of human prostate cancer stem cells reveals a pro-inflammatory phenotype and the importance of extracellular matrix interactions. Genome Biol. 2008;9:R83. [PMC free article] [PubMed]
85. Sferruzzi-Perri AN, Macpherson AM, Roberts CT, Robertson SA. Csf2 null mutation alters placental gene expression and trophoblast glycogen cell and giant cell abundance in mice. Biol Reprod. 2009;81:207–221. [PMC free article] [PubMed]
86. Kepler TB, Sample C, Hudak K, Roach J, Haines A, Walsh A, Ramsburg EA. Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler. BMC Genomics. 2010;11:444. [PMC free article] [PubMed]
87. Engelmann JC, Amann T, Ott-Rötzer B, Nützel M, Reinders Y, Reinders J, Thasler WE, Kristl T, Teufel A, Huber CG, Oefner PJ, Spang R, Hellerbrand C. Causal Modeling of Cancer-Stromal Communication Identifies PAPPA as a Novel Stroma-Secreted Factor Activating NFκB Signaling in Hepatocellular Carcinoma. PLoS Comput Biol. 2015;11:e1004293. [PMC free article] [PubMed]
88. Lissat A, Joerschke M, Shinde DA, Braunschweig T, Meier A, Makowska A, Bortnick R, Henneke P, Herget G, Gorr TA, Kontny U. IL6 secreted by Ewing sarcoma tumor microenvironment confers anti-apoptotic and cell-disseminating paracrine responses in Ewing sarcoma cells. BMC Cancer. 2015;15:552. [PMC free article] [PubMed]
89. Sisler JD, Morgan M, Raje V, Grande RC, Derecka M, Meier J, Cantwell M, Szczepanek K, Korzun WJ, Lesnefsky EJ, Harris TE, Croniger CM, Larner AC. The Signal Transducer and Activator of Transcription 1 (STAT1) Inhibits Mitochondrial Biogenesis in Liver and Fatty Acid Oxidation in Adipocytes. PLoS One. 2015;10:e0144444. [PMC free article] [PubMed]
90. Tijchon E, D van Ingen Schenau, van Opzeeland F, Tirone F, Hoogerbrugge PM, FN Van Leeuwen, Scheijen B. Targeted Deletion of Btg1 and Btg2 Results in Homeotic Transformation of the Axial Skeleton. PLoS One. 2015;10:e0131481. [PMC free article] [PubMed]
91. Voce DJ, Schmitt AM, Uppal A, McNerney ME, Bernal GM, Cahill KE, Wahlstrom JS, Nassiri A, Yu X, Crawley CD, White KP, Onel K, Weichselbaum RR, et al. Nfkb1 is a haploinsufficient DNA damage-specific tumor suppressor. Oncogene. 2015;34:2807–2813. [PMC free article] [PubMed]
92. Krämer OH, Knauer SK, Greiner G, Jandt E, Reichardt S, Gührs KH, Stauber RH, Böhmer FD, Heinzel T. A phosphorylation-acetylation switch regulates STAT1 signaling. Genes Dev. 2009;23:223–235. [PubMed]
93. Weidle UH, Klostermann S, Eggle D, Krüger A. Interleukin 6/interleukin 6 receptor interaction and its role as a therapeutic target for treatment of cachexia and cancer. Cancer Genomics Proteomics. 2010;7:287–302. [PubMed]
94. Gasson J. Molecular physiology of granulocyte-macrophage colony-stimulating factor. Blood. 1991;77:1131–1145. [PubMed]
95. Pollard JW. Role of colony-stimulating factor-1 in reproduction and development. Mol Reprod Dev. 1997;46:54–60. discussion 60–1. [PubMed]
96. Robertson SA, Sjöblom C, Jasper MJ, Norman RJ, Seamark RF. Granulocyte-macrophage colony-stimulating factor promotes glucose transport and blastomere viability in murine preimplantation embryos. Biol Reprod. 2001;64:1206–1215. [PubMed]
97. Wang Y, Cheng Z, Elalieh HZ, Nakamura E, Nguyen MT, Mackem S, Clemens TL, Bikle DD, Chang W. IGF-1R signaling in chondrocytes modulates growth plate development by interacting with the PTHrP/Ihh pathway. J Bone Miner Res. 2011;26:1437–1446. [PMC free article] [PubMed]
98. Choi SW, Jeong DU, Kim JA, Lee B, Joeng KS, Long F, Kim DW. Indian Hedgehog signalling triggers Nkx3.2 protein degradation during chondrocyte maturation. Biochem J. 2012;443:789–798. [PMC free article] [PubMed]
99. Madhala-Levy D, Williams VC, Hughes SM, Reshef R, Halevy O. Cooperation between Shh and IGF-I in promoting myogenic proliferation and differentiation via the MAPK/ERK and PI3K/Akt pathways requires Smo activity. J Cell Physiol. 2012;227:1455–1464. [PMC free article] [PubMed]
100. Webb AE, Brunet A. FOXO transcription factors: key regulators of cellular quality control. Trends Biochem Sci. 2014;39:159–169. [PMC free article] [PubMed]
101. Dluhosova M, Curik N, Vargova J, Jonasova A, Zikmund T, Stopka T. Epigenetic control of SPI1 gene by CTCF and ISWI ATPase SMARCA5. PLoS One. 2014;9:e87448. [PMC free article] [PubMed]
102. Hegde S, Hankey P, Paulson RF. Self-renewal of leukemia stem cells in Friend virus-induced erythroleukemia requires proviral insertional activation of Spi1 and hedgehog signaling but not mutation of p53. Stem Cells. 2012;30:121–130. [PMC free article] [PubMed]
103. Myatt SS, Lam EW. The emerging roles of forkhead box (Fox) proteins in cancer. Nature reviews Cancer. 2007;7:847–859. [PubMed]
104. Holmfeldt P, Pardieck J, Saulsberry AC, Nandakumar SK, Finkelstein D, Gray JT, Persons DA, McKinney-Freeman S. Nfix is a novel regulator of murine hematopoietic stem and progenitor cell survival. Blood. 2013;122:2987–2996. [PubMed]
105. Grabowska MM, Elliott AD, DeGraff DJ, Anderson PD, Anumanthan G, Yamashita H, Sun Q, Friedman DB, Hachey DL, Yu X, Sheehan JH, Ahn JM, Raj GV, et al. NFI transcription factors interact with FOXA1 to regulate prostate-specific gene expression. Mol Endocrinol. 2014;28:949–964. [PubMed]
106. Melchor L, Saucedo-Cuevas LP, Muñoz-Repeto I, Rodríguez-Pinilla SM, Honrado E, Campoverde A, Palacios J, Nathanson KL, García MJ, Benítez J. Comprehensive characterization of the DNA amplification at 13q34 in human breast cancer reveals TFDP1 and CUL4A as likely candidate target genes. Breast Cancer Res. 2009;11:R86. [PMC free article] [PubMed]
107. Moscovich M, LeDoux MS, Xiao J, Rampon GL, Vemula SR, Rodriguez RL, Foote KD, Okun MS. Dystonia, facial dysmorphism, intellectual disability and breast cancer associated with a chromosome 13q34 duplication and overexpression of TFDP1: case report. BMC Med Genet. 2013;14:70. [PMC free article] [PubMed]
108. Engström PG, Steijger T, Sipos B, Grant GR, Kahles A, Rätsch G, Goldman N, Hubbard TJ, Harrow J, Guigó R, Bertone P, RGASP Consortium Systematic evaluation of spliced alignment programs for RNA-seq data. Nat Methods. 2013;10:1185–1191. [PMC free article] [PubMed]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC