Skip to main content
Download PDF
Download ePub
  • Research article
  • Open access
  • Published:

Gene expression profile analyze the molecular mechanism of CXCR7 regulating papillary thyroid carcinoma growth and metastasis

Journal of Experimental & Clinical Cancer Research volume 34, Article number: 16 (2015) Cite this article

Abstract

Background

To detect genetic expression profile alterations after papillary thyroid carcinoma (PTC) cells transfected with chemokine receptor CXCR7 gene by gene microarray, and gain insights into molecular mechanisms of how CXCR7 regulating PTC growth and metastasis.

Methods

The Human OneArray microarray was used for a complete genome-wide transcript profiling of CXCR7 transfected PTCs (K1-CXCR7 cells), defined as experimental group. Non CXCR7 transfected PTCs (K1 cells) were used as control group. Differential analysis for per gene was performed with a random variance model and t test, p values were adjusted to control the false discovery rate. Gene ontology (GO) on differentially expressed genes to identify the biological processes in modulating the progression of papillary thyroid carcinoma. Pathway analysis was used to evaluate the signaling pathway that differentially expressed genes were involved in. In addition, quantitative real-time polymerase chain reaction (q-PCR) and Western blot were used to verify the top differentially expression genes.

Results

Comparative analysis revealed that the expression level of 1149 genes was changed in response to CXCR7 transfection. After unsupervised hierarchical clustering analysis, 270 differentially expressed genes were filtered, of them 156 genes were up-regulated whereas 114 genes were down-regulated in K1-CXCR7 cells. GO enrichment analysis revealed the differentially expressed genes were mainly involved in biopolymer metabolic process, signal transduction and protein metabolism. Pathway enrichment analysis revealed differentially expressed genes were mainly involved in ECM-receptor interaction, Focal adhesion, MAPK signaling pathway and Cytokine-cytokine receptor interaction pathway. More importantly, the expression level of genes closely associated with tumor growth and metastasis was altered significantly in K1-CXCR7 cells, including up-regulated genes FN1, COL1A1, COL4A1, PDGFRB, LTB, CXCL12, MMP-11, MT1-MMP and down-regulated genes ITGA7, and Notch-1.

Conclusions

Gene expression profiling analysis of papillary thyroid carcinoma can further delineate the mechanistic insights on how CXCR7 regulating papillary thyroid carcinoma growth and metastasis. CXCR7 may regulate growth and metastasis of papillary thyroid carcinoma via the activation of PI3K/AKT pathway and its downstream NF-κB signaling, as well as the down-regulation of Notch signaling.

Background

Thyroid carcinoma is the most common endocrine neoplasm, and its incidence has been explosively rising worldwide over the past few decades. Papillary thyroid carcinoma is the most common pathological type of thyroid carcinoma, accounting for at least 70–80% of thyroid carcinoma. Although PTC has a favorable prognosis, certain cases exhibit aggressive clinical characteristics, such as lymph node metastasis.

Chemokines and their receptors play a critical role in tumorigenesis, progression, and metastasis of tumor [1]. The chemokine receptor CXCR7 mediates cellular adhesion, migration, proliferation, and survival by binding its ligands stromal cell-derived factor-1(SDF-1) and Interferon-inducible T cell α-chemoattractant (I-TAC) [2,3]. In recent years, accumulating evidences had demonstrated that expression of CXCR7 played a critical role in tumor cell proliferation, angiogenesis, invasion, and metastasis [4-7]. In our previous study, we have demonstrated that CXCR7 and SDF-1 were over-expressed in PTC tissue compared with peritumoral nonmalignant tissue and thyroid benign lesion tissue, and the expressions of them were positively associated with lymph node metastasis [8]. In addition, we found that knockdown of CXCR7 in PTC cells suppressed cell proliferation, invasion, induced S phase arrest, and promoted apoptosis [9].

To further evaluate the signaling pathways involved in CXCR7 receptor regulated PTC progression, we used gene microarray to detect the altered gene expression in PTC cells transfected with CXCR7 and tried to gain insights into molecular mechanisms of how CXCR7 regulating PTC growth and metastasis.

Methods

Cell lines and culture conditions

The human papillary thyroid carcinoma cell line K1 was purchased from European Collection of Animal Cell Cultures. Stable human CXCR7 cDNA transfected cell line, K1-CXCR7, was established in our laboratory [9]. Both cells were cultured in Dulbecco’s Modified Eagle’s Medium: Ham’s F12:MCDB105 (Sigma-Aldrich, St. Louis, Missouri) containing 10% fetal calf serum (FCS; Sigma-Aldrich) and 2 mmol/L glutamine (Sigma-Aldrich). This study was approved by the Ethics Committee in the Affiliated Shengjing Hospital of China Medical University.

RNA preparation and microarray analysis

The gene microarray analysis was carried out by Phalanx Biotech Group, which included RNA amplification, labeling of probe, hybridization, and data extraction. Briefly, total RNA was extracted from K1-CXCR7 cells as experimental group (O1,O2,O3) and K1 cells as control group (N1, N2, N3) using Trizol reagent (TaKaRa Bio Inc, Japan) according to the manufacturer’s instructions. RNA quantity and purity were assessed by using NanoDrop ND-1000 to measure OD260/280. RNA integrity was ascertained by using Agilent RNA 6000 Nano assay to determine RNA Integrity Number (RIN) values. Gene expression profiling was conducted with the Human OneArray® V6.1 microarray (OneArray, China Taiwan) containing 31741 human genome probes and 938 experimental control probes [10]. After hybridization, arrays were washed, scanned and then gene expression results were extracted by DNA Microarray Scanner G2565B (Agilent Technologies, United States) according to the manufacturer’s instructions. Raw fluorescence intensity values were normalized and log-transformed using GeneSpring GX 10 software (Agilent Technologies, United States).

