Detecting Differential Transcription Factor Activity from ATAC-Seq Data

ONECUT2 acts as a lineage plasticity driver in adenocarcinoma as well as neuroendocrine variants of prostate cancer

Androgen receptor- (AR-) indifference is a mechanism of resistance to hormonal therapy in prostate cancer (PC). Here we demonstrate that ONECUT2 (OC2) activates resistance through multiple drivers associated with adenocarcinoma, stem-like and neuroendocrine (NE) variants. Direct OC2 gene targets include the glucocorticoid receptor (GR; NR3C1) and the NE splicing factor SRRM4, which are key drivers of lineage plasticity. Thus, OC2, despite its previously described NEPC driver function, can indirectly activate a portion of the AR cistrome through epigenetic activation of GR. Mechanisms by which OC2 regulates gene expression include promoter binding, enhancement of genome-wide chromatin accessibility, and super-enhancer reprogramming. Pharmacologic inhibition of OC2 suppresses lineage plasticity reprogramming induced by the AR signaling inhibitor enzalutamide. These results demonstrate that OC2 activation promotes a range of drug resistance mechanisms associated with treatment-emergent lineage variation in PC and support enhanced efforts to therapeutically target OC2 as a means of suppressing treatment-resistant disease.

Viral reprogramming of host transcription initiation

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Computational tools for inferring transcription factor activity

Transcription factors (TFs) are essential players in orchestrating the regulatory landscape in cells. Still, their exact modes of action and dependencies on other regulatory aspects remain elusive. Since TFs act cell type-specific and each TF has its own characteristics, untangling their regulatory interactions from an experimental point of view is laborious and convoluted. Thus, there is an ongoing development of computational tools that estimate transcription factor activity (TFA) from a variety of data modalities, either based on a mapping of TFs to their putative target genes or in a genome-wide, gene-unspecific fashion. These tools can help to gain insights into TF regulation and to prioritize candidates for experimental validation. We want to give an overview of available computational tools that estimate TFA, illustrate examples of their application, debate common result validation strategies, and discuss assumptions and concomitant limitations.

ONECUT2 Activates Diverse Resistance Drivers of Androgen Receptor-Independent Heterogeneity in Prostate Cancer

Androgen receptor- (AR-) indifference is a mechanism of resistance to hormonal therapy in prostate cancer (PC). Here we demonstrate that the HOX/CUT transcription factor ONECUT2 (OC2) activates resistance through multiple drivers associated with adenocarcinoma, stem-like and neuroendocrine (NE) variants. Direct OC2 targets include the glucocorticoid receptor and the NE splicing factor SRRM4, among others. OC2 regulates gene expression by promoter binding, enhancement of chromatin accessibility, and formation of novel super-enhancers. OC2 also activates glucuronidation genes that irreversibly disable androgen, thereby evoking phenotypic heterogeneity indirectly by hormone depletion. Pharmacologic inhibition of OC2 suppresses lineage plasticity reprogramming induced by the AR signaling inhibitor enzalutamide. These results demonstrate that OC2 activation promotes a range of drug resistance mechanisms associated with treatment-emergent lineage variation in PC. Our findings support enhanced efforts to therapeutically target this protein as a means of suppressing treatment-resistant disease.

Loss of SYNCRIP unleashes APOBEC-driven mutagenesis, tumor heterogeneity, and AR-targeted therapy resistance in prostate cancer

