Prediction of recurrence from metabolites and expression of TOP2A and EZH2 in prostate cancer patients treated with radiotherapy.

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.

DNA methylation changes in response to neoadjuvant chemotherapy are associated with breast cancer survival.

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.

MACPET: model-based analysis for ChIA-PET.

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.

Differential Effects of Dietary Supplementation of Krill Meal, Soybean Meal, Butyrate, and Bactocell® on the Gene Expression of Atlantic Salmon Head Kidney.

The head kidney is a key organ that plays a fundamental role in the regulation of the fish immune response and in the maintenance of endocrine homeostasis. Previous studies indicate that the supplementation of exogenous dietary components, such as krill meal (KM), soybean meal (SM), Bactocell® (BA), and butyrate (BU), can have a significant effect on the immune function of the head kidney. The aim of this study was to investigate the differential effect of these four dietary ingredients on the transcriptional profiles of the head kidney of the Atlantic salmon. This study revealed that just a small number of genes were responsive to the feeding regime after a long-term (12 weeks) treatment, and evidenced that the most significant alterations, both in terms of the number of affected genes and magnitude of changes in gene expression, were detectable in the BU- and KM-fed groups compared with controls, while the SM diet had a nearly negligible effect, and BA had no significant effects at all. Most of the differentially expressed genes were involved in the immune response and, in line with data previously obtained from pyloric caeca, major components of the complement system were significantly affected. These alterations were accompanied by an increase in the density of melanomacrophage centers in the KM- and SM-fed group and their reduction in the BU-fed group. While three types of dietary supplements (BU, KM, and SM) were able to produce a significant modulation of some molecular players of the immune system, the butyrate-rich diet was revealed as the one with the most relevant immune-stimulating properties in the head kidney. These preliminary results suggest that further investigations should be aimed towards the elucidation of the potential beneficial effects of butyrate and krill meal supplementation on farmed salmon health and growth performance.

Glutamine to proline conversion is associated with response to glutaminase inhibition in breast cancer.

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-13C] 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. 13C HR MAS MRS analysis of tumor tissue from CB-839-treated and untreated models receiving 13C-labeled glutamine ([5-13C] 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-13C] Pro), alanine ([1-13C] Ala), and glutamate ([1-13C] Glu), which could explain why CB-839 causes tumor growth inhibition in MAS98.06 tumors, but not in MAS98.12 tumors.

Selection response for Streptococcus agalactiae resistance in Nile tilapia Oreochromis niloticus.

The potential of selection to improve resistance to streptococcosis was evaluated in a commercial population of Nile tilapia in Thailand. The base generation (G0) consisted of offspring from 98 sires and 149 dams using a partly nested design. At 60 days post-hatch, 30 fish from each family were injected intraperitoneally with a Streptococcosis agalactiae solution (1 × 109  CFU/ml) and evaluated for 14 days. Disease resistance was recorded as the number of days from challenge until death (DD) and as a binary (BIN) trait (dead/alive) on day 14. Three models were used for genetic analyses: Cox frailty model for DD; animal model for DD; and animal model for BIN. Age at challenge was fitted as a covariate and contemporary group as fixed or random effect, depending on the model. Fish from the 18 most resistant families were selected to produce the first generation (G1). Heritability estimates for G0 were 0.22, 0.14 ± 0.02 and 0.11 ± 0.02 for the Cox, linear DD and linear BIN models, respectively. Selection response indicated that the risk of death decreased to 54%, survival time increased to 3.4 days and survival rate increased to 21%. These results suggest that genetic improvement is possible for this population.

Cholesterol synthesis pathway genes in prostate cancer are transcriptionally downregulated when tissue confounding is minimized.

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.

Phosphatase of regenerating liver 3 (PRL-3) is overexpressed in human prostate cancer tissue and promotes growth and migration.

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.

The Genomic HyperBrowser: an analysis web server for genome-scale data.

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.

Chromatin states reveal functional associations for globally defined transcription start sites in four human cell lines.

