Search

Prostate cancer treatment resistance is a significant challenge facing the field. Genomic and transcriptomic profiling have partially elucidated the mechanisms through which cancer cells escape treatment, but their relation toward the tumor microenvironment (TME) remains elusive. Here we present a comprehensive transcriptomic landscape of the prostate TME at multiple points in the standard treatment timeline employing single-cell RNA-sequencing and spatial transcriptomics data from 120 patients. We identify club-like cells as a key epithelial cell subtype that acts as an interface between the prostate and the immune system. Tissue areas enriched with club-like cells have depleted androgen signaling and upregulated expression of luminal progenitor cell markers. Club-like cells display a senescence-associated secretory phenotype and their presence is linked to increased polymorphonuclear myeloid-derived suppressor cell (PMN-MDSC) activity. Our results indicate that club-like cells are associated with myeloid inflammation previously linked to androgen deprivation therapy resistance, providing a rationale for their therapeutic targeting., Competing Interests: Competing interests C.T.A.P.’s employer may gain commercially from licensing data to Artera AI. G.A. received personal fees, grants, and travel support from Janssen and Astellas Pharma; personal fees or travel support from Pfizer, Novartis/AAA, Bayer Healthcare Pharmaceuticals, AstraZeneca, and Sanofi-Aventis; in addition, G.A.’s former employer, The Institute of Cancer Research, receives royalty income from abiraterone and G.A. receives a share of this income through the Institute’s Rewards to Discoverers Scheme. G.A. has received research funding (institutional) from Janssen, Astellas Pharma, and Novartis. All other authors declare no potential conflicts of interest., (© 2024. The Author(s).)

Prostate tumor heterogeneity is a major obstacle when studying the biological mechanisms of molecular markers. Increased gene expression levels of secreted frizzled-related protein 4 (SFRP4) is a biomarker in aggressive prostate cancer. To understand how SFRP4 relates to prostate cancer we performed comprehensive spatial and multiomics analysis of the same prostate cancer tissue samples. The experimental workflow included spatial transcriptomics, bulk transcriptomics, proteomics, DNA methylomics and tissue staining. SFRP4 mRNA was predominantly located in cancer stroma, produced by fibroblasts and smooth muscle cells, and co-expressed with extracellular matrix components. We also confirmed that higher SFRP4 gene expression is associated with cancer aggressiveness. Gene expression of SFRP4 was affected by gene promotor methylation. Surprisingly, the high mRNA levels did not reflect SFRP4 protein levels, which was much lower. This study contributes previously unknown insights of SFRP4 mRNA in the prostate tumor environment that potentially can improve diagnosis and treatment., (© 2024. The Author(s).)

The DNA damage inducible SOS response in bacteria serves to increase survival of the species at the cost of mutagenesis. The SOS response first initiates error-free repair followed by error-prone repair. Here, we have employed a multi-omics approach to elucidate the temporal coordination of the SOS response. Escherichia coli was grown in batch cultivation in bioreactors to ensure highly controlled conditions, and a low dose of the antibiotic ciprofloxacin was used to activate the SOS response while avoiding extensive cell death. Our results show that expression of genes involved in error-free and error-prone repair were both induced shortly after DNA damage, thus, challenging the established perception that the expression of error-prone repair genes is delayed. By combining transcriptomics and a sub-proteomics approach termed signalomics, we found that the temporal segregation of error-free and error-prone repair is primarily regulated after transcription, supporting the current literature. Furthermore, the heterology index (i.e., the binding affinity of LexA to the SOS box) was correlated to the maximum increase in gene expression and not to the time of induction of SOS genes. Finally, quantification of metabolites revealed increasing pyrimidine pools as a late feature of the SOS response. Our results elucidate how the SOS response is coordinated, showing a rapid transcriptional response and temporal regulation of mutagenesis on the protein and metabolite levels., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Bergum, Singleton, Røst, Bodein, Scott-Boyer, Rye, Droit, Bruheim and Otterlei.)

Background: The dual upregulation of TOP2A and EZH2 gene expression has been proposed as a biomarker for recurrence in prostate cancer patients to be treated with radical prostatectomy. A low tissue level of the metabolite citrate has additionally been connected to aggressive disease and recurrence in this patient group. However, for radiotherapy prostate cancer patients, few prognostic biomarkers have been suggested. The main aim of this study was to use an integrated tissue analysis to evaluate metabolites and expression of TOP2A and EZH2 as predictors for recurrence among radiotherapy patients., Methods: From 90 prostate cancer patients (56 received neoadjuvant hormonal treatment), 172 transrectal ultrasound-guided (TRUS) biopsies were collected prior to radiotherapy. Metabolic profiles were acquired from fresh frozen TRUS biopsies using high resolution-magic angle spinning MRS. Histopathology and immunohistochemistry staining for TOP2A and EZH2 were performed on TRUS biopsies containing cancer cells (n = 65) from 46 patients, where 24 of these patients (n = 31 samples) received hormonal treatment. Eleven radical prostatectomy cohorts of a total of 2059 patients were used for validation in a meta-analysis., Results: Among radiotherapy patients with up to 11 years of follow-up, a low level of citrate was found to predict recurrence, p = 0.001 (C-index = 0.74). Citrate had a higher predictive ability compared with individual clinical variables, highlighting its strength as a potential biomarker for recurrence. The dual upregulation of TOP2A and EZH2 was suggested as a biomarker for recurrence, particularly for patients not receiving neoadjuvant hormonal treatment, p = 0.001 (C-index = 0.84). While citrate was a statistically significant biomarker independent of hormonal treatment status, the current study indicated a potential of glutamine, glutamate and choline as biomarkers for recurrence among patients receiving neoadjuvant hormonal treatment, and glucose among patients not receiving neoadjuvant hormonal treatment., Conclusion: Using an integrated approach, our study shows the potential of citrate and the dual upregulation of TOP2A and EZH2 as biomarkers for recurrence among radiotherapy patients., (© 2022 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.)

