Role of cell-type specific nucleosome positioning in inducible activation of mammalian promoters

Activity-assembled nBAF complex mediates rapid immediate early gene transcription by regulating RNA polymerase II productive elongation

Signal-dependent RNA polymerase II (RNA Pol II) productive elongation is an integral component of gene transcription, including that of immediate early genes (IEGs) induced by neuronal activity. However, it remains unclear how productively elongating RNA Pol II overcomes nucleosomal barriers. Using RNAi, three degraders, and several small-molecule inhibitors, we show that the mammalian switch/sucrose non-fermentable (SWI/SNF) complex of neurons (neuronal BRG1/BRM-associated factor or nBAF) is required for activity-induced transcription of neuronal IEGs, including Arc. The nBAF complex facilitates promoter-proximal RNA Pol II pausing and signal-dependent RNA Pol II recruitment (loading) and, importantly, mediates productive elongation in the gene body via interaction with the elongation complex and elongation-competent RNA Pol II. Mechanistically, RNA Pol II elongation is mediated by activity-induced nBAF assembly (especially ARID1A recruitment) and its ATPase activity. Together, our data demonstrate that the nBAF complex regulates several aspects of RNA Pol II transcription and reveal mechanisms underlying activity-induced RNA Pol II elongation. These findings may offer insights into human maladies etiologically associated with mutational interdiction of BAF functions.

Evolutionary Innovations in Conserved Regulatory Elements Associate With Developmental Genes in Mammals

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Genetic variation within enhancer sequences is an important contributor to phenotypic variation including evolutionary adaptations and human disease. Certain genes and pathways may be more prone to regulatory evolution than others, with different patterns across diverse organisms, but whether such patterns exist has not been investigated at a sufficient scale. To address this question, we identified signatures of accelerated sequence evolution in conserved enhancer elements throughout the mammalian phylogeny at an unprecedented scale. While different genes and pathways were enriched for regulatory evolution in different parts of the tree, we found a striking overall pattern of pleiotropic genes involved in gene regulatory and developmental processes being enriched for accelerated enhancer evolution. These genes were connected to more enhancers than other genes, which was the basis for having an increased amount of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. Detailed study of one acceleration event in an enhancer of HES1 revealed that sequence evolution led to a new activity domain in the developing limb that emerged concurrently with the evolution of digit reduction in hoofed mammals. Our results provide evidence that enhancer evolution has been a frequent contributor to regulatory innovation at conserved developmental signaling genes in mammals.

Signalling Pathways of Inflammation and Cancer in Human Mononuclear Cells: Effect of Nanoparticle Air Pollutants

Fine inhalable particulate matter (PM) triggers an inflammatory response in the airways and activates mononuclear cells, mediators of tissue homeostasis, and tumour-promoting inflammation. We have assessed ex vivo responses of human monocytes and monocyte-derived macrophages to standardised air pollutants: carbon black, urban dust, and nanoparticulate carbon black, focusing on their pro-inflammatory and DNA-damaging properties. None of the PM (100 μg/mL/24 h) was significantly toxic to the cells, aside from inducing oxidative stress, fractional DNA damage, and inhibiting phagocytosis. TNFα was only slightly increased. PM nanoparticles increase the expression and activate DNA-damage-related histone H2A.X as well as pro-inflammatory NF-κB. We have shown that the urban dust stimulates the pathway of DNA damage/repair via the selective post-translational phosphorylation of H2A.X while nanoparticulate carbon black increases inflammation via activation of NF-κB. Moreover, the inflammatory response to lipopolysaccharide was significantly stronger in macrophages pre-exposed to urban dust or nanoparticulate carbon black. Our data show that airborne nanoparticles induce PM-specific, epigenetic alterations in the subsets of cultured mononuclear cells, which may be quantified using binary fluorescence scatterplots. Such changes intercede with inflammatory signalling and highlight important molecular and cell-specific epigenetic mechanisms of tumour-promoting inflammation.

Deep learning with a small dataset predicts chromatin remodelling contribution to winter dormancy of apple axillary buds

