Gene expression signature predicts lymphatic metastasis in squamous cell carcinoma of the oral cavity

Tissue transglutaminase mediates the pro-malignant effects of oncostatin M receptor over-expression in cervical squamous cell carcinoma

Oncostatin M receptor (OSMR) is commonly over-expressed in advanced cervical squamous cell carcinoma (SCC), producing a significantly worse clinical outcome. Cervical SCC cells that over-express OSMR show enhanced responsiveness to the major ligand OSM, which induces multiple pro-malignant effects, including increased cell migration and invasiveness. Here, we show that tissue transglutaminase (TGM2) is an important mediator of the ligand-dependent phenotypic effects of OSMR over-expression in SCC cells. TGM2 expression correlated with disease progression and with OSMR levels in clinical samples of cervical and oral SCC. TGM2 depletion in cervical SCC cells abrogated OSM-induced migration on fibronectin-coated surfaces and invasiveness through extracellular matrix, while ectopic expression of TGM2 increased cell motility and invasiveness. Confocal microscopy and co-immunoprecipitation assays showed that TGM2 interacted with integrin-α5β1 in the presence of fibronectin in cervical SCC cells, with OSM treatment strengthening the interaction. Importantly, integrin-α5β1 and fibronectin were also over-expressed in cervical and oral SCC, where levels correlated with those of OSMR and TGM2. This combined tissue and in vitro study demonstrates for the first time that stimulation of over-expressed OSMR in cervical SCC cells activates TGM2/integrin-α5β1 interactions and induces pro-malignant changes. We conclude that an OSMR/TGM2/integrin-α5β1/fibronectin pathway is of biological significance in cervical SCC and a candidate for therapeutic targeting.

Upregulation of Src homology phosphotyrosyl phosphatase 2 (Shp2) expression in oral cancer and knockdown of Shp2 expression inhibit tumor cell viability and invasion in vitro

OBJECTIVE

This study investigated the clinical significance of Shp2 protein expression in oral squamous cell carcinoma (OSCC) and elucidated its biologic significance in OSCC cells.

STUDY DESIGN

A total of 88 OSCC cases were used to assess Shp2 expression, out of which 70 were for immunohistochemistry and 18 paired tumors vs normal tissues were for Western blot of Shp2 expression. OSCC cells were used to assess the effects of Shp2 knockdown for cell viability, apoptosis, invasion, and protein expressions.

RESULTS

Expression of Shp2 protein was significantly upregulated in OSCC tissues compared with the normal tissues, and Shp2 overexpression was associated with advanced tumor clinical stages and lymph node metastasis ex vivo. Knockdown of Shp2 expression in vitro inhibited OSCC cell viability and invasion but induced apoptosis by regulating expression of the apoptosis-related proteins.

CONCLUSIONS

The data indicated that Shp2 may play an important role in OSCC progression. Further studies will investigate whether a target of Shp2 expression could be a novel therapeutic strategy for clinical control of OSCC.

Unraveling regulatory programs for NF-kappaB, p53 and microRNAs in head and neck squamous cell carcinoma

In head and neck squamous cell carcinoma (HNSCC), mutations of p53 usually coexist with aberrant activation of NF-kappaB (NF-κB), other transcription factors and microRNAs, which promote tumor pathogenesis. However, how these factors and microRNAs interact to globally modulate gene expression and mediate oncogenesis is not fully understood. We devised a novel bioinformatics method to uncover interactive relationships between transcription factors or microRNAs and genes. This approach is based on matrix decomposition modeling under the joint constraints of sparseness and regulator-target connectivity, and able to integrate gene expression profiling and binding data of regulators. We employed this method to infer the gene regulatory networks in HNSCC. We found that the majority of the predicted p53 targets overlapped with those for NF-κB, suggesting that the two transcription factors exert a concerted modulation on regulatory programs in tumor cells. We further investigated the interrelationships of p53 and NF-κB with five additional transcription factors, AP1, CEBPB, EGR1, SP1 and STAT3, and microRNAs mir21 and mir34ac. The resulting gene networks indicate that interactions among NF-κB, p53, and the two miRNAs likely regulate progression of HNSCC. We experimentally validated our findings by determining expression of the predicted NF-κB and p53 target genes by siRNA knock down, and by examining p53 binding activity on promoters of predicted target genes in the tumor cell lines. Our results elucidating the cross-regulations among NF-κB, p53, and microRNAs provide insights into the complex regulatory mechanisms underlying HNSCC, and shows an efficient approach to inferring gene regulatory programs in biological complex systems.

Identification and validation of a multigene predictor of recurrence in primary laryngeal cancer

Purpose

Local recurrence is the major manifestation of treatment failure in patients with operable laryngeal carcinoma. Established clinicopathological factors cannot sufficiently predict patients that are likely to recur after treatment. Additional tools are therefore required to accurately identify patients at high risk for recurrence. This study attempts to identify and independently validate gene expression models, prognostic of disease-free survival (DFS) in operable laryngeal cancer.

Materials and Methods

Using Affymetrix U133A Genechips, we profiled fresh-frozen tumor tissues from 66 patients with laryngeal cancer treated locally with surgery. We applied Cox regression proportional hazards modeling to identify multigene predictors of recurrence. Gene models were then validated in two independent cohorts of 54 and 187 patients (fresh-frozen and formalin-fixed tissue validation sets, respectively).

Results

