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.

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.

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 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.

c-Myc is a universal amplifier of expressed genes in lymphocytes and embryonic stem cells

The c-Myc HLH-bZIP protein has been implicated in physiological or pathological growth, proliferation, apoptosis, metabolism, and differentiation at the cellular, tissue, or organismal levels via regulation of numerous target genes. No principle yet unifies Myc action due partly to an incomplete inventory and functional accounting of Myc's targets. To observe Myc target expression and function in a system where Myc is temporally and physiologically regulated, the transcriptomes and the genome-wide distributions of Myc, RNA polymerase II, and chromatin modifications were compared during lymphocyte activation and in ES cells as well. A remarkably simple rule emerged from this quantitative analysis: Myc is not an on-off specifier of gene activity, but is a nonlinear amplifier of expression, acting universally at active genes, except for immediate early genes that are strongly induced before Myc. This rule of Myc action explains the vast majority of Myc biology observed in literature.

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.

Dual functions of Tet1 in transcriptional regulation in mouse embryonic stem cells

Epigenetic modification of the mammalian genome by DNA methylation (5-methylcytosine) has a profound impact on chromatin structure, gene expression and maintenance of cellular identity. The recent demonstration that members of the Ten-eleven translocation (Tet) family of proteins can convert 5-methylcytosine to 5-hydroxymethylcytosine raised the possibility that Tet proteins are capable of establishing a distinct epigenetic state. We have recently demonstrated that Tet1 is specifically expressed in murine embryonic stem (ES) cells and is required for ES cell maintenance. Using chromatin immunoprecipitation coupled with high-throughput DNA sequencing, here we show in mouse ES cells that Tet1 is preferentially bound to CpG-rich sequences at promoters of both transcriptionally active and Polycomb-repressed genes. Despite an increase in levels of DNA methylation at many Tet1-binding sites, Tet1 depletion does not lead to downregulation of all the Tet1 targets. Interestingly, although Tet1-mediated promoter hypomethylation is required for maintaining the expression of a group of transcriptionally active genes, it is also involved in repression of Polycomb-targeted developmental regulators. Tet1 contributes to silencing of this group of genes by facilitating recruitment of PRC2 to CpG-rich gene promoters. Thus, our study not only establishes a role for Tet1 in modulating DNA methylation levels at CpG-rich promoters, but also reveals a dual function of Tet1 in promoting transcription of pluripotency factors as well as participating in the repression of Polycomb-targeted developmental regulators.

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 analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells

Recent studies have demonstrated that the Ten-eleven translocation (Tet) family proteins can enzymatically convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). While 5mC has been studied extensively, little is known about the distribution and function of 5hmC. Here we present a genome-wide profile of 5hmC in mouse embryonic stem (ES) cells. A combined analysis of global 5hmC distribution and gene expression profile in wild-type and Tet1-depleted ES cells suggests that 5hmC is enriched at both gene bodies of actively transcribed genes and extended promoter regions of Polycomb-repressed developmental regulators. Thus, our study reveals the first genome-wide 5hmC distribution in pluripotent stem cells, and supports its dual function in regulating gene expression.

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.

Genome-wide detection of DNase I hypersensitive sites in single cells and FFPE tissue samples

DNase I hypersensitive sites (DHSs) provide important information on the presence of transcriptional regulatory elements and the state of chromatin in mammalian cells. Conventional DNase sequencing (DNase-seq) for genome-wide DHSs profiling is limited by the requirement of millions of cells. Here we report an ultrasensitive strategy, called single-cell DNase sequencing (scDNase-seq) for detection of genome-wide DHSs in single cells. We show that DHS patterns at the single-cell level are highly reproducible among individual cells. Among different single cells, highly expressed gene promoters and enhancers associated with multiple active histone modifications display constitutive DHS whereas chromatin regions with fewer histone modifications exhibit high variation of DHS. Furthermore, the single-cell DHSs predict enhancers that regulate cell-specific gene expression programs and the cell-to-cell variations of DHS are predictive of gene expression. Finally, we apply scDNase-seq to pools of tumour cells and pools of normal cells, dissected from formalin-fixed paraffin-embedded tissue slides from patients with thyroid cancer, and detect thousands of tumour-specific DHSs. Many of these DHSs are associated with promoters and enhancers critically involved in cancer development. Analysis of the DHS sequences uncovers one mutation (chr18: 52417839G>C) in the tumour cells of a patient with follicular thyroid carcinoma, which affects the binding of the tumour suppressor protein p53 and correlates with decreased expression of its target gene TXNL1. In conclusion, scDNase-seq can reliably detect DHSs in single cells, greatly extending the range of applications of DHS analysis both for basic and for translational research, and may provide critical information for personalized medicine.

