Genomics and proteomics approaches to the study of cancer-stroma interactions

Background: The development and progression of cancer depend on its genetic characteristics as well as on the interactions with its microenvironment. Understanding these interactions may contribute to diagnostic and prognostic evaluations and to the development of new cancer therapies. Aiming to investigate potential mechanisms by which the tumor microenvironment might contribute to a cancer phenotype, we evaluated soluble paracrine factors produced by stromal and neoplastic cells which may influence proliferation and gene and protein expression. Methods: The study was carried out on the epithelial cancer cell line (Hep-2) and fibroblasts isolated from a primary oral cancer. We combined a conditioned-medium technique with subtraction hybridization approach, quantitative PCR and proteomics, in order to evaluate gene and protein expression influenced by soluble paracrine factors produced by stromal and neoplastic cells. Results: We observed that conditioned medium from fibroblast cultures (FCM) inhibited proliferation and induced apoptosis in Hep-2 cells. In neoplastic cells, 41 genes and 5 proteins exhibited changes in expression levels in response to FCM and, in fibroblasts, 17 genes and 2 proteins showed down-regulation in response to conditioned medium from Hep-2 cells (HCM). Nine genes were selected and the expression results of 6 down-regulated genes (ARID4A, CALR, GNB2L1, RNF10, SQSTM1, USP9X) were validated by real time PCR. Conclusions: A significant and common denominator in the results was the potential induction of signaling changes associated with immune or inflammatory response in the absence of a specific protein.

Interaction between polymorphisms of the Human Leukocyte Antigen and HPV-16 Variants on the risk of invasive cervical cancer

Background: Persistent infection with oncogenic types of human papillomavirus (HPV) is the major risk factor for invasive cervical cancer (ICC), and non-European variants of HPV-16 are associated with an increased risk of persistence and ICC. HLA class II polymorphisms are also associated with genetic susceptibility to ICC. Our aim is to verify if these associations are influenced by HPV-16 variability. Methods: We characterized HPV-16 variants by PCR in 107 ICC cases, which were typed for HLA-DQA1, DRB1 and DQB1 genes and compared to 257 controls. We measured the magnitude of associations by logistic regression analysis. Results: European ( E), Asian-American ( AA) and African (Af) variants were identified. Here we show that inverse association between DQB1*05 ( adjusted odds ratio [ OR] = 0.66; 95% confidence interval [CI]: 0.39-1.12]) and HPV-16 positive ICC in our previous report was mostly attributable to AA variant carriers ( OR = 0.27; 95% CI: 0.10-0.75). We observed similar proportions of HLA DRB1*1302 carriers in E-P positive cases and controls, but interestingly, this allele was not found in AA cases ( p = 0.03, Fisher exact test). A positive association with DRB1*15 was observed in both groups of women harboring either E ( OR = 2.99; 95% CI: 1.13-7.86) or AA variants ( OR = 2.34; 95% CI: 1.00-5.46). There was an inverse association between DRB1*04 and ICC among women with HPV-16 carrying the 350T [83L] single nucleotide polymorphism in the E6 gene ( OR = 0.27; 95% CI: 0.08-0.96). An inverse association between DQB1*05 and cases carrying 350G (83V) variants was also found ( OR = 0.37; 95% CI: 0.15-0.89). Conclusion: Our results suggest that the association between HLA polymorphism and risk of ICC might be influenced by the distribution of HPV-16 variants.

Identification of FAM46D as a novel cancer/testis antigen using EST data and serological analysis

Cancer/testis Antigens (CTAs) are immunogenic proteins with a restricted expression pattern in normal tissues and aberrant expression in different types of tumors being considered promising candidates for immunotherapy. We used the alignment between EST sequences and the human genome sequence to identify novel CT genes. By examining the EST tissue composition of known CT clusters we defined parameters for the selection of 1184 EST clusters corresponding to putative CT genes. The expression pattern of 70 CT gene candidates was evaluated by RT-PCR in 21 normal tissues, 17 tumor cell lines and 160 primary tumors. We were able to identify 4 CT genes expressed in different types of tumors. The presence of antibodies against the protein encoded by 1 of these 4 CT genes (FAM46D) was exclusively detected in plasma samples from cancer patients. Due to its restricted expression pattern and immunogenicity FAM46D represents a novel target for cancer immunotherapy. (c) 2009 Elsevier Inc. All rights reserved.

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.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Mitochondrial DNA depletion and its correlation with TFAM, TFB1M, TFB2M and POLG in human diffusely infiltrating astrocytomas

