Cell of origin epigenetic priming determines susceptibility to Tet2 mutation.

Hematopoietic stem cell (HSC) mutations can result in clonal hematopoiesis (CH) with heterogeneous clinical outcomes. Here, we investigate how the cell state preceding Tet2 mutation impacts the pre-malignant phenotype. Using an inducible system for clonal analysis of myeloid progenitors, we find that the epigenetic features of clones at similar differentiation status are highly heterogeneous and functionally respond differently to Tet2 mutation. Cell differentiation stage also influences Tet2 mutation response indicating that the cell of origin's epigenome modulates clone-specific behaviors in CH. Molecular features associated with higher risk outcomes include Sox4 that sensitizes cells to Tet2 inactivation, inducing dedifferentiation, altered metabolism and increasing the in vivo clonal output of mutant cells, as confirmed in primary GMP and HSC models. Our findings validate the hypothesis that epigenetic features can predispose specific clones for dominance, explaining why identical genetic mutations can result in different phenotypes.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements.

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

Kinetic networks identify TWIST2 as a key regulatory node in adipogenesis.

Adipocytes contribute to metabolic disorders such as obesity, diabetes, and atherosclerosis. Prior characterizations of the transcriptional network driving adipogenesis have overlooked transiently acting transcription factors (TFs), genes, and regulatory elements that are essential for proper differentiation. Moreover, traditional gene regulatory networks provide neither mechanistic details about individual regulatory element-gene relationships nor temporal information needed to define a regulatory hierarchy that prioritizes key regulatory factors. To address these shortcomings, we integrate kinetic chromatin accessibility (ATAC-seq) and nascent transcription (PRO-seq) data to generate temporally resolved networks that describe TF binding events and resultant effects on target gene expression. Our data indicate which TF families cooperate with and antagonize each other to regulate adipogenesis. Compartment modeling of RNA polymerase density quantifies how individual TFs mechanistically contribute to distinct steps in transcription. The glucocorticoid receptor activates transcription by inducing RNA polymerase pause release, whereas SP and AP-1 factors affect RNA polymerase initiation. We identify Twist2 as a previously unappreciated effector of adipocyte differentiation. We find that TWIST2 acts as a negative regulator of 3T3-L1 and primary preadipocyte differentiation. We confirm that Twist2 knockout mice have compromised lipid storage within subcutaneous and brown adipose tissue. Previous phenotyping of Twist2 knockout mice and Setleis syndrome Twist2 -/- patients noted deficiencies in subcutaneous adipose tissue. This network inference framework is a powerful and general approach for interpreting complex biological phenomena and can be applied to a wide range of cellular processes.

Functional inference of gene regulation using single-cell multi-omics.

Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human blood cells. Collectively, we generate ~91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scA-TAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

The Drosophila BEAF insulator protein interacts with the polybromo subunit of the PBAP chromatin remodeling complex.

The Drosophila Boundary Element-Associated Factor of 32 kDa (BEAF) binds in promoter regions of a few thousand mostly housekeeping genes. BEAF is implicated in both chromatin domain boundary activity and promoter function, although molecular mechanisms remain elusive. Here, we show that BEAF physically interacts with the polybromo subunit (Pbro) of PBAP, a SWI/SNF-class chromatin remodeling complex. BEAF also shows genetic interactions with Pbro and other PBAP subunits. We examine the effect of this interaction on gene expression and chromatin structure using precision run-on sequencing and micrococcal nuclease sequencing after RNAi-mediated knockdown in cultured S2 cells. Our results are consistent with the interaction playing a subtle role in gene activation. Fewer than 5% of BEAF-associated genes were significantly affected after BEAF knockdown. Most were downregulated, accompanied by fill-in of the promoter nucleosome-depleted region and a slight upstream shift of the +1 nucleosome. Pbro knockdown caused downregulation of several hundred genes and showed a correlation with BEAF knockdown but a better correlation with promoter-proximal GAGA factor binding. Micrococcal nuclease sequencing supports that BEAF binds near housekeeping gene promoters while Pbro is more important at regulated genes. Yet there is a similar general but slight reduction of promoter-proximal pausing by RNA polymerase II and increase in nucleosome-depleted region nucleosome occupancy after knockdown of either protein. We discuss the possibility of redundant factors keeping BEAF-associated promoters active and masking the role of interactions between BEAF and the Pbro subunit of PBAP in S2 cells. We identify Facilitates Chromatin Transcription (FACT) and Nucleosome Remodeling Factor (NURF) as candidate redundant factors.

Deep learning-based enhancement of epigenomics data with AtacWorks.