Quantitative real-time polymerase chain reaction

Total RNA was extracted using Trizol reagent (TaKaRa Bio Inc, Japan) according to the manufacturer’s instructions. Quantitative real-time polymerase chain reaction (q-PCR) analysis was performed on Lightcycler480 (Roche Applied Science) according to the manufacturer’s protocol. GAPDH was used as internal control to normalize mRNA levels. All the experiments were repeated three times. Primer sequences are listed in Table 1.

Table 1 Primers used in this study
Full size table

Western blotting

Cells were washed twice with ice-cold phosphate-buffered saline (PBS) and extracted according to protein extraction protocols. Protein concentrations were determined by the BCA Protein Assay Kit (Beyotime Biotechnology, China). Total protein samples (80 ug/lane) were electrophoresed on sodium dodecyl sulfate-polyacrylamide gel and then transferred to Polyvinylidene Fluoride membrane. After blocking with 5% non-fat dry milk for 2 h, membranes were incubated with primary antibodies overnight at 4°C. The membranes were incubated for 2 h at room temperature with secondary antibody. Antibodies used in this study included the following: rabbit polyclonal anti FN1 (1:400), anti-COL1A1 (1:200), anti-COL4A1 (1:500), anti-CXCL12 (1:400), anti-PDGFRB (1:200), anti-MMP-11 (1:200), anti-MTI-MMP (1:200), anti-ITGA7 (1:200, all from BOSTER, Wuhan, China), rabbit polyclonal anti-LTB (1:500, from Abcam, Cambridge, MA), and rabbit polyclonal anti-Notch1 (1:1000), GAPDH (1:10000, both from Proteintech Group Inc, Chicago, IL), and goat anti-rabbit IgG (1:2000; ZSGB-BIO, China) as secondary antibody. Binding was detected using the enhanced chemiluminescence reagents (Beyotime Biotechnology, China). The ratio between the integated optical density of interest proteins and GAPDH of the same sample was calculated as the relative content of protein detected. All experiments were performed in three times.

Statistical analysis

Raw fluorescence intensity values were normalized and log-transformed. Fold change were calculated by Rosetta Resolver 7.2 with error model adjusted by Amersham Pairwise Ration Builder for signal comparison of sample. In bioinformatic analysis, differentially expressed genes were subjected to hierarchical clustering analysis, principal component analysis, Pathway analysis and Gene Ontology analysis. T-test was applied to analyze differences of measurement and categorical data. The SPSS 19.0 software (SPSS Inc., Chicago, IL) was used. All data were expressed as mean ± standard deviation (SD), and the p value <0.05 was considered statistically significant.

Results

Number of differentially expressed genes

In order to detect the number of genes affected by CXCR7 transfection, we used gene microarray analysis to compare the fluorescence intensity ratio between the experimental and control group. Logarithm of fluorescence intensity ratio was represented by Fold change, and log2 ratios ≥ 1.0 or log2 ratios ≤ −1.0 means two Fold change (Figure 1). Standard selection criteria to further identify differentially expressed genes are as follows: |log2 ratios| ≥ 1 and P < 0.05 (log2 ratios ≥ 1.0 means up-regulated and log2 ratios ≤ −1.0 represents down-regulated) (Figure 2). Our data indicated CXCR7 transfection up-regulated 529 genes and down-regulated 620 genes.

Figure 1
figure 1

Histograms with fold change between the experimental and control groups (O vs N). Fold change represents logarithm of fluorescence signal intensity ratios for differentially expressed genes. And log2 ratios ≥ 1.0 or log2 ratios ≤ −1.0 means two Fold change. The histogram plot shows fold change distribution of all probes excluding control and flagged probes. |Fold change| ≥1 means genes differentially expressed.

Full size image
Figure 2
figure 2

Volcano plot of distribution of differentially expressed genes between the experimental and control groups. The dotted line in red and green represent the cut-off, a measurement of gene expression fold-change on the X-axis versus a measure of statistical significance [−Log10 (P-value)] on the Y-axis. Differentially expressed genes are established at |Fold change| ≥1 and P-value < 0.05 (Blue dots in Figure 2).

Full size image

Hierarchical clustering analysis showed that CXCR7 over-expression modified gene clusters in PTC cells

The correlation of expression profiles between the experimental and control groups was demonstrated by unsupervised hierarchical clustering analysis tree (Figure 3). A subset of differential genes that showed similar properties was selected for clustering analysis. An intensity filter was used to select genes where the difference between the maximum and minimum intensity values exceeds 1200 among all microarrays. For this microarray project, the number of genes clustered was 270, among which 156 genes were up-regulated and 114 were down-regulated.

Figure 3
figure 3

Hierarchical clustering of differentially expressed genes in the experimental and control groups. A hierarchical clustering tree indicates the gene expression patterns similarity of the 270 genes between the experimental and control groups. The expression levels of the 270 genes are showed by different color lump: red, high (up-regulated); black, medium; and green, low (down-regulated).

Full size image

Pathway analysis showed that genes regulated by CXCR7 over-expression are mainly involved in KEGG signaling pathway

