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.

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.

An EBNA3A-Mutated Epstein-Barr Virus Retains the Capacity for Lymphomagenesis in a Cord Blood-Humanized Mouse Model.

Epstein-Barr virus (EBV) causes B cell lymphomas and transforms B cells in vitro The EBV protein EBNA3A collaborates with EBNA3C to repress p16 expression and is required for efficient transformation in vitro An EBNA3A deletion mutant EBV strain was recently reported to establish latency in humanized mice but not cause tumors. Here, we compare the phenotypes of an EBNA3A mutant EBV (Δ3A) and wild-type (WT) EBV in a cord blood-humanized (CBH) mouse model. The hypomorphic Δ3A mutant, in which a stop codon is inserted downstream from the first ATG and the open reading frame is disrupted by a 1-bp insertion, expresses very small amounts of EBNA3A using an alternative ATG at residue 15. Δ3A caused B cell lymphomas at rates similar to their induction by WT EBV but with delayed onset. Δ3A and WT tumors expressed equivalent levels of EBNA2 and p16, but Δ3A tumors in some cases had reduced LMP1. Like the WT EBV tumors, Δ3A lymphomas were oligoclonal/monoclonal, with typically one dominant IGHV gene being expressed. Transcriptome sequencing (RNA-seq) analysis revealed small but consistent gene expression differences involving multiple cellular genes in the WT EBV- versus Δ3A-infected tumors and increased expression of genes associated with T cells, suggesting increased T cell infiltration of tumors. Consistent with an impact of EBNA3A on immune function, we found that the expression of CLEC2D, a receptor that has previously been shown to influence responses of T and NK cells, was markedly diminished in cells infected with EBNA3A mutant virus. Together, these studies suggest that EBNA3A contributes to efficient EBV-induced lymphomagenesis in CBH mice.IMPORTANCE The EBV protein EBNA3A is expressed in latently infected B cells and is important for efficient EBV-induced transformation of B cells in vitro In this study, we used a cord blood-humanized mouse model to compare the phenotypes of an EBNA3A hypomorph mutant virus (Δ3A) and wild-type EBV. The Δ3A virus caused lymphomas with delayed onset compared to the onset of those caused by WT EBV, although the tumors occurred at a similar rate. The WT EBV and EBNA3A mutant tumors expressed similar levels of the EBV protein EBNA2 and cellular protein p16, but in some cases, Δ3A tumors had less LMP1. Our analysis suggested that Δ3A-infected tumors have elevated T cell infiltrates and decreased expression of the CLEC2D receptor, which may point to potential novel roles of EBNA3A in T cell and NK cell responses to EBV-infected tumors.

An EBNA3C-deleted Epstein-Barr virus (EBV) mutant causes B-cell lymphomas with delayed onset in a cord blood-humanized mouse model.

EBV causes human B-cell lymphomas and transforms B cells in vitro. EBNA3C, an EBV protein expressed in latently-infected cells, is required for EBV transformation of B cells in vitro. While EBNA3C undoubtedly plays a key role in allowing EBV to successfully infect B cells, many EBV+ lymphomas do not express this protein, suggesting that cellular mutations and/or signaling pathways may obviate the need for EBNA3C in vivo under certain conditions. EBNA3C collaborates with EBNA3A to repress expression of the CDKN2A-encoded tumor suppressors, p16 and p14, and EBNA3C-deleted EBV transforms B cells containing a p16 germline mutation in vitro. Here we have examined the phenotype of an EBNAC-deleted virus (Δ3C EBV) in a cord blood-humanized mouse model (CBH). We found that the Δ3C virus induced fewer lymphomas (occurring with a delayed onset) in comparison to the wild-type (WT) control virus, although a subset (10/26) of Δ3C-infected CBH mice eventually developed invasive diffuse large B cell lymphomas with type III latency. Both WT and Δ3C viruses induced B-cell lymphomas with restricted B-cell populations and heterogeneous T-cell infiltration. In comparison to WT-infected tumors, Δ3C-infected tumors had greatly increased p16 levels, and RNA-seq analysis revealed a decrease in E2F target gene expression. However, we found that Δ3C-infected tumors expressed c-Myc and cyclin E at similar levels compared to WT-infected tumors, allowing cells to at least partially bypass p16-mediated cell cycle inhibition. The anti-apoptotic proteins, BCL2 and IRF4, were expressed in Δ3C-infected tumors, likely helping cells avoid c-Myc-induced apoptosis. Unexpectedly, Δ3C-infected tumors had increased T-cell infiltration, increased expression of T-cell chemokines (CCL5, CCL20 and CCL22) and enhanced type I interferon response in comparison to WT tumors. Together, these results reveal that EBNA3C contributes to, but is not essential for, EBV-induced lymphomagenesis in CBH mice, and suggest potentially important immunologic roles of EBNA3C in vivo.

