Latent Epstein-Barr virus infection collaborates with Myc over-expression in normal human B cells to induce Burkitt-like Lymphomas in mice.

Epstein-Barr virus (EBV) is an important cause of human lymphomas, including Burkitt lymphoma (BL). EBV+ BLs are driven by Myc translocation and have stringent forms of viral latency that do not express either of the two major EBV oncoproteins, EBNA2 (which mimics Notch signaling) and LMP1 (which activates NF-κB signaling). Suppression of Myc-induced apoptosis, often through mutation of the TP53 (p53) gene or inhibition of pro-apoptotic BCL2L11 (BIM) gene expression, is required for development of Myc-driven BLs. EBV+ BLs contain fewer cellular mutations in apoptotic pathways compared to EBV-negative BLs, suggesting that latent EBV infection inhibits Myc-induced apoptosis. Here we use an EBNA2-deleted EBV virus (ΔEBNA2 EBV) to create the first in vivo model for EBV+ BL-like lymphomas derived from primary human B cells. We show that cord blood B cells infected with both ΔEBNA2 EBV and a Myc-expressing vector proliferate indefinitely on a CD40L/IL21 expressing feeder layer in vitro and cause rapid onset EBV+ BL-like tumors in NSG mice. These LMP1/EBNA2-negative Myc-driven lymphomas have wild type p53 and very low BIM, and express numerous germinal center B cell proteins (including TCF3, BACH2, Myb, CD10, CCDN3, and GCSAM) in the absence of BCL6 expression. Myc-induced activation of Myb mediates expression of many of these BL-associated proteins. We demonstrate that Myc blocks LMP1 expression both by inhibiting expression of cellular factors (STAT3 and Src) that activate LMP1 transcription and by increasing expression of proteins (DNMT3B and UHRF1) known to enhance DNA methylation of the LMP1 promoters in human BLs. These results show that latent EBV infection collaborates with Myc over-expression to induce BL-like human B-cell lymphomas in mice. As NF-κB signaling retards the growth of EBV-negative BLs, Myc-mediated repression of LMP1 may be essential for latent EBV infection and Myc translocation to collaboratively induce human BLs.

Epstein-Barr virus LMP1 protein promotes proliferation and inhibits differentiation of epithelial cells via activation of YAP and TAZ.

Latent Epstein-Barr virus (EBV) infection promotes undifferentiated nasopharyngeal carcinomas (NPCs) in humans, but the mechanism(s) for this effect has been difficult to study because EBV cannot transform normal epithelial cells in vitro and the EBV genome is often lost when NPC cells are grown in culture. Here we show that the latent EBV protein, LMP1 (Latent membrane protein 1), induces cellular proliferation and inhibits spontaneous differentiation of telomerase-immortalized normal oral keratinocytes (NOKs) in growth factor-deficient conditions by increasing the activity of the Hippo pathway effectors, YAP (Yes-associated protein) and TAZ (Transcriptional coactivator with PDZ-binding motif). We demonstrate that LMP1 enhances YAP and TAZ activity in NOKs both by decreasing Hippo pathway-mediated serine phosphorylation of YAP and TAZ and increasing Src kinase-mediated Y357 phosphorylation of YAP. Furthermore, knockdown of YAP and TAZ is sufficient to reduce proliferation and promote differentiation in EBV-infected NOKs. We find that YAP and TAZ are also required for LMP1-induced epithelial-to-mesenchymal transition. Importantly, we demonstrate that ibrutinib (an FDA-approved BTK inhibitor that blocks YAP and TAZ activity through an off-target effect) restores spontaneous differentiation and inhibits proliferation of EBV-infected NOKs at clinically relevant doses. These results suggest that LMP1-induced YAP and TAZ activity contributes to the development of NPC.

Human CD4+ iNKT cell adoptive immunotherapy induces anti-tumour responses against CD1d-negative EBV-driven B lymphoma.