By Pathway analysis, We found that differentially expressed genes were mainly involved in twenty-five KEGG (Kyoto Encyclopedia of Genes and Genomes) classical pathway, especially in ECM-receptor interaction (13 genes), Focal adhesion (18 genes), MAPK signaling pathway (20 genes), Cytokine-cytokine receptor interaction (19 genes), and four BioCarta pathway. Table 2 lists the top 10 significant enrichment pathway terms. And Table 3 lists the top up- and down-regulated genes those were mainly involved in ECM-receptor interaction pathway, Focal adhesion pathway, MAPK signaling pathway and Cytokine-cytokine receptor interaction pathway.

Table 2 The top 10 enrichment pathway terms differentially expressed genes involved in
Full size table
Table 3 The top up- and down-regulated genes involved in the significant enrichment pathway
Full size table

Gene ontology analysis results

Gene ontology (GO) database is organized into three categories describing molecular function (MF), biological process (BP), and cellular component (CC). To analyze function differences represented by differentially expressed genes, we conducted a GO analysis. We found, in MF, differentially expressed genes significantly enriched in KINASE_ACTIVITY (29 genes), ION_BINDING (24 genes) and PHOSPHOTRANSFERASE_ACTIVITY_ALCOHOL_GROUP_AS_ACCEPTOR (30 genes); in BP, the top enrichment terms were SIGNAL_TRANSDUCTION (89 genes) and energy metabolism process, such as BIOPOLYMER_METABOLIC_PROCESS (83 genes), CELLULAR_MACROMOLECULE_METABOLIC_PROCESS (68 genes) and PROTEIN_METABOLIC_PROCESS (76 genes); in CC, the genes involved in the CYTOPLASM (113 genes) and cell MEMBRANE_PART (90 genes) expression altered significantly. The top 10 of enrichment terms were showed in Table 4.

Table 4 The top 10 enrichment gene ontology terms
Full size table

Verification of differentially expression genes

To partly confirm the expression profiles alteration after PTC cells transfected with CXCR7 gene, we selected up-regulated genes of FN1, COL1A1, COL4A1, PDGFRB, CXCL12, LTB, MMP-11, MT-MMP and down-regulated genes of ITGA7 and Notch-1 for further verification. RT-PCR and Western blot analysis revealed that both mRNA and protein expression levels of genes (FN1, COL1A1,COL4A1, PDGFRB, CXCL12, LTB, MMP-11, MT-MMP) were markedly elevated in K1-CXCR7 cells compared with K1 cells, whereas genes ITGA7 and Notch-1 were the opposite (Figure 4 and Table 5).

Figure 4
figure 4

mRNA and protein expression of differentially expression genes in K1 cell and K1-CXCR7 cell. (A) Relative mRNA expression level of differentially expression genes (FN1, COL1A1, COL4A1, PDGFRB, LTB, CXCL12, MMP-11, MT1-MMP, ITGA7 and Notch-1). (B) Upper panel: protein electrophoregram of differentially expression genes. GAPDH protein was used as a loading control. Lower panel: relative protein levels of differentially expression genes. *p < 0.05.

Full size image
Table 5 Expression of mRNA and protein of genes in K1 cells and K1-CXCR7 cells
Full size table

Discussion

In our study, we filtered 270 differentially expressed genes after PTC cells transfected with CXCR7 cDNA by gene microarray. Among them, 156 genes were up-regulated and 144 genes were down-regulated. The best-known effect of thyroid is an increase in basal energy expenditure achieved by acting on protein, carbohydrate and lipid metabolism. As expected, GO enrichment analysis found several high-enrichment terms linked to metabolism, including biopolymer metabolic process (83 genes), cellular macromolecule metabolic process (68 genes), cellular protein metabolic process (76 genes) and protein metabolic process (76 genes). In addition, differentially expressed genes were significantly involved in signal transduction biological process (89 genes). Pathway enrichment analysis found that the most of differentially expressed genes were enriched in ECM-receptor interaction KEGG pathway (13 genes), focal adhesion KEGG pathway (18 genes), MAPK signaling KEGG pathway (20 genes) and Cytokine-cytokine receptor interaction KEGG pathway (19 genes).

Extracellular matrix (ECM) plays an important role in the development and maintenance of tissue and organ architecture and homeostasis. The synthesis and degradation of ECM components (such as type I collagen, type IV collagen and fibronectin) results in its remodeling and promotes the activated endothelial cells (ECs) proliferation, migration and adhesion to ECM, which contributes to angiogenesis [11].

In this study, fibronectin-1 (FN1) was significantly up-regulated (log2 (Ratio) =5.5288). It was involved in ECM-receptor interaction pathway, focal adhesion pathway, pathways in cancer and regulation of actin cytoskeleton pathway. As an important ECM component, FN1 regulates ECs survival, proliferation, adhesion, migration and angiogenesis [12]. Several researches reported that FN1 was over-expressed in PTC, and it may be a useful biomarker to diagnose PTC [13-15]. ECs facilitate tumor angiogenic process through the activation of focal adhesion kinase (FAK) and downstream PI3K/Akt signal pathway, as well as the activation of NF-κB [12,16]. So, we considered that CXCR7 might promote EC adhesion to ECM by up-regulating the expression of FN1, inducing FAK-mediated activation of PI3K/Akt as well as NF-κB pathway, thereby regulating PTC progression.

