Decitabine disrupts EBV genomic epiallele DNA methylation patterns around CTCF binding sites to increase chromatin accessibility and lytic transcription in gastric cancer

Epstein-Barr virus (EBV) is associated with 10% of human gastric carcinomas, which are distinguished by a CpG island methylator phenotype. In gastric carcinoma tumors and cell lines, the EBV genome also exhibits a high degree of 5-methyl cytosine (5mC) marks, which are propagated by host DNA methyltransferases (DNMT) with each cell cycle. Therefore, we sought to determine the effect of DNMT inhibition by the small molecule Decitabine (DCB) on EBV genomic 5mC and chromatin structure in two tumor-derived gastric cancer cell lines, YCCEL1 and SNU719. Decitabine effects on EBV genomic 5mC, chromatin structure, and viral gene expression were profiled by reduced representation bisulfite sequencing, ATAC-seq, and RNA-seq, respectively. Decitabine treatment resulted in global viral genome hypomethylation and a global increase in open chromatin. The most striking finding resulted from analyzing the methylation pattern from single RRBS sequencing reads, showing that the EBV genome contains a heterogeneous pool of epigenetic states, each of which is eroded upon Decitabine treatment. We observed heterogeneous 5mC epiallele patterns around EBV genomic CTCF binding sites and lytic gene transcriptional start sites. These results highlight the importance of 5mC in maintaining EBV genomic chromatin structure and latency. Furthermore, the presence of 5mC epialleles suggests EBV+ gastric cancers harbor transcriptionally distinct EBV episomes, which may exert distinct functional roles in maintaining latency and driving tumorigenesis. IMPORTANCE Epstein-Barr virus (EBV) latency is controlled by epigenetic silencing by DNA methylation [5-methyl cytosine (5mC)], histone modifications,modifications, and chromatin looping. However, how they dictate the transcriptional program in EBV-associated gastric cancers remains incompletely understood. EBV-associated gastric cancer displays a 5mC hypermethylated phenotype. A potential treatment for this cancer subtype is the DNA hypomethylating agent, which induces EBV lytic reactivation and targets hypermethylation of the cellular DNA. In this study, we identified a heterogeneous pool of EBV epialleles within two tumor-derived gastric cancer cell lines that are disrupted with a hypomethylating agent. Stochastic DNA methylation patterning at critical regulatory regions may be an underlying mechanism for spontaneous reactivation. Our results highlight the critical role of epigenetic modulation on EBV latency and life cycle, which is maintained through the interaction between 5mC and the host protein CCCTC-binding factor.