Identification of differentially expressed genes between lung adenocarcinoma and lung squamous cell carcinoma by gene expression profiling.

The present study aimed to identify the differentially expressed genes (DEGs) between lung adenocarcinoma and normal lung tissues, and between lung squamous cell carcinoma and normal lung tissues, with the purpose of identifying potential biomarkers for the treatment of lung cancer. The gene expression profile (GSE6044) was downloaded from the Gene Expression Omnibus database, which included data from 10 lung adenocarcinoma samples, 10 lung squamous cell carcinoma samples, and five matched normal lung tissue samples. After data processing, DEGs were identified using the Student's t‑test adjusted via the Benjamini‑Hochberg method. Subsequently, Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of the DEGs was performed using the Database for Annotation, Visualization and Integrated Discovery, and a global network was constructed. A total of 95 upregulated and 241 downregulated DEGs were detected in lung adenocarcinoma samples, and 204 upregulated and 285 downregulated DEGs were detected in lung squamous cell carcinoma samples, as compared with the normal lung tissue samples. The DEGs in the lung squamous cell carcinoma group were enriched in the following three pathways: Hsa04110, Cell cycle; hsa03030, DNA replication; and hsa03430, mismatch repair. However, the DEGs in the lung adenocarcinoma group were not significantly enriched in any specific pathway. Subsequently, a global network of lung cancer was constructed, which consisted of 341 genes and 1,569 edges, of which the top five genes were HSP90AA1, BCL2, CDK2, KIT and HDAC2. The expression trends of the above genes were different in lung adenocarcinoma and lung squamous cell carcinoma when compared with normal tissues. Therefore, these genes were suggested to be crucial genes for differentiating lung adenocarcinoma and lung squamous cell carcinoma.

Introduction

Lung cancer is a highly prevalent type of cancer (1), and is a leading cause of cancer-associated mortality worldwide (2). In 2012, ~11.9% of cancer cases diagnosed were attributed to lung cancer in Europe (3). In China, lung cancer accounts for 18.5% of cancer cases, and has a mortality rate of 23.1%, ranking it the most prevalent and most lethal among all cancer types (4). Lung cancer is classified as follows: Small cell lung cancer (SCLC) and non-SCLC (NSCLC) (5). Adenocarcinoma, squamous cell carcinoma and large cell cancer are the three main subtypes of NSCLC (6). Following surgical treatment, the majority of cases of NSCLC will relapse within 5 years; therefore, it is of great importance to fully elucidate the mechanism underlying the progression of NSCLC.

Microarray screening has been used to identify differentially expressed genes (DEGs) in cancer samples, in order to better understand the mechanisms underlying cancer development (7). Navab et al (8) identified a prognostic gene-expression signature that contained a subset of 11 genes, which were validated in numerous independent NSCLC gene expression databases. In addition, since lung adenocarcinoma and squamous cell carcinoma exhibit different histology, gene expression levels, particularly the levels of relevant markers such as cytokeratin 5/6, differ between them (9). Therefore, specific genes and microRNAs may be used to distinguish between these two types of cancer (10), and the genetic signatures of these cancer types may differ (11). However, the differences in expression between lung adenocarcinoma and squamous cell carcinoma have yet to be fully elucidated.

The present study used gene expression files downloaded from the Gene Expression Omnibus (GEO) database, compared the DEGs detected between lung adenocarcinoma and squamous cell carcinoma samples, and conducted function and pathway enrichment analyses. In addition, a protein-protein interaction (PPI) network of the DEGs was constructed. Genes that exhibited higher degrees in the networks were selected as the key genes in the two lung cancer groups.

Materials and methods

Microarray data

The GSE6044 gene expression profile was downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo/). The profile included data from 10 lung adenocarcinoma samples, 10 lung squamous cell carcinoma samples, and five matched normal lung tissue samples. The microarray data were based on the GPL201 [HG-Focus] platform (Affymetrix Human HG-Focus Target Array; Affymetrix, Inc., Santa Clara, CA, USA).

