Viral reprogramming of host transcription initiation.

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

SpliceTools, a suite of downstream RNA splicing analysis tools to investigate mechanisms and impact of alternative splicing.

As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS.

A novel role of TRIM28 B box domain in L1 retrotransposition and ORF2p-mediated cDNA synthesis.

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

SON and Its Alternatively Spliced Isoforms Control MLL Complex-Mediated H3K4me3 and Transcription of Leukemia-Associated Genes.

Dysregulation of MLL complex-mediated histone methylation plays a pivotal role in gene expression associated with diseases, but little is known about cellular factors modulating MLL complex activity. Here, we report that SON, previously known as an RNA splicing factor, controls MLL complex-mediated transcriptional initiation. SON binds to DNA near transcription start sites, interacts with menin, and inhibits MLL complex assembly, resulting in decreased H3K4me3 and transcriptional repression. Importantly, alternatively spliced short isoforms of SON are markedly upregulated in acute myeloid leukemia. The short isoforms compete with full-length SON for chromatin occupancy but lack the menin-binding ability, thereby antagonizing full-length SON function in transcriptional repression while not impairing full-length SON-mediated RNA splicing. Furthermore, overexpression of a short isoform of SON enhances replating potential of hematopoietic progenitors. Our findings define SON as a fine-tuner of the MLL-menin interaction and reveal short SON overexpression as a marker indicating aberrant transcriptional initiation in leukemia.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells.

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of 'transitional' L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Organ-, sex- and age-dependent patterns of endogenous L1 mRNA expression at a single locus resolution.

Expression of L1 mRNA, the first step in the L1 copy-and-paste amplification cycle, is a prerequisite for L1-associated genomic instability. We used a reported stringent bioinformatics method to parse L1 mRNA transcripts and measure the level of L1 mRNA expressed in mouse and rat organs at a locus-specific resolution. This analysis determined that mRNA expression of L1 loci in rodents exhibits striking organ specificity with less than 0.8% of loci shared between organs of the same organism. This organ specificity in L1 mRNA expression is preserved in male and female mice and across age groups. We discovered notable differences in L1 mRNA expression between sexes with only 5% of expressed L1 loci shared between male and female mice. Moreover, we report that the levels of total L1 mRNA expression and the number and spectrum of expressed L1 loci fluctuate with age as independent variables, demonstrating different patterns in different organs and sexes. Overall, our comparisons between organs and sexes and across ages ranging from 2 to 22 months establish previously unforeseen dynamic changes in L1 mRNA expression in vivo. These findings establish the beginning of an atlas of endogenous L1 mRNA expression across a broad range of biological variables that will guide future studies.

Transcriptional signatures of Zika virus infection in astrocytes.

Astrocytes are an early and important target of Zika virus (ZIKV) infection in the developing brain, but the impacts of infection on astrocyte function remain controversial. Given that nonhuman primate (NHP) models of ZIKV infection replicate aspects of neurologic disease seen in human infections, we cultured primary astrocytes from the brain tissue of infant rhesus macaques and then infected the cells with Asian or African lineage ZIKV to identify transcriptional patterns associated with infection in these cells. The African lineage virus appeared to have greater infectivity and promote stronger antiviral signaling, but infection by either strain ultimately produced typical virus response patterns. Both viruses induced hypoxic stress, but the Asian lineage strain additionally had an effect on metabolic and lipid biosynthesis pathways. Together, these findings describe an NHP astrocyte model that may be used to assess transcriptional signatures following ZIKV infection.

Screen technical noise in single cell RNA sequencing data.

We proposed a data cleaning pipeline for single cell (SC) RNA-seq data, where we first screen genes (gene-wise screening) followed by screening cell libraries (library-wise screening). Gene-wise screening is based on the expectation that for a gene with a low technical noise, a gene's count in a library will tend to increase with the increase of library size, which was tested using negative binomial regression of gene count (as dependent variable) against library size (as independent variable). Library-wise screening is based on the expectation that across-library correlations for housekeeping (HK) genes is expected to be higher than the correlations for non-housekeeping (NHK) genes in those libraries with low technical noise. We removed those libraries, whose mean pairwise correlation for HK genes is NOT significantly higher than that for NHK genes. We successfully applied the pipeline to two large SC RNA-seq datasets. The pipeline was also developed into an R package.

