Open Access

Lung Cancer Signature Biomarkers: tissue specific semantic similarity based clustering of Digital Differential Display (DDD) data

BMC Research Notes20125:617

DOI: 10.1186/1756-0500-5-617

Received: 3 April 2012

Accepted: 23 October 2012

Published: 2 November 2012

Abstract

Background

The tissue-specific Unigene Sets derived from more than one million expressed sequence tags (ESTs) in the NCBI, GenBank database offers a platform for identifying significantly and differentially expressed tissue-specific genes by in-silico methods. Digital differential display (DDD) rapidly creates transcription profiles based on EST comparisons and numerically calculates, as a fraction of the pool of ESTs, the relative sequence abundance of known and novel genes. However, the process of identifying the most likely tissue for a specific disease in which to search for candidate genes from the pool of differentially expressed genes remains difficult. Therefore, we have used ‘Gene Ontology semantic similarity score’ to measure the GO similarity between gene products of lung tissue-specific candidate genes from control (normal) and disease (cancer) sets. This semantic similarity score matrix based on hierarchical clustering represents in the form of a dendrogram. The dendrogram cluster stability was assessed by multiple bootstrapping. Multiple bootstrapping also computes a p-value for each cluster and corrects the bias of the bootstrap probability.

Results

Subsequent hierarchical clustering by the multiple bootstrapping method (α = 0.95) identified seven clusters. The comparative, as well as subtractive, approach revealed a set of 38 biomarkers comprising four distinct lung cancer signature biomarker clusters (panel 1–4). Further gene enrichment analysis of the four panels revealed that each panel represents a set of lung cancer linked metastasis diagnostic biomarkers (panel 1), chemotherapy/drug resistance biomarkers (panel 2), hypoxia regulated biomarkers (panel 3) and lung extra cellular matrix biomarkers (panel 4).

Conclusions

Expression analysis reveals that hypoxia induced lung cancer related biomarkers (panel 3), HIF and its modulating proteins (TGM2, CSNK1A1, CTNNA1, NAMPT/Visfatin, TNFRSF1A, ETS1, SRC-1, FN1, APLP2, DMBT1/SAG, AIB1 and AZIN1) are significantly down regulated. All down regulated genes in this panel were highly up regulated in most other types of cancers. These panels of proteins may represent signature biomarkers for lung cancer and will aid in lung cancer diagnosis and disease monitoring as well as in the prediction of responses to therapeutics.

Keywords

Digital Differential Display (DDD) Lung tissue cancer Semantic similarity biomarker Clustering analysis Multiple bootstrap

Background