Data processing

The raw CEL data were rectified and standardized by Robust Multichip Average using the Affy package in R language (bioconductor.org/packages/release/bioc/html/affy.html), in order to obtain the corresponding expression data of the probes. Subsequently, redundant probes that did not correspond with an Entrez Gene ID were deleted, whereas the median value was used for probes that corresponded with several Entrez Gene IDs.

DEGs screening

Student's t-test was used to screen for DEGs in the two types of lung cancer samples. The P-values were further adjusted using the false discovery rate (FDR) approach, according to the Benjamini-Hochberg (BH) method (12). Genes with a FDR <0.1 and fold change (FC) >1.5 or <0.67 were considered DEGs between the cancer and normal samples.

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment

In order to identify the pathways associated with the two cancer groups, tools in the Database for Annotation, Visualization and Integrated Discovery (DAVID; david.abcc.ncifcrf.gov) (13,14) were used to screen the pathways enriched for in the DEGs from the cancer samples using the Expression Analysis Systematic Explorer score (a modified Fisher's exact t-test) with BH multiple testing correction (12). A KEGG pathway with a BH-corrected P<0.05 was considered to be significantly enriched.

Global network construction of cancer-associated genes

A total of 14 cancer-associated pathways were downloaded from the KEGG database (www.genome.jp/kegg), including colorectal cancer, pancreatic cancer, glioma, thyroid cancer, acute myeloid leukemia, chronic myeloid leukemia, basal cell carcinoma, melanoma, renal cell carcinoma, bladder cancer, prostate cancer, endometrial cancer, SCLC and NSCLC. Annotated genes from the 14 pathways were used to construct the global network according to their interaction in pathways by the Cytoscape software (www.cytoscape.org). The expression (upregulated or downregulated) of these genes were further annotated by mapping to the present gene analysis results in two types of lung cancer samples in comparison with normal samples. Subsequently, the enriched DEGs were screened to construct the PPI network in which the interaction relationships between proteins were downloaded from the Human Protein Reference Database (HPRD; www.hprd.org) database (15).

Results

DEGs selection

A total of 8,348 genes were obtained from the GSE6044 gene expression profile. Under the screening criteria of FDR<0.1 and FC >1.5 or P<0.67, 95 upregulated DEGs and 241 downregulated DEGs (up vs. down, 1:3) were detected in the lung adenocarcinoma samples, whereas 204 upregulated DEGs and 285 downregulated DEGs (up vs. down, 1:1.5) were detected in the lung squamous cell carcinoma samples, as compared with the normal lung tissue samples. The majority of DEGs were downregulated in both lung cancer groups.

Significantly enriched pathways in the DEGs from the two lung cancer groups

In order to determine the differences between the two types of lung cancer at the level of related pathways, DAVID was used to conduct pathway enrichment. The DEGs in the lung squamous cell carcinoma group were significantly enriched in the following three pathways: Hsa04110, cell cycle (P=3.93E-06), involving genes YWAQ, MCM3, CHEK2, GADD45 G, RAD21 and PCNA; hsa03030, DNA replication (P=8.50E-06), involving genes MCM3, FEN1, RFC4, POLA1, PCNA and MCM2; and hsa03430, mismatch repair (P=0.01309), involving genes RFC4, MSH6, RFC2 and MSH3.

The DEGs may be involved in the lung adenocarcinoma group by participating in metabolism pathways, including Hsa00010, glycolysis/gluconeogenesis (ADH1B, ADH1C, ALDH3A2, LDHA, AKR1A1 and FBP1) and hsa00980, metabolism of xenobiotics by cytochrome P450 (ADH1B, GSTA3, CYP2F1 and ADH6). However, no significant pathways were enriched when using the cut-off value of BH-corrected P<0.05.

Global network of cancer-associated genes

Annotated genes from the 14 cancer-associated pathways downloaded from KEGG were used to construct the global network. There were 341 genes and 1,569 edges (interactions between genes) in the network, as presented in Fig. 1A and B.

Figure 1

Global network of cancer-associated genes. Cancer network of lung adenocarcinoma samples. Red dots represent upregulated genes; green dots represent downregulated genes. The depth of color represents the extent of the fold change.

Two upregulated DEGs and seven downregulated DEGs in the adenocarcinoma group, and eight upregulated DEGs and five downregulated DEGs in the squamous cell carcinoma group were selected in the global network under the FDR<0.1 and FC >1.5 or <0.67 criteria. These DEGs were used to construct the PPI network based on the HRPD database; the network consisted of 21 DEGs and 517 relationships with two common genes, NKX3-1 and MAPK10) in Fig. 2A and B. The degrees (number of interacting partners) of the genes in the PPI networks were calculated (Table I), and the top five genes were identified, as follows: HSP90AA1, BCL2, CDK2, KIT and HDAC2, which had degrees of 81, 64, 63, 50 and 46, respectively (Table II). HSP90AA1 exhibited the highest degree, and was the most downregulated gene detected in the adenocarcinoma group. Furthermore, the expression differences of these 19 DEGs were compared between the two types of lung cancer (Table II). The results demonstrated that some genes exhibited similar expression trends between the groups, such as NKX3-1 and MAPK10, whereas DVL3 was upregulated in squamous cell carcinoma samples; however, its expression was not significantly altered in adenocarcinoma samples. The expression of BCL2 was significantly downregulated in adenocarcinoma samples, but not changed in squamous cell carcinoma samples.