We focused on genes univariately associated with DFS (p<0.01) in the training set. Among several models comprising different numbers of genes, a 30-probe set model demonstrated optimal performance in both the training (log-rank, p<0.001) and 1^st validation (p = 0.010) sets. Specifically, in the 1^st validation set, median DFS as predicted by the 30-probe set model, was 34 and 80 months for high- and low-risk patients, respectively. Hazard ratio (HR) for recurrence in the high-risk group was 3.87 (95% CI 1.28–11.73, Wald's p = 0.017). Testing the expression of selected genes from the above model in the 2^nd validation set, with qPCR, revealed significant associations of single markers, such as ACE2, FLOT1 and PRKD1, with patient DFS. High PRKD1 remained an unfavorable prognostic marker upon multivariate analysis (HR = 2.00, 95% CI 1.28–3.14, p = 0.002) along with positive nodal status.

Conclusions

We have established and validated gene models that can successfully stratify patients with laryngeal cancer, based on their risk for recurrence. It seems worthy to prospectively validate PRKD1 expression as a laryngeal cancer prognostic marker, for routine clinical applications.

Signaling Networks of Activated Oncogenic and Altered Tumor Suppressor Genes in Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC) arises from the upper aerodigestive tract and is the six most common cancers worldwide. HNSCC is associated with high morbidity and mortality, as standard surgery, radiation, and chemotherapy can cause significant disfigurement and only provide 5-year survival rates of ~50-60%. The heterogeneity of HNSCC subsets with different potentials for recurrence and metastasis challenges the traditional pathological classification system, thereby increasing demand for the development of new diagnostic, prognostic, and therapeutic tools based on global molecular signatures of HNSCC. Historically, using classical biological techniques, it has been extremely difficult and time-consuming to survey hundreds or thousands of genes in a given disease. However, the development of high throughput technologies and high-powered computation throughout the last two decades has enabled us to investigate hundreds or thousands of genes simultaneously. Using high throughput technologies, our laboratory has identified the gene signatures and protein networks, which significantly affect HNSCC malignant phenotypes, including TP53/p63/p73 family members, IL-1/TNF-β/NF-κB, PI3K/AKT/mTOR, IL-6/IL-6R/JAK/STAT3, EGFR/MAPK/AP1, HGF/cMET/EGR1, and TGFβ/TGFβR/TAK1/SMAD pathways. This review summarizes the results from high-throughput technological assays conducted on HNSCC samples, including microarray, DNA methylation, miRNA profiling, and protein array, using primarily experimental data and conclusions generated in our own laboratory. The use of bioinformatics and integrated analyses of data sets from different platforms, as well as meta-analysis of large datasets pulled from multiple publicly available studies, provided significantly higher statistical power to extract biologically relevant information. The data suggested that the heterogeneity of HNSCC genotype and phenotype are much more complex than we previously thought. Understanding of global molecular signatures and disease classification for specific subsets of HNSCC will be essential to provide accurate diagnoses for targeted therapy and personalized treatment, which is an important effort toward improving patient outcomes.

EZH2 inhibition: targeting the crossroad of tumor invasion and angiogenesis

Tumor angiogenesis and metastatic spreading are two highly interconnected phenomena, which contribute to cancer-associated deaths. Thus, the identification of novel strategies to target angiogenesis and metastatic spreading is crucial. Polycomb genes are a set of epigenetic effectors, structured in multimeric repressive complexes. EZH2 is the catalytic subunit of Polycomb repressive complex 2 (PRC2), which methylates histone H3 lysine 27, thereby silencing several tumor-suppressor genes. EZH2 is essential for cancer stem cell self-renewal. Interestingly, cancer stem cells are thought to be the seeds of metastatic spreading and are able to differentiate into tumor-associated endothelial cells. Pre-clinical studies showed that EZH2 is able to silence several anti-metastatic genes (e.g., E-cadherin and tissue inhibitors of metalloproteinases), thereby favoring cell invasion and anchorage-independent growth. In addition, EZH2 seems to play a crucial role in the regulation of tumor angiogenesis. High EZH2 expression predicts poor prognosis, high grade, and high stage in several cancer types. Recently, a small molecule inhibitor of PRC2 (DZNeP) demonstrated promising anti-tumor activity, both in vitro and in vivo. Interestingly, DZNeP was able to inhibit cancer cell invasion and tumor angiogenesis in prostate and brain cancers, respectively. At tumor-inhibiting doses, DZNeP is not harmful for non-transformed cells. In the present manuscript, we review current evidence supporting a role of EZH2 in metastatic spreading and tumor angiogenesis. Using Oncomine datasets, we show that DZNeP targets are specifically silenced in some metastatic cancers, and some of them may inhibit angiogenesis. Based on this evidence, we propose the development of EZH2 inhibitors as anti-angiogenic and anti-metastatic therapy.

Network information improves cancer outcome prediction

Disease progression in cancer can vary substantially between patients. Yet, patients often receive the same treatment. Recently, there has been much work on predicting disease progression and patient outcome variables from gene expression in order to personalize treatment options. Despite first diagnostic kits in the market, there are open problems such as the choice of random gene signatures or noisy expression data. One approach to deal with these two problems employs protein-protein interaction networks and ranks genes using the random surfer model of Google's PageRank algorithm. In this work, we created a benchmark dataset collection comprising 25 cancer outcome prediction datasets from literature and systematically evaluated the use of networks and a PageRank derivative, NetRank, for signature identification. We show that the NetRank performs significantly better than classical methods such as fold change or t-test. Despite an order of magnitude difference in network size, a regulatory and protein-protein interaction network perform equally well. Experimental evaluation on cancer outcome prediction in all of the 25 underlying datasets suggests that the network-based methodology identifies highly overlapping signatures over all cancer types, in contrast to classical methods that fail to identify highly common gene sets across the same cancer types. Integration of network information into gene expression analysis allows the identification of more reliable and accurate biomarkers and provides a deeper understanding of processes occurring in cancer development and progression.

Validation of a gene expression signature for assessment of lymph node metastasis in oral squamous cell carcinoma

PURPOSE

Current assessment of lymph node metastasis in patients with head and neck squamous cell carcinoma is not accurate enough to prevent overtreatment. The aim of this study was validation of a gene expression signature for distinguishing metastasizing (N+) from nonmetastasizing (N0) squamous cell carcinoma of the oral cavity (OSCC) and oropharynx (OPSCC) in a large multicenter cohort, using a diagnostic DNA microarray in a Clinical Laboratory Improvement Amendments/International Organization for Standardization-approved laboratory.

METHODS

A multigene signature, previously reported as predictive for the presence of lymph node metastases in OSCC and OPSCC, was first re-evaluated and trained on 94 samples using generic, whole-genome, DNA microarrays. Signature genes were then transferred to a dedicated diagnostic microarray using the same technology platform. Additional samples (n=222) were collected from all head and neck oncologic centers in the Netherlands and analyzed with the diagnostic microarray. Human papillomavirus status was determined by real-time quantitative polymerase chain reaction.

RESULTS

The negative predictive value (NPV) of the diagnostic signature on the entire validation cohort (n=222) was 72%. The signature performed well on the most relevant subset of early-stage (cT1-T2N0) OSCC (n=101), with an NPV of 89%.

CONCLUSION

Combining current clinical assessment with the expression signature would decrease the rate of undetected nodal metastases from 28% to 11% in early-stage OSCC. This should be sufficient to enable clinicians to refrain from elective neck treatment. A new clinical decision model that incorporates the expression signature is therefore proposed for testing in a prospective study, which could substantially improve treatment for this group of patients.

Four-protein signature accurately predicts lymph node metastasis and survival in oral squamous cell carcinoma