BACKGROUND: Deciphering the most common modes by which chromatin regulates transcription, and how this is related to cellular status and processes is an important task for improving our understanding of human cellular biology. The FANTOM5 and ENCODE projects represent two independent large scale efforts to map regulatory and transcriptional features to the human genome. Here we investigate chromatin features around a comprehensive set of transcription start sites in four cell lines by integrating data from these two projects. RESULTS: Transcription start sites can be distinguished by chromatin states defined by specific combinations of both chromatin mark enrichment and the profile shapes of these chromatin marks. The observed patterns can be associated with cellular functions and processes, and they also show association with expression level, location relative to nearby genes, and CpG content. In particular we find a substantial number of repressed inter- and intra-genic transcription start sites enriched for active chromatin marks and Pol II, and these sites are strongly associated with immediate-early response processes and cell signaling. Associations between start sites with similar chromatin patterns are validated by significant correlations in their global expression profiles. CONCLUSIONS: The results confirm the link between chromatin state and cellular function for expressed transcripts, and also indicate that active chromatin states at repressed transcripts may poise transcripts for rapid activation during immune response.

SOS genes are rapidly induced while translesion synthesis polymerase activity is temporally regulated.

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.

Survival, growth and sexual maturation in Atlantic salmon exposed to infectious pancreatic necrosis: a multi-variate mixture model approach.

BACKGROUND: Outbreaks of infectious pancreatic necrosis (IPN) in Atlantic salmon can result in reduced growth rates in a fraction of the surviving fish (runts). Genetic and environmental variation also affects growth rates within different categories of healthy animals and runts, which complicates identification of runts. Mixture models are commonly used to identify the underlying structures in such data, and the aim of this study was to develop Bayesian mixture models for the genetic analysis of health status (runt/healthy) of surviving fish from an IPN outbreak. METHODS: Five statistical models were tested on data consisting of 10 972 fish that died and 3959 survivors with recorded growth data. The most complex models (4 and 5) were multivariate normal-binary mixture models including growth, sexual maturity and field survival traits. Growth rate and liability of sexual maturation were treated as two-component normal mixtures, assuming phenotypes originated from two potentially overlapping distributions, (runt/normal). Runt status was an unobserved binary trait. These models were compared to mixture models with fewer traits (Models 2 and 3) and a classical linear animal model for growth (Model 1). RESULTS: Assuming growth as a mixture trait improved the predictive ability of the statistical model considerably (Model 2 vs. 1). The final models (4 and 5) yielded the following results: estimated (underlying) heritabilities were moderate for growth in healthy fish (0.32 ± 0.04 and 0.35 ± 0.05), runt status (0.39 ± 0.07 and 0.36 ± 0.08) and sexual maturation (0.33 ± 0.05), and high for field survival (0.47 ± 0.03 and 0.48 ± 0.03). Growth in healthy animals, runt status and survival showed consistent favourable genetic associations. Sexual maturation showed an unfavourable non-significant genetic correlation with runt status, but favourable genetic correlations with other traits. The estimated fraction of healthy fish was 81-85%. The estimated breeding values for runt status and (normal) growth were consistent for the most complex models (4 and 5), but showed imperfect correlations with estimated breeding values from the simpler models. CONCLUSIONS: Modelling growth in IPN survivors as a mixture trait improved the predictive ability of the model compared with a classical linear model. The results indicated considerable genetic variation in health status among survivors. Mixture modelling may be useful for the genetic analysis of diseases detected mainly through indicator traits.

A manually curated ChIP-seq benchmark demonstrates room for improvement in current peak-finder programs.

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.

Genome-wide association study for Streptococcus iniae in Nile tilapia (Oreochromis niloticus) identifies a significant QTL for disease resistance.

Streptococcus iniae is a problematic gram-positive bacterium negatively affecting Nile tilapia (Oreochromis niloticus), one of the main aquacultural species produced worldwide. The aim of this study was to identify the genetic architecture of survival to S. iniae and identify single nucleotide polymorphism (SNPs) linked to quantitative trait loci (QTL) related to survival to S. iniae challenge. With this purpose, Nile tilapia from the Spring Genetics breeding program were sent to a controlled S. iniae challenge test where phenotypes were scored as dead for fish that died during challenge test and survivors for the fish alive at the termination of the test. Additionally, fin-clip samples from all fish in the test were collected for DNA extraction. Out of 1904 fish in the challenge test, tissue samples of 321 fish were sent for genotyping using double digest restriction site associated DNA sequencing (ddRADseq). After quality control and filtering, 9,085 SNPs were used to perform a genome-wide association study (GWAS). A significant signal in LG8 was observed indicating association with survival to S. iniae challenge, with SNPs explaining from 12% to 26% of the genetic variance. To demonstrate the usefulness of marker assisted selection (MAS) to selectively breed fish for survival to S. iniae, offspring of breeding candidates classified as "resistant" and "susceptible" based on haplotypes of the four most significant markers were sent to a controlled S. iniae challenge test. At the end of the test, the differences in mortality between the two groups were strikingly different with a final cumulative percent mortality of less than 1% and 73% for offspring from "resistant" and "susceptible" parents, respectively. These results demonstrate that MAS for improved resistance to S. iniae is feasible.

