TOP1 inhibition therapy protects against SARS-CoV-2-induced lethal inflammation.

The ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently affecting millions of lives worldwide. Large retrospective studies indicate that an elevated level of inflammatory cytokines and pro-inflammatory factors are associated with both increased disease severity and mortality. Here, using multidimensional epigenetic, transcriptional, in vitro, and in vivo analyses, we report that topoisomerase 1 (TOP1) inhibition suppresses lethal inflammation induced by SARS-CoV-2. Therapeutic treatment with two doses of topotecan (TPT), an FDA-approved TOP1 inhibitor, suppresses infection-induced inflammation in hamsters. TPT treatment as late as 4 days post-infection reduces morbidity and rescues mortality in a transgenic mouse model. These results support the potential of TOP1 inhibition as an effective host-directed therapy against severe SARS-CoV-2 infection. TPT and its derivatives are inexpensive clinical-grade inhibitors available in most countries. Clinical trials are needed to evaluate the efficacy of repurposing TOP1 inhibitors for severe coronavirus disease 2019 (COVID-19) in humans.

Genomic Decoding of Neuronal Depolarization by Stimulus-Specific NPAS4 Heterodimers.

Cells regulate gene expression in response to salient external stimuli. In neurons, depolarization leads to the expression of inducible transcription factors (ITFs) that direct subsequent gene regulation. Depolarization encodes both a neuron's action potential (AP) output and synaptic inputs, via excitatory postsynaptic potentials (EPSPs). However, it is unclear if distinct types of electrical activity can be transformed by an ITF into distinct modes of genomic regulation. Here, we show that APs and EPSPs in mouse hippocampal neurons trigger two spatially segregated and molecularly distinct induction mechanisms that lead to the expression of the ITF NPAS4. These two pathways culminate in the formation of stimulus-specific NPAS4 heterodimers that exhibit distinct DNA binding patterns. Thus, NPAS4 differentially communicates increases in a neuron's spiking output and synaptic inputs to the nucleus, enabling gene regulation to be tailored to the type of depolarizing activity along the somato-dendritic axis of a neuron.

Transcription Elongation Can Affect Genome 3D Structure.

How transcription affects genome 3D organization is not well understood. We found that during influenza A (IAV) infection, rampant transcription rapidly reorganizes host cell chromatin interactions. These changes occur at the ends of highly transcribed genes, where global inhibition of transcription termination by IAV NS1 protein causes readthrough transcription for hundreds of kilobases. In these readthrough regions, elongating RNA polymerase II disrupts chromatin interactions by inducing cohesin displacement from CTCF sites, leading to locus decompaction. Readthrough transcription into heterochromatin regions switches them from the inert (B) to the permissive (A) chromatin compartment and enables transcription factor binding. Data from non-viral transcription stimuli show that transcription similarly affects cohesin-mediated chromatin contacts within gene bodies. Conversely, inhibition of transcription elongation allows cohesin to accumulate at previously transcribed intragenic CTCF sites and to mediate chromatin looping and compaction. Our data indicate that transcription elongation by RNA polymerase II remodels genome 3D architecture.

Position-dependent function of human sequence-specific transcription factors.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Liver-Derived Signals Sequentially Reprogram Myeloid Enhancers to Initiate and Maintain Kupffer Cell Identity.

Tissue environment plays a powerful role in establishing and maintaining the distinct phenotypes of resident macrophages, but the underlying molecular mechanisms remain poorly understood. Here, we characterized transcriptomic and epigenetic changes in repopulating liver macrophages following acute Kupffer cell depletion as a means to infer signaling pathways and transcription factors that promote Kupffer cell differentiation. We obtained evidence that combinatorial interactions of the Notch ligand DLL4 and transforming growth factor-b (TGF-ß) family ligands produced by sinusoidal endothelial cells and endogenous LXR ligands were required for the induction and maintenance of Kupffer cell identity. DLL4 regulation of the Notch transcriptional effector RBPJ activated poised enhancers to rapidly induce LXRα and other Kupffer cell lineage-determining factors. These factors in turn reprogrammed the repopulating liver macrophage enhancer landscape to converge on that of the original resident Kupffer cells. Collectively, these findings provide a framework for understanding how macrophage progenitor cells acquire tissue-specific phenotypes.

Senataxin suppresses the antiviral transcriptional response and controls viral biogenesis.

