Nour J. Abdulhay

Predicting molecular function directly from DNA sequence remains a grand challenge in computational and molecular biology. Here, we engineer and train bidirectional transformer models to predict the chemical grammar of alternative human mRNA splicing leveraging the largest perturbative full-length RNA dataset to date. By combining high-throughput single-molecule long-read “chemical transcriptomics” in human cells with transformer models, we train AllSplice – a nucleotide foundation model that achieves state-of-the-art prediction of canonical and noncanonical splice junctions across the human transcriptome. We demonstrate improved performance achieved through incorporation of diverse noncanonical splice sites in its training set that were identified through long-read RNA data. Leveraging chemical perturbations and multiple cell types in the data, we fine-tune AllSplice to train ChemSplice – the first predictive model of sequence-dependent and cell-type specific alternative splicing following programmed cellular perturbation. We anticipate the broad application of AllSplice, ChemSplice, and other models fine-tuned on this foundation to myriad areas of RNA therapeutics development.

Nearly all essential nuclear processes act on DNA packaged into arrays of nucleosomes. However, our understanding of how these processes (for example, DNA replication, RNA transcription, chromatin extrusion and nucleosome remodeling) occur on individual chromatin arrays remains unresolved. Here, to address this deficit, we present SAMOSA-ChAAT: a massively multiplex single-molecule footprinting approach to map the primary structure of individual, reconstituted chromatin templates subject to virtually any chromatin-associated reaction. We apply this method to distinguish between competing models for chromatin remodeling by the essential imitation switch (ISWI) ATPase SNF2h: nucleosome-density-dependent spacing versus fixed-linker-length nucleosome clamping. First, we perform in vivo single-molecule nucleosome footprinting in murine embryonic stem cells, to discover that ISWI-catalyzed nucleosome spacing correlates with the underlying nucleosome density of specific epigenomic domains. To establish causality, we apply SAMOSA-ChAAT to quantify the activities of ISWI ATPase SNF2h and its parent complex ACF on reconstituted nucleosomal arrays of varying nucleosome density, at single-molecule resolution. We demonstrate that ISWI remodelers operate as density-dependent, length-sensing nucleosome sliders, whose ability to program DNA accessibility is dictated by single-molecule nucleosome density. We propose that the long-observed, context-specific regulatory effects of ISWI complexes can be explained in part by the sensing of nucleosome density within epigenomic domains. More generally, our approach promises molecule-precise views of the essential processes that shape nuclear physiology.

Our understanding of the beads-on-a-string arrangement of nucleosomes has been built largely on high-resolution sequence-agnostic imaging methods and sequence-resolved bulk biochemical techniques. To bridge the divide between these approaches, we present the single-molecule adenine methylated oligonucleosome sequencing assay (SAMOSA). SAMOSA is a high-throughput single-molecule sequencing method that combines adenine methyltransferase footprinting and single-molecule real-time DNA sequencing to natively and nondestructively measure nucleosome positions on individual chromatin fibres. SAMOSA data allows unbiased classification of single-molecular ‘states’ of nucleosome occupancy on individual chromatin fibres. We leverage this to estimate nucleosome regularity and spacing on single chromatin fibres genome-wide, at predicted transcription factor binding motifs, and across both active and silent human epigenomic domains. Our analyses suggest that chromatin is comprised of a diverse array of both regular and irregular single-molecular oligonucleosome patterns that differ subtly in their relative abundance across epigenomic domains. This irregularity is particularly striking in constitutive heterochromatin, which has typically been viewed as a conformationally static entity. Our proof-of-concept study provides a powerful new methodology for studying nucleosome organization at a previously intractable resolution, and offers up new avenues for modeling and visualizing higher-order chromatin structure. 1-sentence summary: High-throughput single-molecule real-time footprinting of chromatin arrays reveals heterogeneous patterns of oligonucleosome occupancy.

Studies of allelic variation underlying genetic blood disorders have provided important insights into human hematopoiesis. Most often, the identified pathogenic mutations result in loss-of-function or missense changes. However, assessing the pathogenicity of noncoding variants can be challenging. Here, we characterize two unrelated patients with a distinct presentation of dyserythropoietic anemia and other impairments in hematopoiesis associated with an intronic mutation in GATA1 that is 24 nucleotides upstream of the canonical splice acceptor site. Functional studies demonstrate that this single-nucleotide alteration leads to reduced canonical splicing and increased use of an alternative splice acceptor site that causes a partial intron retention event. The resultant altered GATA1 contains a five–amino acid insertion at the C-terminus of the C-terminal zinc finger and has no observable activity. Collectively, our results demonstrate how altered splicing of GATA1, which reduces levels of the normal form of this master transcription factor, can result in distinct changes in human hematopoiesis.

Immunodeficient mouse models have been valuable for studies of human hematopoiesis, but high‐fidelity recapitulation of erythropoiesis in most xenograft recipients remains elusive. Recently developed immunodeficient and Kit mutant mice, however, have provided a suitable background to achieve higher‐level human erythropoiesis after long‐term hematopoietic engraftment. While there has been some characterization of human erythropoiesis in these models, a comprehensive analysis from various human developmental stages has not yet been reported. Here, we have utilized cell surface phenotypes, morphologic analyses, and molecular studies to fully characterize human erythropoiesis from multiple developmental stages in immunodeficient and Kit mutant mouse models following long‐term hematopoietic stem and progenitor cell engraftment. We show that human erythropoiesis in such models demonstrates complete maturation and enucleation, as well as developmentally appropriate globin gene expression. These results provide a framework for future studies to utilize this model system for interrogating disorders affecting human erythropoiesis and for developing improved therapeutic approaches.