Chongzhi Zang

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.

Genome-wide Mapping of HATs and HDACs Reveals Distinct Functions in Active and Inactive Genes

Histone acetyltransferases (HATs) and deacetylases (HDACs) function antagonistically to control histone acetylation. As acetylation is a histone mark for active transcription, HATs have been associated with active and HDACs with inactive genes. We describe here genome-wide mapping of HATs and HDACs binding on chromatin and find that both are found at active genes with acetylated histones. Our data provide evidence that HATs and HDACs are both targeted to transcribed regions of active genes by phosphorylated RNA Pol II. Furthermore, the majority of HDACs in the human genome function to reset chromatin by removing acetylation at active genes. Inactive genes that are primed by MLL-mediated histone H3K4 methylation are subject to a dynamic cycle of acetylation and deacetylation by transient HAT/HDAC binding, preventing Pol II from binding to these genes but poising them for future activation. Silent genes without any H3K4 methylation signal show no evidence of being bound by HDACs.

Global Mapping of H3K4me3 and H3K27me3 Reveals Specificity and Plasticity in Lineage Fate Determination of Differentiating CD4+ T Cells

Multipotential naive CD4(+) T cells differentiate into distinct lineages including T helper 1 (Th1), Th2, Th17, and inducible T regulatory (iTreg) cells. The remarkable diversity of CD4(+) T cells begs the question whether the observed changes reflect terminal differentiation with heritable epigenetic modifications or plasticity in T cell responses. We generated genome-wide histone H3 lysine 4 (H3K4) and lysine 27 (H3K27) trimethylation maps in naive, Th1, Th2, Th17, iTreg, and natural Treg (nTreg) cells. We found that although modifications of signature-cytokine genes (Ifng, Il4, and Il17) partially conform to the expectation of lineage commitment, genes encoding transcription factors like Tbx21 exhibit a broad spectrum of epigenetic states, consistent with our demonstration of T-bet and interferon-gamma induction in nTreg cells. Our data suggest an epigenetic mechanism underlying the specificity and plasticity of effector and regulatory T cells and also provide a framework for understanding complexity of CD4(+) T helper cell differentiation.

A clustering approach for identification of enriched domains from histone modification ChIP-Seq data

MOTIVATION Chromatin states are the key to gene regulation and cell identity. Chromatin immunoprecipitation (ChIP) coupled with high-throughput sequencing (ChIP-Seq) is increasingly being used to map epigenetic states across genomes of diverse species. Chromatin modification profiles are frequently noisy and diffuse, spanning regions ranging from several nucleosomes to large domains of multiple genes. Much of the early work on the identification of ChIP-enriched regions for ChIP-Seq data has focused on identifying localized regions, such as transcription factor binding sites. Bioinformatic tools to identify diffuse domains of ChIP-enriched regions have been lacking. RESULTS Based on the biological observation that histone modifications tend to cluster to form domains, we present a method that identifies spatial clusters of signals unlikely to appear by chance. This method pools together enrichment information from neighboring nucleosomes to increase sensitivity and specificity. By using genomic-scale analysis, as well as the examination of loci with validated epigenetic states, we demonstrate that this method outperforms existing methods in the identification of ChIP-enriched signals for histone modification profiles. We demonstrate the application of this unbiased method in important issues in ChIP-Seq data analysis, such as data normalization for quantitative comparison of levels of epigenetic modifications across cell types and growth conditions. AVAILABILITY http://home.gwu.edu/ approximately wpeng/Software.htm. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Chromatin Signatures in Multipotent Human Hematopoietic Stem Cells Indicate the Fate of Bivalent Genes during Differentiation

Histone modifications have been implicated in stem cell maintenance and differentiation. We have analyzed genome-wide changes in gene expression and histone modifications during differentiation of multipotent human primary hematopoietic stem cells/progenitor cells (HSCs/HPCs) into erythrocyte precursors. Our data indicate that H3K4me1, H3K9me1, and H3K27me1 associate with enhancers of differentiation genes prior to their activation and correlate with basal expression, suggesting that these monomethylations are involved in the maintenance of activation potential required for differentiation. In addition, although the majority of genes associated with both H3K4me3 and H3K27me3 in HSCs/HPCs become silent and lose H3K4me3 after differentiation, those that lose H3K27me3 and become activated after differentiation are associated with increased levels of H2A.Z, H3K4me1, H3K9me1, H4K20me1, and RNA polymerase II in HSCs/HPCs. Thus, our data suggest that gene expression changes during differentiation are programmed by chromatin modifications present at the HSC/HPC stage and provide a resource for enhancer and promoter identification.

Partitioning heritability by functional annotation using genome-wide association summary statistics

Recent work has demonstrated that some functional categories of the genome contribute disproportionately to the heritability of complex diseases. Here we analyze a broad set of functional elements, including cell type–specific elements, to estimate their polygenic contributions to heritability in genome-wide association studies (GWAS) of 17 complex diseases and traits with an average sample size of 73,599. To enable this analysis, we introduce a new method, stratified LD score regression, for partitioning heritability from GWAS summary statistics while accounting for linked markers. This new method is computationally tractable at very large sample sizes and leverages genome-wide information. Our findings include a large enrichment of heritability in conserved regions across many traits, a very large immunological disease–specific enrichment of heritability in FANTOM5 enhancers and many cell type–specific enrichments, including significant enrichment of central nervous system cell types in the heritability of body mass index, age at menarche, educational attainment and smoking behavior.

