Kristian Ovaska

Integrated network analysis platform for protein-protein interactions

There is an increasing demand for network analysis of protein-protein interactions (PPIs). We introduce a web-based protein interaction network analysis platform (PINA), which integrates PPI data from six databases and provides network construction, filtering, analysis and visualization tools. We demonstrated the advantages of PINA by analyzing two human PPI networks; our results suggested a link between LKB1 and TGFβ signaling, and revealed possible competitive interactors of p53 and c-Jun.

Dual role of FoxA1 in androgen receptor binding to chromatin, androgen signalling and prostate cancer

High androgen receptor (AR) level in primary tumour predicts increased prostate cancer‐specific mortality. However, the mechanisms that regulate AR function in prostate cancer are poorly known. We report here a new paradigm for the forkhead protein FoxA1 action in androgen signalling. Besides pioneering the AR pathway, FoxA1 depletion elicited extensive redistribution of AR‐binding sites (ARBs) on LNCaP‐1F5 cell chromatin that was commensurate with changes in androgen‐dependent gene expression signature. We identified three distinct classes of ARBs and androgen‐responsive genes: (i) independent of FoxA1, (ii) pioneered by FoxA1 and (iii) masked by FoxA1 and functional upon FoxA1 depletion. FoxA1 depletion also reprogrammed AR binding in VCaP cells, and glucocorticoid receptor binding and glucocorticoid‐dependent signalling in LNCaP‐1F5 cells. Importantly, FoxA1 protein level in primary prostate tumour had significant association to disease outcome; high FoxA1 level was associated with poor prognosis, whereas low FoxA1 level, even in the presence of high AR expression, predicted good prognosis. The role of FoxA1 in androgen signalling and prostate cancer is distinctly different from that in oestrogen signalling and breast cancer.

Large-scale data integration framework provides a comprehensive view on glioblastoma multiforme.

BackgroundCoordinated efforts to collect large-scale data sets provide a basis for systems level understanding of complex diseases. In order to translate these fragmented and heterogeneous data sets into knowledge and medical benefits, advanced computational methods for data analysis, integration and visualization are needed.MethodsWe introduce a novel data integration framework, Anduril, for translating fragmented large-scale data into testable predictions. The Anduril framework allows rapid integration of heterogeneous data with state-of-the-art computational methods and existing knowledge in bio-databases. Anduril automatically generates thorough summary reports and a website that shows the most relevant features of each gene at a glance, allows sorting of data based on different parameters, and provides direct links to more detailed data on genes, transcripts or genomic regions. Anduril is open-source; all methods and documentation are freely available.ResultsWe have integrated multidimensional molecular and clinical data from 338 subjects having glioblastoma multiforme, one of the deadliest and most poorly understood cancers, using Anduril. The central objective of our approach is to identify genetic loci and genes that have significant survival effect. Our results suggest several novel genetic alterations linked to glioblastoma multiforme progression and, more specifically, reveal Moesin as a novel glioblastoma multiforme-associated gene that has a strong survival effect and whose depletion in vitro significantly inhibited cell proliferation. All analysis results are available as a comprehensive website.ConclusionsOur results demonstrate that integrated analysis and visualization of multidimensional and heterogeneous data by Anduril enables drawing conclusions on functional consequences of large-scale molecular data. Many of the identified genetic loci and genes having significant survival effect have not been reported earlier in the context of glioblastoma multiforme. Thus, in addition to generally applicable novel methodology, our results provide several glioblastoma multiforme candidate genes for further studies.Anduril is available at http://csbi.ltdk.helsinki.fi/anduril/The glioblastoma multiforme analysis results are available at http://csbi.ltdk.helsinki.fi/anduril/tcga-gbm/

FoxA1 Specifies Unique Androgen and Glucocorticoid Receptor Binding Events in Prostate Cancer Cells