Epigenetic changes serve as a cellular memory for cumulative cold recognition in both herbaceous and tree species, including bud dormancy. However, most studies have discussed predicted chromatin structure with respect to histone marks. In the present study, we investigated the structural dynamics of bona fide chromatin to determine how plants recognize prolonged chilling during the initial stage of bud dormancy. The vegetative axillary buds of the 'Fuji' apple, which shows typical low temperature-dependent, but not photoperiod, dormancy induction, were used for the chromatin structure and transcriptional change analyses. The results were integrated using a deep-learning model and interpreted using statistical models, including Bayesian estimation. Although our model was constructed using a small dataset of two time points, chromatin remodelling due to random changes was excluded. The involvement of most nucleosome structural changes in transcriptional changes and the pivotal contribution of cold-driven circadian rhythm-dependent pathways regulated by the mobility of cis-regulatory elements were predicted. These findings may help to develop potential genetic targets for breeding species with less bud dormancy to overcome the effects of short winters during global warming. Our artificial intelligence concept can improve epigenetic analysis using a small dataset, especially in non-model plants with immature genome databases.

Detection and characterization of pancreatic and biliary tract cancers using cell-free DNA fragmentomics

BACKGROUND

Plasma cell-free DNA (cfDNA) fragmentomics has demonstrated significant differentiation power between cancer patients and healthy individuals, but little is known in pancreatic and biliary tract cancers. The aim of this study is to characterize the cfDNA fragmentomics in biliopancreatic cancers and develop an accurate method for cancer detection.

METHODS

One hundred forty-seven patients with biliopancreatic cancers and 71 non-cancer volunteers were enrolled, including 55 patients with cholangiocarcinoma, 30 with gallbladder cancer, and 62 with pancreatic cancer. Low-coverage whole-genome sequencing (median coverage: 2.9 ×) was performed on plasma cfDNA. Three cfDNA fragmentomic features, including fragment size, end motif and nucleosome footprint, were subjected to construct a stacked machine learning model for cancer detection. Integration of carbohydrate antigen 19-9 (CA19-9) was explored to improve model performance.

RESULTS

The stacked model presented robust performance for cancer detection (area under curve (AUC) of 0.978 in the training cohort, and AUC of 0.941 in the validation cohort), and remained consistent even when using extremely low-coverage sequencing depth of 0.5 × (AUC: 0.905). Besides, our method could also help differentiate biliopancreatic cancer subtypes. By integrating the stacked model and CA19-9 to generate the final detection model, a high accuracy in distinguishing biliopancreatic cancers from non-cancer samples with an AUC of 0.995 was achieved.

CONCLUSIONS

Our model demonstrated ultrasensitivity of plasma cfDNA fragementomics in detecting biliopancreatic cancers, fulfilling the unmet accuracy of widely-used serum biomarker CA19-9, and provided an affordable way for accurate noninvasive biliopancreatic cancer screening in clinical practice.

SWI/SNF-dependent genes are defined by their chromatin landscape

SWI/SNF complexes are evolutionarily conserved, ATP-dependent chromatin remodeling machines. Here, we characterize the features of SWI/SNF-dependent genes using BRM014, an inhibitor of the ATPase activity of the complexes. We find that SWI/SNF activity is required to maintain chromatin accessibility and nucleosome occupancy for most enhancers but not for most promoters. SWI/SNF activity is needed for expression of genes with low to medium levels of expression that have promoters with (1) low chromatin accessibility, (2) low levels of active histone marks, (3) high H3K4me1/H3K4me3 ratio, (4) low nucleosomal phasing, and (5) enrichment in TATA-box motifs. These promoters are mostly occupied by the canonical Brahma-related gene 1/Brahma-associated factor (BAF) complex. These genes are surrounded by SWI/SNF-dependent enhancers and mainly encode signal transduction, developmental, and cell identity genes (with almost no housekeeping genes). Machine-learning models trained with different chromatin characteristics of promoters and their surrounding regulatory regions indicate that the chromatin landscape is a determinant for establishing SWI/SNF dependency.

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

Activity-assembled nBAF complex mediates rapid immediate early gene transcription by regulating RNA Polymerase II productive elongation

