High-capacity sample multiplexing for single cell chromatin accessibility profiling.

Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

The NELF pausing checkpoint mediates the functional divergence of Cdk9.

Promoter-proximal pausing by RNA Pol II is a rate-determining step in gene transcription that is hypothesized to be a prominent point at which regulatory factors act. The pausing factor NELF is known to induce and stabilize pausing, but not all kinds of pausing are NELF-mediated. Here, we find that NELF-depleted Drosophila melanogaster cells functionally recapitulate the NELF-independent pausing we previously observed in fission yeast (which lack NELF). Critically, only NELF-mediated pausing establishes a strict requirement for Cdk9 kinase activity for the release of paused Pol II into productive elongation. Upon inhibition of Cdk9, cells with NELF efficiently shutdown gene transcription, while in NELF-depleted cells, defective, non-productive transcription continues unabated. By introducing a strict checkpoint for Cdk9, the evolution of NELF was likely critical to enable increased regulation of Cdk9 in higher eukaryotes, as Cdk9 availability can be restricted to limit gene transcription without inducing wasteful, non-productive transcription.

Transcriptional Activation of Regenerative Hematopoiesis via Vascular Niche Sensing.

Transition between activation and quiescence programs in hematopoietic stem and progenitor cells (HSC/HSPCs) is perceived to be governed intrinsically and by microenvironmental co-adaptation. However, HSC programs dictating both transition and adaptability, remain poorly defined. Single cell multiome analysis divulging differential transcriptional activity between distinct HSPC states, indicated for the exclusive absence of Fli-1 motif from quiescent HSCs. We reveal that Fli-1 activity is essential for HSCs during regenerative hematopoiesis. Fli-1 directs activation programs while manipulating cellular sensory and output machineries, enabling HSPCs co-adoptability with a stimulated vascular niche. During regenerative conditions, Fli-1 presets and enables propagation of niche-derived Notch1 signaling. Constitutively induced Notch1 signaling is sufficient to recuperate functional HSC impairments in the absence of Fli-1. Applying FLI-1 modified-mRNA transduction into lethargic adult human mobilized HSPCs, enables their vigorous niche-mediated expansion along with superior engraftment capacities. Thus, decryption of stem cell activation programs offers valuable insights for immune regenerative medicine.

High-Capacity Sample Multiplexing for Single Cell Chromatin Accessibility Profiling.

Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

RNA polymerase mapping in plants identifies intergenic regulatory elements enriched in causal variants.

Control of gene expression is fundamental at every level of cell function. Promoter-proximal pausing and divergent transcription at promoters and enhancers, which are prominent features in animals, have only been studied in a handful of research experiments in plants. PRO-Seq analysis in cassava (Manihot esculenta) identified peaks of transcriptionally engaged RNA polymerase at both the 5' and 3' end of genes, consistent with paused or slowly moving Polymerase. In addition, we identified divergent transcription at intergenic sites. A full genome search for bi-directional transcription using an algorithm for enhancer detection developed in mammals (dREG) identified many intergenic regulatory element (IRE) candidates. These sites showed distinct patterns of methylation and nucleotide conservation based on genomic evolutionary rate profiling (GERP). SNPs within these IRE candidates explained significantly more variation in fitness and root composition than SNPs in chromosomal segments randomly ascertained from the same intergenic distribution, strongly suggesting a functional importance of these sites. Maize GRO-Seq data showed RNA polymerase occupancy at IREs consistent with patterns in cassava. Furthermore, these IREs in maize significantly overlapped with sites previously identified on the basis of open chromatin, histone marks, and methylation, and were enriched for reported eQTL. Our results suggest that bidirectional transcription can identify intergenic genomic regions in plants that play an important role in transcription regulation and whose identification has the potential to aid crop improvement.

Identifying Transcription Error-Enriched Genomic Loci Using Nuclear Run-on Circular-Sequencing Coupled with Background Error Modeling.

RNA polymerase transcribes certain genomic loci with higher errors rates. These transcription error-enriched genomic loci (TEELs) have implications in disease. Current deep-sequencing methods cannot distinguish TEELs from post-transcriptional modifications, stochastic transcription errors, and technical noise, impeding efforts to elucidate the mechanisms linking TEELs to disease. Here, we describe background error model-coupled precision nuclear run-on circular-sequencing (EmPC-seq) to discern genomic regions enriched for transcription misincorporations. EmPC-seq innovatively combines a nuclear run-on assay for capturing nascent RNA before post-transcriptional modifications, a circular-sequencing step that sequences the same nascent RNA molecules multiple times to improve accuracy, and a statistical model for distinguishing error-enriched regions among stochastic polymerase errors. Applying EmPC-seq to the ribosomal RNA transcriptome, we show that TEELs of RNA polymerase I are not randomly distributed but clustered together, with higher error frequencies at nascent transcript 3' ends. Our study establishes a reliable method of identifying TEELs with nucleotide precision, which can help elucidate their molecular origins.

Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme.