ATAC-seq is a widely-applied assay used to measure genome-wide chromatin accessibility; however, its ability to detect active regulatory regions can depend on the depth of sequencing coverage and the signal-to-noise ratio. Here we introduce AtacWorks, a deep learning toolkit to denoise sequencing coverage and identify regulatory peaks at base-pair resolution from low cell count, low-coverage, or low-quality ATAC-seq data. Models trained by AtacWorks can detect peaks from cell types not seen in the training data, and are generalizable across diverse sample preparations and experimental platforms. We demonstrate that AtacWorks enhances the sensitivity of single-cell experiments by producing results on par with those of conventional methods using ~10 times as many cells, and further show that this framework can be adapted to enable cross-modality inference of protein-DNA interactions. Finally, we establish that AtacWorks can enable new biological discoveries by identifying active regulatory regions associated with lineage priming in rare subpopulations of hematopoietic stem cells.

Pioneer-like factor GAF cooperates with PBAP (SWI/SNF) and NURF (ISWI) to regulate transcription.

Transcriptionally silent genes must be activated throughout development. This requires nucleosomes be removed from promoters and enhancers to allow transcription factor (TF) binding and recruitment of coactivators and RNA polymerase II (Pol II). Specialized pioneer TFs bind nucleosome-wrapped DNA to perform this chromatin opening by mechanisms that remain incompletely understood. Here, we show that GAGA factor (GAF), a Drosophila pioneer-like factor, functions with both SWI/SNF and ISWI family chromatin remodelers to allow recruitment of Pol II and entry to a promoter-proximal paused state, and also to promote Pol II's transition to productive elongation. We found that GAF interacts with PBAP (SWI/SNF) to open chromatin and allow Pol II to be recruited. Importantly, this activity is not dependent on NURF as previously proposed; however, GAF also synergizes with NURF downstream from this process to ensure efficient Pol II pause release and transition to productive elongation, apparently through its role in precisely positioning the +1 nucleosome. These results demonstrate how a single sequence-specific pioneer TF can synergize with remodelers to activate sets of genes. Furthermore, this behavior of remodelers is consistent with findings in yeast and mice, and likely represents general, conserved mechanisms found throughout eukarya.

Epigenomic State Transitions Characterize Tumor Progression in Mouse Lung Adenocarcinoma.

Regulatory networks that maintain functional, differentiated cell states are often dysregulated in tumor development. Here, we use single-cell epigenomics to profile chromatin state transitions in a mouse model of lung adenocarcinoma (LUAD). We identify an epigenomic continuum representing loss of cellular identity and progression toward a metastatic state. We define co-accessible regulatory programs and infer key activating and repressive chromatin regulators of these cell states. Among these co-accessibility programs, we identify a pre-metastatic transition, characterized by activation of RUNX transcription factors, which mediates extracellular matrix remodeling to promote metastasis and is predictive of survival across human LUAD patients. Together, these results demonstrate the power of single-cell epigenomics to identify regulatory programs to uncover mechanisms and key biomarkers of tumor progression.

Inference and effects of barcode multiplets in droplet-based single-cell assays.

A widespread assumption for single-cell analyses specifies that one cell's nucleic acids are predominantly captured by one oligonucleotide barcode. Here, we show that ~13-21% of cell barcodes from the 10x Chromium scATAC-seq assay may have been derived from a droplet with more than one oligonucleotide sequence, which we call "barcode multiplets". We demonstrate that barcode multiplets can be derived from at least two different sources. First, we confirm that approximately 4% of droplets from the 10x platform may contain multiple beads. Additionally, we find that approximately 5% of beads may contain detectable levels of multiple oligonucleotide barcodes. We show that this artifact can confound single-cell analyses, including the interpretation of clonal diversity and proliferation of intra-tumor lymphocytes. Overall, our work provides a conceptual and computational framework to identify and assess the impacts of barcode multiplets in single-cell data.

Defining data-driven primary transcript annotations with primaryTranscriptAnnotation in R.

SUMMARY: Nascent transcript measurements derived from run-on sequencing experiments are critical for the investigation of transcriptional mechanisms and regulatory networks. However, conventional mRNA gene annotations significantly differ from the boundaries of primary transcripts. New primary transcript annotations are needed to accurately interpret run-on data. We developed the primaryTranscriptAnnotation R package to infer the transcriptional start and termination sites of primary transcripts from genomic run-on data. We then used these inferred coordinates to annotate transcriptional units identified de novo. This package provides the novel utility to integrate data-driven primary transcript annotations with transcriptional unit coordinates identified in an unbiased manner. Highlighting the importance of using accurate primary transcript coordinates, we demonstrate that this new methodology increases the detection of differentially expressed transcripts and provides more accurate quantification of RNA polymerase pause indices. AVAILABILITY AND IMPLEMENTATION: https://github.com/WarrenDavidAnderson/genomicsRpackage/tree/master/primaryTranscriptAnnotation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An improved auxin-inducible degron system preserves native protein levels and enables rapid and specific protein depletion.

