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.

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 .

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.

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.

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.