The establishment and maintenance of Epstein-Barr Virus (EBV) latent infection requires distinct viral gene expression programs. These gene expression programs, termed latency types, are determined largely by promoter selection, and controlled through the interplay between cell-type specific transcription factors, chromatin structure, and epigenetic modifications. We used a genome-wide chromatin-immunoprecipitation (ChIP) assay to identify epigenetic modifications that correlate with different latency types. We found that the chromatin insulator protein CTCF binds at several key regulatory nodes in the EBV genome and may compartmentalize epigenetic modifications across the viral genome. Highly enriched CTCF binding sites were identified at the promoter regions upstream of Cp, Wp, EBERs, and Qp. Since Qp is essential for long-term maintenance of viral genomes in type I latency and epithelial cell infections, we focused on the role of CTCF in regulating Qp. Purified CTCF bound approximately 40 bp upstream of the EBNA1 binding sites located at +10 bp relative to the transcriptional initiation site at Qp. Mutagenesis of the CTCF binding site in EBV bacmids resulted in a decrease in the recovery of stable hygromycin-resistant episomes in 293 cells. EBV lacking the Qp CTCF site showed a decrease in Qp transcription initiation and a corresponding increase in Cp and Fp promoter utilization at 8 weeks post-transfection. However, by 16 weeks post-transfection, bacmids lacking CTCF sites had no detectable Qp transcription and showed high levels of histone H3 K9 methylation and CpG DNA methylation at the Qp initiation site. These findings provide direct genetic evidence that CTCF functions as a chromatin insulator that prevents the promiscuous transcription of surrounding genes and blocks the epigenetic silencing of an essential promoter, Qp, during EBV latent infection.

Epstein-Barr Virus (EBV) latent infection is associated with several human malignancies and is a causal agent of lymphoproliferative diseases during immunosuppression. While inhibitors of herpesvirus DNA polymerases, like gancyclovir, reduce EBV lytic cycle infection, these treatments have limited efficacy for treating latent infection. EBNA1 is an EBV-encoded DNA-binding protein required for viral genome maintenance during latent infection.

The evolution of tissue-specific general transcription factors (GTFs), such as testis-specific TBP-related factor 2 (TRF2), enables the spatiotemporal expression of highly specialized genetic programs. Taspase1 is a protease that cleaves nuclear factors MLL1, MLL2, TFIIAα-β, and ALFα-β (TFIIAτ). Here, we demonstrate that Taspase1-mediated processing of TFIIAα-β drives mammalian spermatogenesis. Both Taspase1(-/-) and noncleavable TFIIAα-βnc/nc testes release immature germ cells with impaired transcription of Transition proteins (Tnp) and Protamines (Prm), exhibiting chromatin compaction defects and recapitulating those observed with TRF2(-/-) testes. Although the unprocessed TFIIA still complexes with TRF2, this complex is impaired in targeting and thus activating Tnp1 and Prm1 promoters. The current study presents a paradigm in which a protease (Taspase1) cleaves a ubiquitously expressed GTF (TFIIA) to enable tissue-specific (testis) transcription, meeting the demand for sophisticated regulation of distinct subsets of genes in higher organisms.

CTCF (CCCTC-binding factor) is a highly conserved multifunctional DNA-binding protein with thousands of binding sites genome-wide. Our previous work suggested that differences in CTCF's binding site sequence may affect the regulation of CTCF recruitment and its function. To investigate this possibility, we characterized changes in genome-wide CTCF binding and gene expression during differentiation of mouse embryonic stem cells. After separating CTCF sites into three classes (LowOc, MedOc and HighOc) based on similarity to the consensus motif, we found that developmentally regulated CTCF binding occurs preferentially at LowOc sites, which have lower similarity to the consensus. By measuring the affinity of CTCF for selected sites, we show that sites lost during differentiation are enriched in motifs associated with weaker CTCF binding in vitro. Specifically, enrichment for T at the 18(th) position of the CTCF binding site is associated with regulated binding in the LowOc class and can predictably reduce CTCF affinity for binding sites. Finally, by comparing changes in CTCF binding with changes in gene expression during differentiation, we show that LowOc and HighOc sites are associated with distinct regulatory functions. Our results suggest that the regulatory control of CTCF is dependent in part on specific motifs within its binding site.

