Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet.

BACKGROUND: Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult. RESULTS: Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPß for methylated motifs and the preference of c-Myc for unmethylated E-box motifs. CONCLUSIONS: Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9.

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites. SIGNIFICANCE: The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

An oligodendrocyte silencer element underlies the pathogenic impact of lamin B1 structural variants.

The role of non-coding regulatory elements and how they might contribute to tissue type specificity of disease phenotypes is poorly understood. Autosomal Dominant Leukodystrophy (ADLD) is a fatal, adult-onset, neurological disorder that is characterized by extensive CNS demyelination. Most cases of ADLD are caused by tandem genomic duplications involving the lamin B1 gene ( LMNB1 ) while a small subset are caused by genomic deletions upstream of the gene. Utilizing data from recently identified families that carry LMNB1 gene duplications but do not exhibit demyelination, ADLD patient tissues, CRISPR modified cell lines and mouse models, we have identified a novel silencer element that is lost in ADLD patients and that specifically targets overexpression to oligodendrocytes. This element consists of CTCF binding sites that mediate three-dimensional chromatin looping involving the LMNB1 and the recruitment of the PRC2 repressor complex. Loss of the silencer element in ADLD identifies a previously unknown role for silencer elements in tissue specificity and disease causation.

FACT regulates pluripotency through proximal and distal regulation of gene expression in murine embryonic stem cells.

BACKGROUND: The FACT complex is a conserved histone chaperone with critical roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. RESULTS: To explore the mechanism for FACT function, we generated AID-tagged murine embryonic cell lines for FACT subunit SPT16 and paired depletion with nascent transcription and chromatin accessibility analyses. We also analyzed SPT16 occupancy using CUT&RUN and found that SPT16 localizes to both promoter and enhancer elements, with a strong overlap in binding with OCT4, SOX2, and NANOG. Over a timecourse of SPT16 depletion, nucleosomes invade new loci, including promoters, regions bound by SPT16, OCT4, SOX2, and NANOG, and TSS-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (encoding OCT4), Sox2, Nanog, and enhancer RNAs produced from these genes' associated enhancers are downregulated. CONCLUSIONS: We propose that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency.

Single-Cell Factor Localization on Chromatin using Ultra-Low Input Cleavage Under Targets and Release using Nuclease.

Determining the binding locations of a protein on chromatin is essential for understanding its function and potential regulatory targets. Chromatin Immunoprecipitation (ChIP) has been the gold standard for determining protein localization for over 30 years and is defined by the use of an antibody to pull out the protein of interest from sonicated or enzymatically digested chromatin. More recently, antibody tethering techniques have become popular for assessing protein localization on chromatin due to their increased sensitivity. Cleavage Under Targets & Release Under Nuclease (CUT&RUN) is the genome-wide derivative of Chromatin Immunocleavage (ChIC) and utilizes recombinant Protein A tethered to micrococcal nuclease (pA-MNase) to identify the IgG constant region of the antibody targeting a protein of interest, therefore enabling site-specific cleavage of the DNA flanking the protein of interest. CUT&RUN can be used to profile histone modifications, transcription factors, and other chromatin-binding proteins such as nucleosome remodeling factors. Importantly, CUT&RUN can be used to assess the localization of either euchromatic- or heterochromatic-associated proteins and histone modifications. For these reasons, CUT&RUN is a powerful method for determining the binding profiles of a wide range of proteins. Recently, CUT&RUN has been optimized for transcription factor profiling in low populations of cells and single cells and the optimized protocol has been termed ultra-low input CUT&RUN (uliCUT&RUN). Here, a detailed protocol is presented for single-cell factor profiling using uliCUT&RUN in a manual 96-well format.

Super-resolution imaging reveals the evolution of higher-order chromatin folding in early carcinogenesis.

Genomic DNA is folded into a higher-order structure that regulates transcription and maintains genomic stability. Although progress has been made on understanding biochemical characteristics of epigenetic modifications in cancer, the in-situ higher-order folding of chromatin structure during malignant transformation remains largely unknown. Here, using optimized stochastic optical reconstruction microscopy (STORM) for pathological tissue (PathSTORM), we uncover a gradual decompaction and fragmentation of higher-order chromatin folding throughout all stages of carcinogenesis in multiple tumor types, and prior to tumor formation. Our integrated imaging, genomic, and transcriptomic analyses reveal functional consequences in enhanced transcription activities and impaired genomic stability. We also demonstrate the potential of imaging higher-order chromatin disruption to detect high-risk precursors that cannot be distinguished by conventional pathology. Taken together, our findings reveal gradual decompaction and fragmentation of higher-order chromatin structure as an enabling characteristic in early carcinogenesis to facilitate malignant transformation, which may improve cancer diagnosis, risk stratification, and prevention.