The human helicase senataxin (SETX) has been linked to the neurodegenerative diseases amyotrophic lateral sclerosis (ALS4) and ataxia with oculomotor apraxia (AOA2). Here we identified a role for SETX in controlling the antiviral response. Cells that had undergone depletion of SETX and SETX-deficient cells derived from patients with AOA2 had higher expression of antiviral mediators in response to infection than did wild-type cells. Mechanistically, we propose a model whereby SETX attenuates the activity of RNA polymerase II (RNAPII) at genes stimulated after a virus is sensed and thus controls the magnitude of the host response to pathogens and the biogenesis of various RNA viruses (e.g., influenza A virus and West Nile virus). Our data indicate a potentially causal link among inborn errors in SETX, susceptibility to infection and the development of neurologic disorders.

The selection and function of cell type-specific enhancers.

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

Heterochromatin-Encoded Satellite RNAs Induce Breast Cancer.

Heterochromatic repetitive satellite RNAs are extensively transcribed in a variety of human cancers, including BRCA1 mutant breast cancer. Aberrant expression of satellite RNAs in cultured cells induces the DNA damage response, activates cell cycle checkpoints, and causes defects in chromosome segregation. However, the mechanism by which satellite RNA expression leads to genomic instability is not well understood. Here we provide evidence that increased levels of satellite RNAs in mammary glands induce tumor formation in mice. Using mass spectrometry, we further show that genomic instability induced by satellite RNAs occurs through interactions with BRCA1-associated protein networks required for the stabilization of DNA replication forks. Additionally, de-stabilized replication forks likely promote the formation of RNA-DNA hybrids in cells expressing satellite RNAs. These studies lay the foundation for developing novel therapeutic strategies that block the effects of non-coding satellite RNAs in cancer cells.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation.

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Global changes in the nuclear positioning of genes and intra- and interdomain genomic interactions that orchestrate B cell fate.

The genome is folded into domains located in compartments that are either transcriptionally inert or transcriptionally permissive. Here we used genome-wide strategies to characterize domains during B cell development. Structured interaction matrix analysis showed that occupancy by the architectural protein CTCF was associated mainly with intradomain interactions, whereas sites bound by the histone acetyltransferase p300 or the transcription factors E2A or PU.1 were associated with intra- and interdomain interactions that are developmentally regulated. We identified a spectrum of genes that switched nuclear location during early B cell development. In progenitor cells, the transcriptionally inactive locus encoding early B cell factor (Ebf1) was sequestered at the nuclear lamina, which thereby preserved their multipotency. After development into the pro-B cell stage, Ebf1 and other genes switched compartments to establish new intra- and interdomain interactions associated with a B lineage-specific transcription signature.

Chromosome 10q26-driven age-related macular degeneration is associated with reduced levels of HTRA1 in human retinal pigment epithelium.

Genome-wide association studies have identified the chromosome 10q26 (Chr10) locus, which contains the age-related maculopathy susceptibility 2 (ARMS2) and high temperature requirement A serine peptidase 1 (HTRA1) genes, as the strongest genetic risk factor for age-related macular degeneration (AMD) [L.G. Fritsche et al., Annu. Rev. Genomics Hum. Genet. 15, 151-171, (2014)]. To date, it has been difficult to assign causality to any specific single nucleotide polymorphism (SNP), haplotype, or gene within this region because of high linkage disequilibrium among the disease-associated variants [J. Jakobsdottir et al. Am. J. Hum. Genet. 77, 389-407 (2005); A. Rivera et al. Hum. Mol. Genet. 14, 3227-3236 (2005)]. Here, we show that HTRA1 messenger RNA (mRNA) is reduced in retinal pigment epithelium (RPE) but not in neural retina or choroid tissues derived from human donors with homozygous risk at the 10q26 locus. This tissue-specific decrease is mediated by the presence of a noncoding, cis-regulatory element overlapping the ARMS2 intron, which contains a potential Lhx2 transcription factor binding site that is disrupted by risk variant rs36212733. HtrA1 protein increases with age in the RPE-Bruch's membrane (BM) interface in Chr10 nonrisk donors but fails to increase in donors with homozygous risk at the 10q26 locus. We propose that HtrA1, an extracellular chaperone and serine protease, functions to maintain the optimal integrity of the RPE-BM interface during the aging process and that reduced expression of HTRA1 mRNA and protein in Chr10 risk donors impairs this protective function, leading to increased risk of AMD pathogenesis. HtrA1 augmentation, not inhibition, in high-risk patients should be considered as a potential therapy for AMD.

Human promoters are intrinsically directional.

