Iouri Chepelev

High-resolution profiling of histone methylations in the human genome

Histone modifications are implicated in influencing gene expression. We have generated high-resolution maps for the genome-wide distribution of 20 histone lysine and arginine methylations as well as histone variant H2A.Z, RNA polymerase II, and the insulator binding protein CTCF across the human genome using the Solexa 1G sequencing technology. Typical patterns of histone methylations exhibited at promoters, insulators, enhancers, and transcribed regions are identified. The monomethylations of H3K27, H3K9, H4K20, H3K79, and H2BK5 are all linked to gene activation, whereas trimethylations of H3K27, H3K9, and H3K79 are linked to repression. H2A.Z associates with functional regulatory elements, and CTCF marks boundaries of histone methylation domains. Chromosome banding patterns are correlated with unique patterns of histone modifications. Chromosome breakpoints detected in T cell cancers frequently reside in chromatin regions associated with H3K4 methylations. Our data provide new insights into the function of histone methylation and chromatin organization in genome function.

Intragenic DNA methylation modulates alternative splicing by recruiting MeCP2 to promote exon recognition.

Although the function of DNA methylation in gene promoter regions is well established in transcriptional repression, the function of the evolutionarily conserved widespread distribution of DNA methylation in gene body regions remains incompletely understood. Here, we show that DNA methylation is enriched in included alternatively spliced exons (ASEs), and that inhibition of DNA methylation results in aberrant splicing of ASEs. The methyl-CpG-binding protein MeCP2 is enriched in included ASEs, particularly those that are also highly methylated, and inhibition of DNA methylation disrupts specific targeting of MeCP2 to exons. Interestingly, ablation of MeCP2 results in increased histone acetylation and aberrant ASE-skipping events. We further show that inhibition of histone deacetylase (HDAC) activity leads to exon skipping that shows a highly significant degree of overlap with that caused by MeCP2 knockdown. Together, our data indicate that intragenic DNA methylation operates in exon definition to modulate alternative RNA splicing and can enhance exon recognition via recruitment of the multifunctional protein MeCP2, which thereby maintains local histone hypoacetylation through the subsequent recruitment of HDACs.

Detection of single nucleotide variations in expressed exons of the human genome using RNA-Seq

Whole-genome resequencing is still a costly method to detect genetic mutations that lead to altered forms of proteins and may be associated with disease development. Since the majority of disease-related single nucleotide variations (SNVs) are found in protein-coding regions, we propose to identify SNVs in expressed exons of the human genome using the recently developed RNA-Seq technique. We identify 12 176 and 10 621 SNVs, respectively, in Jurkat T cells and CD4+ T cells from a healthy donor. Interestingly, our data show that one copy of the TAL-1 proto-oncogene has a point mutation in 3′ UTR and only the mutant allele is expressed in Jurkat cells. We provide a comprehensive dataset for further understanding the cancer biology of Jurkat cells. Our results indicate that this is a cost-effective and efficient strategy to systematically identify SNVs in the expressed regions of the human genome.

Characterization of genome-wide enhancer-promoter interactions reveals co-expression of interacting genes and modes of higher order chromatin organization

Recent epigenomic studies have predicted thousands of potential enhancers in the human genome. However, there has not been systematic characterization of target promoters for these potential enhancers. Using H3K4me2 as a mark for active enhancers, we identified genome-wide EP interactions in human CD4+ T cells. Among the 6 520 long-distance chromatin interactions, we identify 2 067 enhancers that interact with 1 619 promoters and enhance their expression. These enhancers exist in accessible chromatin regions and are associated with various histone modifications and polymerase II binding. The promoters with interacting enhancers are expressed at higher levels than those without interacting enhancers, and their expression levels are positively correlated with the number of interacting enhancers. Interestingly, interacting promoters are co-expressed in a tissue-specific manner. We also find that chromosomes are organized into multiple levels of interacting domains. Our results define a global view of EP interactions and provide a data set to further understand mechanisms of enhancer targeting and long-range chromatin organization. The Gene Expression Omnibus accession number for the raw and analyzed chromatin interaction data is GSE32677.

