Guiding the HBO1 complex function through the JADE subunit.

JADE is a core subunit of the HBO1 acetyltransferase complex that regulates developmental and epigenetic programs and promotes gene transcription. Here we describe the mechanism by which JADE facilitates recruitment of the HBO1 complex to chromatin and mediates its enzymatic activity. Structural, genomic and complex assembly in vivo studies show that the PZP (PHD1-zinc-knuckle-PHD2) domain of JADE engages the nucleosome through binding to histone H3 and DNA and is necessary for the association with chromatin targets. Recognition of unmethylated H3K4 by PZP directs enzymatic activity of the complex toward histone H4 acetylation, whereas H3K4 hypermethylation alters histone substrate selectivity. We demonstrate that PZP contributes to leukemogenesis, augmenting transforming activity of the NUP98-JADE2 fusion. Our findings highlight biological consequences and the impact of the intact JADE subunit on genomic recruitment, enzymatic function and pathological activity of the HBO1 complex.

Merging short and stranded long reads improves transcript assembly.

Long-read RNA sequencing has arisen as a counterpart to short-read sequencing, with the potential to capture full-length isoforms, albeit at the cost of lower depth. Yet this potential is not fully realized due to inherent limitations of current long-read assembly methods and underdeveloped approaches to integrate short-read data. Here, we critically compare the existing methods and develop a new integrative approach to characterize a particularly challenging pool of low-abundance long noncoding RNA (lncRNA) transcripts from short- and long-read sequencing in two distinct cell lines. Our analysis reveals severe limitations in each of the sequencing platforms. For short-read assemblies, coverage declines at transcript termini resulting in ambiguous ends, and uneven low coverage results in segmentation of a single transcript into multiple transcripts. Conversely, long-read sequencing libraries lack depth and strand-of-origin information in cDNA-based methods, culminating in erroneous assembly and quantitation of transcripts. We also discover a cDNA synthesis artifact in long-read datasets that markedly impacts the identity and quantitation of assembled transcripts. Towards remediating these problems, we develop a computational pipeline to "strand" long-read cDNA libraries that rectifies inaccurate mapping and assembly of long-read transcripts. Leveraging the strengths of each platform and our computational stranding, we also present and benchmark a hybrid assembly approach that drastically increases the sensitivity and accuracy of full-length transcript assembly on the correct strand and improves detection of biological features of the transcriptome. When applied to a challenging set of under-annotated and cell-type variable lncRNA, our method resolves the segmentation problem of short-read sequencing and the depth problem of long-read sequencing, resulting in the assembly of coherent transcripts with precise 5' and 3' ends. Our workflow can be applied to existing datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.

Specificity Guides Interpretation: On H3K4 Methylation at Enhancers and Broad Promoters.

In 2018, we used internally calibrated chromatin immunoprecipitation (ICeChIP) to find that many of the most commonly used antibodies against H3K4 methylforms had significant off-target binding, which compromised the findings of at least eight literature paradigms that used these antibodies for ChIP-seq (Shah et al., 2018). In many cases, we were able to recapitulate the prior findings in K562 cells with the original, low-quality antibody, only to find that the models did not hold up to scrutiny with highly specific reagents and quantitative calibration. In a recent preprint originally prepared as a Letter to the Editor of Molecular Cell, though they agree with our overarching conclusions, Pekowska and colleagues take issue with analyses presented for two relatively minor points of the paper (Pekowska et al., 2023). We are puzzled by the assertion that these two points constitute the "bulk" of our findings, nor is it clear which components of our "analytical design" they find problematic. We feel their critique, however mild, is misguided.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing.

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Non-canonical H3K79me2-dependent pathways promote the survival of MLL-rearranged leukemia.