The human genome encodes a variety of poorly understood RNA species that remain challenging to identify using existing genomic tools. We developed chromatin run-on and sequencing (ChRO-seq) to map the location of RNA polymerase for almost any input sample, including samples with degraded RNA that are intractable to RNA sequencing. We used ChRO-seq to map nascent transcription in primary human glioblastoma (GBM) brain tumors. Enhancers identified in primary GBMs resemble open chromatin in the normal human brain. Rare enhancers that are activated in malignant tissue drive regulatory programs similar to the developing nervous system. We identified enhancers that regulate groups of genes that are characteristic of each known GBM subtype and transcription factors that drive them. Finally we discovered a core group of transcription factors that control the expression of genes associated with clinical outcomes. This study characterizes the transcriptional landscape of GBM and introduces ChRO-seq as a method to map regulatory programs that contribute to complex diseases.

A Cdk9-PP1 switch regulates the elongation-termination transition of RNA polymerase II.

The end of the RNA polymerase II (Pol II) transcription cycle is strictly regulated to prevent interference between neighbouring genes and to safeguard transcriptome integrity 1 . The accumulation of Pol II downstream of the cleavage and polyadenylation signal can facilitate the recruitment of factors involved in mRNA 3'-end formation and termination 2 , but how this sequence is initiated remains unclear. In a chemical-genetic screen, human protein phosphatase 1 (PP1) isoforms were identified as substrates of positive transcription elongation factor b (P-TEFb), also known as the cyclin-dependent kinase 9 (Cdk9)-cyclin T1 (CycT1) complex 3 . Here we show that Cdk9 and PP1 govern phosphorylation of the conserved elongation factor Spt5 in the fission yeast Schizosaccharomyces pombe. Cdk9 phosphorylates both Spt5 and a negative regulatory site on the PP1 isoform Dis2 4 . Sites targeted by Cdk9 in the Spt5 carboxy-terminal domain can be dephosphorylated by Dis2 in vitro, and dis2 mutations retard Spt5 dephosphorylation after inhibition of Cdk9 in vivo. Chromatin immunoprecipitation and sequencing analysis indicates that Spt5 is dephosphorylated as transcription complexes traverse the cleavage and polyadenylation signal, concomitant with the accumulation of Pol II phosphorylated at residue Ser2 of the carboxy-terminal domain consensus heptad repeat 5 . A conditionally lethal Dis2-inactivating mutation attenuates the drop in Spt5 phosphorylation on chromatin, promotes transcription beyond the normal termination zone (as detected by precision run-on transcription and sequencing 6 ) and is genetically suppressed by the ablation of Cdk9 target sites in Spt5. These results suggest that the transition of Pol II from elongation to termination coincides with a Dis2-dependent reversal of Cdk9 signalling-a switch that is analogous to a Cdk1-PP1 circuit that controls mitotic progression 4 .

Cdk9 regulates a promoter-proximal checkpoint to modulate RNA polymerase II elongation rate in fission yeast.

Post-translational modifications of the transcription elongation complex provide mechanisms to fine-tune gene expression, yet their specific impacts on RNA polymerase II regulation remain difficult to ascertain. Here, in Schizosaccharomyces pombe, we examine the role of Cdk9, and related Mcs6/Cdk7 and Lsk1/Cdk12 kinases, on transcription at base-pair resolution with Precision Run-On sequencing (PRO-seq). Within a minute of Cdk9 inhibition, phosphorylation of Pol II-associated factor, Spt5 is undetectable. The effects of Cdk9 inhibition are more severe than inhibition of Cdk7 and Cdk12, resulting in a shift of Pol II toward the transcription start site (TSS). A time course of Cdk9 inhibition reveals that early transcribing Pol II can escape promoter-proximal regions, but with a severely reduced elongation rate of only ~400 bp/min. Our results in fission yeast suggest the existence of a conserved global regulatory checkpoint that requires Cdk9 kinase activity.

Nascent RNA sequencing reveals a dynamic global transcriptional response at genes and enhancers to the natural medicinal compound celastrol.

