Acid-sensing receptor GPR4 plays a crucial role in lymphatic cancer metastasis.

Cancer tissues exhibit an acidic microenvironment owing to the accumulation of protons and lactic acid produced by cancer and inflammatory cells. To examine the role of an acidic microenvironment in lymphatic cancer metastasis, gene expression profiling was conducted using human dermal lymphatic endothelial cells (HDLECs) treated with a low pH medium. Microarray and gene set enrichment analysis revealed that acid treatment induced the expression of inflammation-related genes in HDLECs, including genes encoding chemokines and adhesion molecules. Acid treatment-induced chemokines C-X3-C motif chemokine ligand 1 (CX3CL1) and C-X-C motif chemokine ligand 6 (CXCL6) autocrinally promoted the growth and tube formation of HDLECs. The expression of vascular cell adhesion molecule 1 (VCAM-1) increased in HDLECs after acid treatment in a time-dependent manner, which, in turn, enhanced their adhesion to melanoma cells. Among various acid-sensing receptors, HDLECs basally expressed G protein-coupled receptor 4 (GPR4), which was augmented under the acidic microenvironment. The induction of chemokines or VCAM-1 under acidic conditions was attenuated by GPR4 knockdown in HDLECs. In addition, lymph node metastases in a mouse melanoma model were suppressed by administering an anti-VCAM-1 antibody or a GPR4 antagonist. These results suggest that an acidic microenvironment modifies the function of lymphatic endothelial cells via GPR4, thereby promoting lymphatic cancer metastasis. Acid-sensing receptors and their downstream molecules might serve as preventive or therapeutic targets in cancer.

Histone methyltransferase SUV420H1/KMT5B contributes to poor prognosis in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has a high rate of recurrence and poor prognosis, even after curative surgery. Multikinase inhibitors have been applied for HCC patients, but their effect has been restricted. This study aims to clarify the clinical impact of SUV420H1/KMT5B, one of the methyltransferases for histone H4 at lysine 20, and elucidate the novel mechanisms of HCC progression. We retrospectively investigated SUV420H1 expression using HCC clinical tissue samples employing immunohistochemical analysis (n = 350). We then performed loss-of-function analysis of SUV420H1 with cell cycle analysis, migration assay, invasion assay and RNA sequence for Gene Ontology (GO) pathway analysis in vitro, and animal experiments with xenograft mice in vivo. The SUV420H1-high-score group (n = 154) had significantly poorer prognosis for both 5-year overall and 2-year/5-year disease-free survival than the SUV420H1-low-score group (n = 196) (p < 0.001 and p < 0.05, respectively). The SUV420H1-high-score group had pathologically larger tumor size, more tumors, poorer differentiation, and more positive vascular invasion than the SUV420H1-low-score group. Multivariate analysis demonstrated that SUV420H1 high score was the poorest independent factor for overall survival. SUV420H1 knockdown could suppress cell cycle from G1 to S phase and cell invasion. GO pathway analysis showed that SUV420H1 contributed to cell proliferation, cell invasion, and/or metastasis. Overexpression of SUV420H1 clinically contributed to poor prognosis in HCC, and the inhibition of SUV420H1 could repress tumor progression and invasion both in vitro and in vivo; thus, further analyses of SUV420H1 are necessary for the discovery of future molecularly targeted drugs.

Suppression of antitumor cytokine IL­24 by PRG4 and PAI­1 may promote myxoid liposarcoma cell survival.

Suppression of the antitumor cytokine interleukin-24 (IL-24) is critical for the survival of myxoid liposarcoma (MLS) cells. It has been previously demonstrated by the authors that an MLS-specific chimeric oncoprotein, translocated in liposarcoma-CCAAT/enhancer-binding protein homologous protein (TLS-CHOP), supresses IL24 mRNA expression via induction of proteoglycan 4 (PRG4) to sustain MLS cell proliferation. However, IL-24 has also been revealed to be suppressed by the ubiquitin-proteasome system in human ovarian and lung cancer cells. Therefore, the aim of the present study was to elucidate the mechanism of IL-24 suppression in MLS cells. The results revealed that the proteasome inhibitor, MG-132, induced cell death in MLS cells in vitro; this effect was reduced following IL-24 knockdown. This indicated that proteasomal degradation of IL-24 may be an important process for MLS cell survival. In addition, it was also previously revealed by the authors that knockdown of plasminogen activator inhibitor-1 (PAI-1), a TLS-CHOP downstream molecule, suppressed the growth of MLS cells, thus instigating the investigation of the effect of PAI-1 on IL-24 expression in MLS cells. Double knockdown of PAI-1 and IL-24 negated the growth-suppressive effect of PAI-1 single knockdown in MLS cells. Interestingly, PAI-1 single knockdown did not increase the mRNA expression of IL24, but it did increase the protein abundance of IL-24, indicating that PAI-1 suppressed IL-24 expression by promoting its proteasomal degradation. Moreover, treatment of MLS cells with a PAI-1 inhibitor, TM5275, induced IL-24 protein expression and apoptosis. Collectively, the results of the present as well as previous studies indicated that IL-24 expression may be suppressed at the transcriptional level by PRG4 and at the protein level by PAI-1 in MLS cells. Accordingly, PAI-1 may represent an effective therapeutic target for MLS treatment.