Mitochondrial activity in cancer cells has been central to cancer research since Otto Warburg first published his thesis on the topic in 1956. Although Warburg proposed that oxidative phosphorylation in the tricarboxylic acid (TCA) cycle was perturbed in cancer, later research has shown that oxidative phosphorylation is activated in most cancers, including prostate cancer (PCa). However, more detailed knowledge on mitochondrial metabolism and metabolic pathways in cancers is still lacking. In this study we expand our previously developed method for analyzing functional homologous proteins (FunHoP), which can provide a more detailed view of metabolic pathways. FunHoP uses results from differential expression analysis of RNA-Seq data to improve pathway analysis. By adding information on subcellular localization based on experimental data and computational predictions we can use FunHoP to differentiate between mitochondrial and non-mitochondrial processes in cancerous and normal prostate cell lines. Our results show that mitochondrial pathways are upregulated in PCa and that splitting metabolic pathways into mitochondrial and non-mitochondrial counterparts using FunHoP adds to the interpretation of the metabolic properties of PCa cells., Competing Interests: The authors have declared that no competing interests exist.

Background: Locally advanced breast cancer is a heterogeneous disease with respect to response to neoadjuvant chemotherapy (NACT) and survival. It is currently not possible to accurately predict who will benefit from the specific types of NACT. DNA methylation is an epigenetic mechanism known to play an important role in regulating gene expression and may serve as a biomarker for treatment response and survival. We investigated the potential role of DNA methylation as a prognostic marker for long-term survival (> 5 years) after NACT in breast cancer., Methods: DNA methylation profiles of pre-treatment (n = 55) and post-treatment (n = 75) biopsies from 83 women with locally advanced breast cancer were investigated using the Illumina HumanMethylation450 BeadChip. The patients received neoadjuvant treatment with epirubicin and/or paclitaxel. Linear mixed models were used to associate DNA methylation to treatment response and survival based on clinical response to NACT (partial response or stable disease) and 5-year survival, respectively. LASSO regression was performed to identify a risk score based on the statistically significant methylation sites and Kaplan-Meier curve analysis was used to estimate survival probabilities using ten years of survival follow-up data. The risk score developed in our discovery cohort was validated in an independent validation cohort consisting of paired pre-treatment and post-treatment biopsies from 85 women with locally advanced breast cancer. Patients included in the validation cohort were treated with either doxorubicin or 5-FU and mitomycin NACT., Results: DNA methylation patterns changed from before to after NACT in 5-year survivors, while no significant changes were observed in non-survivors or related to treatment response. DNA methylation changes included an overall loss of methylation at CpG islands and gain of methylation in non-CpG islands, and these changes affected genes linked to transcription factor activity, cell adhesion and immune functions. A risk score was developed based on four methylation sites which successfully predicted long-term survival in our cohort (p = 0.0034) and in an independent validation cohort (p = 0.049)., Conclusion: Our results demonstrate that DNA methylation patterns in breast tumors change in response to NACT. These changes in DNA methylation show potential as prognostic biomarkers for breast cancer survival., (© 2022. The Author(s).)

High secretion of the metabolites citrate and spermine is a unique hallmark for normal prostate epithelial cells, and is reduced in aggressive prostate cancer. However, the identity of the genes controlling this biological process is mostly unknown. In this study, we have created a gene signature of 150 genes connected to citrate and spermine secretion in the prostate. We have computationally integrated metabolic measurements with multiple transcriptomics datasets from the public domain, including 3826 tissue samples from prostate and prostate cancer. The accuracy of the signature is validated by its unique enrichment in prostate samples and prostate epithelial tissue compartments. The signature highlights genes AZGP1 , ANPEP and metallothioneins with zinc-binding properties not previously studied in the prostate, and the expression of these genes are reduced in more aggressive cancer lesions. However, the absence of signature enrichment in common prostate model systems can make it challenging to study these genes mechanistically., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