MLL-rearranged leukemia depends on H3K79 methylation. Depletion of this transcriptionally activating mark by DOT1L deletion or high concentrations of the inhibitor pinometostat downregulates HOXA9 and MEIS1, and consequently reduces leukemia survival. Yet, some MLL-rearranged leukemias are inexplicably susceptible to low-dose pinometostat, far below concentrations that downregulate this canonical proliferation pathway. In this context, we define alternative proliferation pathways that more directly derive from H3K79me2 loss. By ICeChIP-seq, H3K79me2 is markedly depleted at pinometostat-downregulated and MLL-fusion targets, with paradoxical increases of H3K4me3 and loss of H3K27me3. Although downregulation of polycomb components accounts for some of the proliferation defect, transcriptional downregulation of FLT3 is the major pathway. Loss-of-FLT3-function recapitulates the cytotoxicity and gene expression consequences of low-dose pinometostat, whereas overexpression of constitutively active STAT5A, a target of FLT3-ITD-signaling, largely rescues these defects. This pathway also depends on MLL1, indicating combinations of DOT1L, MLL1 and FLT3 inhibitors should be explored for treating FLT3-mutant leukemia.

Sequence deeper without sequencing more: Bayesian resolution of ambiguously mapped reads.

Next-generation sequencing (NGS) has transformed molecular biology and contributed to many seminal insights into genomic regulation and function. Apart from whole-genome sequencing, an NGS workflow involves alignment of the sequencing reads to the genome of study, after which the resulting alignments can be used for downstream analyses. However, alignment is complicated by the repetitive sequences; many reads align to more than one genomic locus, with 15-30% of the genome not being uniquely mappable by short-read NGS. This problem is typically addressed by discarding reads that do not uniquely map to the genome, but this practice can lead to systematic distortion of the data. Previous studies that developed methods for handling ambiguously mapped reads were often of limited applicability or were computationally intensive, hindering their broader usage. In this work, we present SmartMap: an algorithm that augments industry-standard aligners to enable usage of ambiguously mapped reads by assigning weights to each alignment with Bayesian analysis of the read distribution and alignment quality. SmartMap is computationally efficient, utilizing far fewer weighting iterations than previously thought necessary to process alignments and, as such, analyzing more than a billion alignments of NGS reads in approximately one hour on a desktop PC. By applying SmartMap to peak-type NGS data, including MNase-seq, ChIP-seq, and ATAC-seq in three organisms, we can increase read depth by up to 53% and increase the mapped proportion of the genome by up to 18% compared to analyses utilizing only uniquely mapped reads. We further show that SmartMap enables the analysis of more than 140,000 repetitive elements that could not be analyzed by traditional ChIP-seq workflows, and we utilize this method to gain insight into the epigenetic regulation of different classes of repetitive elements. These data emphasize both the dangers of discarding ambiguously mapped reads and their power for driving biological discovery.

Epigenetic homogeneity in histone methylation underlies sperm programming for embryonic transcription.

Sperm contributes genetic and epigenetic information to the embryo to efficiently support development. However, the mechanism underlying such developmental competence remains elusive. Here, we investigated whether all sperm cells have a common epigenetic configuration that primes transcriptional program for embryonic development. Using calibrated ChIP-seq, we show that remodelling of histones during spermiogenesis results in the retention of methylated histone H3 at the same genomic location in most sperm cell. This homogeneously methylated fraction of histone H3 in the sperm genome is maintained during early embryonic replication. Such methylated histone fraction resisting post-fertilisation reprogramming marks developmental genes whose expression is perturbed upon experimental reduction of histone methylation. A similar homogeneously methylated histone H3 fraction is detected in human sperm. Altogether, we uncover a conserved mechanism of paternal epigenetic information transmission to the embryo through the homogeneous retention of methylated histone in a sperm cells population.

Chromatin-enriched RNAs mark active and repressive cis-regulation: An analysis of nuclear RNA-seq.