Signal-dependent RNA Polymerase II (Pol2) productive elongation is an integral component of gene transcription, including those of immediate early genes (IEGs) induced by neuronal activity. However, it remains unclear how productively elongating Pol2 overcome nucleosomal barriers. Using RNAi, three degraders, and several small molecule inhibitors, we show that the mammalian SWI/SNF complex of neurons (neuronal BAF, or nBAF) is required for activity-induced transcription of neuronal IEGs, including Arc . The nBAF complex facilitates promoter-proximal Pol2 pausing, signal-dependent Pol2 recruitment (loading), and importantly, mediates productive elongation in the gene body via interaction with the elongation complex and elongation-competent Pol2. Mechanistically, Pol2 elongation is mediated by activity-induced nBAF assembly (especially, ARID1A recruitment) and its ATPase activity. Together, our data demonstrate that the nBAF complex regulates several aspects of Pol2 transcription and reveal mechanisms underlying activity-induced Pol2 elongation. These findings may offer insights into human maladies etiologically associated with mutational interdiction of BAF functions.

Comprehensive tissue deconvolution of cell-free DNA by deep learning for disease diagnosis and monitoring

Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

Nucleosome Patterns in Circulating Tumor DNA Reveal Transcriptional Regulation of Advanced Prostate Cancer Phenotypes

Advanced prostate cancers comprise distinct phenotypes, but tumor classification remains clinically challenging. Here, we harnessed circulating tumor DNA (ctDNA) to study tumor phenotypes by ascertaining nucleosome positioning patterns associated with transcription regulation. We sequenced plasma ctDNA whole genomes from patient-derived xenografts representing a spectrum of androgen receptor active (ARPC) and neuroendocrine (NEPC) prostate cancers. Nucleosome patterns associated with transcriptional activity were reflected in ctDNA at regions of genes, promoters, histone modifications, transcription factor binding, and accessible chromatin. We identified the activity of key phenotype-defining transcriptional regulators from ctDNA, including AR, ASCL1, HOXB13, HNF4G, and GATA2. To distinguish NEPC and ARPC in patient plasma samples, we developed prediction models that achieved accuracies of 97% for dominant phenotypes and 87% for mixed clinical phenotypes. Although phenotype classification is typically assessed by IHC or transcriptome profiling from tumor biopsies, we demonstrate that ctDNA provides comparable results with diagnostic advantages for precision oncology.

SIGNIFICANCE

This study provides insights into the dynamics of nucleosome positioning and gene regulation associated with cancer phenotypes that can be ascertained from ctDNA. New methods for classification in phenotype mixtures extend the utility of ctDNA beyond assessments of somatic DNA alterations with important implications for molecular classification and precision oncology. This article is highlighted in the In This Issue feature, p. 517.

Clustered PHD domains in KMT2/MLL proteins are attracted by H3K4me3 and H3 acetylation-rich active promoters and enhancers

Histone lysine-specific methyltransferase 2 (KMT2A-D) proteins, alternatively called mixed lineage leukemia (MLL1-4) proteins, mediate positive transcriptional memory. Acting as the catalytic subunits of human COMPASS-like complexes, KMT2A-D methylate H3K4 at promoters and enhancers. KMT2A-D contain understudied highly conserved triplets and a quartet of plant homeodomains (PHDs). Here, we show that all clustered (multiple) PHDs localize to the well-defined loci of H3K4me3 and H3 acetylation-rich active promoters and enhancers. Surprisingly, we observe little difference in binding pattern between PHDs from promoter-specific KMT2A-B and enhancer-specific KMT2C-D. Fusion of the KMT2A CXXC domain to the PHDs drastically enhances their preference for promoters over enhancers. Hence, the presence of CXXC domains in KMT2A-B, but not KMT2C-D, may explain the promoter/enhancer preferences of the full-length proteins. Importantly, targets of PHDs overlap with KMT2A targets and are enriched in genes involved in the cancer pathways. We also observe that PHDs of KMT2A-D are mutated in cancer, especially within conserved folding motifs (Cys4HisCys2Cys/His). The mutations cause a domain loss-of-function. Taken together, our data suggest that PHDs of KMT2A-D guide the full-length proteins to active promoters and enhancers, and thus play a role in positive transcriptional memory.

