Landscape of Epstein-Barr virus gene expression and perturbations in cancer

Epstein-Barr virus (EBV) is the causative agent for multiple neoplastic diseases of epithelial and lymphocytic origin1-3. The heterogeneity of the viral elements expressed and the mechanisms by which these coding and non-coding genes maintain cancer cell properties in vivo remain elusive4,5. Here we conducted a multi-modal transcriptomic analysis of EBV-associated neoplasms and identified that the ubiquitously expressed RPMS1 non-coding RNAs support cancer cell properties by disruption of the interferon response. Our map of EBV expression shows a variable, but pervasive expression of BNLF2 discerned from the overlapping LMP1 RNA in bulk sequencing data. Using long-read single-molecule sequencing, we identified three new viral elements within the RPMS1 gene. Furthermore, single-cell sequencing datasets allowed for the separation of cancer cells and healthy cells from the same tissue biopsy and the characterization of a microenvironment containing interferon gamma excreted by EBV-stimulated T-lymphocytes. In comparison with healthy epithelium, EBV-transformed cancer cells exhibited increased proliferation and inhibited immune response induced by the RPMS1-encoded microRNAs. Our atlas of EBV expression shows that the EBV-transformed cancer cells express high levels of non-coding RNAs originating from RPMS1 and that the oncogenic properties are maintained by RPMS1 microRNAs. Through bioinformatic disentanglement of single cells from cancer tissues we identified a positive feedback loop where EBV-activated immune cells stimulate cancer cells to proliferate, which in turn undergo viral reactivation and trigger an immune response.

Serologic Profiling Using an Epstein-Barr Virus Mammalian Expression Library Identifies EBNA1 IgA as a Prediagnostic Marker for Nasopharyngeal Carcinoma

PURPOSE

The favorable prognosis of stage I and II nasopharyngeal carcinoma (NPC) has motivated a search for biomarkers for the early detection and risk assessment of Epstein-Barr virus (EBV)-associated NPC. Although EBV seropositivity is ubiquitous among adults, a spike in antibodies against select EBV proteins is a harbinger of NPC. A serologic survey would likely reveal which EBV antibodies could discriminate those at risk of developing NPC.

EXPERIMENTAL DESIGN

Lysates from a new EBV mammalian expression library were used in a denaturing multiplex immunoblot assay to survey antibodies against EBV in sera collected from healthy individuals who later developed NPC (incident cases) in a prospective cohort from Singapore and validated in an independent cohort from Shanghai, P.R. China.

RESULTS

We show that IgA against EBV nuclear antigen 1 (EBNA1) discriminated incident NPC cases from matched controls with 100% sensitivity and 100% specificity up to 4 years before diagnosis in both Singapore and Shanghai cohorts. Incident NPC cases had a greater IgG repertoire against lytic-classified EBV proteins, and the assortment of IgA against EBV proteins detected by the immunoblot assay increased closer to diagnosis.

CONCLUSIONS

Although NPC tumors consistently harbor latent EBV, the observed heightened systemic and mucosal immunity against lytic-classified antigens years prior to clinical diagnosis is consistent with enhanced lytic transcription. We conclude that an expanding EBV mucosal reservoir (which can be latent and/or lytic) is a risk factor for NPC. This presents an opportunity to identify those at risk of developing NPC using IgA against EBNA1 as a biomarker.

FLAME: long-read bioinformatics tool for comprehensive spliceome characterization

