Sreeja Parameswaran

Ebola virus disease (EVD) is characterized by immune dysregulation and damaging hyperinflammation. We aimed to characterize the signaling pathways and regulatory mechanisms dysregulated during EVD. To avoid hyperinflammation, innate immune signaling is regulated by post-translational modifications (PTMs), including protein phosphorylation. Here, we show that the protein tyrosine phosphatase nonreceptor type 13 (PTPN13) negatively regulates Interferon (IFN)-β while also positively regulating the neutrophil chemoattractant CXCL1. Using vectors encoding IRF3 with mutations on phosphorylation sites, we identified Y292 on IRF3 as a PTPN13 target of dephosphorylation. Knockout of PTPN13 increased IRF3 phosphorylation and expression of IFNβ and IFN-stimulated genes (ISGs) following poly(I:C) stimulation. Intriguingly, depletion of PTPN13 during Ebola virus (EBOV) infection resulted in decreased IFNβ and ISG induction at later time points post-infection, which correlated with increased viral titers. We identified PTPN13-mediated dephosphorylation of the viral protein VP35 as one potential mechanism inhibiting virus replication. Additionally, the induction of inflammatory chemokines, including CXCL1, decreased in PTPN13 knockout cells late during EBOV infection. These effects could be explained by increased phosphorylation of the regulatory p85 subunit of PI3K. Dephosphorylation of p85 promotes its degradation, subsequently enhancing PI3K kinase activity and downstream signaling via AKT. Together, our study suggests that PTPN13 is involved in immune regulation and efficient antiviral responses by dephosphorylation of IRF3, EBOV-VP35 and PI3K-p85.

Multiple sclerosis (MS) is an immune-mediated demyelinating disease of the central nervous system affecting 2.8 million people worldwide. Both genetic and environmental factors contribute to MS risk, with Epstein-Barr virus (EBV) infection being an important environmental factor. To better clarify the role of EBV in MS, we examined its impact on gene expression, chromatin accessibility, and transcription factor binding in primary B cells and EBV-transformed B cells derived from patients with MS and healthy controls. RNA-seq and ATAC-seq analyses revealed extensive MS-dependent gene expression and chromatin accessibility differences in EBV-transformed, but not in primary B cells. These changes are largely accounted for by the expression levels of EBNA2, an EBV-encoded transcriptional regulator previously implicated in MS. ChIP-seq analysis revealed that EBNA2 binding with its interacting human partners RBPJ, EBF1, and PU.1 is highly enriched at MS genetic risk loci, with extensive EBNA2 allelic binding and increased enrichment at MS genetic risk loci in MS-derived cells. Our findings demonstrate that enhanced EBNA2 activity in MS alters human gene expression, chromatin accessibility, and transcription factor binding in an MS-dependent manner. Collectively, this study provides new insights into the molecular mechanisms through which EBV, particularly EBNA2, interacts with host genetic risk to contribute to MS pathogenesis.

Transcription factors play key roles in cellular biology. Their genomic binding events are enriched at disease- and trait-associated genetic risk loci in particular cellular contexts. To examine this phenomenon in depth, we constructed a PU.1 (SPI1) binding atlas by uniformly processing 260 PU.1 ChIP-seq datasets spanning many immune cell types. Comparison of ChIP-seq peaks across eight major immune cell types identifies shared and cell type-specific PU.1 binding events. DNA sequence analyses reveals context-specific binding shaped by both canonical PU.1 motifs and motifs from partner transcription factor families. Integration of this atlas with genome-wide association studies of blood cell traits and immune-mediated diseases reveals strong enrichment for PU.1 binding events at genetic risk loci and extensive genotype-dependent PU.1 genomic occupancy. We identify cellular contexts in which PU.1 enrichment is most pronounced, including at autoimmune disease loci within EBV-positive B cells. Together, these results define cellular, infectious, and genetic contexts of PU.1 binding that help connect noncoding variation to human phenotypes.

Human cytomegalovirus (HCMV) infects up to 80% of the world’s population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hr after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon 6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV’s established roles in these processes.