Pol II and its associated epigenetic marks are present at Pol III-transcribed noncoding RNA genes

Epigenetic control is an important aspect of gene regulation. Despite detailed understanding of protein-coding gene expression, the transcription of noncoding RNA genes by RNA polymerase III (Pol III) is less well characterized. Here we profile the epigenetic features of Pol III target genes throughout the human genome. This reveals that the chromatin landscape of Pol III–transcribed genes resembles that of Pol II templates in many ways, although there are also clear differences. Our analysis also uncovered an entirely unexpected phenomenon: namely, that Pol II is present at the majority of genomic loci that are bound by Pol III.

Dynamic regulation of alternative splicing and chromatin structure in Drosophila gonads revealed by RNA-seq

Both transcription and post-transcriptional processes, such as alternative splicing, play crucial roles in controlling developmental programs in metazoans. Recently emerged RNA-seq method has brought our understanding of eukaryotic transcriptomes to a new level, because it can resolve both gene expression level and alternative splicing events simultaneously. To gain a better understanding of cellular differentiation in gonads, we analyzed mRNA profiles from Drosophila testes and ovaries using RNA-seq. We identified a set of genes that have sex-specific isoforms in wild-type (WT) gonads, including several transcription factors. We found that differentiation of sperms from undifferentiated germ cells induced a dramatic downregulation of RNA splicing factors. Our data confirmed that RNA splicing events are significantly more frequent in the undifferentiated cell-enriched bag of marbles (bam) mutant testis, but downregulated upon differentiation in WT testis. Consistent with this, we showed that genes required for meiosis and terminal differentiation in WT testis were mainly regulated at the transcriptional level, but not by alternative splicing. Unexpectedly, we observed an increase in expression of all families of chromatin remodeling factors and histone modifying enzymes in the undifferentiated cell-enriched bam testis. More interestingly, chromatin regulators and histone modifying enzymes with opposite enzymatic activities are coenriched in undifferentiated cells in testis, suggesting that these cells may possess dynamic chromatin architecture. Finally, our data revealed many new features of the Drosophila gonadal transcriptomes, and will lead to a more comprehensive understanding of how differential gene expression and splicing regulate gametogenesis in Drosophila. Our data provided a foundation for the systematic study of gene expression and alternative splicing in many interesting areas of germ cell biology in Drosophila, such as the molecular basis for sexual dimorphism and the regulation of the proliferation vs terminal differentiation programs in germline stem cell lineages. The GEO accession number for the raw and analyzed RNA-seq data is GSE16960.

Prediction of regulatory elements in mammalian genomes using chromatin signatures.

BackgroundRecent genomic scale survey of epigenetic states in the mammalian genomes has shown that promoters and enhancers are correlated with distinct chromatin signatures, providing a pragmatic way for systematic mapping of these regulatory elements in the genome. With rapid accumulation of chromatin modification profiles in the genome of various organisms and cell types, this chromatin based approach promises to uncover many new regulatory elements, but computational methods to effectively extract information from these datasets are still limited.ResultsWe present here a supervised learning method to predict promoters and enhancers based on their unique chromatin modification signatures. We trained Hidden Markov models (HMMs) on the histone modification data for known promoters and enhancers, and then used the trained HMMs to identify promoter or enhancer like sequences in the human genome. Using a simulated annealing (SA) procedure, we searched for the most informative combination and the optimal window size of histone marks.ConclusionCompared with the previous methods, the HMM method can capture the complex patterns of histone modifications particularly from the weak signals. Cross validation and scanning the ENCODE regions showed that our method outperforms the previous profile-based method in mapping promoters and enhancers. We also showed that including more histone marks can further boost the performance of our method. This observation suggests that the HMM is robust and is capable of integrating information from multiple histone marks. To further demonstrate the usefulness of our method, we applied it to analyzing genome wide ChIP-Seq data in three mouse cell lines and correctly predicted active and inactive promoters with positive predictive values of more than 80%. The software is available at http://http:/nash.ucsd.edu/chromatin.tar.gz.