Comprehensive characterization of differentially spliced RNA transcripts with nanopore sequencing is limited by bioinformatics tools that are reliant on existing annotations. We have developed FLAME, a bioinformatics pipeline for alternative splicing analysis of gene-specific or transcriptome-wide long-read sequencing data. FLAME is a Python-based tool aimed at providing comprehensible quantification of full-length splice variants, reliable de novo recognition of splice sites and exons, and representation of consecutive exon connectivity in the form of a weighted adjacency matrix. Notably, this workflow circumvents issues related to inadequate reference annotations and allows for incorporation of short-read sequencing data to improve the confidence of nanopore sequencing reads. In this study, the Epstein-Barr virus long noncoding RNA RPMS1 was used to demonstrate the utility of the pipeline. RPMS1 is ubiquitously expressed in Epstein-Barr virus associated cancer and known to undergo ample differential splicing. To fully resolve the RPMS1 spliceome, we combined gene-specific nanopore sequencing reads from a primary gastric adenocarcinoma and a nasopharyngeal carcinoma cell line with matched publicly available short-read sequencing data sets. All previously reported splice variants, including putative ORFs, were detected using FLAME. In addition, 32 novel exons, including two intron retentions and a cassette exon, were discovered within the RPMS1 gene.

A primary nasopharyngeal three-dimensional air-liquid interface cell culture model of the pseudostratified epithelium reveals differential donor- and cell type-specific susceptibility to Epstein-Barr virus infection

Epstein-Barr virus (EBV) is a ubiquitous γ-herpesvirus with latent and lytic cycles. EBV replicates in the stratified epithelium but the nasopharynx is also composed of pseudostratified epithelium with distinct cell types. Latent infection is associated with nasopharyngeal carcinoma (NPC). Here, we show with nasopharyngeal conditionally reprogrammed cells cultured at the air-liquid interface that pseudostratified epithelial cells are susceptible to EBV infection. Donors varied in susceptibility to de novo EBV infection, but susceptible cultures also displayed differences with respect to pathogenesis. The cultures from one donor yielded lytic infection but cells from two other donors were positive for EBV-encoded EBERs and negative for other lytic infection markers. All cultures stained positive for the pseudostratified markers CK7, MUC5AC, α-tubulin in cilia, and the EBV epithelial cell receptor Ephrin receptor A2. To define EBV transcriptional programs by cell type and to elucidate latent/lytic infection-differential changes, we performed single cell RNA-sequencing on one EBV-infected culture that resulted in alignment with many EBV transcripts. EBV transcripts represented a small portion of the total transcriptome (~0.17%). All cell types in the pseudostratified epithelium had detectable EBV transcripts with suprabasal cells showing the highest number of reads aligning to many EBV genes. Several restriction factors (IRF1, MX1, STAT1, C18orf25) known to limit lytic infection were expressed at lower levels in the lytic subcluster. A third of the differentially-expressed genes in NPC tumors compared to an uninfected pseudostratified ALI culture overlapped with the differentially-expressed genes in the latent subcluster. A third of these commonly perturbed genes were specific to EBV infection and changed in the same direction. Collectively, these findings suggest that the pseudostratified epithelium could harbor EBV infection and that the pseudostratified infection model mirrors many of the transcriptional changes imposed by EBV infection in NPC.

It has been known for over 50 years that EBV infection is associated with NPC. Despite many advances from studies in 2-dimensional cell culture, many aspects of EBV molecular pathogenesis in the nasopharynx remain undefined because the cell types and the biology of the nasopharyngeal epithelium can only be faithfully captured in 3-dimensional cell culture. In the stratified epithelium, cellular differentiation triggers lytic infection but it is not clear to what degree the pseudostratified epithelium is involved. The pseudostratified epithelium is abundant in the lateral wall where the lymphoid-rich fossa of Rosenmüller is located and is a site where NPC tumors most often arises. While the oral epithelium is a site of EBV replication, whether the nasopharyngeal epithelium is a major source of EBV shedding in the nasopharynx is not well defined. Here, we present a 3-dimensional organoid model of the nasopharyngeal pseudostratified epithelium showing that such cells can be infected with EBV in some donor cultures, with examples of both latent and lytic infection. We propose that the cell types of the pseudostratified epithelium should be considered a component of EBV pathogenesis in the nasopharynx and that the difference in donor susceptibility and latent/lytic infection could influence EBV’s fitness in the nasopharynx.