The Epstein Barr virus circRNAome.

Our appreciation for the extent of Epstein Barr virus (EBV) transcriptome complexity continues to grow through findings of EBV encoded microRNAs, new long non-coding RNAs as well as the more recent discovery of over a hundred new polyadenylated lytic transcripts. Here we report an additional layer to the EBV transcriptome through the identification of a repertoire of latent and lytic viral circular RNAs. Utilizing RNase R-sequencing with cell models representing latency types I, II, and III, we identified EBV encoded circular RNAs expressed from the latency Cp promoter involving backsplicing from the W1 and W2 exons to the C1 exon, from the EBNA BamHI U fragment exon, and from the latency long non-coding RPMS1 locus. In addition, we identified circular RNAs expressed during reactivation including backsplicing from exon 8 to exon 2 of the LMP2 gene and a highly expressed circular RNA derived from intra-exonic backsplicing within the BHLF1 gene. While expression of most of these circular RNAs was found to depend on the EBV transcriptional program utilized and the transcription levels of the associated loci, expression of LMP2 exon 8 to exon 2 circular RNA was found to be cell model specific. Altogether we identified over 30 unique EBV circRNAs candidates and we validated and determined the structural features, expression profiles and nuclear/cytoplasmic distributions of several predominant and notable viral circRNAs. Further, we show that two of the EBV circular RNAs derived from the RPMS1 locus are detected in EBV positive clinical stomach cancer specimens. This study increases the known EBV latency and lytic transcriptome repertoires to include viral circular RNAs and it provides an essential foundation and resource for investigations into the functions and roles of this new class of EBV transcripts in EBV biology and diseases.

Maintenance of AP-2-Dependent Functional Activities of Nef Restricts Pathways of Immune Escape from CD8 T Lymphocyte Responses.

Nef-specific CD8+ T lymphocytes (CD8TL) are linked to extraordinary control of primate lentiviral replication, but the mechanisms underlying their efficacy remain largely unknown. The immunodominant, Mamu-B*017:01+-restricted Nef195-203MW9 epitope in SIVmac239 partially overlaps a sorting motif important for interactions with host AP-2 proteins and, hence, downmodulation of several host proteins, including Tetherin (CD317/BST-2), CD28, CD4, SERINC3, and SERINC5. We reasoned that CD8TL-driven evolution in this epitope might compromise Nef's ability to modulate these important molecules. Here, we used deep sequencing of SIV from nine B*017:01+ macaques throughout infection with SIVmac239 to characterize the patterns of viral escape in this epitope and then assayed the impacts of these variants on Nef-mediated modulation of multiple host molecules. Acute variation in multiple Nef195-203MW9 residues significantly compromised Nef's ability to downregulate surface Tetherin, CD4, and CD28 and reduced its ability to prevent SERINC5-mediated reduction in viral infectivity but did not impact downregulation of CD3 or major histocompatibility complex class I, suggesting the selective disruption of immunomodulatory pathways involving Nef AP-2 interactions. Together, our data illuminate a pattern of viral escape dictated by a selective balance to maintain AP-2-mediated downregulation while evading epitope-specific CD8TL responses. These data could shed light on mechanisms of both CD8TL-driven viral control generally and on Mamu-B*017:01-mediated viral control specifically.IMPORTANCE A rare subset of humans infected with HIV-1 and macaques infected with SIV can control the virus without aid of antiviral medications. A common feature of these individuals is the ability to mount unusually effective CD8 T lymphocyte responses against the virus. One of the most formidable aspects of HIV is its ability to evolve to evade immune responses, particularly CD8 T lymphocytes. We show that macaques that target a specific peptide in the SIV Nef protein are capable of better control of the virus and that, as the virus evolves to escape this response, it does so at a cost to specific functions performed by the Nef protein. Our results help show how the virus can be controlled by an immune response, which could help in designing effective vaccines.

A comprehensive approach to expression of L1 loci.