Gene expression analysis in the post genomic era through high throughput genomic studies led to identification of enormous candidate genes related to pathophysiological conditions or altered signal transduction. One such freely available high throughput database is ‘Unigene’ (http://www.ncbi.nlm.nih.gov/Unigene/). The Unigene libraries of interest with varying treatment conditions can be digitally ‘pooled’ and compared to control vs. treatment using Digital Differential Display (DDD). It enables the identification of numerical differences in transcript frequency between the individual or pooled Unigene libraries from the various treatment conditions and multiple cDNA libraries. The frequency of each differentially expressed transcripts and their fold change from the pooled libraries have been calculated using Fisher Exact Test. The prioritisation of DDD identification from differentially expressed candidate genes strictly used relative change in the frequency value and its fold change. Apart from DDD, many web tools are freely available to prioritise candidate genes based on the relative change in gene expression profile [1, 2]. The prioritisation of each tool differs due to their different computational approaches [3]. But the process of identifying the most likely tissue specific disease candidate genes from the pool of differentially expressed genes remained difficult [1].

Recent advances in the systems biology have shown promising results in the elucidation of potential biomarkers of phenotype and clinical relevance, particularly in cancer research sphere [46]. These studies were performed using the predictive integration of gene expression data. Different predictive integration strategies have been developed and were used to study the biological information from public repositories [48]. Amongst such strategies, gene products that are biologically and functionally related would maintain similarity, both in their expression profiles and in the Gene Ontology (GO) annotation [9]. The integration of gene expression data and standardised descriptions of the biological function of gene products were used for the search of candidate prognostic biomarkers and therapeutic targets [1012]. These studies demonstrated that the measure of functional similarity based GO annotations between query genes and the genes of interest can be applied as a complementary predictive feature to characterise gene expression profile. So, we have applied this integrative computational approach to characterise a tissue specific biological data from DDD.

We hypothesised that tissue specific differentially expressed genes can be functionally characterised using their GO semantic similarity score with normal tissue specific genes (query genes). The query genes, in this study, were normal lung tissue specific genes from the Tissue-Specific Genes Database (TiSGeD). The genes of interest were candidate lung cancer genes from DDD [13, 14]. Surprisingly, this approach successfully distinguished 38 signature biomarkers for lung cancer. Thus this suggests that, in principle, this integrated methodology can offer a complementary predictive capability for detecting tissue specific signature biomarkers from the tissue specific differentially expressed data. These tissue specific signature biomarkers may be candidate prognostic biomarkers and therapeutic targets for lung cancer.

Methods

Selection of Human Lung Tissue specific query genes

The normal lung tissue specific genes were collected from TiSGeD (Tissue-specific gene database; http://bioinf.xmu.edu.cn/databases/TiSGeD/index.html). Human adult lung tissue related genes with tissue specificity measure score (SPM) ≥ 0.9 (represents high tissue specificity) were considered. The lung tissue specific “Mouse” and “developmental” genes were omitted.

Collection of Lung Tissue specific differentially Expressed Candidate Genes using DDD

DDD comparisons were made at various tissue stages to elucidate the selective differential expression levels of human lung tissue specific genes for normal (Case 1) and cancerous (Case 2) conditions. In Case 1, the Normal lung tissues (11 tissue libraries) were considered as a ‘Reference’ samples and the remaining normal human tissues (251 tissue libraries) were ‘Query’ samples. In Case 2, the Normal lung tissues (11 tissue libraries) were considered as ‘Reference’ samples and the cancerous human lung tissues (8 tissue libraries) were ‘Query’ samples. These comparisons were designed systematically so as to identify altered Gene expression of varying treatment conditions of ‘Reference’ and ‘Query’ samples. These pair wise comparisons resulted in a relative abundance of ESTs among the contrasting cDNA libraries of digitally ‘pooled’ contracts from Unigene Database.

GO-based similarity assessment

Org.Hs.eg.db package in R-program was used for the computation of Semantic similarity score while the GO-based similarity score was computed based on the three orthogonal gene ontologies generated for Molecular Function (MF), Cellular Component (CC) and Biological process (BP). GOSemSim of R-program was used to calculate semantic similarity between the GO terms and the gene products. In this study, GO terms derived from human annotations were used for calculations. The estimation of between-term similarity was based on the Wang semantic similarity measure [12]. Aggregation of between-term similarities was done with the highest between-term similarity approach, which selectively aggregates maximum between-gene similarity values [9]. Given a pair of gene products, gi and gj, annotated to a set of GO terms, the GO-driven similarity, SIM (gi, gj), is calculated by aggregating the maximum interest similarity values as follows:
S i m g i , g j = 1 i m m a x S i m g i 1 g j + 1 j m m a x S i m g j 1 g i / m + n

where, two sets of GO terms gi = {gi1, gi2, …………., gim} & gj = {gj1, gj2, ……………, gjn} as query and reference sequence. Method max calculates the maximum semantic similarity score over given pairs of GO terms between these two sets, while average calculates the average semantic similarity score over a given pairs of GO terms. The hierarchical clustering of tissue specific, differentially expressed genes in relation with a normal lung tissue is shown in a Dendrogram. In the colour code of heat map, red represents a low semantic similarity below the median level, whereas, the green represents a high semantic similarity above the median level.

Clustering analysis

The clustering analysis was carried out by the program pvclust [15]. It is an add-on package for a statistical software R to perform the bootstrap analysis of clustering and also to assess the uncertainty in hierarchical cluster analysis. The package calculates the approximately unbiased (AU) and bootstrap probability (BP) p-values for each cluster. Stability of the clustering was accessed at 95% probability (α = 0.95).

Results

DDD based prioritisation of lung cancer genes

In order to find the lung tissue specific differentially expressed genes, two Unigene pools (A and B) were constructed (See Additional file 1). For analysis, in the DDD1, we employed the UniGene pool (A) representing 39 human normal tissues excluding normal lung tissue and UniGene pool (B) representing 11 counterpart lung normal tissues were employed for analysis (Table 1). Similarly, in DDD2, UniGene pool (A) representing 8 human lung tumours and UniGene pool (B) representing 11 counterpart lung normal tissues were employed (Table 1). The fold change of normal lung (DDD1) and lung carcinoma candidate genes (DDD2) were calculated based on transcript frequency values. The candidate genes with an expression of at least 2-fold difference were taken into analysis. In DDD1, amongst the total of 519 differentially expressed genes 268 genes were up-regulated (≥2-fold) and 234 genes were down-regulated (≥2-fold). In DDD2, amongst the total of 203 differentially expressed candidate genes, 147 genes (≥2-fold) including 33 unknown were up-regulated (≥2-fold) and 55 genes were down-regulated (≥2-fold). Comparison of DDD1 with DDD2 has revealed that in total 76 genes from DDD1 were differentially expressed in DDD2 (See Additional file 2). From the literature survey, amongst the 76 genes, 18 of them were found to be commonly expressed in all types of cancerous conditions (See Additional file 3) [16]. Excluding these 18 from the 76, the remaining 58 genes were predicted as the lung tissue specific tumour genes (See Additional file 2). The molecular functions of these 58 genes were found to be involved in broad range of cellular functions with majority of the genes playing many different roles like structural, extracellular and intracellular functions. This subtractive approach eliminated most of the commonly expressing genes; for example, housekeeping genes. This approach has also helped to eliminate genes expressing in more than 10 cancerous conditions (See Additional file 4).
Table 1

Different tissue specific Unigene libraries employed in DDD

S. No.

Tissuetypes

No. of ESTlibraries

dbEST identification number

Respective Number of ESTs

1

Adipose

2

10983, 16445

8299, 1646

2

Adrenal gland

4

6791, 6792, 16377, 18302

4582, 1425, 2756, 10026

3

Tongue

4

12982, 18362, 18389, 18479

1116, 7730, 23564, 7420

4

Bladder

1

18307

8220

5

Blood

6

7037, 7038, 8975, 11923, 6824, 9724

1524, 1721, 4215, 6553, 7241, 9352

6

Bone

3

1124, 821, 16433

6209, 1337, 2615

7

Bone marrow

8

6975, 6976, 15949, 15950, 16412, 931, 10409, 10410

2459, 3424, 1379, 2707, 3623, 5336, 1798, 5231

8

Brain

36

1749, 16376, 16390, 18317, 18318, 18352, 18353, 18415, 18466, 14591, 14592, 186, 17380, 742, 16380, 18310, 18311, 16382, 16383, 18322, 1918, 5655, 6811, 6812, 8570, 19377, 14298, 13711, 13053, 536, 7209, 16384, 18319

1600, 1361, 13592, 2579, 15152, 43177, 24735, 23791, 6112, 31569, 15839, 1461, 1334, 18126, 5838, 4298, 3364, 41751, 40272, 2802, 3192, 16612, 7369, 8197, 2790, 3145, 19115, 44785, 1609, 7033, 3931, 4897, 5194, 2565, 25752

9

Eye

26

13915, 17747, 7316, 7315, 10273, 10287, 13901, 19465, 10274, 10281, 10280, 10288, 10279, 10284, 19471, 10285, 10286, 12093, 302, 303, 433, 10966, 10282, 10283, 16572, 13902

3739, 4253, 1836, 1294, 4005, 3543, 2785, 1334, 1595, 1479, 1115, 6010, 1469, 6719, 3043, 8344, 1199, 7816, 9190, 1732, 2174, 4531, 6279, 1185, 6097, 2946

10

Heart

8

15951, 16399, 18354, 18410, 16421, 18503, 16379, 16381

5307, 3284, 2667, 8670, 4000, 8502, 7220, 4698

11

Stomach

12

10299, 10301, 10302, 10305, 10306, 10310, 10311, 10324, 10325, 18488, 18529, 16432

1793, 2409, 1453, 2790, 1692, 1984, 6125, 5913, 1422, 8604, 2574, 3137

12

Testis

5

1752, 16441, 18476, 18517, 19376

6624, 2983, 46964, 44057, 40315

13

Thalamus

3

16437, 18348, 18349

3154, 29651, 23010

14

Thymus

6

16440, 18518, 18519, 18520, 13049, 18375

2365, 1044, 31967, 37541, 3477, 15983

15

Thyroid

3

889, 16408, 7004

1357, 4827, 3342

16

Breast

4

894, 895, 18305, 18475

1786, 6346, 2538, 8256

17

Cartilage

2

8936, 8940

4310, 3858

18

Cervix

2

18425, 18506

2674, 2619

19

Ear

2

371, 18222

12666, 3396

20

Intestine

13

840, 841, 842, 882, 16385, 18350, 18489, 16387, 16400, 17427, 18473, 16425, 18486

1499, 1704, 1759, 11996, 1817, 8191, 2619, 1351, 5536, 8199, 16855, 3403, 2545

21

Epidermis

4

20865, 21612, 7269, 21098

1627, 13186, 10681, 1135

22

Lung

11

10395, 16406, 16413, 16438, 18355, 18363, 18521, 10398, 11912, 18522, 18537

12545, 2448, 6839, 3327, 2565, 16156, 2677, 11510, 15695, 32590, 19278

23

Liver

11

1365, 12531, 12532, 12535, 12549, 12550, 13859, 16392, 18416, 18525, 18893

2302, 1607, 2315, 1136, 2425, 1561, 7537, 6856, 6550, 8424, 31921,

24

Mammary gland

3

6982, 16420, 16436

4561, 3502, 3371

25

Lymph

8

2709, 2710, 2711, 3718, 3719, 3720, 10312, 8613

3434, 3963, 1000, 9867, 1556, 1828, 7590, 1949

26

Teeth

1

12639

1576

27

Medulla

1

9725

9919

28

Muscle

4

530, 16391, 45, 18501

4271, 2154, 2485, 8276

29

Ovary

7

887, 6998, 18421, 18527, 5444, 12637, 12638

2341, 3678, 2300, 2543, 10294, 1050, 1031

30

Nose

2

358, 13908

1702, 24528

31

Placenta

14

13037, 16442, 507, 740, 6999, 10403, 10404, 10405, 10424, 10425, 16422, 17682, 18468, 18484, 13000

11885, 3501, 1370, 1172, 20941, 1862, 1166, 4260, 4517, 1271, 1529, 1032, 16852, 15843, 7344

32

Pancreas

6

16423, 422, 16960, 8840, 3884, 9821

4308, 1799, 13791, 60665, 1234, 17183

33

Prostate

8

888, 17392, 18469, 19880, 16424, 924, 928

1014, 1283, 16483, 41945, 2355, 7609, 1114, 1051

34

Uterus

6

1753, 16443, 18523, 18531, 18544, 12528

1645, 5121, 30124, 2486, 19168, 3510

35

Spleen

2

16431, 18474

2717, 33972

36

Salivary gland

1

16430

2336

37

Kidney

7

16393, 16395, 16410, 18374, 18377, 18524, 16429

5486, 5150, 3565, 17079, 2561, 15730, 1170

37

Pituitary

4

6828, 6829, 13019, 13737

3481, 1444, 7327, 1677

38

immune cells

6

1317, 17555, 17556, 8892, 12072, 12798

11447, 15272, 13685, 6112, 9223, 4548

39

Hair

4

21096, 21100, 21099, 21101

1269, 1288, 1298, 1543

40

Alimentary canal

2

18496, 18418

8371, 2586

41

Cancerous lung

8

1533, 10419, 537, 914, 14132, 14133, 14134, 14135

1377, 4401, 4173, 1850, 4365, 2121, 2847, 1560

Differentially expressed normal lung tissue specific genes (DDD1) were identified from UniGene libraries representing 39 human normal tissues (251 libraries, 1903301 ESTs, S. No. 1–21, 23–40) and counterpart normal lung tissue (11 libraries, 125630 ESTs, S. No. 22 only). Differentially expressed lung tissue specific genes (DDD2) were identified from UniGene libraries representing lung cancer libraries (8 libraries, 22694 ESTs, S. No. 41 only) and counterpart normal lung tissue (11 libraries, 125630 ESTs, S. No. 22 only).

Prediction of Lung tissue specific tumour genes by Semantic similarity score based clustering

To identify lung tissue specific clusters from the 202 genes from DDD2 cancerous condition, firstly they were subjected to similarity clustering analysis using the 47 lung tissue specific genes from TiSGeD (See Additional file 5). Before the semantic similarity clustering analysis, the Unigene ID were converted into Entrez ID. During this process, the 202 genes of DDD2 reduced to 145 and the 47 lung tissue specific genes of TiSGeD were reduced to 28 due to gene duplication. Using GOSemSim package, the similarity correlation matrix was constructed between the 145 predicted lung specific differentially expressed cancer genes from DDD2 and 28 genes from TiSGeD. The differential expression levels of these clustered genes were depicted in the form of a Heat Map (Figure 1). The similarity correlation matrix produced seven gene clusters at 95% confidence level, using the pvclust program (Figure 2). The clusters 1–4 have 14 genes and the clusters 5, 6 and 7 have 36, 74 and 14 genes respectively.
https://static-content.springer.com/image/art%3A10.1186%2F1756-0500-5-617/MediaObjects/13104_2012_Article_1917_Fig1_HTML.jpg
Figure 1

Go semantic similarity score between the set of normal lung tissue specific genes from TiSGeD (28-horizontal, x-axis) and the differentially expressed lung cancer genes from DDD2 (145-vertical, y-axis). The intensity of the color corresponds to the magnitude of the similarity. Red represents low semantic similarity below the median level whereas the green represents high semantic similarity above the median level.

https://static-content.springer.com/image/art%3A10.1186%2F1756-0500-5-617/MediaObjects/13104_2012_Article_1917_Fig2_HTML.jpg
Figure 2

Average correlation distances with hierarchical clustering based on GO semantic similarity score matrix calculated between normal lung tissue specific genes from TiSGed and differentially expressed lung cancer gene from DDD2. Values in red represent AU (Approximately unbiased) p-value and green represents BP (Bootstrap probability) Clusters with AU larger than 95% are highlighted by red rectangle boxes. AU p-value, which is computed by multiscale bootstrap resampling, is a better approximation to unbiased p-value than BP value computer by normal bootstrap resampling.

In the ID conversions from Unigene to Entrez, the 58 lung tissue specific tumour genes were reduced to 38 genes (Table 2). These 38 genes were matched with the 7 clusters. This38 genes formed four panels with the corresponding cluster 4, 5, 6 and 7 respectively. The panels 1–4 contained 2, 9, 21 and 6 genes respectively. This leads to identification of the lung tissue specific clusters of the normal lung tissue specific genes differentially regulated in lung cancer condition.
Table 2

Lung cancer signature biomarker clusters

Clustering of lung tissue specific cancer Biomarkers

DDD1

DDD2

Panel / Cluster of Biomarker

Type / class of Biomarker

GeneSymbol

Transcript frequency

Foldchange

Transcript frequency

Foldchange

   

Pool A (normallung tissues)

Pool B (Normalbody tissuesexcluding lung)

 

Pool A (normallung tissues)

Pool B (lungtissue specificcancers)

 

Panel 1 / Cluster 4

Lung cancer metastasisdiagnostic markers

UCHL1

0

0.0002

-

0

0.0011

+

  

LTF

0.0224

0.0018

+12.44

0.0224

0.0002

-112

Panel 2 / Cluster 5

Chemotherapy/ drugresistance related lungcancer biomarkers

TUBA1B

0

0.0002

-

0

0.0013

+

  

RPSA

0.0001

0.0004

-4

0.0001

0.0027

+27

  

RPL9

0.0002

0.0006

-3

0.0002

0.002

+10

  

TMSB4X

0.0004

0.001

-2.5

0.0004

0.0016

+4

  

COPB1

0.0007

0.0002

+3.5

0.0007

0

-

  

API5

0.0007

0.0003

+2.3

0.0007

0

-

  

NT5C2

0.0008

0.0003

+2.6

0.0008

0

-

  

CPN

0.0009

0.0001

+9

0.0009

0

-

  

PRKAR1A

0.0017

0.0006

+2.83

0.0017

0

-

Panel 3 / Cluster 6

Hypoxia related lungcancer biomarkers

FTL

0.0001

0.0011

-11

0.0001

0.0065

+65

  

COL1A2

0.0001

0.0006

-6

0.0001

0.0023

+23

  

GAPDH

0.0001

0.001

-10

0.0001

0.0011

+11

  

IGKC

0.0002

0.0009

-4.5

0.0002

0.0016

+8

  

ALDOA

0.0002

0.0006

-3

0.0002

0.0014

+7

  

COL1A1

0.0001

0.0004

-4

0.0001

0.0009

+9

  

FN1

0.0025

0.0012

+2.08

0.0025

0.0007

-3.57

  

TGM2

0.0026

0.0008

+3.25

0.0026

0.0007

-3.71

  

FOS

0.0015

0.0007

+2.14

0.0015

0.0002

-7.5

  

CTNNA1

0.0034

0.0008

+4.25

0.0034

0.0003

-11.33

  

FOSB

0.0024

0.0003

+8

0.0024

0.0002

-12

  

APLP2

0.0109

0.0044

+2.47

0.0109

0.0008

-13.63

  

NCOA4

0.0016

0.0005

+3.2

0.0016

0.0001

-16

  

HIF1A

0.0011

0.0004

+2.75

0.0011

0

-

  

AZIN1

0.001

0.0005

+2

0.001

0

-

  

EHF

0.001

0.0001

+10

0.001

0

-

  

TICAM2

0.001

0.0003

+3.33

0.001

0

-

  

NAMPT

0.0008

0.0002

+4

0.0008

0

-

  

TNFRSF1A

0.0008

0.0004

+2

0.0008

0

-

  

DMBT1

0.0008

0.0001

+8

0.0008

0

-

  

CSNK1A1

0.0007

0.0003

+2.33

0.0007

0

-

Panel 4 / Cluster 7

Lung cancer specificextra cellular matrixbiomarkers

TFPI2

0

0.0004

-

0

0.002

+

  

RPL10

0.0001

0.0007

-7

0.0001

0.0014

+14

  

SFTPA1

0.0004

0

+

0.0001

0.0011

+11

  

KIAA1324

0.0013

0.0002

+6.5

0.0013

0

-

  

CRISP3

0.0009

0.0001

+9

0.0009

0

-

  

NET1

0.0007

0.0002

+3.5

0.0007

0

-

Note: (−) represents down regulation and (+) represents up regulation. Except SFTPA1, remaining all the up regulated genes in normal condition observed to be down regulated in cancerous condition and vice versa.

We then analysed the functional significance of each panel as given below.

Analysis of Cluster 4 / Panel 1

The cluster 4 had two-lung cancer related genes ubiquitin thiolesterase (UCHL1) and Lactotransferrin (LTF). In the normal lung (DDD1 data), UCHL1 was down-regulated and LTF was up-regulated (Table 2). This was reversed during the lung cancer condition where UCHL1 up-regulated and the LTF highly down-regulated (Table 2). These two proteins were found to be important in the cancer progression. UCH-L1 up-regulation promoted prostate cancer metastasis through epithelial-to-mesenchymal transition (EMT) induction and LTF expression decreased in lung prostate cancer progression [17, 18]. Both of them were co-expressed in almost six different lung adenocarcinoma cell lines, as evident by mSigDB. This suggested that UCH-L1 and LTF could be novel diagnostic and therapeutic targets for lung cancer metastasis diagnostic markers.

Analysis of Cluster 5 / Panel 2

The cluster 5 was playing the common functional role of immune response and complement activation. The down-regulated RPSA, RPL9, TMSB4X and TUBA1B in normal lung (DDD1) were significantly up-regulated in lung cancer (DDD2) (Table 1). The analysis resulted that all these up regulated genes played the role of tumour cell resistance to the anti-cancer agents. In gastric cancers, the up-regulation of RPSA/LRP contributed to drug resistance via hypoxia-inducible-factor dependent mechanism [19]. Similarly, there was a link between the TMSB4X and TUBA1B and the anti-cancer drug resistance to the drug Paclitaxel (PTX) observed in the cervical and breast/ovarian cancers respectively [20, 21].

In this cluster, NT5C2, API5, CPN, PRKAR1A and COPB1 were fully down-regulated in lung cancer (Table 1). The down regulation of NT5C3 altered the tumour cell sensitivity to cytidine based anti-cancer drugs [22]. The anti-apoptosis gene API5 down-regulation linked to increase in the survival and resistance cancer cells to chemotherapy [23]. To our knowledge, the major copper carrying protein CPN (ceruloplasmin) down regulation link to chemotherapy/drug resistance is not yet studied. But increased level of copper in lewis lung carcinoma cells were related with the development of multi drug resistance [24]. The PRKAR1A down-regulation also linked to multidrug-resistant (MDR) in colon carcinoma cells [25]. The COPB1 was an essential component for the coatomer formation [26]. These coatomers were involved in the drug trafficking pathways and endocytic drug delivery [27]. So, it was expected that the down-regulation of COPB1 might have a role in the chemotherapy which needs to be taken up and studied. We are surprised to find that all these results suggest that the cluster 5 functionally represents a panel of chemotherapy/drug resistance related lung cancer biomarkers.

Analysis of Cluster 6 / Panel 3

In cluster 6, the upregulated FTL (65 fold in our study) and ALDOA (7 fold in our study) were regulated by hypoxia inducible factor (HIF) during lung cancer [2831]. The COL1A1 (23 fold in our study) and GAPDH (11 fold in our study) were regulated by hypoxia [3234]. IGKC (8 fold in our study) up-regulated in lung cancer patients but no literature data was available for its interaction either with HIF or hypoxia [35]. The HIF, TGM2, CSNK1A1, CSNK2A1, CTNNA1, NAMPT)/Visfatin, TNFRSF1A, ETS1 and SRC-1 were down-regulated and proposed as the biomarkers for lung cancer. We found all of them to be interacting with the HIF in cancerous condition [3645]. The down- regulated FN1 and APLP2 showed hypoxia dependent differential regulation [4648]. The DMBT1/SAG interacted HIF-1 was a kind of feedback loop in response to hypoxia. The hypoxia induced HIF-1 to transactivate SAG and the induced SAG then promoted HIF-1alpha ubiquitination and degradation [49]. The FBJ/c-Jun/AP-1 interacted with HIF during hypoxia that controlled the transcriptional regulation of the Cyr61 gene in retinal vascular endothelial cells [50]. The role of AIB1/SRC-3/NCoA during hypoxia condition were exhibited by controlling the expression levels of HIF induced erythropoietin (EPO) gene during hypoxia [42].

However, in this cluster, the AZIN1 and TICAM2 were down-regulated and were lacking direct experimental evidence to support their regulation with HIF or hypoxia during cancer. The following literature analysis suggests their possible regulations either with HIF or hypoxia. The AZIN1 was an inhibitor for the antizyme and both were highly regulated in human cancers and antizyme induced HIF, during increased cellular redox potential [5153]. The TICAM2 physically bridged toll like receptor-4 (TLR4) with TICAM1 and the TLR4 partially regulated by the HIF during adenocarcinoma [54, 55].

All these results suggest that the cluster 6 represents the panel of either HIF or Hypoxia related lung cancer biomarkers.

Analysis of Cluster 7 / Panel 4

In the Cluster 7, there were seven lung biomarkers, mostly encoding for lung tissue specific extra cellular matrix proteins. The epigenetic analysis using Methycancer database (http://methycancer.genomics.org.cn) revealed that amongst the seven, KIAA1324, NET1, NTN3, RPL10 and TFPI2 were epigenetically regulated through DNA methylation. In the remaining two, SFTPA1 was epigenetically regulated [5658]. However, the experimental evidence was lacking the epigenetic related data for CRISP3. However, the Gene card database analysis of CRISP3 showed that the CRISP3 orthlogous gene C-type lectin domain family 18 member A (CLEC18A) epigenetically regulated through DNA methylation (http://www.genecards.org). All these results show that the cluster 7 represented the panel of epigenetically regulated lung cancer specific extra cellular matrix biomarkers.

Discussion

UniGene database using the DDD tool provides us a computational approach to study and understand the lung tissue specific gene expression levels in both disease and normal conditions [59]. Studying their differential expression in disease state (lung cancer) will provide a clue about lung cancer specific candidate genes. However, the candidate identification of the DDD method is relying on the EST frequencies based fold change calculation. In DDD2, the 203 differentially expressed candidate genes (≥2-fold) ranking / prioritisation only based upon fold change did not account for the tissue specific variability of the genes in disease conditions (eg: biomarker identification). To include the tissue specific variability in DDD2 prioritisation, the normal lung tissue specific genes from DDD1 were compared. This approach eliminated most of the house keeping genes from the analysis (gene list reduced from 202 to 76). Further, we detected genes expression selectively altered in the lung cancer by eliminating genes that commonly expressed differentially in more than five tumours (gene list reduced from 76 to 58) (See Additional file 2). Almost all of them have a documented role in the lung cancer (http://www.megabionet.org/bio/hlung). So, these subtractive approaches successfully increase the probability of identifying the lung cancer specific probable candidate biomarkers.

The semantic similarity scores amongst the GO terms and the subsequent hierarchical clustering were calculated using the freely available R-software for lung tissue specific candidate genes from normal and cancer conditions. The analysis of members of individual genes from each cluster revealed the functional significance of each cluster. Out of the seven clusters, our approach identified four functionally important clusters. The four clusters represented metastasis diagnostic markers, chemotherapy/drug resistance related biomarkers, and HIF or Hypoxia induced biomarkers and epigenetically regulated extra cellular matrix biomarkers for lung cancer. This suggests that, especially for lungs tissues, the semantic similarity score amongst GO terms between normal and diseases condition from the same tissue can prioritise biomarkers. But, further study is necessary to extend our hypothesis to other tissues. This subtractive approach integrated with semantic similarity score among GO terms can offer a predictive capability for detecting tissue specific signature biomarkers from the tissue specific differentially expressed data. This approach is also complementary to the network based biomarker prediction approach [60, 61]. Our study is one more example of demonstrating the utility of the Digital differential expression technique.

Our study suggests that amongst the 4 panels, HIF or Hypoxia induced lung cancer biomarkers panel (panel 3) is the most important cluster. Because, in other clusters, most of the identified lung cancer biomarkers follow the same expression pattern (either up or down) in other types cancers like breast, ovarian, cervical etc. However, in our study and literature, the expression pattern of genes down regulated in cluster 6 / panel 3 is distinct from almost all types of other cancers. In panel 3, the expression pattern of the HIF and its modulating proteins are completely different when compared to most of the other types of cancers. For example, in most of the cancerous conditions the HIF level is up-regulated [62]. This up-regulation is expected in cancers due to the acute hypoxic condition exhibited during cancer. In contrast, in lung cancer, the HIF level is completely down-regulated (Table 1).

Therefore, it is evident from our study that the HIF down regulation also affect the expression level of the other HIF modulating lung cancer biomarkers. All the down-regulated genes, in this Panel 3 showed their significant up-regulation in most of many types of cancers (TGM2 [63, 64], CSNK1A1 [65], CTNNA1 [66], NAMPT/Visfatin [67], TNFRSF1A [68], ETS1 [41], SRC-1 [69], FN1 [70], APLP2 [71], DMBT1/SAG [64], AIB1 [72], AZIN1 [72]). Our study further shows that this down-regulation is more than five folds when compared to the normal lungs tissue (Table 1). This fold change level suggests that this fold change seems to be more than enough to detect them in the patient sample. Therefore, this panel of down regulating HIF / hypoxia regulated lung cancer biomarker can help to distinguish lung cancer from other types of cancers.

The identified 38 signature lung cancer specific biomarkers can help to increase the sensitivity and selectivity for early diagnosis of lung cancer.

Conclusion

We could demonstrate that our approach readily predicted lung tissue specific cancer biomarkers from digital differentially expressed lung cancer tissue specific genes. The procedure can easily adapt for the prediction of tissue specific biomarkers from the tissue specific differentially expressed genes. It is necessary to explore the extent to which the proposed approach can be integrated with the prediction of tissue specific biomarkers from tissue specific microarray datasets.

Abbreviations

ESTs: 

Expression Sequence Tags

DDD: 

Digital Differential Display

GO: 

Gene Ontology

TiSGeD: 

Tissue-Specific Genes Database

EMT: 

epithelial-to-mesenchymal transition.

Declarations

Acknowledgements

This work was funded by Defence Research and Development organization (DRDO), India fellowship grants.

Authors’ Affiliations

(1)
Bioinformatics Group, Defence Institute of Physiology and Allied Sciences

References

  1. Zhu M, Zhao S: Candidate gene identification approach: progress and challenges. Int J Biol Sci. 2007, 3 (7): 420-427.PubMedPubMed CentralView ArticleGoogle Scholar
  2. Chen J, Bardes EE, Aronow BJ, Jegga AG: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009, 37 (Web Server issue): W305-W311.PubMedPubMed CentralView ArticleGoogle Scholar
  3. Tranchevent LC, Capdevila FB, Nitsch D, De Moor B, De Causmaecker P, Moreau Y: A guide to web tools to prioritize candidate genes. Brief Bioinform. 2011, 12 (1): 22-32. 10.1093/bib/bbq007.PubMedView ArticleGoogle Scholar
  4. Chuang HY, Hofree M, Ideker T: A decade of systems biology. Annu Rev Cell Dev Biol. 2010, 26: 721-744. 10.1146/annurev-cellbio-100109-104122.PubMedPubMed CentralView ArticleGoogle Scholar
  5. Chen J, Sam L, Huang Y, Lee Y, Li J, Liu Y, Xing HR, Lussier YA: Protein interaction network underpins concordant prognosis among heterogeneous breast cancer signatures. J Biomed Inform. 2010, 43 (3): 385-396. 10.1016/j.jbi.2010.03.009.PubMedPubMed CentralView ArticleGoogle Scholar
  6. Azuaje F: What does systems biology mean for biomarker discovery?. Expert opinion on medical diagnostics. 2010, 4: 1-10. 10.1517/17530050903468709.PubMedView ArticleGoogle Scholar
  7. Taylor IW, Linding R, Warde-Farley D, Liu Y, Pesquita C, Faria D, Bull S, Pawson T, Morris Q, Wrana JL: Dynamic modularity in protein interaction networks predicts breast cancer outcome. Nat Biotechnol. 2009, 27 (2): 199-204. 10.1038/nbt.1522.PubMedView ArticleGoogle Scholar
  8. Azuaje F, Devaux Y, Wagner DR: Coordinated modular functionality and prognostic potential of a heart failure biomarker-driven interaction network. BMC Syst Biol. 2010, 4: 60-10.1186/1752-0509-4-60.PubMedPubMed CentralView ArticleGoogle Scholar
  9. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S: GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics. 2010, 26 (7): 976-978. 10.1093/bioinformatics/btq064.PubMedView ArticleGoogle Scholar
  10. Bolshakova N, Azuaje F, Cunningham P: A knowledge-driven approach to cluster validity assessment. Bioinformatics. 2005, 21 (10): 2546-2547. 10.1093/bioinformatics/bti317.PubMedView ArticleGoogle Scholar
  11. Pesquita C, Faria D, Falcao AO, Lord P, Couto FM: Semantic similarity in biomedical ontologies. PLoS Comput Biol. 2009, 5 (7): e1000443-10.1371/journal.pcbi.1000443.PubMedPubMed CentralView ArticleGoogle Scholar
  12. Wang H, Zheng H, Browne F, Glass DH, Azuaje F: Integration of Gene Ontology-based Similarities for Supporting Analysis of Protein-Protein Interaction Networks. Pattern Recognit Lett. 2010, 31: 2073-2082. 10.1016/j.patrec.2010.04.011.View ArticleGoogle Scholar
  13. Xiao SJ, Zhang C, Zou Q, Ji ZL: TiSGeD: a database for tissue-specific genes. Bioinformatics. 2010, 26 (9): 1273-1275. 10.1093/bioinformatics/btq109.PubMedPubMed CentralView ArticleGoogle Scholar
  14. Liu F, Wang H, Li J: An integrated bioinformatics analysis of mouse testis protein profiles with new understanding. BMB Rep. 2011, 44 (5): 347-351. 10.5483/BMBRep.2011.44.5.347.PubMedView ArticleGoogle Scholar
  15. Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22 (12): 1540-1542. 10.1093/bioinformatics/btl117.PubMedView ArticleGoogle Scholar
  16. Chen S, Zhu B, Yu L: In silico comparison of gene expression levels in ten human tumor types reveals candidate genes associated with carcinogenesis. Cytogenet Genome Res. 2006, 112 (1–2): 53-59.PubMedView ArticleGoogle Scholar
  17. Jang MJ, Baek SH, Kim JH: UCH-L1 promotes cancer metastasis in prostate cancer cells through EMT induction. Cancer Lett. 2011, 302 (2): 128-135. 10.1016/j.canlet.2011.01.006.PubMedView ArticleGoogle Scholar
  18. Shaheduzzaman S, Vishwanath A, Furusato B, Cullen J, Chen Y, Banez L, Nau M, Ravindranath L, Kim KH, Mohammed A: Silencing of Lactotransferrin expression by methylation in prostate cancer progression. Cancer Biol Ther. 2007, 6 (7): 1088-1095. 10.4161/cbt.6.7.4327.PubMedView ArticleGoogle Scholar
  19. Liu L, Sun L, Zhang H, Li Z, Ning X, Shi Y, Guo C, Han S, Wu K, Fan D: Hypoxia-mediated up-regulation of MGr1-Ag/37LRP in gastric cancers occurs via hypoxia-inducible-factor 1-dependent mechanism and contributes to drug resistance. Int J Cancer. 2009, 124 (7): 1707-1715. 10.1002/ijc.24135.PubMedView ArticleGoogle Scholar
  20. Banerjee A: Increased levels of tyrosinated alpha-, beta(III)-, and beta(IV)-tubulin isotypes in paclitaxel-resistant MCF-7 breast cancer cells. Biochem Biophys Res Commun. 2002, 293 (1): 598-601. 10.1016/S0006-291X(02)00269-3.PubMedView ArticleGoogle Scholar
  21. Moon EY, Im YS, Ryu YK, Kang JH: Actin-sequestering protein, thymosin beta-4, is a novel hypoxia responsive regulator. Clin Exp Metastasis. 2010, 27 (8): 601-609. 10.1007/s10585-010-9350-z.PubMedView ArticleGoogle Scholar
  22. Li L, Fridley B, Kalari K, Jenkins G, Batzler A, Safgren S, Hildebrandt M, Ames M, Schaid D, Wang L: Gemcitabine and cytosine arabinoside cytotoxicity: association with lymphoblastoid cell expression. Cancer Res. 2008, 68 (17): 7050-7058. 10.1158/0008-5472.CAN-08-0405.PubMedPubMed CentralView ArticleGoogle Scholar
  23. Goswami S, Wang W, Wyckoff JB, Condeelis JS: Breast cancer cells isolated by chemotaxis from primary tumors show increased survival and resistance to chemotherapy. Cancer Res. 2004, 64 (21): 7664-7667. 10.1158/0008-5472.CAN-04-2027.PubMedView ArticleGoogle Scholar
  24. Majumder S, Dutta P, Choudhuri SK: The role of copper in development of drug resistance in murine carcinoma. Med Chem. 2005, 1 (6): 563-573. 10.2174/157340605774598153.PubMedView ArticleGoogle Scholar
  25. Nesterova MV, Johnson NR, Stewart T, Abrams S, Cho-Chung YS: CpG immunomer DNA enhances antisense protein kinase A RIalpha inhibition of multidrug-resistant colon carcinoma growth in nude mice: molecular basis for combinatorial therapy. Clin Cancer Res. 2005, 11 (16): 5950-5955. 10.1158/1078-0432.CCR-05-0624.PubMedView ArticleGoogle Scholar
  26. Nickel W, Brugger B, Wieland FT: Vesicular transport: the core machinery of COPI recruitment and budding. J Cell Sci. 2002, 115 (Pt 16): 3235-3240.PubMedGoogle Scholar
  27. Watson P, Jones AT, Stephens DJ: Intracellular trafficking pathways and drug delivery: fluorescence imaging of living and fixed cells. Adv Drug Deliv Rev. 2005, 57 (1): 43-61. 10.1016/j.addr.2004.05.003.PubMedView ArticleGoogle Scholar
  28. Kukulj S, Jaganjac M, Boranic M, Krizanac S, Santic Z, Poljak-Blazi M: Altered iron metabolism, inflammation, transferrin receptors, and ferritin expression in non-small-cell lung cancer. Med Oncol. 2010, 27 (2): 268-277. 10.1007/s12032-009-9203-2.PubMedView ArticleGoogle Scholar
  29. Smith TG, Balanos GM, Croft QP, Talbot NP, Dorrington KL, Ratcliffe PJ, Robbins PA: The increase in pulmonary arterial pressure caused by hypoxia depends on iron status. J Physiol. 2008, 586 (Pt 24): 5999-6005.PubMedPubMed CentralView ArticleGoogle Scholar
  30. Rho JH, Roehrl MH, Wang JY: Glycoproteomic analysis of human lung adenocarcinomas using glycoarrays and tandem mass spectrometry: differential expression and glycosylation patterns of vimentin and fetuin A isoforms. Protein J. 2009, 28 (3–4): 148-160.PubMedView ArticleGoogle Scholar
  31. Hamaguchi T, Iizuka N, Tsunedomi R, Hamamoto Y, Miyamoto T, Iida M, Tokuhisa Y, Sakamoto K, Takashima M, Tamesa T: Glycolysis module activated by hypoxia-inducible factor 1alpha is related to the aggressive phenotype of hepatocellular carcinoma. Int J Oncol. 2008, 33 (4): 725-731.PubMedGoogle Scholar
  32. Falanga V, Zhou L, Yufit T: Low oxygen tension stimulates collagen synthesis and COL1A1 transcription through the action of TGF-beta1. J Cell Physiol. 2002, 191 (1): 42-50. 10.1002/jcp.10065.PubMedView ArticleGoogle Scholar
  33. Tokunaga K, Nakamura Y, Sakata K, Fujimori K, Ohkubo M, Sawada K, Sakiyama S: Enhanced expression of a glyceraldehyde-3-phosphate dehydrogenase gene in human lung cancers. Cancer Res. 1987, 47 (21): 5616-5619.PubMedGoogle Scholar
  34. Graven KK, Farber HW: Hypoxia-associated proteins. New Horiz. 1995, 3 (2): 208-218.PubMedGoogle Scholar
  35. Li R, Wang H, Bekele BN, Yin Z, Caraway NP, Katz RL, Stass SA, Jiang F: Identification of putative oncogenes in lung adenocarcinoma by a comprehensive functional genomic approach. Oncogene. 2006, 25 (18): 2628-2635. 10.1038/sj.onc.1209289.PubMedView ArticleGoogle Scholar
  36. Kalousi A, Mylonis I, Politou AS, Chachami G, Paraskeva E, Simos G: Casein kinase 1 regulates human hypoxia-inducible factor HIF-1. J Cell Sci. 2010, 123 (Pt 17): 2976-2986.PubMedView ArticleGoogle Scholar
  37. Bae S-K, Kim S-R, Kim JG, Kim JY, Koo TH, Jang H-O, Yun I, Yoo M-A, Bae M-K: Hypoxic induction of human visfatin gene is directly mediated by hypoxia-inducible factor-1. FEBS Lett. 2006, 580: 4105-4113. 10.1016/j.febslet.2006.06.052.PubMedView ArticleGoogle Scholar
  38. Planque C, Kulasingam V, Smith CR, Reckamp K, Goodglick L, Diamandis EP: Identification of five candidate lung cancer biomarkers by proteomics analysis of conditioned media of four lung cancer cell lines. Mol Cell Proteomics. 2009, 8 (12): 2746-2758. 10.1074/mcp.M900134-MCP200.PubMedPubMed CentralView ArticleGoogle Scholar
  39. van Uden P, Kenneth NS, Webster R, Muller HA, Mudie S, Rocha S: Evolutionary conserved regulation of HIF-1beta by NF-kappaB. PLoS Genet. 2011, 7 (1): e1001285-10.1371/journal.pgen.1001285.PubMedPubMed CentralView ArticleGoogle Scholar
  40. Carrero P, Okamoto K, Coumailleau P, O’Brien S, Tanaka H, Poellinger L: Redox-regulated recruitment of the transcriptional coactivators CREB-binding protein and SRC-1 to hypoxia-inducible factor 1alpha. Mol Cell Biol. 2000, 20 (1): 402-415. 10.1128/MCB.20.1.402-415.2000.PubMedPubMed CentralView ArticleGoogle Scholar
  41. Salnikow K, Aprelikova O, Ivanov S, Tackett S, Kaczmarek M, Karaczyn A, Yee H, Kasprzak KS, Niederhuber J: Regulation of hypoxia-inducible genes by ETS1 transcription factor. Carcinogenesis. 2008, 29 (8): 1493-1499. 10.1093/carcin/bgn088.PubMedPubMed CentralView ArticleGoogle Scholar
  42. Wang F, Zhang R, Wu X, Hankinson O: Roles of coactivators in hypoxic induction of the erythropoietin gene. PLoS One. 2010, 5 (4): e10002-10.1371/journal.pone.0010002.PubMedPubMed CentralView ArticleGoogle Scholar
  43. Wang B, Hasan MK, Alvarado E, Yuan H, Wu H, Chen WY: NAMPT overexpression in prostate cancer and its contribution to tumor cell survival and stress response. Oncogene. 2011, 30 (8): 907-921. 10.1038/onc.2010.468.PubMedView ArticleGoogle Scholar
  44. Filiano AJ, Bailey CD, Tucholski J, Gundemir S, Johnson GV: Transglutaminase 2 protects against ischemic insult, interacts with HIF1beta, and attenuates HIF1 signaling. FASEB J. 2008, 22 (8): 2662-2675. 10.1096/fj.07-097709.PubMedPubMed CentralView ArticleGoogle Scholar
  45. Choi H, Chun YS, Kim TY, Park JW: HIF-2alpha enhances beta-catenin/TCF-driven transcription by interacting with beta-catenin. Cancer Res. 2010, 70 (24): 10101-10111. 10.1158/0008-5472.CAN-10-0505.PubMedView ArticleGoogle Scholar
  46. Arvidsson Y, Andersson E, Bergström A, Andersson MK, Altiparmak G, Illerskog A-C, Ahlman H, Lamazhapova D, Nilsson O: Amyloid precursor-like protein 1 is differentially upregulated in neuroendocrine tumours of the gastrointestinal tract. Endocrinerelated cancer. 2008, 15: 569-581.View ArticleGoogle Scholar
  47. Urtreger AJ, Werbajh SE, Verrecchia F, Mauviel A, Puricelli LI, Kornblihtt AR, Bal de Kier Joffe ED: Fibronectin is distinctly downregulated in murine mammary adenocarcinoma cells with high metastatic potential. Oncol Rep. 2006, 16 (6): 1403-1410.PubMedGoogle Scholar
  48. Xu Y, Shiraishi K, Mori M, Motomiya M: Changes of fibronectin in the right and left ventricles of rats exposed to chronic normobaric hypoxia. Tohoku J Exp Med. 1992, 168 (4): 573-582. 10.1620/tjem.168.573.PubMedView ArticleGoogle Scholar
  49. Tan M, Gu Q, He H, Pamarthy D, Semenza GL, Sun Y: SAG/ROC2/RBX2 is a HIF-1 target gene that promotes HIF-1 alpha ubiquitination and degradation. Oncogene. 2008, 27 (10): 1404-1411. 10.1038/sj.onc.1210780.PubMedView ArticleGoogle Scholar
  50. You J-J, Yang C-M, Chen M-S, Yang C-H: Regulation of Cyr61/CCN1 expression by hypoxia through cooperation of c-Jun/AP-1 and HIF-1α in retinal vascular endothelial cells. Exp Eye Res. 2010, 91: 825-836. 10.1016/j.exer.2010.10.006.PubMedView ArticleGoogle Scholar
  51. Olsen RR, Zetter BR: Evidence of a Role for Antizyme and Antizyme Inhibitor as Regulators of Human Cancer. Mol Cancer Res. 2011, 9 (10): 1285-1293. 10.1158/1541-7786.MCR-11-0178.PubMedPubMed CentralView ArticleGoogle Scholar
  52. Kim JS, Kim TL, Cho EW, Paik SG, Chung HW, Kim IG: Antizyme suppression leads to an increment of the cellular redox potential and an induction of HIF-1alpha: its involvement in resistance to gamma-radiation. IUBMB Life. 2008, 60 (6): 402-409. 10.1002/iub.49.PubMedView ArticleGoogle Scholar
  53. Nicol GR, Han M, Kim J, Birse CE, Brand E, Nguyen A, Mesri M, FitzHugh W, Kaminker P, Moore PA: Use of an immunoaffinity-mass spectrometry-based approach for the quantification of protein biomarkers from serum samples of lung cancer patients. Mol Cell Proteomics. 2008, 7 (10): 1974-1982. 10.1074/mcp.M700476-MCP200.PubMedView ArticleGoogle Scholar
  54. Oshiumi H, Sasai M, Shida K, Fujita T, Matsumoto M, Seya T: TIR-containing adapter molecule (TICAM)-2, a bridging adapter recruiting to toll-like receptor 4 TICAM-1 that induces interferon-beta. J Biol Chem. 2003, 278 (50): 49751-49762. 10.1074/jbc.M305820200.PubMedView ArticleGoogle Scholar
  55. Zhang JJ, Wu HS, Wang L, Tian Y, Zhang JH, Wu HL: Expression and significance of TLR4 and HIF-1alpha in pancreatic ductal adenocarcinoma. World J Gastroenterol. 2010, 16 (23): 2881-2888. 10.3748/wjg.v16.i23.2881.PubMedPubMed CentralView ArticleGoogle Scholar
  56. Warburton D, Olver BE: Coordination of genetic, epigenetic, and environmental factors in lung development, injury, and repair. Chest. 1997, 111 (6 Suppl): 119S-122S.PubMedView ArticleGoogle Scholar
  57. Benlhabib H, Mendelson CR: Epigenetic regulation of surfactant protein A gene (SP-A) expression in fetal lung reveals a critical role for Suv39h methyltransferases during development and hypoxia. Mol Cell Biol. 2011, 31 (10): 1949-1958. 10.1128/MCB.01063-10.PubMedPubMed CentralView ArticleGoogle Scholar
  58. Islam KN, Mendelson CR: Permissive effects of oxygen on cyclic AMP and interleukin-1 stimulation of surfactant protein A gene expression are mediated by epigenetic mechanisms. Mol Cell Biol. 2006, 26 (8): 2901-2912. 10.1128/MCB.26.8.2901-2912.2006.PubMedPubMed CentralView ArticleGoogle Scholar
  59. Scheurle D, DeYoung MP, Binninger DM, Page H, Jahanzeb M, Narayanan R: Cancer gene discovery using digital differential display. Cancer Res. 2000, 60 (15): 4037-4043.PubMedGoogle Scholar
  60. Dudley JT, Butte AJ: Identification of discriminating biomarkers for human disease using integrative network biology. Pac Symp Biocomput. 2009, 22: 27-38.Google Scholar
  61. Frohlich H: Network based consensus gene signatures for biomarker discovery in breast cancer. PLoS One. 2011, 6 (10): e25364-10.1371/journal.pone.0025364.PubMedPubMed CentralView ArticleGoogle Scholar
  62. Talks KL, Turley H, Gatter KC, Maxwell PH, Pugh CW, Ratcliffe PJ, Harris AL: The expression and distribution of the hypoxia-inducible factors HIF-1alpha and HIF-2alpha in normal human tissues, cancers, and tumor-associated macrophages. Am J Pathol. 2000, 157 (2): 411-421. 10.1016/S0002-9440(10)64554-3.PubMedPubMed CentralView ArticleGoogle Scholar
  63. Ibanez de Caceres I, Dulaimi E, Hoffman AM, Al-Saleem T, Uzzo RG, Cairns P: Identification of novel target genes by an epigenetic reactivation screen of renal cancer. Cancer Res. 2006, 66 (10): 5021-5028. 10.1158/0008-5472.CAN-05-3365.PubMedView ArticleGoogle Scholar
  64. Cheung W, Darfler MM, Alvarez H, Hood BL, Conrads TP, Habbe N, Krizman DB, Mollenhauer J, Feldmann G, Maitra A: Application of a global proteomic approach to archival precursor lesions: deleted in malignant brain tumors 1 and tissue transglutaminase 2 are upregulated in pancreatic cancer precursors. Pancreatology. 2008, 8 (6): 608-616. 10.1159/000161012.PubMedPubMed CentralView ArticleGoogle Scholar
  65. Sun D, Zhou M, Kowolik CM, Trisal V, Huang Q, Kernstine KH, Lian F, Shen B: Differential expression patterns of capping protein, protein phosphatase 1, and casein kinase 1 may serve as diagnostic markers for malignant melanoma. Melanoma Res. 2011, 21 (4): 335-343. 10.1097/CMR.0b013e328346b715.PubMedPubMed CentralView ArticleGoogle Scholar
  66. Menke A, Philippi C, Vogelmann R, Seidel B, Lutz MP, Adler G, Wedlich D: Down-regulation of E-cadherin gene expression by collagen type I and type III in pancreatic cancer cell lines. Cancer Res. 2001, 61 (8): 3508-3517.PubMedGoogle Scholar
  67. Bauer L, Venz S, Junker H, Brandt R, Radons J: Nicotinamide phosphoribosyltransferase and prostaglandin H2 synthase 2 are up-regulated in human pancreatic adenocarcinoma cells after stimulation with interleukin-1. Int J Oncol. 2009, 35 (1): 97-107.PubMedGoogle Scholar
  68. Chui YL, Ching AK, Chen S, Yip FP, Rowlands DK, James AE, Lee KK, Chan JY: BRE over-expression promotes growth of hepatocellular carcinoma. Biochem Biophys Res Commun. 2010, 391 (3): 1522-1525. 10.1016/j.bbrc.2009.12.111.PubMedView ArticleGoogle Scholar
  69. Qin L, Chen X, Wu Y, Feng Z, He T, Wang L, Liao L, Xu J: Steroid receptor coactivator-1 upregulates integrin alpha expression to promote breast cancer cell adhesion and migration. Cancer Res. 2011, 71 (5): 1742-1751. 10.1158/0008-5472.CAN-10-3453.PubMedPubMed CentralView ArticleGoogle Scholar
  70. Jinawath N, Vasoontara C, Jinawath A, Fang X, Zhao K, Yap KL, Guo T, Lee CS, Wang W, Balgley BM: Oncoproteomic analysis reveals co-upregulation of RELA and STAT5 in carboplatin resistant ovarian carcinoma. PLoS One. 2010, 5 (6): e11198-10.1371/journal.pone.0011198.PubMedPubMed CentralView ArticleGoogle Scholar
  71. Arvidsson Y, Andersson E, Bergstrom A, Andersson MK, Altiparmak G, Illerskog AC, Ahlman H, Lamazhapova D, Nilsson O: Amyloid precursor-like protein 1 is differentially upregulated in neuroendocrine tumours of the gastrointestinal tract. Endocr Relat Cancer. 2008, 15 (2): 569-581. 10.1677/ERC-07-0145.PubMedView ArticleGoogle Scholar
  72. Luo JH, Xie D, Liu MZ, Chen W, Liu YD, Wu GQ, Kung HF, Zeng YX, Guan XY: Protein expression and amplification of AIB1 in human urothelial carcinoma of the bladder and overexpression of AIB1 is a new independent prognostic marker of patient survival. Int J Cancer. 2008, 122 (11): 2554-2561. 10.1002/ijc.23399.PubMedView ArticleGoogle Scholar

Copyright

© Srivastava et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement