SUMMARY 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), a FDA-approved Top1 inhibitor, suppresses infection-induced inflammation in hamsters. TPT treatment as late as four 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 COVID-19 in humans.
Abstract The conserved Runt-related (RUNX) transcription factor family are well-known master regulators of developmental and regenerative processes. Runx1 and Runx2 are both expressed in satellite cells (SC) and skeletal myotubes. Conditional deletion of Runx1 in adult SC negatively impacted self-renewal and impaired skeletal muscle maintenance. Runx1- deficient SC retain Runx2 expression but cannot support muscle regeneration in response to injury. To determine the unique molecular functions of Runx1 that cannot be compensated by Runx2 we deleted Runx1 in C2C12 that retain Runx2 expression and established that myoblasts differentiation was blocked in vitro due in part to ectopic expression of Mef2c, a target repressed by Runx1 . Structure-function analysis demonstrated that the Ets-interacting MID/EID region of Runx1, absent from Runx2, is critical to regulating myoblasts proliferation, differentiation, and fusion. Analysis of in-house and published ChIP-seq datasets from Runx1 (T-cells, muscle) versus Runx2 (preosteoblasts) dependent tissue identified enrichment for a Ets:Runx composite site in Runx1 -dependent tissues. Comparing ATACseq datasets from WT and Runx1KO C2C12 cells showed that the Ets:Runx composite motif was enriched in peaks open exclusively in WT cells compared to peaks unique to Runx1KO cells. Thus, engagement of a set of targets by the RUNX1/ETS complex define the non-redundant functions of Runx1 .
Atopic dermatitis (AD) is one of the most common skin disorders among children. Disease etiology involves genetic and environmental factors, with 29 independent AD risk loci enriched for risk allele-dependent gene expression in the skin and CD4+ T cell compartments. We investigated the potential epigenetic mechanisms responsible for the genetic susceptibility of CD4+ T cells. To understand the differences in gene regulatory activity in peripheral blood T cells in AD, we measured chromatin accessibility (an assay based on transposase-accessible chromatin sequencing, ATAC-seq), nuclear factor kappa B subunit 1 (NFKB1) binding (chromatin immunoprecipitation with sequencing, ChIP-seq), and gene expression levels (RNA-seq) in stimulated CD4+ T cells from subjects with active moderate-to-severe AD, as well as in age-matched non-allergic controls. Open chromatin regions in stimulated CD4+ T cells were highly enriched for AD genetic risk variants, with almost half of the AD risk loci overlapping AD-dependent ATAC-seq peaks. AD-specific open chromatin regions were strongly enriched for NF-κB DNA-binding motifs. ChIP-seq identified hundreds of NFKB1-occupied genomic loci that were AD- or control-specific. As expected, the AD-specific ChIP-seq peaks were strongly enriched for NF-κB DNA-binding motifs. Surprisingly, control-specific NFKB1 ChIP-seq peaks were not enriched for NFKB1 motifs, but instead contained motifs for other classes of human transcription factors, suggesting a mechanism involving altered indirect NFKB1 binding. Using DNA sequencing data, we identified 63 instances of altered genotype-dependent chromatin accessibility at 36 AD risk variant loci (30% of AD risk loci) that might lead to genotype-dependent gene expression. Based on these findings, we propose that CD4+ T cells respond to stimulation in an AD-specific manner, resulting in disease- and genotype-dependent chromatin accessibility alterations involving NFKB1 binding.
Myeloid lineage cells use TLRs to recognize and respond to diverse microbial ligands. Although unique transcription factors dictate the outcome of specific TLR signaling, whether lineage-specific differences exist to further modulate the quality of TLR-induced inflammation remains unclear. Comprehensive analysis of global gene transcription in human monocytes, monocyte-derived macrophages, and monocyte-derived dendritic cells stimulated with various TLR ligands identifies multiple lineage-specific, TLR-responsive gene programs. Monocytes are hyperresponsive to TLR7/8 stimulation that correlates with the higher expression of the receptors. While macrophages and monocytes express similar levels of TLR4, macrophages, but not monocytes, upregulate interferon-stimulated genes (ISGs) in response to TLR4 stimulation. We find that TLR4 signaling in macrophages uniquely engages transcription factor IRF1, which facilitates the opening of ISG loci for transcription. This study provides a critical mechanistic basis for lineage-specific TLR responses and uncovers IRF1 as a master regulator for the ISG transcriptional program in human macrophages.
Abstract A status-quo in targeted cancer therapy is that out of the thousands of somatic alterations found in a cancer cell, alterations only in driver genes determine therapeutic strategy. Despite unimpressive results of some driver-targeted therapies, and given that the majority of genomic alterations in cancer are not ‘drivers’, but ‘passengers’ / bystander alterations, it remains underappreciated whether targeting built-in vulnerabilities imposed by passenger gene alterations may provide therapeutic value. The tumor suppressor PTEN undergoes widespread functional inactivation including deletion in human cancer. PTEN deletion occurs frequently in GBM, sometimes as part of the 10q loss or chromosome 10 monosomy. We discovered that several genes including Stearoyl Co-A Desaturase, SCD (10q24.31), located 12 MB telomeric to PTEN is frequently co-deleted hemizygously and unintentionally in PTEN-deleted cancers. Strikingly, in a subset of GBM, SCD was also epigenetically silenced. A combination of SCD deletion and methylation resulted in two molecular subgroups – one that expressed SCD, and another that showed little or no detectable SCD. SCD, is an integral membrane protein of the endoplasmic reticulum and converts saturated fatty acids to monounsaturated fatty acids (Oleic and palmitoleic acids) that are critical for membrane fluidity and function, and thus SCD is generally overexpressed in most cancers. We show that SCD expressing lines are highly sensitive to multiple SCD inhibitors, while non-expressors are resistant. Despite modest BBB penetration, one SCD inhibitor was remarkably efficient in blocking intracranial tumor growth. SCD is an oxygen-dependent enzyme. We show that SCD retains significant enzymatic activity even in highly hypoxic conditions. Finally, through RNAseq, functional proteomics and ATACseq, we demonstrate an evolutionarily conserved mechanism of acquired resistance to SCD inhibitor through drug-induced acute phase signaling response in multiple SCD expressing cancers.
Abstract Acute lymphoblastic leukemia (ALL) is the most common mortal cancer in children. Chronic lymphocytic leukemia (CLL) is the most prevalent form of adult leukemia in western countries. B lymphocytes dominate the origin of both diseases. We have found a possible role for Epstein-Barr virus (EBV) in the origins of both ALL and CLL. We applied our strategy (Nat Genet 50:699, 2018) to determine whether the binding of EBV transcription factors (TFs) was concentrate at the 84 known risk loci for CLL and 16 loci for ALL. We evaluated 52 virally encoded TF ChIP-seq (chromatin immunoprecipitation with DNA sequencing) datasets and complemented this analysis with the results from 1535 human TF ChIP-seq datasets. We found that Epstein-Barr nuclear antigen leader protein (EBNALP), EBNA3C and EBNA2 were concentrated in the CLL loci by 3.7, 3.7 and 3.5-fold with p=4.89*10−19, p=2.74*10−11 and p=1.07*10−8, respectively. The viral (n=3) and human TFs (n=40) cluster together in an optimal subset of ~15 of the 84 known loci in CLL at p<10−6. Eighty percent of the most highly associated viral and human TF ChIP-seq datasets were collected from EBV transformed B cell lines in the Latency III program of viral expression, for which EBNALP, EBNA3C and EBNA2 are viral gene products. In ALL 3 of 16 loci were occupied by EBNA3ABC by a 17.7-fold enrichment with p=2.0*10−10. Among human TFs only c-MYC reached statistical significance binding 2 non-overlapping risk loci with 28.4-fold enrichment with p=2.9*10−10. These results nominate EBV for a role in the pathogenesis of CLL and ALL by a mechanism operating in transformed B cells through the EBV Latency III program of viral expression.
Abstract Innate lymphoid cells (ILCs) are rare tissue-resident “helper” lymphocytes that do not express diversified antigen receptors. Type 3 ILCs (ILC3s) are an important class of these cells enriched in the respiratory and intestinal mucosa, where they regulate inflammation and mucosal homeostasis. To gain insight into the cis-regulatory circuitries underlying ILC3 function, we used high-resolution Capture Hi-C to profile promoter-anchored chromosomal contacts in human primary ILC3s. Combining significant interaction detection with the Activity-By-Contact approach adapted to Capture Hi-C, we reveal a multitude of contacts between promoters and distal regulatory elements and obtain evidence for distinct regulatory wiring of alternative promoters. We find that promoter-interacting regions in ILC3s are enriched for genetic variants associated with multiple immune diseases. Focusing on Crohn’s disease (CD), in which ILC3s are established mediators, we devised a Bayesian approach that incorporates multivariate fine-mapping to link CD-associated genetic variants with putative target genes. We identify known and previously unimplicated genes in conferring genetic risk of CD through activity in ILC3s. This includes the C LN3 gene that is mutated in most cases of the neurodegenerative disorder Batten disease. Using Cln3 mutant mice, we show that CLN3 is a putative negative regulator of IL-17 production in an inflammatory subset of ILC3s. This finding suggests a functional role for CLN3 in ILC3 biology, with mechanistic implications for Crohn’s and Batten diseases.
The etiology of steroid-sensitive nephrotic syndrome (SSNS) is not well understood. Genetic studies have established common single nucleotide polymorphisms (SNPs) that are associated with increased SSNS disease risk. We review previous genetic association studies of SSNS and nominate particular transcriptional regulators and immune cells as potential key players in the etiology of this disease.A list of SNPs associated with SSNS was compiled from published genome wide association and candidate gene studies. The Regulatory Element Locus Intersection (RELI) tool was used to calculate the enrichment of the overlap between disease risk SNPs and the genomic coordinates of data from a collection of >10,000 chromatin immunoprecipitation sequencing experiments.After linkage disequilibrium expansion of the previously reported tag associated SNPs, we identified 192 genetic variants at 8 independent risk loci. Using the Regulatory Element Locus Intersection algorithm, we identified transcriptional regulators with enriched binding at SSNS risk loci (10-05 < Pcorrected < 10-124), including ZNF530, CIITA, CD74, RFX5, and ZNF425. Many of these regulators have well-described roles in the immune response. RNA polymerase II binding in B cells also demonstrated enriched binding at SSNS risk loci (10-37
Abstract Objective IBD is a chronic inflammatory disorder of the GI tract with complex etiology that involves both genetic variants and environmental factors. Several viruses and bacteria have been suggested as potential causes of the disease. The clinical course of IBD appears to be influenced by EBV tissue infection. With a prevalence of >90% in the adult human population, EBV infection is nearly ubiquitous. We reported that the EBV transcription co-factor EBNA2 was concentrated in the risk loci of IBD relative to the remainder of the genome (Pc=1.24E-13, (RR) = 3.33, Nature Genetics 50:699, 2018). The objective of this study was to determine whether new data and studies made available since the data used for our 2018 study further supported the possibility that EBV was related to some IBD as a potential etiologic agent. Methods Genomic approaches offer previously unavailable perspectives for understanding disease mechanisms. The previous analysis was based upon the 9 viral transcription factor and co-factor (TF) datasets available in human cells in 2015.Of the 52 now available there are 29 from EBV. The number of established IBD loci qualifying for analysis has expanded from 112 in 2015 to 324 available now. We applied simulation analysis to assess statistical significance of TF associations with IBD risk loci using RELI (Regulatory Element Locus Intersection) (Nat Genet 50:699, 2018). RELI tests the probability of the observed association by the count of intersections of the same variant from two sources. We extracted 700 SNPs in European ancestry from 93 published GWAS and 350 additional candidate gene studies for IBD. All variants with linkage disequilibrium r2>0.8 of the best variant were included, while loci with r2>0.2 with a more highly significant locus were excluded. Results Among the 52 viral TF datasets, 139 of the 324 IBD loci were occupied by EBNALP (RR=2.08793, Pc=3.7321E-23), 113 by EBNA3C (RR=2.35113, Pc=2.84678E-22), 88 loci by EBNA2 (RR=3.2, Pc=7.8E-32) and 28 loci by EBNA3ABC, which used an antibody that did not distinguish between EBNA3 subtypes. Analysis of 11,535 human TF ChIP-seq datasets produced 1578 with strong associations (p<E-6) with IBD risk loci. The locus intersections of the three EBV gene products optimally clustered at 150 IBD loci (46%) with a set of 26 human TFs (p<2E-308), thereby nominating these loci as potentially having EBV-dependent mechanisms that alter genetic risk for IBD. The human TFs were enriched for (ChIP-seq) data obtained from EBV-infected B lymphocyte cell line (OR = 10.12425329, Chi-squared = 241.2925605) Moreover, the human TFs that participate in the super-enhancers that form in the latency III EBV infected and transformed B cell are also concentrated among the statistically significant associations (p<E-6) with the IBD risk loci (OR) = 2.9448, P = 10E-7). Conclusion These preliminary findings nominate EBV for a role in the pathogenesis and etiology of IBD by mechanisms operating in transformed B cells through the latency III EBV program of viral gene expression.