Tumor mutational burden and heterogeneity has been suggested to fuel resistance to many targeted therapies. The cytosine deaminase APOBEC proteins have been implicated in the mutational signatures of more than 70% of human cancers. However, the mechanism underlying how cancer cells hijack the APOBEC mediated mutagenesis machinery to promote tumor heterogeneity, and thereby foster therapy resistance remains unclear. We identify SYNCRIP as an endogenous molecular brake which suppresses APOBEC-driven mutagenesis in prostate cancer (PCa). Overactivated APOBEC3B, in SYNCRIP-deficient PCa cells, is a key mutator, representing the molecular source of driver mutations in some frequently mutated genes in PCa, including FOXA1, EP300. Functional screening identifies eight crucial drivers for androgen receptor (AR)-targeted therapy resistance in PCa that are mutated by APOBEC3B: BRD7, CBX8, EP300, FOXA1, HDAC5, HSF4, STAT3, and AR. These results uncover a cell-intrinsic mechanism that unleashes APOBEC-driven mutagenesis, which plays a significant role in conferring AR-targeted therapy resistance in PCa.

The time-resolved genomic impact of Wnt/β-catenin signaling

Wnt signaling orchestrates gene expression via its effector, β-catenin. However, it is unknown whether β-catenin binds its target genomic regions simultaneously and how this impacts chromatin dynamics to modulate cell behavior. Using a combination of time-resolved CUT&RUN against β-catenin, ATAC-seq, and perturbation assays in different cell types, we show that Wnt/β-catenin physical targets are tissue-specific, β-catenin "moves" on different loci over time, and its association to DNA accompanies changing chromatin accessibility landscapes that determine cell behavior. In particular, Wnt/β-catenin progressively shapes the chromatin of human embryonic stem cells (hESCs) as they undergo mesodermal differentiation, a behavior that we define as "plastic." In HEK293T cells, on the other hand, Wnt/β-catenin drives a transient chromatin opening, followed by re-establishment of the pre-stimulation state, a response that we define as "elastic." Future experiments shall assess whether other cell communication mechanisms, in addition to Wnt signaling, are ruled by time, cellular idiosyncrasies, and chromatin constraints. A record of this paper's transparent peer review process is included in the supplemental information.

Chi3l1 Is a Modulator of Glioma Stem Cell States and a Therapeutic Target in Glioblastoma

Chitinase 3-like 1 (Chi3l1) is a secreted protein that is highly expressed in glioblastoma. Here, we show that Chi3l1 alters the state of glioma stem cells (GSC) to support tumor growth. Exposure of patient-derived GSCs to Chi3l1 reduced the frequency of CD133+SOX2+ cells and increased the CD44+Chi3l1+ cells. Chi3l1 bound to CD44 and induced phosphorylation and nuclear translocation of β-catenin, Akt, and STAT3. Single-cell RNA sequencing and RNA velocity following incubation of GSCs with Chi3l1 showed significant changes in GSC state dynamics driving GSCs towards a mesenchymal expression profile and reducing transition probabilities towards terminal cellular states. ATAC-seq revealed that Chi3l1 increases accessibility of promoters containing a Myc-associated zinc finger protein (MAZ) transcription factor footprint. Inhibition of MAZ downregulated a set of genes with high expression in cellular clusters that exhibit significant cell state transitions after treatment with Chi3l1, and MAZ deficiency rescued the Chi3L-induced increase of GSC self-renewal. Finally, targeting Chi3l1 in vivo with a blocking antibody inhibited tumor growth and increased the probability of survival. Overall, this work suggests that Chi3l1 interacts with CD44 on the surface of GSCs to induce Akt/β-catenin signaling and MAZ transcriptional activity, which in turn upregulates CD44 expression in a pro-mesenchymal feed-forward loop. The role of Chi3l1 in regulating cellular plasticity confers a targetable vulnerability to glioblastoma.

SIGNIFICANCE

Chi3l1 is a modulator of glioma stem cell states that can be targeted to promote differentiation and suppress growth of glioblastoma.

Microscale combinatorial stimulation of human myeloid cells reveals inflammatory priming by viral ligands