Invariant natural killer T (iNKT) cells are a conserved population of innate T lymphocytes that are uniquely suitable as off-the-shelf cellular immunotherapies due to their lack of alloreactivity. Two major subpopulations of human iNKT cells have been delineated, a CD4- subset that has a TH1/cytolytic profile, and a CD4+ subset that appears polyfunctional and can produce both regulatory and immunostimulatory cytokines. Whether these two subsets differ in anti-tumour effects is not known. Using live cell imaging, we found that CD4- iNKT cells limited growth of CD1d+ Epstein-Barr virus (EBV)-infected B-lymphoblastoid spheroids in vitro, whereas CD4+ iNKT cells showed little or no direct anti-tumour activity. However, the effects of the two subsets were reversed when we tested them as adoptive immunotherapies in vivo using a xenograft model of EBV-driven human B cell lymphoma. We found that EBV-infected B cells down-regulated CD1d in vivo, and administering CD4- iNKT cells had no discernable impact on tumour mass. In contrast, xenotransplanted mice bearing lymphomas showed rapid reduction in tumour mass after administering CD4+ iNKT cells. Immunotherapeutic CD4+ iNKT cells trafficked to both spleen and tumour and were associated with subsequently enhanced responses of xenotransplanted human T cells against EBV. CD4+ iNKT cells also had adjuvant-like effects on monocyte-derived DCs and promoted antigen-dependent responses of human T cells in vitro. These results show that allogeneic CD4+ iNKT cellular immunotherapy leads to marked anti-tumour activity through indirect pathways that do not require tumour cell CD1d expression and that are associated with enhanced activity of antigen-specific T cells.

Rta is the principal activator of Epstein-Barr virus epithelial lytic transcription.

The transition from latent Epstein-Barr virus (EBV) infection to lytic viral replication is mediated by the viral transcription factors Rta and Zta. Although both are required for virion production, dissecting the specific roles played by Rta and Zta is challenging because they induce each other's expression. To circumvent this, we constructed an EBV mutant deleted for the genes encoding Rta and Zta (BRLF1 and BZLF1, respectively) in the Akata strain BACmid. This mutant, termed EBVΔRZ, was used to infect several epithelial cell lines, including telomerase-immortalized normal oral keratinocytes, a highly physiologic model of EBV epithelial cell infection. Using RNA-seq, we determined the gene expression induced by each viral transactivator. Surprisingly, Zta alone only induced expression of the lytic origin transcripts BHLF1 and LF3. In contrast, Rta activated the majority of EBV early gene transcripts. As expected, Zta and Rta were both required for expression of late gene transcripts. Zta also cooperated with Rta to enhance a subset of early gene transcripts (Rtasynergy transcripts) that Zta was unable to activate when expressed alone. Interestingly, Rta and Zta each cooperatively enhanced the other's binding to EBV early gene promoters, but this effect was not restricted to promoters where synergy was observed. We demonstrate that Zta did not affect Rtasynergy transcript stability, but increased Rtasynergy gene transcription despite having no effect on their transcription when expressed alone. Our results suggest that, at least in epithelial cells, Rta is the dominant transactivator and that Zta functions primarily to support DNA replication and co-activate a subset of early promoters with Rta. This closely parallels the arrangement in KSHV where ORF50 (Rta homolog) is the principal activator of lytic transcription and K8 (Zta homolog) is required for DNA replication at oriLyt.

Reduced IRF4 expression promotes lytic phenotype in Type 2 EBV-infected B cells.

