# Single-cell profiling of HDAC inhibitor-induced EBV lytic heterogeneity defines abortive and refractory states in B lymphoblasts

**Authors:** Lauren E. Haynes, Ashley P. Barry, Micah A. Luftig, John Karijolich, John Karijolich, John Karijolich

PMC · DOI: 10.1371/journal.ppat.1013610 · PLOS Pathogens · 2026-03-23

## TL;DR

The study uses single-cell profiling to explore how HDAC inhibitors affect EBV lytic reactivation in cancer cells, revealing abortive and refractory states.

## Contribution

The study identifies CD137 signaling as a novel factor influencing abortive lytic reactivation in EBV-positive malignancies.

## Key findings

- HDAC inhibitors induce heterogeneous lytic reactivation in EBV-positive cell lines.
- Abortive lytic cells upregulate NF-κB and immune signaling genes like CD137/CD137L.
- CD137 receptor is functionally involved in preventing successful lytic reactivation.

## Abstract

Epstein-Barr virus (EBV) is associated with multiple malignancies including Burkitt lymphoma (BL), Hodgkin’s lymphomas, nasopharyngeal carcinomas (NPC), and gastric cancers. Canonically, EBV positive tumors display latent gene expression programs that are difficult to target pharmacologically. To overcome this hurdle, lytic reactivation therapies have been developed based on HDAC inhibition with limited mechanistic studies. We therefore characterized the impact of pan-HDAC inhibitor, panobinostat, and class I HDAC inhibitor, nanatinostat, on the growth, survival, and lytic reactivation of four EBV-positive cell lines: P3HR1-ZHT BL, Jijoye BL, IBL-1 immunoblastic lymphoma, and de novo infection derived lymphoblastoid cell lines (LCL). All lines were sensitive, enabling us to define ranges of sensitivity within which to use single cell approaches to assess early EBV lytic gene expression, cell cycle state, and apoptosis. We observed that each EBV-positive model of malignancy responded uniquely to the same HDAC inhibitors and that lytic reactivation was successful in only a small percentage of the cell population. To elucidate the potential role of host factors in preventing successful lytic reactivation, we performed single-cell RNA sequencing on the P3HR1-ZHT BL line treated with the HDAC inhibitor panobinostat. We observed that abortive lytic cells, or cells that do not successfully progress through the lytic cycle, upregulated genes downstream of NF-κB activity. Additionally, genes involved in immune signaling including the CD137/CD137L signaling axis, were upregulated in abortive lytic cells. Functional validation through a Cas9-RNP approach revealed that the CD137 receptor is indeed involved in preventing successful lytic reactivation. These data have important implications for how we approach oncolytic therapies for EBV-associated malignancies.

Epstein-Barr virus (EBV) is an extremely prevalent human herpesvirus that is associated with a variety of cancers and autoimmune diseases. EBV establishes latent infection in the host and, under various circumstances, can reactivate the lytic cycle to produce more infectious particles. In the context of EBV-associated malignancies, the virus is most often maintained in a latent state, which makes it difficult to target with pharmaceuticals. To develop more viral targeted strategies, kick and kill regimens have been investigated. This therapy involves reactivating the virus with an HDAC inhibitor followed by treatment with an antiviral drug. It is well established that reactivating EBV with pharmaceuticals is often inefficient and leads to heterogeneous responses, including an abortive lytic trajectory. To better characterize the overall effect of two classes of HDAC inhibitors in various models of EBV-associated malignancies, we utilized single-cell techniques to capture various responses to stimuli. Consistent with prior studies, HDAC inhibition led to both successful and abortive lytic populations. Single-cell RNA sequencing provided evidence of upregulated immune signaling pathways in this abortive lytic population. Downstream functional validation revealed that CD137 signaling is important for this abortive lytic phenotype. This study provides in depth characterization of lytic reactivation with a biologically relevant stimulus.

## Linked entities

- **Genes:** TNFRSF9 (TNF receptor superfamily member 9) [NCBI Gene 3604], NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790]
- **Proteins:** TNFRSF9 (TNF receptor superfamily member 9), TNFSF9 (TNF superfamily member 9)
- **Chemicals:** panobinostat (PubChem CID 6918837), nanatinostat (PubChem CID 49857317)
- **Diseases:** Burkitt lymphoma (MONDO:0007243)

## Full-text entities

