PubMed Health⌕ Search

Biomedical subjects

Itamar Simon

Publications and source records attributed to Itamar Simon.

14 recordsLinked to original sources

Improved spike-in normalization clarifies the relationship between active histone modifications and transcription.

Spike-in normalization enables quantitative analysis of chromatin immunoprecipitation sequencing (ChIP-seq) signal. Here we introduce a robust dual spike-in normalization approach for ChIP-seq (ChIP-wrangler), optimize parameters and verify its accuracy in quantifying changes in ChIP-seq signal and detecting technical artifacts. We use ChIP-wrangler to revisit recent claims that active histone marks depend on transcription. We show that acute depletion of RNA polymerase II (RNAPII) has a modest impact on H3K27ac levels, with only 6% of peaks significantly changing after RNAPII depletion, indicating that histone acetylation maintenance is not entirely dependent on ongoing transcription. Promoters and enhancers are differentially affected, with 82% of decreasing acetylation peaks located at promoter-distal elements with enhancer-related motifs. ChIP-wrangler provides increased rigor and 'guardrails' for successful spike-in normalization and, as applied here, refines the understanding of crosstalk between RNAPII activity and transcription-associated histone marks.

Histones↗

Reconstructing dynamic regulatory maps.

Even simple organisms have the ability to respond to internal and external stimuli. This response is carried out by a dynamic network of protein-DNA interactions that allows the specific regulation of genes needed for the response. We have developed a novel computational method that uses an input-output hidden Markov model to model these regulatory networks while taking into account their dynamic nature. Our method works by identifying bifurcation points, places in the time series where the expression of a subset of genes diverges from the rest of the genes. These points are annotated with the transcription factors regulating these transitions resulting in a unified temporal map. Applying our method to study yeast response to stress, we derive dynamic models that are able to recover many of the known aspects of these responses. Predictions made by our method have been experimentally validated leading to new roles for Ino4 and Gcn4 in controlling yeast response to stress. The temporal cascade of factors reveals common pathways and highlights differences between master and secondary factors in the utilization of network motifs and in condition-specific regulation.

Algorithms↗

Polycomb-mediated methylation on Lys27 of histone H3 pre-marks genes for de novo methylation in cancer.

Many genes associated with CpG islands undergo de novo methylation in cancer. Studies have suggested that the pattern of this modification may be partially determined by an instructive mechanism that recognizes specifically marked regions of the genome. Using chromatin immunoprecipitation analysis, here we show that genes methylated in cancer cells are specifically packaged with nucleosomes containing histone H3 trimethylated on Lys27. This chromatin mark is established on these unmethylated CpG island genes early in development and then maintained in differentiated cell types by the presence of an EZH2-containing Polycomb complex. In cancer cells, as opposed to normal cells, the presence of this complex brings about the recruitment of DNA methyl transferases, leading to de novo methylation. These results suggest that tumor-specific targeting of de novo methylation is pre-programmed by an established epigenetic system that normally has a role in marking embryonic genes for repression.

Caco-2 Cells↗

An analysis of intra array repeats: the good, the bad and the non informative.

BACKGROUND: On most common microarray platforms many genes are represented by multiple probes. Although this is quite common no one has systematically explored the concordance between probes mapped to the same gene. RESULTS: Here we present an analysis of all the cases of multiple probe sets measuring the same gene on the Affymetrix U133a GeneChip and found that although in the majority of cases both measurements tend to agree there are a significant number of cases in which the two measurements differ from each other. In these cases the measurements can not be simply averaged but rather should be handled individually. CONCLUSION: Our analysis allows us to provide a comprehensive list of the correlation between all pairs of probe sets that are mapped to the same gene and thus allows microarray users to sort out the cases that deserve further analysis. Comparison between the set of highly correlated pairs and the set of pairs that tend to differ from each other reveals potential factors that may affect it.

Data Interpretation, Statistical↗

Evidence for an instructive mechanism of de novo methylation in cancer cells.

DNA methylation has a role in the regulation of gene expression during normal mammalian development but can also mediate epigenetic silencing of CpG island genes in cancer and other diseases. Many individual genes (including tumor suppressors) have been shown to undergo de novo methylation in specific tumor types, but the biological logic inherent in this process is not understood. To decipher this mechanism, we have adopted a new approach for detecting CpG island DNA methylation that can be used together with microarray technology. Genome-wide analysis by this technique demonstrated that tumor-specific methylated genes belong to distinct functional categories, have common sequence motifs in their promoters and are found in clusters on chromosomes. In addition, many are already repressed in normal cells. These results are consistent with the hypothesis that cancer-related de novo methylation may come about through an instructive mechanism.

Animals↗

MILANO--custom annotation of microarray results using automatic literature searches.

