Identification of MEF2B, EBF1, and IL6R as Direct Gene Targets of Epstein-Barr Virus (EBV) Nuclear Antigen 1 Critical for EBV-Infected B-Lymphocyte Survival.

UNLABELLED: Epstein-Barr virus (EBV) nuclear antigen 1 (EBNA1) is the EBV-encoded nuclear antigen and sequence-specific DNA binding protein required for viral origin binding and episome maintenance during latency. EBNA1 can also bind to numerous sites in the cellular genome and can provide a host cell survival function, but it is not yet known how EBNA1 sequence-specific binding is responsible for host cell survival. Here, we integrate EBNA1 chromatin immunoprecipitation sequencing (ChIP-Seq) with transcriptome sequencing (RNA-Seq) after EBNA1 depletion to identify cellular genes directly regulated by EBNA1 that are also essential for B-cell survival. We first compared EBNA1 ChIP-Seq patterns in four different EBV-positive cell types, including Burkitt lymphoma (BL) cells, nasopharyngeal carcinoma (NPC) cells, and lymphoblastoid cell lines (LCLs). EBNA1 binds to ~1,000 sites that are mostly invariant among cell types and share a consensus recognition motif. We found that a large subset of EBNA1 binding sites are located proximal to transcription start sites and correlate genome-wide with transcription activity. EBNA1 bound to genes of high significance for B-cell growth and function, including MEF2B, IL6R, and EBF1. EBNA1 depletion from latently infected LCLs results in the loss of cell proliferation and the loss of gene expression for some EBNA1-bound genes, including MEF2B, EBF1, and IL6R. Depletion of MEF2B, EBF1, or IL6R partially phenocopies EBNA1 depletion by decreasing the cell growth and viability of cells latently infected with EBV. These findings suggest that EBNA1 binds to a large cohort of cellular genes important for cell viability and implicates EBNA1 as a critical regulator of transcription of host cell genes important for enhanced survival of latently infected cells. IMPORTANCE: Epstein-Barr virus (EBV) latent infection is responsible for a variety of lymphoid and epithelial cell malignancies. EBNA1 is the EBV-encoded nuclear antigen that is consistently expressed in all EBV-associated cancers. EBNA1 is known to provide a host cell survival function, but the mechanism is not known. EBNA1 is a sequence-specific binding protein important for viral genome maintenance during latency. Here, by integrating ChIP-Seq and RNA-Seq, we demonstrate that EBNA1 binds directly to the promoter regulatory regions and upregulates the transcription of host genes that are important for the survival of EBV-infected cells. Identification of EBNA1 target genes provides potential new targets for therapeutic intervention in EBV-associated disease.

A role for CTCF and cohesin in subtelomere chromatin organization, TERRA transcription, and telomere end protection.

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.

Cohesins repress Kaposi's sarcoma-associated herpesvirus immediate early gene transcription during latency.

Chromatin-organizing factors such as CTCF and cohesins have been implicated in the control of complex viral regulatory programs. We investigated the role of CTCF and cohesins in the control of the switch from latency to the lytic cycle for Kaposi's sarcoma-associated herpesvirus (KSHV). We found that cohesin subunits but not CTCF are required for the repression of KSHV immediate early gene transcription. Depletion of the cohesin subunits Rad21, SMC1, and SMC3 resulted in lytic cycle gene transcription and viral DNA replication. In contrast, depletion of CTCF failed to induce lytic transcription or DNA replication. Chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq) revealed that cohesins and CTCF bound to several sites within the immediate early control region for ORF50 and to more distal 5' sites that also regulate the divergently transcribed ORF45-ORF46-ORF47 gene cluster. Rad21 depletion led to a robust increase in ORF45, ORF46, ORF47, and ORF50 transcripts, with similar kinetics to that observed with chemical induction by sodium butyrate. During latency, the chromatin between the ORF45 and ORF50 transcription start sites was enriched in histone H3K4me3, with elevated H3K9ac at the ORF45 promoter and elevated H3K27me3 at the ORF50 promoter. A paused form of RNA polymerase II (Pol II) was loosely associated with the ORF45 promoter region during latency but was converted to an active elongating form upon reactivation induced by Rad21 depletion. Butyrate treatment caused a rapid dissociation of cohesins and loss of CTCF binding at the immediate early gene locus, suggesting that cohesins may be a direct target of butyrate-mediated lytic induction. Our findings implicate cohesins as a major repressor of KSHV lytic gene activation and show that they function coordinately with CTCF to regulate the switch between latent and lytic gene activity.

Identification of host-chromosome binding sites and candidate gene targets for Kaposi's sarcoma-associated herpesvirus LANA.