Summary: We systematically investigate interactions between Epstein-Barr virus (EBV) transcriptional regulators (vTRs) and the human genome. Starting with 16 known and candidate vTRs, we identify nine whose introduction into human cells results in substantial alterations to host gene expression. Genome-scale determination of vTR genomic binding events and alterations to chromatin accessibility reveals a detailed map of EBV’s functional interactions with the human genome, including >100,000 vTR binding events impacting almost a quarter of all human genes. BMRF1 emerges as a potent regulator, impacting >7,000 genes and altering >37,000 chromatin regions. Our results provide new evidence that EBV RTA interacts with and stabilizes the binding of human RBPJ. Network analysis reveals that many human genes are targeted by multiple EBV vTRs, highlighting the vast coordinated impact of EBV on human gene expression. This study provides a valuable, extensive resource for examining EBV-induced alterations to human gene regulation, with data available on multiple platforms.

There are many well-established relationships between pathogens and human disease, but far fewer when focusing on non-communicable diseases (NCDs). We leverage data from The UK Biobank and TriNetX to perform a systematic survey across 20 pathogens and 426 diseases, primarily NCDs. To this end, we assess the association between disease status and infection history proxies. We identify 206 pathogen-disease pairs that replicate in both cohorts. We replicate many established relationships, including Helicobacter pylori with several gastroenterological diseases and connections between Epstein-Barr virus with multiple sclerosis and lupus. Overall, our approach identified evidence of association for 15 pathogens and 96 distinct diseases, including a currently controversial link between human cytomegalovirus (CMV) and ulcerative colitis (UC). We validate this connection through two orthogonal analyses, revealing increased CMV gene expression in UC patients and enrichment for UC genetic risk signal near human genes that have altered expression upon CMV infection. Collectively, these results form a foundation for future investigations into mechanistic roles played by pathogens in NCDs. All results are easily accessible on our website, https://tf.cchmc.org/pathogen-disease .

Runt-related transcription factor 1 (Runx1) can act as both an activator and a repressor. Here we show that CRISPR-mediated deletion of Runx1 in an embryonic kidney-derived cell (mK4) results in large-scale genome-wide changes to chromatin accessibility and gene expression. Open chromatin regions near down-regulated loci are enriched for Runx sites, remain bound by Runx2, but lose chromatin accessibility and expression in Runx1 knockout cells. Unexpectedly, regions near upregulated genes are depleted of Runx sites and are instead enriched for Zeb transcription factor binding sites. Re-expressing Zeb2 in Runx1 knockout cells restores suppression. These data confirm that Runx1 activity is uniquely needed to maintain open chromatin at many loci, and demonstrate that genome-scale derepression is an indirect consequence of losing Runx1-dependent Zeb expression.

The interplay between environmental and genetic factors plays a key role in the development of many autoimmune diseases. In particular, the Epstein-Barr virus (EBV) is an established contributor to multiple sclerosis, lupus, and other disorders. Previously, we demonstrated that the EBV nuclear antigen 2 (EBNA2) transactivating protein occupies up to half of the risk loci for a set of seven autoimmune disorders. To further examine the mechanistic roles played by EBNA2 at these loci on a genome-wide scale, we globally examined gene expression, chromatin accessibility, chromatin looping, and EBNA2 binding, in a B cell line that was 1) uninfected, 2) infected with a strain of EBV lacking EBNA2, or 3) infected with a strain that expresses EBNA2. We identified >400 EBNA2-dependent differentially expressed human genes and >5,000 EBNA2 binding events in the human genome. ATAC-seq analysis revealed >2,000 regions in the human genome with EBNA2-dependent chromatin accessibility, and HiChIP-seq data revealed >1,700 regions where EBNA2 altered chromatin looping interactions. Importantly, autoimmune genetic risk loci were highly enriched at the sites of these EBNA2-dependent chromatin-altering events. We present examples of autoimmune risk genotype-dependent EBNA2 events, nominating genetic risk mechanisms for autoimmune risk loci such as ZMIZ1 . Taken together, our results reveal important interactions between host genetic variation and EBNA2-driven disease mechanisms. Further, our study highlights a critical role for EBNA2 in rewiring human gene regulatory programs through rearrangement of the chromatin landscape and nominates these interactions as components of genetic mechanisms that influence the risk of multiple autoimmune diseases.

Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3,073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we constructed a massively parallel reporter assay (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into the Epstein-Barr virus-transformed B cell line GM12878 revealed 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Comparison of MPRA results in GM12878 and Jurkat T cell lines highlights shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around allelic variants identifies one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second, larger class of TFs that bind allelically without direct alteration of their motif by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.