Suresh Cuddapah

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.

Combinatorial patterns of histone acetylations and methylations in the human genome

Histones are characterized by numerous posttranslational modifications that influence gene transcription. However, because of the lack of global distribution data in higher eukaryotic systems, the extent to which gene-specific combinatorial patterns of histone modifications exist remains to be determined. Here, we report the patterns derived from the analysis of 39 histone modifications in human CD4+ T cells. Our data indicate that a large number of patterns are associated with promoters and enhancers. In particular, we identify a common modification module consisting of 17 modifications detected at 3,286 promoters. These modifications tend to colocalize in the genome and correlate with each other at an individual nucleosome level. Genes associated with this module tend to have higher expression, and addition of more modifications to this module is associated with further increased expression. Our data suggest that these histone modifications may act cooperatively to prepare chromatin for transcriptional activation.

Dynamic Regulation of Nucleosome Positioning in the Human Genome

The positioning of nucleosomes with respect to DNA plays an important role in regulating transcription. However, nucleosome mapping has been performed for only limited genomic regions in humans. We have generated genome-wide maps of nucleosome positions in both resting and activated human CD4+ T cells by direct sequencing of nucleosome ends using the Solexa high-throughput sequencing technique. We find that nucleosome phasing relative to the transcription start sites is directly correlated to RNA polymerase II (Pol II) binding. Furthermore, the first nucleosome downstream of a start site exhibits differential positioning in active and silent genes. TCR signaling induces extensive nucleosome reorganization in promoters and enhancers to allow transcriptional activation or repression. Our results suggest that H2A.Z-containing and modified nucleosomes are preferentially lost from the -1 nucleosome position. Our data provide a comprehensive view of the nucleosome landscape and its dynamic regulation in the human genome.

Global analysis of the insulator binding protein CTCF in chromatin barrier regions reveals demarcation of active and repressive domains

Insulators are DNA elements that prevent inappropriate interactions between the neighboring regions of the genome. They can be functionally classified as either enhancer blockers or domain barriers. CTCF (CCCTC-binding factor) is the only known major insulator-binding protein in the vertebrates and has been shown to bind many enhancer-blocking elements. However, it is not clear whether it plays a role in chromatin domain barriers between active and repressive domains. Here, we used ChIP-seq to map the genome-wide binding sites of CTCF in three cell types and identified significant binding of CTCF to the boundaries of repressive chromatin domains marked by H3K27me3. Although we find an extensive overlapping of CTCF-binding sites across the three cell types, its association with the domain boundaries is cell-type-specific. We further show that the nucleosomes flanking CTCF-binding sites are well positioned. Interestingly, we found a complementary pattern between the repressive H3K27me3 and the active H2AK5ac regions, which are separated by CTCF. Our data indicate that CTCF may play important roles in the barrier activity of insulators, and this study provides a resource for further investigation of the CTCF function in organizing chromatin in the human genome.

Genome-wide identification of in vivo protein–DNA binding sites from ChIP-Seq data

ChIP-Seq, which combines chromatin immunoprecipitation (ChIP) with ultra high-throughput massively parallel sequencing, is increasingly being used for mapping protein–DNA interactions in-vivo on a genome scale. Typically, short sequence reads from ChIP-Seq are mapped to a reference genome for further analysis. Although genomic regions enriched with mapped reads could be inferred as approximate binding regions, short read lengths (∼25–50 nt) pose challenges for determining the exact binding sites within these regions. Here, we present SISSRs (Site Identification from Short Sequence Reads), a novel algorithm for precise identification of binding sites from short reads generated from ChIP-Seq experiments. The sensitivity and specificity of SISSRs are demonstrated by applying it on ChIP-Seq data for three widely studied and well-characterized human transcription factors: CTCF (CCCTC-binding factor), NRSF (neuron-restrictive silencer factor) and STAT1 (signal transducer and activator of transcription protein 1). We identified 26 814, 5813 and 73 956 binding sites for CTCF, NRSF and STAT1 proteins, respectively, which is 32, 299 and 78% more than that inferred previously for the respective proteins. Motif analysis revealed that an overwhelming majority of the identified binding sites contained the previously established consensus binding sequence for the respective proteins, thus attesting for SISSRs’ accuracy. SISSRs’ sensitivity and precision facilitated further analyses of ChIP-Seq data revealing interesting insights, which we believe will serve as guidance for designing ChIP-Seq experiments to map in vivo protein–DNA interactions. We also show that tag densities at the binding sites are a good indicator of protein–DNA binding affinity, which could be used to distinguish and characterize strong and weak binding sites. Using tag density as an indicator of DNA-binding affinity, we have identified core residues within the NRSF and CTCF binding sites that are critical for a stronger DNA binding.

The genomic landscape of histone modifications in human T cells