Long noncoding RNAs (lncRNAs) localize in the cell nucleus and influence gene expression through a variety of molecular mechanisms. Chromatin-enriched RNAs (cheRNAs) are a unique class of lncRNAs that are tightly bound to chromatin and putatively function to locally cis-activate gene transcription. CheRNAs can be identified by biochemical fractionation of nuclear RNA followed by RNA sequencing, but until now, a rigorous analytic pipeline for nuclear RNA-seq has been lacking. In this study, we survey four computational strategies for nuclear RNA-seq data analysis and develop a new pipeline, Tuxedo-ch, which outperforms other approaches. Tuxedo-ch assembles a more complete transcriptome and identifies cheRNA with higher accuracy than other approaches. We used Tuxedo-ch to analyze benchmark datasets of K562 cells and further characterize the genomic features of intergenic cheRNA (icheRNA) and their similarity to enhancer RNAs (eRNAs). We quantify the transcriptional correlation of icheRNA and adjacent genes and show that icheRNA is more positively associated with neighboring gene expression than eRNA or cap analysis of gene expression (CAGE) signals. We also explore two novel genomic associations of cheRNA, which indicate that cheRNAs may function to promote or repress gene expression in a context-dependent manner. IcheRNA loci with significant levels of H3K9me3 modifications are associated with active enhancers, consistent with the hypothesis that enhancers are derived from ancient mobile elements. In contrast, antisense cheRNA (as-cheRNA) may play a role in local gene repression, possibly through local RNA:DNA:DNA triple-helix formation.

Native internally calibrated chromatin immunoprecipitation for quantitative studies of histone post-translational modifications.

Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) has served as the central method for the study of histone modifications for the past decade. In ChIP-seq analyses, antibodies selectively capture nucleosomes bearing a modification of interest and the associated DNA is then mapped to the genome to determine the distribution of the mark. This approach has several important drawbacks: (i) ChIP interpretation necessitates the assumption of perfect antibody specificity, despite growing evidence that this is often not the case. (ii) Common methods for evaluating antibody specificity in other formats have little or no bearing on specificity within a ChIP experiment. (iii) Uncalibrated ChIP is reported as relative enrichment, which is biologically meaningless outside the experimental reference frame defined by a discrete immunoprecipitation (IP), thus preventing facile comparison across experimental conditions or modifications. (iv) Differential library amplification and loading onto next-generation sequencers, as well as computational normalization, can further compromise quantitative relationships that may exist between samples. Consequently, the researcher is presented with a series of potential pitfalls and is blind to nearly all of them. Here we provide a detailed protocol for internally calibrated ChIP (ICeChIP), a method we recently developed to resolve these problems by spike-in of defined nucleosomal standards within a ChIP procedure. This protocol is optimized for specificity and quantitative power, allowing for measurement of antibody specificity and absolute measurement of histone modification density (HMD) at genomic loci on a biologically meaningful scale enabling unambiguous comparisons. We provide guidance on optimal conditions for next-generation sequencing (NGS) and instructions for data analysis. This protocol takes between 17 and 18 h, excluding time for sequencing or bioinformatic analysis. The ICeChIP procedure enables accurate measurement of histone post-translational modifications (PTMs) genome-wide in mammalian cells as well as Drosophila melanogaster and Caenorhabditis elegans, indicating suitability for use in eukaryotic cells more broadly.

A Mutation in Histone H2B Represents a New Class of Oncogenic Driver.