An analysis modality for vascular structures combining tissue-clearing technology and topological data analysis.

The blood and lymphatic vasculature networks are not yet fully understood even in mouse because of the inherent limitations of imaging systems and quantification methods. This study aims to evaluate the usefulness of the tissue-clearing technology for visualizing blood and lymphatic vessels in adult mouse. Clear, unobstructed brain/body imaging cocktails and computational analysis (CUBIC) enables us to capture the high-resolution 3D images of organ- or area-specific vascular structures. To evaluate these 3D structural images, signals are first classified from the original captured images by machine learning at pixel base. Then, these classified target signals are subjected to topological data analysis and non-homogeneous Poisson process model to extract geometric features. Consequently, the structural difference of vasculatures is successfully evaluated in mouse disease models. In conclusion, this study demonstrates the utility of CUBIC for analysis of vascular structures and presents its feasibility as an analysis modality in combination with 3D images and mathematical frameworks.

Bone Morphogenetic Protein Signaling in Cancer; Some Topics in the Recent 10 Years.

Bone morphogenetic proteins (BMPs), members of the transforming growth factor-ß (TGF-ß) family, are multifunctional cytokines. BMPs have a broad range of functions, and abnormalities in BMP signaling pathways are involved in cancer progression. BMPs activate the proliferation of certain cancer cells. Malignant phenotypes of cancer cells, such as increased motility, invasiveness, and stemness, are enhanced by BMPs. Simultaneously, BMPs act on various cellular components and regulate angiogenesis in the tumor microenvironment. Thus, BMPs function as pro-tumorigenic factors in various types of cancer. However, similar to TGF-ß, which shows both positive and negative effects on tumorigenesis, BMPs also act as tumor suppressors in other types of cancers. In this article, we review important findings published in the recent decade and summarize the pro-oncogenic functions of BMPs and their underlying mechanisms. The current status of BMP-targeted therapies for cancers is also discussed.

An in vivo orthotopic serial passaging model for a metastatic renal cancer study.

We developed an in vivo serial passaging model for renal cancer with orthotopic renal subcapsular inoculation. We detail the procedures for renal subcapsular inoculation of cancer cells in mice, followed by in vivo and exvivo bioluminescence imaging, tumor-bearing kidney dissociation, and in vivo passaging. This protocol is capable of reproducing the coevolution between cancer cells and the primary tumor microenvironment. It enables us to unveil the systemic dynamics of metastasis and develop a therapeutic strategy for metastatic renal cancer. For complete details on the use and execution of this protocol, please refer to Nishida et al. (2020).

Genome-wide analysis of DNA methylation identifies the apoptosis-related gene UQCRH as a tumor suppressor in renal cancer.

DNA hypermethylation is frequently observed in clear cell renal cell carcinoma (ccRCC) and correlates with poor clinical outcomes. However, the detailed function of DNA hypermethylation in ccRCC has not been fully uncovered. Here, we show the role of DNA methylation in ccRCC progression through the identification of a target(s) of DNA methyltransferases (DNMT). Our preclinical model of ccRCC using the serial orthotopic inoculation model showed the upregulation of DNMT3B in advanced ccRCC. Pretreatment of advanced ccRCC cells with 5-aza-deoxycytidine, a DNMT inhibitor, attenuated the formation of primary tumors through the induction of apoptosis. DNA methylated sites were analyzed genome-wide using methylation array in reference to RNA-sequencing data. The gene encoding ubiquinol cytochrome c reductase hinge protein (UQCRH), one of the components of mitochondrial complex III, was extracted as a methylation target in advanced ccRCC. Immunohistochemical analysis revealed that the expression of UQCRH in human ccRCC tissues was lower than normal adjacent tissues. Silencing of UQCRH attenuated the cytochrome c release in response to apoptotic stimuli and resulted in enhancement of primary tumor formation in vivo, implying the tumor-suppressive role of UQCRH. Moreover, 5-aza-deoxycytidine enhanced the therapeutic efficiency of mammalian target of rapamycin inhibitor everolimus in vivo. These findings suggested that the DNMT3B-induced methylation of UQCRH may contribute to renal cancer progression and implicated clinical significance of DNMT inhibitor as a therapeutic option for ccRCC.

