Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease

Profiling of transcriptional and epigenetic changes during directed endothelial differentiation of human embryonic stem cells identifies FOXA2 as a marker of early mesoderm commitment

Introduction: Differentiation of vascular endothelial cells (ECs) in clinically relevant numbers for injection into ischaemic areas could offer therapeutic potential in the treatment of cardiovascular conditions, including myocardial infarction, peripheral vascular disease and stroke. While we and others have demonstrated successful generation of functional endothelial-like cells from human embryonic stem cells (hESCs), little is understood regarding the complex transcriptional and epigenetic changes that occur during differentiation, in particular during early commitment to a mesodermal lineage. Methods: We performed the first gene expression microarray study of hESCs undergoing directed differentiation to ECs using a monolayer-based, feeder-free and serum-free protocol. Microarray results were confirmed by quantitative RT-PCR and immunocytochemistry, and chromatin immunoprecipitation (ChIP)-PCR analysis was utilised to determine the bivalent status of differentially expressed genes. Results: We identified 22 transcription factors specific to early mesoderm commitment. Among these factors, FOXA2 was observed to be the most significantly differentially expressed at the hESC–EC day 2 timepoint. ChIP-PCR analysis revealed that the FOXA2 transcription start site is bivalently marked with histone modifications for both gene activation (H3K4me3) and repression (H3K27me3) in hESCs, suggesting the transcription factor may be a key regulator of hESC differentiation. Conclusion: This enhanced knowledge of the lineage commitment process will help improve the design of directed differentiation protocols, increasing the yield of endothelial-like cells for regenerative medicine therapies in cardiovascular disease

A signal processing approach for enriched region detection in RNA polymerase II ChIP-seq data

International audienc

Q&A: ChIP-seq technologies and the study of gene regulation

10.1186/1741-7007-8-56BMC Biology85

Repression of Esophageal Neoplasia and Inflammatory Signaling by Anti-miR-31 Delivery In Vivo.

BACKGROUND: Overexpression of microRNA-31 (miR-31) is implicated in the pathogenesis of esophageal squamous cell carcinoma (ESCC), a deadly disease associated with dietary zinc deficiency. Using a rat model that recapitulates features of human ESCC, the mechanism whereby Zn regulates miR-31 expression to promote ESCC is examined. METHODS: To inhibit in vivo esophageal miR-31 overexpression in Zn-deficient rats (n = 12-20 per group), locked nucleic acid-modified anti-miR-31 oligonucleotides were administered over five weeks. miR-31 expression was determined by northern blotting, quantitative polymerase chain reaction, and in situ hybridization. Physiological miR-31 targets were identified by microarray analysis and verified by luciferase reporter assay. Cellular proliferation, apoptosis, and expression of inflammation genes were determined by immunoblotting, caspase assays, and immunohistochemistry. The miR-31 promoter in Zn-deficient esophagus was identified by ChIP-seq using an antibody for histone mark H3K4me3. Data were analyzed with t test and analysis of variance. All statistical tests were two-sided. RESULTS: In vivo, anti-miR-31 reduced miR-31 overexpression (P = .002) and suppressed the esophageal preneoplasia in Zn-deficient rats. At the same time, the miR-31 target Stk40 was derepressed, thereby inhibiting the STK40-NF-κΒ-controlled inflammatory pathway, with resultant decreased cellular proliferation and activated apoptosis (caspase 3/7 activities, fold change = 10.7, P = .005). This same connection between miR-31 overexpression and STK40/NF-κΒ expression was also documented in human ESCC cell lines. In Zn-deficient esophagus, the miR-31 promoter region and NF-κΒ binding site were activated. Zn replenishment restored the regulation of this genomic region and a normal esophageal phenotype. CONCLUSIONS: The data define the in vivo signaling pathway underlying interaction of Zn deficiency and miR-31 overexpression in esophageal neoplasia and provide a mechanistic rationale for miR-31 as a therapeutic target for ESCC

The Genomic HyperBrowser: inferential genomics at the sequence level

The immense increase in the generation of genomic scale data poses an unmet analytical challenge, due to a lack of established methodology with the required flexibility and power. We propose a first principled approach to statistical analysis of sequence-level genomic information. We provide a growing collection of generic biological investigations that query pairwise relations between tracks, represented as mathematical objects, along the genome. The Genomic HyperBrowser implements the approach and is available at http://hyperbrowser.uio.no

Inferring nucleosome positions with their histone mark annotation from ChIP data

MOTIVATION: The nucleosome is the basic repeating unit of chromatin. It contains two copies each of the four core histones H2A, H2B, H3 and H4 and about 147 bp of DNA. The residues of the histone proteins are subject to numerous post-translational modifications, such as methylation or acetylation. Chromatin immunoprecipitiation followed by sequencing (ChIP-seq) is a technique that provides genome-wide occupancy data of these modified histone proteins, and it requires appropriate computational methods. RESULTS: We present NucHunter, an algorithm that uses the data from ChIP-seq experiments directed against many histone modifications to infer positioned nucleosomes. NucHunter annotates each of these nucleosomes with the intensities of the histone modifications. We demonstrate that these annotations can be used to infer nucleosomal states with distinct correlations to underlying genomic features and chromatin-related processes, such as transcriptional start sites, enhancers, elongation by RNA polymerase II and chromatin-mediated repression. Thus, NucHunter is a versatile tool that can be used to predict positioned nucleosomes from a panel of histone modification ChIP-seq experiments and infer distinct histone modification patterns associated to different chromatin states. AVAILABILITY: The software is available at http://epigen.molgen.mpg.de/nuchunter/. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online

lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs.

Although recent studies have indicated roles of long non-coding RNAs (lncRNAs) in physiological aspects of cell-type determination and tissue homeostasis, their potential involvement in regulated gene transcription programs remains rather poorly understood. The androgen receptor regulates a large repertoire of genes central to the identity and behaviour of prostate cancer cells, and functions in a ligand-independent fashion in many prostate cancers when they become hormone refractory after initial androgen deprivation therapy. Here we report that two lncRNAs highly overexpressed in aggressive prostate cancer, PRNCR1 (also known as PCAT8) and PCGEM1, bind successively to the androgen receptor and strongly enhance both ligand-dependent and ligand-independent androgen-receptor-mediated gene activation programs and proliferation in prostate cancer cells. Binding of PRNCR1 to the carboxy-terminally acetylated androgen receptor on enhancers and its association with DOT1L appear to be required for recruitment of the second lncRNA, PCGEM1, to the androgen receptor amino terminus that is methylated by DOT1L. Unexpectedly, recognition of specific protein marks by PCGEM1-recruited pygopus 2 PHD domain enhances selective looping of androgen-receptor-bound enhancers to target gene promoters in these cells. In 'resistant' prostate cancer cells, these overexpressed lncRNAs can interact with, and are required for, the robust activation of both truncated and full-length androgen receptor, causing ligand-independent activation of the androgen receptor transcriptional program and cell proliferation. Conditionally expressed short hairpin RNA targeting these lncRNAs in castration-resistant prostate cancer cell lines strongly suppressed tumour xenograft growth in vivo. Together, these results indicate that these overexpressed lncRNAs can potentially serve as a required component of castration-resistance in prostatic tumours