Applying genetic technologies to combat infectious diseases in aquaculture.

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

The Triform algorithm: improved sensitivity and specificity in ChIP-Seq peak finding.

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.

Cell-type specificity of ChIP-predicted transcription factor binding sites.

BACKGROUND: Context-dependent transcription factor (TF) binding is one reason for differences in gene expression patterns between different cellular states. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identifies genome-wide TF binding sites for one particular context-the cells used in the experiment. But can such ChIP-seq data predict TF binding in other cellular contexts and is it possible to distinguish context-dependent from ubiquitous TF binding? RESULTS: We compared ChIP-seq data on TF binding for multiple TFs in two different cell types and found that on average only a third of ChIP-seq peak regions are common to both cell types. Expectedly, common peaks occur more frequently in certain genomic contexts, such as CpG-rich promoters, whereas chromatin differences characterize cell-type specific TF binding. We also find, however, that genotype differences between the cell types can explain differences in binding. Moreover, ChIP-seq signal intensity and peak clustering are the strongest predictors of common peaks. Compared with strong peaks located in regions containing peaks for multiple transcription factors, weak and isolated peaks are less common between the cell types and are less associated with data that indicate regulatory activity. CONCLUSIONS: Together, the results suggest that experimental noise is prevalent among weak peaks, whereas strong and clustered peaks represent high-confidence binding events that often occur in other cellular contexts. Nevertheless, 30-40% of the strongest and most clustered peaks show context-dependent regulation. We show that by combining signal intensity with additional data-ranging from context independent information such as binding site conservation and position weight matrix scores to context dependent chromatin structure-we can predict whether a ChIP-seq peak is likely to be present in other cellular contexts.

Clustered ChIP-Seq-defined transcription factor binding sites and histone modifications map distinct classes of regulatory elements.

BACKGROUND: Transcription factor binding to DNA requires both an appropriate binding element and suitably open chromatin, which together help to define regulatory elements within the genome. Current methods of identifying regulatory elements, such as promoters or enhancers, typically rely on sequence conservation, existing gene annotations or specific marks, such as histone modifications and p300 binding methods, each of which has its own biases. RESULTS: Herein we show that an approach based on clustering of transcription factor peaks from high-throughput sequencing coupled with chromatin immunoprecipitation (Chip-Seq) can be used to evaluate markers for regulatory elements. We used 67 data sets for 54 unique transcription factors distributed over two cell lines to create regulatory element clusters. By integrating the clusters from our approach with histone modifications and data for open chromatin, we identified general methylation of lysine 4 on histone H3 (H3K4me) as the most specific marker for transcription factor clusters. Clusters mapping to annotated genes showed distinct patterns in cluster composition related to gene expression and histone modifications. Clusters mapping to intergenic regions fall into two groups either directly involved in transcription, including miRNAs and long noncoding RNAs, or facilitating transcription by long-range interactions. The latter clusters were specifically enriched with H3K4me1, but less with acetylation of lysine 27 on histone 3 or p300 binding. CONCLUSION: By integrating genomewide data of transcription factor binding and chromatin structure and using our data-driven approach, we pinpointed the chromatin marks that best explain transcription factor association with different regulatory elements. Our results also indicate that a modest selection of transcription factors may be sufficient to map most regulatory elements in the human genome.

FunHoP analysis reveals upregulation of mitochondrial genes in prostate cancer.

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.

The genes controlling normal function of citrate and spermine secretion are lost in aggressive prostate cancer and prostate model systems.

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.