L1 elements represent the only currently active, autonomous retrotransposon in the human genome, and they make major contributions to human genetic instability. The vast majority of the 500 000 L1 elements in the genome are defective, and only a relatively few can contribute to the retrotransposition process. However, there is currently no comprehensive approach to identify the specific loci that are actively transcribed separate from the excess of L1-related sequences that are co-transcribed within genes. We have developed RNA-Seq procedures, as well as a 1200 bp 5΄ RACE product coupled with PACBio sequencing that can identify the specific L1 loci that contribute most of the L1-related RNA reads. At least 99% of L1-related sequences found in RNA do not arise from the L1 promoter, instead representing pieces of L1 incorporated in other cellular RNAs. In any given cell type a relatively few active L1 loci contribute to the 'authentic' L1 transcripts that arise from the L1 promoter, with significantly different loci seen expressed in different tissues.

TAF4b Regulates Oocyte-Specific Genes Essential for Meiosis.

TAF4b is a gonadal-enriched subunit of the general transcription factor TFIID that is implicated in promoting healthy ovarian aging and female fertility in mice and humans. To further explore the potential mechanism of TAF4b in promoting ovarian follicle development, we analyzed global gene expression at multiple time points in the human fetal ovary. This computational analysis revealed coordinate expression of human TAF4B and critical regulators and effectors of meiosis I including SYCP3, YBX2, STAG3, and DAZL. To address the functional relevance of this analysis, we turned to the embryonic Taf4b-deficient mouse ovary where, for the first time, we demonstrate, severe deficits in prophase I progression as well as asynapsis in Taf4b-deficient oocytes. Accordingly, TAF4b occupies the proximal promoters of many essential meiosis and oogenesis regulators, including Stra8, Dazl, Figla, and Nobox, and is required for their proper expression. These data reveal a novel TAF4b function in regulating a meiotic gene expression program in early mouse oogenesis, and support the existence of a highly conserved TAF4b-dependent gene regulatory network promoting early oocyte development in both mice and women.

p53 Enables metabolic fitness and self-renewal of nephron progenitor cells.

Contrary to its classic role in restraining cell proliferation, we demonstrate here a divergent function of p53 in the maintenance of self-renewal of the nephron progenitor pool in the embryonic mouse kidney. Nephron endowment is regulated by progenitor availability and differentiation potential. Conditional deletion of p53 in nephron progenitor cells (Six2Cre(+);p53(fl/fl)) induces progressive depletion of Cited1(+)/Six2(+) self-renewing progenitors and loss of cap mesenchyme (CM) integrity. The Six2(p53-null) CM is disorganized, with interspersed stromal cells and an absence of a distinct CM-epithelia and CM-stroma interface. Impaired cell adhesion and epithelialization are indicated by decreased E-cadherin and NCAM expression and by ineffective differentiation in response to Wnt induction. The Six2Cre(+);p53(fl/fl) cap has 30% fewer Six2(GFP(+)) cells. Apoptotic index is unchanged, whereas proliferation index is significantly reduced in accordance with cell cycle analysis showing disproportionately fewer Six2Cre(+);p53(fl/fl) cells in the S and G2/M phases compared with Six2Cre(+);p53(+/+) cells. Mutant kidneys are hypoplastic with fewer generations of nascent nephrons. A significant increase in mean arterial pressure is observed in early adulthood in both germline and conditional Six2(p53-null) mice, linking p53-mediated defects in kidney development to hypertension. RNA-Seq analyses of FACS-isolated wild-type and Six2(GFP(+)) CM cells revealed that the top downregulated genes in Six2Cre(+);p53(fl/fl) CM belong to glucose metabolism and adhesion and/or migration pathways. Mutant cells exhibit a ∼ 50% decrease in ATP levels and a 30% decrease in levels of reactive oxygen species, indicating energy metabolism dysfunction. In summary, our data indicate a novel role for p53 in enabling the metabolic fitness and self-renewal of nephron progenitors.

Reduced mitochondrial activity in colonocytes facilitates AMPKα2-dependent inflammation.

