Genome-wide Location Data Research Articles

A genome-wide CRISPR-Cas9 knockout screen identifies essential and growth-restricting genes in human trophoblast stem cells

The recent derivation of human trophoblast stem cells (hTSCs) provides a scalable in vitro model system of human placental development, but the molecular regulators of hTSC identity have not been systematically explored thus far. Here, we utilize a genome-wide CRISPR-Cas9 knockout screen to comprehensively identify essential and growth-restricting genes in hTSCs. By cross-referencing our data to those from similar genetic screens performed in other cell types, as well as gene expression data from early human embryos, we define hTSC-specific and -enriched regulators. These include both well-established and previously uncharacterized trophoblast regulators, such as ARID3A, GATA2, and TEAD1 (essential), and GCM1, PTPN14, and TET2 (growth-restricting). Integrated analysis of chromatin accessibility, gene expression, and genome-wide location data reveals that the transcription factor TEAD1 regulates the expression of many trophoblast regulators in hTSCs. In the absence of TEAD1, hTSCs fail to complete faithful differentiation into extravillous trophoblast (EVT) cells and instead show a bias towards syncytiotrophoblast (STB) differentiation, thus indicating that this transcription factor safeguards the bipotent lineage potential of hTSCs. Overall, our study provides a valuable resource for dissecting the molecular regulation of human placental development and diseases.

BENIN: Biologically enhanced network inference.

Gene regulatory network inference is one of the central problems in computational biology. We need models that integrate the variety of data available in order to use their complementarity information to overcome the issues of noisy and limited data. BENIN: Biologically Enhanced Network INference is our proposal to integrate data and infer more accurate networks. BENIN is a general framework that jointly considers different types of prior knowledge with expression datasets to improve the network inference. The method states the network inference as a feature selection problem and uses a popular penalized regression method, the Elastic net, combined with bootstrap resampling to solve it. BENIN significantly outperforms the state-of-the-art methods on the simulated data from the DREAM 4 challenge when combining genome-wide location data, knockout gene expression data, and time series expression data.

Epigenomic Co-localization and Co-evolution Reveal a Key Role for 5hmC as a Communication Hub in the Chromatin Network of ESCs.

Epigenetic communication through histone and cytosine modifications is essential for gene regulation and cell identity. Here, we propose a framework that is based on a chromatin communication model to get insight on the function of epigenetic modifications in ESCs. The epigenetic communication network was inferred from genome-wide location data plus extensive manual annotation. Notably, we found that 5-hydroxymethylcytosine (5hmC) is the most-influential hub of this network, connecting DNA demethylation to nucleosome remodeling complexes and to key transcription factors of pluripotency. Moreover, an evolutionary analysis revealed a central role of 5hmC in the co-evolution of chromatin-related proteins. Further analysis of regions where 5hmC co-localizes with specific interactors shows that each interaction points to chromatin remodeling, stemness, differentiation, or metabolism. Our results highlight the importance of cytosine modifications in the epigenetic communication of ESCs.

Parameter Estimation for Gene Regulatory Networks from Microarray Data: Cold Shock Response in Saccharomyces cerevisiae.

We investigated the dynamics of a gene regulatory network controlling the cold shock response in budding yeast, Saccharomyces cerevisiae. The medium-scale network, derived from published genome-wide location data, consists of 21 transcription factors that regulate one another through 31 directed edges. The expression levels of the individual transcription factors were modeled using mass balance ordinary differential equations with a sigmoidal production function. Each equation includes a production rate, a degradation rate, weights that denote the magnitude and type of influence of the connected transcription factors (activation or repression), and a threshold of expression. The inverse problem of determining model parameters from observed data is our primary interest. We fit the differential equation model to published microarray data using a penalized nonlinear least squares approach. Model predictions fit the experimental data well, within the 95 % confidence interval. Tests of the model using randomized initial guesses and model-generated data also lend confidence to the fit. The results have revealed activation and repression relationships between the transcription factors. Sensitivity analysis indicates that the model is most sensitive to changes in the production rate parameters, weights, and thresholds of Yap1, Rox1, and Yap6, which form a densely connected core in the network. The modeling results newly suggest that Rap1, Fhl1, Msn4, Rph1, and Hsf1 play an important role in regulating the early response to cold shock in yeast. Our results demonstrate that estimation for a large number of parameters can be successfully performed for nonlinear dynamic gene regulatory networks using sparse, noisy microarray data.Electronic supplementary materialThe online version of this article (doi:10.1007/s11538-015-0092-6) contains supplementary material, which is available to authorized users.