Discrete Roles of STAT4 and STAT6 Transcription Factors in Tuning Epigenetic Modifications and Transcription during T Helper Cell Differentiation

Signal transducer and activator of transcription 4 (STAT4) and STAT6 are key factors in the specification of helper T cells; however, their direct roles in driving differentiation are not well understood. Using chromatin immunoprecipitation and massive parallel sequencing, we quantitated the full complement of STAT-bound genes, concurrently assessing global STAT-dependent epigenetic modifications and gene transcription by using cells from cognate STAT-deficient mice. STAT4 and STAT6 each bound over 4000 genes with distinct binding motifs. Both played critical roles in maintaining chromatin configuration and transcription of a core subset of genes through the combination of different epigenetic patterns. Globally, STAT4 had a more dominant role in promoting active epigenetic marks, whereas STAT6 had a more prominent role in antagonizing repressive marks. Clusters of genes negatively regulated by STATs were also identified, highlighting previously unappreciated repressive roles of STATs. Therefore, STAT4 and STAT6 play wide regulatory roles in T helper cell specification.

Genome-wide analysis of histone methylation reveals chromatin state-based regulation of gene transcription and function of memory CD8+ T cells.

Memory lymphocytes are characterized by their ability to exhibit a rapid response to the recall antigen, in which differential transcription is important, yet the underlying mechanism is not understood. We report here a genome-wide analysis of histone methylation on two histone H3 lysine residues (H3K4me3 and H3K27me3) and gene expression profiles in naive and memory CD8(+) T cells. We found that specific correlation exists between gene expression and the amounts of H3K4me3 (positive correlation) and H3K27me3 (negative correlation) across the gene body. These correlations displayed four distinct modes (repressive, active, poised, and bivalent), reflecting different functions of these genes in CD8(+) T cells. Furthermore, a permissive chromatin state of each gene was established by a combination of different histone modifications. Our findings reveal a complex regulation by histone methylation in differential gene expression and suggest that histone methylation may be responsible for memory CD8(+) T cell function.

PTIP Promotes Chromatin Changes Critical for Immunoglobulin Class Switch Recombination

PTIP in Immunoglobulin Switching One of the hallmarks of humoral immunity is the ability of immunoglobulins (Ig) to undergo class switch recombination (CSR). Through genetic recombining of the Ig heavy chain, Igs maintain their antigen specificity but gain the ability to interact with different cell surface receptors required for successful pathogen clearance. CSR requires transcription at the Ig heavy chain locus to initiate genetic rearrangement. Changes in chromatin accessibility are thought to promote CSR-associated transcription. Daniel et al. (p. 917, published online 29 July; see the Perspective by Singh and Demarco) now show that trimethylation of histone 3 at lysine 4 (H3K4me3) controls the accessibility of the Ig heavy chain locus to CSR and that PTIP (Pax interaction with transcription-activation domain protein-1), a component of the histone methylase complex, is required for this modification. Mouse PTIP-deficient B cells exhibited impaired CSR. PTIP was required both for the recruitment RNA Polymerase II and for subsequent chromatin remodeling, including histone acetylation, which occurs during CSR. Largely independent of its function in transcription initiation in CSR, PTIP also associated with double-stranded DNA breaks during CSR and promoted genome stability. These dual functions of PTIP may be important for the precise coordination of chromatin accessibility and recombination required during CSR. A factor that regulates chromatin accessibility and recombination during immunoglobulin rearrangements is identified. Programmed genetic rearrangements in lymphocytes require transcription at antigen receptor genes to promote accessibility for initiating double-strand break (DSB) formation critical for DNA recombination and repair. Here, we showed that activated B cells deficient in the PTIP component of the MLL3 (mixed-lineage leukemia 3)–MLL4 complex display impaired trimethylation of histone 3 at lysine 4 (H3K4me3) and transcription initiation of downstream switch regions at the immunoglobulin heavy-chain (Igh) locus, leading to defective immunoglobulin class switching. We also showed that PTIP accumulation at DSBs contributes to class switch recombination (CSR) and genome stability independently of Igh switch transcription. These results demonstrate that PTIP promotes specific chromatin changes that control the accessibility of the Igh locus to CSR and suggest a nonredundant role for the MLL3-MLL4 complex in altering antibody effector function.

Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification.

Sequencing of DNase I hypersensitive sites (DNase-seq) is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize performance of DNase-seq. Sequencing short fragments of 50–100 base pairs (bp) that accumulate in long internucleosome linker regions was more efficient for identifying transcription factor binding sites compared to sequencing longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy via generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNase I cutting bias, we found a strong effect that varied over more than two orders of magnitude. This indicates that the nucleotide-resolution cleavage patterns at many transcription factor binding sites are derived from intrinsic DNase I cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNase I hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.