Intestinal inflammation is associated with low levels of mucosal ATP, highlighting the importance of mitochondrial function associated with ATP production in the pathophysiology of the disease. In the inflamed colon of humans and mice, we found decreased levels of mitochondrial complex cytochrome c oxidase I/IV and lower ATP levels. Thus, we generated colonic ρ0 cells with reduced mitochondrial function linked to ATP production by selective depletion of mitochondrial DNA. In these cells, RNA sequencing revealed a substantial number of differentially expressed transcripts, among which 240 belonged to inflammatory pathways activated in human inflamed colon and TNF-α-treated cells (false discovery rate < 0.05). TNF-α treatment of colonic ρ0 cells augmented IL-8 expression by 9-fold (P < 0.01) via NF-κB compared to TNF-α-treated control. Moreover, reduced mitochondrial function facilitated TNF-α-mediated NF-κB luciferase promoter activity as a result of lowered inhibitory IκBα (nuclear factor of κ light polypeptide gene enhancer in B-cell inhibitor, α), leading to elevated NF-κB. In cells with reduced mitochondrial function, TNF-α facilitated AMPKα2 activation by 8-fold (P < 0.01), which was involved in NF-κB-dependent IL-8 expression. Last, in human and mouse colon, anti-TNF-α treatment restored reduced mitochondria-dependent inflammation. We propose that selective targeting of this novel mechanism provides new treatment opportunities for intestinal inflammation.-Heller, S., Penrose, H. M., Cable, C., Biswas, D., Nakhoul, H., Baddoo, M., Flemington, E., Crawford, S. E., Savkovic, S. D. Reduced mitochondrial activity in colonocytes facilitates AMPKα2-dependent inflammation.

Global transcript structure resolution of high gene density genomes through multi-platform data integration.

Annotation of herpesvirus genomes has traditionally been undertaken through the detection of open reading frames and other genomic motifs, supplemented with sequencing of individual cDNAs. Second generation sequencing and high-density microarray studies have revealed vastly greater herpesvirus transcriptome complexity than is captured by existing annotation. The pervasive nature of overlapping transcription throughout herpesvirus genomes, however, poses substantial problems in resolving transcript structures using these methods alone. We present an approach that combines the unique attributes of Pacific Biosciences Iso-Seq long-read, Illumina short-read and deepCAGE (Cap Analysis of Gene Expression) sequencing to globally resolve polyadenylated isoform structures in replicating Epstein-Barr virus (EBV). Our method, Transcriptome Resolution through Integration of Multi-platform Data (TRIMD), identifies nearly 300 novel EBV transcripts, quadrupling the size of the annotated viral transcriptome. These findings illustrate an array of mechanisms through which EBV achieves functional diversity in its relatively small, compact genome including programmed alternative splicing (e.g. across the IR1 repeats), alternative promoter usage by LMP2 and other latency-associated transcripts, intergenic splicing at the BZLF2 locus, and antisense transcription and pervasive readthrough transcription throughout the genome.

High-fat diet induced leptin and Wnt expression: RNA-sequencing and pathway analysis of mouse colonic tissue and tumors.

Obesity, an immense epidemic affecting approximately half a billion adults, has doubled in prevalence in the last several decades. Epidemiological data support that obesity, due to intake of a high-fat, western diet, increases the risk of colon cancer; however, the mechanisms underlying this risk remain unclear. Here, utilizing next generation RNA sequencing, we aimed to determine the high-fat diet (HFD) mediated expression profile in mouse colon and the azoxymethane/dextran sulfate sodium model of colon cancer. Mice on HFD had significantly higher colonic inflammation, tumor burden, and a number of differentially expressed transcripts compared to mice on regular diet (RD). We identified 721 transcripts differentially expressed in mouse HFD colon that were in a shared pattern with colonic tumors (RD and HFD). Importantly, in mouse colon, HFD stimulated an expression signature strikingly similar to human colon cancer, especially those with inflammatory microsatellite instability. Furthermore, pathway analysis of these transcripts demonstrated their association with active inflammation and colon cancer signaling, with leptin and Wnt as the top two transcripts elevated in mouse HFD colon shared with tumors. Moreover, in mouse colon, HFD-stimulated tumorigenic Wnt pathway activation was further validated by upregulation of ß-catenin transcriptional targets. Finally, in human colon cancer, upregulation of leptin pathway members was shown with a large network of dysregulated transcripts being linked with worse overall survival.