Humans are infected with two types of EBV (Type 1 (T1) and Type 2 (T2)) that differ substantially in their EBNA2 and EBNA 3A/B/C latency proteins and have different phenotypes in B cells. T1 EBV transforms B cells more efficiently than T2 EBV in vitro, and T2 EBV-infected B cells are more lytic. We previously showed that both increased NFATc1/c2 activity, and an NFAT-binding motif within the BZLF1 immediate-early promoter variant (Zp-V3) contained in all T2 strains, contribute to lytic infection in T2 EBV-infected B cells. Here we compare cellular and viral gene expression in early-passage lymphoblastoid cell lines (LCLs) infected with either T1 or T2 EBV strains. Using bulk RNA-seq, we show that T2 LCLs are readily distinguishable from T1 LCLs, with approximately 600 differentially expressed cellular genes. Gene Set Enrichment Analysis (GSEA) suggests that T2 LCLs have increased B-cell receptor (BCR) signaling, NFAT activation, and enhanced expression of epithelial-mesenchymal-transition-associated genes. T2 LCLs also have decreased RNA and protein expression of a cellular gene required for survival of T1 LCLs, IRF4. In addition to its essential role in plasma cell differentiation, IRF4 decreases BCR signaling. Knock-down of IRF4 in a T1 LCL (infected with the Zp-V3-containing Akata strain) induced lytic reactivation whereas over-expression of IRF4 in Burkitt lymphoma cells inhibited both NFATc1 and NFATc2 expression and lytic EBV reactivation. Single-cell RNA-seq confirmed that T2 LCLs have many more lytic cells compared to T1 LCLs and showed that lytically infected cells have both increased NFATc1, and decreased IRF4, compared to latently infected cells. These studies reveal numerous differences in cellular gene expression in B cells infected with T1 versus T2 EBV and suggest that decreased IRF4 contributes to both the latent and lytic phenotypes in cells with T2 EBV.

Type 1 and Type 2 Epstein-Barr viruses induce proliferation, and inhibit differentiation, in infected telomerase-immortalized normal oral keratinocytes.

Differentiated epithelial cells are an important source of infectious EBV virions in human saliva, and latent Epstein-Barr virus (EBV) infection is strongly associated with the epithelial cell tumor, nasopharyngeal carcinoma (NPC). However, it has been difficult to model how EBV contributes to NPC, since EBV has not been shown to enhance proliferation of epithelial cells in monolayer culture in vitro and is not stably maintained in epithelial cells without antibiotic selection. In addition, although there are two major types of EBV (type 1 (T1) and type 2 (T2)), it is currently unknown whether T1 and T2 EBV behave differently in epithelial cells. Here we inserted a G418 resistance gene into the T2 EBV strain, AG876, allowing us to compare the phenotypes of T1 Akata virus versus T2 AG876 virus in a telomerase-immortalized normal oral keratinocyte cell line (NOKs) using a variety of different methods, including RNA-seq analysis, proliferation assays, immunoblot analyses, and air-liquid interface culture. We show that both T1 Akata virus infection and T2 AG876 virus infection of NOKs induce cellular proliferation, and inhibit spontaneous differentiation, in comparison to the uninfected cells when cells are grown without supplemental growth factors in monolayer culture. T1 EBV and T2 EBV also have a similar ability to induce epithelial-to-mesenchymal (EMT) transition and activate canonical and non-canonical NF-κB signaling in infected NOKs. In contrast to our recent results in EBV-infected lymphoblastoid cells (in which T2 EBV infection is much more lytic than T1 EBV infection), we find that NOKs infected with T1 and T2 EBV respond similarly to lytic inducing agents such as TPA treatment or differentiation. These results suggest that T1 and T2 EBV have similar phenotypes in infected epithelial cells, with both EBV types enhancing cellular proliferation and inhibiting differentiation when growth factors are limiting.

Accurate Quantification of Overlapping Herpesvirus Transcripts from RNA Sequencing Data.

