Mapping Genome-wide Accessible Chromatin in Primary Human T Lymphocytes by ATAC-Seq

Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method used for the identification of open (accessible) regions ...

Libraries and Information System of Bar-Ilan University

Bar-Ilan University

Higher-order genome organization shows tissue-specific patterns. However, functional relevance and the mechanisms shaping the genome architecture are poorly understood. Here we report a profound shift from promiscuous to highly selective genome organization that accompanies the effector lineage choice of differentiating T cells. As multipotent naive cells receive antigenic signals and commit to a T helper (Th) pathway, the genome-wide contacts of a lineage-specific cytokine locus are preferentially enriched for functionally relevant genes. Despite the establishment of divergent interactomes and global reprogramming of transcription in Th1 versus Th2, the overall expression status of the contact genes is surprisingly similar between the two lineages. Importantly, during differentiation, the genomic contacts are retained and strengthened precisely at DNA binding sites of the specific lineage-determining STAT transcription factor. In cells from the specific STAT knock-out mouse, the signature cytokine locus is unable to shed the promiscuous contacts established in the naive T cells, indicating the importance of genomic STAT binding. Altogether, the global aggregation of STAT binding loci from genic and nongenic regions highlights a new role for differentiation-promoting transcription factors in direct specification of higher-order nuclear architecture through interacting with regulatory regions. Such subnuclear environments have significant implications for efficient functioning of the mature effector lymphocytes.

Spatial congregation of STAT binding directs selective nuclear architecture during T-cell functional differentiation.

The "CTCF code" hypothesis posits that CTCF pleiotropic functions are driven by recognition of diverse sequences through combinatorial use of its 11 zinc fingers (ZFs). This model, however, is supported by in vitro binding studies of a limited number of sequences. To study CTCF multivalency in vivo, we define ZF binding requirements at ∼50,000 genomic sites in primary lymphocytes. We find that CTCF reads sequence diversity through ZF clustering. ZFs 4-7 anchor CTCF to ∼80% of targets containing the core motif. Nonconserved flanking sequences are recognized by ZFs 1-2 and ZFs 8-11 clusters, which also stabilize CTCF broadly. Alternatively, ZFs 9-11 associate with a second phylogenetically conserved upstream motif at ∼15% of its sites. Individually, ZFs increase overall binding and chromatin residence time. Unexpectedly, we also uncovered a conserved downstream DNA motif that destabilizes CTCF occupancy. Thus, CTCF associates with a wide array of DNA modules via combinatorial clustering of its 11 ZFs.

A genome-wide map of CTCF multivalency redefines the CTCF code.

A key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation.

The cell nucleus is a busy and organized organelle. In this megalopolis made of billions of nucleotides, protein factors find their target loci to exert nuclear functions such as transcription and replication. Remarkably, despite the lack of internal membrane barrier, the interlinked and tightly regulated nuclear processes occur in spatially organized fashion. These processes can lead to double-strand breaks (DSBs) that compromise the integrity of the genome. Moreover, in some cells like lymphocytes, DNA damage is also targeted within the context of immunoglobulin gene recombination. If not repaired correctly, DSBs can cause chromosomal rearrangements, including translocations which are etiological in numerous tumors. Therefore, the chromosomal locations of DSBs, as well as their spatial positioning, are important contributors to formation of chromosomal translocations at specific genomic loci. To obtain a mechanistic understanding of chromosomal translocations these parameters should be accounted for in a global and integrative fashion. In this review we will discuss recent findings addressing how genome architecture, DNA damage, and repair contribute to the genesis of chromosomal translocations.

3D view of chromosomes, DNA damage, and translocations.

Appropriate DNA double-strand break (DSB) repair factor choice is essential for ensuring accurate repair outcome and genomic integrity. The factors that regulate this process remain poorly understood. Here, we identify two repressive chromatin components, the macrohistone variant macroH2A1 and the H3K9 methyltransferase and tumor suppressor PRDM2, which together direct the choice between the antagonistic DSB repair mediators BRCA1 and 53BP1. The macroH2A1/PRDM2 module mediates an unexpected shift from accessible to condensed chromatin that requires the ataxia telangiectasia mutated (ATM)-dependent accumulation of both proteins at DSBs in order to promote DSB-flanking H3K9 dimethylation. Remarkably, loss of macroH2A1 or PRDM2, as well as experimentally induced chromatin decondensation, impairs the retention of BRCA1, but not 53BP1, at DSBs. As a result, macroH2A1 and/or PRDM2 depletion causes epistatic defects in DSB end resection, homology-directed repair, and the resistance to poly(ADP-ribose) polymerase (PARP) inhibition-all hallmarks of BRCA1-deficient tumors. Together, these findings identify dynamic, DSB-associated chromatin reorganization as a critical modulator of BRCA1-dependent genome maintenance.