Measurement variation determines the gene network topology reconstructed from experimental data: a case study of the yeast cyclin network

Inference of the topology of gene regulatory networks from experimental data is one of the primary challenges of systems biology. In an example of a genetic network of cyclins in the yeast cell cycle, we analyzed static genome-wide location data together with microarray kinetic measurements using a recurrent neural network-based model of gene expression and a newly developed, unbiased algorithm based on evolutionary programming principles. The modeling and simulation of gene expression dynamics identified cyclin genetic networks that were active during the cell cycle. We document that because there is inherent experimental variation, it is not possible to identify a single genetic network, only a set of equivalent networks with the same probability of occurrence. Analysis of these networks showed that each target gene was controlled by only a few regulators and that the control was robust. These results led to the reformulation of the cyclin genetic network in the yeast cell cycle as previously published. The analysis shows that with the methodologies that are currently available, it is not possible to predict only one genetic network; rather, we must work with the hypothesis of multiple, equivalent networks. Chromatin immunoprecipitation (ChIP)-on-chip experiments are not sufficient to predict the functional networks that are active during an investigated process. Such predictions must be considered as only potential, and their actual realization during particular cellular processes must be identified by incorporating both kinetic and other types of data.

Investigating transcriptional regulation: From analysis of complex networks to discovery of cis-regulatory elements

Regulation of gene expression at the transcriptional level is a fundamental mechanism that is well conserved in all cellular systems. Due to advances in large-scale experimental analyses, we now have a wealth of information on gene regulation such as mRNA expression level across multiple conditions, genome-wide location data of transcription factors and data on transcription factor binding sites. This knowledge can be used to reconstruct transcriptional regulatory networks. Such networks are usually represented as directed graphs where regulatory interactions are depicted as directed edges from the transcription factor nodes to the target gene nodes. This abstract representation allows us to apply graph theory to study transcriptional regulation at global and local levels, to predict regulatory motifs and regulatory modules such as regulons and to compare the regulatory network of different genomes. Here we review some of the available computational methodologies for studying transcriptional regulatory networks as well as their interpretation.

Integrative genomics analysis reveals silencing of β-adrenergic signaling by polycomb in prostate cancer: Yu J, Cao Q, Mehra R, Laxman B, Yu J, Tomlins SA, Creighton CJ, Dhanasekaran SM, Shen R, Chen G, Morris DS, Marquez VE, Shah RB, Ghosh D, Varambally S, Chinnaiyan AM, Department of Pathology, Michigan Center for Translational Pathology, University of Michigan Medical School, Ann Arbor,

The Polycomb group (PcG) protein EZH2 possesses oncogenic properties for which the underlying mechanism is unclear. We integrated in vitro cell line, in vivo tumor profiling, and genome-wide location data to nominate key targets of EZH2. One of the candidates identified was ADRB2 (Adrenergic Receptor, β-2), a critical mediator of beta-adrenergic signaling. EZH2 is recruited to the ADRB2 promoter and represses ADRB2 expression. ADRB2 inhibition confers cell invasion and transforms benign prostate epithelial cells, whereas ADRB2 overexpression counteracts EZH2-mediated oncogenesis. ADRB2 is underexpressed in metastatic prostate cancer, and clinically localized tumors that express lower levels of ADRB2 exhibit a poor prognosis. Taken together, we demonstrate the power of integrating multiple diverse genomic data to decipher targets of disease-related genes.

Computational approaches to study transcriptional regulation

In recent years, a number of technical and experimental advances have allowed us to obtain an unprecedented amount of information about living systems on a genomic scale. Although the complete genomes of many organisms are available due to the progress made in sequencing technology, the challenge to understand how the individual genes are regulated within the cell remains. Here, I provide an overview of current computational methods to investigate transcriptional regulation. I will first discuss how representing protein-DNA interactions as a network provides us with a conceptual framework to understand the organization of regulatory interactions in an organism. I will then describe methods to predict transcription factors and cis-regulatory elements using information such as sequence, structure and evolutionary conservation. Finally, I will discuss approaches to infer genome-scale transcriptional regulatory networks using experimentally characterized interactions from model organisms and by reverse-engineering regulatory interactions that makes use of gene expression data and genomewide location data. The methods summarized here can be exploited to discover previously uncharacterized transcriptional pathways in organisms whose genome sequence is known. In addition, such a framework and approach can be invaluable to investigate transcriptional regulation in complex microbial communities such as the human gut flora or populations of emerging pathogens. Apart from these medical applications, the concepts and methods discussed can be used to understand the combinatorial logic of transcriptional regulation and can be exploited in biotechnological applications, such as in synthetic biology experiments aimed at engineering regulatory circuits for various purposes.