Maternal genomic imprints are established during oogenesis. Histone deacetylases (HDACs) 1 and 2 are required for oocyte development in mouse, but their role in genomic imprinting is unknown. We find that Hdac1:Hdac2(-/-) double-mutant growing oocytes exhibit global DNA hypomethylation and fail to establish imprinting marks for Igf2r, Peg3, and Srnpn. Global hypomethylation correlates with increased retrotransposon expression and double-strand DNA breaks. Nuclear-associated DNMT3A2 is reduced in double-mutant oocytes, and injecting these oocytes with Hdac2 partially restores DNMT3A2 nuclear staining. DNMT3A2 co-immunoprecipitates with HDAC2 in mouse embryonic stem cells. Partial loss of nuclear DNMT3A2 and HDAC2 occurs in Sin3a(-/-) oocytes, which exhibit decreased DNA methylation of imprinting control regions for Igf2r and Srnpn, but not Peg3. These results suggest seminal roles of HDAC1/2 in establishing maternal genomic imprints and maintaining genomic integrity in oocytes mediated in part through a SIN3A complex that interacts with DNMT3A2.

Exposure to endocrine disruptors is associated with developmental defects. One compound of concern, to which humans are widely exposed, is bisphenol A (BPA). In model organisms, BPA exposure is linked to metabolic disorders, infertility, cancer, and behavior anomalies. Recently, BPA exposure has been linked to DNA methylation changes, indicating that epigenetic mechanisms may be relevant. We investigated effects of exposure on genomic imprinting in the mouse as imprinted genes are regulated by differential DNA methylation and aberrant imprinting disrupts fetal, placental, and postnatal development. Through allele-specific and quantitative real-time PCR analysis, we demonstrated that maternal BPA exposure during late stages of oocyte development and early stages of embryonic development significantly disrupted imprinted gene expression in embryonic day (E) 9.5 and 12.5 embryos and placentas. The affected genes included Snrpn, Ube3a, Igf2, Kcnq1ot1, Cdkn1c, and Ascl2; mutations and aberrant regulation of these genes are associated with imprinting disorders in humans. Furthermore, the majority of affected genes were expressed abnormally in the placenta. DNA methylation studies showed that BPA exposure significantly altered the methylation levels of differentially methylated regions (DMRs) including the Snrpn imprinting control region (ICR) and Igf2 DMR1. Moreover, exposure significantly reduced genome-wide methylation levels in the placenta, but not the embryo. Histological and immunohistochemical examinations revealed that these epigenetic defects were associated with abnormal placental development. In contrast to this early exposure paradigm, exposure outside of the epigenetic reprogramming window did not cause significant imprinting perturbations. Our data suggest that early exposure to common environmental compounds has the potential to disrupt fetal and postnatal health through epigenetic changes in the embryo and abnormal development of the placenta.

The contribution of human subtelomeric DNA and chromatin organization to telomere integrity and chromosome end protection is not yet understood in molecular detail. Here, we show by ChIP-Seq that most human subtelomeres contain a CTCF- and cohesin-binding site within ∼1-2 kb of the TTAGGG repeat tract and adjacent to a CpG-islands implicated in TERRA transcription control. ChIP-Seq also revealed that RNA polymerase II (RNAPII) was enriched at sites adjacent to the CTCF sites and extending towards the telomere repeat tracts. Mutation of CTCF-binding sites in plasmid-borne promoters reduced transcriptional activity in an orientation-dependent manner. Depletion of CTCF by shRNA led to a decrease in TERRA transcription, and a loss of cohesin and RNAPII binding to the subtelomeres. Depletion of either CTCF or cohesin subunit Rad21 caused telomere-induced DNA damage foci (TIF) formation, and destabilized TRF1 and TRF2 binding to the TTAGGG proximal subtelomere DNA. These findings indicate that CTCF and cohesin are integral components of most human subtelomeres, and important for the regulation of TERRA transcription and telomere end protection.