The presence of lymph node (LN) metastasis significantly affects the survival of patients with oral squamous cell carcinoma (OSCC). Successful detection and removal of positive LNs are crucial in the treatment of this disease. Current evaluation methods still have their limitations in detecting the presence of tumor cells in the LNs, where up to a third of clinically diagnosed metastasis-negative (N0) patients actually have metastasis-positive LNs in the neck. We developed a molecular signature in the primary tumor that could predict LN metastasis in OSCC. A total of 211 cores from 55 individuals were included in the study. Eleven proteins were evaluated using immunohistochemical analysis in a tissue microarray. Of the 11 biomarkers evaluated using receiver operating curve analysis, epidermal growth factor receptor (EGFR), v-erb-b2 erythroblastic leukemia viral oncogene homolog 2 (HER-2/neu), laminin, gamma 2 (LAMC2), and ras homolog family member C (RHOC) were found to be significantly associated with the presence of LN metastasis. Unsupervised hierarchical clustering-demonstrated expression patterns of these 4 proteins could be used to differentiate specimens that have positive LN metastasis from those that are negative for LN metastasis. Collectively, EGFR, HER-2/neu, LAMC2, and RHOC have a specificity of 87.5% and a sensitivity of 70%, with a prognostic accuracy of 83.4% for LN metastasis. We also demonstrated that the LN signature could independently predict disease-specific survival (P = .036). The 4-protein LN signature validated in an independent set of samples strongly suggests that it could reliably distinguish patients with LN metastasis from those who were metastasis-free and therefore could be a prognostic tool for the management of patients with OSCC.

Prognostic biological features in neck dissection specimens

The superior prognostic value offered by routine histopathological staging of neck dissections, as compared to clinical staging using palpation and modern imaging techniques, is well established in the literature concerning the management of squamous cell carcinoma of the head and neck. In this review, we discuss the definitions and criteria used in standardised routine histopathological reporting and explore additional potential nodal prognostic features. In addition, we critically appraise the value of immunohistochemistry, histochemistry, molecular and other non-morphological techniques and suggest tumour and host features that merit further investigations.

A microRNA activity map of human mesenchymal tumors: connections to oncogenic pathways; an integrative transcriptomic study

Background

MicroRNAs (miRNAs) are nucleic acid regulators of many human mRNAs, and are associated with many tumorigenic processes. miRNA expression levels have been used in profiling studies, but some evidence suggests that expression levels do not fully capture miRNA regulatory activity. In this study we integrate multiple gene expression datasets to determine miRNA activity patterns associated with cancer phenotypes and oncogenic pathways in mesenchymal tumors – a very heterogeneous class of malignancies.

Results

Using a computational method, we identified differentially activated miRNAs between 77 normal tissue specimens and 135 sarcomas and we validated many of these findings with microarray interrogation of an independent, paraffin-based cohort of 18 tumors. We also showed that miRNA activity is imperfectly correlated with miRNA expression levels. Using next-generation miRNA sequencing we identified potential base sequence alterations which may explain differential activity. We then analyzed miRNA activity changes related to the RAS-pathway and found 21 miRNAs that switch from silenced to activated status in parallel with RAS activation. Importantly, nearly half of these 21 miRNAs were predicted to regulate integral parts of the miRNA processing machinery, and our gene expression analysis revealed significant reductions of these transcripts in RAS-active tumors. These results suggest an association between RAS signaling and miRNA processing in which miRNAs may attenuate their own biogenesis.

Conclusions

Our study represents the first gene expression-based investigation of miRNA regulatory activity in human sarcomas, and our findings indicate that miRNA activity patterns derived from integrated transcriptomic data are reproducible and biologically informative in cancer. We identified an association between RAS signaling and miRNA processing, and demonstrated sequence alterations as plausible causes for differential miRNA activity. Finally, our study highlights the value of systems level integrative miRNA/mRNA assessment with high-throughput genomic data, and the applicability of paraffin-tissue-derived RNA for validation of novel findings.

Possibility of the use of public microarray database for identifying significant genes associated with oral squamous cell carcinoma

There are lots of studies attempting to identify the expression changes in oral squamous cell carcinoma. Most studies include insufficient samples to apply statistical methods for detecting significant gene sets. This study combined two small microarray datasets from a public database and identified significant genes associated with the progress of oral squamous cell carcinoma. There were different expression scales between the two datasets, even though these datasets were generated under the same platforms - Affymetrix U133A gene chips. We discretized gene expressions of the two datasets by adjusting the differences between the datasets for detecting the more reliable information. From the combination of the two datasets, we detected 51 significant genes that were upregulated in oral squamous cell carcinoma. Most of them were published in previous studies as cancer-related genes. From these selected genes, significant genetic pathways associated with expression changes were identified. By combining several datasets from the public database, sufficient samples can be obtained for detecting reliable information. Most of the selected genes were known as cancer-related genes, including oral squamous cell carcinoma. Several unknown genes can be biologically evaluated in further studies.

Effect of siRNA-mediated downregulation of VEGF in Tca8113 cells on the activity of monocyte-derived dendritic cells

Vascular endothelial growth factor (VEGF) is a tumor angiogenesis factor that is important in immune regulation. In our previous study, we found that VEGF expression in the peripheral blood and neoplasm nest from patients with oral squamous cell carcinoma (OSCC) was positively correlated with the course of disease, while an inverse correlation between VEGF expression and dendritic cells (DCs) was identified in the peripheral blood. Therefore, in the present study, we investigated whether inhibition of human VEGF in the human tongue carcinoma cell line Tca8113 had effects on the activity of monocyte-derived DCs. We knocked down the expression of human VEGF in Tca8113 cells using the small interfering RNA (siRNA) technique. Tca8113 cells pre-transfected with siRNA targeting VEGF were co-cultured with monocyte‑derived immature and mature DCs. Cell proliferation was evaluated by a WST-8 assay. Cell apoptosis, cell cycle and cell phenotypes were determined by flow cytometry. The data revealed that downregulation of the human VEGF significantly inhibited the proliferation of Tca8113 cells and increased apoptosis. Inhibition of human VEGF arrested the cell cycle of Tca8113 cells at the G0/G1 phase. Our results showed that the co-culture of DCs with Tca8113 cells markedly inhibited the expression of the mature markers of DCs including HLA-DR, CD80, CD86, CD40 and CD1a, as well as the immature marker CD83, while inhibition of human VEGF in Tca8113 cells significantly reversed these effects. Therefore, human VEGF in Tca8113 cells may not only regulate the cell proliferation and apoptosis of oral squamous cell carcinoma cells, but may also inhibit DC maturation.