By examination of the cancer genomics database, we identified a new set of mutations in core histones that frequently recur in cancer patient samples and are predicted to disrupt nucleosome stability. In support of this idea, we characterized a glutamate to lysine mutation of histone H2B at amino acid 76 (H2B-E76K), found particularly in bladder and head and neck cancers, that disrupts the interaction between H2B and H4. Although H2B-E76K forms dimers with H2A, it does not form stable histone octamers with H3 and H4 in vitro, and when reconstituted with DNA forms unstable nucleosomes with increased sensitivity to nuclease. Expression of the equivalent H2B mutant in yeast restricted growth at high temperature and led to defective nucleosome-mediated gene repression. Significantly, H2B-E76K expression in the normal mammary epithelial cell line MCF10A increased cellular proliferation, cooperated with mutant PIK3CA to promote colony formation, and caused a significant drift in gene expression and fundamental changes in chromatin accessibility, particularly at gene regulatory elements. Taken together, these data demonstrate that mutations in the globular domains of core histones may give rise to an oncogenic program due to nucleosome dysfunction and deregulation of gene expression. SIGNIFICANCE: Mutations in the core histones frequently occur in cancer and represent a new mechanism of epigenetic dysfunction that involves destabilization of the nucleosome, deregulation of chromatin accessibility, and alteration of gene expression to drive cellular transformation.See related commentary by Sarthy and Henikoff, p. 1346.This article is highlighted in the In This Issue feature, p. 1325.

Examining the Roles of H3K4 Methylation States with Systematically Characterized Antibodies.

Histone post-translational modifications (PTMs) are important genomic regulators often studied by chromatin immunoprecipitation (ChIP), whereby their locations and relative abundance are inferred by antibody capture of nucleosomes and associated DNA. However, the specificity of antibodies within these experiments has not been systematically studied. Here, we use histone peptide arrays and internally calibrated ChIP (ICeChIP) to characterize 52 commercial antibodies purported to distinguish the H3K4 methylforms (me1, me2, and me3, with each ascribed distinct biological functions). We find that many widely used antibodies poorly distinguish the methylforms and that high- and low-specificity reagents can yield dramatically different biological interpretations, resulting in substantial divergence from the literature for numerous H3K4 methylform paradigms. Using ICeChIP, we also discern quantitative relationships between enhancer H3K4 methylation and promoter transcriptional output and can measure global PTM abundance changes. Our results illustrate how poor antibody specificity contributes to the "reproducibility crisis," demonstrating the need for rigorous, platform-appropriate validation.

Regulated Capture of Vκ Gene Topologically Associating Domains by Transcription Factories.

Expression of vast repertoires of antigen receptors by lymphocytes, with each cell expressing a single receptor, requires stochastic activation of individual variable (V) genes for transcription and recombination. How this occurs remains unknown. Using single-cell RNA sequencing (scRNA-seq) and allelic variation, we show that individual pre-B cells monoallelically transcribe divergent arrays of Vκ genes, thereby opening stochastic repertoires for subsequent Vκ-Jκ recombination. Transcription occurs upon translocation of Vκ genes to RNA polymerase II arrayed on the nuclear matrix in transcription factories. Transcription is anchored by CTCF-bound sites or E2A-loaded Vκ promotors and continues over large genomic distances delimited only by topological associating domains (TADs). Prior to their monoallelic activation, Vκ loci are transcriptionally repressed by cyclin D3, which prevents capture of Vκ gene containing TADs by transcription factories. Cyclin D3 also represses protocadherin, olfactory, and other monoallelically expressed genes, suggesting a widely deployed mechanism for coupling monoallelic gene activation with cell cycle exit.

Targeting the MTF2-MDM2 Axis Sensitizes Refractory Acute Myeloid Leukemia to Chemotherapy.