Divergent transcription, in which reverse-oriented transcripts occur upstream of eukaryotic promoters in regions devoid of annotated genes, has been suggested to be a general property of active promoters. Here we show that the human basal RNA polymerase II transcriptional machinery and core promoter are inherently unidirectional and that reverse-oriented transcripts originate from their own cognate reverse-directed core promoters. In vitro transcription analysis and mapping of nascent transcripts in HeLa cells revealed that sequences at reverse start sites are similar to those of their forward counterparts. The use of DNase I accessibility to define proximal promoter borders revealed that about half of promoters are unidirectional and that unidirectional promoters are depleted at their upstream edges of reverse core promoter sequences and their associated chromatin features. Divergent transcription is thus not an inherent property of the transcription process but rather the consequence of the presence of both forward- and reverse-directed core promoters.

Identification and dynamic quantification of regulatory elements using total RNA.

The spatial and temporal regulation of transcription initiation is pivotal for controlling gene expression. Here, we introduce capped-small RNA-seq (csRNA-seq), which uses total RNA as starting material to detect transcription start sites (TSSs) of both stable and unstable RNAs at single-nucleotide resolution. csRNA-seq is highly sensitive to acute changes in transcription and identifies an order of magnitude more regulated transcripts than does RNA-seq. Interrogating tissues from species across the eukaryotic kingdoms identified unstable transcripts resembling enhancer RNAs, pri-miRNAs, antisense transcripts, and promoter upstream transcripts in multicellular animals, plants, and fungi spanning 1.6 billion years of evolution. Integration of epigenomic data from these organisms revealed that histone H3 trimethylation (H3K4me3) was largely confined to TSSs of stable transcripts, whereas H3K27ac marked nucleosomes downstream from all active TSSs, suggesting an ancient role for posttranslational histone modifications in transcription. Our findings show that total RNA is sufficient to identify transcribed regulatory elements and capture the dynamics of initiated stable and unstable transcripts at single-nucleotide resolution in eukaryotes.

A global network of transcription factors, involving E2A, EBF1 and Foxo1, that orchestrates B cell fate.

It is now established that the transcription factors E2A, EBF1 and Foxo1 have critical roles in B cell development. Here we show that E2A and EBF1 bound regulatory elements present in the Foxo1 locus. E2A and EBF1, as well as E2A and Foxo1, in turn, were wired together by a vast spectrum of cis-regulatory sequences. These associations were dynamic during developmental progression. Occupancy by the E2A isoform E47 directly resulted in greater abundance, as well as a pattern of monomethylation of histone H3 at lysine 4 (H3K4) across putative enhancer regions. Finally, we divided the pro-B cell epigenome into clusters of loci with occupancy by E2A, EBF and Foxo1. From this analysis we constructed a global network consisting of transcriptional regulators, signaling and survival factors that we propose orchestrates B cell fate.

Genome evolution in the allotetraploid frog Xenopus laevis.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Paradoxical association of TET loss of function with genome-wide DNA hypomethylation.

Cancer genomes are characterized by focal increases in DNA methylation, co-occurring with widespread hypomethylation. Here, we show that TET loss of function results in a similar genomic footprint. Both 5hmC in wild-type (WT) genomes and DNA hypermethylation in TET-deficient genomes are largely confined to the active euchromatic compartment, consistent with the known functions of TET proteins in DNA demethylation and the known distribution of 5hmC at transcribed genes and active enhancers. In contrast, an unexpected DNA hypomethylation noted in multiple TET-deficient genomes is primarily observed in the heterochromatin compartment. In a mouse model of T cell lymphoma driven by TET deficiency (Tet2/3 DKO T cells), genomic analysis of malignant T cells revealed DNA hypomethylation in the heterochromatic genomic compartment, as well as reactivation of repeat elements and enrichment for single-nucleotide alterations, primarily in heterochromatic regions of the genome. Moreover, hematopoietic stem/precursor cells (HSPCs) doubly deficient for Tet2 and Dnmt3a displayed greater losses of DNA methylation than HSPCs singly deficient for Tet2 or Dnmt3a alone, potentially explaining the unexpected synergy between DNMT3A and TET2 mutations in myeloid and lymphoid malignancies. Tet1-deficient cells showed decreased localization of DNMT3A in the heterochromatin compartment compared with WT cells, pointing to a functional interaction between TET and DNMT proteins and providing a potential explanation for the hypomethylation observed in TET-deficient genomes. Our data suggest that TET loss of function may at least partially underlie the characteristic pattern of global hypomethylation coupled to regional hypermethylation observed in diverse cancer genomes, and highlight the potential contribution of heterochromatin hypomethylation to oncogenesis.

Identification of Lineage-Specific Transcription Factors That Prevent Activation of Hepatic Stellate Cells and Promote Fibrosis Resolution.

