Morten Beck Rye

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–400u2009bp, which is sufficient to identify unique peaks and binding sites. Based on these results, we propose a meta-approach that gives improved peak definitions.

Identification of serum microRNA profiles in colon cancer.

Background:microRNAs (miRNAs) exist in blood in an apparently stable form. We have explored whether serum miRNAs can be used as non-invasive early biomarkers of colon cancer.Methods:Serum samples from 30 patients with colon cancer stage IV and 10 healthy controls were examined for the expression of 375 cancer-relevant miRNAs. Based on the miRNA profile in this study, 34 selected miRNAs were measured in serum from 40 patients with stage I–II colon cancer and from 10 additional controls.Results:Twenty miRNAs were differentially expressed in serum from stage IV patients compared with controls (P<0.01). Unsupervised clustering revealed four subgroups; one corresponding mostly to the control group and the three others to the patient groups. Of the 34 miRNAs measured in the follow-up study of stage I–II patients, 21 showed concordant expression between stage IV and stage I–II patient. Based on the profiles of these 21 miRNAs, a supervised linear regression analysis (Partial Least Squares Regression) was performed. Using this model we correctly assigned stage I–II colon cancer patients based on miRNA profiles of stage IV patients.Conclusion:Serum miRNA expression profiling may be utilised in early detection of colon 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.

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

BackgroundTranscription 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.ResultsHerein 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.ConclusionBy 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.

An improved pixel-based approach for analyzing images in two-dimensional gel electrophoresis.

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.

A ChIP-Seq Benchmark Shows That Sequence Conservation Mainly Improves Detection of Strong Transcription Factor Binding Sites

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.

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

BackgroundContext-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?ResultsWe 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.ConclusionsTogether, 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.

A multivariate spot filtering model for two-dimensional gel electrophoresis.

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.

Gene signatures ESC, MYC and ERG-fusion are early markers of a potentially dangerous subtype of prostate cancer

BackgroundGood 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.MethodsWe 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.ResultsAn 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 MYCu2009+u2009rich 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.ConclusionOur 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.

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

BackgroundChromatin 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.ResultsWe 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.ConclusionsTriform 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.