Heterogenous expression of endoglin marks advanced renal cancer with distinct tumor microenvironment fitness.

Intratumoral heterogeneity, including in clear cell renal cell carcinoma, is a potential cause of drug resistance and metastatic cancer progression. We specified the heterogeneous population marked by endoglin (also known as CD105) in a preclinical model of clear cell renal cell carcinoma progression. Highly malignant derivatives of human clear cell renal cell carcinoma OS-RC-2 cells were established as OS5Ks by serial orthotopic inoculation in our previous study. Expression of both ENG (encoding endoglin) mRNA and protein were heterogeneously upregulated in OS5Ks, and the endoglin-positive (ENG+ ) population exhibited growth dependency on endoglin in anchorage-independent cultures. Despite the function of endoglin as a type III receptor, transforming growth factor ß and bone morphogenetic protein-9 signaling were unlikely to contribute to the proliferative phenotype. Although endoglin has been proposed as a marker for renal cancer-initiating cells, the OS5K-3 ENG+ population did not enrich other reported cancer-initiating cell markers or differentiate into the ENG- population. Mouse tumor inoculation models revealed that the tumor-forming capabilities of OS5K-3 ENG+ and ENG- cells in vivo were highly dependent on the microenvironment, with the renal microenvironment most preferable to ENG+ cells. In conclusion, the renal microenvironment, rather than the hypothesized ENG+ cell-centered hierarchy, maintains cellular heterogeneity in clear cell renal cell carcinoma. Therefore, the effect of the microenvironment should be considered when evaluating the proliferative capability of renal cancer cells in the experimental settings.

Visualization of the cancer cell cycle by tissue-clearing technology using the Fucci reporter system.

Tissue-clearing technology is an emerging imaging technique currently utilized not only in neuroscience research but also in cancer research. In our previous reports, tissue-clearing methods were used for the detection of metastatic tumors. Here, we showed that the cell cycles of primary and metastatic tumors were visualized by tissue-clearing methods using a reporter system. First, we established cancer cell lines stably expressing fluorescent ubiquitination-based cell cycle indicator (Fucci) reporter with widely used cancer cell lines A549 and 4T1. Fluorescence patterns of the Fucci reporter were investigated in various tumor inoculation models in mice. Interestingly, fluorescence patterns of the Fucci reporter of tumor colonies were different between various organs, and even among colonies in the same organs. The effects of antitumor drugs were also evaluated using these Fucci reporter cells. Of the three antitumor drugs studied, 5-fluorouracil treatment on 4T1-Fucci cells resulted in characteristic fluorescent patterns by the induction of G2 /M arrest both in vitro and in vivo. Thus, the combination of a tissue-clearing method with the Fucci reporter is useful for analyzing the mechanisms of cancer metastasis and drug resistance.

EHF suppresses cancer progression by inhibiting ETS1-mediated ZEB expression.

ETS homologous factor (EHF) belongs to the epithelium-specific subfamily of the E26 transformation-specific (ETS) transcription factor family. Currently, little is known about EHF's function in cancer. We previously reported that ETS1 induces expression of the ZEB family proteins ZEB1/δEF1 and ZEB2/SIP1, which are key regulators of the epithelial-mesenchymal transition (EMT), by activating the ZEB1 promoters. We have found that EHF gene produces two transcript variants, namely a long form variant that includes exon 1 (EHF-LF) and a short form variant that excludes exon 1 (EHF-SF). Only EHF-SF abrogates ETS1-mediated activation of the ZEB1 promoter by promoting degradation of ETS1 proteins, thereby inhibiting the EMT phenotypes of cancer cells. Most importantly, we identified a novel point mutation within the conserved ETS domain of EHF, and found that EHF mutations abolish its original function while causing the EHF protein to act as a potential dominant negative, thereby enhancing metastasis in vivo. Therefore, we suggest that EHF acts as an anti-EMT factor by inhibiting the expression of ZEBs, and that EHF mutations exacerbate cancer progression.