CD44 expression in oro-pharyngeal carcinoma tissues and cell lines

Expression of CD44, a transmembrane hyaluronan-binding glycoprotein, is variably considered to have prognostic significance for different cancers, including oral squamous cell carcinoma. Although unclear at present, tissue-specific expression of particular isoforms of CD44 might underlie the different outcomes in currently available studies. We mined public transcriptomics databases for gene expression data on CD44, and analyzed normal, immortalized and tumour-derived human cell lines for splice variants of CD44 at both the transcript and protein levels. Bioinformatics readouts, from a total of more than 15,000 analyses, implied an increased CD44 expression in head and neck cancer, including increased expression levels relative to many normal and tumor tissue types. Also, meta-analysis of over 260 cell lines and over 4,000 tissue specimens of diverse origins indicated lower CD44 expression levels in cell lines compared to tissue. With minor exceptions, reverse transcribed polymerase chain reaction identified expression of the four main isoforms of CD44 in normal oral keratinocytes, transformed lines termed DT and HaCaT, and a series of paired primary and metastasis-derived cell lines from oral or pharyngeal carcinomas termed HN4/HN12, HN22/HN8 and HN30/HN31. Immunocytochemistry, Western blotting and flow cytometric assessments all confirmed the isoform expression pattern at the protein level. Overall, bioinformatic processing of large numbers of global gene expression analyses demonstrated elevated CD44 expression in head and neck cancer relative to other cancer types, and that the application of standard cell culture protocols might decrease CD44 expression. Additionally, the results show that the many variant CD44 exons are not fundamentally deregulated in a diverse range of cultured normal and transformed keratinocyte lines.

Differentiation-promoting culture of competent and noncompetent keratinocytes identifies biomarkers for head and neck cancer

Aberrant contact-inhibited proliferation and differentiation induction couple with tumor severity, albeit with an imprecise association with prognosis. Assessment of contact inhibition and differentiation-promoting culture in this study of normal and immortalized oral keratinocytes (NOK and SVpgC2a, respectively) demonstrated elevated cloning ability and saturation density in the immortalized versus normal state, including consistent absence of differentiated morphological features. Transcriptomic analysis implicated 48 gene ontology categories, 8 molecular networks, and 10 key regulator genes in confluency-induced differentiation of NOK, all of which remained nonregulated in SVpgC2a. The SVpgC2a versus NOK transcriptome enriched 52 gene ontology categories altogether, 18 molecular networks, and 39 key regulator genes, several of which were associated with epithelial-mesenchymal transition. Assessment of the previously described gene sets relative to training data sets of head and neck squamous cell carcinoma samples, one including data on tumor differentiation and patient outcome and one present in the Human Gene Expression Map, identified four genes with association to poor survival (COX7A1, MFAP5, MPDU1, and POLD1). This gene set predicted poor outcome in an independent data set of 71 head and neck squamous cell carcinomas. The present study defines, for the first time to our knowledge, the broad gene spectrum that couples to induction, and loss, of oral keratinocyte differentiation. Bioinformatics assessments of the results relative to clinical data generated novel differentiation-related tumor biomarkers relevant to patient outcome.

Consistent Differential Expression Pattern (CDEP) on microarray to identify genes related to metastatic behavior

Background

To utilize the large volume of gene expression information generated from different microarray experiments, several meta-analysis techniques have been developed. Despite these efforts, there remain significant challenges to effectively increasing the statistical power and decreasing the Type I error rate while pooling the heterogeneous datasets from public resources. The objective of this study is to develop a novel meta-analysis approach, Consistent Differential Expression Pattern (CDEP), to identify genes with common differential expression patterns across different datasets.

Results

We combined False Discovery Rate (FDR) estimation and the non-parametric RankProd approach to estimate the Type I error rate in each microarray dataset of the meta-analysis. These Type I error rates from all datasets were then used to identify genes with common differential expression patterns. Our simulation study showed that CDEP achieved higher statistical power and maintained low Type I error rate when compared with two recently proposed meta-analysis approaches. We applied CDEP to analyze microarray data from different laboratories that compared transcription profiles between metastatic and primary cancer of different types. Many genes identified as differentially expressed consistently across different cancer types are in pathways related to metastatic behavior, such as ECM-receptor interaction, focal adhesion, and blood vessel development. We also identified novel genes such as AMIGO2, Gem, and CXCL11 that have not been shown to associate with, but may play roles in, metastasis.

Conclusions

CDEP is a flexible approach that borrows information from each dataset in a meta-analysis in order to identify genes being differentially expressed consistently. We have shown that CDEP can gain higher statistical power than other existing approaches under a variety of settings considered in the simulation study, suggesting its robustness and insensitivity to data variation commonly associated with microarray experiments.

Availability: CDEP is implemented in R and freely available at: http://genomebioinfo.musc.edu/CDEP/

Contact: zhengw@musc.edu

Identification of PRDX4 and P4HA2 as metastasis-associated proteins in oral cavity squamous cell carcinoma by comparative tissue proteomics of microdissected specimens using iTRAQ technology

Cervical lymph node metastasis represents the major prognosticator for oral cavity squamous cell carcinoma (OSCC). Here, we used an iTRAQ-based quantitative proteomic approach to identify proteins that are differentially expressed between microdissected primary and metastatic OSCC tumors. The selected candidates were examined in tissue sections via immunohistochemistry, and their roles in OSCC cell function investigated using RNA interference. Seventy-four differentially expressed proteins in nodal metastases, including PRDX4 and P4HA2, were identified. Immunohistochemical analysis revealed significantly higher levels of PRDX4 and P4HA2 in tumor cells than adjacent non-tumor epithelia (P < 0.0001 and P < 0.0001, respectively), and even higher expression in the 31 metastatic tumors of lymph nodes, compared to the corresponding primary tumors (P = 0.060 and P = 0.002, respectively). Overexpression of PRDX4 and P4HA2 was significantly associated with positive pN status (P = 0.048 and P = 0.021, respectively). PRDX4 overexpression was a significant prognostic factor for disease-specific survival in both univariate and multivariate analyses (P = 0.034 and P = 0.032, respectively). Additionally, cell migration and invasiveness were attenuated in OEC-M1 cells upon in vitro knockdown of PRDX4 and P4HA2 with specific interfering RNA. Novel metastasis-related prognostic markers for OSCC could be identified by our approach.

Discovering biological connections between experimental conditions based on common patterns of differential gene expression