Viruses suppress immune recognition through diverse mechanisms. Epstein-Barr Virus (EBV) establishes latent infection in memory B-lymphocytes and B-cell malignancies where it impacts B-cell immune function. We show here that EBV primary infection of naïve B-cells results in a robust down-regulation of HLA genes. We found that the viral encoded transcriptional regulatory factor EBNA2 bound to multiple regulatory regions in the HLA locus. Conditional expression of EBNA2 correlated with the down regulation of HLA class II transcription. EBNA2 down-regulation of HLA transcription was found to be dependent on CIITA, the major transcriptional activator of HLA class II gene transcription. We identified a major EBNA2 binding site downstream of the CIITA gene and upstream of DEXI, a dexamethasone inducible gene that is oriented head-to-head with CIITA gene transcripts. CRISPR/Cas9 deletion of the EBNA2 site upstream of DEXI attenuated CIITA transcriptional repression. EBNA2 caused an increase in DEXI transcription and a graded change in histone modifications with activation mark H3K27ac near the DEXI locus, and a loss of activation marks at the CIITA locus. A prominent CTCF binding site between CIITA and DEXI enhancers was mutated and further diminished the effects of EBNA2 on CIITA. Analysis of HiC data indicate that DEXI and CIITA enhancers are situated in different chromosome topological associated domains (TADs). These findings suggest that EBNA2 down regulates HLA-II genes through the down regulation of CIITA, and that this down regulation is an indirect consequence of EBNA2 enhancer formation at a neighboring TAD. We propose that enhancer competition between these neighboring chromosome domains represents a novel mechanism for gene regulation demonstrated by EBNA2.

Subpopulations of B-lymphocytes traffic to different sites and organs to provide diverse and tissue-specific functions. Here, we provide evidence that epigenetic differences confer a neuroinvasive phenotype. An EBV+ B cell lymphoma cell line (M14) with low frequency trafficking to the CNS was neuroadapted to generate a highly neuroinvasive B-cell population (MUN14). MUN14 B cells efficiently infiltrated the CNS within one week and produced neurological pathologies. We compared the gene expression profiles of viral and cellular genes using RNA-Seq and identified one viral (EBNA1) and several cellular gene candidates, including secreted phosphoprotein 1/osteopontin (SPP1/OPN), neuron navigator 3 (NAV3), CXCR4, and germinal center-associated signaling and motility protein (GCSAM) that were selectively upregulated in MUN14. ATAC-Seq and ChIP-qPCR revealed that these gene expression changes correlated with epigenetic changes at gene regulatory elements. The neuroinvasive phenotype could be attenuated with a neutralizing antibody to OPN, confirming the functional role of this protein in trafficking EBV+ B cells to the CNS. These studies indicate that B-cell trafficking to the CNS can be acquired by epigenetic adaptations and provide a new model to study B-cell neuroinvasion associated CNS lymphoma and autoimmune disease of the CNS, including multiple sclerosis (MS).

Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiologic agent of Kaposi's sarcoma (KS). It is endemic in a number of sub-Saharan African countries with infection rate of >50%. The high prevalence of HIV-1 coupled with late presentation of advanced cancer staging make KS the leading cancer in the region with poor prognosis and high mortality. Disease markers and cellular functions associated with KS tumorigenesis remain ill-defined. Several studies have attempted to investigate changes of the gene profile with in vitro infection of monoculture models, which are not likely to reflect the cellular complexity of the in vivo lesion environment. Our approach is to characterize and compare the gene expression profile in KS lesions versus non-cancer tissues from the same individual. Such comparisons could identify pathways critical for KS formation and maintenance. This is the first study that utilized high throughput RNA-seq to characterize the viral and cellular transcriptome in tumor and non-cancer biopsies of African epidemic KS patients. These patients were treated anti-retroviral therapy with undetectable HIV-1 plasma viral load. We found remarkable variability in the viral transcriptome among these patients, with viral latency and immune modulation genes most abundantly expressed. The presence of KSHV also significantly affected the cellular transcriptome profile. Specifically, genes involved in lipid and glucose metabolism disorder pathways were substantially affected. Moreover, infiltration of immune cells into the tumor did not prevent KS formation, suggesting some functional deficits of these cells. Lastly, we found only minimal overlaps between our in vivo cellular transcriptome dataset with those from in vitro studies, reflecting the limitation of in vitro models in representing tumor lesions. These findings could lead to the identification of diagnostic and therapeutic markers for KS, and will provide bases for further mechanistic studies on the functions of both viral and cellular genes that are involved.

