Abstract
Polyadenylated nuclear (PAN) RNA is a long noncoding transcript involved in Kaposi's sarcoma-associated herpesvirus (KSHV) lytic reactivation and regulation of cellular and viral gene expression. We have previously shown that PAN RNA has dynamic secondary structure and protein binding profiles that can be influenced by epitranscriptomic modifications. N6-methyladenosine (m6A) is one of the most abundant chemical signatures found in viral RNA genomes and virus-encoded RNAs. Here, we combined antibody-independent next-generation mapping with direct RNA sequencing to address the epitranscriptomic status of PAN RNA in KSHV infected cells. We showed that PAN m6A status is dynamic, reaching the highest number of modifications at the late lytic stages of KSHV infection. Using a newly developed method, termed selenium-modified deoxythymidine triphosphate (SedTTP)-reverse transcription (RT) and ligation assisted PCR analysis of m6A (SLAP), we gained insight into the fraction of modification at identified sites. By applying comprehensive proteomic approaches, we identified writers and erasers that regulate the m6A status of PAN, and readers that can convey PAN m6A phenotypic effects. We verified the temporal and spatial subcellular availability of the methylome components for PAN modification by performing confocal microscopy analysis. Additionally, the RNA biochemical probing (SHAPE-MaP) outlined local and global structural alterations invoked by m6A in the context of full-length PAN RNA. This work represents the first comprehensive overview of the dynamic interplay that takes place between the cellular epitranscriptomic machinery and a specific viral RNA in the context of KSHV infected cells.