Cells sense a wide variety of signals and respond by adopting complex transcriptional states. Most single-cell profiling is carried out today at cellular baseline, blind to cells' potential spectrum of functional responses. Exploring the space of cellular responses experimentally requires access to a large combinatorial perturbation space. Single-cell genomics coupled with multiplexing techniques provide a useful tool for characterizing cell states across several experimental conditions. However, current multiplexing strategies require programmatic handling of many samples in macroscale arrayed formats, precluding their application in large-scale combinatorial analysis. Here, we introduce StimDrop, a method that combines antibody-based cell barcoding with parallel droplet processing to automatically formulate cell population × stimulus combinations in a microfluidic device. We applied StimDrop to profile the effects of 512 sequential stimulation conditions on human dendritic cells. Our results demonstrate that priming with viral ligands potentiates hyperinflammatory responses to a second stimulus, and show transcriptional signatures consistent with this phenomenon in myeloid cells of patients with severe COVID-19.

Targeting KDM4 for treating PAX3-FOXO1-driven alveolar rhabdomyosarcoma

Chimeric transcription factors drive lineage-specific oncogenesis but are notoriously difficult to target. Alveolar rhabdomyosarcoma (RMS) is an aggressive childhood soft tissue sarcoma transformed by the pathognomonic Paired Box 3-Forkhead Box O1 (PAX3-FOXO1) fusion protein, which governs a core regulatory circuitry transcription factor network. Here, we show that the histone lysine demethylase 4B (KDM4B) is a therapeutic vulnerability for PAX3-FOXO1⁺ RMS. Genetic and pharmacologic inhibition of KDM4B substantially delayed tumor growth. Suppression of KDM4 proteins inhibited the expression of core oncogenic transcription factors and caused epigenetic alterations of PAX3-FOXO1-governed superenhancers. Combining KDM4 inhibition with cytotoxic chemotherapy led to tumor regression in preclinical PAX3-FOXO1⁺ RMS subcutaneous xenograft models. In summary, we identified a targetable mechanism required for maintenance of the PAX3-FOXO1-related transcription factor network, which may translate to a therapeutic approach for fusion-positive RMS.

Epigenetic Regulation of Endothelial Cell Lineages During Zebrafish Development-New Insights From Technical Advances

Epigenetic regulation is integral in orchestrating the spatiotemporal regulation of gene expression which underlies tissue development. The emergence of new tools to assess genome-wide epigenetic modifications has enabled significant advances in the field of vascular biology in zebrafish. Zebrafish represents a powerful model to investigate the activity of cis-regulatory elements in vivo by combining technologies such as ATAC-seq, ChIP-seq and CUT&Tag with the generation of transgenic lines and live imaging to validate the activity of these regulatory elements. Recently, this approach led to the identification and characterization of key enhancers of important vascular genes, such as gata2a, notch1b and dll4. In this review we will discuss how the latest technologies in epigenetics are being used in the zebrafish to determine chromatin states and assess the function of the cis-regulatory sequences that shape the zebrafish vascular network.

Macrophage programming is regulated by a cooperative interaction between fatty acid binding protein 5 and peroxisome proliferator-activated receptor γ

Resolution of inflammation is an active process that is tightly regulated to achieve repair and tissue homeostasis. In the absence of resolution, persistent inflammation underlies the pathogenesis of chronic lung disease such as chronic obstructive pulmonary disease (COPD) with recurrent exacerbations. Over the course of inflammation, macrophage programming transitions from pro-inflammatory to pro-resolving, which is in part regulated by the nuclear receptor Peroxisome Proliferator-Activated Receptor γ (PPARγ). Our previous work demonstrated an association between Fatty Acid Binding Protein 5 (FABP5) expression and PPARγ activity in peripheral blood mononuclear cells of healthy and COPD patients. However, a role for FABP5 in macrophage programming has not been examined. Here, using a combination of in vitro and in vivo approaches, we demonstrate that FABP5 is necessary for PPARγ activation. In turn, PPARγ acts directly to increase FABP5 expression in primary human alveolar macrophages. We further illustrate that lack of FABP5 expression promotes a pro-inflammatory macrophage programming with increased secretion of pro-inflammatory cytokines and increased chromatin accessibility for pro-inflammatory transcription factors (e.g., NF-κB and MAPK). And finally, real-time cell metabolic analysis using the Seahorse technology shows an inhibition of oxidative phosphorylation in FABP5-deficient macrophages. Taken together, our data indicate that FABP5 and PPARγ reciprocally regulate each other's expression and function, consistent with a novel positive feedback loop between the two factors that mediates macrophage pro-resolving programming. Our studies highlight the importance of defining targets and regulatory mechanisms that control the resolution of inflammation and may serve to inform novel interventional strategies directed towards COPD.