Whole-organ analysis of TGF-ß-mediated remodelling of the tumour microenvironment by tissue clearing.

Tissue clearing is one of the most powerful strategies for a comprehensive analysis of disease progression. Here, we established an integrated pipeline that combines tissue clearing, 3D imaging, and machine learning and applied to a mouse tumour model of experimental lung metastasis using human lung adenocarcinoma A549 cells. This pipeline provided the spatial information of the tumour microenvironment. We further explored the role of transforming growth factor-ß (TGF-ß) in cancer metastasis. TGF-ß-stimulated cancer cells enhanced metastatic colonization of unstimulated-cancer cells in vivo when both cells were mixed. RNA-sequencing analysis showed that expression of the genes related to coagulation and inflammation were up-regulated in TGF-ß-stimulated cancer cells. Further, whole-organ analysis revealed accumulation of platelets or macrophages with TGF-ß-stimulated cancer cells, suggesting that TGF-ß might promote remodelling of the tumour microenvironment, enhancing the colonization of cancer cells. Hence, our integrated pipeline for 3D profiling will help the understanding of the tumour microenvironment.

Neurotensin receptor 1 signaling promotes pancreatic cancer progression.

Pancreatic cancer is one of the cancers with the poorest prognosis, with a 5-year survival rate of approximately 5-10%. Thus, it is urgent to identify molecular targets for the treatment of pancreatic cancer. Using serial transplantations in a mouse pancreatic orthotopic inoculation model, we previously produced highly malignant pancreatic cancer sublines with increased tumor-forming abilities in vivo. Here, we used these sublines to screen molecular targets for the treatment of pancreatic cancer. Among the genes with increased expression levels in the sublines, we focused on those encoding cell surface receptors that may be involved in the interactions between cancer cells and the tumor microenvironment. Based on our previous RNA-sequence analysis, we found increased expression levels of neurotensin (NTS) receptor 1 (NTSR1) in highly malignant pancreatic cancer sublines. Furthermore, re-analysis of clinical databases revealed that the expression level of NTSR1 was increased in advanced pancreatic cancer and that high NTSR1 levels were correlated with a poor prognosis. Overexpression of NTSR1 in human pancreatic cancer cells Panc-1 and SUIT-2 accelerated their tumorigenic and metastatic abilities in vivo. In addition, RNA-sequence analysis showed that MAPK and NF-κB signaling pathways were activated upon NTS stimulation in highly malignant cancer sublines and also revealed many new target genes for NTS in pancreatic cancer cells. NTS stimulation increased the expression of MMP-9 and other pro-inflammatory cytokines and chemokines in pancreatic cancer cells. Moreover, the treatment with SR48692, a selective NTSR1 antagonist, suppressed the activation of the MAPK and NF-κB signaling pathways and induction of target genes in pancreatic cancer cells in vitro, while the administration of SR48692 attenuated the tumorigenicity of pancreatic cancer cells in vivo. These findings suggest that NTSR1 may be a prognostic marker and a molecular target for pancreatic cancer treatment.

Protocol for Imaging and Analysis of Mouse Tumor Models with CUBIC Tissue Clearing.

Tissue-clearing technologies have developed rapidly in the past decade, especially for use in neuroscience research. We previously reported that CUBIC, which is one tissue-clearing method, is useful for applications in cancer research. CUBIC cancer analysis can be used to detect cancer metastasis with single-cell resolution at whole mouse body/organ level. This approach can also analyze the tumor characteristics with high-quality 3D images. Here, we describe a detailed CUBIC cancer protocol from tissue clearing, capturing 3D images and post-immunohistochemistry. For complete details on the use and execution of this protocol, please refer to Kubota et al. (2017).

Epigenetic remodelling shapes inflammatory renal cancer and neutrophil-dependent metastasis.

Advanced clear cell renal cell carcinoma (ccRCC) frequently causes systemic inflammation. Recent studies have shown that cancer cells reshape the immune landscape by secreting cytokines or chemokines. This phenotype, called cancer-cell-intrinsic inflammation, triggers a metastatic cascade. Here, we identified the functional role and regulatory mechanism of inflammation driven by advanced ccRCC cells. The inflammatory nature of advanced ccRCC was recapitulated in a preclinical model of ccRCC. Amplification of cancer-cell-intrinsic inflammation during ccRCC progression triggered neutrophil-dependent lung metastasis. Massive expression of inflammation-related genes was transcriptionally activated by epigenetic remodelling through mechanisms such as DNA demethylation and super-enhancer formation. A bromodomain and extra-terminal motif inhibitor synchronously suppressed C-X-C-type chemokines in ccRCC cells and decreased neutrophil-dependent lung metastasis. Overall, our findings provide insight into the nature of inflammatory ccRCC, which triggers metastatic cascades, and suggest a potential therapeutic strategy.

