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.

Recognition of a mononucleosomal histone modification pattern by BPTF via multivalent interactions.

Little is known about how combinations of histone marks are interpreted at the level of nucleosomes. The second PHD finger of human BPTF is known to specifically recognize histone H3 when methylated on lysine 4 (H3K4me2/3). Here, we examine how additional heterotypic modifications influence BPTF binding. Using peptide surrogates, three acetyllysine ligands are indentified for a PHD-adjacent bromodomain in BPTF via systematic screening and biophysical characterization. Although the bromodomain displays limited discrimination among the three possible acetyllysines at the peptide level, marked selectivity is observed for only one of these sites, H4K16ac, in combination with H3K4me3 at the mononucleosome level. In support, these two histone marks constitute a unique trans-histone modification pattern that unambiguously resides within a single nucleosomal unit in human cells, and this module colocalizes with these marks in the genome. Together, our data call attention to nucleosomal patterning of covalent marks in dictating critical chromatin associations.

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.

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.

Calibrating ChIP-Seq with Nucleosomal Internal Standards to Measure Histone Modification Density Genome Wide.

Chromatin immunoprecipitation (ChIP) serves as a central experimental technique in epigenetics research, yet there are serious drawbacks: it is a relative measurement, which untethered to any external scale obscures fair comparison among experiments; it employs antibody reagents that have differing affinities and specificities for target epitopes that vary in abundance; and it is frequently not reproducible. To address these problems, we developed Internal Standard Calibrated ChIP (ICeChIP), wherein a native chromatin sample is spiked with nucleosomes reconstituted from recombinant and semisynthetic histones on barcoded DNA prior to immunoprecipitation. ICeChIP measures local histone modification densities on a biologically meaningful scale, enabling unbiased trans-experimental comparisons, and reveals unique insight into the nature of bivalent domains. This technology provides in situ assessment of the immunoprecipitation step, accommodating for many experimental pitfalls as well as providing a critical examination of untested assumptions inherent to conventional ChIP.

An Interactive Database for the Assessment of Histone Antibody Specificity.

Access to high-quality antibodies is a necessity for the study of histones and their posttranslational modifications (PTMs). Here we debut the Histone Antibody Specificity Database (http://www.histoneantibodies.com), an online and expanding resource cataloging the behavior of widely used, commercially available histone antibodies by peptide microarray. This interactive web portal provides a critical resource to the biological research community that routinely uses these antibodies as detection reagents for a wide range of applications.

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.

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.

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.

Multivalent engagement of chromatin modifications by linked binding modules.

Various chemical modifications on histones and regions of associated DNA play crucial roles in genome management by binding specific factors that, in turn, serve to alter the structural properties of chromatin. These so-called effector proteins have typically been studied with the biochemist's paring knife--the capacity to recognize specific chromatin modifications has been mapped to an increasing number of domains that frequently appear in the nuclear subset of the proteome, often present in large, multisubunit complexes that bristle with modification-dependent binding potential. We propose that multivalent interactions on a single histone tail and beyond may have a significant, if not dominant, role in chromatin transactions.

The United States of histone ubiquitylation and methylation.

In this issue of Molecular Cell, Wu et al. (2011) reveal that ubiquitylation of histone 2B lysine 34 stimulates histone methyltransferase activity on nucleosomes, a finding with implications for the general mechanism by which monoubiquitylation may influence subsequent modification activities.

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.

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.

Multiple interactions recruit MLL1 and MLL1 fusion proteins to the HOXA9 locus in leukemogenesis.

MLL1 fusion proteins activate HoxA9 gene expression and cause aggressive leukemias that respond poorly to treatment, but how they recognize and stably bind to HoxA9 is not clearly understood. In a systematic analysis of MLL1 domain recruitment activity, we identified an essential MLL1 recruitment domain that includes the CXXC domain and PHD fingers and is controlled by direct interactions with the PAF elongation complex and H3K4Me2/3. MLL1 fusion proteins lack the PHD fingers and require prebinding of a wild-type MLL1 complex and CXXC domain recognition of DNA for stable HoxA9 association. Together, these results suggest that specific recruitment of MLL1 requires multiple interactions and is a precondition for stable recruitment of MLL1 fusion proteins to HoxA9 in leukemogenesis. Since wild-type MLL1 and oncogenic MLL1 fusion proteins have overlapping yet distinct recruitment mechanisms, this creates a window of opportunity that could be exploited for the development of targeted therapies.

Recombinant antibodies to histone post-translational modifications.

Variability in the quality of antibodies to histone post-translational modifications (PTMs) is a widely recognized hindrance in epigenetics research. Here, we produced recombinant antibodies to the trimethylated lysine residues of histone H3 with high specificity and affinity and no lot-to-lot variation. These recombinant antibodies performed well in common epigenetics applications, and enabled us to identify positive and negative correlations among histone PTMs.

Traceless semisynthesis of a set of histone 3 species bearing specific lysine methylation marks.

Considerable mechanistic insight into the function of histone post-translational modifications and the enzymes that install and remove them derives from in vitro experiments with modified histones, often embedded in nucleosomes. We report the first semisyntheses of native-like histone 3 (H3) bearing tri- and dimethyllysines at position 79 and trimethyllysine at position 36, as well as more facile and traceless semisyntheses of K9 and K27 trimethylated species. These semisyntheses are practical on a multi-milligram scale and can also generate H3 with combinations of marks. Each of these modifications has distinct functional consequences, although the pathways by which H3K36me3 and H3K79me2/3 act have not been entirely mapped. To this end, we demonstrated that our semisynthetic histones, when reconstituted into nucleosomes, are valuable affinity reagents for unbiased binding partner discovery and compare them to their methyllysine analogue (MLA) counterparts at the nucleosome level.

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.

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.

Crystal structure of Staphylococcus aureus tRNA adenosine deaminase TadA in complex with RNA.

Bacterial tRNA adenosine deaminases (TadAs) catalyze the hydrolytic deamination of adenosine to inosine at the wobble position of tRNA(Arg2), a process that enables this single tRNA to recognize three different arginine codons in mRNA. In addition, inosine is also introduced at the wobble position of multiple eukaryotic tRNAs. The genes encoding these deaminases are essential in bacteria and yeast, demonstrating the importance of their biological activity. Here we report the crystallization and structure determination to 2.0 A of Staphylococcus aureus TadA bound to the anticodon stem-loop of tRNA(Arg2) bearing nebularine, a non-hydrolyzable adenosine analog, at the wobble position. The cocrystal structure reveals the basis for both sequence and structure specificity in the interactions of TadA with RNA, and it additionally provides insight into the active site architecture that promotes efficient hydrolytic deamination.

Regulation of MLL1 H3K4 methyltransferase activity by its core components.

Histone H3 Lys4 (H3K4) methylation is a prevalent mark associated with transcription activation. A common feature of several H3K4 methyltransferase complexes is the presence of three structural components (RbBP5, Ash2L and WDR5) and a catalytic subunit containing a SET domain. Here we report the first biochemical reconstitution of a functional four-component mixed-lineage leukemia protein-1 (MLL1) core complex. This reconstitution, combined with in vivo assays, allows direct analysis of the contribution of each component to MLL1 enzymatic activity and their roles in transcriptional regulation. Moreover, taking clues from a crystal structure analysis, we demonstrate that WDR5 mediates interactions of the MLL1 catalytic unit both with the common structural platform and with the histone substrate. Mechanistic insights gained from this study can be generalized to the whole family of SET1-like histone methyltransferases in mammals.