Directional gene expression and antisense transcripts in sexual and asexual stages of Plasmodium falciparum

BACKGROUND

It has been shown that nearly a quarter of the initial predicted gene models in the Plasmodium falciparum genome contain errors. Although there have been efforts to obtain complete cDNA sequences to correct the errors, the coverage of cDNA sequences on the predicted genes is still incomplete, and many gene models for those expressed in sexual or mosquito stages have not been validated. Antisense transcripts have widely been reported in P. falciparum; however, the extent and pattern of antisense transcripts in different developmental stages remain largely unknown.

RESULTS

We have sequenced seven bidirectional libraries from ring, early and late trophozoite, schizont, gametocyte II, gametocyte V, and ookinete, and four strand-specific libraries from late trophozoite, schizont, gametocyte II, and gametocyte V of the 3D7 parasites. Alignment of the cDNA sequences to the 3D7 reference genome revealed stage-specific antisense transcripts and novel intron-exon splicing junctions. Sequencing of strand-specific cDNA libraries suggested that more genes are expressed in one direction in gametocyte than in schizont. Alternatively spliced genes, antisense transcripts, and stage-specific expressed genes were also characterized.

CONCLUSIONS

It is necessary to continue to sequence cDNA from different developmental stages, particularly those of non-erythrocytic stages. The presence of antisense transcripts in some gametocyte and ookinete genes suggests that these antisense RNA may play an important role in gene expression regulation and parasite development. Future gene expression studies should make use of directional cDNA libraries. Antisense transcripts may partly explain the observed discrepancy between levels of mRNA and protein expression.

H2A.Z facilitates access of active and repressive complexes to chromatin in embryonic stem cell self-renewal and differentiation

Chromatin modifications have been implicated in the self-renewal and differentiation of embryonic stem cells (ESCs). However, the function of histone variant H2A.Z in ESCs remains unclear. We show that H2A.Z is highly enriched at promoters and enhancers and is required for both efficient self-renewal and differentiation of murine ESCs. H2A.Z deposition leads to an abnormal nucleosome structure, decreased nucleosome occupancy, and increased chromatin accessibility. In self-renewing ESCs, knockdown of H2A.Z compromises OCT4 binding to its target genes and leads to decreased binding of MLL complexes to active genes and of PRC2 complex to repressed genes. During differentiation of ESCs, inhibition of H2A.Z also compromises RA-induced RARα binding, activation of differentiation markers, and the repression of pluripotency genes. We propose that H2A.Z mediates such contrasting activities by acting as a general facilitator that generates access for a variety of complexes, both activating and repressive.

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.

Overall survival in the OlympiA phase III trial of adjuvant olaparib in patients with germline pathogenic variants in BRCA1/2 and high risk, early breast cancer

BACKGROUND

The randomized, double-blind OlympiA trial compared 1 year of the oral poly(adenosine diphosphate-ribose) polymerase inhibitor, olaparib, to matching placebo as adjuvant therapy for patients with pathogenic or likely pathogenic variants in germline BRCA1 or BRCA2 (gBRCA1/2pv) and high-risk, human epidermal growth factor receptor 2-negative, early breast cancer (EBC). The first pre-specified interim analysis (IA) previously demonstrated statistically significant improvement in invasive disease-free survival (IDFS) and distant disease-free survival (DDFS). The olaparib group had fewer deaths than the placebo group, but the difference did not reach statistical significance for overall survival (OS). We now report the pre-specified second IA of OS with updates of IDFS, DDFS, and safety.

PATIENTS AND METHODS

One thousand eight hundred and thirty-six patients were randomly assigned to olaparib or placebo following (neo)adjuvant chemotherapy, surgery, and radiation therapy if indicated. Endocrine therapy was given concurrently with study medication for hormone receptor-positive cancers. Statistical significance for OS at this IA required P < 0.015.

RESULTS

With a median follow-up of 3.5 years, the second IA of OS demonstrated significant improvement in the olaparib group relative to the placebo group [hazard ratio 0.68; 98.5% confidence interval (CI) 0.47-0.97; P = 0.009]. Four-year OS was 89.8% in the olaparib group and 86.4% in the placebo group (Δ 3.4%, 95% CI -0.1% to 6.8%). Four-year IDFS for the olaparib group versus placebo group was 82.7% versus 75.4% (Δ 7.3%, 95% CI 3.0% to 11.5%) and 4-year DDFS was 86.5% versus 79.1% (Δ 7.4%, 95% CI 3.6% to 11.3%), respectively. Subset analyses for OS, IDFS, and DDFS demonstrated benefit across major subgroups. No new safety signals were identified including no new cases of acute myeloid leukemia or myelodysplastic syndrome.

