Biomarkers of human skin cells identified using DermArray DNA arrays and new bioinformatics methods

Biomarker genes of human skin-derived cells were identified by new simple bioinformatic methods and DNA microarray analysis utilizing in vitro cultures of normal neonatal human epidermal keratinocytes, melanocytes, and dermal fibroblasts. A survey of 4405 human cDNAs was performed using DermArray DNA microarrays. Biomarkers were rank ordered by "likelihood ratio" algorithms and stringent selection criteria that have general applicability for analyzing a minimum of three RNA samples. Signature biomarker genes (up-regulated in one cell type) and anti-signature biomarker genes (down-regulated in one cell type) were determined for the three major skin cell types. Many of the signature genes are known biomarkers for these cell types. In addition, 17 signature genes were identified as ESTs, and 22 anti-signature biomarkers were discovered. Quantitative RT-PCR was used to verify nine signature biomarker genes. A total of 158 biomarkers of normal human skin cells were identified, many of which may be valuable in diagnostic applications and as molecular targets for drug discovery and therapeutic intervention.     

Keratinocyte growth factor induces gene expression signature associated with suppression of malignant phenotype of cutaneous squamous carcinoma cells

Keratinocyte growth factor (KGF, fibroblast growth factor-7) is a fibroblast-derived mitogen, which stimulates proliferation of epithelial cells. The expression of KGF by dermal fibroblasts is induced following injury and it promotes wound repair. However, the role of KGF in cutaneous carcinogenesis and cancer progression is not known. We have examined the role of KGF in progression of squamous cell carcinoma (SCC) of the skin. The expression of KGF receptor (KGFR) mRNA was lower in cutaneous SCCs (n = 6) than in normal skin samples (n = 6). Expression of KGFR mRNA was detected in 6 out of 8 cutaneous SCC cell lines and the levels were downregulated by 24-h treatment with KGF. KGF did not stimulate SCC cell proliferation, but it reduced invasion of SCC cells through collagen. Gene expression profiling of three cutaneous SCC cell lines treated with KGF for 24 h revealed a specific gene expression signature characterized by upregulation of a set of genes specifically downregulated in SCC cells compared to normal epidermal keratinocytes, including genes with tumor suppressing properties (SPRY4, DUSP4, DUSP6, LRIG1, PHLDA1). KGF also induced downregulation of a set of genes specifically upregulated in SCC cells compared to normal keratinocytes, including genes associated with tumor progression (MMP13, MATN2, CXCL10, and IGFBP3). Downregulation of MMP-13 and KGFR expression in SCC cells and HaCaT cells was mediated via ERK1/2. Activation of ERK1/2 in HaCaT cells and tumorigenic Ha-ras-transformed HaCaT cells resulted in downregulation of MMP-13 and KGFR expression. These results provide evidence, that KGF does not promote progression of cutaneous SCC, but rather suppresses the malignant phenotype of cutaneous SCC cells by regulating the expression of several genes differentially expressed in SCC cells, as compared to normal keratinocytes.     

DNA microarrays and likelihood ratio bioinformatic methods: discovery of human melanocyte biomarkers

In this article, some of the advantages and limitations of DNA microarray technologies for gene expression profiling are summarized. As a model experiment, DermArray DNA microarrays were utilized to identify potential biomarkers of cultured normal human melanocytes in two different experimental comparisons. In the first case, melanocyte RNA was compared with vastly dissimilar non-melanocytic RNA samples of normal skin keratinocytes and fibroblasts. In the second case, melanocyte RNA was compared with a primary cutaneous melanoma line (MS7) and a metastatic melanoma cell line (SKMel-28). The alternative approaches provide dramatically different lists of 'normal melanocyte' biomarkers. The most robust biomarkers were identified using principal component analysis bioinformatic methods related to likelihood ratios. Only three of 25 robust biomarkers in the melanocyte-proximal study (i.e. melanocytes vs. melanoma cells) were coincidentally identified in the melanocyte-distal study (i.e. melanocytes vs. non-melanocytic cells). Selected up-regulated biomarkers of melanocytes (i.e. TRP-1, melan-A/MART-1, silver/Pmel17, and nidogen-2) were validated by qRT-PCR. Some of the melanocytic biomarkers identified here may be useful in molecular diagnostics, as potential molecular targets for drug discovery, and for understanding the biochemistry of melanocytic cells.     

A method to improve selection of molecular targets by circumventing the ADME pharmacokinetic system utilizing PharmArray DNA microarrays