The forkhead protein FoxA1 has functions other than a pioneer factor, in that its depletion brings about a significant redistribution in the androgen receptor (AR) and glucocorticoid receptor (GR) cistromes. In this study, we found a novel function for FoxA1 in defining the cell-type specificity of AR- and GR-binding events in a distinct fashion, namely, for AR in LNCaP-1F5 cells and for GR in VCaP cells. We also found different, cell-type and receptor-specific compilations of cis-elements enriched adjacent to the AR- and GR-binding sites. The AR pathway is central in prostate cancer biology, but the role of GR is poorly known. We find that AR and GR cistromes and transcription programs exhibit significant overlap, and GR regulates a large number of genes considered to be AR pathway-specific. This raises questions about the role of GR in maintaining the AR pathway under androgen-deprived conditions in castration-resistant prostate cancer patients. However, in the presence of androgen, ligand-occupied GR acts as a partial antiandrogen and attenuates the AR-dependent transcription program. .

Fast Gene Ontology based clustering for microarray experiments

BackgroundAnalysis of a microarray experiment often results in a list of hundreds of disease-associated genes. In order to suggest common biological processes and functions for these genes, Gene Ontology annotations with statistical testing are widely used. However, these analyses can produce a very large number of significantly altered biological processes. Thus, it is often challenging to interpret GO results and identify novel testable biological hypotheses.ResultsWe present fast software for advanced gene annotation using semantic similarity for Gene Ontology terms combined with clustering and heat map visualisation. The methodology allows rapid identification of genes sharing the same Gene Ontology cluster.ConclusionOur R based semantic similarity open-source package has a speed advantage of over 2000-fold compared to existing implementations. From the resulting hierarchical clustering dendrogram genes sharing a GO term can be identified, and their differences in the gene expression patterns can be seen from the heat map. These methods facilitate advanced annotation of genes resulting from data analysis.

Transmembrane Prostatic Acid Phosphatase (TMPAP) Interacts with Snapin and Deficient Mice Develop Prostate Adenocarcinoma

The molecular mechanisms underlying prostate carcinogenesis are poorly understood. Prostatic acid phosphatase (PAP), a prostatic epithelial secretion marker, has been linked to prostate cancer since the 1930s. However, the contribution of PAP to the disease remains controversial. We have previously cloned and described two isoforms of this protein, a secretory (sPAP) and a transmembrane type-I (TMPAP). The goal in this work was to understand the physiological function of TMPAP in the prostate. We conducted histological, ultra-structural and genome-wide analyses of the prostate of our PAP-deficient mouse model (PAP−/−) with C57BL/6J background. The PAP−/− mouse prostate showed the development of slow-growing non-metastatic prostate adenocarcinoma. In order to find out the mechanism behind, we identified PAP-interacting proteins byyeast two-hybrid assays and a clear result was obtained for the interaction of PAP with snapin, a SNARE-associated protein which binds Snap25 facilitating the vesicular membrane fusion process. We confirmed this interaction by co-localization studies in TMPAP-transfected LNCaP cells (TMPAP/LNCaP cells) and in vivo FRET analyses in transient transfected LNCaP cells. The differential gene expression analyses revealed the dysregulation of the same genes known to be related to synaptic vesicular traffic. Both TMPAP and snapin were detected in isolated exosomes. Our results suggest that TMPAP is involved in endo-/exocytosis and disturbed vesicular traffic is a hallmark of prostate adenocarcinoma.

Neuron navigator 3 alterations in nervous system tumors associate with tumor malignancy grade and prognosis

Copy number changes or reduced expression of the Neuron navigator 3 (NAV3) gene occurs in neuroblastomas and malignancies of epithelial or lymphoid origin. To elucidate whether NAV3 has a role in the tumorigenesis of nervous system tumors in general, we studied central and peripheral nervous system tumors for NAV3 copy number changes. In search for common tumorigenic denominators, we analyzed 113 central and peripheral nervous system tumors, including glial tumors (grades I–IV gliomas), medulloblastomas, and neuroblastomas. NAV3 copy number changes were studied by fluorescence in situ hybridization and correlated to survival analyses. To identify target genes of NAV3 deletion, NAV3 was silenced by siRNA in glioblastoma cell lines and gene expression profiles were analyzed by Agilent 4×44k dual‐color microarrays. Selected upregulations were confirmed by immunohistochemistry and quantitative polymerase chain reaction. We found NAV3 amplifications to dominate in neuronally differentiated tumors, whereas glial tumors showed almost equal proportions of NAV3 deletion and amplification. However, Grade IV gliomas had more frequent NAV3 deletions than grades I–III gliomas. Silencing of NAV3 in glioma cell lines led to the upregulation of receptor genes associated with gonadotropin‐releasing hormone and Jak‐Stat signaling pathways. Kaplan–Meier analysis of the entire clinical tumor material showed association between NAV3 amplifications and favorable prognosis, as well as NAV3 deletions and unfavorable prognosis. With Cox regression model, a hazard ratio of 0.51 was observed for NAV3 amplifications and 1.36 for NAV3 deletions. We conclude that NAV3 may be a potential new prognostic biomarker and a potential therapeutic target.