Figure 2

Global network of cancer-associated genes. Cancer network of lung squamous cell carcinoma samples. Red dots represent upregulated genes; green dots represent downregulated genes. The depth of color represents the extent of the fold change.

Table I

Degrees of the top 10 genes in the protein-protein interaction network.

Symbol	Degree
HSP90AA1	81
BCL2	64
CDK2	63
KIT	50
HDAC2	46
CDK4	36
FGFR3	30
TGFBR2	28
MAP2K1	27
RAD51	23

Table II

Comparison of 19 differentially expressed genes between the two types of lung cancer.

Symbol	Gene ID	Adenocarcinoma			Squamous cell carcinoma
		Dir	FC	BH-P	Dir	FC	BH-P
FGFR3	2261	−1	0.120	0.006	−1	0.400	0.368
EPAS1	2034	−1	0.503	0.195	−1	0.290	0.040
KIT	3815	−1	0.591	0.565	−1	0.310	0.012
NKX3-1	4824	−1	0.338	0.020	−1	0.348	0.029
TGFBR2	7048	1	1.210	0.769	−1	0.435	0.093
WNT4	54361	−1	0.452	0.017	−1	0.591	0.116
HSP90AA1	3320	−1	0.479	0.027	−1	0.771	0.352
MAPK10	5602	−1	0.510	0.032	−1	0.577	0.064
BCL2	596	−1	0.546	0.056	−1	0.877	0.908
HDAC2	3066	1	1.147	0.806	1	1.751	0.037
CDK2	1017	1	1.310	0.153	1	1.809	0.056
RAD51	5888	1	1.205	0.657	1	1.812	0.048
MAP2K1	5604	1	1.451	0.288	1	1.840	0.060
ITGA3	3675	1	2.120	0.049	1	1.008	0.981
DVL3	1857	−1	0.942	0.912	1	2.279	0.066
CDK4	1019	1	1.644	0.236	1	2.799	0.004
CKS1B	1163	1	1.427	0.754	1	3.061	0.040
DAPK1	1612	1	3.083	0.095	1	1.148	0.749
FZD7	8324	1	1.325	0.792	1	4.668	0.068

FC, fold change; BH-P, Benjamini-Hochberg method modified P-value; dir, direction; −1, downregulated expression; 1, upregulated expression; significance was defined as BH-P<0.1 and FC > 1.5 or < 0.67.

Discussion