UBR7 acts as a histone chaperone for post-nucleosomal histone H3

Histone chaperones modulate the stability of histones beginning from histone synthesis, through incorporation into DNA, and during recycling during transcription and replication. Following histone removal from DNA, chaperones regulate histone storage and degradation. Here, we demonstrate that UBR7 is a histone H3.1 chaperone that modulates the supply of pre-existing post-nucleosomal histone complexes. We demonstrate that UBR7 binds to post-nucleosomal H3K4me3 and H3K9me3 histones via its UBR box and PHD. UBR7 binds to the non-nucleosomal histone chaperone NASP. In the absence of UBR7, the pool of NASP-bound post-nucleosomal histones accumulate and chromatin is depleted of H3K4me3-modified histones. We propose that the interaction of UBR7 with NASP and histones opposes the histone storage functions of NASP and that UBR7 promotes reincorporation of post-nucleosomal H3 complexes.

Deconvolution of multiplexed transcriptional responses to wood smoke particles defines rapid aryl hydrocarbon receptor signaling dynamics

The heterogeneity of respirable particulates and compounds complicates our understanding of transcriptional responses to air pollution. Here, we address this by applying precision nuclear run-on sequencing and the assay for transposase-accessible chromatin sequencing to measure nascent transcription and chromatin accessibility in airway epithelial cells after wood smoke particle (WSP) exposure. We used transcription factor enrichment analysis to identify temporally distinct roles for ternary response factor-serum response factor complexes, the aryl hydrocarbon receptor (AHR), and NFκB in regulating transcriptional changes induced by WSP. Transcription of canonical targets of the AHR, such as CYP1A1 and AHRR, was robustly increased after just 30 min of WSP exposure, and we discovered novel AHR-regulated pathways and targets including the DNA methyltransferase, DNMT3L. Transcription of these genes and associated enhancers rapidly returned to near baseline by 120 min after exposure. The kinetics of AHR- and NFκB-regulated responses to WSP were distinguishable based on the timing of both transcriptional responses and chromatin remodeling, with induction of several cytokines implicated in maintaining NFκB-mediated responses through 120 min of exposure. In aggregate, our data establish a direct and primary role for AHR in mediating airway epithelial responses to WSP and identify crosstalk between AHR and NFκB signaling in controlling proinflammatory gene expression. This work also defines an integrated genomics-based strategy for deconvoluting multiplexed transcriptional responses to heterogeneous environmental exposures.

Chromatin remodeling due to degradation of citrate carrier impairs osteogenesis of aged mesenchymal stem cells

Aging is accompanied by a general decline in the function of many cellular pathways. However, whether these are causally or functionally interconnected remains elusive. Here, we study the effect of mitochondrial-nuclear communication on stem cell aging. We show that aged mesenchymal stem cells exhibit reduced chromatin accessibility and lower histone acetylation, particularly on promoters and enhancers of osteogenic genes. The reduced histone acetylation is due to impaired export of mitochondrial acetyl-CoA, owing to the lower levels of citrate carrier (CiC). We demonstrate that aged cells showed enhanced lysosomal degradation of CiC, which is mediated via mitochondrial-derived vesicles. Strikingly, restoring cytosolic acetyl-CoA levels either by exogenous CiC expression or via acetate supplementation, remodels the chromatin landscape and rescues the osteogenesis defects of aged mesenchymal stem cells. Collectively, our results establish a tight, age-dependent connection between mitochondrial quality control, chromatin and stem cell fate, which are linked together by CiC.

