- Research article
- Open Access
Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA
BMC Genomics volume 14, Article number: 543 (2013)
Epstein-Barr virus (EBV) is a human herpesvirus implicated in cancer and autoimmune disorders. Little is known concerning the roles of RNA structure in this important human pathogen. This study provides the first comprehensive genome-wide survey of RNA and RNA structure in EBV.
Novel EBV RNAs and RNA structures were identified by computational modeling and RNA-Seq analyses of EBV. Scans of the genomic sequences of four EBV strains (EBV-1, EBV-2, GD1, and GD2) and of the closely related Macacine herpesvirus 4 using the RNAz program discovered 265 regions with high probability of forming conserved RNA structures. Secondary structure models are proposed for these regions based on a combination of free energy minimization and comparative sequence analysis. The analysis of RNA-Seq data uncovered the first observation of a stable intronic sequence RNA (sisRNA) in EBV. The abundance of this sisRNA rivals that of the well-known and highly expressed EBV-encoded non-coding RNAs (EBERs).
This work identifies regions of the EBV genome likely to generate functional RNAs and RNA structures, provides structural models for these regions, and discusses potential functions suggested by the modeled structures. Enhanced understanding of the EBV transcriptome will guide future experimental analyses of the discovered RNAs and RNA structures.
Epstein-Barr virus (EBV) is widely disseminated in the human population. Upwards of 95% of the adult human population is infected with EBV . EBV is implicated in a number of different cancers, including Hodgkin’s disease , nasopharyngeal carcinoma , hepatocellular carcinomas , lymphoepithelioma-like carcinomas , some breast cancers , and in autoimmune disorders such as Sjögren’s syndrome , dermatomyositis , lupus , rheumatoid arthritis , and multiple sclerosis . EBV was the first cancer-associated virus to be discovered when in 1964  it was isolated from tumors occurring in African children (Burkitt’s lymphoma ). Despite intense investigation for more than 50 years, the precise roles played by the virus in these diseases remains to be elucidated.
The ~170000 bp genome of EBV is a linear double-stranded (ds)DNA that circularizes to form an episome in the host cell nucleus. Infection occurs by entry of the EBV virion into human host epithelial cells and initially proceeds via an aggressive lytic phase. The virus migrates to B cells where it causes persistent lifelong infections marked by extended periods of latency with interspersed lytic reactivation . EBV latency proceeds via three distinct programs, each expressing a different set of coding and non-coding viral gene products. Viral latent gene products rewire B cells to evade the host immune system and propagate the virus . In a manner not yet fully understood, this rewiring increases the tumorigenic potential of EBV-infected cells, including making infected cells resistant to apoptotic pathways that would otherwise kill cancerous cells .
The EBV transcriptome is complex, consisting of multiple pre-mRNA transcripts that undergo extensive alternative splicing events and can, at various states of infection, yield different sets of gene products . Two well-studied EBV-encoded ncRNAs (EBER1 and EBER2) are known to be expressed throughout infection , as are subsets of as many as 50 different EBV-encoded microRNAs (miRNAs) [19, 20]. Although the precise functions of the EBERs have remained obscure, they are the most highly expressed RNAs in EBV-infected cells (~107 copies per cell ). Roles in cancer have been proposed  and the EBERs are capable of inducing tumors in immune-suppressed mice . Both EBERs form ribonucleoprotein complexes [23–26], including host proteins, and adopt well-defined secondary structures . Except for the EBERs, RNA structure is a relatively understudied aspect of EBV virology; an analysis of RNA secondary structure in EBV could enable significant advances in understanding this important human pathogen.
RNA secondary structure plays important roles in many viruses, mediating functions as diverse as translation initiation , catalysis , viral genome packaging , alternative splicing [31, 32], interactions with the host innate immune system , RNA lifetime/stability [34–36], and regulation of gene expression . The discovery of RNA structure is an active area of research [38–40] and is aided by the availability of sophisticated non-coding (nc)RNA-discovery algorithms [41–43]. A successful strategy for discovering conserved RNA structures is implemented in the program RNAz . In this approach, homologous sequences are aligned and divided into fragments. For each fragment, calculated parameters that include measures of the thermodynamic stability (z-score ) and conservation of secondary structure are used to make predictions based on a support vector machine (SVM) trained on known ncRNAs, which generates a probability classifier (p-class, scale of 0.0 to 1.0) for the fragment containing conserved secondary structure. The RNAz approach has been used to discover RNA secondary structures in a number of viruses [38, 46–48].
Results and discussion
RefSeq genomes (~170 kbp long) from four EBV strains (NC_007605.1, NC_009334.1, AY961628.3, HQ020558.1) and one from the closely-related (Figure 1) Macacine herpesvirus 4 (MHV4, NC_006146.1; also known as rhesus lymphocryptovirus) were aligned using the MAFFT alignment program [49, 50]. The average pairwise sequence identity (APSI) of the alignment is 84.8%; the four EBV strains were close in sequence (APSI of 95.6%), while MHV4 had an APSI of 68.6% versus EBV.
Across the entire EBV genome the z-score average is -0.19 and the average p-class is 0.13. RNAz discovered 265 regions that map to known transcripts from the EBV genome (named locus 1 to 265 in Additional files 1 and 2) and have strong evidence of locally conserved and stable RNA structure. These regions (average length ~200 nt) are comprised of overlapping windows with high predicted probability of generating conserved RNA structure. The average p-class for the set of predictions was 0.91 with minimum and maximum values of 0.61 and 1.00, respectively (Additional file 1). The thermodynamic z-score is an important component in determining the p-class. This parameter reflects the thermodynamic stability of the native sequence predicted structure versus that of random sequence  (e.g. measures greater than expected folding stability). For predicted structured regions in EBV RNA, the average z-score was -2.38 with minimum and maximum values of -1.18 and -9.71, respectively. This analysis, however, does not account for long-range RNA interactions that may occur outside the overlapped prediction windows.
Identification of known EBV ncRNAs
The genomic region specifying EBER1 occurs within locus 12 (Additional file 1). The structure model derived for EBER1 (Figure 2A) is almost identical to the model in the Rfam database (ID# RF01789). A single base pair in the Rfam model, between nt 6760 and 6770, is shifted to form a pair between 6646 and 6770, and the predicted closing AU pair of this hairpin is open in the Rfam model (Figure 2A). Additionally, C6649 is in a bulge-loop in the Rfam model, but the insertion of a G residue across from this nt allows a GC pair to form. EBER2, surprisingly, does not occur in an RNAz-predicted locus. The region containing EBER2 was not calculated to have a favorable z-score; however, this is most likely due to the inability of the primary sequence alignment to capture the conservation of structure for this RNA. EBER2 is less conserved in sequence than EBER1 (50.3% sequence conservation vs. 78.3% for EBER1) making primary sequence alignment difficult in this region. Common helices in EBER2 were not properly aligned, thus reducing the predicted free energy of the consensus structure, which resulted in an unfavorable z-score.
All but two of the EBV miRNAs annotated in the miRBase database  occur in RNAz-predicted structured regions. Each of the mature miRNA sequences, when mapped to predicted structure models, falls within a hairpin RNA structure that resembles a canonical pre-miRNA (Additional files 1 and 3). For example, Figure 2B shows three hairpins in locus 209, which include the BART3, BART4, and BART1 miRNAs. Each of these hairpins places the 5P and 3P mature sequences in canonical structural contexts: the mature miRNAs are offset and imperfectly base pair with each other. This is also true for the other 36 mature miRNA sequences that fall within predicted structured regions (Additional file 3).
Finally, a reported EBV-encoded viral small nucleolar RNA (v-snoRNA; ) also falls within an RNAz-predicted structured region (locus 235 in Additional file 1). The v-snoRNA corresponds to nt 154352 to 154416 (Additional file 3) and, in isolation, has a predicted folding free energy and z-score of -17.8 kcal/mol and -0.81, respectively (Figure 2C).
In addition to the regions corresponding to EBER1, the v-snoRNA, and 14 loci specifying EBV pre-miRNAs, there are 249 additional EBV regions with predicted conserved RNA structure (Additional file 1).
Structure within introns
Sixty predicted structured regions overlap introns. Introns are a fertile source of structured ncRNAs [54, 55] that play roles in a number of important biological processes: snoRNAs , miRNAs , and piRNAs  are expressed from introns. RNA structure within introns can alter accessibility or distance between functional elements and thus regulate splicing .
The most extensive regions of predicted structure in introns occur within a repetitive region of the EBV genome (the W repeats). The number of W repeats can vary, and the optimal number appears to be five to eight . The W repeats are transcribed at the 5′ end of the long (~100000 nt) primary transcript encoding the Epstein-Barr nuclear antigen (EBNA) latent proteins (Figure 3). EBNA transcripts are particularly interesting because they have high nuclear-to-cytoplasmic abundance, suggesting non-coding nuclear functions .
Six EBNA proteins are expressed from mRNAs produced by extensive alternative splicing. Transcription of the EBNA RNA can occur from an upstream promoter, the C promoter (Cp), or a downstream W promoter (Wp) . Early in latency III Wp is utilized, where splicing of W0 to W1 generates the start codon for EBNA-LP production ; later in latency III there is a switch to Cp where the non-coding C2 exon is joined to W1 to create the EBNA-LP start codon (Figure 3). Within the W repeats two short coding exons, W1 and W2, which are 66 and 132 nt long, respectively, are joined to form the EBNA-LP protein open reading frame (ORF) by removal of a long (2791 nt) and short (81 nt) intron. RNAz predicts widespread conserved and stable RNA structure covering 49% of the long W repeat intron. There are five predicted regions, which are identical in sequence in the W repeats (Figure 3): loci 18, 23, 28, 33, 38, 43, 48, and 53 are identical; loci 19, 24, 29, 34, 39, 44, 49, and 54 are identical; loci 20, 25, 30, 35, 40, 45, and 50 are identical; loci 21, 26, 31, 36, 41, 46, and 51 are identical; and finally, loci 22, 27, 32, 37, 42, 47, and 52 are identical. Loci 18, 19, and 20, although they are multiplied in the W repeats, are transcribed only when Cp is used to initiate transcription (Figure 3). Cp occurs just upstream of locus 17 and the C1 non-coding exon is within locus 17 (Figure 3).The first copy of Wp overlaps locus 21, which also contains W0 (Figure 3).
Within loci 18, 23, 28, 33, 38, 43, 48, and 53 is a highly unusual modeled RNA structure, a very long (586 nt) tetraloop hairpin (Figure 4A). The INFERNAL package [64, 65] was used to search for additional homologous structures. Three other lymphocryptovirus species [MHV4, Callitrichine herpesvirus 3 (CaHV3), and Cercopithecine herpesvirus 12 (CeHV12, also known as papiine herpesvirus)], members of the herpesvirus taxonomical group to which EBV belongs (Figure 1), possess homologous structures. These discovered hairpins occur within the long introns of repeats homologous to the W repeats of EBV, appearing six times in MHV4 (NC_006146.1 nt 12738 to 13335, 15820 to 16416, 18901 to 19497, 21982 to 22578, 25063 to 25659, and 28144 to 28740) and 4 times in CaHV3 (NC_004367.1 in the negative strand nt 120959 to 121293, 123871 to 124205, 126783 to 127117, and 129695 to 130029). Three hairpin homologs are found within the partial sequence of CeHV12 (AF200364.1 nt 6498 to 7112, 9542 to 10156, 12586 to 13200).
The APSI of these long hairpin sequences is 66.7%, with the CaHV3 sequence being most distant (APSI 45.1% versus EBV). The MHV4 and CeHV12 hairpins are most similar to EBV in sequence (APSI of 80.5% and 74.5%, respectively) and structure (Figure 4B and C and Additional file 4). The EBV, MHV4, and CeHV12 hairpins are each capped with a tetraloop closed by a GC pair (Figure 4A, B, and C) and only the first nt in the loop varies (UUGG, CUGG, and GUGG, respectively). The more distant CaHV3 hairpin, in contrast, is capped with a seven-nt loop (Figure 4D). Perhaps the identities of the capping nucleotides are of lesser importance to the function of this long hairpin than the overall conserved structure, which is marked by extensively base paired stem regions interspersed with internal (or bulged) loops (Figure 4). Internal loops may adopt interesting 2D and 3D structures; for example, the three-by-three internal loop that occurs 20 bps from the hairpin loop in EBV [Figure 4A, also conserved in MHV4 (Figure 4B)], may fold into a “3 RRs” motif  where the opposing internal loop nucleotides form consecutive sheared purine-purine pairs.
Additionally, global thermodynamic z-scores were calculated for each hairpin (Figure 4): all four have very negative z-scores (ranging from -6.3 to -5.9). Thus, the structures generated by these sequences are more than six standard deviations more thermodynamically stable than random sequences with the same dinucleotide content. This indicates that evolution is acting to preserve base pairing within each sequence, as well as the striking structural homology between these viral hairpins.
Very long RNA hairpins that are structurally similar to the EBV W repeat hairpin have been described in humans  and in Caenorhabditis elegans. These structures are excellent substrates for adenosine deaminase acting on RNA (ADAR) editing enzymes, which convert adenosine to inosine (A-to-I editing). Interestingly, a previous study of A-to-I editing in EBV  found evidence for an editing site that maps to a region of the genome that we predict to be part of the long hairpin (Figure 4A). ADARs act on structured RNAs and have preferences for adenosines in particular sequence contexts , allowing the prediction of potential A-to-I editing sites and of the extent of editing. The reported editing site is predicted to be moderately edited (12.9% by ADAR2, Additional file 4) and is modeled to be in a CA internal loop (Figure 4A). There are additional sites predicted to be strongly edited in the EBV W repeat hairpin and in the homologous structures found in MHV4, CeHV12, and CaHV3 (Additional file 4). Few of these strongly predicted sites (seven), however, have apparent conservation (Additional file 4), which is unsurprising, as selection of ADAR sites may be partially stochastic .
RNA editing in introns, typically within Alu inserts , has been proposed to regulate RNA splicing . Inosine can form wobble pairs with G and C residues and can contribute to RNA-protein interactions: thus A-to-I editing can alter important RNA-RNA and RNA-protein interactions . Editing can alter splicing regulatory elements or create/destroy splice sites . Perhaps editing of this long intronic hairpin also plays a role in splicing or in other functions post-splicing.
Discovery of an EBV-encoded sisRNA
A small RNA library (30-200 nt size range) was constructed using nuclear RNA isolated from cultured human B lymphocytes stably infected with EBV (BJAB-B1 cells ). These cells exhibit a latency III program of gene expression, the most transcriptionally active program, which also produces latency I and II transcripts. Nuclear RNAs include not only EBER1 and EBER2  but also intronic sequences that are excised during the splicing process. Unsurprisingly, the abundant EBER1 and EBER2 are present in the library and had the highest number of aligned RNA-Seq reads (Figure 5 and Table 1). A significant peak (~7.8% EBER1 abundance) corresponding to the v-snoRNA (Figure 5 and Table 1) was also observed. There was a surprisingly large peak originating from the W repeat region (Figure 5). This peak corresponds to reads that cover the entirety of the small intron (81 nt) that separates the W1 and W2 ORFs (Figure 3). Next to the EBERs, this was the largest peak observed. Comparing length-normalized reads between the intron and EBER1, the intron is estimated to be present at 21% the level of EBER1 or roughly on a par with EBER2, which is estimated to be 25% as abundant as EBER1 (Table 1). The presence of the intron in total nuclear RNA samples was confirmed by Northern blot (Figure 6A) and RT-PCR (Figure 6B). RT-qPCR then revealed that the intron is ~9-fold enriched in nuclear versus cytoplasmic RNA (Figure 6C).
Stable intronic sequence (sis)RNAs have been described in Xenopus oocytes  and, interestingly, in herpes simplex virus-1 (HSV-1, an α-herpesvirus)  and in cytomegalovirus (a β- herpesvirus) . In these other herpesviruses a stable intron, known as the latency-associated transcript (LAT) in HSV-1, is implicated in maintenance of latency . The EBV sisRNA (ebv-sisRNA-1) is the first of this class of ncRNAs to be described in EBV. The ebv-sisRNA-1 differs from the LAT in two major respects: the LAT is much larger (>2 kb), and the functional form of the LAT is believed to be the lariat-intron splicing intermediate, whereas the ebv-sisRNA-1 is likely to be a linear molecule. For both ends of ebv-sisRNA-1 to have been free for ligation in the small RNA library construction step, the 5′ end of the released intron could not have been sequestered in a lariat structure. This suggests that ebv-sisRNA-1 has been acted upon by the debranching enzyme, which also linearizes intron-derived snoRNAs during the maturation of this important class of ncRNAs .
Introns homologous to ebv-sisRNA-1 were found in other lymphocryptoviruses (Figure 6D). The ebv-sisRNA-1 sequence and these RNAs have an APSI of 83.3%, with the CaHV3 sequence being most distant (57.6% APSI). In all analyzed EBV strains (EBV-1, EBV-2, GD1, and GD2), the 81-nt ebv-sisRNA-1 sequence is 100% conserved, whereas the APSI of the genome sequences of these strains is 95.6%. Within the CA-rich region an 11-nt sequence, from positions 49 to 59, is 100% conserved in the analyzed lymphocryptoviruses (Figure 6D). The distribution of dinucleotide frequencies in the ebv-sisRNA-1 sequence is skewed, with CA being the most frequent, followed by AC and CC. The presence of CA-rich regions in the sisRNA is interesting as such sequences are able to bind hnRNP L protein and modulate splicing .
Although ebv-sisRNA-1 did not overlap an RNAz-predicted region, the secondary structure for this RNA was modeled on the basis of the eight sequences shown in Figure 6E. Two small hairpins are predicted (Figure 6E): an upstream hairpin that places the U-rich motif in a loop and a weak (3 base pair) hairpin in the CA-rich region. A compensatory mutation (double point-mutation that preserves base pairing) is observed in the more stable upstream hairpin (Figure 6E, a conserved GC pair is converted to a UA pair in CaHV3), which suggests that evolution is acting to preserve this structure. It is interesting to find evidence of a conserved hairpin that places the U-rich motif in a loop, as U-rich sequences often bind proteins  and hairpins can pre-organize RNA sequence motifs for RNA-protein interactions . The CA-rich region is likely unstructured and free to associate with other interactors.
Additional putative sisRNAs in EBV
Other than the W repeat sisRNAs, EBV possesses 18 additional short introns (Table 2) interspersed throughout its genome. These short introns range in conservation from 82.9% to 100.0% APSI (average of 97.2%) and are predicted to fold into structures with a range of thermodynamic stabilities (from -0.4 to -59.6 kcal/mol). Four sequences have significantly negative z-scores [< -1.0, from nt 11466 to 11610, 80362 to 80446, 87014 to 87087, and 105929 to 106056 (Table 2)] and one intron has a moderately favorable z-score (-0.54) from nt 1496 to 1573. These short introns also overlap with RNAz-predicted structured regions (Table 2).
Particularly interesting is the sequence from nt 105929 to 106056, which has the most favorable predicted folding free energy and z-score of any EBV intron (Table 2). This intron occurs within the lytically-expressed BBLF2/3 gene, which generates part of a helix-primase complex essential for replication of the EBV genome . The modeled secondary structure of the intron is shown in Figure 7A; it forms a tetraloop hairpin that includes all of the intron sequence. This intron structure is contained within RNAz-predicted locus 152, where the basal stem of the hairpin extends into coding regions by 25 nt (Additional file 3, nt 66692 to 66844 in the reverse genome sequence) and the splice donor and acceptor are base paired within helices. A similar hairpin structure, identified by INFERNAL, occurs in the MHV4 genome (Figure 7B). Rather than being joined together via splicing, MHV4 BBLF2/3 is translated from two ORFs that produce BBLF2 and BBLF3 as separate proteins. This same intronic region in EBV partially codes for BBLF2 in MHV4 (the translation stop codon occurs within the UUAA tetraloop) and is also present in the 5′ UTR of BBLF3 (Figure 7B). It is remarkable that the hairpin structure has been maintained evolutionarily between EBV and MHV4 despite different apparent roles in protein expression: in EBV it brings splice sites closer together in space and may be processed out to form an independent structured RNA, whereas in MHV4 it presents a stop codon in a loop and is not processed out of the mRNA.
Four EBV short introns have significant numbers of mapped reads in the small RNA library (Table 1); however, they are lower in abundance than ebv-sisRNA-1. In the EBV latent membrane protein (LMP-2B) gene, for example, there are three short introns with significant numbers of aligned read peaks (nt 459 to 539, nt 789 to 870, and nt 952 to 1025 in Table 1). Indeed, the intron encompassing nt 952 to 1025 has an abundance of ~5.7% relative to EBER1, which is similar to the value observed for the v-snoRNA (7.8%, Table 2). The short introns in EBV, including those that are expressed during latency III and those with few or no mapped reads, may represent a family of sisRNAs. It will be interesting to learn if they are differentially expressed in other viral latency programs or during the lytic phase. Short introns are an attractive source of small ncRNAs, as they are released as by-products of pre-mRNA splicing.
Structure overlapping splice sites
The importance of RNA structure in the regulation of splicing has been highlighted in two excellent reviews [59, 87]. Structure can be inhibitory, by sequestering splice sites  or regulatory elements , or it can enhance splicing by presenting sites in an accessible conformation  or by bringing splice sites into closer proximity with each other . Conformational switching between accessible and inaccessible structures may also regulate splicing .
In EBV, there are seven predicted structured regions that include splice sites (Additional file 1). Of these, three include a splice donor or acceptor, while the remaining regions include both sites. For example, in the LMP-2B gene, locus 1 overlaps a splice donor site (after nt 788, Additional file 3), which is buried within the stem of a tetraloop hairpin; towards the 3′ end of this structure (nt 827 to 833) is a putative intronic splicing enhancer (Figure 8A). Locus 2 spans an entire intron and includes both the donor and acceptor sites (after nt 1495 and 1573, respectively, in Additional file 3). Additionally, there is a putative exonic splicing enhancer at the beginning of locus 2 (nt 1475 to 1482). The modeled structure for locus 2 (Figure 8B) organizes this region into a complex multibranched structure where the splice donor and acceptor sites are each in helices: the donor at the end of the basal stem and the acceptor within a hairpin loop. Interestingly, each splice site occurs directly after or near a predicted UU/UU two-by-two internal loop motif, in which UU mismatches maintain stable [93, 94] helical 3D conformations . Internal loops and non-canonical pairs can present complex patterns of chemical functional groups in protein-accessible regions of a folded RNA , which may influence the binding of splicing factors. Additionally, the 3D folding of this structure may bring splice sites close together .
Structure overlapping coding regions
Many (150) predictions overlap EBV coding regions, in some cases including a significant portion of an ORF. In the lytically-expressed EBV FGAM-synthase gene, for example, predictions cover 40% the coding nt (locus 3 to locus 9; Additional file 1). More dramatically, the BGLF3 gene overlaps four predicted structured regions (locus 163 to locus 166; Additional file 1), which span 84% of the coding region including part of the start codon. Interestingly, many viral genomes appear to be enriched for RNA secondary structure [97, 98], including within coding regions [38, 99, 100]. RNA structure in coding regions can have a variety of functions , such as altering reading frames via frame-shifting pseudoknots  or by inducing ribosomal pausing , which can impact numerous co-translational processes [103, 104]. The identification of conserved and stable RNA secondary structures within EBV coding regions suggests that they may play important roles in this virus.
Twenty-one predicted structured regions include translation start or stop codons. Locus 263, for example, includes the start codon for the latent membrane protein 2A (LMP-2A), which is involved in the regulation of host B-cell signaling and is implicated in EBV-associated tumorigenesis . The first 63 nt of locus 263 fall within the 5′ UTR of the LMP-2A mRNA (Figure 8C). The start codon is predicted to occur within the stem helix of a multibranch loop. RNA structure can play a role in regulating translation through modulating the accessibility of codons. Strong RNA structure at stop codons, for example, can stimulate stop codon read-through , whereas modulation of start codon accessibility can affect translational efficiency . Perhaps the stable structure (z-score -1.48, Figure 8C) that encompasses the LMP-2A start codon modulates the expression of this important gene.
Bioinformatics analyses of the EBV genome revealed 265 regions with putative conserved RNA structure. The computational screen successfully identified most of the known EBV ncRNAs: EBER1, a virally-encoded snoRNA, and 42 miRNAs. Other predicted structured regions of the EBV genome occur at or near sites that suggest potential functions for RNA structures. A very long hairpin structure in the W repeat region that is conserved in lymphocryptoviruses contains at least one A-to-I editing site, with possible implications for splicing. Additional regions span splice sites or start and stop codons. These structure predictions may prove useful since RNA and RNA structures are attractive targets for chemotherapy [108–111]. The modeled RNA structures are rich in internal loops, e.g. the UU/UU internal loops near two splice sites in LMP-2B, which are potentially drugable structural motifs . RNA splicing, in particular, is an attractive target for RNA therapeutics [113, 114].
RNA-Seq data identified several putative sisRNAs, including one in the W repeat region (ebv-sisRNA-1) that is abundant and localizes in the nucleus. These results lay the groundwork for future experiments to further characterize the proposed structures and to understand the functional significance of the discovered RNAs and RNA structures in EBV. All structural data are available in the Additional files. In addition, data for the W repeat long hairpin and ebv-sisRNA-1 will be deposited in the Rfam database.
Structure prediction and modeling
EBV RefSeq genomes from four different viral strains, EBV-1 (NC_007605.1), EBV-2 (NC_009334.1), GD1 (AY961628.3), and GD2 (HQ020558.1), plus one closely-related genome from Macacine herpesvirus 4 (NC_006146.1), were obtained from the NCBI nt database. Genome sequences were aligned using the MAFFT program [49, 50] implementing the FFT-NS-1 alignment strategy. The alignment was processed into 120-nt windows with a 10-nt step-size using the Perl script rnazWindow.pl provided with the RNAz program. Predictions of windows with likely conserved and stable RNA structure were made for both the plus and minus strands with the RNAz 2.1 [42, 44] program using a dinucleotide randomization model. The RNAz output was processed to combine windows with higher probability of structure (p-class > 0.5) into loci. These loci were further filtered to include only those that had at least one window with p-class > 0.9. When predictions had ambiguous strand orientations, e.g. occurred in a genome location that could be transcribed from either strand, the program RNAstrand  was used to determine if structure was more likely to be in the forward or reverse strand.
From each predicted locus the EBV-2 sequence was extracted and used in a BLAST  search against the NCBI nt database. The returned homologous sequences were compiled and each locus was aligned using the MAFFT program implementing the Q-INS-i method. Initial models for these alignments were made using RNAalifold . Alignments and predicted consensus structures were converted to Stockholm format alignments and used to build covariance models with the INFERNAL package . These covariance models were used to search for potentially homologous sequences/structures in a database of all available herpesvirus RefSeq genomes. Returned sequences were added to the alignments and structural predictions were made using RNAalifold as well as the program TurboFold  followed by manual model refinement. Resulting models were used to perform additional iterative INFERNAL searches and structure revisions until no further model improvement could be made (finding the structure that was best conserved in the alignment).
2D renderings of modeled RNAs were generated using the PseudoViewer program  and manually processed using the open-source graphics program Inkscape and ImageJ [120, 121]. The Gibbs free energy of folding at 37°C [122, 123] for models was calculated using the efn2 method implemented in the RNAstructure package . The z-scores for sequences were calculated by generating sets of 200 dinucleotide randomizations using the SIMMONICS program  and then calculating energies with RNAfold . The difference in native free energy versus the randomized set was normalized by the standard deviation to give the z-score.
Small RNA-Seq analysis
EBV-positive BJAB-B1 cells were maintained at 37°C and 5% CO2 in RPMI media supplemented with FBS (10%) and penicillin . Cells that had been expanded from 1 ml frozen stock to 200 ml culture volume were grown for an additional 48 hrs and were then pelleted by centrifugation at 1000 × g for 5 min at 4°C. Pellets were washed 3× with ice-cold PBS buffer (plus DTT). The pellet, on ice, was resuspended in cell disruption buffer  (10 mM KCl, 1.5 mM MgCl2, 20 mM Tris-Cl, 1 mM DTT, 0.1% Triton-X) and moved to a chilled Dounce homogenizer (Kontes, type B). Homogenization proceeded on ice until mostly free nuclei were observed under microscope. The lysate was transferred to centrifuge bottles and spun at 1500 × g for 5 min at 4°C. The supernatant, containing the cytoplasmic fraction, was removed to a fresh tube and a 3× volume of Trizol reagent was added. The nuclear fraction containing pellet was dissolved in 3 ml of Trizol reagent and both fractions were flash frozen in a dry-ice/ethanol bath and stored at -80°C overnight. Each sample was chloroform extracted and then precipitated at room temperature for 10 min using isopropanol. RNA was pelleted by centrifugation at 10000 × g for 10 min at 4°C. Pellets were resuspended in water and treated with RNase-free DNase I according to manufacturer’s protocol (Promega). RNA was further purified/recovered by RNeasy RNA purification kit (Qiagen). Samples were treated again using DNase I to remove any traces of genomic DNA contamination and recovered using the RNeasy kit. RNA integrity was checked using an Agilent BioAnalyzer: the RNA integrity number for the undigested samples was equal to 10, indicating that sample annealing did not degrade the RNA. RNA samples were used to build libraries using a small RNA library kit (Illumina), where specific adapters were ligated to the 5′ and 3′ ends, before barcoding and cDNA synthesis. The cDNA was fragmented and size-selected to contain fragments in the range of 30 to 200 bp. RNA-Seq data, 75 bp reads, were acquired on an Illumina HiSeq instrument run in paired-end mode. Reads were aligned to the EBV-2 RefSeq genome (NC_009334.1) using the Bowtie 2 program  (data given in Additional file 5).
Analysis of ebv-sisRNA-1 expression
20 μg of total RNA from EBV-positive BJAB-B1 and EBV-negative (but otherwise isogenic) BJAB cells were fractionated by denaturing 10% PAGE. RNA was transferred to an Amersham Hybond-N + positively charged nylon membrane then cross-linked to the membrane using UV radiation. The membrane was incubated at 37°C for 1 hr in pre-hybridization buffer (GE) and washed 3× in double-distilled water (ddH2O). The membrane was then incubated overnight at 37°C with 5′-radiolabeled DNA probe complementary to nts 37 to 66 of the 106 nt human U6 snRNA (NR_046494.1), probe sequence: 5′-GCAGGGGCCATGCTAATCTTCTCTGTATCG-3′. The blot was washed 3× in wash buffer (0.005% SDS and 3× SSC), wrapped in plastic film, exposed to a phosphor screen for 16 hrs and imaged on a Storm phosphorimager. The blot was rehybridized with 5′-radiolabeled DNA probe complementary to nt 19 to 46 of ebv-sisRNA-1, probe sequence: 5′- TGTGGTGGAGTGTTGGGCTTAGCAGAAA -3′.
For the RT-PCR and RT-qPCR analyses, 2 μg nuclear or cytoplasmic RNA were used to generate cDNA with the High Capacity cDNA synthesis kit (Applied Biosystems). 5 μl 100× diluted cDNA was used for PCR with EBER1 primers complementary to the transcript’s 5′ and 3′ ends: FWD 5′-AGGACCTACGCTGCCCTAGA-3′ REV 5′-AAAACATGCGGACCACCAGCTGG-3′, and ebv-sisRNA-1 primers complementary to the 5′ end and nt 52 to 71: FWD 5′-GTAAGTGGACTTTAATTTTTTCTGCTAAGCCC-3′ REV 5′-TGGGTGTGTGTAGTGTGTGC-3′. Template was combined with 5 μl 10× PCR buffer (NEB), 1 μl 10 mM dNTPs, 1 μl each of 10 μM primer and 0.25 μl Taq (NEB) and ddH2O to a final volume of 50 μl. The thermocycler protocol was: (i) initial denaturation at 90°C for 10 min, (ii) melt at 90°C for 15 s, (iii) anneal at 53°C for 30 s, (iv) extend at 60°C for 30 s, (v) cycle back to step ii 40×. Amplification products were analyzed by denaturing 10% PAGE and stained with ethidium bromide. Amplicon sizes were determined by comparison to a 10 bp ladder (Invitrogen).
For RT-qPCR a master mix was made of FastStart Universal SYBR Green Master Mix (Roche), ddH2O and each primer (final concentration of 300 nM). 45 μl aliquots were added to a 96-well plate before adding 5 μl 100× dilute cDNA (three biological replicates, each run with three technical replicates, including no template and reverse transcriptase controls). In addition to the ebv-sisRNA-1 primers, amplification of nt 8 to 104 of the nuclear-enriched U2 snRNA (NR_002716.3) used primers: FWD 5′-CTCGGCCTTTTGGCTAAGAT-3′ REV 5′-TATTCCATCTCCCTGCTCCA-3′; and nt 994 to 1071 of the cytoplasmically-enriched 18S rRNA (NR_003286.2) used primers: FWD 5′-CGAAAGCATTTGCCAAGAAT-3′ REV 5′-GCATCGTTTATGGTCGGAAC-3′. The same thermocycler protocol used for PCR was used for the qPCR run, and data was collected on an Applied Biosystems StepOnePlus Real-Time PCR system. Data were exported to Microsoft Excel and analyzed using the Delta Ct method to find the ratio of nuclear to cytoplasmic abundance.
Average pairwise sequence identity
Callitrichine herpesvirus 3
Cercopithecine herpesvirus 12
Macacine herpesvirus 4
Stable intronic sequence RNA
Small nucleolar RNA
Epstein-Barr nuclear antigen
Open reading frame
Adenosine deaminase acting on RNA
Latency associated transcript
Epstein-Barr virus encoded RNA.
Epstein-Barr virus and infectious mononucleosis.http://www.cdc.gov/ncidod/diseases/ebv.htm,
Flavell KJ, Murray PG: Hodgkin’s disease and the Epstein-Barr virus. Mol Pathol: MP. 2000, 53 (5): 262-269. 10.1136/mp.53.5.262.
Burgos JS: Involvement of the Epstein-Barr virus in the nasopharyngeal carcinoma pathogenesis. Med Oncol. 2005, 22 (2): 113-121. 10.1385/MO:22:2:113.
Akhter S, Liu H, Prabhu R, DeLucca C, Bastian F, Garry RF, Schwartz M, Thung SN, Dash S: Epstein-Barr virus and human hepatocellular carcinoma. Cancer Letters. 2003, 192 (1): 49-57. 10.1016/S0304-3835(02)00695-X.
Begin LR, Eskandari J, Joncas J, Panasci L: Epstein-Barr virus related lymphoepithelioma-like carcinoma of lung. J Surg Oncol. 1987, 36 (4): 280-283. 10.1002/jso.2930360413.
Hippocrate A, Oussaief L, Joab I: Possible role of EBV in breast cancer and other unusually EBV-associated cancers. Cancer Letters. 2011, 305 (2): 144-149. 10.1016/j.canlet.2010.11.007.
Wen S, Shimizu N, Yoshiyama H, Mizugaki Y, Shinozaki F, Takada K: Association of Epstein-Barr virus (EBV) with Sjogren’s syndrome: differential EBV expression between epithelial cells and lymphocytes in salivary glands. Am J Pathol. 1996, 149 (5): 1511-1517.
Chen DY, Chen YM, Lan JL, Chen HH, Hsieh CW, Wey SJ, Lu JJ: Polymyositis/dermatomyositis and nasopharyngeal carcinoma: the Epstein-Barr virus connection?. J Clin Virol. 2010, 49 (4): 290-295. 10.1016/j.jcv.2010.08.015.
Parks CG, Cooper GS, Hudson LL, Dooley MA, Treadwell EL, St Clair EW, Gilkeson GS, Pandey JP: Association of Epstein-Barr virus with systemic lupus erythematosus: effect modification by race, age, and cytotoxic T lymphocyte-associated antigen 4 genotype. Arthritis Rheum. 2005, 52 (4): 1148-1159. 10.1002/art.20997.
Alspaugh MA, Jensen FC, Rabin H, Tan EM: Lymphocytes transformed by Epstein-Barr virus. Induction of nuclear antigen reactive with antibody in rheumatoid arthritis. J Exp Med. 1978, 147 (4): 1018-1027. 10.1084/jem.147.4.1018.
Levin LI, Munger KL, Rubertone MV, Peck CA, Lennette ET, Spiegelman D, Ascherio A: Multiple sclerosis and Epstein-Barr virus. JAMA: J Am Med Assoc. 2003, 289 (12): 1533-1536. 10.1001/jama.289.12.1533.
Epstein MA, Achong BG, Barr YM: Virus particles in cultured lymphoblasts from Burkitt’s lymphoma. Lancet. 1964, 1 (7335): 702-703.
Burkitt D: A sarcoma involving the jaws in African children. Brit J Surg. 1958, 46 (197): 218-223. 10.1002/bjs.18004619704.
Young LS, Dawson CW, Eliopoulos AG: The expression and function of Epstein-Barr virus encoded latent genes. Mol Pathol: MP. 2000, 53 (5): 238-247. 10.1136/mp.53.5.238.
Ning S: Innate immune modulation in EBV infection. Herpesviridae. 2011, 2 (1): 1-10.1186/2042-4280-2-1.
Takada K: Role of Epstein-Barr virus in Burkitt’s lymphoma. Curr Top Microbiol Immunol. 2001, 258: 141-151. 10.1007/978-3-642-56515-1_9.
Young LS, Arrand JR, Murray PG: EBV gene expression and regulation. Human Herpes viruses: biology, therapy, and immunoprophylaxis. Edited by: Arvin A, Campadelli-Fiume G, Mocarski E, Moore PS, Roizman B, Whitley R, Yamanishi K. 2007, Cambridge: Cambridge University Press
Lerner MR, Andrews NC, Miller G, Steitz JA: Two small RNAs encoded by Epstein-Barr virus and complexed with protein are precipitated by antibodies from patients with systemic lupus erythematosus. Proc Natl Acad Sci USA. 1981, 78 (2): 805-809. 10.1073/pnas.78.2.805.
Chen SJ, Chen GH, Chen YH, Liu CY, Chang KP, Chang YS, Chen HC: Characterization of Epstein-Barr virus miRNAome in nasopharyngeal carcinoma by deep sequencing. PloS One. 2010, 5 (9): e12745-10.1371/journal.pone.0012745.
Riley KJ, Rabinowitz GS, Yario TA, Luna JM, Darnell RB, Steitz JA: EBV and human microRNAs co-target oncogenic and apoptotic viral and human genes during latency. EMBO J. 2012, 31 (9): 2207-2221. 10.1038/emboj.2012.63.
Yoshizaki T, Endo K, Ren Q, Wakisaka N, Murono S, Kondo S, Sato H, Furukawa M: Oncogenic role of Epstein-Barr virus-encoded small RNAs (EBERs) in nasopharyngeal carcinoma. Auris, Nasus, Larynx. 2007, 34 (1): 73-78. 10.1016/j.anl.2006.09.025.
Yamamoto N, Takizawa T, Iwanaga Y, Shimizu N: Malignant transformation of B lymphoma cell line BJAB by Epstein-Barr virus-encoded small RNAs. FEBS Letters. 2000, 484 (2): 153-158. 10.1016/S0014-5793(00)02145-1.
Conrad NK, Fok V, Cazalla D, Borah S, Steitz JA: The challenge of viral snRNPs. Cold Spring Harb Symp Quant Biol. 2006, 71: 377-384. 10.1101/sqb.2006.71.057.
Fok V, Mitton-Fry RM, Grech A, Steitz JA: Multiple domains of EBER 1, an Epstein-Barr virus noncoding RNA, recruit human ribosomal protein L22. RNA. 2006, 12 (5): 872-882. 10.1261/rna.2339606.
Toczyski DP, Matera AG, Ward DC, Steitz JA: The Epstein-Barr virus (EBV) small RNA EBER1 binds and relocalizes ribosomal protein L22 in EBV-infected human B lymphocytes. Proc Natl Acad Sci USA. 1994, 91 (8): 3463-3467. 10.1073/pnas.91.8.3463.
Lee N, Pimienta G, Steitz JA: AUF1/hnRNP D is a novel protein partner of the EBER1 noncoding RNA of Epstein-Barr virus. RNA. 2012, 18 (11): 2073-2082. 10.1261/rna.034900.112.
Glickman JN, Howe JG, Steitz JA: Structural analyses of EBER1 and EBER2 ribonucleoprotein particles present in Epstein-Barr virus-infected cells. J Virol. 1988, 62 (3): 902-911.
Pfingsten JS, Costantino DA, Kieft JS: Structural basis for ribosome recruitment and manipulation by a viral IRES RNA. Science. 2006, 314 (5804): 1450-1454. 10.1126/science.1133281.
Butcher SE: Structure and function of the small ribozymes. Curr Opin Struct Biol. 2001, 11 (3): 315-320. 10.1016/S0959-440X(00)00207-4.
Johnson SF, Telesnitsky A: Retroviral RNA dimerization and packaging: the what, how, when, where, and why. PLoS Pathog. 2010, 6 (10): e1001007-10.1371/journal.ppat.1001007.
Abbink TE, Berkhout B: RNA structure modulates splicing efficiency at the human immunodeficiency virus type 1 major splice donor. J Virol. 2008, 82 (6): 3090-3098. 10.1128/JVI.01479-07.
Moss WN, Dela-Moss LI, Priore SF, Turner DH: The influenza A segment 7 mRNA 3′ splices site pseudoknot/hairpin family. RNA Biol. 2012, 9 (11): 1305-1310.
Malathi K, Saito T, Crochet N, Barton DJ, Gale M, Silverman RH: RNase L releases a small RNA from HCV RNA that refolds into a potent PAMP. RNA. 2010, 16 (11): 2108-2119. 10.1261/rna.2244210.
Conrad NK, Shu MD, Uyhazi KE, Steitz JA: Mutational analysis of a viral RNA element that counteracts rapid RNA decay by interaction with the polyadenylate tail. Proc Natl Acad Sci USA. 2007, 104 (25): 10412-10417. 10.1073/pnas.0704187104.
Mitton-Fry RM, DeGregorio SJ, Wang J, Steitz TA, Steitz JA: Poly (A) tail recognition by a viral RNA element through assembly of a triple helix. Science. 2010, 330 (6008): 1244-1247. 10.1126/science.1195858.
Tycowski KT, Shu MD, Borah S, Shi M, Steitz JA: Conservation of a triple-helix-forming RNA stability element in noncoding and genomic RNAs of diverse viruses. Cell Rep. 2012, 2 (1): 26-32. 10.1016/j.celrep.2012.05.020.
Giedroc DP, Cornish PV: Frameshifting RNA pseudoknots: structure and mechanism. Virus Res. 2009, 139 (2): 193-208. 10.1016/j.virusres.2008.06.008.
Moss WN, Priore SF, Turner DH: Identification of potential conserved RNA secondary structure throughout influenza A coding regions. RNA. 2011, 17 (6): 991-1011. 10.1261/rna.2619511.
Hofacker IL, Fekete M, Flamm C, Huynen MA, Rauscher S, Stolorz PE, Stadler PF: Automatic detection of conserved RNA structure elements in complete RNA virus genomes. Nucleic Acids Res. 1998, 26 (16): 3825-3836. 10.1093/nar/26.16.3825.
Tuplin A, Wood J, Evans DJ, Patel AH, Simmonds P: Thermodynamic and phylogenetic prediction of RNA secondary structures in the coding region of hepatitis C virus. RNA. 2002, 8 (6): 824-841. 10.1017/S1355838202554066.
Mathews DH, Moss WN, Turner DH: Folding and finding RNA secondary structure. Cold Spring Harb Perspect Biol. 2010, 2 (12): a003665-10.1101/cshperspect.a003665.
Washietl S, Hofacker IL, Stadler PF: Fast and reliable prediction of noncoding RNAs. Proc Natl Acad Sci USA. 2005, 102 (7): 2454-2459. 10.1073/pnas.0409169102.
Schroeder SJ: Advances in RNA structure prediction from sequence: new tools for generating hypotheses about viral RNA structure-function relationships. J Virol. 2009, 83 (13): 6326-6334. 10.1128/JVI.00251-09.
Gruber AR, Findeiss S, Washietl S, Hofacker IL, Stadler PF: RNAz 2.0: improved noncoding RNA detection. Pac Symp Biocomp. 2010, 15: 69-79.
Clote P, Ferre F, Kranakis E, Krizanc D: Structural RNA has lower folding energy than random RNA of the same dinucleotide frequency. RNA. 2005, 11 (5): 578-591. 10.1261/rna.7220505.
Hofacker IL, Stadler PF, Stocsits RR: Conserved RNA secondary structures in viral genomes: a survey. Bioinformatics. 2004, 20 (10): 1495-1499. 10.1093/bioinformatics/bth108.
Lodeiro MF, Filomatori CV, Gamarnik AV: Structural and functional studies of the promoter element for dengue virus RNA replication. J Virol. 2009, 83 (2): 993-1008. 10.1128/JVI.01647-08.
Cordey S, Gerlach D, Junier T, Zdobnov EM, Kaiser L, Tapparel C: The cis-acting replication elements define human enterovirus and rhinovirus species. RNA. 2008, 14 (8): 1568-1578. 10.1261/rna.1031408.
Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33 (2): 511-518. 10.1093/nar/gki198.
Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30 (14): 3059-3066. 10.1093/nar/gkf436.
Ehlers B, Spiess K, Leendertz F, Peeters M, Boesch C, Gatherer D, McGeoch DJ: Lymphocryptovirus phylogeny and the origins of Epstein-Barr virus. J Gen Virol. 2010, 91 (Pt 3): 630-642.
Hutzinger R, Feederle R, Mrazek J, Schiefermeier N, Balwierz PJ, Zavolan M, Polacek N, Delecluse HJ, Huttenhofer A: Expression and processing of a small nucleolar RNA from the Epstein-Barr virus genome. PLoS Pathog. 2009, 5 (8): e1000547-10.1371/journal.ppat.1000547.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 (Database issue): D154-D158.
Rearick D, Prakash A, McSweeny A, Shepard SS, Fedorova L, Fedorov A: Critical association of ncRNA with introns. Nucleic Acids Res. 2011, 39 (6): 2357-2366. 10.1093/nar/gkq1080.
St Laurent G, Shtokalo D, Tackett MR, Yang Z, Eremina T, Wahlestedt C, Urcuqui-Inchima S, Seilheimer B, McCaffrey TA, Kapranov P: Intronic RNAs constitute the major fraction of the non-coding RNA in mammalian cells. BMC Genomics. 2012, 13: 504-10.1186/1471-2164-13-504.
Dieci G, Preti M, Montanini B: Eukaryotic snoRNAs: a paradigm for gene expression flexibility. Genomics. 2009, 94 (2): 83-88. 10.1016/j.ygeno.2009.05.002.
Lin SL, Kim H, Ying SY: Intron-mediated RNA interference and microRNA (miRNA). Front Biosci. 2008, 13: 2216-2230. 10.2741/2836.
Beyret E, Liu N, Lin H: piRNA biogenesis during adult spermatogenesis in mice is independent of the ping-pong mechanism. Cell Res. 2012, 22 (10): 1429-1439. 10.1038/cr.2012.120.
Warf MB, Berglund JA: Role of RNA structure in regulating pre-mRNA splicing. Trends Biochem Sci. 2010, 35 (3): 169-178. 10.1016/j.tibs.2009.10.004.
Tierney RJ, Kao KY, Nagra JK, Rickinson AB: Epstein-Barr virus BamHI W repeat number limits EBNA2/EBNA-LP coexpression in newly infected B cells and the efficiency of B-cell transformation: a rationale for the multiple W repeats in wild-type virus strains. J Virol. 2011, 85 (23): 12362-12375. 10.1128/JVI.06059-11.
Concha M, Wang X, Cao S, Baddoo M, Fewell C, Lin Z, Hulme W, Hedges D, McBride J, Flemington EK: Identification of new viral genes and transcript isoforms during Epstein-Barr virus reactivation using RNA-Seq. J Virol. 2012, 86 (3): 1458-1467. 10.1128/JVI.06537-11.
Woisetschlaeger M, Yandava CN, Furmanski LA, Strominger JL, Speck SH: Promoter switching in Epstein-Barr virus during the initial stages of infection of B lymphocytes. Proc Natl Acad Sci USA. 1990, 87 (5): 1725-1729. 10.1073/pnas.87.5.1725.
Rogers RP, Woisetschlaeger M, Speck SH: Alternative splicing dictates translational start in Epstein-Barr virus transcripts. EMBO J. 1990, 9 (7): 2273-2277.
Nawrocki EP, Kolbe DL, Eddy SR: Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009, 25 (10): 1335-1337. 10.1093/bioinformatics/btp157.
Gardner PP: The use of covariance models to annotate RNAs in whole genomes. Brief Funct Genomic Proteomic. 2009, 8 (6): 444-450. 10.1093/bfgp/elp042.
Arvey A, Tempera I, Tsai K, Chen HS, Tikhmyanova N, Klichinsky M, Leslie C, Lieberman PM: An atlas of the Epstein-Barr virus transcriptome and epigenome reveals host-virus regulatory interactions. Cell Host Microbe. 2012, 12 (2): 233-245. 10.1016/j.chom.2012.06.008.
Lerman YV, Kennedy SD, Shankar N, Parisien M, Major F, Turner DH: NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions. RNA. 2011, 17 (9): 1664-1677. 10.1261/rna.2641911.
Morse DP, Aruscavage PJ, Bass BL: RNA hairpins in noncoding regions of human brain and Caenorhabditis elegans mRNA are edited by adenosine deaminases that act on RNA. Proc Natl Acad Sci USA. 2002, 99 (12): 7906-7911. 10.1073/pnas.112704299.
Morse DP, Bass BL: Long RNA hairpins that contain inosine are present in Caenorhabditis elegans poly (A) + RNA. Proc Natl Acad Sci USA. 1999, 96 (11): 6048-6053. 10.1073/pnas.96.11.6048.
Eggington JM, Greene T, Bass BL: Predicting sites of ADAR editing in double-stranded RNA. Nat Commun. 2011, 2: 319-
Kim DD, Kim TT, Walsh T, Kobayashi Y, Matise TC, Buyske S, Gabriel A: Widespread RNA editing of embedded alu elements in the human transcriptome. Genome Res. 2004, 14 (9): 1719-1725. 10.1101/gr.2855504.
Lev-Maor G, Ram O, Kim E, Sela N, Goren A, Levanon EY, Ast G: Intronic Alus influence alternative splicing. PLoS Genet. 2008, 4 (9): e1000204-10.1371/journal.pgen.1000204.
Bond CS, Fox AH: Paraspeckles: nuclear bodies built on long noncoding RNA. J Cell Biol. 2009, 186 (5): 637-644. 10.1083/jcb.200906113.
Fresen KO, Hausen H: Establishment of EBNA-expressing cell lines by infection of Epstein-Barr virus (EBV)-genome-negative human lymphoma cells with different EBV strains. Int J Cancer J Int du Cancer. 1976, 17 (2): 161-166. 10.1002/ijc.2910170203.
Howe JG, Steitz JA: Localization of Epstein-Barr virus-encoded small RNAs by in situ hybridization. Proc Natl Acad Sci USA. 1986, 83 (23): 9006-9010. 10.1073/pnas.83.23.9006.
Thorvaldsdottir H, Robinson JT, Mesirov JP: Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013, 14 (2): 178-192. 10.1093/bib/bbs017.
Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP: Integrative genomics viewer. Nat Biotechnol. 2011, 29 (1): 24-26. 10.1038/nbt.1754.
Gardner EJ, Nizami ZF, Talbot CC, Gall JG: Stable intronic sequence RNA (sisRNA), a new class of noncoding RNA from the oocyte nucleus of Xenopus tropicalis. Genes Dev. 2012, 26 (22): 2550-2559. 10.1101/gad.202184.112.
Farrell MJ, Dobson AT, Feldman LT: Herpes simplex virus latency-associated transcript is a stable intron. Proc Natl Acad Sci USA. 1991, 88 (3): 790-794. 10.1073/pnas.88.3.790.
Kulesza CA, Shenk T: Murine cytomegalovirus encodes a stable intron that facilitates persistent replication in the mouse. Proc Natl Acad Sci USA. 2006, 103 (48): 18302-18307. 10.1073/pnas.0608718103.
Hesselberth JR: Lives that introns lead after splicing. WIREs RNA. 2013, 10.1002/wrna.1187.
Petfalski E, Dandekar T, Henry Y, Tollervey D: Processing of the precursors to small nucleolar RNAs and rRNAs requires common components. Mol Cell Biol. 1998, 18 (3): 1181-1189.
Hui J, Hung LH, Heiner M, Schreiner S, Neumuller N, Reither G, Haas SA, Bindereif A: Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing. EMBO J. 2005, 24 (11): 1988-1998. 10.1038/sj.emboj.7600677.
You Y, Chen CY, Shyu AB: U-rich sequence-binding proteins (URBPs) interacting with a 20-nucleotide U-rich sequence in the 3′ untranslated region of c-fos mRNA may be involved in the first step of c-fos mRNA degradation. Mol Cell Biol. 1992, 12 (7): 2931-2940.
Alian A, DeGiovanni A, Griner SL, Finer-Moore JS, Stroud RM: Crystal structure of an RluF-RNA complex: a base-pair rearrangement is the key to selectivity of RluF for U2604 of the ribosome. J Mol Biol. 2009, 388 (4): 785-800. 10.1016/j.jmb.2009.03.029.
Liao G, Huang J, Fixman ED, Hayward SD: The Epstein-Barr virus replication protein BBLF2/3 provides an origin-tethering function through interaction with the zinc finger DNA binding protein ZBRK1 and the KAP-1 corepressor. J Virol. 2005, 79 (1): 245-256. 10.1128/JVI.79.1.245-256.2005.
Buratti E, Baralle FE: Influence of RNA secondary structure on the pre-mRNA splicing process. Mol Cell Biol. 2004, 24 (24): 10505-10514. 10.1128/MCB.24.24.10505-10514.2004.
Blanchette M, Chabot B: A highly stable duplex structure sequesters the 5′ splice site region of hnRNP A1 alternative exon 7B. RNA. 1997, 3 (4): 405-419.
Sirand-Pugnet P, Durosay P, Clouet d’Orval BC, Brody E, Marie J: Beta-Tropomyosin pre-mRNA folding around a muscle-specific exon interferes with several steps of spliceosome assembly. J Mol Biol. 1995, 251 (5): 591-602. 10.1006/jmbi.1995.0458.
Buratti E, Muro AF, Giombi M, Gherbassi D, Iaconcig A, Baralle FE: RNA folding affects the recruitment of SR proteins by mouse and human polypurinic enhancer elements in the fibronectin EDA exon. Mol Cell Biol. 2004, 24 (3): 1387-1400. 10.1128/MCB.24.3.1387-1400.2004.
Libri D, Stutz F, McCarthy T, Rosbash M: RNA structural patterns and splicing: molecular basis for an RNA-based enhancer. RNA. 1995, 1 (4): 425-436.
Moss WN, Dela-Moss LI, Kierzek E, Kierzek R, Priore SF, Turner DH: The 3′ splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin. PloS One. 2012, 7 (6): e38323-10.1371/journal.pone.0038323.
SantaLucia J, Kierzek R, Turner DH: Stabilities of consecutive A.C, C.C, G.G, U.C, and U.U mismatches in RNA internal loops: Evidence for stable hydrogen-bonded U.U and C.C. + pairs. Biochemistry. 1991, 30 (33): 8242-8251. 10.1021/bi00247a021.
Wu M, McDowell JA, Turner DH: A periodic table of symmetric tandem mismatches in RNA. Biochemistry. 1995, 34 (10): 3204-3211. 10.1021/bi00010a009.
Sashital DG, Venditti V, Angers CG, Cornilescu G, Butcher SE: Structure and thermodynamics of a conserved U2 snRNA domain from yeast and human. RNA. 2007, 13 (3): 328-338. 10.1261/rna.418407.
Hermann T, Westhof E: Non-Watson-Crick base pairs in RNA-protein recognition. Chem Biol. 1999, 6 (12): R335-R343. 10.1016/S1074-5521(00)80003-4.
Simmonds P, Tuplin A, Evans DJ: Detection of genome-scale ordered RNA structure (GORS) in genomes of positive-stranded RNA viruses: implications for virus evolution and host persistence. RNA. 2004, 10 (9): 1337-1351. 10.1261/rna.7640104.
Priore SF, Moss WN, Turner DH: influenza A virus coding regions exhibit host-specific global ordered RNA structure. PloS One. 2012, 7 (4): e35989-10.1371/journal.pone.0035989.
Cuceanu NM, Tuplin A, Simmonds P: Evolutionarily conserved RNA secondary structures in coding and non-coding sequences at the 3′ end of the hepatitis G virus/GB-virus C genome. J Gen Virol. 2001, 82 (Pt 4): 713-722.
Clyde K, Harris E: RNA secondary structure in the coding region of dengue virus type 2 directs translation start codon selection and is required for viral replication. J Virol. 2006, 80 (5): 2170-2182. 10.1128/JVI.80.5.2170-2182.2006.
Mao Y, Li Q, Wang W, Liang P, Tao S: Number variation of high stability regions is correlated with gene functions. Genome Biol Evol. 2013, 5 (3): 484-493. 10.1093/gbe/evt020.
Wen JD, Lancaster L, Hodges C, Zeri AC, Yoshimura SH, Noller HF, Bustamante C, Tinoco I: Following translation by single ribosomes one codon at a time. Nature. 2008, 452 (7187): 598-603. 10.1038/nature06716.
Cabrita LD, Dobson CM, Christodoulou J: Protein folding on the ribosome. Curr Opin Struct Biol. 2010, 20 (1): 33-45. 10.1016/j.sbi.2010.01.005.
Siller E, DeZwaan DC, Anderson JF, Freeman BC, Barral JM: Slowing bacterial translation speed enhances eukaryotic protein folding efficiency. J Mol Biol. 2010, 396 (5): 1310-1318. 10.1016/j.jmb.2009.12.042.
Fukuda M, Longnecker R: Epstein-Barr virus latent membrane protein 2A mediates transformation through constitutive activation of the Ras/PI3-K/Akt pathway. J Virol. 2007, 81 (17): 9299-9306. 10.1128/JVI.00537-07.
Firth AE, Wills NM, Gesteland RF, Atkins JF: Stimulation of stop codon readthrough: frequent presence of an extended 3′ RNA structural element. Nucleic Acids Res. 2011, 39 (15): 6679-6691. 10.1093/nar/gkr224.
Gu W, Zhou T, Wilke CO: A universal trend of reduced mRNA stability near the translation-initiation site in prokaryotes and eukaryotes. PLoS Comput Biol. 2010, 6 (2): e1000664-10.1371/journal.pcbi.1000664.
Das A, Bhadra K, Suresh Kumar G: Targeting RNA by small molecules: comparative structural and thermodynamic aspects of aristololactam-beta-D-glucoside and daunomycin binding to tRNA (phe). PloS One. 2011, 6 (8): e23186-10.1371/journal.pone.0023186.
Stelzer AC, Frank AT, Kratz JD, Swanson MD, Gonzalez-Hernandez MJ, Lee J, Andricioaei I, Markovitz DM, Al-Hashimi HM: Discovery of selective bioactive small molecules by targeting an RNA dynamic ensemble. Nat Chem Biol. 2011, 7 (8): 553-559. 10.1038/nchembio.596.
Pushechnikov A, Lee MM, Childs-Disney JL, Sobczak K, French JM, Thornton CA, Disney MD: Rational design of ligands targeting triplet repeating transcripts that cause RNA dominant disease: application to myotonic muscular dystrophy type 1 and spinocerebellar ataxia type 3. J Am Chem Soc. 2009, 131 (28): 9767-9779. 10.1021/ja9020149.
Childs JL, Disney MD, Turner DH: Oligonucleotide directed misfolding of RNA inhibits Candida albicans group I intron splicing. Proc Natl Acad Sci USA. 2002, 99 (17): 11091-11096. 10.1073/pnas.172391199.
Childs-Disney JL, Wu M, Pushechnikov A, Aminova O, Disney MD: A small molecule microarray platform to select RNA internal loop-ligand interactions. ACS Chem Biol. 2007, 2 (11): 745-754. 10.1021/cb700174r.
Havens MA, Duelli DM, Hastings ML: Targeting RNA splicing for disease therapy. WIREs RNA. 2013, 4: 247-266. 10.1002/wrna.1158.
Zhou J, Zheng X, Shen H: Targeting RNA-splicing for SMA treatment. Mol Cells. 2012, 33 (3): 223-228. 10.1007/s10059-012-0005-6.
Reiche K, Stadler PF: RNAstrand: reading direction of structured RNAs in multiple sequence alignments. Algor Mol Biol: AMB. 2007, 2: 6-10.1186/1748-7188-2-6.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Bernhart SH, Hofacker IL, Will S, Gruber AR, Stadler PF: RNAalifold: improved consensus structure prediction for RNA alignments. BMC Bioinformatics. 2008, 9: 474-10.1186/1471-2105-9-474.
Harmanci AO, Sharma G, Mathews DH: TurboFold: iterative probabilistic estimation of secondary structures for multiple RNA sequences. BMC Bioinformatics. 2011, 12: 108-10.1186/1471-2105-12-108.
Han K, Lee Y, Kim W: PseudoViewer: automatic visualization of RNA pseudoknots. Bioinformatics. 2002, 18 (Suppl 1): S321-S328. 10.1093/bioinformatics/18.suppl_1.S321.
Schneider CA, Rasband WS, Eliceiri KW: NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012, 9 (7): 671-675. 10.1038/nmeth.2089.
Girish V, Vijayalakshmi A: Affordable image analysis using NIH Image/ImageJ. Indian J Cancer. 2004, 41 (1): 47-
Turner DH, Mathews DH: NNDB: the nearest neighbor parameter database for predicting stability of nucleic acid secondary structure. Nucleic Acids Res. 2010, 38 (Database issue): D280-D282.
Mathews DH, Disney MD, Childs JL, Schroeder SJ, Zuker M, Turner DH: Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. Proc Natl Acad Sci USA. 2004, 101 (19): 7287-7292. 10.1073/pnas.0401799101.
Reuter JS, Mathews DH: RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics. 2010, 11: 129-10.1186/1471-2105-11-129.
Simmonds P: SSE: a nucleotide and amino acid sequence analysis platform. BMC Res Notes. 2012, 5: 50-10.1186/1756-0500-5-50.
Gruber AR, Lorenz R, Bernhart SH, Neubock R, Hofacker IL: The Vienna RNA websuite. Nucleic Acids Res. 2008, 36 (Web Server issue): W70-W74.
Rio DC, Ares M, Hannon GJ, Nilsen TW: Preparation of cytoplasmic and nuclear RNA from tissue culture cells. Cold Spring Harb Protoc. 2010, 2010 (6): 10.1101/pdb.prot5441.
Langmead B, Salzberg SL: Fast gapped-read alignment with bowtie 2. Nat Methods. 2012, 9 (4): 357-359. 10.1038/nmeth.1923.
We thank Dr. Genaro Pimienta for advice on RNA-Seq data analysis; and Drs. Anna Vilborg and Nara Lee for help with experimental analyses of ebv-sisRNA-1. We also thank Professor Douglas H. Turner (University of Rochester) and Dr. Kazimierz Tycowski for their helpful comments on the manuscript. This work was supported by grant CA16038 from the NIH. JAS is an investigator of the Howard Hughes Medical Institute.
The authors declare that they have no competing interests.
WNM: design and execution of experiments, data analysis/interpretation, and manuscript preparation. JAS: data analysis/interpretation, manuscript preparation, and final approval of the version to be published. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Table showing RNAz-predicted regions with likely stable and conserved RNA structure. Annotated regions (e.g. splice sites and miRNAs) are highlighted in different colors, which match color-annotated structure models provided in Additional file 3 with start and stop codons in green and red, respectively; splice sites in blue; miRNAs in orange; the W repeat hairpin in yellow; EBER1 in purple; and the v-snoRNA in brown. The columns for the 5′ and 3′ nt are the first and last nt of the EBV genome (NC_009334.1) that comprise the predicted structured region. The 5′ and 3′ nt-Rev numbers indicate the location of structure that likely occurs in RNA transcribed from the reverse genome sequence (corresponding to the structure models shown in the REV_STR worksheet in Additional file 3). (XLSX 31 KB)
Additional file 3:Spreadsheet showing structure models for predicted EBV structured regions. The sequence shown is for the EBV-2 strain (NC_009334.1). In EBV, transcription can occur from either DNA strand and results for the forward and reverse sense genome sequences are presented in separate worksheets (FWD_STR and REV_STR). RNA structures are in “dot-bracket” notation, with paired sites indicated with matched brackets and unpaired with dots. Colored regions correspond to annotated RNAs in Additional file 1. (XLSX 4 MB)
Additional file 4:Structure models and predictions of possible A-to-I editing sites (using the Inosine Predict program) for lymphocryptovirus repeat long hairpin RNAs. Each virus has a separate listing with nt position numbered for just the hairpin, followed by the sequence, the secondary structure in dot-bracket notation, the calculated percent editing for four ADAR specificities, and the maximum predicted value of all ADAR specificities. A reported A-to-I editing site  is colored red. (XLSX 131 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Moss, W.N., Steitz, J.A. Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA. BMC Genomics 14, 543 (2013). https://doi.org/10.1186/1471-2164-14-543
- Epstein-Barr virus (EBV)
- RNA structure
- Non-coding RNA (ncRNA)
- W repeat
- RNA editing