The molecular basis for the formation of functional, higher-ordered macro-molecular domains is not completely known. The Kaposi's Sarcoma-Associated Herpesvirus (KSHV) genome forms a super-molecular domain structure during latent infection that is strictly dependent on the DNA binding of the viral nuclear antigen LANA to the viral terminal repeats (TR). LANA is known to form oligomeric structures that have been implicated in viral episome maintenance. In this study, we show that the LANA oligomerization interface is required for the formation of higher-order nuclear bodies that partially colocalize with DAXX, EZH2, H3K27me3, and ORC2 but not with PML. These nuclear bodies assemble at the periphery of condensed cellular chromosomes during mitotic cell division. We demonstrate that the LANA oligomerization interface contributes to the cooperative DNA binding at the viral TR and the recruitment of ORC to the viral episome. Oligomerization mutants failed to auto-regulate LANA/ORF73 transcription, and this correlated with the loss of a chromosome conformational DNA-loop between the TR and LANA promoter. Viral genomes with LANA oligomerization mutants were subject to genome rearrangements including the loss of subgenomic DNA. Our data suggests that LANA oligomerization drives stable binding to the TR and formation of an epigenetically stable chromatin architecture resulting in higher-order LANA nuclear bodies important for viral genome integrity and long-term episome persistence.

Epstein-Barr Virus (EBV) latent infection is a causative co-factor for endemic Nasopharyngeal Carcinoma (NPC). NPC-associated variants have been identified in EBV-encoded nuclear antigen EBNA1. Here, we solve the X-ray crystal structure of an NPC-derived EBNA1 DNA binding domain (DBD) and show that variant amino acids are found on the surface away from the DNA binding interface. We show that NPC-derived EBNA1 is compromised for DNA replication and episome maintenance functions. Recombinant virus containing the NPC EBNA1 DBD are impaired in their ability to immortalize primary B-lymphocytes and suppress lytic transcription during early stages of B-cell infection. We identify Survivin as a host protein deficiently bound by the NPC variant of EBNA1 and show that Survivin depletion compromises EBV episome maintenance in multiple cell types. We propose that endemic variants of EBNA1 play a significant role in EBV-driven carcinogenesis by altering key regulatory interactions that destabilize latent infection.

Liquid-liquid phase separation (LLPS) can drive formation of diverse and essential macromolecular structures, including those specified by viruses. Kaposi's Sarcoma-Associated Herpesvirus (KSHV) genomes associate with the viral encoded Latency-Associated Nuclear Antigen (LANA) to form stable nuclear bodies (NBs) during latent infection. Here, we show that LANA-NB formation and KSHV genome conformation involves LLPS. Using LLPS disrupting solvents, we show that LANA-NBs are partially disrupted, while DAXX and PML foci are highly resistant. LLPS disruption altered the LANA-dependent KSHV chromosome conformation but did not stimulate lytic reactivation. We found that LANA-NBs undergo major morphological transformation during KSHV lytic reactivation to form LANA-associated replication compartments encompassing KSHV DNA. DAXX colocalizes with the LANA-NBs during latency but is evicted from the LANA-associated lytic replication compartments. These findings indicate the LANA-NBs are dynamic super-molecular nuclear structures that partly depend on LLPS and undergo morphological transitions corresponding to the different modes of viral replication.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that persists as a multicopy episome in proliferating host cells. Episome maintenance is strictly dependent on EBNA1, a sequence-specific DNA-binding protein with no known enzymatic activities. Here, we show that EBNA1 forms a cell cycle-dependent DNA crosslink with the EBV origin of plasmid replication oriP. EBNA1 tyrosine 518 (Y518) is essential for crosslinking to oriP and functionally required for episome maintenance and generation of EBV-transformed lymphoblastoid cell lines (LCLs). Mechanistically, Y518 is required for replication fork termination at oriP in vivo and for formation of SDS-resistant complexes in vitro. EBNA1-DNA crosslinking corresponds to single-strand endonuclease activity specific to DNA structures enriched at replication-termination sites, such as 4-way junctions. These findings reveal that EBNA1 forms tyrosine-dependent DNA-protein crosslinks and single-strand cleavage at oriP required for replication termination and viral episome maintenance.