Factoring local sequence composition in motif significance analysis.

We recently introduced a biologically realistic and reliable significance analysis of the output of a popular class of motif finders. In this paper we further improve our significance analysis by incorporating local base composition information. Relying on realistic biological data simulation, as well as on FDR analysis applied to real data, we show that our method is significantly better than the increasingly popular practice of using the normal approximation to estimate the significance of a finder's output. Finally we turn to leveraging our reliable significance analysis to improve the actual motif finding task. Specifically, endowing a variant of the Gibbs Sampler with our improved significance analysis we demonstrate that de novo finders can perform better than has been perceived. Significantly, our new variant outperforms all the finders reviewed in a recently published comprehensive analysis of the Harbison genome-wide binding location data. Interestingly, many of these finders incorporate additional information such as nucleosome positioning and the significance of binding data.

Statistical methods to infer cooperative binding among transcription factors in Saccharomyces cerevisiae

Transcription factors regulate transcription in prokaryotes and eukaryotes by binding to specific DNA sequences in the regulatory regions of the genes. This regulation usually occurs in a coordinated manner involving multiple transcription factors. Genome-wide location data, also called ChIP-chip data, have enabled researchers to infer the binding sites for individual regulatory proteins. However, current methods to infer binding sites, such as simple thresholding based on p-values, are not optimal for a number of study objectives like combinatorial regulation, leading to potential loss of information. Hence, there is a need to develop more efficient statistical methods for analyzing such data. We propose to use log-linear models to study cooperative binding among transcription factors and have developed an Expectation-Maximization algorithm for statistical inferences. Our method is advantageous over simple thresholding methods both based on simulation and real data studies. We apply our method to infer the cooperative network of 204 regulators in Rich Medium and a subset of them in four different environmental conditions. Our results indicate that the cooperative network is condition specific; for a set of regulators, the network structure changes under different environmental conditions. Our program is available at http://bioinformatics.med.yale.edu/TFcooperativity.

Integrative Genomics Analysis Reveals Silencing of β-Adrenergic Signaling by Polycomb in Prostate Cancer

The Polycomb group (PcG) protein EZH2 possesses oncogenic properties for which the underlying mechanism is unclear. We integrated in vitro cell line, in vivo tumor profiling, and genome-wide location data to nominate key targets of EZH2. One of the candidates identified was ADRB2 (Adrenergic Receptor, Beta-2), a critical mediator of beta-adrenergic signaling. EZH2 is recruited to the ADRB2 promoter and represses ADRB2 expression. ADRB2 inhibition confers cell invasion and transforms benign prostate epithelial cells, whereas ADRB2 overexpression counteracts EZH2-mediated oncogenesis. ADRB2 is underexpressed in metastatic prostate cancer, and clinically localized tumors that express lower levels of ADRB2 exhibit a poor prognosis. Taken together, we demonstrate the power of integrating multiple diverse genomic data to decipher targets of disease-related genes.

Incorporating prior information via shrinkage: a combined analysis of genome‐wide location data and gene expression data

Transcriptional control is a critical step in regulation of gene expression. Understanding such a control on a genomic level involves deciphering the mechanisms and structures of regulatory programmes and networks. A difficulty arises due to the weak signal and high noise in various sources of data while most current approaches are limited to analysis of a single source of data. A natural alternative is to improve statistical efficiency and power by a combined analysis of multiple sources of data. Here we propose a shrinkage method to combine genome-wide location data and gene expression data to detect the binding sites or target genes of a transcription factor. Specifically, a prior 'non-target' gene list is generated by analysing the expression data, and then this information is incorporated into the subsequent binding data analysis via a shrinkage method. There is a Bayesian justification for this shrinkage method. Both simulated and real data were used to evaluate the proposed method and compare it with analysing binding data alone. In simulation studies, the proposed method gives higher sensitivity and lower false discovery rate (FDR) in detecting the target genes. In real data example, the proposed method can reduce the estimated FDR and increase the power to detect the previously known target genes of a broad transcription regulator, leucine responsive regulatory protein (Lrp) in Escherichia coli. This method can also be used to incorporate other information, such as gene ontology (GO), to microarray data analysis to detect differentially expressed genes.

Identification of transcription factor cooperativity via stochastic system model