Herpesviruses employ extensive bidirectional transcription of overlapping genes to overcome length constraints on their gene product repertoire. As a consequence, many lytic transcripts cannot be measured individually by reverse transcription-quantitative PCR (RT-qPCR) or conventional RNA sequencing (RNA-seq) analysis. A. G. Bruce, S. Barcy, T. DiMaio, E. Gan, et al. (Pathogens 6:11, 2017, https://doi.org/10.3390/pathogens6010011) have proposed an approximation method using unique coding sequences (UCDS) to estimate lytic gene abundance from Kaposi's sarcoma-associated herpesvirus (KSHV) RNA-seq data. Although UCDS has been widely employed, its accuracy, to our knowledge, has never been rigorously validated for any herpesvirus. In this study, we use cap analysis of gene expression sequencing (CAGE-seq) as a gold-standard to determine the accuracy of UCDS for estimating Epstein-Barr virus (EBV) lytic gene expression levels from RNA-seq data. We also introduce the Unique TranScript (UTS) method, which, like UCDS, estimates transcript abundance from changes in mean RNA-seq read depth. UTS is distinguished by its use of empirically determined 5' and 3' transcript ends rather than coding sequence annotations. Compared to conventional read assignment, both UCDS and UTS improved the accuracy of quantitation of overlapping genes, with UTS giving the most-accurate results. The UTS method discards fewer reads and may be advantageous for experiments with less sequencing depth. UTS is compatible with any aligner and, unlike isoform-aware alignment methods, can be implemented on a laptop computer. Our findings demonstrate that the accuracy achieved by complex and expensive techniques such as CAGE-seq can be approximated using conventional short-read RNA-seq data when read assignment methods address transcript overlap. Although our study focuses on EBV transcription, the UTS method should be applicable across all herpesviruses as well as to other genomes with extensively overlapping transcriptomes. IMPORTANCE Many viruses employ extensively overlapping transcript structures. This complexity makes it difficult to quantify gene expression by using conventional methods, including RNA-seq. Although high-throughput techniques that overcome these limitations exist, they are complex, expensive, and scarce in the herpesvirus literature relative to short-read RNA-seq. Here, using Epstein-Barr virus (EBV) as a model, we demonstrate that conventional RNA-seq analysis methods fail to accurately quantify the abundances of many overlapping transcripts. We further show that the previously described Unique CoDing Sequence (UCDS) method and our Unique TranScript (UTS) method greatly improve the accuracy of EBV lytic gene measurements obtained from RNA-seq data. The UTS method has the advantages of discarding fewer reads and being implementable on a laptop computer. Although this study focuses on EBV, the UCDS and UTS methods should be applicable across herpesviruses and for other viruses that make extensive use of overlapping transcription.

Epstein-Barr virus nuclear antigen 3C (EBNA3C) interacts with the metabolism sensing C-terminal binding protein (CtBP) repressor to upregulate host genes.

Epstein-Barr virus (EBV) infection is associated with the development of specific types of lymphoma and some epithelial cancers. EBV infection of resting B-lymphocytes in vitro drives them to proliferate as lymphoblastoid cell lines (LCLs) and serves as a model for studying EBV lymphomagenesis. EBV nuclear antigen 3C (EBNA3C) is one of the genes required for LCL growth and previous work has suggested that suppression of the CDKN2A encoded tumor suppressor p16INK4A and possibly p14ARF is central to EBNA3C's role in this growth transformation. To directly assess whether loss of p16 and/or p14 was sufficient to explain EBNA3C growth effects, we used CRISPR/Cas9 to disrupt specific CDKN2A exons in EBV transformed LCLs. Disruption of p16 specific exon 1α and the p16/p14 shared exon 2 were each sufficient to restore growth in the absence of EBNA3C. Using EBNA3C conditional LCLs knocked out for either exon 1α or 2, we identified EBNA3C induced and repressed genes. By trans-complementing with EBNA3C mutants, we determined specific genes that require EBNA3C interaction with RBPJ or CtBP for their regulation. Unexpectedly, interaction with the CtBP repressor was required not only for repression, but also for EBNA3C induction of many host genes. Contrary to previously proposed models, we found that EBNA3C does not recruit CtBP to the promoters of these genes. Instead, our results suggest that CtBP is bound to these promoters in the absence of EBNA3C and that EBNA3C interaction with CtBP interferes with the repressive function of CtBP, leading to EBNA3C mediated upregulation.

Hippo signaling effectors YAP and TAZ induce Epstein-Barr Virus (EBV) lytic reactivation through TEADs in epithelial cells.

The Epstein-Barr virus (EBV) human herpesvirus is associated with B-cell and epithelial-cell malignancies, and both the latent and lytic forms of viral infection contribute to the development of EBV-associated tumors. Here we show that the Hippo signaling effectors, YAP and TAZ, promote lytic EBV reactivation in epithelial cells. The transcriptional co-activators YAP/TAZ (which are inhibited by Hippo signaling) interact with DNA-binding proteins, particularly TEADs, to induce transcription. We demonstrate that depletion of either YAP or TAZ inhibits the ability of phorbol ester (TPA) treatment, cellular differentiation or the EBV BRLF1 immediate-early (IE) protein to induce lytic EBV reactivation in oral keratinocytes, and show that over-expression of constitutively active forms of YAP and TAZ reactivate lytic EBV infection in conjunction with TEAD family members. Mechanistically, we find that YAP and TAZ interact with, and activate, the EBV BZLF1 immediate-early promoter. Furthermore, we demonstrate that YAP, TAZ, and TEAD family members are expressed at much higher levels in epithelial cell lines in comparison to B-cell lines, and find that EBV infection of oral keratinocytes increases the level of activated (dephosphorylated) YAP and TAZ. Finally, we have discovered that lysophosphatidic acid (LPA), a known YAP/TAZ activator that plays an important role in inflammation, induces EBV lytic reactivation in epithelial cells through a YAP/TAZ dependent mechanism. Together these results establish that YAP/TAZ are powerful inducers of the lytic form of EBV infection and suggest that the ability of EBV to enter latency in B cells at least partially reflects the extremely low levels of YAP/TAZ and TEADs in this cell type.

ΔNp63α promotes Epstein-Barr virus latency in undifferentiated epithelial cells.

Epstein-Barr virus (EBV) is a human herpesvirus that causes infectious mononucleosis and contributes to both B-cell and epithelial-cell malignancies. EBV-infected epithelial cell tumors, including nasopharyngeal carcinoma (NPC), are largely composed of latently infected cells, but the mechanism(s) maintaining viral latency are poorly understood. Expression of the EBV BZLF1 (Z) and BRLF1 (R) encoded immediate-early (IE) proteins induces lytic infection, and these IE proteins activate each other's promoters. ΔNp63α (a p53 family member) is required for proliferation and survival of basal epithelial cells and is over-expressed in NPC tumors. Here we show that ΔNp63α promotes EBV latency by inhibiting activation of the BZLF1 IE promoter (Zp). Furthermore, we find that another p63 gene splice variant, TAp63α, which is expressed in some Burkitt and diffuse large B cell lymphomas, also represses EBV lytic reactivation. We demonstrate that ΔNp63α inhibits the Z promoter indirectly by preventing the ability of other transcription factors, including the viral IE R protein and the cellular KLF4 protein, to activate Zp. Mechanistically, we show that ΔNp63α promotes viral latency in undifferentiated epithelial cells both by enhancing expression of a known Zp repressor protein, c-myc, and by decreasing cellular p38 kinase activity. Furthermore, we find that the ability of cis-platinum chemotherapy to degrade ΔNp63α contributes to the lytic-inducing effect of this agent in EBV-infected epithelial cells. Together these findings demonstrate that the loss of ΔNp63α expression, in conjunction with enhanced expression of differentiation-dependent transcription factors such as BLIMP1 and KLF4, induces lytic EBV reactivation during normal epithelial cell differentiation. Conversely, expression of ΔNp63α in undifferentiated nasopharyngeal carcinoma cells and TAp63α in Burkitt lymphoma promotes EBV latency in these malignancies.