In the present study, DEGs in lung adenocarcinoma and lung squamous cell carcinoma samples were screened by comparing gene expression between the cancer samples and normal lung tissue samples. The majority of DEGs were revealed to be downregulated. Compared with lung adenocarcinoma, there were more DEGs in the lung squamous cell carcinoma group. The disheveled homolog DVL3 was upregulated in squamous cell carcinoma and was downregulated in lung adenocarcinoma. Disheveled family proteins have been reported to be associated with the development of cancer, as the cytoplasmic mediators to activate the WNT/β-catenin signaling pathway (16). In addition, positive DVL3 expression has previously been detected in NSCLC samples (17). DVL3 has also been identified as a candidate driver for genomic amplification of chromosome 3q26-29 in lung squamous cell carcinoma of the lung. Knocking down of DVL3 protein can lead to cell growth inhibition (18) Therefore, DVL3 may be considered a reliable biomarker to distinguish between these two types of lung cancer. Furthermore, BCL2 was downregulated in lung adenocarcinoma; however, it was not detected as a DEG in lung squamous cell carcinoma. It has previously been reported that the BCL2 protein acts as an inhibitor of apoptosis, and its family members, Bcl-2 and Bcl-xL, are overexpressed in numerous types of tumor cell, including NSCLC (19). Targeted inhibition of BCL2 overcomes the resistance of NSCLC cell lines to chemotherapy (20). However, recent meta-analysis studies suggest BCL2-negative expression is associated with poor prognosis (21,22), indicating its cancer inhibition roles, which was also confirmed by analysis of lung adenocarcinoma samples in the present study.

A PPI network analysis demonstrated that the top five genes were HSP90AA1, BCL2, CDK2, KIT and HDAC2, indicating these genes may be important for lung cancer. These results were consistent with those of previous studies although their roles require further elucidation. The HSP90AA1 gene (generally referred to as HSP90) is located on chromosome 14q32.2 (23). HSP90 is an emerging focus of cancer therapy due to its ability to simultaneously inhibit several signaling pathways (24). Sequist et al (25) reported that a potent inhibitor of HSP90, IPI-504, exerts clinical activity in patients with NSCLC, particularly among patients with anaplastic lymphoma kinase gene rearrangements. However, HSP90AA1 was observed to be downregulated in the current study, indicating its underlying tumor suppressor effect, which was also demonstrated by Li et al (26) who identified the expression of HSP90AA1 was lower in nasopharyngeal carcinoma radioresistant cells than that in sensitive cells.

CDK2 is a member of the cyclin-dependent kinases, and inhibition of CDK2 has been reported to induce anaphase catastrophe and lead to apoptosis in NSCLC (27). As expected, CDK2 was also upregulated in the lung cancer samples in the present study, particularly in squamous cell carcinoma. Furthermore, Abrams et al (28) evaluated the activity of the indolinone kinase inhibitor SU11248 against the receptor tyrosine kinase KIT in SCLC; the results suggested that SU11248 may have clinical potential for the treatment of SCLC via direct KIT-mediated antitumor activity. Nevertheless, KIT was downregulated in squamous cell carcinoma in the present study, which may be attributed to the lower malignancy in NSCLC compared with SCLC. Histone deacetylases (HDACs) have been identified as therapeutic targets due to their regulatory function on DNA structure and organization (29). HDAC inhibitors (HDACIs) represent a novel class of anticancer agents. The HDACIs LBH589 (29), trichostatin A (30) and suberoylanilide hydroxamic acid (31) exert profound anti-growth activity against lung cancer cells. In addition, the other DEGs identified in the present study may be considered potential biomarkers for lung cancer therapy.

In conclusion, the present study identifies various crucial genes, including HSP90AA1, BCL2, CDK2, KIT and HDAC2, that display different expression trends in lung adenocarcinoma and lung squamous cell carcinoma via DEG analysis compared with normal lung tissue and global network construction. These genes may be potential biomarkers for differentiating these two types of lung cancer in clinic.