Aging is associated with impaired angiogenesis, but normal microvascular network structure, in the rat mesentery.

A big problem associated with aging is thought to be impaired microvascular growth or angiogenesis. However, to link the evidence for impaired angiogenesis to microvascular dysfunction in aged tissues, we must compare adult vs. aged microvascular networks in unstimulated scenarios. The objective of this study was to test the hypothesis that aged microvascular networks are characterized by both fewer vessels and the impaired ability to undergo angiogenesis. Mesentery tissues from adult (9-mo) and aged (24-mo) male Fischer 344 rats were harvested and immunolabeled for platelet/endothelial cell adhesion molecule (an endothelial cell marker) according to two scenarios: unstimulated and stimulated. For unstimulated groups, tissues harvested from adult and aged rats were compared. For stimulated groups, tissues were harvested 3 or 10 days after compound 48/80-induced mast cell degranulation stimulation. Unstimulated aged microvascular networks displayed larger mean vascular area per tissue area compared with the unstimulated adult networks. The lack of a decrease in vessel density was supported at the gene expression level with RNA-Seq analysis and with comparison of vessel densities in soleus muscle. Following stimulation, capillary sprouting and vessel density were impaired in aged networks at 3 and 10 days, respectively. Our results suggest that aging associated with impaired angiogenesis mechanisms might not influence normal microvascular function, since unstimulated aged microvascular networks can display a "normal adult-like" vessel density and architecture. NEW & NOTEWORTHY: Using a multidimensional approach, we present evidence supporting that aged microvascular networks display vessel density and patterning similar to adult networks despite also being characterized by a decreased capacity to undergo angiogenesis. Thus, vessel loss is not necessarily a characteristic of aging.

Latent Expression of the Epstein-Barr Virus (EBV)-Encoded Major Histocompatibility Complex Class I TAP Inhibitor, BNLF2a, in EBV-Positive Gastric Carcinomas.

The Epstein-Barr virus (EBV) BNLF2a gene product provides immune evasion properties to infected cells through inhibition of transporter associated with antigen processing (TAP)-mediated transport of antigen peptides. Although BNLF2a is considered to be a lytic gene, we demonstrate that it is expressed in nearly half of the EBV-associated gastric carcinomas analyzed. Further, we show that BNLF2a expression is dissociated from lytic gene expression. BNLF2a is therefore expressed in this latency setting, potentially helping protect the infected tumor cells from immunosurveillance.

High-throughput RNA sequencing-based virome analysis of 50 lymphoma cell lines from the Cancer Cell Line Encyclopedia project.

