Control of inducible gene expression links cohesin to hematopoietic progenitor self-renewal and differentiation.

Cohesin is important for 3D genome organization. Nevertheless, even the complete removal of cohesin has surprisingly little impact on steady-state gene transcription and enhancer activity. Here we show that cohesin is required for the core transcriptional response of primary macrophages to microbial signals, and for inducible enhancer activity that underpins inflammatory gene expression. Consistent with a role for inflammatory signals in promoting myeloid differentiation of hematopoietic stem and progenitor cells (HPSCs), cohesin mutations in HSPCs led to reduced inflammatory gene expression and increased resistance to differentiation-inducing inflammatory stimuli. These findings uncover an unexpected dependence of inducible gene expression on cohesin, link cohesin with myeloid differentiation, and may help explain the prevalence of cohesin mutations in human acute myeloid leukemia.

Competition between PRC2.1 and 2.2 subcomplexes regulates PRC2 chromatin occupancy in human stem cells.

Polycomb repressive complex 2 (PRC2) silences expression of developmental transcription factors in pluripotent stem cells by methylating lysine 27 on histone H3. Two mutually exclusive subcomplexes, PRC2.1 and PRC2.2, are defined by the set of accessory proteins bound to the core PRC2 subunits. Here we introduce separation-of-function mutations into the SUZ12 subunit of PRC2 to drive it into a PRC2.1 or 2.2 subcomplex in human induced pluripotent stem cells (iPSCs). We find that PRC2.2 occupies polycomb target genes at low levels and that homeobox transcription factors are upregulated when this complex is exclusively present. In contrast with previous studies, we find that chromatin occupancy of PRC2 increases drastically when it is forced to form PRC2.1. Additionally, several cancer-associated mutations also coerce formation of PRC2.1. We suggest that PRC2 chromatin occupancy can be altered in the context of disease or development by tuning the ratio of PRC2.1 to PRC2.2.

HIF1A employs CDK8-mediator to stimulate RNAPII elongation in response to hypoxia.

The transcription factor HIF1A is a key mediator of the cellular response to hypoxia. Despite the importance of HIF1A in homeostasis and various pathologies, little is known about how it regulates RNA polymerase II (RNAPII). We report here that HIF1A employs a specific variant of the Mediator complex to stimulate RNAPII elongation. The Mediator-associated kinase CDK8, but not the paralog CDK19, is required for induction of many HIF1A target genes. HIF1A induces binding of CDK8-Mediator and the super elongation complex (SEC), containing AFF4 and CDK9, to alleviate RNAPII pausing. CDK8 is dispensable for HIF1A chromatin binding and histone acetylation, but it is essential for binding of SEC and RNAPII elongation. Global analysis of active RNAPII reveals that hypoxia-inducible genes are paused and active prior to their induction. Our results provide a mechanistic link between HIF1A and CDK8, two potent oncogenes, in the cellular response to hypoxia.

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.

Selective inhibition of CDK7 reveals high-confidence targets and new models for TFIIH function in transcription.

CDK7 associates with the 10-subunit TFIIH complex and regulates transcription by phosphorylating the C-terminal domain (CTD) of RNA polymerase II (RNAPII). Few additional CDK7 substrates are known. Here, using the covalent inhibitor SY-351 and quantitative phosphoproteomics, we identified CDK7 kinase substrates in human cells. Among hundreds of high-confidence targets, the vast majority are unique to CDK7 (i.e., distinct from other transcription-associated kinases), with a subset that suggest novel cellular functions. Transcription-associated factors were predominant CDK7 substrates, including SF3B1, U2AF2, and other splicing components. Accordingly, widespread and diverse splicing defects, such as alternative exon inclusion and intron retention, were characterized in CDK7-inhibited cells. Combined with biochemical assays, we establish that CDK7 directly activates other transcription-associated kinases CDK9, CDK12, and CDK13, invoking a "master regulator" role in transcription. We further demonstrate that TFIIH restricts CDK7 kinase function to the RNAPII CTD, whereas other substrates (e.g., SPT5 and SF3B1) are phosphorylated by the three-subunit CDK-activating kinase (CAK; CCNH, MAT1, and CDK7). These results suggest new models for CDK7 function in transcription and implicate CAK dissociation from TFIIH as essential for kinase activation. This straightforward regulatory strategy ensures CDK7 activation is spatially and temporally linked to transcription, and may apply toward other transcription-associated kinases.