Most studies of responses to transcriptional stimuli measure changes in cellular mRNA concentrations. By sequencing nascent RNA instead, it is possible to detect changes in transcription in minutes rather than hours and thereby distinguish primary from secondary responses to regulatory signals. Here, we describe the use of PRO-seq to characterize the immediate transcriptional response in human cells to celastrol, a compound derived from traditional Chinese medicine that has potent anti-inflammatory, tumor-inhibitory, and obesity-controlling effects. Celastrol is known to elicit a cellular stress response resembling the response to heat shock, but the transcriptional basis of this response remains unclear. Our analysis of PRO-seq data for K562 cells reveals dramatic transcriptional effects soon after celastrol treatment at a broad collection of both coding and noncoding transcription units. This transcriptional response occurred in two major waves, one within 10 min, and a second 40-60 min after treatment. Transcriptional activity was generally repressed by celastrol, but one distinct group of genes, enriched for roles in the heat shock response, displayed strong activation. Using a regression approach, we identified key transcription factors that appear to drive these transcriptional responses, including members of the E2F and RFX families. We also found sequence-based evidence that particular transcription factors drive the activation of enhancers. We observed increased polymerase pausing at both genes and enhancers, suggesting that pause release may be widely inhibited during the celastrol response. Our study demonstrates that a careful analysis of PRO-seq time-course data can disentangle key aspects of a complex transcriptional response, and it provides new insights into the activity of a powerful pharmacological agent.

Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq).

We provide a protocol for precision nuclear run-on sequencing (PRO-seq) and its variant, PRO-cap, which map the location of active RNA polymerases (PRO-seq) or transcription start sites (TSSs) (PRO-cap) genome-wide at high resolution. The density of RNA polymerases at a particular genomic locus directly reflects the level of nascent transcription at that region. Nuclei are isolated from cells and, under nuclear run-on conditions, transcriptionally engaged RNA polymerases incorporate one or, at most, a few biotin-labeled nucleotide triphosphates (biotin-NTPs) into the 3' end of nascent RNA. The biotin-labeled nascent RNA is used to prepare sequencing libraries, which are sequenced from the 3' end to provide high-resolution positional information for the RNA polymerases. PRO-seq provides much higher sensitivity than ChIP-seq, and it generates a much larger fraction of usable sequence reads than ChIP-seq or NET-seq (native elongating transcript sequencing). Similarly to NET-seq, PRO-seq maps the RNA polymerase at up to base-pair resolution with strand specificity, but unlike NET-seq it does not require immunoprecipitation. With the protocol provided here, PRO-seq (or PRO-cap) libraries for high-throughput sequencing can be generated in 4-5 working days. The method has been applied to human, mouse, Drosophila melanogaster and Caenorhabditis elegans cells and, with slight modifications, to yeast.

Divergence of a conserved elongation factor and transcription regulation in budding and fission yeast.

Complex regulation of gene expression in mammals has evolved from simpler eukaryotic systems, yet the mechanistic features of this evolution remain elusive. Here, we compared the transcriptional landscapes of the distantly related budding and fission yeast. We adapted the Precision Run-On sequencing (PRO-seq) approach to map the positions of RNA polymerase active sites genome-wide in Schizosaccharomyces pombe and Saccharomyces cerevisiae. Additionally, we mapped preferred sites of transcription initiation in each organism using PRO-cap. Unexpectedly, we identify a pause in early elongation, specific to S. pombe, that requires the conserved elongation factor subunit Spt4 and resembles promoter-proximal pausing in metazoans. PRO-seq profiles in strains lacking Spt4 reveal globally elevated levels of transcribing RNA Polymerase II (Pol II) within genes in both species. Messenger RNA abundance, however, does not reflect the increases in Pol II density, indicating a global reduction in elongation rate. Together, our results provide the first base-pair resolution map of transcription elongation in S. pombe and identify divergent roles for Spt4 in controlling elongation in budding and fission yeast.

PPP1CC2 can form a kinase/phosphatase complex with the testis-specific proteins TSSK1 and TSKS in the mouse testis.

The mouse protein phosphatase gene Ppp1cc is essential for male fertility, with mutants displaying a failure in spermatogenesis including a widespread loss of post-meiotic germ cells and abnormalities in the mitochondrial sheath. This phenotype is hypothesized to be responsible for the loss of the testis-specific isoform PPP1CC2. To identify PPP1CC2-interacting proteins with a function in spermatogenesis, we carried out GST pull-down assays in mouse testis lysates. Amongst the identified candidate interactors was the testis-specific protein kinase TSSK1, which is also essential for male fertility. Subsequent interaction experiments confirmed the capability of PPP1CC2 to form a complex with TSSK1 mediated by the direct interaction of each with the kinase substrate protein TSKS. Interaction between PPP1CC2 and TSKS is mediated through an RVxF docking motif on the TSKS surface. Phosphoproteomic analysis of the mouse testis identified a novel serine phosphorylation site within the TSKS RVxF motif that appears to negatively regulate binding to PPP1CC2. Immunohistochemical analysis of TSSK1 and TSKS in the Ppp1cc mutant testis showed reduced accumulation to distinct cytoplasmic foci and other abnormalities in their distribution consistent with the loss of germ cells and seminiferous tubule disorganization observed in the Ppp1cc mutant phenotype. A comparison of Ppp1cc and Tssk1/2 knockout phenotypes via electron microscopy revealed similar abnormalities in the morphology of the mitochondrial sheath. These data demonstrate a novel kinase/phosphatase complex in the testis that could play a critical role in the completion of spermatogenesis.