BACKGROUND: High-throughput genomic research tools are becoming standard in the biologist's toolbox. After processing the genomic data with one of the many available statistical algorithms to identify statistically significant genes, these genes need to be further analyzed for biological significance in light of all the existing knowledge. Literature mining--the process of representing literature data in a fashion that is easy to relate to genomic data--is one solution to this problem. RESULTS: We present a web-based tool, MILANO (Microarray Literature-based Annotation), that allows annotation of lists of genes derived from microarray results by user defined terms. Our annotation strategy is based on counting the number of literature co-occurrences of each gene on the list with a user defined term. This strategy allows the customization of the annotation procedure and thus overcomes one of the major limitations of the functional annotations usually provided with microarray results. MILANO expands the gene names to include all their informative synonyms while filtering out gene symbols that are likely to be less informative as literature searching terms. MILANO supports searching two literature databases: GeneRIF and Medline (through PubMed), allowing retrieval of both quick and comprehensive results. We demonstrate MILANO's ability to improve microarray analysis by analyzing a list of 150 genes that were affected by p53 overproduction. This analysis reveals that MILANO enables immediate identification of known p53 target genes on this list and assists in sorting the list into genes known to be involved in p53 related pathways, apoptosis and cell cycle arrest. CONCLUSIONS: MILANO provides a useful tool for the automatic custom annotation of microarray results which is based on all the available literature. MILANO has two major advances over similar tools: the ability to expand gene names to include all their informative synonyms while removing synonyms that are not informative and access to the GeneRIF database which provides short summaries of curated articles relevant to known genes. MILANO is available at http://milano.md.huji.ac.il.

Algorithms↗

Combined static and dynamic analysis for determining the quality of time-series expression profiles.

Expression profiling of time-series experiments is widely used to study biological systems. However, determining the quality of the resulting profiles remains a fundamental problem. Because of inadequate sampling rates, the effect of arrest-and-release methods and loss of synchronization, the measurements obtained from a series of time points may not accurately represent the underlying expression profiles. To solve this, we propose an approach that combines time-series and static (average) expression data analysis--for each gene, we determine whether its temporal expression profile can be reconciled with its static expression levels. We show that by combining synchronized and unsynchronized human cell cycle data, we can identify many cycling genes that are missed when using only time-series data. The algorithm also correctly distinguishes cycling genes from genes that specifically react to an environmental stimulus even if they share similar temporal expression profiles. Experimental validation of these results shows the utility of this analytical approach for determining the accuracy of gene expression patterns.

Algorithms↗

Deconvolving cell cycle expression data with complementary information.

MOTIVATION: In the study of many systems, cells are first synchronized so that a large population of cells exhibit similar behavior. While synchronization can usually be achieved for a short duration, after a while cells begin to lose their synchronization. Synchronization loss is a continuous process and so the observed value in a population of cells for a gene at time t is actually a convolution of its values in an interval around t. Deconvolving the observed values from a mixed population will allow us to obtain better models for these systems and to accurately detect the genes that participate in these systems. RESULTS: We present an algorithm which combines budding index and gene expression data to deconvolve expression profiles. Using the budding index data we first fit a synchronization loss model for the cell cycle system. Our deconvolution algorithm uses this loss model and can also use information from co-expressed genes, making it more robust against noise and missing values. Using expression and budding data for yeast we show that our algorithm is able to reconstruct a more accurate representation when compared with the observed values. In addition, using the deconvolved profiles we are able to correctly identify 15% more cycling genes when compared to a set identified using the observed values. AVAILABILITY: Matlab implementation can be downloaded from the supporting website http://www.cs.cmu.edu/~zivbj/decon/decon.html

Biological Clocks↗

Comparing the continuous representation of time-series expression profiles to identify differentially expressed genes.

We present a general algorithm to detect genes differentially expressed between two nonhomogeneous time-series data sets. As increasing amounts of high-throughput biological data become available, a major challenge in genomic and computational biology is to develop methods for comparing data from different experimental sources. Time-series whole-genome expression data are a particularly valuable source of information because they can describe an unfolding biological process such as the cell cycle or immune response. However, comparisons of time-series expression data sets are hindered by biological and experimental inconsistencies such as differences in sampling rate, variations in the timing of biological processes, and the lack of repeats. Our algorithm overcomes these difficulties by using a continuous representation for time-series data and combining a noise model for individual samples with a global difference measure. We introduce a corresponding statistical method for computing the significance of this differential expression measure. We used our algorithm to compare cell-cycle-dependent gene expression in wild-type and knockout yeast strains. Our algorithm identified a set of 56 differentially expressed genes, and these results were validated by using independent protein-DNA-binding data. Unlike previous methods, our algorithm was also able to identify 22 non-cell-cycle-regulated genes as differentially expressed. This set of genes is significantly correlated in a set of independent expression experiments, suggesting additional roles for the transcription factors Fkh1 and Fkh2 in controlling cellular activity in yeast.