CONCLUSION

With 3.5 years of median follow-up, OlympiA demonstrates statistically significant improvement in OS with adjuvant olaparib compared with placebo for gBRCA1/2pv-associated EBC and maintained improvements in the previously reported, statistically significant endpoints of IDFS and DDFS with no new safety signals.

BAF complexes facilitate decatenation of DNA by topoisomerase IIα

Recent exon-sequencing studies of human tumours have revealed that subunits of BAF (mammalian SWI/SNF) complexes are mutated in more than 20% of all human malignancies, but the mechanisms involved in tumour suppression are unclear. BAF chromatin-remodelling complexes are polymorphic assemblies that use energy provided by ATP hydrolysis to regulate transcription through the control of chromatin structure and the placement of Polycomb repressive complex 2 (PRC2) across the genome. Several proteins dedicated to this multisubunit complex, including BRG1 (also known as SMARCA4) and BAF250a (also known as ARID1A), are mutated at frequencies similar to those of recognized tumour suppressors. In particular, the core ATPase BRG1 is mutated in 5-10% of childhood medulloblastomas and more than 15% of Burkitt's lymphomas. Here we show a previously unknown function of BAF complexes in decatenating newly replicated sister chromatids, a requirement for proper chromosome segregation during mitosis. We find that deletion of Brg1 in mouse cells, as well as the expression of BRG1 point mutants identified in human tumours, leads to anaphase bridge formation (in which sister chromatids are linked by catenated strands of DNA) and a G2/M-phase block characteristic of the decatenation checkpoint. Endogenous BAF complexes interact directly with endogenous topoisomerase IIα (TOP2A) through BAF250a and are required for the binding of TOP2A to approximately 12,000 sites across the genome. Our results demonstrate that TOP2A chromatin binding is dependent on the ATPase activity of BRG1, which is compromised in oncogenic BRG1 mutants. These studies indicate that the ability of TOP2A to prevent DNA entanglement at mitosis requires BAF complexes and suggest that this activity contributes to the role of BAF subunits as tumour suppressors.

SMARCB1 is required for widespread BAF complex-mediated activation of enhancers and bivalent promoters

Perturbations to mammalian SWI/SNF (BAF) complexes contribute to over 20% of human cancers, with driving roles first identified in malignant rhabdoid tumor (MRT), an aggressive pediatric cancer characterized by biallelic inactivation of the core BAF complex subunit SMARCB1 (BAF47). However, the mechanism by which this alteration contributes to tumorigenesis remains poorly understood. We find that BAF47 loss destabilizes BAF complexes on chromatin, absent significant changes in intra-complex integrity. Rescue of BAF47 in BAF47-deficient sarcoma cell lines results in increased genome-wide BAF complex occupancy, facilitating widespread enhancer activation and opposition of polycomb-mediated repression at bivalent promoters. We demonstrate differential regulation by BAF and PBAF complexes at enhancers and promoters, respectively, suggesting distinct functions of each complex which are perturbed upon BAF47 loss. Our results demonstrate collaborative mechanisms of mSWI/SNF-mediated gene activation, identifying functions that are coopted or abated to drive human cancers and developmental disorders.

PBAF chromatin-remodeling complex requires a novel specificity subunit, BAF200, to regulate expression of selective interferon-responsive genes

PBAF and BAF are two chromatin-remodeling complexes of the SWI/SNF family essential for mammalian transcription and development. Although these complexes share eight identical subunits, only PBAF can facilitate transcriptional activation by nuclear receptors in vitro. Here we show that these complexes have selectivity in mediating transcription of different interferon-responsive genes. The selectivity by PBAF requires a novel subunit, BAF200, but not the previously described PBAF-specificity subunit, BAF180 (Polybromo). Our study provides in vivo evidence that PBAF and BAF regulate expression of distinct genes, and suggests that BAF200 plays a key role in PBAF function.

RNA Polymerase II Regulates Topoisomerase 1 Activity to Favor Efficient Transcription

We report a mechanism through which the transcription machinery directly controls topoisomerase 1 (TOP1) activity to adjust DNA topology throughout the transcription cycle. By comparing TOP1 occupancy using chromatin immunoprecipitation sequencing (ChIP-seq) versus TOP1 activity using topoisomerase 1 sequencing (TOP1-seq), a method reported here to map catalytically engaged TOP1, TOP1 bound at promoters was discovered to become fully active only after pause-release. This transition coupled the phosphorylation of the carboxyl-terminal-domain (CTD) of RNA polymerase II (RNAPII) with stimulation of TOP1 above its basal rate, enhancing its processivity. TOP1 stimulation is strongly dependent on the kinase activity of BRD4, a protein that phosphorylates Ser2-CTD and regulates RNAPII pause-release. Thus the coordinated action of BRD4 and TOP1 overcame the torsional stress opposing transcription as RNAPII commenced elongation but preserved negative supercoiling that assists promoter melting at start sites. This nexus between transcription and DNA topology promises to elicit new strategies to intercept pathological gene expression.