Uncovering Divergence in Gene Expression Regulation in the Adaptation of Yeast to Nitrogen Scarcity

Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.

Inferring TF activities and activity regulators from gene expression data with constraints from TF perturbation data

Motivation

The activity of a transcription factor (TF) in a sample of cells is the extent to which it is exerting its regulatory potential. Many methods of inferring TF activity from gene expression data have been described, but due to the lack of appropriate large-scale datasets, systematic and objective validation has not been possible until now.

Results

We systematically evaluate and optimize the approach to TF activity inference in which a gene expression matrix is factored into a condition-independent matrix of control strengths and a condition-dependent matrix of TF activity levels. We find that expression data in which the activities of individual TFs have been perturbed are both necessary and sufficient for obtaining good performance. To a considerable extent, control strengths inferred using expression data from one growth condition carry over to other conditions, so the control strength matrices derived here can be used by others. Finally, we apply these methods to gain insight into the upstream factors that regulate the activities of yeast TFs Gcr2, Gln3, Gcn4 and Msn2.

Availability and implementation

Evaluation code and data are available at https://doi.org/10.5281/zenodo.4050573.

Supplementary information

Supplementary data are available at Bioinformatics online.

Prenet: Predictive network from ATAC-SEQ data

Assays for transposase-accessible chromatin sequencing (ATAC-seq) provides an innovative approach to study chromatin status in multiple cell types. Moreover, it is also possible to efficiently extract differentially accessible chromatin (DACs) regions by using state-of-the-art algorithms (e.g. DESeq2) to predict gene activity in specific samples. Furthermore, it has recently been shown that small dips in sequencing peaks can be attributed to the binding of transcription factors. These dips, also known as footprints, can be used to identify trans-regulating interactions leading to gene expression. Current protocols used to identify footprints (e.g. pyDNAse and HINT-ATAC) have shown limitations resulting in the discovery of many false positive footprints. We generated a novel approach to identify genuine footprints within any given ATAC-seq dataset. Herein, we developed a new pipeline embedding DACs together with bona fide footprints resulting in the generation of a Predictive gene regulatory Network (PreNet) simply from ATAC-seq data. We further demonstrated that PreNet can be used to unveil meaningful molecular regulatory pathways in a given cell type.

Metabolic effects of air pollution exposure and reversibility

Air pollution involving particulate matter smaller than 2.5 μm in size (PM2.5) is the world's leading environmental risk factor contributing to mortality through cardiometabolic pathways. In this study, we modeled early life exposure using chow-fed C57BL/6J male mice that were exposed to real-world inhaled, concentrated PM2.5 (~10 times ambient levels/~60-120 μg/m3) or filtered air over a 14-week period. We investigated the effects of PM2.5 on phenotype, the transcriptome, and chromatin accessibility and compared these with the effects of a prototypical high-fat diet (HFD) as well as cessation of exposure on phenotype reversibility. Exposure to PM2.5 impaired glucose and insulin tolerance and reduced energy expenditure and 18FDG-PET uptake in brown adipose tissue. Multiple differentially expressed gene clusters in pathways involving metabolism and circadian rhythm were noted in insulin-responsive tissues. Although the magnitude of transcriptional change detected with PM2.5 exposure was lower than that observed with a HFD, the degree of alteration in chromatin accessibility after PM2.5 exposure was significant. The novel chromatin remodeler SMARCA5 (SWI/SNF complex) was regulated in response to PM2.5 exposure, the cessation of which was associated with a reversal of insulin resistance and restoration of chromatin accessibility and nucleosome positioning near transcription start sites, as well as a reversal of exposure-induced changes in the transcriptome, including SMARCA5. These changes indicate pliable epigenetic control mechanisms following cessation of exposure.

The MUC5B-associated variant rs35705950 resides within an enhancer subject to lineage- and disease-dependent epigenetic remodeling