A cancer-associated Epstein-Barr virus BZLF1 promoter variant enhances lytic infection.

Latent Epstein-Barr virus (EBV) infection contributes to both B-cell and epithelial-cell malignancies. However, whether lytic EBV infection also contributes to tumors is unclear, although the association between malaria infection and Burkitt lymphomas (BLs) may involve excessive lytic EBV replication. A particular variant of the viral promoter (Zp) that controls lytic EBV reactivation is over-represented, relative to its frequency in non-malignant tissue, in EBV-positive nasopharyngeal carcinomas and AIDS-related lymphomas. To date, no functional differences between the prototype Zp (Zp-P) and the cancer-associated variant (Zp-V3) have been identified. Here we show that a single nucleotide difference between the Zp-V3 and Zp-P promoters creates a binding site for the cellular transcription factor, NFATc1, in the Zp-V3 (but not Zp-P) variant, and greatly enhances Zp activity and lytic viral reactivation in response to NFATc1-inducing stimuli such as B-cell receptor activation and ionomycin. Furthermore, we demonstrate that restoring this NFATc1-motif to the Zp-P variant in the context of the intact EBV B95.8 strain genome greatly enhances lytic viral reactivation in response to the NFATc1-activating agent, ionomycin, and this effect is blocked by the NFAT inhibitory agent, cyclosporine, as well as NFATc1 siRNA. We also show that the Zp-V3 variant is over-represented in EBV-positive BLs and gastric cancers, and in EBV-transformed B-cell lines derived from EBV-infected breast milk of Kenyan mothers that had malaria during pregnancy. These results demonstrate that the Zp-V3 enhances EBV lytic reactivation to physiologically-relevant stimuli, and suggest that increased lytic infection may contribute to the increased prevalence of this variant in EBV-associated malignancies.

CAGE-seq analysis of Epstein-Barr virus lytic gene transcription: 3 kinetic classes from 2 mechanisms.

Epstein-Barr virus (EBV) lytic replication proceeds through an ordered cascade of gene expression that integrates lytic DNA amplification and late gene transcription. We and others previously demonstrated that 6 EBV proteins that have orthologs in ß- and γ-, but not in α-herpesviruses, mediate late gene transcription in a lytic DNA replication-dependent manner. We proposed a model in which the ßγ gene-encoded viral pre-initiation complex (vPIC) mediates transcription from newly replicated viral DNA. While this model explains the dependence of late gene transcription on lytic DNA replication, it does not account for this dependence in α-herpesviruses nor for recent reports that some EBV late genes are transcribed independently of vPIC. To rigorously define which transcription start sites (TSS) are dependent on viral lytic DNA replication or the ßγ complex, we performed Cap Analysis of Gene Expression (CAGE)-seq on cells infected with wildtype EBV or EBV mutants defective for DNA replication, ßγ function, or lacking an origin of lytic replication (OriLyt). This approach identified 16 true-late, 32 early, and 16 TSS that are active at low levels early and are further upregulated in a DNA replication-dependent manner (leaky late). Almost all late gene transcription is vPIC-dependent, with BCRF1 (vIL10), BDLF2, and BDLF3 transcripts being notable exceptions. We present evidence that leaky late transcription is not due to a distinct mechanism, but results from superimposition of the early and late transcription mechanisms at the same promoter. Our results represent the most comprehensive characterization of EBV lytic gene expression kinetics reported to date and suggest that most, but not all EBV late genes are vPIC-dependent.

Epstein-Barr virus ensures B cell survival by uniquely modulating apoptosis at early and late times after infection.

Latent Epstein-Barr virus (EBV) infection is causally linked to several human cancers. EBV expresses viral oncogenes that promote cell growth and inhibit the apoptotic response to uncontrolled proliferation. The EBV oncoprotein LMP1 constitutively activates NFκB and is critical for survival of EBV-immortalized B cells. However, during early infection EBV induces rapid B cell proliferation with low levels of LMP1 and little apoptosis. Therefore, we sought to define the mechanism of survival in the absence of LMP1/NFκB early after infection. We used BH3 profiling to query mitochondrial regulation of apoptosis and defined a transition from uninfected B cells (BCL-2) to early-infected (MCL-1/BCL-2) and immortalized cells (BFL-1). This dynamic change in B cell survival mechanisms is unique to virus-infected cells and relies on regulation of MCL-1 mitochondrial localization and BFL-1 transcription by the viral EBNA3A protein. This study defines a new role for EBNA3A in the suppression of apoptosis with implications for EBV lymphomagenesis.