CTCF-Mediated Enhancer-Promoter Interaction Is a Critical Regulator of Cell-to-Cell Variation of Gene Expression

Recent studies indicate that even a homogeneous population of cells display heterogeneity in gene expression and response to environmental stimuli. Although promoter structure critically influences the cell-to-cell variation of gene expression in bacteria and lower eukaryotes, it remains unclear what controls the gene expression noise in mammals. Here we report that CTCF decreases cell-to-cell variation of expression by stabilizing enhancer-promoter interaction. We show that CTCF binding sites are interwoven with enhancers within topologically associated domains (TADs) and a positive correlation is found between CTCF binding and the activity of the associated enhancers. Deletion of CTCF sites compromises enhancer-promoter interactions. Using single-cell flow cytometry and single-molecule RNA-FISH assays, we demonstrate that knocking down of CTCF or deletion of a CTCF binding site results in increased cell-to-cell variation of gene expression, indicating that long-range promoter-enhancer interaction mediated by CTCF plays important roles in controlling the cell-to-cell variation of gene expression in mammalian cells.

High-resolution genome-wide mapping of histone modifications

The expression patterns of eukaryotic genomes are controlled by their chromatin structure, consisting of nucleosome subunits in which DNA of approximately 146 bp is wrapped around a core of 8 histone molecules. Post-translational histone modifications play an essential role in modifying chromatin structure. Here we apply a combination of SAGE and chromatin immunoprecipitation (ChIP) protocols to determine the distribution of hyperacetylated histones H3 and H4 in the Saccharomyces cerevisiae genome. We call this approach genome-wide mapping technique (GMAT). Using GMAT, we find that the highest acetylation levels are detected in the 5' end of a gene's coding region, but not in the promoter. Furthermore, we show that the histone acetyltransferase, GCN5p, regulates H3 acetylation in the promoter and 5' end of the coding regions. These findings indicate that GMAT should find valuable applications in mapping target sites of chromatin-modifying enzymes.

The transcription factors T-bet and Runx are required for the ontogeny of pathogenic interferon-γ-producing T helper 17 cells

T helper 17 (Th17) cells can give rise to interleukin-17A (IL-17A)- and interferon (IFN)-γ-double-producing cells that are implicated in development of autoimmune diseases. However, the molecular mechanisms that govern generation of IFN-γ-producing Th17 cells are unclear. We found that coexpression of the Th1 and Th17 cell master transcription factors, T-bet and retinoid-related orphan receptor gamma-t (RORγt), respectively, did not generate Th cells with robust IL-17 and IFN-γ expression. Instead, development of IFN-γ-producing Th17 cells required T-bet and Runx1 or Runx3. IL-12-stimulated Th17 cells upregulated Runx1, which bound to the Ifng locus in a T-bet-dependent manner. Reciprocally, T-bet or Runx1 deficiency or inhibition of Runx transcriptional activity impaired the development of IFN-γ-producing Th17 cells during experimental autoimmune encephalomyelitis, which correlated with substantially ameliorated disease course. Thus, our studies identify a critical role for T-bet and Runx transcription factors in the generation of pathogenic IFN-γ-producing Th17 cells.

Critical Role of STAT5 transcription factor tetramerization for cytokine responses and normal immune function

Cytokine-activated STAT proteins dimerize and bind to high-affinity motifs, and N-terminal domain-mediated oligomerization of dimers allows tetramer formation and binding to low-affinity tandem motifs, but the functions of dimers versus tetramers are unknown. We generated Stat5a-Stat5b double knockin (DKI) N-domain mutant mice in which STAT5 proteins form dimers but not tetramers, identified cytokine-regulated genes whose expression required STAT5 tetramers, and defined dimer versus tetramer consensus motifs. Whereas Stat5-deficient mice exhibited perinatal lethality, DKI mice were viable; thus, STAT5 dimers were sufficient for survival. Nevertheless, STAT5 DKI mice had fewer CD4(+)CD25(+) T cells, NK cells, and CD8(+) T cells, with impaired cytokine-induced and homeostatic proliferation of CD8(+) T cells. Moreover, DKI CD8(+) T cell proliferation after viral infection was diminished and DKI Treg cells did not efficiently control colitis. Thus, tetramerization of STAT5 is critical for cytokine responses and normal immune function, establishing a critical role for STAT5 tetramerization in vivo.