A macrohistone variant links dynamic chromatin compaction to BRCA1-dependent genome maintenance.

The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA(+) enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

B cell super-enhancers and regulatory clusters recruit AID tumorigenic activity.

Although physiological steroid levels are often pulsatile (ultradian), the genomic effects of this pulsatility are poorly understood. By utilizing glucocorticoid receptor (GR) signaling as a model system, we uncovered striking spatiotemporal relationships between receptor loading, lifetimes of the DNase I hypersensitivity sites (DHSs), long-range interactions, and gene regulation. We found that hormone-induced DHSs were enriched within ± 50 kb of GR-responsive genes and displayed a broad spectrum of lifetimes upon hormone withdrawal. These lifetimes dictate the strength of the DHS interactions with gene targets and contribute to gene regulation from a distance. Our results demonstrate that pulsatile and constant hormone stimulations induce unique, treatment-specific patterns of gene and regulatory element activation. These modes of activation have implications for corticosteroid function in vivo and for steroid therapies in various clinical settings.

Dynamics of chromatin accessibility and long-range interactions in response to glucocorticoid pulsing.

Three-dimensional (3-D) genome organization in the nuclear space affects various genomic functions. Circular chromosome conformation capture (4C-seq) is a powerful technique that allows researchers to measure long-range chromosomal interactions with a locus of interest across the entire genome. This method relies on enzymatic cleavage of cross-linked chromatin and consecutive ligation to create ligation junctions between physically adjacent loci, followed by PCR amplification of locus-specific associating loci. The enzymes used must meet 4C standards because variations in their efficiency and performance may affect the quality of the obtained data. Here we systematically compare the efficiency and reliability of different T4 DNA ligases and PCR DNA polymerases, assessing the most critical and technically challenging steps in 4C. The results of this analysis enable the use of cost-effective enzymes with superior specificity and efficiency for 4C and save time in screening for appropriate primers. This information provides users with flexibility in their experimental design and guidelines for adapting and testing any enzyme of choice for obtaining standardized results.

Comparative analysis of T4 DNA ligases and DNA polymerases used in chromosome conformation capture assays.

Polycomb group (PcG) proteins are evolutionarily conserved chromatin modifiers that regulate developmental pathways in plants. PcGs form nuclear multi-subunit Polycomb Repressive Complexes (PRCs). The PRC2 complex mediates gene repression via methylation of lysine 27 on histone H3, which consequently leads to chromatin condensation. In Arabidopsis thaliana, several PRC2 complexes with different compositions were identified, each controlling a particular developmental program.The core subunit FIE is crucial for PRC2 function throughout the plant life cycle, yet accurate information on its spatial and temporal localization was absent. This study focused on identifying FIE accumulation patterns, using microscopy and biochemical approaches. Analysing endogenous FIE and transgenic gFIE-green fluorescent protein fusion protein (gFIE-GFP) showed that FIE accumulates in the nuclei of every cell type examined. Interestingly, gFIE-GFP, as well as the endogenous FIE, also localized to the cytoplasm in all examined tissues. In both vegetative and reproductive organs, FIE formed cytoplasmic high-molecular-mass complexes, in parallel to the nuclear PRC2 complexes. Moreover, size-exclusion chromatography and bimolecular fluorescence complementation assays indicated that in inflorescences FIE formed a cytoplasmic complex with MEA, a PRC2 histone methyltransferase subunit. In contrast, CLF and SWN histone methyltransferases were strictly nuclear. Presence of PRC2 subunits in cytoplasmic complexes has not been previously described in plants. Our findings are in agreement with accumulating evidence demonstrating cytoplasmic localization and function of PcGs in metazoa. The cytosolic accumulation of PRC2 components in plants supports the model that PcGs have alternative non-nuclear functions that go beyond chromatin methylation.

FIE, a nuclear PRC2 protein, forms cytoplasmic complexes in Arabidopsis thaliana.