Epigenetic liquid biopsies for minimal residual disease, what's around the corner?

Liquid biopsy assays for minimal residual disease (MRD) are used to monitor and inform oncological treatment and predict the risk of relapse in cancer patients. To-date, most MRD assay development has focused on targeting somatic mutations. However, epigenetic changes are more frequent and universal than genetic alterations in cancer and circulating tumor DNA (ctDNA) retains much of these changes. Here, we review the epigenetic signals that can be used to detect MRD, including DNA methylation alterations and fragmentation patterns that differentiate ctDNA from noncancerous circulating cell-free DNA (ccfDNA). We then summarize the current state of MRD monitoring; highlight the advantages of epigenetics over genetics-based approaches; and discuss the emerging paradigm of assaying both genetic and epigenetic targets to monitor treatment response, detect disease recurrence, and inform adjuvant therapy.

Histone H1 binding to nucleosome arrays depends on linker DNA length and trajectory

Throughout the genome, nucleosomes often form regular arrays that differ in nucleosome repeat length (NRL), occupancy of linker histone H1 and transcriptional activity. Here, we report cryo-EM structures of human H1-containing tetranucleosome arrays with four physiologically relevant NRLs. The structures show a zig-zag arrangement of nucleosomes, with nucleosomes 1 and 3 forming a stack. H1 binding to stacked nucleosomes depends on the NRL, whereas H1 always binds to the non-stacked nucleosomes 2 and 4. Short NRLs lead to altered trajectories of linker DNA, and these altered trajectories sterically impair H1 binding to the stacked nucleosomes in our structures. As the NRL increases, linker DNA trajectories relax, enabling H1 contacts and binding. Our results provide an explanation for why arrays with short NRLs are depleted of H1 and suited for transcription, whereas arrays with long NRLs show full H1 occupancy and can form transcriptionally silent heterochromatin regions.

Cryo-EM structures of human H1-containing tetranucleosome arrays with distinct, physiological nucleosome repeat lengths reveal that nucleosomes assume a zig-zag arrangement and H1 binds to stacked nucleosomes with longer linker DNA.

Nucleosome footprinting in plasma cell-free DNA for the pre-surgical diagnosis of ovarian cancer

Fragmentation patterns of plasma cell-free DNA (cfDNA) are known to reflect nucleosome positions of cell types contributing to cfDNA. Based on cfDNA fragmentation patterns, the deviation in nucleosome footprints was quantified between diagnosed ovarian cancer patients and healthy individuals. Multinomial modeling was subsequently applied to capture these deviations in a per sample nucleosome footprint score. Validation was performed in 271 cfDNAs pre-surgically collected from women with an adnexal mass. We confirmed that nucleosome scores were elevated in invasive carcinoma patients, but not in patients with benign or borderline disease. Combining nucleosome scores with chromosomal instability scores assessed in the same cfDNA improved prediction of malignancy. Nucleosome scores were, however, more reliable to predict non-high-grade serous ovarian tumors, which are characterized by low chromosomal instability. These data highlight that compared to chromosomal instability, nucleosome footprinting provides a complementary and more generic read-out for pre-surgical diagnosis of invasive disease in women with adnexal masses.

Differences in RNA polymerase II complexes and their interactions with surrounding chromatin on human and cytomegalovirus genomes

Interactions of the RNA polymerase II (Pol II) preinitiation complex (PIC) and paused early elongation complexes with the first downstream (+1) nucleosome are thought to be functionally important. However, current methods are limited for investigating these relationships, both for cellular chromatin and the human cytomegalovirus (HCMV) genome. Digestion with human DNA fragmentation factor (DFF) before immunoprecipitation (DFF-ChIP) precisely revealed both similarities and major differences in PICs driven by TBP on the host genome in comparison with PICs driven by TBP or the viral-specific, late initiation factor UL87 on the viral genome. Host PICs and paused Pol II complexes are frequently found in contact with the +1 nucleosome and paused Pol II can also be found in a complex involved in the initial invasion of the +1 nucleosome. In contrast, viral transcription complexes have very limited nucleosomal interactions, reflecting a relative lack of chromatinization of transcriptionally active regions of HCMV genomes.