Comparative analysis of TTF-1 binding DNA regions in small-cell lung cancer and non-small-cell lung cancer.

Thyroid transcription factor-1 (TTF-1, encoded by the NKX2-1 gene) is highly expressed in small-cell lung carcinoma (SCLC) and lung adenocarcinoma (LADC), but how its functional roles differ between SCLC and LADC remains to be elucidated. Here, we compared the genome-wide distributions of TTF-1 binding regions and the transcriptional programs regulated by TTF-1 between NCI-H209 (H209), a human SCLC cell line, and NCI-H441 (H441), a human LADC cell line, using chromatin immunoprecipitation-sequencing (ChIP-seq) and RNA-sequencing (RNA-seq). TTF-1 binding regions in H209 and H441 cells differed by 75.0% and E-box motifs were highly enriched exclusively in the TTF-1 binding regions of H209 cells. Transcriptome profiling revealed that TTF-1 is involved in neuroendocrine differentiation in H209 cells. We report that TTF-1 and achaete-scute homolog 1 (ASCL1, also known as ASH1, an E-box binding basic helix-loop-helix transcription factor, and a lineage-survival oncogene of SCLC) are coexpressed and bound to adjacent sites on target genes expressed in SCLC, and cooperatively regulate transcription. Furthermore, TTF-1 regulated expression of the Bcl-2 gene family and showed antiapoptotic function in SCLC. Our findings suggest that TTF-1 promotes SCLC growth and contributes to neuroendocrine and antiapoptotic gene expression by partly coordinating with ASCL1.

c-Ski accelerates renal cancer progression by attenuating transforming growth factor ß signaling.

Although transforming growth factor beta (TGF-ß) is known to be involved in the pathogenesis and progression of many cancers, its role in renal cancer has not been fully investigated. In the present study, we examined the role of TGF-ß in clear cell renal carcinoma (ccRCC) progression in vitro and in vivo. First, expression levels of TGF-ß signaling pathway components were examined. Microarray and immunohistochemical analyses showed that the expression of c-Ski, a transcriptional corepressor of Smad-dependent TGF-ß and bone morphogenetic protein (BMP) signaling, was higher in ccRCC tissues than in normal renal tissues. Next, a functional analysis of c-Ski effects was carried out. Bioluminescence imaging of renal orthotopic tumor models demonstrated that overexpression of c-Ski in human ccRCC cells promoted in vivo tumor formation. Enhancement of tumor formation was also reproduced by the introduction of a dominant-negative mutant TGF-ß type II receptor into ccRCC cells. In contrast, introduction of the BMP signaling inhibitor Noggin failed to accelerate tumor formation, suggesting that the tumor-promoting effect of c-Ski depends on the inhibition of TGF-ß signaling rather than of BMP signaling. Finally, the molecular mechanism of the tumor-suppressive role of TGF-ß was assessed. Although TGF-ß signaling did not affect tumor angiogenesis, apoptosis of ccRCC cells was induced by TGF-ß. Taken together, these findings suggest that c-Ski suppresses TGF-ß signaling in ccRCC cells, which, in turn, attenuates the tumor-suppressive effect of TGF-ß.

Efficacy of an orally active small-molecule inhibitor of RANKL in bone metastasis.

Bone is one of the preferred sites for the metastasis of malignant tumours, such as breast cancer, lung cancer and malignant melanoma. Tumour cells colonizing bone have the capacity to induce the expression of receptor activator of nuclear factor-κB ligand (RANKL), which promotes osteoclast differentiation and activation. Tumour-induced osteoclastic bone resorption leads to a vicious cycle between tumours and bone cells that fuels osteolytic tumour growth, causing bone pain and hypercalcaemia. Furthermore, RANKL contributes to bone metastasis by acting as a chemoattractant to bone for tumour cells that express its receptor, RANK. Thus inhibition of the RANKL-RANK pathway is a promising treatment for bone metastasis, and a human monoclonal anti-RANKL antibody, denosumab, has been used in the clinic. However, orally available drugs targeting RANKL must be developed to increase the therapeutic benefits to patients. Here we report the efficacy of the small-molecule RANKL inhibitor AS2676293 in treating bone metastasis using mouse models. Oral administration of AS2676293 markedly inhibited bone metastasis of human breast cancer cells MDA-MB-231-5a-D-Luc2 as well as tumour-induced osteolysis. AS2676293 suppressed RANKL-mediated tumour migration in the transwell assay and inhibited bone metastasis of the murine cell line B16F10, which is known not to trigger osteoclast activation. Based on the results from this study, RANKL inhibition with a small-molecule compound constitutes a promising therapeutic strategy for treating bone metastasis by inhibiting both osteoclastic bone resorption and tumour migration to bone.