Transcriptional Responses to IFN-γ Require Mediator Kinase-Dependent Pause Release and Mechanistically Distinct CDK8 and CDK19 Functions.

Transcriptional responses to external stimuli remain poorly understood. Using global nuclear run-on followed by sequencing (GRO-seq) and precision nuclear run-on sequencing (PRO-seq), we show that CDK8 kinase activity promotes RNA polymerase II pause release in response to interferon-γ (IFN-γ), a universal cytokine involved in immunity and tumor surveillance. The Mediator kinase module contains CDK8 or CDK19, which are presumed to be functionally redundant. We implemented cortistatin A, chemical genetics, transcriptomics, and other methods to decouple their function while assessing enzymatic versus structural roles. Unexpectedly, CDK8 and CDK19 regulated different gene sets via distinct mechanisms. CDK8-dependent regulation required its kinase activity, whereas CDK19 governed IFN-γ responses through its scaffolding function (i.e., it was kinase independent). Accordingly, CDK8, not CDK19, phosphorylates the STAT1 transcription factor (TF) during IFN-γ stimulation, and CDK8 kinase inhibition blocked activation of JAK-STAT pathway TFs. Cytokines such as IFN-γ rapidly mobilize TFs to "reprogram" cellular transcription; our results implicate CDK8 and CDK19 as essential for this transcriptional reprogramming.

Macrophages derived from human induced pluripotent stem cells (iPSCs) serve as a high-fidelity cellular model for investigating HIV-1, dengue, and influenza viruses.

Macrophages are important target cells for diverse viruses and thus represent a valuable system for studying virus biology. Isolation of primary human macrophages is done by culture of dissociated tissues or from differentiated blood monocytes, but these methods are both time consuming and result in low numbers of recovered macrophages. Here, we explore whether macrophages derived from human induced pluripotent stem cells (iPSCs)-which proliferate indefinitely and potentially provide unlimited starting material-could serve as a faithful model system for studying virus biology. Human iPSC-derived monocytes were differentiated into macrophages and then infected with HIV-1, dengue virus, or influenza virus as model human viruses. We show that iPSC-derived macrophages support the replication of these viruses with kinetics and phenotypes similar to human blood monocyte-derived macrophages. These iPSC-derived macrophages were virtually indistinguishable from human blood monocyte-derived macrophages based on surface marker expression (flow cytometry), transcriptomics (RNA sequencing), and chromatin accessibility profiling. iPSC lines were additionally generated from non-human primate (chimpanzee) fibroblasts. When challenged with dengue virus, human and chimpanzee iPSC-derived macrophages show differential susceptibility to infection, thus providing a valuable resource for studying the species-tropism of viruses. We also show that blood- and iPSC-derived macrophages both restrict influenza virus at a late stage of the virus lifecycle. Collectively, our results substantiate iPSC-derived macrophages as an alternative to blood monocyte-derived macrophages for the study of virus biology. IMPORTANCE: Macrophages have complex relationships with viruses: while macrophages aid in the removal of pathogenic viruses from the body, macrophages are also manipulated by some viruses to serve as vessels for viral replication, dissemination, and long-term persistence. Here, we show that iPSC-derived macrophages are an excellent model that can be exploited in virology.

Internal and external normalization of nascent RNA sequencing run-on experiments.