Transcription factor binding sites are known to co-occur in the same gene owing to cooperativity of the transcription factors (TFs) that bind to them. Genome-wide location data can help us understand how an individual TF regulates its target gene. Nevertheless, how TFs cooperate to regulate their target genes still needs further study. In this study, genome-wide location data and expression profiles are integrated to reveal how TFs cooperate to regulate their target genes from the stochastic system perspective. Based on a stochastic dynamic model, a new measurement of TF cooperativity is developed according to the regulatory abilities of cooperative TF pairs and the number of their occurrences. Our method is employed to the yeast cell cycle and reveals successfully many cooperative TF pairs confirmed by previous experiments, e.g. Swi4-Swi6 in G1/S phase and Ndd1-Fkh2 in G2/M phase. Other TF pairs with potential cooperativity mentioned in our results can provide new directions for future experiments. Finally, a cooperative TF network of cell cycle is constructed from significant cooperative TF pairs. bschen@ee.nthu.edu.tw http://www.ee.nthu.edu.tw/~bschen/cooperativity/

Incorporating Biological Information as a Prior in an Empirical Bayes Approach to Analyzing Microarray Data

Currently the practice of using existing biological knowledge in analyzing high throughput genomic and proteomic data is mainly for the purpose of validations. Here we take a different approach of incorporating biological knowledge into statistical analysis to improve statistical power and efficiency. Specifically, we consider how to fuse biological information into a mixture model to analyze microarray data. In contrast to a standard mixture model where it is assumed that all the genes come from the same (marginal) distribution, including an equal prior probability of having an event, such as having differential expression or being bound by a transcription factor (TF), our proposed mixture model allows the genes in different groups to have different distributions while the grouping of the genes reflects biological information. Using a list of about 800 putative cell cycle-regulated genes as prior biological knowledge, we analyze a genome-wide location data to detect binding sites of TF Fkh1. We find that our proposal improves over the standard approach, resulting in reduced false discovery rates (FDR), and hence it is a useful alternative to the current practice.

Learning module networks from genome-wide location and expression data

We develop a systematic algorithm for discovering network of regulatory modules, which identifies regulatory modules and their regulation program by integrating genome-wide location and expression data. Unlike previous approaches [Eisen, M.B., Spellman, P.T., Brown, P.O. and Botstein, D. (1998) Proc. Natl. Acad. Sci. USA 95, 14863–14868; Tavazoie, S., Hughes, J.D., Campbell, M.J., Cho, R.J. and Church, G.M. (1999) Nat. Genet. 22, 281–285; Ihmels, J., Friedlander, G., Bergmann, S., Sarig, O., Ziv, Y. and Barkai, N. (2002) Nat. Genet. 31, 370–377; Segal, E., Shapira, M., Regev, A., Pe’er, D., Botstein, D., Koller, D. and Friedman, N. (2003) Nat. Genet. 34, 166–176] that relied primarily on gene expression data, our algorithm regards the regulator binding data as prior knowledge that provide direct evidence of physical regulatory interactions. We applied the method to a Saccharomyces cerevisiae genome-wide location data [Lee, T.I., Rinaldi, N.J., Robert, F., Odom, D.T., Bar-Joseph, Z., Gerber, G.K., Hannett, N.M., Harbison, C.T., Thompson, C.M., Simon, I., Zeitlinger, J., Jennings, E.G., Murray, H.L. Gordon, D.B., Ren, B., Wyrick, J.J., Tagne, J.B., Volkert, T.L., Fraenkel, E., Gifford, D.K. and Young, R.A. (2002) Science 298, 799–804] for 106 DNA-binding transcription factors and 250 gene expression experiments under the conditions from the cell cycle to responses to various stress conditions. The results show that our method is able to identify functionally coherent modules and their proper regulators. Supplementary materials are available at http://compbio.sibsnet.org/projects/module-network/.

A Mixed Model Approach to Identify Yeast Transcriptional Regulatory Motifs via Microarray Experiments

A genome-wide location analysis method has been introduced as a means to simultaneously study protein-DNA binding interactions for a large number of genes on a microarray platform. Identification of interactions between transcription factors (TF) and genes provide insight into the mechanisms that regulate a variety of cellular responses. Drawing proper inferences from the experimental data is key to finding statistically significant TF-gene binding interactions. We describe how the analysis and interpretation of genome-wide location data can be fit into a traditional statistical modeling framework that considers the data across all arrays and formulizes appropriate hypothesis tests. The approach is illustrated with data from a yeast transcription factor binding experiment that illustrates how identified TF-gene interactions can enhance initial exploration of transcriptional regulatory networks. Examples of five kinds of transcriptional regulatory structure are also demonstrated. Some stark differences with previously published results are explored.