To understand the molecular basis that supports the dynamic gene expression programs unique to T cells, we investigated the genomic landscape of activating histone modifications, including histone H3 K9/K14 diacetylation (H3K9acK14ac), H3 K4 trimethylation (H3K4me3), and the repressive histone modification H3 K27 trimethylation (H3K27me3) in primary human T cells. We show that H3K9acK14ac and H3K4me3 are associated with active genes required for T cell function and development, whereas H3K27me3 is associated with silent genes that are involved in development in other cell types. Unexpectedly, we find that 3,330 gene promoters are associated with all of these histone modifications. The gene expression levels are correlated with both the absolute and relative levels of the activating H3K4me3 and the repressive H3K27me3 modifications. Our data reveal that rapidly inducible genes are associated with the H3 acetylation and H3K4me3 modifications, suggesting they assume a chromatin structure poised for activation. In addition, we identified a subpopulation of chromatin regions that are associated with high levels of H3K4me3 and H3K27me3 but low levels of H3K9acK14ac. Therefore, these regions have a distinctive chromatin modification pattern and thus may represent a distinct class of chromatin domains.

Chromatin poises miRNA- and protein-coding genes for expression

Chromatin modifications have been implicated in the regulation of gene expression. While association of certain modifications with expressed or silent genes has been established, it remains unclear how changes in chromatin environment relate to changes in gene expression. In this article, we used ChIP-seq (chromatin immunoprecipitation with massively parallel sequencing) to analyze the genome-wide changes in chromatin modifications during activation of total human CD4(+) T cells by T-cell receptor (TCR) signaling. Surprisingly, we found that the chromatin modification patterns at many induced and silenced genes are relatively stable during the short-term activation of resting T cells. Active chromatin modifications were already in place for a majority of inducible protein-coding genes, even while the genes were silent in resting cells. Similarly, genes that were silenced upon T-cell activation retained positive chromatin modifications even after being silenced. To investigate if these observations are also valid for miRNA-coding genes, we systematically identified promoters for known miRNA genes using epigenetic marks and profiled their expression patterns using deep sequencing. We found that chromatin modifications can poise miRNA-coding genes as well. Our data suggest that miRNA- and protein-coding genes share similar mechanisms of regulation by chromatin modifications, which poise inducible genes for activation in response to environmental stimuli.

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.

Acute depletion of Tet1-dependent 5-hydroxymethylcytosine levels impairs LIF/Stat3 signaling and results in loss of embryonic stem cell identity

The TET family of FE(II) and 2-oxoglutarate-dependent enzymes (Tet1/2/3) promote DNA demethylation by converting 5-methylcytosine to 5-hydroxymethylcytosine (5hmC), which they further oxidize into 5-formylcytosine and 5-carboxylcytosine. Tet1 is robustly expressed in mouse embryonic stem cells (mESCs) and has been implicated in mESC maintenance. Here we demonstrate that, unlike genetic deletion, RNAi-mediated depletion of Tet1 in mESCs led to a significant reduction in 5hmC and loss of mESC identity. The differentiation phenotype due to Tet1 depletion positively correlated with the extent of 5hmC loss. Meta-analyses of genomic data sets suggested interaction between Tet1 and leukemia inhibitory factor (LIF) signaling. LIF signaling is known to promote self-renewal and pluripotency in mESCs partly by opposing MAPK/ERK-mediated differentiation. Withdrawal of LIF leads to differentiation of mESCs. We discovered that Tet1 depletion impaired LIF-dependent Stat3-mediated gene activation by affecting Stat3s ability to bind to its target sites on chromatin. Nanog overexpression or inhibition of MAPK/ERK signaling, both known to maintain mESCs in the absence of LIF, rescued Tet1 depletion, further supporting the dependence of LIF/Stat3 signaling on Tet1. These data support the conclusion that analysis of mESCs in the hours/days immediately following efficient Tet1 depletion reveals Tet1s normal physiological role in maintaining the pluripotent state that may be subject to homeostatic compensation in genetic models.

Epigenomics of T cell activation, differentiation and memory

Activation of T cells is an essential step in the immunological response to infection. Although activation of naïve T cells results in proliferation and slow differentiation into cytokine-producing effector cells, antigen engagement with memory cells leads to cytokine production immediately. Even though the cell surface signaling events are similar in both the cases, the outcome is different, suggesting that distinct regulatory mechanisms may exist downstream of the activation signals. Recent advances in the understanding of global epigenetic patterns in T cells have resulted in the appreciation of the role of epigenetic mechanisms in processes such as activation and differentiation. In this review we discuss recent data suggesting that naïve T cell activation, differentiation, and lineage commitment result in epigenetic changes and a fine balance between different histone modifications is required. On the other hand, memory T cells are poised and do not require epigenetic changes for short-term activation.