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

ChIP-seq analysis reveals distinct H3K27me3 profiles that correlate with transcriptional activity

Transcriptional control is dependent on a vast network of epigenetic modifications. One epigenetic mark of particular interest is tri-methylation of lysine 27 on histone H3 (H3K27me3), which is catalysed and maintained by Polycomb Repressive Complex 2 (PRC2). Although this histone mark is studied widely, the precise relationship between its local pattern of enrichment and regulation of gene expression is currently unclear. We have used ChIP-seq to generate genome-wide maps of H3K27me3 enrichment, and have identified three enrichment profiles with distinct regulatory consequences. First, a broad domain of H3K27me3 enrichment across the body of genes corresponds to the canonical view of H3K27me3 as inhibitory to transcription. Second, a peak of enrichment around the transcription start site (TSS) is commonly associated with ‘bivalent’ genes, where H3K4me3 also marks the TSS. Finally and most surprisingly, we identified an enrichment profile with a peak in the promoter of genes that is associated with active transcription. Genes with each of these three profiles were found in different proportions in each of the cell types studied. The data analysis techniques developed here will be useful for the identification of common enrichment profiles for other histone modifications that have important consequences for transcriptional regulation

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>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.</p

Identification of an enhancer that increases miR-200b~200a~429 gene expression in breast cancer cells

The miR-200b~200a~429 gene cluster is a key regulator of EMT and cancer metastasis, however the transcription-based mechanisms controlling its expression during this process are not well understood. We have analyzed the miR-200b~200a~429 locus for epigenetic modifications in breast epithelial and mesenchymal cell lines using chromatin immunoprecipitation assays and DNA methylation analysis. We discovered a novel enhancer located approximately 5.1kb upstream of the miR-200b~200a~429 transcriptional start site. This region was associated with the active enhancer chromatin signature comprising H3K4me1, H3K27ac, RNA polymerase II and CpG dinucleotide hypomethylation. Luciferase reporter assays revealed the upstream enhancer stimulated the transcription of the miR-200b~200a~429 minimal promoter region approximately 27-fold in breast epithelial cells. Furthermore, we found that a region of the enhancer was transcribed, producing a short, GC-rich, mainly nuclear, non-polyadenylated RNA transcript designated miR-200b eRNA. Over-expression of miR-200b eRNA had little effect on miR-200b~200a~429 promoter activity and its production did not correlate with miR-200b~200a~429 gene expression. While additional investigations of miR-200b eRNA function will be necessary, it is possible that miR-200b eRNA may be involved in the regulation of miR-200b~200a~429 gene expression and silencing. Taken together, these findings reveal the presence of a novel enhancer, which contributes to miR-200b~200a~429 transcriptional regulation in epithelial cells.Joanne L. Attema, Andrew G. Bert, Yat-Yuen Lim, Natasha Kolesnikoff, David M. Lawrence, Katherine A. Pillman, Eric Smith, Paul A. Drew, Yeesim Khew-Goodall, Frances Shannon, Gregory J. Goodal

PcG methylation of the HIST1 cluster defines an epigenetic marker of acute myeloid leukemia

International audienceno abstrac

Genetic recombination is directed away from functional genomic elements in mice.

International audienceGenetic recombination occurs during meiosis, the key developmental programme of gametogenesis. Recombination in mammals has been recently linked to the activity of a histone H3 methyltransferase, PR domain containing 9 (PRDM9), the product of the only known speciation-associated gene in mammals. PRDM9 is thought to determine the preferred recombination sites--recombination hotspots--through sequence-specific binding of its highly polymorphic multi-Zn-finger domain. Nevertheless, Prdm9 knockout mice are proficient at initiating recombination. Here we map and analyse the genome-wide distribution of recombination initiation sites in Prdm9 knockout mice and in two mouse strains with different Prdm9 alleles and their F(1) hybrid. We show that PRDM9 determines the positions of practically all hotspots in the mouse genome, with the exception of the pseudo-autosomal region (PAR)--the only area of the genome that undergoes recombination in 100% of cells. Surprisingly, hotspots are still observed in Prdm9 knockout mice, and as in wild type, these hotspots are found at H3 lysine 4 (H3K4) trimethylation marks. However, in the absence of PRDM9, most recombination is initiated at promoters and at other sites of PRDM9-independent H3K4 trimethylation. Such sites are rarely targeted in wild-type mice, indicating an unexpected role of the PRDM9 protein in sequestering the recombination machinery away from gene-promoter regions and other functional genomic elements