DNA microarrays may be used to identify potential molecular targets for drug discovery. Yet, DNA microarray experiments provide massive amounts of data. To limit the choice of potential molecular targets, it may be desirable to eliminate genes coincidentally up-regulated in tissues implicated in absorption, distribution, metabolism, and excretion (ADME) pharmacokinetics. DNA microarray experiments were performed to demonstrate a gene-exclusion approach using as an example RNA samples of neural origin, i.e., a human neuroblastoma cell line (SK-N-SH) and brain tissue, as the intended hypothetical site(s) of drug action. Biomarkers were identified using PharmArray DNA microarrays. The lists of neuroblastoma and neural biomarkers were constrained by limiting selection to the subset of genes that were not highly expressed in three transformed cell lines from liver, colon, and kidney (HepG2, Caco-2, and 786-O, respectively) that are routinely used as representatives of the ADME system during in vitro pharmacology and toxicology experiments. Principal component analysis methods with likelihood ratio-related bioinformatic tools were utilized to identify robust potential biomarker genes for the three ADME-related cell lines, neuroblastoma, and normal brain. Biomarkers of each sample were identified and selected genes were validated by qRT-PCR. Hundreds of biomarkers of the three ADME-related cell types, representing hepatocytes, kidney epithelium, and gastrointestinal tract, may now be used as a valuable database to restrict selection of biomarkers as potential molecular targets from the intended samples (e.g., neuroblastoma in this work). In addition to biomarker discovery per se, this demonstration suggests that our model method may be viable to help restrict gene lists during selection of potential molecular targets for subsequent drug discovery.     

BRAF and RKIP are significantly decreased in cutaneous squamous cell carcinoma

Background. Actinic keratosis (AK) is a well-established pre-cancerous skin lesion that has the potential to progress to squamous cell carcinoma (SCC). However, little is known about the implication of BRAF and RKIP expression, or about the incidence of BRAF mutations in the formation of these cutaneous diseases. The RAS oncogene has been proposed to significantly contribute to skin cancer development. Moreover, numerous BRAF mutations have been detected in melanoma biopsy specimens and cell lines.

Objectives. This study aimed to measure the mRNA levels of the genes BRAF and RKIP in AK, as well as their possible implication in the progress of AK to SCC. All biopsy specimens were also screened for BRAF mutations within exons 11 and 15.

Patients and methods. Expression levels of the genes BRAF and RKIP were examined in 16 AKs and 12 SCCs by RT-qPCR. A novel allele-specific qPCR method, in combination with direct DNA sequencing, was performed in order to inspect the frequency of the V600E mutation in exon 15, as well as to examine the mutation status of the gene within exon 11.

Results. Significant down-regulation was noted for both genes in SCC, compared to normal tissue (for BRAF, P=0.002; for RKIP, P     

Identification of a potential ovarian cancer stem cell gene expression profile from advanced stage papillary serous ovarian cancer

Identification of gene expression profiles of cancer stem cells may have significant implications in the understanding of tumor biology and for the design of novel treatments targeted toward these cells. Here we report a potential ovarian cancer stem cell gene expression profile from isolated side population of fresh ascites obtained from women with high-grade advanced stage papillary serous ovarian adenocarcinoma. Affymetrix U133 Plus 2.0 microarrays were used to interrogate the differentially expressed genes between side population (SP) and main population (MP), and the results were analyzed by paired T-test using BRB-ArrayTools. We identified 138 up-regulated and 302 down-regulated genes that were differentially expressed between all 10 SP/MP pairs. Microarray data was validated using qRT-PCR and17/19 (89.5%) genes showed robust correlations between microarray and qRT-PCR expression data. The Pathway Studio analysis identified several genes involved in cell survival, differentiation, proliferation, and apoptosis which are unique to SP cells and a mechanism for the activation of Notch signaling is identified. To validate these findings, we have identified and isolated SP cells enriched for cancer stem cells from human ovarian cancer cell lines. The SP populations were having a higher colony forming efficiency in comparison to its MP counterpart and also capable of sustained expansion and differentiation in to SP and MP phenotypes. 50,000 SP cells produced tumor in nude mice whereas the same number of MP cells failed to give any tumor at 8 weeks after injection. The SP cells demonstrated a dose dependent sensitivity to specific γ-secretase inhibitors implicating the role of Notch signaling pathway in SP cell survival. Further the generated SP gene list was found to be enriched in recurrent ovarian cancer tumors.     

Kallikrein-5 promotes cleavage of desmoglein-1 and loss of cell-cell cohesion in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) ranks among the top 8 causes of cancer death worldwide, with only a 60% 5-year survival rate, highlighting the need for discovery of novel biomarkers and therapeutic targets. We have previously reported that expression of a panel of serine proteinase kallikreins (KLK 5, 7, 8, and 10) is correlated with formation of more aggressive OSCC tumors in a murine orthotopic OSCC model and is elevated in human OSCC. Current studies focus on understanding the potential role of KLK5 in OSCC progression. In initial studies, KLK levels in malignant OSCC cells (SCC25) were compared with cells from normal oral mucosa (OKF/6) and pre-malignant oral keratinocytes (pp126) using qPCR. A marked elevation of all KLKs was observed in aggressive SCC25 cells relative to OKF/6 cells. In normal skin, KLKs are involved in desquamation during epidermal differentiation via proteolytic cleavage of the desmosomal cadherin component desmoglein 1 (Dsg1). As loss of cell-cell cohesion is prevalent in tumor metastasis, Dsg1 integrity was evaluated. Results show that SCC25 cells exhibit cleavage of Dsg1, which is blocked by proteinase inhibitor treatment as well as by siRNA silencing of KLK5 expression. Furthermore, cell-cell aggregation assays demonstrate that silencing of KLK5 enforces cell-cell adhesion; conversely, overexpression of KLK5 in normal oral mucosal cells (OKF/6) enhances cell dispersal. These data suggest that KLK5 may promote metastatic dissemination of OSCC by promoting loss of junctional integrity through cleavage of desmoglein 1.     