CCCTC-binding factor (CTCF) is an organizer of higher-order chromatin structure and regulates gene expression. Genetic studies have implicated mutations in CTCF in intellectual disabilities. However, the role of CTCF-mediated chromatin structure in learning and memory is unclear. We show that depletion of CTCF in postmitotic neurons, or depletion in the hippocampus of adult mice through viral-mediated knockout, induces deficits in learning and memory. These deficits in learning and memory at the beginning of adulthood are correlated with impaired long-term potentiation and reduced spine density, with no changes in basal synaptic transmission and dendritic morphogenesis and arborization. Cognitive disabilities are associated with downregulation of cadherin and learning-related genes. In addition, CTCF knockdown attenuates fear-conditioning-induced hippocampal gene expression of key learning genes and loss of long-range interactions at the BDNF and Arc loci. This study thus suggests that CTCF-dependent gene expression regulation and genomic organization are regulators of learning and memory.

Neuronal CTCF Is Necessary for Basal and Experience-Dependent Gene Regulation, Memory Formation, and Genomic Structure of BDNF and Arc.

Glucocorticoids and their receptor (GR) have been an important area of research because of their pleiotropic physiological functions and extensive use in the clinic. In addition, the association between GR and glucocorticoids, which is highly specific, leads to rapid nuclear translocation where GR associates with chromatin to regulate gene transcription. This simplified model system has been instrumental for studying the complexity of transcription regulation processes occurring at chromatin. In this review we discuss our current understanding of GR action that has been enhanced by recent developments in genome wide measurements of chromatin accessibility, histone marks, chromatin remodeling and 3D chromatin structure in various cell types responding to glucocorticoids.

Genomic effects of glucocorticoids.

The three dimensional folding of mammalian genomes is cell type specific and difficult to alter suggesting that it is an important component of gene regulation. However, given the multitude of chromatin-associating factors, the mechanisms driving the colocalization of active chromosomal domains and the role of this organization in regulating the transcription program in adipocytes are not clear. Analysis of genome-wide chromosomal associations revealed cell type-specific spatial clustering of adipogenic genes in 3T3-L1 cells. Time course analysis demonstrated that the adipogenic 'hub', sampled by PPARγ and Lpin1, undergoes orchestrated reorganization during adipogenesis. Coupling the dynamics of genome architecture with multiple chromatin datasets indicated that among all the transcription factors (TFs) tested, RXR is central to genome reorganization at the beginning of adipogenesis. Interestingly, at the end of differentiation, the adipogenic hub was shifted to an H3K27me3-repressive environment in conjunction with attenuation of gene transcription. We propose a stage-specific hierarchy for the activity of TFs contributing to the establishment of an adipogenic genome architecture that brings together the adipogenic genetic program. In addition, the repositioning of this network in a H3K27me3-rich environment at the end of differentiation may contribute to the stabilization of gene transcription levels and reduce the developmental plasticity of these specialized cells.

Hierarchical role for transcription factors and chromatin structure in genome organization along adipogenesis.

Ectopic expression of lineage master regulators induces transdifferentiation. Whether cell fate transitions can be induced during various developmental stages has not been systemically examined. Here we discover that amongst different developmental stages, mouse embryonic stem cells (mESCs) are resistant to cell fate conversion induced by the melanocyte lineage master regulator MITF. By generating a transgenic system we exhibit that in mESCs, the pluripotency master regulator Oct4, counteracts pro-differentiation induced by Mitf by physical interference with MITF transcriptional activity. We further demonstrate that mESCs must be released from Oct4-maintained pluripotency prior to ectopically induced differentiation. Moreover, Oct4 induction in various differentiated cells represses their lineage identity in vivo. Alongside, chromatin architecture combined with ChIP-seq analysis suggest that Oct4 competes with various lineage master regulators for binding promoters and enhancers. Our analysis reveals pluripotency and transdifferentiation regulatory principles and could open new opportunities in the field of regenerative medicine.

Ofir Hakim

The Mina and Everard Goodman Faculty of Life Sciences

Libraries and Information System of Bar-Ilan University

Mapping Genome-wide Accessible Chromatin in Primary Human T Lymphocytes by ATAC-Seq

Ofir Hakim

.css-o250i3{font-size:18px;font-family:Open Sans,sans-serif;line-height:21.6px;}The Mina and Everard Goodman Faculty of Life Sciences

Libraries and Information System of Bar-Ilan University

Mapping Genome-wide Accessible Chromatin in Primary Human T Lymphocytes by ATAC-Seq

The Mina and Everard Goodman Faculty of Life Sciences