Background

Identifying similarities between patterns of differential gene expression provides an opportunity to identify similarities between the experimental and biological conditions that give rise to these gene expression alterations. The growing volume of gene expression data in open data repositories such as the NCBI Gene Expression Omnibus (GEO) presents an opportunity to identify these gene expression similarities on a large scale across a diverse collection of datasets. We have developed a fast, pattern-based computational approach, named openSESAME (Search of Expression Signatures Across Many Experiments), that identifies datasets enriched in samples that display coordinate differential expression of a query signature. Importantly, openSESAME performs this search without prior knowledge of the phenotypic or experimental groups in the datasets being searched. This allows openSESAME to identify perturbations of gene expression that are due to phenotypic attributes that may not have been described in the sample annotation included in the repository.

To demonstrate the utility of openSESAME, we used gene expression signatures of two biological perturbations to query a set of 75,164 human expression profiles that were generated using Affymetrix microarrays and deposited in GEO. The first query, using a signature of estradiol treatment, identified experiments in which estrogen signaling was perturbed and also identified differences in estrogen signaling between estrogen receptor-positive and -negative breast cancers. The second query, which used a signature of silencing of the transcription factor p63 (a key regulator of epidermal differentiation), identified datasets related to stratified squamous epithelia or epidermal diseases such as melanoma.

Conclusions

openSESAME is a tool for leveraging the growing body of publicly available microarray data to discover relationships between different biological states based on common patterns of differential gene expression. These relationships may serve to generate hypotheses about the causes and consequences of specific patterns of observed differential gene expression. To encourage others to explore the utility of this approach, we have made a website for performing openSESAME queries freely available at http://opensesame.bu.edu.

Deregulation of Snai2 is associated with metastasis and poor prognosis in tongue squamous cell carcinoma

The members of the Snail superfamily of zinc-finger transcription factors, including Snai1 and Snai2, are involved in essential biological processes, such as epithelial-mesenchymal transition (EMT). Although Snai1 has been investigated in a number of cancers, our knowledge on Snai2 and its role(s) in squamous cell carcinoma of oral tongue (SCCOT) is limited. In this study, we confirmed the previous observation that over-expression of Snai2 is a frequent event in SCCOT. We further demonstrated that Snai2 over-expression is associated with lymph node metastasis in two independent SCCOT patient cohorts (total n = 129). Statistical analysis revealed that Snai2 over-expression was correlated with reduced overall survival. Furthermore, over-expression of Snai2 was correlated with reduced E-cadherin expression and enhanced Vimentin expression, suggesting a functional role of Snai2 in EMT. These observations were confirmed in vitro, in which knockdown of Snai2 induced a switch from a mesenchymal-like morphology to an epithelial-like morphology in SCCOT cell lines, and suppressed the cell invasion and migration. In contrast, ectopic transfection of Snai2 led to enhanced cell invasion and migration. Furthermore, Snai2 knockdown attenuated TGFβ1-induced EMT in SCCOT cell lines. Taken together, these data suggest that Snai2 plays major roles in EMT and the progression of SCCOT and may serve as a therapeutic target for patients at risk of metastasis.

Identification of guanylate-binding protein 1 as a potential oral cancer marker involved in cell invasion using omics-based analysis

Oral cavity squamous cell carcinoma (OSCC) is a devastating disease that accounts for 3% of all cancer cases diagnosed annually. OSCC is usually diagnosed at advanced clinical stages, resulting in poor outcomes. To identify effective biomarkers for improved OSCC diagnosis and/or management, we simultaneously analyzed the OSCC cell secretome and tissue transcriptome. Among the 19 candidates isolated, guanylate-binding protein 1 (GBP1) was selected for further validation using serum samples from OSCC patients and healthy controls. Notably, the serum level of GBP1 was higher in OSCC patients, compared to that in healthy controls. Immunohistochemical analysis further revealed GBP1 overexpression in OSCC tissues, compared with adjacent noncancerous epithelia. Importantly, the higher GBP1 level in OSCC tissue was associated with higher overall pathological stage, positive perineural invasion, and poorer prognosis. Moreover, GBP1 modulated the migration and invasion of OSCC cells in vitro. Our results collectively indicate that integrated analysis of the cancer secretome and transcriptome is a feasible strategy for the efficient identification of novel OSCC markers.

Meta-analysis of multiple microarray datasets reveals a common gene signature of metastasis in solid tumors

Background

Metastasis is the number one cause of cancer deaths. Expression microarrays have been widely used to study metastasis in various types of cancer. We hypothesize that a meta-analysis of publicly available gene expression datasets in various tumor types can identify a signature of metastasis that is common to multiple tumor types. This common signature of metastasis may help us to understand the shared steps in the metastatic process and identify useful biomarkers that could predict metastatic risk.

Methods

We identified 18 publicly available gene expression datasets in the Oncomine database comparing distant metastases to primary tumors in various solid tumors which met our eligibility criteria. We performed a meta-analysis using a modified permutation counting method in order to obtain a common gene signature of metastasis. We then validated this signature in independent datasets using gene set expression comparison analysis with the LS-statistic.

Results

A common metastatic signature of 79 genes was identified in the metastatic lesions compared with primaries with a False Discovery Proportion of less than 0.1. Interestingly, all the genes in the signature, except one, were significantly down-regulated, suggesting that overcoming metastatic suppression may be a key feature common to all metastatic tumors. Pathway analysis of the significant genes showed that the genes were involved in known metastasis-associated pathways, such as integrin signaling, calcium signaling, and VEGF signaling. To validate the signature, we used an additional six expression datasets that were not used in the discovery study. Our results showed that the signature was significantly enriched in four validation sets with p-values less than 0.05.

Conclusions

We have modified a previously published meta-analysis method and identified a common metastatic signature by comparing primary tumors versus metastases in various tumor types. This approach, as well as the gene signature identified, provides important insights to the common metastatic process and a foundation for future discoveries that could have broad application, such as drug discovery, metastasis prediction, and mechanistic studies.

Novel p63 target genes involved in paracrine signaling and keratinocyte differentiation