Expression profiling in squamous carcinoma cells reveals pleiotropic effects of vitamin D3 analog EB1089 signaling on cell proliferation,differentiation, and immune system regulation

The active form of vitamin D3, 1alpha,25-dihydroxyvitamin D3 [1,25-(OH)2D3] is key mediator of calcium homeostasis and is a component of the complex homeostatic system of the skin. 1,25-(OH)2D3 regulates cellular differentiation and proliferation and has broad potential as an anticancer agent. Oligonucleotide microarrays were used to assess profiles of target gene regulation at several points over a 48 h period by the low calcemic 1,25-(OH)2D3 analog EB1089 in human SCC25 head and neck squamous carcinoma cells. One hundred fifty-two targets were identified, composed of 89 up- and 63 down-regulated genes distributed in multiple profiles of regulation. Results are consistent with EB1089 driving SCC25 cells toward a less malignant phenotype, similar to that of basal keratinocytes. Targets identified control inter- and intra-cellular signaling, G protein-coupled receptor function, intracellular redox balance, cell adhesion, and extracellular matrix composition, cell cycle progression, steroid metabolism, and more than 20 genes modulating immune system function. The data indicate that EB1089 performs three key functions of a cancer chemoprevention agent; it is antiproliferative, it induces cellular differentiation, and has potential genoprotective effects. While no evidence was found for gene-specific differences in efficacy of 1,25-(OH)2D3 and EB1089, gene regulation by 1,25-(OH)2D3 was generally more transient. Treatment of cells with 1,25-(OH)2D3 and the cytochrome P450 inhibitor ketoconazole produced profiles of regulation essentially identical to those observed with EB1089 alone, indicating that the more sustained regulation by EB1089 was due to its resistance to inactivation by induced 24-hydroxylase activity. This suggests that differences in action of the two compounds arise more from their relative sensitivities to metabolism than from differing effects on VDR function.     

REG4 Is Highly Expressed in Mucinous Ovarian Cancer: A Potential Novel Serum Biomarker

Preoperative diagnostics of ovarian neoplasms rely on ultrasound imaging and the serum biomarkers CA125 and HE4. However, these markers may be elevated in non-neoplastic conditions and may fail to identify most non-serous epithelial cancer subtypes. The objective of this study was to identify histotype-specific serum biomarkers for mucinous ovarian cancer. The candidate genes with mucinous histotype specific expression profile were identified from publicly available gene-expression databases and further in silico data mining was performed utilizing the MediSapiens database. Candidate biomarker validation was done using qRT-PCR, western blotting and immunohistochemical staining of tumor tissue microarrays. The expression level of the candidate gene in serum was compared to the serum CA125 and HE4 levels in a patient cohort of prospectively collected advanced ovarian cancer. Database searches identified REG4 as a potential biomarker with specificity for the mucinous ovarian cancer subtype. The specific expression within epithelial ovarian tumors was further confirmed by mRNA analysis. Immunohistochemical staining of ovarian tumor tissue arrays showed distinctive cytoplasmic expression pattern only in mucinous carcinomas and suggested differential expression between benign and malignant mucinous neoplasms. Finally, an ELISA based serum biomarker assay demonstrated increased expression only in patients with mucinous ovarian cancer. This study identifies REG4 as a potential serum biomarker for histotype-specific detection of mucinous ovarian cancer and suggests serum REG4 measurement as a non-invasive diagnostic tool for postoperative follow-up of patients with mucinous ovarian cancer.     