Soluble RANKL is physiologically dispensable but accelerates tumour metastasis to bone.

Receptor activator of NF-κB ligand (RANKL) is a multifunctional cytokine known to affect immune and skeletal systems, as well as oncogenesis and metastasis1-4. RANKL is synthesized as a membrane-bound molecule, and cleaved into its soluble form by proteases5-7. As the soluble form of RANKL does not contribute greatly to bone remodelling or ovariectomy-induced bone loss8, whether soluble RANKL has a role in pathological settings remains unclear. Here we show that soluble RANKL promotes the formation of tumour metastases in bone. Mice that selectively lack soluble RANKL (Tnfsf11ΔS/ΔS)5-7,9 have normal bone homoeostasis and develop a normal immune system but display markedly reduced numbers of bone metastases after intracardiac injection of RANK-expressing melanoma and breast cancer cells. Deletion of soluble RANKL does not affect osteoclast numbers in metastatic lesions or tumour metastasis to non-skeletal tissues. Therefore, soluble RANKL is dispensable for physiological regulation of bone and immune systems, but has a distinct and pivotal role in the promotion of bone metastases.

Intracellular and extracellular TGF-ß signaling in cancer: some recent topics.

Transforming growth factor (TGF)-ß regulates a wide variety of cellular responses, including cell growth arrest, apoptosis, cell differentiation, motility, invasion, extracellular matrix production, tissue fibrosis, angiogenesis, and immune function. Although tumor-suppressive roles of TGF-ß have been extensively studied and well-characterized in many cancers, especially at early stages, accumulating evidence has revealed the critical roles of TGF-ß as a pro-tumorigenic factor in various types of cancer. This review will focus on recent findings regarding epithelial-mesenchymal transition (EMT) induced by TGF-ß, in relation to crosstalk with some other signaling pathways, and the roles of TGF-ß in lung and pancreatic cancers, in which TGF-ß has been shown to be involved in cancer progression. Recent findings also strongly suggested that targeting TGF-ß signaling using specific inhibitors may be useful for the treatment of some cancers. TGF-ß plays a pivotal role in the differentiation and function of regulatory T cells (Tregs). TGF-ß is produced as latent high molecular weight complexes, and the latent TGF-ß complex expressed on the surface of Tregs contains glycoprotein A repetitions predominant (GARP, also known as leucine-rich repeat containing 32 or LRRC32). Inhibition of the TGF-ß activities through regulation of the latent TGF-ß complex activation will be discussed.

Pancreatic tumor microenvironment confers highly malignant properties on pancreatic cancer cells.

Tumor microenvironment plays a pivotal role in cancer progression; however, little is known regarding how differences in the microenvironment affect characteristics of cancer cells. Here, we investigated the effects of tumor microenvironment on cancer cells by using mouse tumor models. After three cycles of inoculation and extraction of human pancreatic cancer cells, including SUIT-2 and Panc-1 cells, from tumors, distinct cancer cell lines were established: 3P cells from the pancreas obtained using the orthotopic tumor model and 3sc cells from subcutaneous tissue obtained using the subcutaneous tumor model. On re-inoculation of these cells, the 3sc cells and, more prominently, the 3P cells, exhibited higher tumorigenic activity than the parental cells. The 3P cells specifically exhibited low E-cadherin expression and high invasiveness, suggesting that they were endowed with the highest malignant characteristics. RNA-sequence analysis demonstrated that distinct signaling pathways were activated in each cell line and that the 3P cells acquired a cancer stem cell-like phenotype. Among cancer stem cell-related genes, those specifically expressed in the 3P cells, including NES, may be potential new targets for cancer therapy. The mechanisms underlying the development of highly malignant cancer cell lines were investigated. Individual cell clones within the parental cells varied in tumor-forming ability, indicating the presence of cellular heterogeneity. Moreover, the tumor-forming ability and the gene expression profile of each cell clone were altered after serial orthotopic inoculations. The present study thus suggests that both selection and education processes by tumor microenvironment are involved in the development of highly malignant cancer cells.