LANA is essential for tethering the Kaposi's sarcoma-associated herpesvirus (KSHV) genome to metaphase chromosomes and for modulating host-cell gene expression, but the binding sites in the host-chromosome remain unknown. Here, we use LANA-specific chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) to identify LANA binding sites in the viral and host-cell genomes of a latently infected pleural effusion lymphoma cell line BCBL1. LANA bound with high occupancy to the KSHV genome terminal repeats (TR) and to a few minor binding sites in the KSHV genome, including the LANA promoter region. We identified 256 putative LANA binding site peaks with P < 0.01 and overlap in two independent ChIP-Seq experiments. We validated several of the high-occupancy binding sites by conventional ChIP assays and quantitative PCR. Candidate cellular LANA binding motifs were identified and assayed for binding to purified recombinant LANA protein in vitro but bound with low affinity compared to the viral TR binding site. More than half of the LANA binding sites (170/256) could be mapped to within 2.5 kb of a cellular gene transcript. Pathways and Gene Ontogeny (GO) analysis revealed that LANA binds to genes within the p53 and tumor necrosis factor (TNF) regulatory network. Further analysis revealed partial overlap of LANA and STAT1 binding sites in several gamma interferon (IFN-γ)-regulated genes. We show that ectopic expression of LANA can downmodulate IFN-γ-mediated activation of a subset of genes, including the TAP1 peptide transporter and proteasome subunit beta type 9 (PSMB9), both of which are required for class I antigen presentation. Our data provide a potential mechanism through which LANA may regulate several host cell pathways by direct binding to gene regulatory elements.

Annotation of gene promoters by integrative data-mining of ChIP-seq Pol-II enrichment data.

BACKGROUND: Use of alternative gene promoters that drive widespread cell-type, tissue-type or developmental gene regulation in mammalian genomes is a common phenomenon. Chromatin immunoprecipitation methods coupled with DNA microarray (ChIP-chip) or massive parallel sequencing (ChIP-seq) are enabling genome-wide identification of active promoters in different cellular conditions using antibodies against Pol-II. However, these methods produce enrichment not only near the gene promoters but also inside the genes and other genomic regions due to the non-specificity of the antibodies used in ChIP. Further, the use of these methods is limited by their high cost and strong dependence on cellular type and context. METHODS: We trained and tested different state-of-art ensemble and meta classification methods for identification of Pol-II enriched promoter and Pol-II enriched non-promoter sequences, each of length 500 bp. The classification models were trained and tested on a bench-mark dataset, using a set of 39 different feature variables that are based on chromatin modification signatures and various DNA sequence features. The best performing model was applied on seven published ChIP-seq Pol-II datasets to provide genome wide annotation of mouse gene promoters. RESULTS: We present a novel algorithm based on supervised learning methods to discriminate promoter associated Pol-II enrichment from enrichment elsewhere in the genome in ChIP-chip/seq profiles. We accumulated a dataset of 11,773 promoter and 46,167 non-promoter sequences, each of length 500 bp, generated from RNA Pol-II ChIP-seq data of five tissues (Brain, Kidney, Liver, Lung and Spleen). We evaluated the classification models in building the best predictor and found that Bagging and Random Forest based approaches give the best accuracy. We implemented the algorithm on seven different published ChIP-seq datasets to provide a comprehensive set of promoter annotations for both protein-coding and non-coding genes in the mouse genome. The resulting annotations contain 13,413 (4,747) protein-coding (non-coding) genes with single promoters and 9,929 (1,858) protein-coding (non-coding) genes with two or more alternative promoters, and a significant number of unassigned novel promoters. CONCLUSION: Our new algorithm can successfully predict the promoters from the genome wide profile of Pol-II bound regions. In addition, our algorithm performs significantly better than existing promoter prediction methods and can be applied for genome-wide predictions of Pol-II promoters.

Genome-wide analysis of host-chromosome binding sites for Epstein-Barr Virus Nuclear Antigen 1 (EBNA1).

The Epstein-Barr Virus (EBV) Nuclear Antigen 1 (EBNA1) protein is required for the establishment of EBV latent infection in proliferating B-lymphocytes. EBNA1 is a multifunctional DNA-binding protein that stimulates DNA replication at the viral origin of plasmid replication (OriP), regulates transcription of viral and cellular genes, and tethers the viral episome to the cellular chromosome. EBNA1 also provides a survival function to B-lymphocytes, potentially through its ability to alter cellular gene expression. To better understand these various functions of EBNA1, we performed a genome-wide analysis of the viral and cellular DNA sites associated with EBNA1 protein in a latently infected Burkitt lymphoma B-cell line. Chromatin-immunoprecipitation (ChIP) combined with massively parallel deep-sequencing (ChIP-Seq) was used to identify cellular sites bound by EBNA1. Sites identified by ChIP-Seq were validated by conventional real-time PCR, and ChIP-Seq provided quantitative, high-resolution detection of the known EBNA1 binding sites on the EBV genome at OriP and Qp. We identified at least one cluster of unusually high-affinity EBNA1 binding sites on chromosome 11, between the divergent FAM55 D and FAM55B genes. A consensus for all cellular EBNA1 binding sites is distinct from those derived from the known viral binding sites, suggesting that some of these sites are indirectly bound by EBNA1. EBNA1 also bound close to the transcriptional start sites of a large number of cellular genes, including HDAC3, CDC7, and MAP3K1, which we show are positively regulated by EBNA1. EBNA1 binding sites were enriched in some repetitive elements, especially LINE 1 retrotransposons, and had weak correlations with histone modifications and ORC binding. We conclude that EBNA1 can interact with a large number of cellular genes and chromosomal loci in latently infected cells, but that these sites are likely to represent a complex ensemble of direct and indirect EBNA1 binding sites.