Identification of head and neck squamous cell carcinoma biomarker candidates through proteomic analysis of cancer cell secretome

Protein biomarker discovery for early detection of head and neck squamous cell carcinoma (HNSCC) is a crucial unmet need to improve patient outcomes. Mass spectrometry-based proteomics has emerged as a promising tool for identification of biomarkers in different cancer types. Proteins secreted from cancer cells can serve as potential biomarkers for early diagnosis. In the current study, we have used isobaric tag for relative and absolute quantitation (iTRAQ) labeling methodology coupled with high resolution mass spectrometry to identify and quantitate secreted proteins from a panel of head and neck carcinoma cell lines. In all, we identified 2,472 proteins, of which 225 proteins were secreted at higher or lower abundance in HNSCC-derived cell lines. Of these, 148 were present in higher abundance and 77 were present in lower abundance in the cancer-cell derived secretome. We detected a higher abundance of some previously known markers for HNSCC including insulin like growth factor binding protein 3, IGFBP3 (11-fold) and opioid growth factor receptor, OGFR (10-fold) demonstrating the validity of our approach. We also identified several novel secreted proteins in HNSCC including olfactomedin-4, OLFM4 (12-fold) and hepatocyte growth factor activator, HGFA (5-fold). IHC-based validation was conducted in HNSCC using tissue microarrays which revealed overexpression of IGFBP3 and OLFM4 in 70% and 75% of the tested cases, respectively. Our study illustrates quantitative proteomics of secretome as a robust approach for identification of potential HNSCC biomarkers. This article is part of a Special Issue entitled: An Updated Secretome.     

Genetic tumor archeology: microdissection and genetic heterogeneity in squamous and basal cell carcinoma

Carcinogenesis is a multi-step series of somatic genetic events. The complexity of this multi-hit process makes it difficult to determine each single event and the definitive outcome of such events. To investigate the genetic alterations in cancer-related genes, sensitive and reliable detection methods are of major importance for generating relevant results. Another critical issue is the quality of starting material which largely affects the outcome of the analysis. Microdissection of cells defined under the microscope ensures a selection of representative material for subsequent genetic analysis. Skin cancer provides an advantageous model for studying the development of cancer. Detectable lesions occur early during tumor progression, facilitating molecular analysis of the cell populations from both preneoplastic and neoplastic lesions. Alterations of the p53 tumor suppressor gene are very common in non-melanoma skin cancer, and dysregulation of p53 pathways appear to be an early event in the tumor development. A high frequency of epidermal p53 clones has been detected in chronically sun-exposed skin. The abundance of clones containing p53 mutated keratinocytes adjacent to basal cell (BCC) and squamous cell carcinoma (SCC) suggests a role in human skin carcinogenesis. Studies using p53 mutations as a clonality marker have suggested a direct link between actinic keratosis, SCC in situ and invasive SCC. Microdissection-based studies have also shown that different parts of individual BCC tumors can share a common p53 mutation yet differ with respect to additional alterations within the p53 gene, consistent with subclonal development within tumors. Here, we present examples of using well-defined cell populations, including single cells, from complex tissue in combination with molecular tools to reveal features involved in skin carcinogenesis.     

Identification of proteins secreted by head and neck cancer cell lines using LC-MS/MS: Strategy for discovery of candidate serological biomarkers