- **Genes:** STAT1 (signal transducer and activator of transcription 1) [NCBI Gene 6772] {aka CANDF7, IMD31A, IMD31B, IMD31C, ISGF-3, STAT91}, SP1 (Sp1 transcription factor) [NCBI Gene 6667], ATM (ATM serine/threonine kinase) [NCBI Gene 472] {aka AT1, ATA, ATC, ATD, ATDC, ATE}, CASP1 (caspase 1) [NCBI Gene 834] {aka ICE, IL1BC, P45}, VCAM1 (vascular cell adhesion molecule 1) [NCBI Gene 7412] {aka CD106, INCAM-100}, TNFRSF9 (TNF receptor superfamily member 9) [NCBI Gene 3604] {aka 4-1BB, CD137, CDw137, ILA, IMD109}, IL10 (interleukin 10) [NCBI Gene 3586] {aka CSIF, GVHDS, IL-10, IL10A, TGIF}, TNFRSF1A (TNF receptor superfamily member 1A) [NCBI Gene 7132] {aka CD120a, FPF, TBP1, TNF-R, TNF-R-I, TNF-R55}, HDAC11 (histone deacetylase 11) [NCBI Gene 79885] {aka HD11}, HDAC9 (histone deacetylase 9) [NCBI Gene 9734] {aka HD7, HD7b, HD9, HDAC, HDAC7B, HDAC9B}, HDAC2 (histone deacetylase 2) [NCBI Gene 3066] {aka HD2, KDAC2, RPD3, YAF1}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, ICAM1 (intercellular adhesion molecule 1) [NCBI Gene 3383] {aka BB2, CD54, P3.58}, TNFSF9 (TNF superfamily member 9) [NCBI Gene 8744] {aka 4-1BB-L, CD137L, TNLG5A}, PRKCD (protein kinase C delta) [NCBI Gene 5580] {aka ALPS3, CVID9, MAY1, PKCD, nPKC-delta}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, ESR1 (estrogen receptor 1) [NCBI Gene 2099] {aka ER, ESR, ESRA, ESTRR, Era, NR3A1}, CTCF (CCCTC-binding factor) [NCBI Gene 10664] {aka CFAP108, FAP108, MRD21}, RBFOX2 (RNA binding fox-1 homolog 2) [NCBI Gene 23543] {aka FOX2, Fox-2, HNRBP2, HRNBP2, RBM9, RTA}, NLRP3 (NLR family pyrin domain containing 3) [NCBI Gene 114548] {aka AGTAVPRL, AII, AVP, C1orf7, CIAS1, CLR1.1}, KAT5 (lysine acetyltransferase 5) [NCBI Gene 10524] {aka ESA1, HTATIP, HTATIP1, NEDFASB, PLIP, TIP}, TSN (translin) [NCBI Gene 7247] {aka BCLF-1, C3PO, RCHF1, REHF-1, TBRBP, TRSLN}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, CASP3 (caspase 3) [NCBI Gene 836] {aka CPP32, CPP32B, SCA-1}, CD46 (CD46 molecule) [NCBI Gene 4179] {aka AHUS2, MCP, MIC10, TLX, TRA2.10}, TXNIP (thioredoxin interacting protein) [NCBI Gene 10628] {aka ARRDC6, EST01027, HHCPA78, THIF, VDUP1}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}
- **Diseases:** immunoblastic lymphoma (MESH:D016400), ABC-DLBCL (MESH:D016403), Hodgkin lymphomas (MESH:D006689), B cell malignancies (MESH:D016393), inflammatory (MESH:D007249), EBV associated (MESH:D020031), AIDS (MESH:D000163), gastric cancer (MESH:D013274), LCL (MESH:D002292), infectious mononucleosis (MESH:D007244), non-Hodgkin lymphomas (MESH:D008228), cytotoxicity (MESH:D064420), cancer (MESH:D009369), BL (MESH:D002051), NPC (MESH:D000077274), NK/T cell lymphomas (MESH:D016399), autoimmune diseases (MESH:D001327), AIDS immunoblastic lymphoma (MESH:D016483), infected (MESH:D007239), multiple myeloma (MESH:D009101)
- **Chemicals:** phorbol (MESH:C033085), Nanatinostat (MESH:C555225), water (MESH:D014867), Valganciclovir (MESH:D000077562), romidepsin (MESH:C087123), 4-hydroxytamoxifen (MESH:C016601), FITC (MESH:D016650), CO2 (MESH:D002245), arginine butyrate (MESH:C033018), Neon (MESH:D009356), vorinostat (MESH:D000077337), sodium butyrate (MESH:D020148), BrdU (MESH:D001973), Panobinostat (MESH:D000077767), PBS (MESH:D007854), GCV (MESH:D015774), BD (MESH:C028491), DMSO (MESH:D004121), CellTiter Glo (-), 7-AAD (MESH:C025942), PAA (MESH:D010746), belinostat (MESH:C487081), Pan (MESH:C041728)
- **Species:** human gammaherpesvirus 4 (Epstein Barr virus, no rank) [taxon 10376], Homo sapiens (human, species) [taxon 9606], herpesvirus [taxon 39059], Mus musculus (house mouse, species) [taxon 10090]
- **Cell lines:** P3HR1 — Homo sapiens (Human), EBV-related Burkitt lymphoma, Cancer cell line (CVCL_2676), Akata BL — Homo sapiens (Human), Burkitt lymphoma, Cancer cell line (CVCL_M562), BL — Homo sapiens (Human), EBV-related Burkitt lymphoma, Cancer cell line (CVCL_6354), IBL-1 — Homo sapiens (Human), AIDS-related immunoblastic lymphoma, Cancer cell line (CVCL_9638), BL41-B958 — Saguinus oedipus (Cotton-top tamarin), Transformed cell line (CVCL_1953), Jijoye — Homo sapiens (Human), EBV-related Burkitt lymphoma, Cancer cell line (CVCL_1317), BL41 — Homo sapiens (Human), Burkitt lymphoma, Cancer cell line (CVCL_1087)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC13029708/full.md

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13029708/full.md

## References

88 references — full list in the complete paper: https://tomesphere.com/paper/PMC13029708/full.md

---
Source: https://tomesphere.com/paper/PMC13029708