COL1A1 (collagen, type I, alpha 1) was also obviously up-regulated (log2 (Ratio) =4.3820), involved in ECM-receptor interaction pathway and focal adhesion pathway in our study, as well as COL4A1 (collagen, type IV, alpha 1). As is well-known, COL1A1 is expressed in many tumor cells and tumor-associated stromal cells. Several researches have demonstrated that COL1A1 plays an important role in angiogenesis and desmoplasia, and the over-expression of COL1A1 was associated with invasion process in PTC [17,18]. What’s more, it was reported that PI3K/AKT and NF-κB pathway were involved in regulating the expression of COL1A1 [19,20]. And COL4A1, as an essential component of ECM, plays an important role in angiogenesis and tumor progression [21]. These suggested that PI3K/AKT and NF-κB pathway may mediate the induction of collagen (COL1A1 and COL4A1) by CXCR7 in PTC cells.

Chemokines and their receptors are best known for their ability to mediate the direct migration of leukocytes in the immune system. As we all known, a close connection exists between inflammation and cancer. Recently, accumulating evidence suggests that alterations in immune function may play a crucial role in PTC initiation. And a history of autoimmune disease, always resulting in tissue destruction and inflammation, has also been associated with the increased risk of PTC [22-24]. As shown in our study, several up-regulated differentially expressed genes were significantly involved in Cytokine-cytokine receptor interaction pathway, including PDGFRB (log2 (Ratio) =3.7321), LTB (log2 (Ratio) =2.9402), CXCL12 (log2 (Ratio) =2.7037). They are all important pro-inflammatory molecules, which play important roles in regulating tumor cells proliferation, invasion, metastasis, angiogenesis and apoptosis. Platelet-derived growth factors (PDGFs), a family of peptides of growth factors, bind to their receptors (PDGFR-α and -β) and stimulate growth and diffusion of cancer cells. It was reported that PDGFR-α and -β were expressed in PTC but not found in normal thyroid tissues [25]. And PDGFR-β can induce the transcription and secretion of vascular endothelial growth factor (VEGF), which plays a critical role in tumor growth, angiogenesis and metastasis [26]. LTB, lymphotoxin beta (TNF superfamily, member 3), is a multifunctional pro-inflammatory cytokine that belongs to the tumor necrosis factor (TNF) superfamily. As we all known, TNF is involved in the regulation of multiple biological processes, including cell proliferation, differentiation and apoptosis. It was reported that an autocrine inflammatory cytokine network of inflammatory cytokine TNF-α, angiogenic factor VEGF, and chemokine CXCL12 existed in ovarian cancer microenvironment and stimulated tumor neoangiogenesis [27]. What’s more, TNF-αinduced the production of VEGF and CXCL12, and VEGF also induced CXCL12 [28,29]. Meanwhile, CXCL12 and VEGF were synergized in facilitating angiogenesis of ovarian cancer [30]. In addition, anti-apoptotic and pro-inflammatory abilities of TNF-αwere resulted from the induction of NF-κB [31,32]. Indeed, NF-κB is one of the main TNF-α downstream effectors, and functions as a transcription factor to regulate genes associated with inflammatory, anti-apoptotic and cell proliferation. Madge LA et al. [33] also demonstrated that lymphotoxin-beta receptor (LTβR) ligands LIGHT and LTα1β2 activate NF-κB pathway, and then up-regulated chemokine CXCL12. CXCL12, also named stromal cell-derived factor 1 (SDF-1), functions as a chemotactic factor for many cells, including T cells, pre-B-cells, monocytes, dendritic cells, and hematopoietic progenitor cells [34,35]. It was also involved in tumor cells migration, invasion and metastasis [36,37]. What’s important, we have demonstrated a significantly higher expression level of SDF-1 in PTC tissue compared with peritumoral nonmalignant tissue and thyroid benign lesion tissue, and the expression of SDF-1 was closely associated with lymph node metastasis of PTC [8]. In addition, our research results suggested that SDF-1 binding to its receptor CXCR7 regulated the directed migration and invasion of PTC cells [9]. So, it is likely that CXCR7 interaction with its ligand CXCL12, induced expression alterations of pro-inflammation genes through the induction of NF-κB signaling, therefore regulated anti-apoptotic and proliferation of PTC cells.

Matrix metalloproteinases (MMPs) mediate tumor cell invasion and metastasis by degrading of ECM. Maeta H et al. [38] demonstrated that the expression of MMPs (MMP-2 and MMP-9), as well as their inhibitors (TIMP-1 and TIMP-2) were enhanced in aggressive PTC. Wani N et al. [39] reported that CXCR7 promoted the metastasis of breast cancer by up-regulating the activity of MMPs. In our study, MMP-11 and membrane type 1 matrix metalloproteinase (MT1-MMP) were both up-regulated (log2 (Ratio) =1.5055, log2 (Ratio) =1.0861), which demonstrated that CXCR7 might promote the secretion of MMP-11 and MT1-MMP by up-regulating genes expression at transcriptional level. Although MMP-11 can’t degrade any ECM component, it is also associated with tumor progression and poor prognosis. MMP-11 was found negative in normal thyroid tissue and thyroid follicular cells but significantly expressed in PTC tissues and thyroid carcinoma cell lines [18,40]. Focal degradation of ECM is the key step in the invasion of cancer cells, MT1-MMP degrades ECM by activating proMMP-2, which promote cancer cells invasion [41]. Nakamura H et al. [42] reported that MT1-MMP expression was correlated with the activation of proMMP-2 and lymph node metastasis of PTC, which suggested expression of proMMP-2 and MT1-MMP-mediated activation played a vital role in the lymph node metastasis of PTC. These results indicated that chemokine receptor CXCR7-induced the transcriptional activation of MMPs promoted the secretion of MMPs proteins and resulted in the degradation of ECM, which induced PTC invasion and lymph node metastasis.