In search of blood-based biomarkers that would enhance the ability to diagnose head and neck/oral squamous cell carcinoma (HNOSCC) in early stages or predict its prognosis, we analyzed the HNOSCC secretome (ensemble of proteins secreted and/or shed from the tumor cells) for potential biomarkers using proteomic technologies. LC-MS/MS was used to identify proteins in the conditioned media of four HNOSCC cell lines (SCC4, HSC2, SCC38, and AMOSIII); 140 unique proteins were identified on the basis of 5% global false discovery rate, 122 of which were secretory proteins, with 29 being previously reported to be overexpressed in HNOSCC in comparison to normal head and neck tissues. Of these, five proteins including α-enolase, peptidyl prolyl isomerase A/cyclophilin A, 14-3-3 ζ, heterogeneous ribonucleoprotein K, and 14-3-3 σ were detected in the sera of HNOSCC patients by Western blot analysis. Our study provides the evidence that analysis of head and neck cancer cells' secretome is a viable strategy for identifying candidate serological biomarkers for HNOSCC. In future, these biomarkers may be useful in predicting the likelihood of transformation of oral pre-malignant lesions, prognosis of HNOSCC patients and evaluate response to therapy using minimally invasive tests.     

Proteome analysis of human pancreatic cancer cell lines with highly liver metastatic potential by antibody microarray

Antibody microarrays have been successfully used to determine relative abundance of key proteins in various cancers and other diseases. We have previously showed liver metastatic-related genes between the metastatic pancreatic cancer line (SW1990HM) and its parental line (SW1990). In this study, we searched for potential markers for metastatic progression using antibody microarrays. The SpringBio Antibody Microarrays were used to analysis the different proteomes between SW1990HM and SW1990 cells. A standard ≥2.0-fold cutoff value was used to determine differentially expressed proteins and Western blotting analysis further confirmed the results. Antibody microarrays revealed that 40 proteins were reproducibly altered more than 2-fold between the selected variant and its parental counterpart; 14 of the proteins were up-regulated, and 26 were down-regulated. Most of the up-regulated proteins (7/14) play a role in tumor signal transduction, while a number of down-regulated proteins (10/26) function in cell differentiation; this might be crucial for pancreatic cancer metastasis. Four dysregulated proteins were validated by western blotting in the cell lines. Interestingly, the up-regulation of Glucagon and down-regulation of Prolactin were further confirmed in the culture supernatants by western blotting. These proteomic data are valuable for understanding pancreatic cancer metastasis and searching for potential markers of metastatic progression.     

A miR‐34a‐SIRT6 axis in the squamous cell differentiation network

Squamous cell carcinomas (SCCs) are highly heterogeneous tumours, resulting from deranged expression of genes involved in squamous cell differentiation. Here we report that microRNA‐34a (miR‐34a) functions as a novel node in the squamous cell differentiation network, with SIRT6 as a critical target. miR‐34a expression increases with keratinocyte differentiation, while it is suppressed in skin and oral SCCs, SCC cell lines, and aberrantly differentiating primary human keratinocytes (HKCs). Expression of this miRNA is restored in SCC cells, in parallel with differentiation, by reversion of genomic DNA methylation or wild‐type p53 expression. In normal HKCs, the pro‐differentiation effects of increased p53 activity or UVB exposure are miR‐34a‐dependent, and increased miR‐34a levels are sufficient to induce differentiation of these cells both in vitro and in vivo. SIRT6, a sirtuin family member not previously connected with miR‐34a function, is a direct target of this miRNA in HKCs, and SIRT6 down‐modulation is sufficient to reproduce the miR‐34a pro‐differentiation effects. The findings are of likely biological significance, as SIRT6 is oppositely expressed to miR‐34a in normal keratinocytes and keratinocyte‐derived tumours.     

Systematic Identification and Characterization of Novel Human Skin-Associated Genes Encoding Membrane and Secreted Proteins

Through bioinformatics analyses of a human gene expression database representing 105 different tissues and cell types, we identified 687 skin-associated genes that are selectively and highly expressed in human skin. Over 50 of these represent uncharacterized genes not previously associated with skin and include a subset that encode novel secreted and plasma membrane proteins. The high levels of skin-associated expression for eight of these novel therapeutic target genes were confirmed by semi-quantitative real time PCR, western blot and immunohistochemical analyses of normal skin and skin-derived cell lines. Four of these are expressed specifically by epidermal keratinocytes; two that encode G-protein-coupled receptors (GPR87 and GPR115), and two that encode secreted proteins (WFDC5 and SERPINB7). Further analyses using cytokine-activated and terminally differentiated human primary keratinocytes or a panel of common inflammatory, autoimmune or malignant skin diseases revealed distinct patterns of regulation as well as disease associations that point to important roles in cutaneous homeostasis and disease. Some of these novel uncharacterized skin genes may represent potential biomarkers or drug targets for the development of future diagnostics or therapeutics.