Proteogenomic characterization and mapping of nucleosomes decoded by Brd and HP1 proteins

BackgroundHistone post-translational modifications (PTMs) constitute a branch of epigenetic mechanisms that can control the expression of eukaryotic genes in a heritable manner. Recent studies have identified several PTM-binding proteins containing diverse specialized domains whose recognition of specific PTM sites leads to gene activation or repression. Here, we present a high-throughput proteogenomic platform designed to characterize the nucleosomal make-up of chromatin enriched with a set of histone PTM binding proteins known as histone PTM readers. We support our findings with gene expression data correlating to PTM distribution.ResultsWe isolated human mononucleosomes bound by the bromodomain-containing proteins Brd2, Brd3 and Brd4, and by the chromodomain-containing heterochromatin proteins HP1β and HP1α. Histone PTMs were quantified by mass spectrometry (ChIP-qMS), and their associated DNAs were mapped using deep sequencing. Our results reveal that Brd- and HP1-bound nucleosomes are enriched in histone PTMs consistent with actively transcribed euchromatin and silent heterochromatin, respectively. Data collected using RNA-Seq show that Brd-bound sites correlate with highly expressed genes. In particular, Brd3 and Brd4 are most enriched on nucleosomes located within HOX gene clusters, whose expression is reduced upon Brd4 depletion by short hairpin RNA.ConclusionsProteogenomic mapping of histone PTM readers, alongside the characterization of their local chromatin environments and transcriptional information, should prove useful for determining how histone PTMs are bound by these readers and how they contribute to distinct transcriptional states.

Cell Fate Determination Factor Dachshund Reprograms Breast Cancer Stem Cell Function

The cell fate determination factor Dachshund was cloned as a dominant inhibitor of the hyperactive epidermal growth factor receptor ellipse. The expression of Dachshund is lost in human breast cancer associated with poor prognosis. Breast tumor-initiating cells (TIC) may contribute to tumor progression and therapy resistance. Here, endogenous DACH1 was reduced in breast cancer cell lines with high expression of TIC markers and in patient samples of the basal breast cancer phenotype. Re-expression of DACH1 reduced new tumor formation in serial transplantations in vivo, reduced mammosphere formation, and reduced the proportion of CD44high/CD24low breast tumor cells. Conversely, lentiviral shRNA to DACH1 increased the breast (B)TIC population. Genome-wide expression studies of mammary tumors demonstrated DACH1 repressed a molecular signature associated with stem cells (SOX2, Nanog, and KLF4) and genome-wide ChIP-seq analysis identified DACH1 binding to the promoter of the Nanog, KLF4, and Lin28 genes. KLF4/c-Myc and Oct4/Sox2 antagonized DACH1 repression of BTIC. Mechanistic studies demonstrated DACH1 directly repressed the Nanog and Sox2 promoters via a conserved domain. Endogenous DACH1 regulates BTIC in vitro and in vivo.

Attenuation of Forkhead signaling by the retinal determination factor DACH1

The Drosophila Dachshund (Dac) gene, cloned as a dominant inhibitor of the hyperactive growth factor mutant ellipse, encodes a key component of the retinal determination gene network that governs cell fate. Herein, cyclic amplification and selection of targets identified a DACH1 DNA-binding sequence that resembles the FOX (Forkhead box–containing protein) binding site. Genome-wide in silico promoter analysis of DACH1 binding sites identified gene clusters populating cellular pathways associated with the cell cycle and growth factor signaling. ChIP coupled with high-throughput sequencing mapped DACH1 binding sites to corresponding gene clusters predicted in silico and identified as weight matrix resembling the cyclic amplification and selection of targets–defined sequence. DACH1 antagonized FOXM1 target gene expression, promoter occupancy in the context of local chromatin, and contact-independent growth. Attenuation of FOX function by the cell fate determination pathway has broad implications given the diverse role of FOX proteins in cellular biology and tumorigenesis.