Differentiation-Dependent LMP1 Expression Is Required for Efficient Lytic Epstein-Barr Virus Reactivation in Epithelial Cells.

Epstein-Barr virus (EBV)-associated diseases of epithelial cells, including tumors that have latent infection, such as nasopharyngeal carcinoma (NPC), and oral hairy leukoplakia (OHL) lesions that have lytic infection, frequently express the viral latent membrane protein 1 (LMP1). In lytically infected cells, LMP1 expression is activated by the BRLF1 (R) immediate early (IE) protein. However, the mechanisms by which LMP1 expression is normally regulated in epithelial cells remain poorly understood, and its potential roles in regulating lytic reactivation in epithelial cells are as yet unexplored. We previously showed that the differentiation-dependent cellular transcription factors KLF4 and BLIMP1 induce lytic EBV reactivation in epithelial cells by synergistically activating the two EBV immediate early promoters (Zp and Rp). Here we show that epithelial cell differentiation also induces LMP1 expression. We demonstrate that KLF4 and BLIMP1 cooperatively induce the expression of LMP1, even in the absence of the EBV IE proteins BZLF1 (Z) and R, via activation of the two LMP1 promoters. Furthermore, we found that differentiation of NOKs-Akata cells by either methylcellulose suspension or organotypic culture induces LMP1 expression prior to Z and R expression. We show that LMP1 enhances the lytic infection-inducing effects of epithelial cell differentiation, as well as 12-O-tetradecanoylphorbol-13-acetate (TPA) and sodium butyrate treatment, in EBV-infected epithelial cells by increasing expression of the Z and R proteins. Our results suggest that differentiation of epithelial cells activates a feed-forward loop in which KLF4 and BLIMP1 first activate LMP1 expression and then cooperate with LMP1 to activate Z and R expression.IMPORTANCE The EBV protein LMP1 is expressed in EBV-associated epithelial cell diseases, regardless of whether these diseases are due to lytic infection (such as oral hairy leukoplakia) or latent infection (such as nasopharyngeal carcinoma). However, surprisingly little is known about how LMP1 expression is regulated in epithelial cells, and there are conflicting reports about whether it plays any role in regulating viral lytic reactivation. In this study, we show that epithelial cell differentiation induces LMP1 expression by increasing expression of two cellular transcription factors (KLF4 and BLIMP1) which cooperatively activate the two LMP1 promoters. We also demonstrate that LMP1 promotes efficient lytic reactivation in EBV-infected epithelial cells by enhancing expression of the Z and R proteins. Thus, in EBV-infected epithelial cells, LMP1 expression is promoted by differentiation and positively regulates lytic viral reactivation.

An Epstein-Barr Virus-Encoded Protein Complex Requires an Origin of Lytic Replication In Cis to Mediate Late Gene Transcription.

Epstein-Barr virus lytic replication is accomplished by an intricate cascade of gene expression that integrates viral DNA replication and structural protein synthesis. Most genes encoding structural proteins exhibit "true" late kinetics-their expression is strictly dependent on lytic DNA replication. Recently, the EBV BcRF1 gene was reported to encode a TATA box binding protein homolog, which preferentially recognizes the TATT sequence found in true late gene promoters. BcRF1 is one of seven EBV genes with homologs found in other ß- and γ-, but not in α-herpesviruses. Using EBV BACmids, we systematically disrupted each of these "ßγ" genes. We found that six of them, including BcRF1, exhibited an identical phenotype: intact viral DNA replication with loss of late gene expression. The proteins encoded by these six genes have been found by other investigators to form a viral protein complex that is essential for activation of TATT-containing reporters in EBV-negative 293 cells. Unexpectedly, in EBV infected 293 cells, we found that TATT reporter activation was weak and non-specific unless an EBV origin of lytic replication (OriLyt) was present in cis. Using two different replication-defective EBV genomes, we demonstrated that OriLyt-mediated DNA replication is required in cis for TATT reporter activation and for late gene expression from the EBV genome. We further demonstrate by fluorescence in situ hybridization that the late BcLF1 mRNA localizes to EBV DNA replication factories. These findings support a model in which EBV true late genes are only transcribed from newly replicated viral genomes.