Mitochondrial DNA (mtDNA) alterations and their clinical and pathological implications have been analyzed in several human malignancies. A marked decrease in mtDNA copy number along with the increase in malignancy was observed in diffusely infiltrating astrocytomas (24 WHO grade II, 18 grade III, and 78 grade IV or GBM) compared to non-neoplastic brain tissues, being mostly depleted in GBM. Although high relative gene expression levels of mtDNA replication regulators (mitochondrial polymerase catalytic subunit (POLG), transcription factors A (TFAM), B1 (TFB1M) and B2 (TFB2M)) were detected, it cannot successfully revert the mtDNA depletion observed in our samples. On the other hand, a strong correlation among the expression levels of mitochondrial transcription factors corroborates with the TFAM role in the direct control of TFB1M and TFB2M during initiation of mtDNA replication. POLG expression was related to decreased mtDNA copy number, and its overexpression associated with TFAM expression levels also have an impact on long-term survival among GBM patients, interpreted as a potential predictive factor for better prognosis. (C) 2010 Elsevier B.V. and Mitochondria Research Society. All rights reserved.

Bioinformatics construction of the human cell surfaceome

Cell surface proteins are excellent targets for diagnostic and therapeutic interventions. By using bioinformatics tools, we generated a catalog of 3,702 transmembrane proteins located at the surface of human cells (human cell surfaceome). We explored the genetic diversity of the human cell surfaceome at different levels, including the distribution of polymorphisms, conservation among eukaryotic species, and patterns of gene expression. By integrating expression information from a variety of sources, we were able to identify surfaceome genes with a restricted expression in normal tissues and/or differential expression in tumors, important characteristics for putative tumor targets. A high-throughput and efficient quantitative real-time PCR approach was used to validate 593 surfaceome genes selected on the basis of their expression pattern in normal and tumor samples. A number of candidates were identified as potential diagnostic and therapeutic targets for colorectal tumors and glioblastoma. Several candidate genes were also identified as coding for cell surface cancer/testis antigens. The human cell surfaceome will serve as a reference for further studies aimed at characterizing tumor targets at the surface of human cells.

Annexin A1 subcellular expression in laryngeal squamous cell carcinoma

Annexin A1 (ANXA1) is a soluble cytoplasmic protein, moving to membranes when calcium levels are elevated. ANXA1 has also been shown to move to the nucleus or outside the cells, depending on tyrosine-kinase signalling, thus interfering in cytoskeletal organization and cell differentiation, mostly in inflammatory and neoplastic processes. The aim was to investigate subcellular patterns of immunohistochemical expression of ANXA1 in neoplastic and non-neoplastic samples from patients with laryngeal squamous cell carcinomas (LSCC), to elucidate the role of ANXA1 in laryngeal carcinogenesis. Serial analysis of gene expression experiments detected reduced expression of ANXA1 gene in LSCC compared with the corresponding non-neoplastic margins. Quantitative polymerase chain reaction confirmed ANXA1 low expression in 15 LSCC and eight matched normal samples. Thus, we investigated subcellular patterns of immunohistochemical expression of ANXA1 in 241 paraffin-embedded samples from 95 patients with LSCC. The results showed ANXA1 down-regulation in dysplastic, tumourous and metastatic lesions and provided evidence for the progressive migration of ANXA1 from the nucleus towards the membrane during laryngeal tumorigenesis. ANXA1 dysregulation was observed early in laryngeal carcinogenesis, in intra-epithelial neoplasms; it was not found related to prognostic parameters, such as nodal metastases.

Association of EGFR c.2073A > T polymorphism with decreased risk of diffusely infiltrating astrocytoma in a Brazilian case-control study

Epidermal growth factor receptor (EGFR) gene overexpression has been implicated in the development of many types of tumors, including glioblastomas, the most frequent diffusely infiltrating astrocytomas. However, little is known about the influence of the polymorphisms of EGFR on EGFR production and/or activity, possibly modulating the susceptibility to astrocytomas. This study aimed to examine the association of two EGFR promoter polymorphisms (c.-191C > A and c.-216G > T) and the c.2073A > T polymorphism located in exon 16 with susceptibility to astrocytomas, EGFR gene expression and survival in a case-control study of 193 astrocytoma patients and 200 cancer-free controls. We found that the variant TT genotype of the EGFR c.2073A > T polymorphism was associated with a significantly decreased risk of astrocytoma when compared with the AA genotype [sex- and age-adjusted odds ratio 0.51, 95% confidence interval 0.26-0.98]. No association of the two promoter EGFR polymorphisms (or combinations of these polymorphisms) and risk of astrocytomas, EGFR expression or survival was found. Our findings suggest that modulation of the EGFR c.2073A > T polymorphism could play a role in future therapeutic approaches to astrocytoma. (Int J Biol Markers 2008; 23: 140-6)

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

Proteomic analysis of low- to high-grade astrocytomas reveals an alteration of the expression level of raf kinase inhibitor protein and nucleophosmin