The transcription factor p63 is required for proper epidermal barrier formation and maintenance. Herein, we used chromatin immunoprecipitation coupled with DNA sequencing to identify novel p63 target genes involved in normal human epidermal keratinocyte (NHEKs) growth and differentiation. We identified over 2000 genomic sites bound by p63, of which 82 were also transcriptionally regulated by p63 in NHEKs. Through the discovery of interleukin-1-α as a p63 target gene, we identified that p63 is a regulator of epithelial–mesenchymal crosstalk. Further, three-dimensional organotypic co-cultures revealed TCF7L1, another novel p63 target gene, as a regulator of epidermal proliferation and differentiation, providing a mechanism by which p63 maintains the proliferative potential of basal epidermal cells. The discovery of new target genes links p63 to diverse signaling pathways required for epidermal development, including regulation of paracrine signaling to proliferative potential. Further mechanistic insight into p63 regulation of epidermal cell growth and differentiation is provided by the identification of a number of novel p63 target genes in this study.

KiSS1 mediates platinum sensitivity and metastasis suppression in head and neck squamous cell carcinoma

Although surgery and radiotherapy have been the standard treatment modalities for head and neck squamous cell carcinoma (HNSCC), the integration of cisplatin (CDDP)-based therapy has led to improvements in local and regional control of disease for patients. However, many trials show that only 10-20% of patients benefit from this treatment intensification, which can result in profound treatment-associated morbidity and mortality. Moreover, the marginal survival improvement suggests that CDDP resistance is an innate characteristic of HNSCC. To elucidate the biological mechanisms underpinning CDDP resistance in HNSCC, we utilized an experimental model of CDDP resistance in this disease. We first observed significant enhancements in local tumor growth and metastasis, as well as adverse survival, in CDDP-resistant (CR) tumors compared with sensitive tumors. To elucidate the molecular mechanisms of this phenotype, we undertook a systems biology-based approach utilizing high-throughput PCR arrays, and we identified a significant suppression of KiSS1 mRNA and protein expression in the CR cells, but no significant regions of genomic loss with array comparative genomic hybridization. Genetic suppression of KiSS1 in CDDP-sensitive cell lines rendered them CR, an observation that was mechanistically linked to alterations in glutathione S-transferase-π expression and function. We next confirmed that, in human HNSCC tumors, loss of KiSS1 expression was associated with metastatic human HNSCC tumors compared with non-metastatic tumors. Genetic reconstitution of KiSS1 in CR cells abrogated cellular migration and induced CDDP sensitivity. To confirm these findings in a murine model, either CR or KiSS1-transfected CR cells were studied in an orthotopic model of HNSCC, or survival studies revealed significant improvement in survival of the mice bearing CR-KiSS1 tumors. Mechanistically, alterations in apoptotic pathways and CDDP metabolism contributed to KiSS1-associated chemotherapy sensitization. These studies provided further direct evidence for the role of KiSS1 loss in biologically aggressive HNSCC and suggest potential targets for therapy in CR cancers.

The effect of a DNA repair gene on cellular invasiveness: XRCC3 over-expression in breast cancer cells

Over-expression of DNA repair genes has been associated with resistance to radiation and DNA-damage induced by chemotherapeutic agents such as cisplatin. More recently, based on the analysis of genome expression profiling, it was proposed that over-expression of DNA repair genes enhances the invasive behaviour of tumour cells. In this study we present experimental evidence utilizing functional assays to test this hypothesis. We assessed the effect of the DNA repair proteins known as X-ray complementing protein 3 (XRCC3) and RAD51, to the invasive behavior of the MCF-7 luminal epithelial-like and BT20 basal-like triple negative human breast cancer cell lines. We report that stable or transient over-expression of XRCC3 but not RAD51 increased invasiveness in both cell lines in vitro. Moreover, XRCC3 over-expressing MCF-7 cells also showed a higher tumorigenesis in vivo and this phenotype was associated with increased activity of the metalloproteinase MMP-9 and the expression of known modulators of cell-cell adhesion and metastasis such as CD44, ID-1, DDR1 and TFF1. Our results suggest that in addition to its' role in facilitating repair of DNA damage, XRCC3 affects invasiveness of breast cancer cell lines and the expression of genes associated with cell adhesion and invasion.

Is gene expression profiling of head and neck cancers ready for the clinic?

One of the major challenges in the head and neck oncology clinic is the need to identify biomarkers and/or gene expression signatures that complement, strengthen and increase the sensitivity and specificity of the current clinicopathologic analyses. Microarray analysis of head and neck tumors has demonstrated that the combined influence of many genes or biomarkers can make superior identifiers and/or predictors of tumor behavior and patient outcome. Here, an update of the recent literature on the prognostic and predictive value of microarrays for patients with head and neck squamous cell carcinomas is presented. Microarray technology has the potential for improved decision-making and corroboration within the clinical setting. However, further integration, standardization, validation and research are required before the use of microarray analysis is ready for routine clinical management of head and neck squamous cell carcinomas.

Gene Expression Signatures of Lymph Node Metastasis in Oral Cancer: Molecular Characteristics and Clinical Significances

Even though lymph node metastasis accounts for the vast majority of cancer death in patients with oral cancer (OC), the molecular mechanisms of lymph node metastasis remain elusive. Genome-wide microarray analyses and functional studies in vitro and in vivo, along with detailed clinical observations, have identified a number of molecules that may contribute to lymph node metastasis. These include lymphangionenic cytokines, cell adhesion molecules, basement membrane-interacting molecules, matrix enzymes and relevant downstream signaling pathways. However, defined gene signatures from different studies are highly variable, which hinders their translation to clinically relevant applications. To date, none of the identified signatures or molecular biomarkers has been successfully implemented as a diagnostic or prognostic tool applicable to routine clinical practice. In this review, we will first introduce the significance of lymph node metastasis in OC, and clinical/experimental evidences that support the underlying molecular mechanisms. We will then provide a comprehensive review and integrative analysis of the existing gene expression studies that aim to identify the metastasis-related signatures in OC. Finally, the remaining challenges will be discussed and our insights on future directions will be provided.

Current potential and limitations of molecular diagnostic methods in head and neck cancer

Traditional diagnostic methods such as clinical assessment, histopathological examination and imaging techniques are limited in their capacity to provide information on prognosis and treatment choice of head and neck cancer. In recent years, molecular techniques have been developed that enabled us to get more insight into the molecular biological cellular pathways underlying tumor progression and metastasis. Correlation of these molecular changes with clinical events has been explored. However, consistently useful markers have not been identified yet, although many promising developments are in progress. It may be expected that in the near future, molecular markers will be useful for clinical purposes. In this paper, an overview will be given of the several molecular techniques that may have potential to be introduced in clinical practice in the management of head and neck squamous cell carcinoma.

Deregulation of manganese superoxide dismutase (SOD2) expression and lymph node metastasis in tongue squamous cell carcinoma

Background