Algorithms↗

Systematic analysis of essential yeast TAFs in genome-wide transcription and preinitiation complex assembly.

The general transcription factor TFIID is composed of the TATA box binding protein (TBP) and a set of conserved TBP-associated factors (TAFs). Here we report the completion of genome-wide expression profiling analyses of yeast strains bearing temperature-sensitive mutations in each of the 13 essential TAFs. The percentage of the yeast genome dependent on each TAF ranges from 3% (TAF2) to 59-61% (TAF9). Approximately 84% of yeast genes are dependent upon one or more TAFs and 16% of yeast genes are TAF independent. In addition, this complete analysis defines three distinct classes of yeast promoters whose transcriptional requirements for TAFs differ substantially. Using this collection of temperature-sensitive mutants, we show that in all cases the transcriptional dependence for a TAF can be explained by a requirement for TBP recruitment and assembly of the preinitiation complex (PIC). Unexpectedly, these assembly experiments reveal that TAF11 and TAF13 appear to provide the critical functional contacts with TBP during PIC assembly. Collectively, our results confirm and extend the proposal that individual TAFs have selective transcriptional roles and distinct functions.

Base Sequence↗

Program-specific distribution of a transcription factor dependent on partner transcription factor and MAPK signaling.

Specialized gene expression programs are induced by signaling pathways that act on transcription factors. Whether these transcription factors can function in multiple developmental programs through a global switch in promoter selection is not known. We have used genome-wide location analysis to show that the yeast Ste12 transcription factor, which regulates mating and filamentous growth, is bound to distinct program-specific target genes dependent on the developmental condition. This condition-dependent distribution of Ste12 requires concurrent binding of the transcription factor Tec1 during filamentation and is differentially regulated by the MAP kinases Fus3 and Kss1. Program-specific distribution across the genome may be a general mechanism by which transcription factors regulate distinct gene expression programs in response to signaling.

DNA-Binding Proteins↗

Continuous representations of time-series gene expression data.

We present algorithms for time-series gene expression analysis that permit the principled estimation of unobserved time points, clustering, and dataset alignment. Each expression profile is modeled as a cubic spline (piecewise polynomial) that is estimated from the observed data and every time point influences the overall smooth expression curve. We constrain the spline coefficients of genes in the same class to have similar expression patterns, while also allowing for gene specific parameters. We show that unobserved time points can be reconstructed using our method with 10-15% less error when compared to previous best methods. Our clustering algorithm operates directly on the continuous representations of gene expression profiles, and we demonstrate that this is particularly effective when applied to nonuniformly sampled data. Our continuous alignment algorithm also avoids difficulties encountered by discrete approaches. In particular, our method allows for control of the number of degrees of freedom of the warp through the specification of parameterized functions, which helps to avoid overfitting. We demonstrate that our algorithm produces stable low-error alignments on real expression data and further show a specific application to yeast knock-out data that produces biologically meaningful results.

Algorithms↗

Transcriptional regulatory networks in Saccharomyces cerevisiae.

We have determined how most of the transcriptional regulators encoded in the eukaryote Saccharomyces cerevisiae associate with genes across the genome in living cells. Just as maps of metabolic networks describe the potential pathways that may be used by a cell to accomplish metabolic processes, this network of regulator-gene interactions describes potential pathways yeast cells can use to regulate global gene expression programs. We use this information to identify network motifs, the simplest units of network architecture, and demonstrate that an automated process can use motifs to assemble a transcriptional regulatory network structure. Our results reveal that eukaryotic cellular functions are highly connected through networks of transcriptional regulators that regulate other transcriptional regulators.

Algorithms↗

The genome-wide localization of Rsc9, a component of the RSC chromatin-remodeling complex, changes in response to stress.

The cellular response to environmental changes includes widespread modifications in gene expression. Here we report the identification and characterization of Rsc9, a member of the RSC chromatin-remodeling complex in yeast. The genome-wide localization of Rsc9 indicated a relationship between genes targeted by Rsc9 and genes regulated by stress; treatment with hydrogen peroxide or rapamycin, which inhibits TOR signaling, resulted in genome-wide changes in Rsc9 occupancy. We further show that Rsc9 is involved in both repression and activation of mRNAs regulated by TOR as well as the synthesis of rRNA. Our results illustrate the response of a chromatin-remodeling factor to signaling cascades and suggest that changes in the activity of chromatin-remodeling factors are reflected in changes in their localization in the genome.

Amino Acid Sequence↗