Multiple myeloma (MM) is an incurable hematologic malignancy resulting from the clonal expansion of plasma cells. MM cells are interacting with components of the bone marrow microenvironment such as cytokines to survive and proliferate. Phosphatase of regenerating liver (PRL)-3, a cytokine-induced oncogenic phosphatase, is highly expressed in myeloma patients and is a mediator of metabolic reprogramming of cancer cells. To find novel pathways and genes regulated by PRL-3, we characterized the global transcriptional response to PRL-3 overexpression in two MM cell lines. We used pathway enrichment analysis to identify pathways regulated by PRL-3. We further confirmed the hits from the enrichment analysis with in vitro experiments and investigated their function. We found that PRL-3 induced expression of genes belonging to the type 1 interferon (IFN-I) signaling pathway due to activation of signal transducer and activator of transcription (STAT) 1 and STAT2. This activation was independent of autocrine IFN-I secretion. The increase in STAT1 and STAT2 did not result in any of the common consequences of increased IFN-I or STAT1 signaling in cancer. Knockdown of STAT1/2 did not affect the viability of the cells, but decreased PRL-3-induced glycolysis. Interestingly, glucose metabolism contributed to the activation of STAT1 and STAT2 and expression of IFN-I-stimulated genes in PRL-3-overexpressing cells. In summary, we describe a novel signaling circuit where the key IFN-I-activated transcription factors STAT1 and STAT2 are important drivers of the increase in glycolysis induced by PRL-3. Subsequently, increased glycolysis regulates the IFN-I-stimulated genes by augmenting the activation of STAT1/2., (© 2021 The Authors. The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Cytoscape is often used for visualization and analysis of metabolic pathways. For example, based on KEGG data, a reader for KEGG Markup Language (KGML) is used to load files into Cytoscape. However, although multiple genes can be responsible for the same reaction, the KGML-reader KEGGScape only presents the first listed gene in a network node for a given reaction. This can lead to incorrect interpretations of the pathways. Our new method, FunHoP, shows all possible genes in each node, making the pathways more complete. FunHoP collapses all genes in a node into one measurement using read counts from RNA-seq. Assuming that activity for an enzymatic reaction mainly depends upon the gene with the highest number of reads, and weighting the reads on gene length and ratio, a new expression value is calculated for the node as a whole. Differential expression at node level is then applied to the networks. Using prostate cancer as model, we integrate RNA-seq data from two patient cohorts with metabolism data from literature. Here we show that FunHoP gives more consistent pathways that are easier to interpret biologically. Code and documentation for running FunHoP can be found at https://github.com/kjerstirise/FunHoP., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Background: Prostate cancer tissues are inherently heterogeneous, which presents a challenge for metabolic profiling using traditional bulk analysis methods that produce an averaged profile. The aim of this study was therefore to spatially detect metabolites and lipids on prostate tissue sections by using mass spectrometry imaging (MSI), a method that facilitates molecular imaging of heterogeneous tissue sections, which can subsequently be related to the histology of the same section., Methods: Here, we simultaneously obtained metabolic and lipidomic profiles in different prostate tissue types using matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MSI. Both positive and negative ion mode were applied to analyze consecutive sections from 45 fresh-frozen human prostate tissue samples (N = 15 patients). Mass identification was performed with tandem MS., Results: Pairwise comparisons of cancer, non-cancer epithelium, and stroma revealed several metabolic differences between the tissue types. We detected increased levels of metabolites crucial for lipid metabolism in cancer, including metabolites involved in the carnitine shuttle, which facilitates fatty acid oxidation, and building blocks needed for lipid synthesis. Metabolites associated with healthy prostate functions, including citrate, aspartate, zinc, and spermine had lower levels in cancer compared to non-cancer epithelium. Profiling of stroma revealed higher levels of important energy metabolites, such as ADP, ATP, and glucose, and higher levels of the antioxidant taurine compared to cancer and non-cancer epithelium., Conclusions: This study shows that specific tissue compartments within prostate cancer samples have distinct metabolic profiles and pinpoint the advantage of methodology providing spatial information compared to bulk analysis. We identified several differential metabolites and lipids that have potential to be developed further as diagnostic and prognostic biomarkers for prostate cancer. Spatial and rapid detection of cancer-related analytes showcases MALDI-TOF MSI as a promising and innovative diagnostic tool for the clinic.

We present model-based analysis for ChIA-PET (MACPET), which analyzes paired-end read sequences provided by ChIA-PET for finding binding sites of a protein of interest. MACPET uses information from both tags of each PET and searches for binding sites in a two-dimensional space, while taking into account different noise levels in different genomic regions. MACPET shows favorable results compared with MACS in terms of motif occurrence and spatial resolution. Furthermore, significant binding sites discovered by MACPET are involved in a higher number of significant three-dimensional interactions than those discovered by MACS. MACPET is freely available on Bioconductor. ChIA-PET; MACPET; Model-based clustering; Paired-end tags; Peak-calling algorithm., (© The Author 2019. Published by Oxford University Press.)

Background: Prostate cancer (PCa) has the highest incidence rates of cancers in men in western countries. Unlike several other types of cancer, PCa has few genetic drivers, which has led researchers to look for additional epigenetic and transcriptomic contributors to PCa development and progression. Especially datasets on DNA methylation, the most commonly studied epigenetic marker, have recently been measured and analysed in several PCa patient cohorts. DNA methylation is most commonly associated with downregulation of gene expression. However, positive associations of DNA methylation to gene expression have also been reported, suggesting a more diverse mechanism of epigenetic regulation. Such additional complexity could have important implications for understanding prostate cancer development but has not been studied at a genome-wide scale., Results: In this study, we have compared three sets of genome-wide single-site DNA methylation data from 870 PCa and normal tissue samples with multi-cohort gene expression data from 1117 samples, including 532 samples where DNA methylation and gene expression have been measured on the exact same samples. Genes were classified according to their corresponding methylation and expression profiles. A large group of hypermethylated genes was robustly associated with increased gene expression (UPUP group) in all three methylation datasets. These genes demonstrated distinct patterns of correlation between DNA methylation and gene expression compared to the genes showing the canonical negative association between methylation and expression (UPDOWN group). This indicates a more diversified role of DNA methylation in regulating gene expression than previously appreciated. Moreover, UPUP and UPDOWN genes were associated with different compartments - UPUP genes were related to the structures in nucleus, while UPDOWN genes were linked to extracellular features., Conclusion: We identified a robust association between hypermethylation and upregulation of gene expression when comparing samples from prostate cancer and normal tissue. These results challenge the classical view where DNA methylation is always associated with suppression of gene expression, which underlines the importance of considering corresponding expression data when assessing the downstream regulatory effect of DNA methylation.

Sequencing technologies have changed not only our approaches to classical genetics, but also the field of epigenetics. Specific methods allow scientists to identify novel genome-wide epigenetic patterns of DNA methylation down to single-nucleotide resolution. DNA methylation is the most researched epigenetic mark involved in various processes in the human cell, including gene regulation and development of diseases, such as cancer. Increasing numbers of DNA methylation sequencing datasets from human genome are produced using various platforms-from methylated DNA precipitation to the whole genome bisulfite sequencing. Many of those datasets are fully accessible for repeated analyses. Sequencing experiments have become routine in laboratories around the world, while analysis of outcoming data is still a challenge among the majority of scientists, since in many cases it requires advanced computational skills. Even though various tools are being created and published, guidelines for their selection are often not clear, especially to non-bioinformaticians with limited experience in computational analyses. Separate tools are often used for individual steps in the analysis, and these can be challenging to manage and integrate. However, in some instances, tools are combined into pipelines that are capable to complete all the essential steps to achieve the result. In the case of DNA methylation sequencing analysis, the goal of such pipeline is to map sequencing reads, calculate methylation levels, and distinguish differentially methylated positions and/or regions. The objective of this review is to describe basic principles and steps in the analysis of DNA methylation sequencing data that in particular have been used for mammalian genomes, and more importantly to present and discuss the most pronounced computational pipelines that can be used to analyze such data. We aim to provide a good starting point for scientists with limited experience in computational analyses of DNA methylation and hydroxymethylation data, and recommend a few tools that are powerful, but still easy enough to use for their own data analysis.

Introduction: Glutaminase inhibitors target cancer cells by blocking the conversion of glutamine to glutamate, thereby potentially interfering with anaplerosis and synthesis of amino acids and glutathione. The drug CB-839 has shown promising effects in preclinical experiments and is currently undergoing clinical trials in several human malignancies, including triple-negative breast cancer (TNBC). However, response to glutaminase inhibitors is variable and there is a need for identification of predictive response biomarkers. The aim of this study was to determine how glutamine is utilized in two patient-derived xenograft (PDX) models of breast cancer representing luminal-like/ER+ (MAS98.06) and basal-like/triple-negative (MAS98.12) breast cancer and to explore the metabolic effects of CB-839 treatment., Experimental: MAS98.06 and MAS98.12 PDX mice received CB-839 (200 mg/kg) or drug vehicle two times daily p.o. for up to 28 days (n = 5 per group), and the effect on tumor growth was evaluated. Expression of 60 genes and seven glutaminolysis key enzymes were determined using gene expression microarray analysis and immunohistochemistry (IHC), respectively, in untreated tumors. Uptake and conversion of glutamine were determined in the PDX models using HR MAS MRS after i.v. infusion of [5- 13 C] glutamine when the models had received CB-839 (200 mg/kg) or vehicle for 2 days (n = 5 per group)., Results: Tumor growth measurements showed that CB-839 significantly inhibited tumor growth in MAS98.06 tumors, but not in MAS98.12 tumors. Gene expression and IHC analysis indicated a higher proline synthesis from glutamine in untreated MAS98.06 tumors. This was confirmed by HR MAS MRS of untreated tumors demonstrating that MAS98.06 used glutamine to produce proline, glutamate, and alanine, and MAS98.12 to produce glutamate and lactate. In both models, treatment with CB-839 resulted in accumulation of glutamine. In addition, CB-839 caused depletion of alanine, proline, and glutamate ([1-13C] glutamate) in the MAS98.06 model., Conclusion: Our findings indicate that TNBCs may not be universally sensitive to glutaminase inhibitors. The major difference in the metabolic fate of glutamine between responding MAS98.06 xenografts and non-responding MAS98.12 xenografts is the utilization of glutamine for production of proline. We therefore suggest that addiction to proline synthesis from glutamine is associated with response to CB-839 in breast cancer. The effect of glutaminase inhibition in two breast cancer patient-derived xenograft (PDX) models. 13 C HR MAS MRS analysis of tumor tissue from CB-839-treated and untreated models receiving 13 C-labeled glutamine ([5- 13 C] Gln) shows that the glutaminase inhibitor CB-839 is causing an accumulation of glutamine (arrow up) in two PDX models representing luminal-like breast cancer (MAS98.06) and basal-like breast cancer (MAS98.12). In MAS98.06 tumors, CB-839 is in addition causing depletion of proline ([5- 13 C] Pro), alanine ([1- 13 C] Ala), and glutamate ([1- 13 C] Glu), which could explain why CB-839 causes tumor growth inhibition in MAS98.06 tumors, but not in MAS98.12 tumors.

The Norwegian e-Infrastructure for Life Sciences (NeLS) has been developed by ELIXIR Norway to provide its users with a system enabling data storage, sharing, and analysis in a project-oriented fashion. The system is available through easy-to-use web interfaces, including the Galaxy workbench for data analysis and workflow execution. Users confident with a command-line interface and programming may also access it through Secure Shell (SSH) and application programming interfaces (APIs).  NeLS has been in production since 2015, with training and support provided by the help desk of ELIXIR Norway. Through collaboration with NorSeq, the national consortium for high-throughput sequencing, an integrated service is offered so that sequencing data generated in a research project is provided to the involved researchers through NeLS. Sensitive data, such as individual genomic sequencing data, are handled using the TSD (Services for Sensitive Data) platform provided by Sigma2 and the University of Oslo. NeLS integrates national e-infrastructure storage and computing resources, and is also integrated with the SEEK platform in order to store large data files produced by experiments described in SEEK.   In this article, we outline the architecture of NeLS and discuss possible directions for further development., Competing Interests: No competing interests were disclosed.

Background: The relationship between cholesterol and prostate cancer has been extensively studied for decades, where high levels of cellular cholesterol are generally associated with cancer progression and less favorable outcomes. However, the role of in vivo cellular cholesterol synthesis in this process is unclear, and data on the transcriptional activity of cholesterol synthesis pathway genes in tissue from prostate cancer patients are inconsistent., Methods: A common problem with cancer tissue data from patient cohorts is the presence of heterogeneous tissue which confounds molecular analysis of the samples. In this study we present a general method to minimize systematic confounding from stroma tissue in any prostate cancer cohort comparing prostate cancer and normal samples. In particular we use samples assessed by histopathology to identify genes enriched and depleted in prostate stroma. These genes are then used to assess stroma content in tissue samples from other prostate cancer cohorts where no histopathology is available. Differential expression analysis is performed by comparing cancer and normal samples where the average stroma content has been balanced between the sample groups. In total we analyzed seven patient cohorts with prostate cancer consisting of 1713 prostate cancer and 230 normal tissue samples., Results: When stroma confounding was minimized, differential gene expression analysis over all cohorts showed robust and consistent downregulation of nearly all genes in the cholesterol synthesis pathway. Additional Gene Ontology analysis also identified cholesterol synthesis as the most significantly altered metabolic pathway in prostate cancer at the transcriptional level., Conclusion: The surprising observation that cholesterol synthesis genes are downregulated in prostate cancer is important for our understanding of how prostate cancer cells regulate cholesterol levels in vivo. Moreover, we show that tissue heterogeneity explains the lack of consistency in previous expression analysis of cholesterol synthesis genes in prostate cancer.

Docetaxel is the chemotherapeutic choice for metastatic hormone-refractory prostate cancer, however, it only marginally improves the survival rate. The purpose of the present study was to examine if a peptide targeting the cellular scaffold protein PCNA could improve docetaxel's efficacy. We found that docetaxel given in combination with a cell penetrating peptide containing the AlkB homolog 2 PCNA interacting motif (APIM-peptide), reduced the prostate volume and limited prostate cancer regrowth in vivo in the immunocompetent transgenic adenocarcinoma model of prostate cancer (TRAMP). In accordance with this, we found that the APIM-peptide enhanced the efficacy of docetaxel in vitro . Gene expression analysis on prostate cancer cell lines indicated that the combination of docetaxel and APIM-peptide alters expression of genes involved in cellular signaling, apoptosis, and prostate cancer development. These changes were not detected in single agent treated cells. Our results suggest that targeting PCNA and thereby affecting multiple cellular pathways simultaneously has the potential to improve docetaxel therapy of advanced prostate cancer., Competing Interests: CONFLICTS OF INTEREST APIM Therapeutics is a spin-off company of the Norwegian University of Science and Technology, and has co-funded this study. Professor Marit Otterlei is an inventor, minority shareholder and CSO of this company. Patent application no: PCT/GB2009/000489 “New PCNA interacting motif”, filed on February 20, 2009. There are no further patents, products or development or marked products to declare. The other authors declare no conflicts of interest.

Increased knowledge of the molecular differences between indolent and aggressive prostate cancer is needed for improved risk stratification and treatment selection. Secreted frizzled-related protein 4 (SFRP4) is a modulator of the cancer-associated Wnt pathway, and previously suggested as a potential marker for prostate cancer aggressiveness. In this study, we investigated and validated the association between SFRP4 gene expression and aggressiveness in nine independent cohorts (n = 2157). By differential expression and combined meta-analysis of all cohorts, we detected significantly higher SFRP4 expression in cancer compared with normal samples, and in high (3-5) compared with low (1-2) Grade Group samples. SFRP4 expression was a significant predictor of biochemical recurrence in six of seven cohorts and in the overall analysis, and was a significant predictor of metastatic event in one cohort. In our study cohort, where metabolic information was available, SFRP4 expression correlated significantly with the concentrations of citrate and spermine, two previously suggested biomarkers for aggressive prostate cancer. SFRP4 immunohistochemistry in an independent cohort (n = 33) was not associated with aggressiveness. To conclude, high SFRP4 gene expression is associated with high Grade Group and recurrent prostate cancer after surgery. Future studies investigating the mechanistic and clinical usefulness of SFRP4 in prostate cancer are warranted.

Activation of the Canonical Wnt pathway (CWP) has been linked to advanced and metastatic prostate cancer, whereas the Wnt5a-induced non-canonical Wnt pathway (NCWP) has been associated with both good and poor prognosis. A newly discovered NCWP, Wnt5/Fzd2, has been shown to induce epithelial-to-mesenchymal transition (EMT) in cancers, but has not been investigated in prostate cancer. The aim of this study was to investigate if the CWP and NCWP, in combination with EMT, are associated with metabolic alterations, aggressive disease and biochemical recurrence in prostate cancer. An initial analysis was performed using integrated transcriptomics, ex vivo and in vivo metabolomics, and histopathology of prostatectomy samples (n=129), combined with at least five-year follow-up. This analysis detected increased activation of NCWP through Wnt5a/ Fzd2 as the most common mode of Wnt activation in prostate cancer. This activation was associated with increased expression of EMT markers and higher Gleason score. The transcriptional association between NCWP and EMT was confirmed in five other publicly available patient cohorts (1519 samples in total). A novel gene expression signature of concordant activation of NCWP and EMT (NCWP-EMT) was developed, and this signature was significantly associated with metastasis and shown to be a significant predictor of biochemical recurrence. The NCWP-EMT signature was also associated with decreased concentrations of the metabolites citrate and spermine, which have previously been linked to aggressive prostate cancer. Our results demonstrate the importance of NCWP and EMT in prostate cancer aggressiveness, suggest a novel gene signature for improved risk stratification, and give new molecular insight.

TMPRSS2-ERG has been proposed to be a prognostic marker for prostate cancer. The aim of this study was to identify changes in metabolism, genes and biochemical recurrence related to TMPRSS2-ERG by using an integrated approach, combining metabolomics, transcriptomics, histopathology and clinical data in a cohort of 129 human prostate samples (41 patients). Metabolic analyses revealed lower concentrations of citrate and spermine comparing ERGhigh to ERGlow samples, suggesting an increased cancer aggressiveness of ERGhigh compared to ERGlow. These results could be validated in a separate cohort, consisting of 40 samples (40 patients), and magnetic resonance spectroscopy imaging (MRSI) indicated an in vivo translational potential. Alterations of gene expression levels associated with key enzymes in the metabolism of citrate and polyamines were in consistence with the metabolic results. Furthermore, the metabolic alterations between ERGhigh and ERGlow were more pronounced in low Gleason samples than in high Gleason samples, suggesting it as a potential tool for risk stratification. However, no significant difference in biochemical recurrence was detected, although a trend towards significance was detected for low Gleason samples. Using an integrated approach, this study suggests TMPRSS2-ERG as a potential risk stratification tool for inclusion of active surveillance patients., Competing Interests: None.

Molecular analysis of patient tissue samples is essential to characterize the in vivo variability in human cancers which are not accessible in cell-lines or animal models. This applies particularly to studies of tumor metabolism. The challenge is, however, the complex mixture of various tissue types within each sample, such as benign epithelium, stroma and cancer tissue, which can introduce systematic biases when cancers are compared to normal samples. In this study we apply a simple strategy to remove such biases using sample selections where the average content of stroma tissue is balanced between the sample groups. The strategy is applied to a prostate cancer patient cohort where data from MR spectroscopy and gene expression have been collected from and integrated on the exact same tissue samples. We reveal in vivo changes in cancer-relevant metabolic pathways which are otherwise hidden in the data due to tissue confounding. In particular, lowered levels of putrescine are connected to increased expression of SRM, reduced levels of citrate are attributed to upregulation of genes promoting fatty acid synthesis, and increased succinate levels coincide with reduced expression of SUCLA2 and SDHD. In addition, the strategy also highlights important metabolic differences between the stroma, epithelium and prostate cancer. These results show that important in vivo metabolic features of cancer can be revealed from patient data only if the heterogeneous tissue composition is properly accounted for in the analysis.

Background: PRL-3 is a phosphatase implicated in oncogenesis in multiple cancers. In some cancers, notably carcinomas, PRL-3 is also associated with inferior prognosis and increased metastatic potential. In this study we investigated the expression of PRL-3 mRNA in fresh-frozen samples from patients undergoing radical prostatectomy because of prostate cancer (PC) and the biological function of PRL-3 in prostate cancer cells., Methods: Samples from 41 radical prostatectomy specimens (168 samples in total) divided into low (Gleason score ≤ 6), intermediate (Gleason score = 7) and high (Gleason score ≥ 8) risk were analyzed with gene expression profiling and compared to normal prostate tissue. PRL-3 was identified as a gene with differential expression between healthy and cancerous tissue in these analyses. We used the prostate cancer cell lines PC3 and DU145 and a small molecular inhibitor of PRL-3 to investigate whether PRL-3 had a functional role in cancer. Relative ATP-measurement and thymidine incorporation were used to assess the effect of PRL-3 on growth of the cancer cells. We performed an in vitro scratch assay to investigate the involvement of PRL-3 in migration. Immunohistochemistry was used to identify PRL-3 protein in prostate cancer primary tumor and corresponding lymph node metastases., Results: Compared to normal prostate tissue, the prostate cancer tissue expressed a significantly higher level of PRL-3. We found PRL-3 to be present in both PC3 and DU145, and that inhibition of PRL-3 led to growth arrest and apoptosis in these two cell lines. Inhibition of PRL-3 led to reduced migration of the PC3 cells. Immunohistochemistry showed PRL-3 expression in both primary tumor and corresponding lymph node metastases., Conclusions: PRL-3 mRNA was expressed to a greater extent in prostate cancer tissue compared to normal prostate tissue. PRL-3 protein was expressed in both prostate cancer primary tumor and corresponding lymph node metastases. The results from our in vitro assays suggest that PRL-3 promotes growth and migration in prostate cancer. In conclusion, these results imply that PRL-3 has a role in the pathogenesis of prostate cancer.

Background: Good prognostic tools for predicting disease progression in early stage prostate cancer (PCa) are still missing. Detection of molecular subtypes, for instance by using microarray gene technology, can give new prognostic information which can assist personalized treatment planning. The detection of new subtypes with validation across additional and larger patient cohorts is important for bringing a potential prognostic tool into the clinic., Methods: We used fresh frozen prostatectomy tissue of high molecular quality to further explore four molecular subtype signatures of PCa based on Gene Set Enrichment Analysis (GSEA) of 15 selected gene sets published in a previous study. For this analysis we used a statistical test of dependent correlations to compare reference signatures to signatures in new normal and PCa samples, and also explore signatures within and between sample subgroups in the new samples., Results: An important finding was the consistent signatures observed for samples from the same patient independent of Gleason score. This proves that the signatures are robust and can surpass a normally high tumor heterogeneity within each patient. Our data did not distinguish between four different subtypes of PCa as previously published, but rather highlighted two groups of samples which could be related to good and poor prognosis based on survival data from the previous study.The poor prognosis group highlighted a set of samples characterized by enrichment of ESC, ERG-fusion and MYC + rich signatures in patients diagnosed with low Gleason score,. The other group consisted of PCa samples showing good prognosis as well as normal samples. Accounting for sample composition (the amount of benign structures such as stroma and epithelial cells in addition to the cancer component) was important to improve subtype assignments and should also be considered in future studies., Conclusion: Our study validates a previous molecular subtyping of PCa in a new patient cohort, and identifies a subgroup of PCa samples highly interesting for detecting high risk PCa at an early stage. The importance of taking sample tissue composition into account when assigning subtype is emphasized.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

The immense increase in availability of genomic scale datasets, such as those provided by the ENCODE and Roadmap Epigenomics projects, presents unprecedented opportunities for individual researchers to pose novel falsifiable biological questions. With this opportunity, however, researchers are faced with the challenge of how to best analyze and interpret their genome-scale datasets. A powerful way of representing genome-scale data is as feature-specific coordinates relative to reference genome assemblies, i.e. as genomic tracks. The Genomic HyperBrowser (http://hyperbrowser.uio.no) is an open-ended web server for the analysis of genomic track data. Through the provision of several highly customizable components for processing and statistical analysis of genomic tracks, the HyperBrowser opens for a range of genomic investigations, related to, e.g., gene regulation, disease association or epigenetic modifications of the genome.

Background: Chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-Seq) is the most frequently used method to identify the binding sites of transcription factors. Active binding sites can be seen as peaks in enrichment profiles when the sequencing reads are mapped to a reference genome. However, the profiles are normally noisy, making it challenging to identify all significantly enriched regions in a reliable way and with an acceptable false discovery rate., Results: We present the Triform algorithm, an improved approach to automatic peak finding in ChIP-Seq enrichment profiles for transcription factors. The method uses model-free statistics to identify peak-like distributions of sequencing reads, taking advantage of improved peak definition in combination with known characteristics of ChIP-Seq data., Conclusions: Triform outperforms several existing methods in the identification of representative peak profiles in curated benchmark data sets. We also show that Triform in many cases is able to identify peaks that are more consistent with biological function, compared with other methods. Finally, we show that Triform can be used to generate novel information on transcription factor binding in repeat regions, which represents a particular challenge in many ChIP-Seq experiments. The Triform algorithm has been implemented in R, and is available via http://tare.medisin.ntnu.no/triform.

Background: Transcription factors are important controllers of gene expression and mapping transcription factor binding sites (TFBS) is key to inferring transcription factor regulatory networks. Several methods for predicting TFBS exist, but there are no standard genome-wide datasets on which to assess the performance of these prediction methods. Also, it is believed that information about sequence conservation across different genomes can generally improve accuracy of motif-based predictors, but it is not clear under what circumstances use of conservation is most beneficial., Results: Here we use published ChIP-seq data and an improved peak detection method to create comprehensive benchmark datasets for prediction methods which use known descriptors or binding motifs to detect TFBS in genomic sequences. We use this benchmark to assess the performance of five different prediction methods and find that the methods that use information about sequence conservation generally perform better than simpler motif-scanning methods. The difference is greater on high-affinity peaks and when using short and information-poor motifs. However, if the motifs are specific and information-rich, we find that simple motif-scanning methods can perform better than conservation-based methods., Conclusions: Our benchmark provides a comprehensive test that can be used to rank the relative performance of transcription factor binding site prediction methods. Moreover, our results show that, contrary to previous reports, sequence conservation is better suited for predicting strong than weak transcription factor binding sites.

Chromatin immunoprecipitation (ChIP) followed by high throughput sequencing (ChIP-seq) is rapidly becoming the method of choice for discovering cell-specific transcription factor binding locations genome wide. By aligning sequenced tags to the genome, binding locations appear as peaks in the tag profile. Several programs have been designed to identify such peaks, but program evaluation has been difficult due to the lack of benchmark data sets. We have created benchmark data sets for three transcription factors by manually evaluating a selection of potential binding regions that cover typical variation in peak size and appearance. Performance of five programs on this benchmark showed, first, that external control or background data was essential to limit the number of false positive peaks from the programs. However, >80% of these peaks could be manually filtered out by visual inspection alone, without using additional background data, showing that peak shape information is not fully exploited in the evaluated programs. Second, none of the programs returned peak-regions that corresponded to the actual resolution in ChIP-seq data. Our results showed that ChIP-seq peaks should be narrowed down to 100-400 bp, which is sufficient to identify unique peaks and binding sites. Based on these results, we propose a meta-approach that gives improved peak definitions.

Image segmentation plays an important role in the automatic analysis of protein spots in 2-DE. Using image segments representing protein spots, the amount of protein in each segment can be quantified, and corresponding segments can be matched and compared for multiple gels. However, the common presence of image segments caused by noise and unwanted artefacts highly disturb the analysis and comparison of the gels. Common sources of such artefacts are cracks in the gel surface, fingerprints, dust and other pollutions. It would be advantageous to remove these unwanted artefacts during or after the segmentation procedure. To achieve this task a multivariate spot filtering model is developed using image segments from a gel segmentation. Parameters in the model are based on texture, shape and intensity measurements of the image segments. The model successfully managed to separate segments caused by noise, artefacts and cracks from image segments representing true protein spots. The classification method used is discriminant partial least squares regression.

The benefits of defining common spot boundaries when several gels from 2-DE are compared and analyzed have lately been stressed by both commercial software producers and users of this software. Though the importance of common spot boundaries is clearly stated, few reports exist that target this issue explicitly. In this study a method for defining common spots boundaries is developed, called the spot density method. The method consists of the following steps: segmentation and spot identification on each individual gel, transferring the spot-center coordinates for all gels onto a single new gel, collecting spot centers clustered together in the new gel and finally assigning pixels and new spot boundaries based on the spots in each cluster. The method is compared to a synthetic gel approach, and validated by visual inspection of three representative areas in the gels. The gel images need to be aligned prior to segmentation and spot identification, but the method can be used regardless of the choice of segmentation procedure. This makes the method an easy extension to existing methods for spot identification and matching. Conclusions based on the visual inspection are that the spot density method identifies partly overlapping spots and low-intensity spots better than the synthetic gel approach.

An improved pixel-based approach for analyzing 2-DE images is presented. The key feature of the method is to create a mask based on all gels in the experiment using image morphology, followed by multivariate analysis on the pixel level. The method reduces the impact of noise and background by identifying regions in the image where protein spots are present, but make no assumption on individual spot boundaries for isolated spots. This makes it possible to detect significant changes in complex regions, and visualize these changes over multiple gels in an easy way. False missing values and spot volumes caused by imposing erroneous spot boundaries are thus circumvented. The approach presented gives improved pixel-based information from the gels, and is also an alternative to existing methods for data-reduction, significance testing and visualization of 2-DE data. Results are compared with software using a common spot boundary approach on an experiment consisting of 35 full size gel images. Gel alignment is required before analysis.