Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

SVM-T-RFE: a novel gene selection algorithm for identifying metastasis-related genes in colorectal cancer using gene expression profiles

Although metastasis is the principal cause of death cause for colorectal cancer (CRC) patients, the molecular mechanisms underlying CRC metastasis are still not fully understood. In an attempt to identify metastasis-related genes in CRC, we obtained gene expression profiles of 55 early stage primary CRCs, 56 late stage primary CRCs, and 34 metastatic CRCs from the expression project in Oncology (http://www.intgen.org/expo/). We developed a novel gene selection algorithm (SVM-T-RFE), which extends support vector machine recursive feature elimination (SVM-RFE) algorithm by incorporating T-statistic. We achieved highest classification accuracy (100%) with smaller gene subsets (10 and 6, respectively), when classifying between early and late stage primary CRCs, as well as between metastatic CRCs and late stage primary CRCs. We also compared the performance of SVM-T-RFE and SVM-RFE gene selection algorithms on another large-scale CRC dataset and the five public microarray datasets. SVM-T-RFE bestowed SVM-RFE algorithm in identifying more differentially expressed genes, and achieving highest prediction accuracy using equal or smaller number of selected genes. A fraction of selected genes have been reported to be associated with CRC development or metastasis.     

Comprehensive gene and microRNA expression profiling reveals the crucial role of hsa-let-7i and its target genes in colorectal cancer metastasis

Accumulating evidence has demonstrated that miRNAs play important roles in the occurrence and development of colorectal cancer (CRC). However, whether miRNAs are associated with the metastasis of CRC remains largely unexplored. The aim of the current study is to profile miRNAs in different CRC metastatic cell lines to identify the biomarkers in CRC metastasis. Gene and miRNA expression profiling was performed to analyze the global expression of mRNAs and miRNAs in the four human CRC cell lines (LoVo, SW480, HT29 and Caco-2) with different potential of metastasis. Expression patterns of mRNAs and miRNAs were altered in different CRC cell lines. By developing an integrated bioinformatics analysis of gene and miRNA expression patterns, hsa-let-7i was identified to show the highest degree in the microRNA-GO-network and microRNA-Gene-network. The expression level of hsa-let-7i was further validated by qRT-PCR in CRC cells. In addition, the targets of hsa-let-7i were predicted by two programs TargetScan and PicTar, and target genes were validated by expression profiling in the most epresentative LoVo and Caco-2 cell lines. Eight genes including TRIM41, SOX13, SLC25A4, SEMA4F, RPUSD2, PLEKHG6, CCND2, and BTBD3 were identified as hsa-let-7i targets. Our data showed the power of comprehensive gene and miRNA expression profiling and the application of bioinformatics tools in the identification of novel biomarkers in CRC metastasis.     

Network cluster analysis of protein–protein interaction network identified biomarker for early onset colorectal cancer

Colorectal cancer (CRC) is a major cause of morbidity and mortality throughout the world. However, the genetic alterations and molecular mechanism of the early onset CRCs are not fully investigated. The present study aimed to characterize early onset CRC by analyzing its gene expression compared with normal controls and to identify network-based biomarkers of early onset CRC. The gene expression profiles of early onset CRC were downloaded from Gene Expression Omnibus and the differentially expressed genes (DEGs) in CRC patients were identified. Then, a protein–protein interaction (PPI) network was constructed and the clusters in PPI were analyzed by ClusterONE. Furthermore, the gene ontology functional analysis and pathway enrichment analysis were conducted to the modules in PPI network. A systems biology approach integrating microarray data and PPI was further applied to construct a PPI network in CRC. Total 631 DEGs were identified from the early onset CRC compared to healthy controls. These genes were found to be involved in several biological processes, including cell communication, cell proliferation, cell shape and apoptosis. Five functional modules which may play important roles in the initiation of early onset CRC were identified from the PPI network. Functional annotation revealed that these five modules were involved in the pathways of signal transduction, carcinogenesis and metastasis. The hub nodes of these five modules, CDC42, TEX11, QKI, CAV1 and FN1, may serve as the biomarkers of early onset CRC and have the potential to be targets for therapeutic intervention. However, further investigations are still needed to confirm our findings.     

Identification of IMPA2 as the hub gene associated with colorectal cancer and liver metastasis by integrated bioinformatics analysis

Background and ObjectivesColorectal cancer (CRC) is one of the most common malignant tumors worldwide with high incidence and mortality rate, while colorectal liver metastasis (CRLM) is one of the major causes of cancer-related deaths. Therefore, the present study aims to identify the hub gene associated with CRC carcinogenesis and liver metastasis, and then explore its diagnostic and prognostic value as well as the potential regulation mechanism.MethodsThe overlapping differential co-expression genes among CRC, CRLM, and normal tissues were explored on the GSE49355 and GSE81582 datasets from the Gene Expression Omnibus (GEO) database by integrated bioinformatics analysis. Then, the hub prognostic genes were selected from the overlapping genes by univariate Cox proportional hazard analysis and online database Gene Expression Profiling Interactive Analysis 2 (GEPIA2). Subsequently, the clinical value of the hub genes was evaluated in the TCGA and GSE39582 cohorts. Finally, the underlying mechanisms of the hub gene regulating CRC carcinogenesis and metastasis were explored by Gene function annotation and DNA methylation analysis.ResultsInositol mono-phosphatase 2 (IMPA2) was identified as the hub gene associated with CRC carcinogenesis and liver metastasis. IMPA2 had an excellent diagnostic efficiency, and its expression was significantly decreased in CRC and liver metastasis samples, being positively correlated with poor prognosis. Moreover, its low expression was associated with AJCC stage III+IV, T4, N1+2, and M1. In addition, our results revealed that the potential mechanisms used by IMPA2 to mediate CRC carcinogenesis and metastasis could be associated with lipid metabolism and epithelial mesenchymal transition (EMT). Finally, IMPA2 expression could be regulated by DNA methylation.ConclusionsIMPA2 was identified and reported for the first time as a hub gene biomarker in the diagnosis and prognosis of CRC, which could regulate CRC carcinogenesis and liver metastasis through the regulation of lipid metabolism, EMT, and DNA methylation.     

Silencing of RHEB inhibits cell proliferation and promotes apoptosis in colorectal cancer cells via inhibition of the mTOR signaling pathway

Colorectal cancer (CRC) is commonly known as one of the most prominent reasons for cancer-related death in China. Ras homolog enriched in brain (RHEB) and the mammalian target activity of rapamycin (mTOR) signaling pathway were found correlated with CRC, but their specific interaction in CRC was still to be investigated. Therefore, we explored whether RHEB gene silencing affected the cell proliferation, differentiation, and apoptosis by directly targeting the mTOR signaling pathway in cells previously harvested from CRC patients. A microarray analysis was subsequently conducted to investigate the relationship between RHEB and mTOR. Eighty-three adjacent normal tissues and CRC tissues were selected. Immunohistochemistry was carried out to detect the positive expression rates of RHEB and Ki-67 in the CRC tissues. Cells were then transfected with different siRNAs to investigate the potential effects RHEB would have on CRC progression. The expressions of RHEB, 4EBP1, ribosomal protein S6 kinase (p70S6K), proliferating cell nuclear antigen (PCNA), B cell lymphoma 2 (bcl-2), and bcl-2-associated X protein (bax) were determined and then the cell cycle, cell proliferation, and apoptotic rate were also measured. We identified RHEB and mTOR as upregulated genes in CRC. Cells treated with RHEB silencing showed a decreased extent of mTOR, p70S6K, 4EBP1 phosphorylation and expression of RHEB, Ki-67, mTOR, p70S6K, 4EBP1, bcl-2, and PCNA as well as decreased activity of cell proliferation and differentiation; although, the expression of bax was evidently higher. Collectively, our data propose the idea that RHEB gene silencing might repress cell proliferation and differentiation while accelerating apoptosis via inactivating the mTOR signaling pathway.     

Matrix metalloproteinase 2 promotes cell growth and invasion in colorectal cancer

Colorectal cancer (CRC) is the third leading cause of cancer-related death in the western world. In this study, we evaluated the expression of matrix metalloproteinase 2 gene (MMP2) in CRC and analyzed its correlation with clinicopathological features. We found that the expression of MMP2 was significantly higher in CRC tissues than in the colorectal tissues. In addition, high levels of MMP2 protein were positively correlated with the status of tumor size, lymph node metastasis, distant metastasis, Dukes' stage, and tumor invasion. Moreover, patients with higher MMP2 levels had markedly shorter overall survivals than those with low MMP2 levels. Multivariate analysis results suggested that the level of MMP2 expression is an independent prognostic indicator for the survival of patients with CRC. Silencing MMP2 expression in CRC cell lines with lentiviral-mediated shRNA markedly suppressed cell proliferation, colony formation, and invasion. Furthermore, we observed that vascular endothelial growth factor (VEGF) and membrane type 1 (MT1)-MMP protein levels were decreased in MMP2-down-regulated colorectal cells. Therefore, our study demonstrated that MMP2 is an important factor related to carcinogenesis and metastasis of CRC, and MMP2 promotes CRC cell growth and invasion by up-regulating VEGF and MT1-MMP expression, which makes this pathway a potential target for cancer treatment.     

Human embryonic stem cells and metastatic colorectal cancer cells shared the common endogenous human microRNA‐26b

The increase in proliferation and the lack of differentiation of cancer cells resemble what occur in the embryonic stem cells during physiological process of embryogenesis. There are also striking similarities in the behaviour between the invasive placental cells and invasive cancer cells. In the present study, microarrays were used to analyse the global expression of microRNAs in a human embryonic stem cell line (i.e. HUES‐17) and four colorectal cancer (CRC) cell lines (i.e. LoVo, SW480, HT29 and Caco‐2) with different metastatic potentialities. Only the expression of miR‐26b was significant decreased in HUES‐17s and LoVo cells, compared with other three cell lines (P < 0.01). The quantitative real‐time PCR analysis confirmed the results of the microarray analysis. Overexpression of miR‐26b expression by miR‐26 mimics transfection and led to the significant suppression of the cell growth and the induction of apoptosis in LoVo cells in vitro, and the inhibition of tumour growth in vivo. Moreover, the potential targets of miR‐26b was predicted by using bioinformatics, and then the predicted target genes were further validated by comparing gene expression profiles between LoVo and NCM460 cell lines. Four genes (TAF12, PTP4A1, CHFR and ALS2CR2) with intersection were found to be the targets of miR‐26b. MetaCore network analysis further showed that the regulatory pathways of miR‐26b were significantly associated with the invasiveness and metastasis of CRC cells. These data suggest that miR‐26b might serve as a novel prognostic factor and a potential therapeutic target for CRC.     

Ezrin expression and cell survival regulation in colorectal cancer

Colorectal cancer (CRC) is the second largest cause of cancer deaths in the United States. A key barrier that prevents better outcomes for this type of cancer as well as other solid tumors is the lack of effective therapies against the metastatic disease. Thus there is an urgent need to fill this gap in cancer therapy. We utilized a 2D-DIGE proteomics approach to identify and characterize proteins that are differentially regulated between primary colon tumor and liver metastatic deposits of the IGF1R-dependent GEO human CRC xenograft, orthotopically implanted in athymic nude mice that may serve as potential therapeutic targets against CRC metastasis. We observed increased expression of ezrin in liver metastasis in comparison to the primary colonic tumor. Increased ezrin expression was further confirmed by western blot and microarray analyses. Ezrin, a cytoskeletal protein belonging to Ezrin–Radixin–Moesin (ERM) family plays important roles in cell motility, invasion and metastasis. However, its exact function in colorectal cancer is not well characterized. Establishment of advanced GEO cell lines with enhanced liver-metastasizing ability showed a significant increase in ezrin expression in liver metastasis. Increased phosphorylation of ezrin at the T567 site (termed here as p-ezrin T567) was observed in liver metastasis. IHC studies of human CRC patient specimens showed an increased expression of p-ezrin T567 in liver metastasis compared to the primary tumors of the same patient. Ezrin modulation by siRNA, inhibitors and T567A/D point mutations significantly downregulated inhibitors of apoptosis (IAP) proteins XIAP and survivin that have been linked to increased aberrant cell survival and metastasis and increased cell death. Inhibition of the IGF1R signaling pathway by humanized recombinant IGF1R monoclonal antibody MK-0646 in athymic mouse subcutaneous xenografts resulted in inhibition of p-ezrin T567 indicating ezrin signaling is downstream of the IGF1R signaling pathway. We identified increased expression of p-ezrin T567 in CRC liver metastasis in both orthotopically implanted GEO tumors as well as human patient specimens. We report for the first time that p-ezrin T567 is downstream of the IGF1R signaling and demonstrate that ezrin regulates cell survival through survivin/XIAP modulation.     

Identification of potential genes/proteins regulated by Tiam1 in colorectal cancer by microarray analysis and proteome analysis

Tiam1 (T-cell lymphoma invasion and metastasis-inducing protein 1), a guanine nucleotide exchange factor that activates Rac, is a colorectal cancer metastasis-related gene. In this study, we aimed to better understand the mechanism underlying Tiam1-mediated metastasis. We applied gene microarray and proteome analysis and compared expression of genes and proteins in a stable Tiam1-silencing colorectal cancer cell line and in a control cell line. Our analysis identified three genes, high-mobility group box1 (HMGB1), annexin IV (ANXA4) and phosphoglycerate mutase 1 (PGAM1) that were associated with Tiam1. Analysis of these proteins, which may be directly or indirectly regulated by Tiam1, may provide insight into the role and mechanism of Tiam1 in colorectal cancer metastasis.     

Systematic analysis of genes involved in oral cancer metastasis to lymph nodes

Oral cancer remains a deadly disease worldwide. Lymph node metastasis and invasion is one of the causes of death from oral cancer. Elucidating the mechanism of oral cancer lymph node metastasis and identifying critical regulatory genes are important for the treatment of this disease. This study aimed to identify differentially expressed genes (gene signature) and pathways that contribute to oral cancer metastasis to lymph nodes. The GSE70604-associated study compared gene profiles in lymph nodes with metastasis of oral cancer to those of normal lymph nodes. The GSE2280-associated study compared gene profiles in primary tumor of oral cancer with lymph node metastasis to those in tumors without lymph node metastasis. There are 28 common differentially expressed genes (DEGs) showing consistent changes in both datasets in overlapping analysis. GO biological process and KEGG pathway analysis of these 28 DEGs identified the gene signature CCND1, JUN and SPP1, which are categorized as key regulatory genes involved in the focal adhesion pathway. Silencing expression of CCND1, JUN and SPP1 in the human oral cancer cell line OECM-1 confirmed that those genes play essential roles in oral cancer cell invasion. Analysis of clinical samples of oral cancer found a strong correlation of these genes with short survival, especially JUN expression associated with metastasis. Our study identified a unique gene signature – CCND1, JUN and SPP1 – which may be involved in oral cancer lymph node metastasis.     

Layered Signaling Regulatory Networks Analysis of Gene Expression Involved in Malignant Tumorigenesis of Non-Resolving Ulcerative Colitis via Integration of Cross-Study Microarray Profiles

Background

Ulcerative colitis (UC) was the most frequently diagnosed inflammatory bowel disease (IBD) and closely linked to colorectal carcinogenesis. By far, the underlying mechanisms associated with the disease are still unclear. With the increasing accumulation of microarray gene expression profiles, it is profitable to gain a systematic perspective based on gene regulatory networks to better elucidate the roles of genes associated with disorders. However, a major challenge for microarray data analysis is the integration of multiple-studies generated by different groups.

Methodology/Principal Findings

In this study, firstly, we modeled a signaling regulatory network associated with colorectal cancer (CRC) initiation via integration of cross-study microarray expression data sets using Empirical Bayes (EB) algorithm. Secondly, a manually curated human cancer signaling map was established via comprehensive retrieval of the publicly available repositories. Finally, the co-differently-expressed genes were manually curated to portray the layered signaling regulatory networks.

Results

Overall, the remodeled signaling regulatory networks were separated into four major layers including extracellular, membrane, cytoplasm and nucleus, which led to the identification of five core biological processes and four signaling pathways associated with colorectal carcinogenesis. As a result, our biological interpretation highlighted the importance of EGF/EGFR signaling pathway, EPO signaling pathway, T cell signal transduction and members of the BCR signaling pathway, which were responsible for the malignant transition of CRC from the benign UC to the aggressive one.

Conclusions

The present study illustrated a standardized normalization approach for cross-study microarray expression data sets. Our model for signaling networks construction was based on the experimentally-supported interaction and microarray co-expression modeling. Pathway-based signaling regulatory networks analysis sketched a directive insight into colorectal carcinogenesis, which was of significant importance to monitor disease progression and improve therapeutic interventions.     

Genome-wide integrative analysis revealed a correlation between lengths of copy number segments and corresponding gene expression profile

Microarray analysis has been applied to comprehensively reveal the abnormalities of DNA copy number (CN) and gene expression in human cancer research during the last decade. These analyses have individually contributed to identify the genes associated with carcinogenesis, progression, metastasis of tumor cells and poor prognosis of cancer patients. However, it is known that the correlation between profiles of CN and gene expression does not highly correlate. Factors which determine the degree of correlation remain largely unexplained. To investigate one such factor, we performed trend analyses between the lengths of CN segments and corresponding gene expression profiles from microarray data in hepatocellular carcinoma (HCC) and colorectal carcinoma (CRC). Significant correlations were observed in CN gain of HCC and CRC (p<0.05). The trend of the CN loss showed a significant correlation in HCC although there was no correlation between the length of CN loss segments and gene expression in CRC. Our findings suggest that the influence of CN on gene expression highly depends on the length of CN region, especially in the case of CN gain. To the best of our knowledge, this is the first study describing the correlation between lengths of CNA segments and expression profiles of corresponding genes.     

Temporal expression and functional analysis of long non‐coding RNAs in colorectal cancer initiation

Long non‐coding RNAs (lncRNAs) have potential applications in clinical diagnosis and targeted cancer therapies. However, the expression profile of lncRNAs in colorectal cancer (CRC) initiation is still unclear. In this study, the expression profiles of lncRNAs and mRNAs were determined by microarray at specific tumour stages in an AOM/DSS‐induced primary colon cancer model. The temporal expression of lncRNAs was analysed by K‐means clustering. Additionally, weighted correlation network analysis (WGCNA) and gene ontology analysis were performed to construct co‐expression networks and establish functions of the identified lncRNAs and mRNAs. Our results suggested that 4307 lncRNAs and 5798 mRNAs are deregulated during CRC initiation. These differential expression genes (DEGs) exhibited a clear correlation with the differential stage of tumour initiation. WGCNA results suggested that a series of hub lncRNAs are involved in regulating cell stemness, colon inflammation, oxidative stress response and cell death at each stage. Among them, lncRNA H19 was up‐regulated in colon tumours and correlated with poor patient prognosis. Collectively, we have been the first to demonstrate the temporal expression and function of lncRNAs in CRC initiation. These results provide novel diagnosis and therapy targets for CRC.     

乳腺癌转移相关基因与功能识别的可重复性

基于基因表达谱识别乳腺癌转移相关差异表达基因及其功能时,由于基因表达在个体间的变异相对较高而样本量相对较少,由不同研究识别的差异表达基因的可重复性较低。本文基于两套乳腺癌转移基因表达谱,评价两组差异表达基因及其所富集的功能的可重复性。结果显示:在两套表达谱中识别的差异表达基因的表达改变方向高度一致并具有显著的表达相关性;基于两组差异表达基因识别的转移相关功能在两套表达谱中高度可重复,主要涉及细胞分裂、细胞周期、DNA复制、染色体分离、磷酸肌醇介导信号转导和DNA损伤刺激应答等。     

The aberrant cancer metabolic gene carbohydrate sulfotransferase 11 promotes non-small cell lung cancer cell metastasis via dysregulation of ceruloplasmin and intracellular iron balance

Aberrant metabolism has been proposed as one of the emerging hallmarks of cancer. However, the interplay between metabolic disorders and cancer metastasis remains to be defined. To explore the sophisticated metabolic processes during metastatic progression, we analyzed differentially expressed metabolic genes during the epithelial-mesenchymal transition (EMT) of lung cancer cells and defined the EMT-associated metabolic gene signature in lung adenocarcinoma patients. We found that the glycosaminoglycan (GAG)-chondroitin sulfate (CS) biosynthesis pathway was upregulated in the mesenchymal state of lung cancer and associated with poor prognosis. Notably, carbohydrate sulfotransferase 11 (CHST11), a crucial CS biosynthetic enzyme, was confirmed as a poor prognosis marker in non-small cell lung cancer (NSCLC) by immunohistochemical analysis. Moreover, forced CHST11 expression promoted invasion and metastasis, which was abolished by depleting the final product of CS biosynthesis by chondroitinase ABC treatment or active-domain negative CHST11. In vivo metastasis mouse models showed that CHST11 increased lung colonies number and sulfated mucosubstance expression. Furthermore, microarray analysis revealed ceruloplasmin (CP), which facilitated iron metabolism, was the downstream effector of CHST11. CP was upregulated by CHST11 through interferon-γ signaling pathway stimulation and related to unfavorable prognosis. Both forced CP expression and long-term iron treatment increased invasion and lung colony formation. Furthermore, we found 3-AP, an iron chelator, hampered the CHST11-induced metastasis. Our findings implicate that the novel CHST11-CP-iron axis enhances EMT and may serve as a new therapeutic target to treat NSCLC patients.     

Colorectal cancer cell intrinsic fibroblast activation protein alpha binds to Enolase1 and activates NF-κB pathway to promote metastasis

Fibroblast activation protein alpha (FAP) is a marker of cancer-associated fibroblast, which is also expressed in cancer epithelial cells. However, the role of FAP in colorectal cancer (CRC) cells remains to be elucidated. Here we investigate the expression pattern of FAP in CRC tissues and cells to prove that FAP is upregulated in CRC cells. Loss- of and gain-of-function assays identified FAP promotes migration and invasion instead of an effect on cell proliferation. Microarray assays are adopted to identify the different expressed genes after FAP knockdown and gene set enrichment analysis (GSEA) is used to exploit the involved signaling pathway. Our works reveal FAP exerts a function dependent on NF-κB signaling pathway and FAP expression is associated with NF-κB signaling pathway in clinical samples. Our work shows FAP is secreted by CRC cells and soluble FAP could promote metastasis. To investigate the mechanism of FAP influencing the NF-κB signaling pathway, LC/MS is performed to identify the proteins interacting with FAP. We find that FAP binds to ENO1 and activates NF-κB signaling pathway dependent on ENO1. Blocking ENO1 could partially reverse the pro-metastatic effect mediated by FAP. We also provide evidences that both FAP and ENO1 are associated with CRC stages, and high levels of FAP and ENO1 predict a poor survival in CRC patients. In summary, our work could provide a novel mechanism of FAP in CRC cells and a potential strategy for treatment of metastatic CRC.Subject terms: Cancer immunotherapy, Non-coding RNAs     

Alteration in gene expression profile and biological behavior in human lung cancer cell line NL9980 by nm23-H1 gene silencing

Lung cancer is the leading cause of cancer death in both men and women. Tumor metastasis is an essential aspect of lung cancer progression. nm23-H1 is a metastasis suppressor gene. The molecular mechanism by which nm23-H1 suppresses the metastasis is still unclear. Here, we compared the gene expression profile of human large cell lung cancer cell line NL9980 by nm23-H1 gene silencing with that of negative control cells to comprehensively investigate nm23-H1-mediated changes in gene expression of NL9980 cells. Microarray assay revealed that expression of 733-known genes (1.9%, 733/38,500) were altered in response to nm23-H1 gene silencing, including 466 upregulated genes and 267 downregulated. real-time PCR assay of the expression changes indicated that 81.82% (45/55) of verified genes were consistent with that observed in microarray assay. The upregulated genes included MMP-1, -2, SNAI2, CXCL1, 2, 3, PAI-2, while the downregulated genes included cystatin B, TIMP-2, E-cadherin, centrin-2, all of which have been associated with tumor metastasis. Furthermore, we confirmed by Western blot that the expression of MMP-1 and -2 were significantly increased while that of cystatin B was dramatically decreased in NL9980-nm23-H1 silencing cells. The NL9980-nm23-H1 silencing cells exhibited significantly more S phase growth and invasive ability. Thus, silencing of nm23-H1 gene caused metastasis-related gene expression changes in lung cancer cells. The knockdown of nm23-H1 expression may change the lung cancer cells to a more invasive phenotype through alteration in the expression of a set of genes.