Rapid perturbation of protein function permits the ability to define primary molecular responses while avoiding downstream cumulative effects of protein dysregulation. The auxin-inducible degron (AID) system was developed as a tool to achieve rapid and inducible protein degradation in nonplant systems. However, tagging proteins at their endogenous loci results in chronic auxin-independent degradation by the proteasome. To correct this deficiency, we expressed the auxin response transcription factor (ARF) in an improved inducible degron system. ARF is absent from previously engineered AID systems but is a critical component of native auxin signaling. In plants, ARF directly interacts with AID in the absence of auxin, and we found that expression of the ARF PB1 (Phox and Bem1) domain suppresses constitutive degradation of AID-tagged proteins. Moreover, the rate of auxin-induced AID degradation is substantially faster in the ARF-AID system. To test the ARF-AID system in a quantitative and sensitive manner, we measured genome-wide changes in nascent transcription after rapidly depleting the ZNF143 transcription factor. Transcriptional profiling indicates that ZNF143 activates transcription in cis and regulates promoter-proximal paused RNA polymerase density. Rapidly inducible degradation systems that preserve the target protein's native expression levels and patterns will revolutionize the study of biological systems by enabling specific and temporally defined protein dysregulation.

Droplet-based combinatorial indexing for massive-scale single-cell chromatin accessibility.

Recent technical advancements have facilitated the mapping of epigenomes at single-cell resolution; however, the throughput and quality of these methods have limited their widespread adoption. Here we describe a high-quality (105 nuclear fragments per cell) droplet-microfluidics-based method for single-cell profiling of chromatin accessibility. We use this approach, named 'droplet single-cell assay for transposase-accessible chromatin using sequencing' (dscATAC-seq), to assay 46,653 cells for the unbiased discovery of cell types and regulatory elements in adult mouse brain. We further increase the throughput of this platform by combining it with combinatorial indexing (dsciATAC-seq), enabling single-cell studies at a massive scale. We demonstrate the utility of this approach by measuring chromatin accessibility across 136,463 resting and stimulated human bone marrow-derived cells to reveal changes in the cis- and trans-regulatory landscape across cell types and under stimulatory conditions at single-cell resolution. Altogether, we describe a total of 510,123 single-cell profiles, demonstrating the scalability and flexibility of this droplet-based platform.

Molecular mechanisms driving transcriptional stress responses.

Proteotoxic stress, that is, stress caused by protein misfolding and aggregation, triggers the rapid and global reprogramming of transcription at genes and enhancers. Genome-wide assays that track transcriptionally engaged RNA polymerase II (Pol II) at nucleotide resolution have provided key insights into the underlying molecular mechanisms that regulate transcriptional responses to stress. In addition, recent kinetic analyses of transcriptional control under heat stress have shown how cells 'prewire' and rapidly execute genome-wide changes in transcription while concurrently becoming poised for recovery. The regulation of Pol II at genes and enhancers in response to heat stress is coupled to chromatin modification and compartmentalization, as well as to co-transcriptional RNA processing. These mechanistic features seem to apply broadly to other coordinated genome-regulatory responses.

Transcription factors GAF and HSF act at distinct regulatory steps to modulate stress-induced gene activation.

The coordinated regulation of gene expression at the transcriptional level is fundamental to development and homeostasis. Inducible systems are invaluable when studying transcription because the regulatory process can be triggered instantaneously, allowing the tracking of ordered mechanistic events. Here, we use precision run-on sequencing (PRO-seq) to examine the genome-wide heat shock (HS) response in Drosophila and the function of two key transcription factors on the immediate transcription activation or repression of all genes regulated by HS. We identify the primary HS response genes and the rate-limiting steps in the transcription cycle that GAGA-associated factor (GAF) and HS factor (HSF) regulate. We demonstrate that GAF acts upstream of promoter-proximally paused RNA polymerase II (Pol II) formation (likely at the step of chromatin opening) and that GAF-facilitated Pol II pausing is critical for HS activation. In contrast, HSF is dispensable for establishing or maintaining Pol II pausing but is critical for the release of paused Pol II into the gene body at a subset of highly activated genes. Additionally, HSF has no detectable role in the rapid HS repression of thousands of genes.

Mammalian Heat Shock Response and Mechanisms Underlying Its Genome-wide Transcriptional Regulation.

The heat shock response (HSR) is critical for survival of all organisms. However, its scope, extent, and the molecular mechanism of regulation are poorly understood. Here we show that the genome-wide transcriptional response to heat shock in mammals is rapid and dynamic and results in induction of several hundred and repression of several thousand genes. Heat shock factor 1 (HSF1), the "master regulator" of the HSR, controls only a fraction of heat shock-induced genes and does so by increasing RNA polymerase II release from promoter-proximal pause. Notably, HSF2 does not compensate for the lack of HSF1. However, serum response factor appears to transiently induce cytoskeletal genes independently of HSF1. The pervasive repression of transcription is predominantly HSF1-independent and is mediated through reduction of RNA polymerase II pause release. Overall, mammalian cells orchestrate rapid, dynamic, and extensive changes in transcription upon heat shock that are largely modulated at pause release, and HSF1 plays a limited and specialized role.