Deep sequencing has revealed that epigenetic modifiers are the most mutated genes in acute myeloid leukemia (AML). Thus, elucidating epigenetic dysregulation in AML is crucial to understand disease mechanisms. Here, we demonstrate that metal response element binding transcription factor 2/polycomblike 2 (MTF2/PCL2) plays a fundamental role in the polycomb repressive complex 2 (PRC2) and that its loss elicits an altered epigenetic state underlying refractory AML. Unbiased systems analyses identified the loss of MTF2-PRC2 repression of MDM2 as central to, and therefore a biomarker for, refractory AML. Thus, immature MTF2-deficient CD34+CD38- cells overexpress MDM2, thereby inhibiting p53 that leads to chemoresistance due to defects in cell-cycle regulation and apoptosis. Targeting this dysregulated signaling pathway by MTF2 overexpression or MDM2 inhibitors sensitized refractory patient leukemic cells to induction chemotherapeutics and prevented relapse in AML patient-derived xenograft mice. Therefore, we have uncovered a direct epigenetic mechanism by which MTF2 functions as a tumor suppressor required for AML chemotherapeutic sensitivity and identified a potential therapeutic strategy to treat refractory AML.Significance: MTF2 deficiency predicts refractory AML at diagnosis. MTF2 represses MDM2 in hematopoietic cells and its loss in AML results in chemoresistance. Inhibiting p53 degradation by overexpressing MTF2 in vitro or by using MDM2 inhibitors in vivo sensitizes MTF2-deficient refractory AML cells to a standard induction-chemotherapy regimen. Cancer Discov; 8(11); 1376-89. ©2018 AACR. See related commentary by Duy and Melnick, p. 1348 This article is highlighted in the In This Issue feature, p. 1333.

Facile target validation in an animal model with intracellularly expressed monobodies.

Rapidly determining the biological effect of perturbing a site within a potential drug target could guide drug discovery efforts, but it remains challenging. Here, we describe a facile target validation approach that exploits monobodies, small synthetic binding proteins that can be fully functionally expressed in cells. We developed a potent and selective monobody to WDR5, a core component of the mixed lineage leukemia (MLL) methyltransferase complex. The monobody bound to the MLL interaction site of WDR5, the same binding site for small-molecule inhibitors whose efficacy has been demonstrated in cells but not in animals. As a genetically encoded reagent, the monobody inhibited proliferation of an MLL-AF9 cell line in vitro, suppressed its leukemogenesis and conferred a survival benefit in an in vivo mouse leukemia model. The capacity of this approach to readily bridge biochemical, structural, cellular characterization and tests in animal models may accelerate discovery and validation of druggable sites.

Quantitative and Structural Assessment of Histone Methyllysine Analogue Engagement by Cognate Binding Proteins Reveals Affinity Decrements Relative to Those of Native Counterparts.

Methyllysine analogues (MLAs), furnished by aminoethylation of engineered cysteine residues, are widely used surrogates of histone methyllysine and are considered to be effective proxies for studying these epigenetic marks in vitro. Here we report the first structure of a trimethyllysine MLA histone in complex with a protein binding partner, quantify the thermodynamic distinctions between MLAs and their native methyllysine counterparts, and demonstrate that these differences can compromise qualitative interpretations of binding at the nucleosome level. Quantitative measurements with two methyllysine binding protein modules reveal substantial affinity losses for the MLA peptides versus the corresponding native methyllysine species in both cases, although the thermodynamic underpinnings are distinct. MLA and methyllysine adopt distinct conformational geometries when in complex with the BPTF PHD finger, a well-established H3K4me3 binding partner. In this case, an ∼13-fold Kd difference at the peptide level translates to nucleosomal affinities for MLA analogues that fall outside of the detectable range in a pull-down format, whereas the methyllysine species installed by native chemical ligation demonstrates robust binding. Thus, despite their facile production and commercial availability, there is a significant caveat of potentially altered binding affinity when MLAs are used in place of native methyllysine residues.

Transcription-factor-dependent enhancer transcription defines a gene regulatory network for cardiac rhythm.