In experiments with significant perturbations to transcription, nascent RNA sequencing protocols are dependent on external spike-ins for reliable normalization. Unlike in RNA-seq, these spike-ins are not standardized and, in many cases, depend on a run-on reaction that is assumed to have constant efficiency across samples. To assess the validity of this assumption, we analyze a large number of published nascent RNA spike-ins to quantify their variability across existing normalization methods. Furthermore, we develop a new biologically-informed Bayesian model to estimate the error in spike-in based normalization estimates, which we term Virtual Spike-In (VSI). We apply this method both to published external spike-ins as well as using reads at the [Formula: see text] end of long genes, building on prior work from Mahat (Mol Cell 62(1):63-78, 2016. https://doi.org/10.1016/j.molcel.2016.02.025 ) and Vihervaara (Nat Commun 8(1):255, 2017. https://doi.org/10.1038/s41467-017-00151-0 ). We find that spike-ins in existing nascent RNA experiments are typically under sequenced, with high variability between samples. Furthermore, we show that these high variability estimates can have significant downstream effects on analysis, complicating biological interpretations of results.

The Δ40p53 isoform inhibits p53-dependent eRNA transcription and enables regulation by signal-specific transcription factors during p53 activation.

The naturally occurring Δ40p53 isoform heterotetramerizes with wild-type p53 (WTp53) to regulate development, aging, and stress responses. How Δ40p53 alters WTp53 function remains enigmatic because their co-expression causes tetramer heterogeneity. We circumvented this issue with a well-tested strategy that expressed Δ40p53:WTp53 as a single transcript, ensuring a 2:2 tetramer stoichiometry. Human MCF10A cell lines expressing Δ40p53:WTp53, WTp53, or WTp53:WTp53 (as controls) from the native TP53 locus were examined with transcriptomics (precision nuclear run-on sequencing [PRO-seq] and RNA sequencing [RNA-seq]), metabolomics, and other methods. Δ40p53:WTp53 was transcriptionally active, and, although phenotypically similar to WTp53 under normal conditions, it failed to induce growth arrest upon Nutlin-induced p53 activation. This occurred via Δ40p53:WTp53-dependent inhibition of enhancer RNA (eRNA) transcription and subsequent failure to induce mRNA biogenesis, despite similar genomic occupancy to WTp53. A different stimulus (5-fluorouracil [5FU]) also showed Δ40p53:WTp53-specific changes in mRNA induction; however, other transcription factors (TFs; e.g., E2F2) could then drive the response, yielding similar outcomes vs. WTp53. Our results establish that Δ40p53 tempers WTp53 function to enable compensatory responses by other stimulus-specific TFs. Such modulation of WTp53 activity may be an essential physiological function for Δ40p53. Moreover, Δ40p53:WTp53 functional distinctions uncovered herein suggest an eRNA requirement for mRNA biogenesis and that human p53 evolved as a tetramer to support eRNA transcription.

Just 2% of SARS-CoV-2-positive individuals carry 90% of the virus circulating in communities.

We analyze data from the fall 2020 pandemic response efforts at the University of Colorado Boulder, where more than 72,500 saliva samples were tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using qRT-PCR. All samples were collected from individuals who reported no symptoms associated with COVID-19 on the day of collection. From these, 1,405 positive cases were identified. The distribution of viral loads within these asymptomatic individuals was indistinguishable from what has been previously observed in symptomatic individuals. Regardless of symptomatic status, ∼50% of individuals who test positive for SARS-CoV-2 seem to be in noninfectious phases of the disease, based on having low viral loads in a range from which live virus has rarely been isolated. We find that, at any given time, just 2% of individuals carry 90% of the virions circulating within communities, serving as viral "supercarriers" and possibly also superspreaders.

Transcription dosage compensation does not occur in Down syndrome.

BACKGROUND: The increase in DNA copy number in Down syndrome (DS; caused by trisomy 21) has led to the DNA dosage hypothesis, which posits that the level of gene expression is proportional to the gene's DNA copy number. Yet many reports have suggested that a proportion of chromosome 21 genes are dosage compensated back towards typical expression levels (1.0×). In contrast, other reports suggest that dosage compensation is not a common mechanism of gene regulation in trisomy 21, providing support to the DNA dosage hypothesis. RESULTS: In our work, we use both simulated and real data to dissect the elements of differential expression analysis that can lead to the appearance of dosage compensation, even when compensation is demonstrably absent. Using lymphoblastoid cell lines derived from a family with an individual with Down syndrome, we demonstrate that dosage compensation is nearly absent at both nascent transcription (GRO-seq) and steady-state RNA (RNA-seq) levels. Furthermore, we link the limited apparent dosage compensation to expected allelic variation in transcription levels. CONCLUSIONS: Transcription dosage compensation does not occur in Down syndrome. Simulated data containing no dosage compensation can appear to have dosage compensation when analyzed via standard methods. Moreover, some chromosome 21 genes that appear to be dosage compensated are consistent with allele specific expression.

Integrated genomics approaches identify transcriptional mediators and epigenetic responses to Afghan desert particulate matter in small airway epithelial cells.

Military Deployment to Southwest Asia and Afghanistan and exposure to toxic airborne particulates have been associated with an increased risk of developing respiratory disease, collectively termed deployment-related respiratory diseases (DRRDs). Our knowledge about how particulates mediate respiratory disease is limited, precluding the appropriate recognition or management. Central to this limitation is the lack of understanding of how exposures translate into dysregulated cell identity with dysregulated transcriptional programs. The small airway epithelium is involved in both the pathobiology of DRRD and fine particulate matter deposition. To characterize small airway epithelial cell epigenetic and transcriptional responses to Afghan desert particulate matter (APM) and investigate the functional interactions of transcription factors that mediate these responses, we applied two genomics assays, the assay for transposase accessible chromatin with sequencing (ATAC-seq) and Precision Run-on sequencing (PRO-seq). We identified activity changes in a series of transcriptional pathways as candidate regulators of susceptibility to subsequent insults, including signal-dependent pathways, such as loss of cytochrome P450 or P53/P63, and lineage-determining transcription factors, such as GRHL2 loss or TEAD3 activation. We further demonstrated that TEAD3 activation was unique to APM exposure despite similar inflammatory responses when compared with wood smoke particle exposure and that P53/P63 program loss was uniquely positioned at the intersection of signal-dependent and lineage-determining transcriptional programs. Our results establish the utility of an integrated genomics approach in characterizing responses to exposures and identifying genomic targets for the advanced investigation of the pathogenesis of DRRD.

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.

Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries.

BACKGROUND: A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries. RESULTS: Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ' end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription. CONCLUSIONS: Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

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.

Characterization of Global Gene Expression, Regulation of Metal Ions, and Infection Outcomes in Immune-Competent 129S6 Mouse Macrophages.

Nutritional immunity involves cellular and physiological responses to invading pathogens, such as limiting iron, increasing exposure to bactericidal copper, and altering zinc to restrict the growth of pathogens. Here, we examine infection of bone marrow-derived macrophages from 129S6/SvEvTac mice by Salmonella enterica serovar Typhimurium. The 129S6/SvEvTac mice possess a functional Slc11a1 (Nramp-1), a phagosomal transporter of divalent cations that plays an important role in modulating metal availability to the pathogen. We carried out global RNA sequencing upon treatment with live or heat-killed Salmonella at 2 h and 18 h postinfection and observed widespread changes in metal transport, metal-dependent genes, and metal homeostasis genes, suggesting significant remodeling of iron, copper, and zinc availability by host cells. Changes in host cell gene expression suggest infection increases cytosolic zinc while simultaneously limiting zinc within the phagosome. Using a genetically encoded sensor, we demonstrate that cytosolic labile zinc increases 45-fold at 12 h postinfection. Further, manipulation of zinc in the medium alters bacterial clearance and replication, with zinc depletion inhibiting both processes. Comparing the transcriptomic changes to published data on infection of C57BL/6 macrophages revealed notable differences in metal regulation and the global immune response. Our results reveal that 129S6 macrophages represent a distinct model system compared to C57BL/6 macrophages. Further, our results indicate that manipulation of zinc at the host-pathogen interface is more nuanced than that of iron or copper. The 129S6 macrophages leverage intricate means of manipulating zinc availability and distribution to limit the pathogen's access to zinc, while simultaneously ensuring sufficient zinc to support the immune response.

Enhancer RNA profiling predicts transcription factor activity.

Transcription factors (TFs) exert their regulatory influence through the binding of enhancers, resulting in coordination of gene expression programs. Active enhancers are often characterized by the presence of short, unstable transcripts termed enhancer RNAs (eRNAs). While their function remains unclear, we demonstrate that eRNAs are a powerful readout of TF activity. We infer sites of eRNA origination across hundreds of publicly available nascent transcription data sets and show that eRNAs initiate from sites of TF binding. By quantifying the colocalization of TF binding motif instances and eRNA origins, we derive a simple statistic capable of inferring TF activity. In doing so, we uncover dozens of previously unexplored links between diverse stimuli and the TFs they affect.

Genome-wide dose-dependent inhibition of histone deacetylases studies reveal their roles in enhancer remodeling and suppression of oncogenic super-enhancers.

Histone deacetylase inhibitors (HDACIs) are known to alter gene expression by both up- and down-regulation of protein-coding genes in normal and cancer cells. However, the exact regulatory mechanisms of action remain uncharacterized. Here we investigated genome wide dose-dependent epigenetic and transcriptome changes in response to HDACI largazole in a transformed and a non-transformed cell line. Exposure to low nanomolar largazole concentrations (

A trans-acting Variant within the Transcription Factor RIM101 Interacts with Genetic Background to Determine its Regulatory Capacity.

Most genetic variants associated with disease occur within regulatory regions of the genome, underscoring the importance of defining the mechanisms underlying differences in regulation of gene expression between individuals. We discovered a pair of co-regulated, divergently oriented transcripts, AQY2 and ncFRE6, that are expressed in one strain of Saccharomyces cerevisiae, ∑1278b, but not in another, S288c. By combining classical genetics techniques with high-throughput sequencing, we identified a trans-acting single nucleotide polymorphism within the transcription factor RIM101 that causes the background-dependent expression of both transcripts. Subsequent RNA-seq experiments revealed that RIM101 regulates many more targets in S288c than in ∑1278b and that deletion of RIM101 in both backgrounds abrogates the majority of differential expression between the strains. Strikingly, only three transcripts undergo a significant change in expression after swapping RIM101 alleles between backgrounds, implying that the differences in the RIM101 allele lead to a remarkably focused transcriptional response. However, hundreds of RIM101-dependent targets undergo a subtle but consistent shift in expression in the S288c RIM101-swapped strain, but not its ∑1278b counterpart. We conclude that ∑1278b may harbor a variant(s) that buffers against widespread transcriptional dysregulation upon introduction of a non-native RIM101 allele, emphasizing the importance of accounting for genetic background when assessing the impact of a regulatory variant.

The Influence of Polyploidy on the Evolution of Yeast Grown in a Sub-Optimal Carbon Source.

Polyploidization events have occurred during the evolution of many fungi, plant, and animal species and are thought to contribute to speciation and tumorigenesis, however little is known about how ploidy level contributes to adaptation at the molecular level. Here we integrate whole genome sequencing, RNA expression analysis, and relative fitness of ∼100 evolved clones at three ploidy levels. Independent haploid, diploid, and tetraploid populations were grown in a low carbon environment for 250 generations. We demonstrate that the key adaptive mutation in the evolved clones is predicted by a gene expression signature of just five genes. All of the adaptive mutations identified encompass a narrow set of genes, however the tetraploid clones gain a broader spectrum of adaptive mutations than haploid or diploid clones. While many of the adaptive mutations occur in genes that encode proteins with known roles in glucose sensing and transport, we discover mutations in genes with no canonical role in carbon utilization (IPT1 and MOT3), as well as identify novel dominant mutations in glucose signal transducers thought to only accumulate recessive mutations in carbon limited environments (MTH1 and RGT1). We conclude that polyploid cells explore more genotypic and phenotypic space than lower ploidy cells. Our study provides strong evidence for the beneficial role of polyploidization events that occur during the evolution of many species and during tumorigenesis.