The G/T transversion rs35705950, located approximately 3 kb upstream of the MUC5B start site, is the cardinal risk factor for idiopathic pulmonary fibrosis (IPF). Here, we investigate the function and chromatin structure of this –3 kb region and provide evidence that it functions as a classically defined enhancer subject to epigenetic programming. We use nascent transcript analysis to show that RNA polymerase II loads within 10 bp of the G/T transversion site, definitively establishing enhancer function for the region. By integrating Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) analysis of fresh and cultured human airway epithelial cells with nuclease sensitivity data, we demonstrate that this region is in accessible chromatin that affects the expression of MUC5B. Through applying paired single-nucleus RNA- and ATAC-seq to frozen tissue from IPF lungs, we extend these findings directly to disease, with results indicating that epigenetic programming of the –3 kb enhancer in IPF occurs in both MUC5B-expressing and nonexpressing lineages. In aggregate, our results indicate that the MUC5B-associated variant rs35705950 resides within an enhancer that is subject to epigenetic remodeling and contributes to pathologic misexpression in IPF.

Chromatin accessibility in canine stromal cells and its implications for canine somatic cell reprogramming

Naturally occurring disease in pet dogs is an untapped and unique resource for stem cell-based regenerative medicine translational research, given the many similarities and complexity such disease shares with their human counterparts. Canine-specific regulators of somatic cell reprogramming and pluripotency maintenance are poorly understood. While retroviral delivery of the four Yamanaka factors successfully reprogrammed canine embryonic fibroblasts, adult stromal cells remained resistant to reprogramming in spite of effective viral transduction and transgene expression. We hypothesized that adult stromal cells fail to reprogram due to an epigenetic barrier. Here, we performed assay for transposase-accessible chromatin using sequencing (ATAC-seq) on canine stromal and pluripotent stem cells, analyzing 51 samples in total, and establishing the global landscape of chromatin accessibility before and after reprogramming to induced pluripotent stem cells (iPSC). We also studied adult stromal cells that do not yield iPSC colonies to identify potential reprogramming barriers. ATAC-seq analysis identified distinct cell type clustering patterns and chromatin remodeling during embryonic fibroblast reprogramming. Compared with embryonic fibroblasts, adult stromal cells had a chromatin accessibility landscape that reflects phenotypic differentiation and somatic cell-fate stability. We ultimately identified 76 candidate genes and several transcription factor binding motifs that may be impeding somatic cell reprogramming to iPSC, and could be targeted for inhibition or activation, in order to improve the process in canines. These results provide a vast resource for better understanding of pluripotency regulators in dogs and provide an unbiased rationale for novel canine-specific reprogramming approaches.

Integrative transcriptome and chromatin landscape analysis reveals distinct epigenetic regulations in human memory B cells

Memory B cells (MBCs) are long-lived and produce high-affinity, generally, class-switched antibodies. Here, we use a multiparameter approach involving CD27 to segregate naïve B cells (NBC), IgD⁺ unswitched (unsw)MBCs and IgG⁺ or IgA⁺ class-switched (sw)MBCs from humans of different age, sex and race. Conserved antibody variable gene expression indicates that MBCs emerge through unbiased selection from NBCs. Integrative analyses of mRNAs, miRNAs, lncRNAs, chromatin accessibility and cis-regulatory elements uncover a core mRNA-ncRNA transcriptional signature shared by IgG⁺ and IgA⁺ swMBCs and distinct from NBCs, while unswMBCs display a transitional transcriptome. Some swMBC transcriptional signature loci are accessible but not expressed in NBCs. Profiling miRNAs reveals downregulated MIR181, and concomitantly upregulated MIR181 target genes such as RASSF6, TOX, TRERF1, TRPV3 and RORα, in swMBCs. Finally, lncRNAs differentially expressed in swMBCs cluster proximal to the IgH chain locus on chromosome 14. Our findings thus provide new insights into MBC transcriptional programs and epigenetic regulation, opening new investigative avenues on these critical cell elements in human health and disease.