KSHV-associated cancers have poor prognoses and lack therapeutics that selectively target viral gene functions. We developed a screening campaign to identify known drugs that could be repurposed for the treatment of KSHV-associated cancers. We focused on primary effusion lymphoma (PEL), which has particularly poor treatment outcomes. We developed a luciferase reporter assay to test the ability of drugs to inhibit DNA binding of the KSHV LANA DNA binding domain (DBD). In parallel, we screened drugs for selective inhibition of a KSHV+ PEL cells. While potent hits were identified in each assay, only one hit, Mubritinib, was found to score in both assays. Mubritinib caused PEL cells to undergo cell cycle arrest with accumulation of sub-G1 population and Annexin V. Mubritinib inhibited LANA binding to KSHV terminal repeat (TR) DNA in KSHV+ PEL cells, but did not lead to KSHV lytic cycle reactivation. Mubritinib was originally identified as a receptor tyrosine kinase (RTK) inhibitor selective for HER2/ErbB2. But recent studies have revealed that Mubritinib can also inhibit the electron transport chain (ETC) complex at nanomolar concentrations. We found that other related ETC complex inhibitors (Rotenone and Deguelin) exhibited PEL cell growth inhibition while RTK inhibitors failed. Seahorse analysis demonstrated that Mubritinib selectively inhibits the maximal oxygen consumption (OCR) in PEL cells and metabolomics revealed changes in ATP/ADP and ATP/AMP ratios. These findings indicate that PEL cells are selectively sensitive to ETC complex inhibitors and provide a rationale for repurposing Mubritinib for selective treatment of PEL.

The telomeric shelterin protein telomeric repeat-binding factor 2 (TRF2) recruits origin recognition complex (ORC) proteins, the foundational building blocks of DNA replication origins, to telomeres. We seek to determine whether TRF2-recruited ORC proteins give rise to functional origins in telomere repeat tracts. We find that reduction of telomeric recruitment of ORC2 by expression of an ORC interaction-defective TRF2 mutant significantly reduces telomeric initiation events in human cells. This reduction in initiation events is accompanied by telomere repeat loss, telomere aberrations and dysfunction. We demonstrate that telomeric origins are activated by induced replication stress to provide a key rescue mechanism for completing compromised telomere replication. Importantly, our studies also indicate that the chromatin remodeler SNF2H promotes telomeric initiation events by providing access for ORC2. Collectively, our findings reveal that active recruitment of ORC by TRF2 leads to formation of functional origins, providing an important mechanism for avoiding telomere dysfunction and rescuing challenged telomere replication.