Lymph node metastasis is a critical event in the progression of tongue squamous cell carcinoma (TSCC). The identification of biomarkers associated with the metastatic process would provide critical prognostic information to facilitate clinical decision making. Previous studies showed that deregulation of manganese superoxide dismutase (SOD2) expression is a frequent event in TSCC and may be associated with enhanced cell invasion. The purpose of this study is to further evaluate whether the expression level of SOD2 is correlated with the metastatic status in TSCC patients.

Methods

We first examined the SOD2 expression at mRNA level on 53 TSCC and 22 normal control samples based on pooled-analysis of existing microarray datasets. To confirm our observations, we examined the expression of SOD2 at protein level on an additional TSCC patient cohort (n = 100), as well as 31 premalignant dysplasias, 15 normal tongue mucosa, and 32 lymph node metastatic diseases by immunohistochemistry (IHC).

Results

The SOD2 mRNA level in primary TSCC tissue is reversely correlated with lymph node metastasis in the first TSCC patient cohort. The SOD2 protein level in primary TSCC tissue is also reversely correlated with lymph node metastasis in the second TSCC patient cohort. Deregulation of SOD2 expression is a common event in TSCC and appears to be associated with disease progression. Statistical analysis revealed that the reduced SOD2 expression in primary tumor tissue is associated with lymph node metastasis in both TSCC patient cohorts examined.

Conclusions

Our study suggested that the deregulation of SOD2 in TSCC has potential predictive values for lymph node metastasis, and may serve as a therapeutic target for patients at risk of metastasis.

Hypoxic conversion of SMAD7 function from an inhibitor into a promoter of cell invasion

Smad7 is an inhibitor of the transforming growth factor-beta-activated signaling pathway. Under well-oxygenated conditions, Smad7 is a potent inhibitor of carcinoma cell invasion. Paradoxically, however, the expression of Smad7 is upregulated across several cancers and may promote cancer progression. Hypoxia, which is frequently met in solid tumors, is an enhancer of carcinoma cell invasion and cancer progression. Here, we report that hypoxia activates the expression of Smad7 in a hypoxia-inducible factor- and von Hippel-Lindau protein-dependent manner. As expected, in normoxia, the forced expression of Smad7 inhibited carcinoma cell invasion. In contrast with the normoxic condition, the inhibitory effect of Smad7 was lost under hypoxia. The block in carcinoma cell invasion by forced expression of Smad7 was released by hypoxia in two invasive carcinoma cell lines. Moreover, the noninvasive HaCaT keratinocytes become invasive upon simultaneous hypoxia and transforming growth factor-beta stimulus. The hypoxia-activated invasion was attenuated by inhibiting Smad7 expression by short interfering RNA. Finally, the increased Smad7 expression in human carcinomas correlated with hypoxic gene expression. The data provide evidence that hypoxia could convert Smad7 function from an invasion inhibitor into an activator of invasion. Furthermore, they might shed light as to why increased Smad7 expression is detected in cancers.

Altered phenotype of lymphatic endothelial cells induced by highly metastatic OTSCC cells contributed to the lymphatic metastasis of OTSCC cells

The lymphatic endothelial cell (LEC) is an interactive surface for cancer cells. This article aims to explore cancer cell-induced changes of LEC, and study the tumor-lymphatic endothelium interaction. Here, LECs were co-cultured with highly and poorly metastatic tongue cancer cells. The differences in biologic behaviors and gene expression profiles between them were examined. The results showed that LECs induced by highly metastatic cancer cells displayed abnormal biologic behaviors, and could secrete chemokines to promote the migration of cancer cells. Therefore, biologic properties and functional status of LECs in oral tongue squamous cell carcinoma (OTSCC) might be a positive factor in lymphatic dissemination.

TrkB induces EMT and has a key role in invasion of head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a significant public health problem, accounting for over 5% of all cancer-related deaths, and these deaths primarily result from metastatic disease. The molecular processes involved in HNSCC pathogenesis and progression are poorly understood, and here we present experimental evidence for a direct role of the cell surface receptor tyrosine kinase, TrkB, in HNSCC tumor progression. Using immunohistochemical analysis and transcriptional profiling of archival HNSCC tumor specimens, we found that TrkB and its secreted ligand, brain-derived neurotrophic factor (BDNF), are expresses in greater than 50% of human HNSCC tumors, but not in normal upper aerodigestive tract (UADT) epithelia. Studies with HNSCC cell lines reveal that in vitro stimulation with BDNF, the ligand for TrkB, upregulates the migration and invasion of HNSCC cells, and both transient and stable suppressions of TrkB result in significant abrogation of constitutive and ligand-mediated migration and invasion. Furthermore, enforced overexpression of TrkB results in altered expression of molecular mediators of epithelial-to-mesenchymal transition (EMT), including downregulation of E-cadherin and upregulation of Twist. Using an in vivo mouse model of HNSCC, we were able to show that downregulation of TrkB suppresses tumor growth. These results directly implicate TrkB in EMT and the invasive behavior of HNSCC, and correlate with the in vivo overexpression of TrkB in human HNSCC. Taken together, these data suggest that the TrkB receptor may be a critical component in the multi-step tumor progression of HNSCC, and may be an attractive target for much needed new therapies for this disease.

Fatty-acid-binding protein 5 promotes cell proliferation and invasion in oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is one of the most malignant neoplasms worldwide, and the molecular mechanism of oral tumorigenesis is still unclear. Fatty-acid-binding protein 5 (FABP5) was found in our previous study to be upregulated in oral squamous cell carcinomas by proteomic analysis. The implications of FABP5 overexpression in oral cancer progression have not yet been elucidated.

MATERIALS AND METHODS

In this study, the recombinant adeno-associated virus vectors were used to deliver and increase the expression of FABP5 in human OSCC cell lines. U6 promoter-driven short-hairpin RNA (shRNA)-triggered RNA interference was used to block FABP5 gene expression. Transwell Matrigel invasion assay, MTS cell proliferation assay, reverse transcription-polymerase chain reaction, Western blot, and gelatin zymography analysis were used to investigate the effects of FABP5 on cell invasion, growth, and matrix metalloproteinase (MMP) production.

RESULTS

Overexpression of FABP5 in oral cancer cells increased cell proliferation and invasiveness by increasing the expression of MMP-9. Silencing FABP5 with shRNA significantly suppressed cell proliferation, MMP-9 activities, and invasiveness.

CONCLUSION

Our study provides the first evidence that FABP5 expression modulated MMP-9 production and the invasive behavior of oral cancer cells and suggests that FABP5 may provide novel targets for therapeutic intervention.