In addition, ITGA7 (integrin, alpha 7), as a cell adhesion molecule, contributes to the interaction between ECM and cells by binding to its ligand integrin β and is involved in multiple biological processes, such as human tissue development, tissue differentiation, and immune responses. Recently, increasing evidences have demonstrated ITGA7 as a tumor suppressor gene in many human malignant neoplasms [43-45], including prostate cancer, liver cancer, glioblastoma multiforme and leiomyosarcoma. Han YC et al. [44] reported that integrin-link kinase interaction with miniature chromosome maintenance 7 (MCM7) and MCM7 phosphorylation may be a critical event in ITGA7 signaling pathway leading to tumor suppression. Tan LZ et al. [45] reported a possible mechanisms of ITGA7-mediated tumor cell growth suppression is that IGTA7 interacts with tissue inhibitor of metalloproteinase 3 (TIMP3), results in the relocation of NF-κB from nucleus to cytoplasm, and down-regulates cyclin D1. These processes led to a remarkable suppression of cell growth. In our study, IGTA7 gene was also significantly down-regulated (log2 (Ratio) = −1.99689) in K1-CXCR7 cell, which suggested that CXCR7 might promote the progression of PTC cell by inhibiting the expression of IGTA7, but functions of ITGA7 in PTC need to be further studied.

Notch1, as a multifunction transmembrane receptor, plays a key role in metazoan development. It also participates in the maintenance of tissue homeostasis by regulating cell proliferation, differentiation and apoptosis [46,47]. Notch1 functions as either a tumor suppressor gene or an oncogene in many human carcinomas, which is cell type-specific. In thyroid cancers, the role of Notch1 signaling is tumor histological differentiation dependent. Notch1 signaling was significantly down-regulated in human anaplastic thyroid carcinoma compared with normal thyroid cells, over-expression of Notch1 reduced cancer cells growth and restored differentiation [48]. It was reported that Notch1 acted as a tumor suppressor in medullary thyroid carcinoma (MTC), and activation of Notch1 inhibited growth of MTC cells and induced apoptosis of MTC cells [49-51]. However, Notch1 as an oncogene or a tumor suppressor gene in PTC remains controversial. Zhang J et al. [52] showed that Notch1 expression was higher in PTC compared with normal thyroid tissues, and higher expression levels of Notch1 was closely associated with lymphatic metastasis and poor prognosis of PTC [53,54]. However, Xiao X et al. [55] showed that Notch1 expression was minimal in papillary and follicular thyroid cancer cells, and activation of Notch1 inhibited growth and proliferation of thyroid cancer cells. Our data indicated that Notch-1 gene was significantly down-regulated (log2 (Ratio) = −1.3895) in K1-CXCR7 cell, which suggested Notch-1 gene might have the effect of tumor suppression in PTC. Several studies have shown that Notch1 signaling might be involved in regulating tumor progression through interacting with multiple signaling pathways, including PI3K/AKT, NF-κB and Wnt pathway [56-58], but further studies were needed to identify the functions and mechanisms of Notch1 signaling in PTC.

Other than the genes mentioned above, there are some genes found to be regulated by CXCR7 in our study but have few reports on, such as deltex homolog (DTX3), matrix-remodeling associated 7 (MXRA7). In our study, expressions of DTX3 and MXRA7 were dramatically increased at transcriptional level (log2 (Ratio) =3.2667, log2 (Ratio) =2.3092). It was reported that DTX3 may function as an ubiquitin ligase protein, which regulate the Notch pathway via some ubiquitin ligase activity [59]. Veiga-Castelli LC et al. [60] demonstrated MXRA7 was over-expressed in ectopic endometrium, and it may result in tissue remodeling according to Gene Ontology analysis. However, the exact functions of these genes are not clear. Further studies on these genes are needed to demonstrate their function and correlation with PTC progression.

Conclusions

In summary, here we demonstrate that CXCR7 mediates the transcriptional expression of multiple signaling molecules, including FN1, COL1A1, COL4A1, PDGFRB, MMP-11, MT1-MMP, LTB, CXCL12, ITGA7 and Notch-1. These signal molecules are involved in PI3K/AKT, NF-κB, Notch signal pathway, which may be associated with papillary thyroid carcinoma growth and metastasis. Chemokine receptor CXCR7 may promote PTC growth and metastasis via the activation of PI3K/AKT pathway and its downstream NF-κB signaling, as well as the down-regulation of Notch1 signaling. In addition, several genes with unknown function were found by gene microarray, such as DTX3 and MXRA7, which need to be further studied.