Proteomic approaches have been useful for the identification of aberrantly expressed proteins in complex diseases such as cancer. These proteins are not only potential disease biomarkers, but also targets for therapy. The aim of this study was to identify differentially expressed proteins in diffuse astrocytoma grade II, anaplastic astrocytoma grade III and glioblastoma multiforme grade IV in human tumor samples and in non-neoplastic brain tissue as control using 2-DE and MS. Tumor and control brain tissue dissection was guided by histological hematoxylin/eosin tissue sections to provide more than 90% of tumor cells and astrocytes. Six proteins were detected as up-regulated in higher grade astrocytomas and the most important finding was nucleophosmin (NPM) (p < 0.05), whereas four proteins were down-regulated, among them raf kinase inhibitor protein (RKIP) (p < 0.05). We report here for the first time the alteration of NPM and RKIP expression in brain cancer. Our focus on these proteins was due to the fact that they are involved in the PI3K/AKT/mTOR and RAS/RAF/MAPK pathways, known for their contribution to the development and progression of gliomas. The proteomic data for NPM and RKIP were confirmed by Western blot, quantitative real-time PCR and immunohistochemistry. Due to the participation of NPM and RKIP in uncontrolled proliferation and evasion of apoptosis, these proteins are likely targets for drug development.

Alternative splicing and genetic diversity: silencers are more frequently modified by SNVs associated with alternative exon/intron borders

With the availability of a large amount of genomic data it is expected that the influence of single nucleotide variations (SNVs) in many biological phenomena will be elucidated. Here, we approached the problem of how SNVs affect alternative splicing. First, we observed that SNVs and exonic splicing regulators (ESRs) independently show a biased distribution in alternative exons. More importantly, SNVs map more frequently in ESRs located in alternative exons than in ESRs located in constitutive exons. By looking at SNVs associated with alternative exon/intron borders (by their common presence in the same cDNA molecule), we observed that a specific type of ESR, the exonic splicing silencers (ESSs), are more frequently modified by SNVs. Our results establish a clear association between genetic diversity and alternative splicing involving ESSs.

Gene expression profiling reveals molecular marker candidates of laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma is very common in head and neck cancer, with high mortality rates and poor prognosis. In this study, we compared expression profiles of clinical samples from 13 larynx tumors and 10 non-neoplastic larynx tissues using a custom-built cDNA microarray containing 331 probes for 284 genes previously identified by informatics analysis of EST databases as markers of head and neck tumors. Thirty-five genes showed statistically significant differences (SNR >= 11.01, p <= 0.001) in the expression between tumor and non-tumor larynx tissue samples. Functional annotation indicated that these genes are involved in cellular processes relevant to the cancer phenotype, such as apoptosis, cell cycle, DNA repair, proteolysis, protease inhibition, signal transduction and transcriptional regulation. Six of the identified transcripts map to intronic regions of protein-coding genes and may comprise non-annotated exons or as yet uncharacterized long ncRNAs with a regulatory role in the gene expression program of larynx tissue. The differential expression of 10 of these genes (ADCY6, AES, AL2SCR3, CRR9, CSTB, DUSP1, MAP3K5, PLAT, UBL1 and ZNF706) was independently confirmed by quantitative real-time RT-PCR. Among these, the CSTB gene product has cysteine protease inhibitor activity that has been associated with an antimetastatic function. Interestingly, CSTB showed a low expression in the tumor samples analyzed (p<0.0001). The set of genes identified here contribute to a better understanding of the molecular basis of larynx cancer, and provide candidate markers for improving diagnosis, prognosis and treatment of this carcinoma.

SOX2 is an amplified lineage-survival oncogene in lung and esophageal squamous cell carcinomas

Lineage-survival oncogenes are activated by somatic DNA alterations in cancers arising from the cell lineages in which these genes play a role in normal development(1,2). Here we show that a peak of genomic amplification on chromosome 3q26.33 found in squamous cell carcinomas (SCCs) of the lung and esophagus contains the transcription factor gene SOX2, which is mutated in hereditary human esophageal malformations(3), is necessary for normal esophageal squamous development(4), promotes differentiation and proliferation of basal tracheal cells(5) and cooperates in induction of pluripotent stem cells(6-8). SOX2 expression is required for proliferation and anchorage-independent growth of lung and esophageal cell lines, as shown by RNA interference experiments. Furthermore, ectopic expression of SOX2 here cooperated with FOXE1 or FGFR2 to transform immortalized tracheobronchial epithelial cells. SOX2-driven tumors show expression of markers of both squamous differentiation and pluripotency. These characteristics identify SOX2 as a lineage-survival oncogene in lung and esophageal SCC.

The Genetic Landscape of the Childhood Cancer Medulloblastoma

Medulloblastoma (MB) is the most common malignant brain tumor of children. To identify the genetic alterations in this tumor type, we searched for copy number alterations using high-density microarrays and sequenced all known protein-coding genes and microRNA genes using Sanger sequencing in a set of 22 MBs. We found that, on average, each tumor had 11 gene alterations, fewer by a factor of 5 to 10 than in the adult solid tumors that have been sequenced to date. In addition to alterations in the Hedgehog and Wnt pathways, our analysis led to the discovery of genes not previously known to be altered in MBs. Most notably, inactivating mutations of the histone-lysine N-methyltransferase genes MLL2 or MLL3 were identified in 16% of MB patients. These results demonstrate key differences between the genetic landscapes of adult and childhood cancers, highlight dysregulation of developmental pathways as an important mechanism underlying MBs, and identify a role for a specific type of histone methylation in human tumorigenesis.