Genomic Region Operation Kit for Flexible Processing of Deep Sequencing Data

Computational analysis of data produced in deep sequencing (DS) experiments is challenging due to large data volumes and requirements for flexible analysis approaches. Here, we present a mathematical formalism based on set algebra for frequently performed operations in DS data analysis to facilitate translation of biomedical research questions to language amenable for computational analysis. With the help of this formalism, we implemented the Genomic Region Operation Kit (GROK), which supports various DS-related operations such as preprocessing, filtering, file conversion, and sample comparison. GROK provides high-level interfaces for R, Python, Lua, and command line, as well as an extension C++ API. It supports major genomic file formats and allows storing custom genomic regions in efficient data structures such as red-black trees and SQL databases. To demonstrate the utility of GROK, we have characterized the roles of two major transcription factors (TFs) in prostate cancer using data from 10 DS experiments. GROK is freely available with a user guide from http://csbi.ltdk.helsinki.fi/grok/.

NAV3 copy number changes and target genes in basal and squamous cell cancers

Abstract:  The neuron navigator 3 (NAV3) gene on chromosome 12q21 encodes a microtubule plus end tracking protein and belongs to the navigator family of cytoskeletal regulators. Loss of heterozygosity on 12q has previously been suggested to be associated with poor prognosis in cancers of epithelial origin. In this study, we characterized copy number changes of NAV3 in 24 basal cell cancers (BCCs), eight squamous cell cancers (SCCs) and eight non‐malignant inflammatory skin lesions by fluorescent in situ hybridization. To identify genes affected by NAV3, we used oligo siRNA gene silencing and gene microarrays to analyse gene expression profiles at several time points post‐transfection in primary human keratinocytes. We found NAV3 copy number loss and decreased protein expression in 21% of the BCCs and 25% of the SCCs. In the nodular/superficial BCC subgroup, low‐level NAV3 amplification was also observed. NAV3 aberrations were independent of the known chromosome 6 amplifications in BCC. Chromosome 12 polysomy was detected in 33% and 25% of the invasive type of BCC and SCC, respectively. Silencing of NAV3 in primary human keratinocytes revealed 22 differentially expressed genes, mostly related to inflammation. The most relevant of these were validated with qPCR or immunohistochemistry. This pilot study suggests that NAV3 is a novel cancer‐associated gene that contributes to the pathogenesis of a subgroup of BCC and SCC.

Integrative Analysis of Deep Sequencing Data Identifies Estrogen Receptor Early Response Genes and Links ATAD3B to Poor Survival in Breast Cancer

Identification of responsive genes to an extra-cellular cue enables characterization of pathophysiologically crucial biological processes. Deep sequencing technologies provide a powerful means to identify responsive genes, which creates a need for computational methods able to analyze dynamic and multi-level deep sequencing data. To answer this need we introduce here a data-driven algorithm, SPINLONG, which is designed to search for genes that match the user-defined hypotheses or models. SPINLONG is applicable to various experimental setups measuring several molecular markers in parallel. To demonstrate the SPINLONG approach, we analyzed ChIP-seq data reporting PolII, estrogen receptor (), H3K4me3 and H2A.Z occupancy at five time points in the MCF-7 breast cancer cell line after estradiol stimulus. We obtained 777 early responsive genes and compared the biological functions of the genes having binding within 20 kb of the transcription start site (TSS) to genes without such binding site. Our results show that the non-genomic action of via the MAPK pathway, instead of direct binding, may be responsible for early cell responses to activation. Our results also indicate that the responsive genes triggered by the genomic pathway are transcribed faster than those without binding sites. The survival analysis of the 777 responsive genes with 150 primary breast cancer tumors and in two independent validation cohorts indicated the ATAD3B gene, which does not have binding site within 20 kb of its TSS, to be significantly associated with poor patient survival.