In sub-Saharan Africa, endemic Kaposi's sarcoma (EnKS) is still prevalent despite high incidence of epidemic Kaposi's sarcoma (EpKS) resulting from the on-going HIV-1 epidemic. While KSHV is clearly the etiologic agent of KS, the mechanisms underlying KS development are not fully understood. For example, HIV-1 co-infection and concomitant immune dysfunction have been associated with EpKS development. However, the direct or indirect role(s) of HIV-1, and therefore of immune suppression, in EpKS remains unclear. How, or whether, EpKS is mechanistically distinct from EnKS is unknown. Thus, the absence of HIV-1 co-infection in EnKS provides a unique control for investigating and deciphering whether HIV-1 plays a direct or indirect role in the EpKS tumor microenvironment. We hypothesized that HIV-1 co-infection would induce transcriptome changes that differentiate EpKS from EnKS, thereby defining the direct intra-tumor role of HIV-1 in KS. Comparison of ART-treated and -naïve patients would further define the impact of ART on the KS transcriptome. We utilized RNA-seq followed by multiparameter bioinformatics analysis to compare transcriptomes from KS lesions to uninvolved control skin. We provide the first transcriptomic comparison of EpKS versus EnKS, ART-treated vs-naïve EpKS and male vs female EpKS to define the roles of HIV-1 co-infection, the impact of ART, and gender on KS gene expression profiles. Our findings suggest that ART-use and gender have minimal impact on transcriptome profiles of KS lesions. Gene expression profiles strongly correlated between EpKS and EnKS patients (Spearman r = 0.83, p<10-10). A subset of genes involved in tumorigenesis and inflammation/immune responses showed higher magnitude, but not unique dysregulation in EnKS compared to EpKS. While gender and ART had no detectable contribution, the trend toward higher magnitude of gene dysregulation in EnKS coupled with the absence of HIV-1 transcripts in EpKS may suggest an indirect or systemic effect of HIV-1 to promote KS tumorigenesis.

Kaposi's Sarcoma (KS) is a heterogenous, multifocal vascular malignancy caused by the human herpesvirus 8 (HHV8), also known as Kaposi's Sarcoma-Associated Herpesvirus (KSHV). Here, we show that KS lesions express iNOS/NOS2 broadly throughout KS lesions, with enrichment in LANA positive spindle cells. The iNOS byproduct 3-nitrotyrosine is also enriched in LANA positive tumor cells and colocalizes with a fraction of LANA-nuclear bodies. We show that iNOS is highly expressed in the L1T3/mSLK tumor model of KS. iNOS expression correlated with KSHV lytic cycle gene expression, which was elevated in late-stage tumors (>4 weeks) but to a lesser degree in early stage (1 week) xenografts. Further, we show that L1T3/mSLK tumor growth is sensitive to an inhibitor of nitric oxide, L-NMMA. L-NMMA treatment reduced KSHV gene expression and perturbed cellular gene pathways relating to oxidative phosphorylation and mitochondrial dysfunction. These finding suggest that iNOS is expressed in KSHV infected endothelial-transformed tumor cells in KS, that iNOS expression depends on tumor microenvironment stress conditions, and that iNOS enzymatic activity contributes to KS tumor growth.

DNA methylation erasure is required for mammalian primordial germ cell reprogramming. TET enzymes iteratively oxidize 5-methylcytosine to generate 5-hyroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxycytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during germline reprogramming remains unresolved due to the lack of genetic models that decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 ( Tet1-HxD ) and TET1 that stalls oxidation at 5hmC ( Tet1-V ). Tet1 -/- , Tet1 V/V , and Tet1 HxD/HxD sperm methylomes show that TET1 V and TET1 HxD rescue most Tet1 -/- hypermethylated regions, demonstrating the importance of TET1’s extra-catalytic functions. Imprinted regions, in contrast, require iterative oxidation. We further reveal a broader class of hypermethylated regions in sperm of Tet1 mutant mice that are excluded from de novo methylation during male germline development and depend on TET oxidation for reprogramming. Our study underscores the link between TET1-mediated demethylation during reprogramming and sperm methylome patterning.

About 15% of human cancer cases are attributed to viral infections. To date, virus expression in tumor tissues has been mostly studied by aligning tumor RNA sequencing reads to databases of known viruses. To allow identification of divergent viruses and rapid characterization of the tumor virome, we develop viRNAtrap, an alignment-free pipeline to identify viral reads and assemble viral contigs. We utilize viRNAtrap, which is based on a deep learning model trained to discriminate viral RNAseq reads, to explore viral expression in cancers and apply it to 14 cancer types from The Cancer Genome Atlas (TCGA). Using viRNAtrap, we uncover expression of unexpected and divergent viruses that have not previously been implicated in cancer and disclose human endogenous viruses whose expression is associated with poor overall survival. The viRNAtrap pipeline provides a way forward to study viral infections associated with different clinical conditions.