Human memory B cells differentiate from naïve B cells and can express different immunoglobulin (Ig) isotypes resulted from class-switch recombination. Here the authors describe, using transcriptional and epigenetic data from human memory B cells and integrated multi-omics analyses, the differentiation regulation and trajectory of IgG⁺, IgA⁺ and IgD⁺ memory B cells.

ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation

While footprinting analysis of ATAC-seq data can theoretically enable investigation of transcription factor (TF) binding, the lack of a computational tool able to conduct different levels of footprinting analysis has so-far hindered the widespread application of this method. Here we present TOBIAS, a comprehensive, accurate, and fast footprinting framework enabling genome-wide investigation of TF binding dynamics for hundreds of TFs simultaneously. We validate TOBIAS using paired ATAC-seq and ChIP-seq data, and find that TOBIAS outperforms existing methods for bias correction and footprinting. As a proof-of-concept, we illustrate how TOBIAS can unveil complex TF dynamics during zygotic genome activation in both humans and mice, and propose how zygotic Dux activates cascades of TFs, binds to repeat elements and induces expression of novel genetic elements.

Analytical Approaches for ATAC-seq Data Analysis

ATAC-seq, the assay for transposase-accessible chromatin using sequencing, is a quick and efficient approach to investigating the chromatin accessibility landscape. Investigating chromatin accessibility has broad utility for answering many biological questions, such as mapping nucleosomes, identifying transcription factor binding sites, and measuring differential activity of DNA regulatory elements. Because the ATAC-seq protocol is both simple and relatively inexpensive, there has been a rapid increase in the availability of chromatin accessibility data. Furthermore, advances in ATAC-seq protocols are rapidly extending its breadth to additional experimental conditions, cell types, and species. Accompanying the increase in data, there has also been an explosion of new tools and analytical approaches for analyzing it. Here, we explain the fundamentals of ATAC-seq data processing, summarize common analysis approaches, and review computational tools to provide recommendations for different research questions. This primer provides a starting point and a reference for analysis of ATAC-seq data. © 2020 Wiley Periodicals LLC.

Loss of CHD1 Promotes Heterogeneous Mechanisms of Resistance to AR-Targeted Therapy via Chromatin Dysregulation

Metastatic prostate cancer is characterized by recurrent genomic copy number alterations that are presumed to contribute to resistance to hormone therapy. We identified CHD1 loss as a cause of antiandrogen resistance in an in vivo small hairpin RNA (shRNA) screen of 730 genes deleted in prostate cancer. ATAC-seq and RNA-seq analyses showed that CHD1 loss resulted in global changes in open and closed chromatin with associated transcriptomic changes. Integrative analysis of this data, together with CRISPR-based functional screening, identified four transcription factors (NR3C1, POU3F2, NR2F1, and TBX2) that contribute to antiandrogen resistance, with associated activation of non-luminal lineage programs. Thus, CHD1 loss results in chromatin dysregulation, thereby establishing a state of transcriptional plasticity that enables the emergence of antiandrogen resistance through heterogeneous mechanisms.

TFIID Enables RNA Polymerase II Promoter-Proximal Pausing

RNA polymerase II (RNAPII) transcription is governed by the pre-initiation complex (PIC), which contains TFIIA, TFIIB, TFIID, TFIIE, TFIIF, TFIIH, RNAPII, and Mediator. After initiation, RNAPII enzymes pause after transcribing less than 100 bases; precisely how RNAPII pausing is enforced and regulated remains unclear. To address specific mechanistic questions, we reconstituted human RNAPII promoter-proximal pausing in vitro, entirely with purified factors (no extracts). As expected, NELF and DSIF increased pausing, and P-TEFb promoted pause release. Unexpectedly, the PIC alone was sufficient to reconstitute pausing, suggesting RNAPII pausing is an inherent PIC function. In agreement, pausing was lost upon replacement of the TFIID complex with TATA-binding protein (TBP), and PRO-seq experiments revealed widespread disruption of RNAPII pausing upon acute depletion (t = 60 min) of TFIID subunits in human or Drosophila cells. These results establish a TFIID requirement for RNAPII pausing and suggest pause regulatory factors may function directly or indirectly through TFIID.