UNLABELLED: Using high-throughput RNA sequencing data from 50 common lymphoma cell culture models from the Cancer Cell Line Encyclopedia project, we performed an unbiased global interrogation for the presence of a panel of 740 viruses and strains known to infect human and other mammalian cells. This led to the findings of previously identified infections by Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV), and human T-lymphotropic virus type 1 (HTLV-1). In addition, we also found a previously unreported infection of one cell line (DEL) with a murine leukemia virus. High expression of murine leukemia virus (MuLV) transcripts was observed in DEL cells, and we identified four transcriptionally active integration sites, one being in the TNFRSF6B gene. We also found low levels of MuLV reads in a number of other cell lines and provided evidence suggesting cross-contamination during sequencing. Analysis of HTLV-1 integrations in two cell lines, HuT 102 and MJ, identified 14 and 66 transcriptionally active integration sites with potentially activating integrations in immune regulatory genes, including interleukin-15 (IL-15), IL-6ST, STAT5B, HIVEP1, and IL-9R. Although KSHV and EBV do not typically integrate into the genome, we investigated a previously identified integration of EBV into the BACH2 locus in Raji cells. This analysis identified a BACH2 disruption mechanism involving splice donor sequestration. Through viral gene expression analysis, we detected expression of stable intronic RNAs from the EBV BamHI W repeats that may be part of long transcripts spanning the repeat region. We also observed transcripts at the EBV vIL-10 locus exclusively in the Hodgkin's lymphoma cell line, Hs 611.T, the expression of which were uncoupled from other lytic genes. Assessment of the KSHV viral transcriptome in BCP-1 cells showed expression of the viral immune regulators, K2/vIL-6, K4/vIL-8-like vCCL1, and K5/E2-ubiquitin ligase 1 that was significantly higher than expression of the latency-associated nuclear antigen. Together, this investigation sheds light into the virus composition across these lymphoma model systems and provides insights into common viral mechanistic principles. IMPORTANCE: Viruses cause cancer in humans. In lymphomas the Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV) and human T-lymphotropic virus type 1 are major contributors to oncogenesis. We assessed virus-host interactions using a high throughput sequencing method that facilitates the discovery of new virus-host associations and the investigation into how the viruses alter their host environment. We found a previously unknown murine leukemia virus infection in one cell line. We identified cellular genes, including cytokine regulators, that are disrupted by virus integration, and we determined mechanisms through which virus integration causes deregulation of cellular gene expression. Investigation into the KSHV transcriptome in the BCP-1 cell line revealed high-level expression of immune signaling genes. EBV transcriptome analysis showed expression of vIL-10 transcripts in a Hodgkin's lymphoma that was uncoupled from lytic genes. These findings illustrate unique mechanisms of viral gene regulation and to the importance of virus-mediated host immune signaling in lymphomas.

New Noncoding Lytic Transcripts Derived from the Epstein-Barr Virus Latency Origin of Replication, oriP, Are Hyperedited, Bind the Paraspeckle Protein, NONO/p54nrb, and Support Viral Lytic Transcription.

UNLABELLED: We have previously shown that the Epstein-Barr virus (EBV) likely encodes hundreds of viral long noncoding RNAs (vlncRNAs) that are expressed during reactivation. Here we show that the EBV latency origin of replication (oriP) is transcribed bi-directionally during reactivation and that both leftward (oriPtLs) and rightward (oriPtRs) transcripts are largely localized in the nucleus. While the oriPtLs are most likely noncoding, at least some of the oriPtRs contain the BCRF1/vIL10 open reading frame. Nonetheless, oriPtR transcripts with long 5' untranslated regions may partially serve noncoding functions. Both oriPtL and oriPtR transcripts are expressed with late kinetics, and their expression is inhibited by phosphonoacetic acid. RNA sequencing (RNA-seq) analysis showed that oriPtLs and oriPtRs exhibited extensive "hyperediting" at their Family of Repeat (FR) regions. RNA secondary structure prediction revealed that the FR region of both oriPtLs and oriPtRs may form large evolutionarily conserved and thermodynamically stable hairpins. The double-stranded RNA-binding protein and RNA-editing enzyme ADAR was found to bind to oriPtLs, likely facilitating editing of the FR hairpin. Further, the multifunctional paraspeckle protein, NONO, was found to bind to oriPt transcripts, suggesting that oriPts interact with the paraspeckle-based innate antiviral immune pathway. Knockdown and ectopic expression of oriPtLs showed that it contributes to global viral lytic gene expression and viral DNA replication. Together, these results show that these new vlncRNAs interact with cellular innate immune pathways and that they help facilitate progression of the viral lytic cascade. IMPORTANCE: Recent studies have revealed that the complexity of lytic herpesviral transcriptomes is significantly greater than previously appreciated with hundreds of viral long noncoding RNAs (vlncRNAs) being recently discovered. Work on cellular lncRNAs over the past several years has just begun to give us an initial appreciation for the array of functions they play in complex formation and regulatory processes in the cell. The newly identified herpesvirus lncRNAs are similarly likely to play a variety of different functions, although these functions are likely tailored to specific needs of the viral infection cycles. Here we describe novel transcripts derived from the EBV latency origin of replication. We show that they are hyperedited, that they interact with a relatively newly appreciated antiviral pathway, and that they play a role in facilitating viral lytic gene expression. These investigations are a starting point to unraveling the complex arena of vlncRNA function in herpesvirus lytic replication.