High-throughput binding characterization of RNA aptamer selections using a microplate-based multiplex microcolumn device.

We describe a versatile 96-well microplate-based device that utilizes affinity microcolumn chromatography to complement downstream plate-based processing in aptamer selections. This device is reconfigurable and is able to operate in serial and/or parallel mode with up to 96 microcolumns. We demonstrate the utility of this device by simultaneously performing characterizations of target binding using five RNA aptamers and a random library. This was accomplished through 96 total selection tests. Three sets of selections tested the effects of target concentration on aptamer binding compared to the random RNA library using aptamers to the proteins green fluorescent protein (GFP), human heat shock factor 1 (hHSF1), and negative elongation factor E (NELF-E). For all three targets, we found significant effects consistent with steric hindrance with optimum enrichments at predictable target concentrations. In a fourth selection set, we tested the partitioning efficiency and binding specificity of our three proteins' aptamers, as well as two suspected background binding sequences, to eight targets running serially. The targets included an empty microcolumn, three affinity resins, three specific proteins, and a non-specific protein control. The aptamers showed significant enrichments only on their intended targets. Specifically, the hHSF1 and NELF-E aptamers enriched over 200-fold on their protein targets, and the GFP aptamer enriched 750-fold. By utilizing our device's plate-based format with other complementary plate-based systems for all downstream biochemical processes and analysis, high-throughput selections, characterizations, and optimization were performed to significantly reduce the time and cost for completing large-scale aptamer selections.

Multiplexed microcolumn-based process for efficient selection of RNA aptamers.

We describe a reusable microcolumn and process for the efficient discovery of nucleic acid aptamers for multiple target molecules. The design of our device requires only microliter volumes of affinity chromatography resin-a condition that maximizes the enrichment of target-binding sequences over non-target-binding (i.e., background) sequences. Furthermore, the modular design of the device accommodates a multiplex aptamer selection protocol. We optimized the selection process performance using microcolumns filled with green fluorescent protein (GFP)-immobilized resin and monitoring, over a wide range of experimental conditions, the enrichment of a known GFP-binding RNA aptamer (GFPapt) against a random RNA aptamer library. We validated the multiplex approach by monitoring the enrichment of GFPapt in de novo selection experiments with GFP and other protein preparations. After only three rounds of selection, the cumulative GFPapt enrichment on the GFP-loaded resin was greater than 10(8) with no enrichment for the other nonspecific targets. We used this optimized protocol to perform a multiplex selection to two human heat shock factor (hHSF) proteins, hHSF1 and hHSF2. High-throughput sequencing was used to identify aptamers for each protein that were preferentially enriched in just three selection rounds, which were confirmed and isolated after five rounds. Gel-shift and fluorescence polarization assays showed that each aptamer binds with high-affinity (KD < 20 nM) to the respective targets. The combination of our microcolumns with a multiplex approach and high-throughput sequencing enables the selection of aptamers to multiple targets in a high-throughput and efficient manner.

Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production.

Bioethanol is a biofuel produced mainly from the fermentation of carbohydrates derived from agricultural feedstocks by the yeast Saccharomyces cerevisiae. One of the most widely adopted strains is PE-2, a heterothallic diploid naturally adapted to the sugar cane fermentation process used in Brazil. Here we report the molecular genetic analysis of a PE-2 derived diploid (JAY270), and the complete genome sequence of a haploid derivative (JAY291). The JAY270 genome is highly heterozygous (approximately 2 SNPs/kb) and has several structural polymorphisms between homologous chromosomes. These chromosomal rearrangements are confined to the peripheral regions of the chromosomes, with breakpoints within repetitive DNA sequences. Despite its complex karyotype, this diploid, when sporulated, had a high frequency of viable spores. Hybrid diploids formed by outcrossing with the laboratory strain S288c also displayed good spore viability. Thus, the rearrangements that exist near the ends of chromosomes do not impair meiosis, as they do not span regions that contain essential genes. This observation is consistent with a model in which the peripheral regions of chromosomes represent plastic domains of the genome that are free to recombine ectopically and experiment with alternative structures. We also explored features of the JAY270 and JAY291 genomes that help explain their high adaptation to industrial environments, exhibiting desirable phenotypes such as high ethanol and cell mass production and high temperature and oxidative stress tolerance. The genomic manipulation of such strains could enable the creation of a new generation of industrial organisms, ideally suited for use as delivery vehicles for future bioenergy technologies.