The noncoding genome is pervasively transcribed. Noncoding RNAs (ncRNAs) generated from enhancers have been proposed as a general facet of enhancer function and some have been shown to be required for enhancer activity. Here we examine the transcription-factor-(TF)-dependence of ncRNA expression to define enhancers and enhancer-associated ncRNAs that are involved in a TF-dependent regulatory network. TBX5, a cardiac TF, regulates a network of cardiac channel genes to maintain cardiac rhythm. We deep sequenced wildtype and Tbx5-mutant mouse atria, identifying ~2600 novel Tbx5-dependent ncRNAs. Tbx5-dependent ncRNAs were enriched for tissue-specific marks of active enhancers genome-wide. Tbx5-dependent ncRNAs emanated from regions that are enriched for TBX5-binding and that demonstrated Tbx5-dependent enhancer activity. Tbx5-dependent ncRNA transcription provided a quantitative metric of Tbx5-dependent enhancer activity, correlating with target gene expression. We identified RACER, a novel Tbx5-dependent long noncoding RNA (lncRNA) required for the expression of the calcium-handling gene Ryr2. We illustrate that TF-dependent enhancer transcription can illuminate components of TF-dependent gene regulatory networks.

Chromatin-enriched lncRNAs can act as cell-type specific activators of proximal gene transcription.

We recently described a new class of long noncoding RNAs (lncRNAs) that are distinguished by especially tight chromatin association and whose presence is strongly correlated to expression of nearby genes. Here, we examine the cis-enhancer mechanism of this class of chromatin-enriched RNA (cheRNA) across multiple human cell lines. cheRNAs are largely cell type specific and provide the most reliable chromatin signature to predict cis-gene transcription in every human cell type examined. Targeted depletion of three cheRNAs decreases expression of their neighboring genes, indicating potential co-activator function, and single-molecule fluorescence in situ hybridization (smFISH) of one cheRNA-distal target gene pair suggests a spatial overlap consistent with a role in chromosome looping. Additionally, the cheRNA HIDALGO stimulates the fetal hemoglobin subunit gamma 1 (HBG1) gene during erythroid differentiation by promoting contacts to a downstream enhancer. Our results suggest that multiple cheRNAs activate proximal lineage-specific gene transcription.

Antigen clasping by two antigen-binding sites of an exceptionally specific antibody for histone methylation.

Antibodies have a well-established modular architecture wherein the antigen-binding site residing in the antigen-binding fragment (Fab or Fv) is an autonomous and complete unit for antigen recognition. Here, we describe antibodies departing from this paradigm. We developed recombinant antibodies to trimethylated lysine residues on histone H3, important epigenetic marks and challenging targets for molecular recognition. Quantitative characterization demonstrated their exquisite specificity and high affinity, and they performed well in common epigenetics applications. Surprisingly, crystal structures and biophysical analyses revealed that two antigen-binding sites of these antibodies form a head-to-head dimer and cooperatively recognize the antigen in the dimer interface. This "antigen clasping" produced an expansive interface where trimethylated Lys bound to an unusually extensive aromatic cage in one Fab and the histone N terminus to a pocket in the other, thereby rationalizing the high specificity. A long-neck antibody format with a long linker between the antigen-binding module and the Fc region facilitated antigen clasping and achieved both high specificity and high potency. Antigen clasping substantially expands the paradigm of antibody-antigen recognition and suggests a strategy for developing extremely specific antibodies.

Nuclear Fractionation Reveals Thousands of Chromatin-Tethered Noncoding RNAs Adjacent to Active Genes.

A number of long noncoding RNAs (lncRNAs) have been reported to regulate transcription via recruitment of chromatin modifiers or bridging distal enhancer elements to gene promoters. However, the generality of these modes of regulation and the mechanisms of chromatin attachment for thousands of unstudied human lncRNAs remain unclear. To address these questions, we performed stringent nuclear fractionation coupled to RNA sequencing. We provide genome-wide identification of human chromatin-associated lncRNAs and demonstrate tethering of RNA to chromatin by RNAPII is a pervasive mechanism of attachment. We also uncovered thousands of chromatin-enriched RNAs (cheRNAs) that share molecular properties with known lncRNAs. Although distinct from eRNAs derived from active prototypical enhancers, the production of cheRNAs is strongly correlated with the expression of neighboring protein-coding genes. This work provides an updated framework for nuclear RNA organization that includes a large chromatin-associated transcript population correlated with active genes and may prove useful in de novo enhancer annotation.