Here the authors digested chromatin with DNA fragmentation factor (DFF) prior to chromatin immunoprecipitation (DFF-ChIP) to depict transcription complex interactions with neighboring nucleosomes in cells. Applying this method to human cytomegalovirus (HMCV)-infected cells, they find that the viral genome is underchromatinized, leading to fewer transcription complex interactions with nucleosomes.

The MuvB complex binds and stabilizes nucleosomes downstream of the transcription start site of cell-cycle dependent genes

The chromatin architecture in promoters is thought to regulate gene expression, but it remains uncertain how most transcription factors (TFs) impact nucleosome position. The MuvB TF complex regulates cell-cycle dependent gene-expression and is critical for differentiation and proliferation during development and cancer. MuvB can both positively and negatively regulate expression, but the structure of MuvB and its biochemical function are poorly understood. Here we determine the overall architecture of MuvB assembly and the crystal structure of a subcomplex critical for MuvB function in gene repression. We find that the MuvB subunits LIN9 and LIN37 function as scaffolding proteins that arrange the other subunits LIN52, LIN54 and RBAP48 for TF, DNA, and histone binding, respectively. Biochemical and structural data demonstrate that MuvB binds nucleosomes through an interface that is distinct from LIN54-DNA consensus site recognition and that MuvB increases nucleosome occupancy in a reconstituted promoter. We find in arrested cells that MuvB primarily associates with a tightly positioned +1 nucleosome near the transcription start site (TSS) of MuvB-regulated genes. These results support a model that MuvB binds and stabilizes nucleosomes just downstream of the TSS on its target promoters to repress gene expression.

The MuvB protein complex regulates genes that are differentially expressed through the cell cycle, yet its precise molecular function has remained unclear. Here the authors reveal MuvB associates with the nucleosome adjacent to the transcription start site of cell-cycle genes and that the tight positioning of this nucleosome correlates with MuvB-dependent gene repression.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.

A Comprehensive Toolbox to Analyze Enhancer-Promoter Functions

Knowledge in gene transcription and chromatin regulation has been intensely studied for decades, but thanks to next-generation sequencing (NGS) techniques there has been a major leap forward in the last few years. Historically, identification of specific enhancer elements has led to the identification of master transcription factors (TFs) in the 1990s. Genetic and biochemical experiments have identified the key regulators controlling RNA polymerase II (RNAPII) transcription and structurally analyses have elucidated detailed mechanisms. NGS and the development of chromatin immunoprecipitation (ChIP) have accelerated the gain of knowledge in the recent years. By now, we have a dazzling wealth of techniques that are currently used to put gene expression into a genome-wide context. This book is an attempt to assemble useful protocols for many researchers within and nearby research areas. In general, these innovative techniques focus on enhancer and promoter studies. The techniques should also be of interest for related fields such as DNA repair and replication.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

MNase Profiling of Promoter Chromatin in Salmonella typhimurium-Stimulated GM12878 Cells Reveals Dynamic and Response-Specific Nucleosome Architecture

The nucleosome is the primary unit of chromatin structure and commonly imputed as a regulator of nuclear events, although the exact mechanisms remain unclear. Recent studies have shown that certain nucleosomes can have different sensitivities to micrococcal nuclease (MNase) digestion, resulting in the release of populations of nucleosomes dependent on the concentration of MNase. Mapping MNase sensitivity of nucleosomes at transcription start sites genome-wide reveals an important functional nucleosome organization that correlates with gene expression levels and transcription factor binding. In order to understand nucleosome distribution and sensitivity dynamics during a robust genome response, we mapped nucleosome position and sensitivity using multiple concentrations of MNase. We used the innate immune response as a model system to understand chromatin-mediated regulation. Herein we demonstrate that stimulation of a human lymphoblastoid cell line (GM12878) with heat-killed Salmonella typhimurium (HKST) results in changes in nucleosome sensitivity to MNase. We show that the HKST response alters the sensitivity of -1 nucleosomes at highly expressed promoters. Finally, we correlate the increased sensitivity with response-specific transcription factor binding. These results indicate that nucleosome sensitivity dynamics reflect the cellular response to HKST and pave the way for further studies that will deepen our understanding of the specificity of genome response.