Oral cancer cells with different potential of lymphatic metastasis displayed distinct biologic behaviors and gene expression profiles

OBJECTIVE

Oral squamous cell carcinoma (OSCC) often spreads from the primary tumor to regional lymph nodes in the early stage. Better understanding of the biology of lymphatic spread of oral cancer cells is important for improving the survival rate of cancer patients.

METHODS

We established the cell line LNMTca8113 by repeated injections in foot pads of nude mice, which had a much higher lymphatic metastasis rate than its parental cell line Tca8113. Then, we compared the biologic behaviors of cancer cells between them. Moreover, microarray-based expression profiles between them were also compared, and a panel of differential genes was validated using real-time-PCR.

RESULTS

In contrast to Tca8113 cells, LNMTca8113 cells were more proliferative and resistant to apoptosis in the absence of serum, and had enhanced ability of inducing capillary-like structures. Moreover, microarray-based expression profiles between them identified 1341 genes involved in cell cycle, cell adhesion, lymphangiogenesis, regulation of apoptosis, and so on. Some genes dedicating to the metastatic potential, including JAM2, TNC, CTSC, LAMB1, VEGFC, HAPLN1, ACPP, GDF9 and FGF11, were upregulated in LNMTca8113 cells.

CONCLUSION

These results suggested that LNMTca8113 and Tca8113 cells were proper models for lymphatic metastasis study because there were differences in biologic behaviors and metastasis-related genes between them. Additionally, the differentially expressed gene profiles in cancer progression may be helpful in exploring therapeutic targets and provide the foundation for further functional validation of these specific candidate genes for OSCC.

Prediction of disease severity in patients with early rheumatoid arthritis by gene expression profiling

In order to test the ability of peripheral blood gene expression profiles to predict future disease severity in patients with early rheumatoid arthritis (RA), a group of 17 patients (1 ± 0.2 years disease duration) was evaluated at baseline for gene expression profiles. Disease status was evaluated after a mean of 5 years using an index combining pain, global and recoded MHAQ scores. Unsupervised and supervised algorithms identified "predictor genes" whose combined expression levels correlated with follow-up disease severity scores. Unsupervised clustering algorithms separated patients into two branches. The only significant difference between these two groups was the disease severity score; demographic variables and medication usage were not different. Supervised T-Test analysis identified 19 "predictor genes" of future disease severity. Results were validated in an independent cohort of subjects of established RA with using Support Vector Machines and K-Nearest-Neighbor Classification. Our study demonstrates that peripheral blood gene expression profiles may be a useful tool to predict future disease severity in patients with early and established RA.

Exploring phenotype-associated modules in an oral cavity tumor using an integrated framework

MOTIVATION

Like most human diseases, tumors are complex traits, the genesis and development of which recruit a number of genes and several important biological processes. As proteins involved in common processes tend to be centralized in the same local area of protein-protein interaction networks, here a novel framework has been developed to identify which areas of the networks are most relevant to a phenotype.

RESULTS

These areas termed 'coherent modules' can be regarded as gene sets dynamically defined in the networks. Compared with previous analogous approaches, one critical feature of our method is the optimization of coherent modules for two distinct aspects balanced by tuning a parameter in the framework. First, we seek the low coupling between coherent modules and then maximize the intrinsic similarity within a module. The framework has good expansibility, with classical expression data analysis methods generalized as particular cases. This coherent module approach was applied to an oral cavity tumor dataset with 18 significant coherent modules identified, which could indicate the presence of lymph node metastasis. Further examination shows that most of the modules are responsible for comparatively independent biological processes. Our framework is helpful for the prognosis of tumors and offers a new perspective for tumor research.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of oral cancer cell-induced changes in gene expression profile of lymphatic endothelial cell

Understading functional properties of tumor-derived lymphatic endothelial cells (TLEC) are relevant for blocking lymphatic metastasis. The changes of lymphatic endothelial cells (LEC) cocultured with oral cancer cells in a vitro model were examined. TLEC, in contrast to LEC, were more proliferative and have enhanced ability of lymphangiogenesis and anti-apoptosis. Gene microarrays revealed that 677 unique genes had two-fold or higher change between the two groups. Differential expressions of selected genes were confirmed by real-time PCR. Our results indicate that TLEC display abnormal characteristics and are distinct at the molecular level. Manipulation of TLEC is encouraging for therapy of lymphatic metastasis.

Oral tongue cancer gene expression profiling: Identification of novel potential prognosticators by oligonucleotide microarray analysis

Background

The present study is aimed at identifying potential candidate genes as prognostic markers in human oral tongue squamous cell carcinoma (SCC) by large scale gene expression profiling.

Methods

The gene expression profile of patients (n=37) with oral tongue SCC were analyzed using Affymetrix HG_U95Av2 high-density oligonucleotide arrays. Patients (n=20) from which there were available tumor and matched normal mucosa were grouped into stage (early vs. late) and nodal disease (node positive vs. node negative) subgroups and genes differentially expressed in tumor vs. normal and between the subgroups were identified. Three genes, GLUT3, HSAL2, and PACE4, were selected for their potential biological significance in a larger cohort of 49 patients via quantitative real-time RT-PCR.

Results

Hierarchical clustering analyses failed to show significant segregation of patients. In patients (n=20) with available tumor and matched normal mucosa, 77 genes were found to be differentially expressed (P< 0.05) in the tongue tumor samples compared to their matched normal controls. Among the 45 over-expressed genes, MMP-1 encoding interstitial collagenase showed the highest level of increase (average: 34.18 folds). Using the criterion of two-fold or greater as overexpression, 30.6%, 24.5% and 26.5% of patients showed high levels of GLUT3, HSAL2 and PACE4, respectively. Univariate analyses demonstrated that GLUT3 over-expression correlated with depth of invasion (P<0.0001), tumor size (P=0.024), pathological stage (P=0.009) and recurrence (P=0.038). HSAL2 was positively associated with depth of invasion (P=0.015) and advanced T stage (P=0.047). In survival studies, only GLUT3 showed a prognostic value with disease-free (P=0.049), relapse-free (P=0.002) and overall survival (P=0.003). PACE4 mRNA expression failed to show correlation with any of the relevant parameters.

Conclusion

The characterization of genes identified to be significant predictors of prognosis by oligonucleotide microarray and further validation by real-time RT-PCR offers a powerful strategy for identification of novel targets for prognostication and treatment of oral tongue carcinoma.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes.

METHODS

Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries.

RESULTS

Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues.

CONCLUSION

To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.