From reads to insight: a hitchhiker's guide to ATAC-seq data analysis

Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) is widely used in studying chromatin biology, but a comprehensive review of the analysis tools has not been completed yet. Here, we discuss the major steps in ATAC-seq data analysis, including pre-analysis (quality check and alignment), core analysis (peak calling), and advanced analysis (peak differential analysis and annotation, motif enrichment, footprinting, and nucleosome position analysis). We also review the reconstruction of transcriptional regulatory networks with multiomics data and highlight the current challenges of each step. Finally, we describe the potential of single-cell ATAC-seq and highlight the necessity of developing ATAC-seq specific analysis tools to obtain biologically meaningful insights.

Computational approaches for characterizing the tumor immune microenvironment

Recent advances in high-throughput molecular profiling technologies and multiplexed imaging platforms have revolutionized our ability to characterize the tumor immune microenvironment. As a result, studies of tumor-associated immune cells increasingly involve complex data sets that require sophisticated methods of computational analysis. In this review, we present an overview of key assays and related bioinformatics tools for analyzing the tumor-associated immune system in bulk tissues and at the single-cell level. In parallel, we describe how data science strategies and novel technologies have advanced tumor immunology and opened the door for new opportunities to exploit host immunity to improve cancer clinical outcomes.

Nascent transcript analysis of glucocorticoid crosstalk with TNF defines primary and cooperative inflammatory repression

The glucocorticoid receptor (NR3C1, also known as GR) binds to specific DNA sequences and directly induces transcription of anti-inflammatory genes that contribute to cytokine repression, frequently in cooperation with NF-kB. Whether inflammatory repression also occurs through local interactions between GR and inflammatory gene regulatory elements has been controversial. Here, using global run-on sequencing (GRO-seq) in human airway epithelial cells, we show that glucocorticoid signaling represses transcription within 10 min. Many repressed regulatory regions reside within “hyper-ChIPable” genomic regions that are subject to dynamic, yet nonspecific, interactions with some antibodies. When this artifact was accounted for, we determined that transcriptional repression does not require local GR occupancy. Instead, widespread transcriptional induction through canonical GR binding sites is associated with reciprocal repression of distal TNF-regulated enhancers through a chromatin-dependent process, as evidenced by chromatin accessibility and motif displacement analysis. Simultaneously, transcriptional induction of key anti-inflammatory effectors is decoupled from primary repression through cooperation between GR and NF-kB at a subset of regulatory regions. Thus, glucocorticoids exert bimodal restraints on inflammation characterized by rapid primary transcriptional repression without local GR occupancy and secondary anti-inflammatory effects resulting from transcriptional cooperation between GR and NF-kB.

A stable mode of bookmarking by TBP recruits RNA polymerase II to mitotic chromosomes

Maintenance of transcription programs is challenged during mitosis when chromatin becomes condensed and transcription is silenced. How do the daughter cells re-establish the original transcription program? Here, we report that the TATA-binding protein (TBP), a key component of the core transcriptional machinery, remains bound globally to active promoters in mouse embryonic stem cells during mitosis. Using live-cell single-molecule imaging, we observed that TBP mitotic binding is highly stable, with an average residence time of minutes, in stark contrast to typical TFs with residence times of seconds. To test the functional effect of mitotic TBP binding, we used a drug-inducible degron system and found that TBP promotes the association of RNA Polymerase II with mitotic chromosomes, and facilitates transcriptional reactivation following mitosis. These results suggest that the core transcriptional machinery promotes efficient transcription maintenance globally.