References

  1. Vandercappellen J, Van Damme J, Struyf S. The role of CXC chemokines and their receptors in cancer. Cancer Lett. 2008;267:226–44.

    Article  CAS  PubMed  Google Scholar 

  2. Boldajipour B, Mahabaleshwar H, Kardash E, Reichman-Fried M, Blaser H, Minina S, et al. Control of chemokine-guided cell migration by ligand sequestration. Cell. 2008;132:463–73.

    Article  CAS  PubMed  Google Scholar 

  3. Burns JM, Summers BC, Wang Y, Melikian A, Berahovich R, Miao Z, et al. A novel chemokine receptor for SDF-1 and I-TAC involved in cell survival, cell adhesion, and tumor development. J Exp Med. 2006;203:2201–13.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Sun X, Cheng G, Hao M, Zheng J, Zhou X, Zhang J, et al. CXCL12/CXCR4/CXCR7 chemokine axis and cancer progression. Cancer Metastasis Rev. 2010;29:709–22.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Hao M, Zheng J, Hou K, Wang J, Chen X, Lu X, et al. Role of chemokine receptor CXCR7 in bladder cancer progression. Biochem Pharmacol. 2012;84:204–14.

    Article  CAS  PubMed  Google Scholar 

  6. Zheng K, Li HY, Su XL, Wang XY, Tian T, Li F, et al. Chemokine receptor CXCR7 regulates the invasion, angiogenesis and tumor growth of human hepatocellular carcinoma cells. J Exp Clin Cancer Res. 2010;29:31.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Singh RK, Lokeshwar BL. The IL-8-regulated chemokine receptor CXCR7 stimulates EGFR signaling to promote prostate cancer growth. Cancer Res. 2011;71:3268–77.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Liu Z, Sun DX, Teng XY, Xu WX, Meng XP, Wang BS. Expression of stromal cellderived factor 1 and CXCR7 in papillary thyroid carcinoma. Endocr Pathol. 2012;23:247–53.

    Article  CAS  PubMed  Google Scholar 

  9. Liu Z, Yang L, Teng X, Zhang H, Guan H. The involvement of CXCR7 in modulating the progression of papillary thyroid carcinoma. J Surg Res. 2014;191:379–88.

    Article  CAS  PubMed  Google Scholar 

  10. The Human OneArray® Service. [http://www.onearray.com.cn/Services/Genome_Std.php]

  11. Wang D, Anderson JC, Gladson CL. The role of the extracellular matrix in angiogenesis in malignant glioma tumors. Brain Pathol. 2005;15:318–26.

    Article  CAS  PubMed  Google Scholar 

  12. Paik JY, Ko BH, Jung KH, Lee KH. Fibronectin stimulates endothelial cell 18F-FDG uptake through focal adhesion kinase-mediated phosphatidylinositol 3-kinase/Akt signaling. J Nucl Med. 2009;50:618–24.

    Article  CAS  PubMed  Google Scholar 

  13. Huang Y, Prasad M, Lemon WJ, Hampel H, Wright FA, Kornacker K, et al. Gene expression in papillary thyroid carcinoma reveals highly consistent profiles. Proc Natl Acad Sci U S A. 2001;98:15044–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Huang Y, Prasad M, Lemon WJ, Hampel H, Wright FA, Kornacker K, et al. Hepatocyte growth factor receptor, matrix metalloproteinase-11, tissue inhibitor of metalloproteinase-1, and fibronectin are up-regulated in papillary thyroid carcinoma: a cDNA and tissue microarray study. Clin Cancer Res. 2003;9:68–75.

    Google Scholar 

  15. Huang Y, Prasad M, Lemon WJ, Hampel H, Wright FA, Kornacker K, et al. Hashimoto’s thyroiditis with papillary thyroid carcinoma (PTC)-like nuclear alterations express molecular markers of PTC. Histopathology. 2004;45:39–46.

    Article  PubMed  Google Scholar 

  16. Wang YH, Dong YY, Wang WM, Xie XY, Wang ZM, Chen RX, et al. Vascular endothelial cells facilitated HCC invasion and metastasis through the Akt and NF-κB pathways induced by paracrine cytokines. J Exp Clin Cancer Res. 2013;32:51.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. St Croix B, Rago C, Velculescu V, Traverso G, Romans KE, Montgomery E, et al. Genes expressed in human tumor endothelium. Science. 2000;289:1197–202.

    Article  CAS  PubMed  Google Scholar 

  18. Lee KY, Huang SM, Li S, Kim JM. Identification of differentially expressed genes in papillary thyroid cancers. Yonsei Med J. 2009;50:60–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Ricupero DA, Poliks CF, Rishikof DC, Cuttle KA, Kuang PP, Goldstein RH. Phosphatidylinositol 3-kinase-dependent stabilization of alpha1(I) collagen mRNA in human lung fibroblasts. Am J Physiol Cell Physiol. 2001;281:C99–105.

    CAS  PubMed  Google Scholar 

  20. Chetty A, Cao GJ, Nielsen HC. Insulin-like Growth Factor-I signaling mechanisms, type I collagen and alpha smooth muscle actin in human fetal lung fibroblasts. Pediatr Res. 2006;60:389–94.

    Article  CAS  PubMed  Google Scholar 

  21. Kalluri R. Basement membranes: structure, assembly and role in tumour angiogenesis. Nat Rev Cancer. 2003;3:422–33.

    Article  CAS  PubMed  Google Scholar 

  22. Bozec A, Lassalle S, Hofman V, Ilie M, Santini J, Hofman P. The thyroid gland: a crossroad in inflammation-induced carcinoma? An ongoing debate with new therapeutic potential. Curr Med Chem. 2010;17:3449–61.

    Article  CAS  PubMed  Google Scholar 

  23. Guarino V, Castellone MD, Avilla E, Melillo RM. Thyroid cancer and inflammation. Mol Cell Endocrinol. 2010;321:94–102.

    Article  CAS  PubMed  Google Scholar 

  24. Weng MY, Huang YT, Liu MF, Lu TH. Incidence of cancer in a nationwide population cohort of 7852 patients with primary Sjogren’s syndrome in Taiwan. Ann Rheum Dis. 2012;71:524–7.

    Article  PubMed  Google Scholar 

  25. Zhang J, Wang P, Dykstra M, Gelebart P, Williams D, Ingham R, et al. Platelet-derived growth factor receptor-α promotes lymphatic metastases in papillary thyroid cancer. J Pathol. 2012;228:241–50.

    Article  CAS  PubMed  Google Scholar 

  26. Gong L, Chen P, Liu X, Han Y, Zhou Y, Zhang W, et al. Expressions of D2-40, CK19, galectin-3, VEGF and EGFR in papillary thyroid carcinoma. Gland Surg. 2012;1:25–32.

    PubMed Central  PubMed  Google Scholar 

  27. Kulbe H, Thompson R, Wilson JL, Robinson S, Hagemann T, Fatah R, et al. The inflammatory cytokine tumor necrosis factor-alpha generates an autocrine tumor-promoting network in epithelial ovarian cancer cells. Cancer Res. 2007;67:585–92.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Galbán S, Fan J, Martindale JL, Cheadle C, Hoffman B, Woods MP, et al. von Hippel-Lindau protein-mediated repression of tumor necrosis factor alpha translation revealed through use of cDNA arrays. Mol Cell Biol. 2003;23:2316–28.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Kaplan RN, Riba RD, Zacharoulis S, Bramley AH, Vincent L, Costa C, et al. VEGFR1-positive haematopoietic bone marrow progenitors initiate the pre-metastatic niche. Nature. 2005;438:820–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Kryczek I, Lange A, Mottram P, Alvarez X, Cheng P, Hogan M, et al. CXCL12 and vascular endothelial growth factor synergistically induce neoangiogenesis in human ovarian cancers. Cancer Res. 2005;65:465–72.

    CAS  PubMed  Google Scholar 

  31. Balkwill F. Tumor necrosis factor or tumor promoting factor? Cytokine Growth Factor Rev. 2002;13:135–41.

    Article  CAS  PubMed  Google Scholar 

  32. Karin M. Nuclear factor-kappaB in cancer development and progression. Nature. 2006;441:431–6.

    Article  CAS  PubMed  Google Scholar 

  33. Madge LA, Kluger MS, Orange JS, May MJ. Lymphotoxin-alpha 1 beta 2 and LIGHT induce classical and noncanonical NF-kappa B-dependent proinflammatory gene expression in vascular endothelial cells. J Immunol. 2008;180:3467–77.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Luster AD. Chemokines–chemotactic cytokines that mediate inflammation. N Engl J Med. 1998;338:436–45.

    Article  CAS  PubMed  Google Scholar 

  35. Christopherson K, Hromas R. Chemokine regulation of normal and pathologic immune responses. Stem Cells. 2001;19:388–96.

    Article  CAS  PubMed  Google Scholar 

  36. Sutton A, Friand V, Brulé-Donneger S, Chaigneau T, Ziol M, Sainte-Catherine O, et al. Stromal cell-derived factor-1/chemokine (C-X-C motif) ligand 12 stimulates human hepatoma cell growth, migration, and invasion. Mol Cancer Res. 2007;5:21–33.

    Article  CAS  PubMed  Google Scholar 

  37. Zhang S, Qi L, Li M, Zhang D, Xu S, Wang N, et al. Chemokine CXCL12 and its receptor CXCR4 expression are associated with perineural invasion of prostate cancer. J Exp Clin Cancer Res. 2008;27:62.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Maeta H, Ohgi S, Terada T. Protein expression of matrix metalloproteinases 2 and 9 and tissue inhibitors of metalloproteinase 1 and 2 in papillary thyroid carcinomas. Virchows Arch. 2001;438:121–8.

    Article  CAS  PubMed  Google Scholar 

  39. Wani N, Nasser MW, Ahirwar DK, Zhao H, Miao Z, Shilo K, et al. C-X-C motif chemokine 12/C-X-C chemokine receptor type 7 signaling regulates breast cancer growth and metastasis by modulating the tumor microenvironment. Breast Cancer Res. 2014;16:R54.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Baldini E, Toller M, Graziano FM, Russo FP, Pepe M, Biordi L, et al. Expression of matrix metalloproteinases and their specific inhibitors in normal and different human thyroid tumor cell lines. Thyroid. 2004;14:881–8.

    Article  CAS  PubMed  Google Scholar 

  41. Hoshino D, Koshikawa N, Suzuki T, Quaranta V, Weaver AM, Seiki M, et al. Establishment and validation of computational model for MT1-MMP dependent ECM degradation and intervention strategies. PLoS Comput Biol. 2012;8:e1002479.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Nakamura H, Ueno H, Yamashita K, Shimada T, Yamamoto E, Noguchi M, et al. Enhanced production and activation of progelatinase A mediated by membrane-type 1 matrix metalloproteinase in human papillary thyroid carcinomas. Cancer Res. 1999;59:467–73.

    CAS  PubMed  Google Scholar 

  43. Ren B, Yu YP, Tseng GC, Wu C, Chen K, Rao UN, et al. Analysis of integrin alpha7 mutations in prostate cancer, liver cancer, glioblastoma multiforme, and leiomyosarcoma. J Natl Cancer Inst. 2007;99:868–80.

    Article  CAS  PubMed  Google Scholar 

  44. Han YC, Yu YP, Nelson J, Wu C, Wang H, Michalopoulos GK, et al. Interaction of integrin-linked kinase and miniature chromosome maintenance 7-mediating integrin {alpha}7 induced cell growth suppression. Cancer Res. 2010;70:4375–84.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Tan LZ, Song Y, Nelson J, Yu YP, Luo JH. Integrin α7 binds tissue inhibitor of metalloproteinase 3 to suppress growth of prostate cancer cells. Am J Pathol. 2013;183:831–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Tschaharganeh DF, Chen X, Latzko P, Malz M, Gaida MM, Felix K, et al. Yes-associated protein up-regulates jagged-1 and activates the NOTCH pathway in human hepatocellular carcinoma. Gastroenterology. 2013;144:1530–42.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Ercan C, Vermeulen JF, Hoefnagel L, Bult P, van der Groep P, van der Wall E, et al. HIF-1α and NOTCH signaling in ductal and lobular carcinomas of the breast. Cell Oncol (Dordr). 2012;35:435–42.

    Article  CAS  Google Scholar 

  48. Ferretti E, Tosi E, Po A, Scipioni A, Morisi R, Espinola MS, et al. Notch signaling is involved in expression of thyrocyte differentiation markers and is down-regulated in thyroid tumors. J Clin Endocrinol Metab. 2008;93:4080–7.

    Article  CAS  PubMed  Google Scholar 

  49. Ning L, Greenblatt DY, Kunnimalaiyaan M, Chen H. Suberoyl bis-hydroxamic acid activates Notch-1 signaling and induces apoptosis in medullary thyroid carcinoma cells. Oncologist. 2008;13:98–104.

    Article  CAS  PubMed  Google Scholar 

  50. Kunnimalaiyaan M, Chen H. Tumor suppressor role of Notch-1 signaling in neuroendocrine tumors. Oncologist. 2007;12:535–42.

    Article  CAS  PubMed  Google Scholar 

  51. Kunnimalaiyaan M, Vaccaro AM, Ndiaye MA, Chen H. Over-expression of the NOTCH1 intracellular domain inhibits cell proliferation and alters the neuroendocrine phenotype of medullary thyroid cancer cells. J Biol Chem. 2006;281:39819–30.

    Article  CAS  PubMed  Google Scholar 

  52. Zhang J, Wang Y, Li D, Jing S. Notch and TGF-β/Smad3 pathways are involved in the interaction between cancer cells and cancer-associated fibroblasts in papillary thyroid carcinoma. Tumour Biol. 2014;35:379–85.

    Article  CAS  PubMed  Google Scholar 

  53. Geers C, Colin IM, Gérard AC. Delta-like 4/Notch pathway is differentially regulated in benign and malignant thyroid tissues. Thyroid. 2011;21:1323–30.

    Article  CAS  PubMed  Google Scholar 

  54. Park HS, Jung CK, Lee SH, Chae BJ, Lim DJ, Park WC, et al. Notch1 receptor as a marker of lymph node metastases in papillary thyroid cancer. Cancer Sci. 2012;103:305–9.

    Article  CAS  PubMed  Google Scholar 

  55. Xiao X, Ning L, Chen H. Notch1 mediates growth suppression of papillary and follicular thyroid cancer cells by histone deacetylase inhibitors. Mol Cancer Ther. 2009;8:350–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Song J, Park S, Kim M, Shin I. Down-regulation of Notch-dependent transcription by Akt in vitro. FEBS Lett. 2008;582:1693–9.

    Article  CAS  PubMed  Google Scholar 

  57. Shahi P, Seethammagari MR, Valdez JM, Xin L, Spencer DM. Wnt and Notch pathways have interrelated opposing roles on prostate progenitor cell proliferation and differentiation. Stem Cells. 2011;29:678–88.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Su C, Chen Z, Luo H, Su Y, Liu W, Cai L, et al. Different patterns of NF-κB and Notch1 signaling contribute to tumor-inducedlymphangiogenesis of esophageal squamous cell carcinoma. J Exp Clin Cancer Res. 2011;22(30):85.

    Article  Google Scholar 

  59. Takeyama K, Aguiar RC, Gu L, He C, Freeman GJ, Kutok JL, et al. The BAL-binding protein BBAP and related Deltex family members exhibit ubiquitin-protein isopeptide ligase activity. J Biol Chem. 2003;278:21930–7.

    Article  CAS  PubMed  Google Scholar 

  60. Veiga-Castelli LC, Silva JC, Meola J, Ferriani RA, Yoshimoto M, Santos SA, et al. Genomic alterations detected by comparative genomic hybridization in ovarian endometriomas. Braz J Med Biol Res. 2010;43:799–805.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This project was supported by the National Natural Science Foundation of China (NO.81072182) and Science Foundation of Liaoning Province of China (NO.2013021100).

Author information

Author notes

    Authors and Affiliations

    1. Department General Surgery, Affiliated Shenjing Hospital, China Medical University, No.36 Sanhao Street, Shenyang, 110004, China

      Hengwei Zhang, Xuyong Teng, Zhangyi Liu, Lei Zhang & Zhen Liu

    Authors
    1. Hengwei Zhang
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Xuyong Teng
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Zhangyi Liu
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Lei Zhang
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Zhen Liu
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding author

    Correspondence to Zhen Liu.

    Additional information

    Competing interests

    The authors declare that they have no competing interests.

    Authors’ contributions

    Zhen L designed and supervised the experiments. HZ performed the experiments, and drafted the manuscript. XT, LZ and Zhangyi L participated in the experiments, and performed the statistical analysis. All authors read and approved the final manuscript.

    Hengwei Zhang, Xuyong Teng, Zhangyi Liu and Lei Zhang contributed equally to this work.

    Rights and permissions

    Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

    The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

    To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

    The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Zhang, H., Teng, X., Liu, Z. et al. Gene expression profile analyze the molecular mechanism of CXCR7 regulating papillary thyroid carcinoma growth and metastasis. J Exp Clin Cancer Res 34, 16 (2015). https://doi.org/10.1186/s13046-015-0132-y

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s13046-015-0132-y

    Keywords

    Journal of Experimental & Clinical Cancer Research

    ISSN: 1756-9966

    Contact us