BACKGROUND & AIMS: Development of liver fibrosis is associated with activation of quiescent hepatic stellate cells (HSCs) into collagen type I-producing myofibroblasts (activated HSCs). Cessation of liver injury often results in fibrosis resolution and inactivation of activated HSCs/myofibroblasts into a quiescent-like state (inactivated HSCs). We aimed to identify molecular features of phenotypes of HSCs from mice and humans. METHODS: We performed studies with LratCre, Ets1-floxed, Nf1-floxed, Pparγ-floxed, Gata6-floxed, Rag2-/-γc-/-, and C57/Bl6 (control) mice. Some mice were given carbon tetrachloride (CCl4) to induce liver fibrosis, with or without a peroxisome proliferator-activated receptor-γ (PPARγ) agonist. Livers from mice were analyzed by immunohistochemistry. Quiescent, activated, and inactivated HSCs were isolated from livers of Col1α1YFP mice and analyzed by chromatin immunoprecipitation and sequencing. Human HSCs were isolated from livers denied for transplantation. We compared changes in gene expression patterns and epigenetic modifications (histone H3 lysine 4 dimethylation and histone H3 lysine 27 acetylation) in primary mouse and human HSCs. Transcription factors were knocked down with small hairpin RNAs in mouse HSCs. RESULTS: Motif enrichment identified E26 transcription-specific transcription factors (ETS) 1, ETS2, GATA4, GATA6, interferon regulatory factor (IRF) 1, and IRF2 transcription factors as regulators of the mouse and human HSC lineage. Small hairpin RNA-knockdown of these transcription factors resulted in increased expression of genes that promote fibrogenesis and inflammation, and loss of HSC phenotype. Disruption of Gata6 or Ets1, or Nf1 or Pparγ (which are regulated by ETS1), increased the severity of CCl4-induced liver fibrosis in mice compared to control mice. Only mice with disruption of Gata6 or Pparγ had defects in fibrosis resolution after CCl4 administration was stopped, associated with persistent activation of HSCs. Administration of a PPARγ agonist accelerated regression of liver fibrosis after CCl4 administration in control mice but not in mice with disruption of Pparγ. CONCLUSIONS: Phenotypes of HSCs from humans and mice are regulated by transcription factors, including ETS1, ETS2, GATA4, GATA6, IRF1, and IRF2. Activated mouse and human HSCs can revert to a quiescent-like, inactivated phenotype. We found GATA6 and PPARγ to be required for inactivation of human HSCs and regression of liver fibrosis in mice.

Remodeling of the enhancer landscape during macrophage activation is coupled to enhancer transcription.

Recent studies suggest a hierarchical model in which lineage-determining factors act in a collaborative manner to select and prime cell-specific enhancers, thereby enabling signal-dependent transcription factors to bind and function in a cell-type-specific manner. Consistent with this model, TLR4 signaling primarily regulates macrophage gene expression through a pre-existing enhancer landscape. However, TLR4 signaling also induces priming of ∼3,000 enhancer-like regions de novo, enabling visualization of intermediates in enhancer selection and activation. Unexpectedly, we find that enhancer transcription precedes local mono- and dimethylation of histone H3 lysine 4 (H3K4me1/2). H3K4 methylation at de novo enhancers is primarily dependent on the histone methyltransferases Mll1, Mll2/4, and Mll3 and is significantly reduced by inhibition of RNA polymerase II elongation. Collectively, these findings suggest an essential role of enhancer transcription in H3K4me1/2 deposition at de novo enhancers that is independent of potential functions of the resulting eRNA transcripts.

Deconvolution of pro- and antiviral genomic responses in Zika virus-infected and bystander macrophages.

Genome-wide investigations of host-pathogen interactions are often limited by analyses of mixed populations of infected and uninfected cells, which lower sensitivity and accuracy. To overcome these obstacles and identify key mechanisms by which Zika virus (ZIKV) manipulates host responses, we developed a system that enables simultaneous characterization of genome-wide transcriptional and epigenetic changes in ZIKV-infected and neighboring uninfected primary human macrophages. We demonstrate that transcriptional responses in ZIKV-infected macrophages differed radically from those in uninfected neighbors and that studying the cell population as a whole produces misleading results. Notably, the uninfected population of macrophages exhibits the most rapid and extensive changes in gene expression, related to type I IFN signaling. In contrast, infected macrophages exhibit a delayed and attenuated transcriptional response distinguished by preferential expression of IFNB1 at late time points. Biochemical and genomic studies of infected macrophages indicate that ZIKV infection causes both a targeted defect in the type I IFN response due to degradation of STAT2 and reduces RNA polymerase II protein levels and DNA occupancy, particularly at genes required for macrophage identity. Simultaneous evaluation of transcriptomic and epigenetic features of infected and uninfected macrophages thereby reveals the coincident evolution of dominant proviral or antiviral mechanisms, respectively, that determine the outcome of ZIKV exposure.