Human regulatory proteins associate with non-coding RNAs from the EBV IR1 region
BMC Research Notes volume 11, Article number: 139 (2018)
The function of Epstein–Barr virus (EBV) stable intronic sequence (sis)RNAs, non-coding RNAs transcribed from a region required for EBV-mediated cellular transformation, remain unknown. To better understand the function of ebv-sisRNA-1 and ebv-sisRNA-2 from the internal repeat (IR)1 region of EBV, we used a combination of bioinformatics and biochemistry to identify associated RNA binding proteins. The findings reported here are part of ongoing studies to determine the functions of non-coding RNAs from the IR1 region of EBV.
Human regulatory proteins HNRNPA1 (heterogeneous nuclear ribonucleoprotein A1), HNRNPC, HNRNPL, HuR (human antigen R), and protein LIN28A (lin-28 homolog A) were predicted to bind ebv-sisRNA-1 and/or ebv-sisRNA-2; FUS (fused in sarcoma) was predicted to associate with ebv-sisRNA-2. Protein interactions were validated using a combination of RNA immunoprecipitation and biotin pulldown assays. Both sisRNAs also precipitated with HNRNPD and NONO (non-POU domain-containing octamer-binding protein). Interestingly, each of these interacting proteins also precipitated non-spliced non-coding RNA sequences transcribed from the IR1 region. Our findings suggest interesting roles for sisRNAs (through their interactions with regulatory proteins) and provide further evidence for the existence of non-spliced stable non-coding RNAs.
Stable intronic sequences (sis)RNAs and long non-coding (lnc)RNAs are non-coding RNAs (nc)RNAs known to perform many functions that can regulate gene expression [1, 2]. Aberrant expression of lncRNAs is thought to contribute to disease states such as cancer . The impact of the relatively new category of sisRNAs are not well understood but could include roles in development or cancer [1, 4]. Rarely (if ever), is the function of these RNAs independent of the proteins that bind them. Numerous RNA binding proteins have been identified and many of their important functions have been discussed elsewhere in detail , including their implication in disease states such as cancer .
Two abundant sisRNAs expressed from the internal repeat (IR)1 region, containing the W-repeats, of the Epstein–Barr virus (EBV; human herpes virus 4, HHV4) genome were identified previously . Their transcription is linked to that of the Epstein–Barr nuclear antigen leader protein (EBNA-LP), and occurs during cancer-associated latency III infection as well as in a rare latency program found in Burkitts lymphomas: Wp-restricted . Efficient transformation of B cells requires EBNA-LP expression . Indeed, the connection between EBV and cancer is well established [10, 11], but the role of the IR1 region stable introns in cellular transformation is not. Recently, a mutation in ebv-sisRNA-1 that inhibited EBV-mediated transformation of B cells was reported . Though the mutation did not occur in isolation, this maintains the idea that the ncRNA from this region is important for tumorigenesis. Here, we extend the understanding of the non-coding RNAs, ebv-sisRNA-1 and ebv-sisRNA-2, through the demonstration of their association with a number of human regulatory proteins. Further, we provide evidence that these intronic sequences are more likely to be found embedded within non-spliced transcripts that retain both exons of the IR1-region.
Descriptions of cell lines and transductions, biotin precipitations and mass spectrometry, bioinformatics and data processing, RNA immunoprecipitation, RNA extraction, and PCR are in Additional file 1.
Results and discussion
Initial identification and validation of protein interactors of ebv-sisRNA-1
The majority of ebv-sisRNA-1 is predicted to be unstructured (Fig. 1a), providing easier access to single-stranded RNA binding proteins. To identify proteins that bind ebv-sisRNA-1, biotin precipitation assays were performed using BJAB cell lysates (EBV negative) and either a wild-type or deletion mutant of the conserved CA-rich region (ΔCA). This region was speculated previously to be important for heterogeneous nuclear ribonucleoprotein (HNRNP)L binding . Silver staining revealed losses of several protein bands, prominently at both 60–65 kDa (size of HNRNPL) and 35–40 kDa (most intense band for wt) compared to ΔCA mutation (Fig. 1b). Mass spectrometry of these regions identified HNRNPL and HNRNPD/AUF1 as the most probable interactors (Additional file 2). HNRNPL binding is fully consistent with reported CA-rich consensus binding sequences [13, 14]. HNRNPD/AUF1 typically binds AU-rich elements (AREs) of RNA  but has also been shown to interact with HNRNPL . RNA immunoprecipitation (RIP) for HNRNPL and HNRNPD/AUF1 on both Raji and BJAB-B1 cells, which harbor type 1 and type 2 EBV genomes, respectively, revealed an association with ebv-sisRNA-1 (Fig. 1c; Additional file 3). Ebv-sisRNA-1 is highly abundant and expressed at the same time as EBNA-LP , and recent genetic evidence suggests that EBNA-LP and possibly ebv-sisRNA-1 may be important for EBV-mediated B cell transformation ; thus, its interactions with two key regulatory proteins have interesting implications.
HNRNPD/AUF1 is best known for its function in regulating mRNA abundance (e.g., binding to ARE-containing 3′ untranslated regions (UTRs) and destabilizing transcripts) . Overexpression of HNRNPD/AUF1 can promote oncogenesis and control translation [18, 19]. Interestingly, HNRNPD/AUF1 binds directly to EBER1 , another EBV ncRNA that is highly expressed throughout all stages of latent/lytic infection. The proposed function of this EBER1 interaction is to perturb the normal homeostasis between HNRNPD/AUF1 and target mRNAs in EBV-infected cells. Given that both ebv-sisRNA-1 and EBER1 are abundantly present during the initial patterning of B-cell gene expression (e.g., latency III), it is enticing to speculate that ebv-sisRNA-1 and EBER1 may coordinately contribute to disease-promoting functions of HNRNPD/AUF1.
HNRNPL has a proposed role in various cancers and is known to affect the expression (e.g., alternative splicing) of a wide array of human genes [21,22,23,24]. HNRNPL shuttles between the cytoplasm and nucleus; thus, its interaction with a viral ncRNA that is abundantly present in the nucleus, suggests a potential role in disrupting HNRNPL localization/interactions. This could have a vast impact on the alternative splicing of many host genes; as well as, interestingly, in the regulation of the EBNA-LP gene from which ebv-sisRNA-1 is derived: e.g., through a competition between EBNA-LP pre-mRNA and the accumulating sisRNA that forms stable ribonucleoproteins (RNPs) with HNRNPL. These possibilities are being explored.
Identification of additional interactors of ebv-sisRNA-1 and -2
To better understand the protein binding repertoire of ebv-sisRNA-1 and ebv-sisRNA-2, we used the program RBPmap  to predict human host RNA binding proteins (RBPs) that might directly interact. Both type 1 (NC_007605) and type 2 (NC_009334) EBV sequences were analyzed using a sequence fragment of the IR1 W-repeat region from each genome, which contains exon W1, ebv-sisRNA-1, exon W2, and ebv-sisRNA-2 (Fig. 2a, W1s1W2s2). RBPmap compares these sequences to known consensus binding sites of human and mouse proteins and classifies/ranks hits based on their fit to the consensus motif. 81 unique proteins were identified that had significant binding site consensus throughout the region of interest (Fig. 2b; see Additional file 4 for binding region detailed information). The putative interacting proteins are known to be involved in biological processes of RNA processing, splicing, and gene expression (Fig. 2c; Additional files 3 and 5).
Consistent with our findings above, ebv-sisRNA-1 contained predicted HNRNPL binding regions (Fig. 2b). Ebv-sisRNA-2 did as well, and this interaction was verified by RIP (Fig. 2d). Similar to ebv-sisRNA-1, ebv-sisRNA-2 precipitated with HNRNPD/AUF1. Multiple interaction sites for several RBPs throughout the sisRNAs suggest similar functions for potential regulation of EBNA-LP and sisRNA processing, as well as downstream effects on cellular homeostasis of RBP interactions with endogenously regulated host genes, which remains to be explored.
Validation of select interactors
Eight other predicted proteins (Fig. 2b, blue; Additional file 6) were tested for their ability to co-precipitate either ebv-sisRNA-1 (HuR/ELAVL1, IGF2BP2, IGF2BP3), ebv-sisRNA-2 (FUS, HNRNPA1, LIN28, MBNL1), or both (HNRNPC). Quantitative (q)PCR for both sisRNAs was performed for all RIPs (Fig. 3a; Additional file 3). HuR, HNRNPA1, LIN28, and HNRNPC precipitated both ebv-sisRNAs with at least twofold enrichment over normal IgG in two cell lines. FUS bound ebv-sisRNA-2 in both cell lines, but only precipitated ebv-sisRNA-1 in BJAB-B1 cells. IGF2BP2, IGF2BP3, and MBNL1 did not bind under these experimental conditions.
Two RBPs not predicted by RBPmap were also tested. To determine if sisRNAs were potential targets of micro (mi)RNAs (e.g., acting as “sponges”), an antibody that targets multiple Argonaute proteins was used. Neither sisRNA was enriched in our experiments. This result challenges a previous in silico predicted interaction  of ebv-sisRNA-1 and miR-142-3p (previously proposed to be targeted by another herpes virus ncRNA [26, 27]). Lytic EBV ncRNAs originating from the OriP region contain long hairpin stems similar in size/structure to ebv-sisRNA-2 [7, 26] and were found to interact with the paraspeckle-associated protein p54nrb/NONO . The OriP ncRNA was found to promote transcription/replication of EBV and to modulate host innate immunity through its interactions with paraspeckles. Likewise, NONO interacts with EBER2  and may bridge an interaction with host transcription factor PAX5, which supports EBV replication through recruitment by EBER2 to the viral genome . Similar to the OriP and EBER2 ncRNAs, NONO bound both ebv-sisRNAs by RIP. The common interaction of these viral ncRNAs with NONO suggests potential overlapping functions (e.g., replication or immune modulation) that need to be tested. Additionally, localization of the sisRNAs (and other viral ncRNAs) to paraspeckles remains to be shown.
These results show binding—either direct or indirect—of the IR1 sisRNAs to HNRNPL, HNRNPD/AUF1, HuR, HNRNPA1, LIN28, NONO, and HNRNPC (Fig. 3a). Intriguingly, each of these proteins are implicated in pro-cancer-related phenotypes [6, 31]. Binding of these regulatory proteins to ebv-sisRNAs may play roles in the emergence of these phenotypes in EBV-associated diseases, which merits further investigation.
Identification of a junction-spanning IR1-derived ncRNA that co-precipitates with RBPs
Intrigued that both ebv-sisRNAs associated with HuR, HNRNPA1, FUS, and LIN28, when only one or the other was predicted to bind, a primer set that spanned the ebv-sisRNA-1 to W2 splice junction was tested (Figs. 2b, amplicon 3; 3b; Additional file 3). This RNA, containing an intact splice junction, precipitated with the majority of interacting proteins, suggesting that these RBPs also bound non-spliced RNA from the W-repeat region. To test for the relative abundance of spliced to non-spliced ebv-sisRNA-1, RT-PCR was performed using a forward primer targeting exon W1 and a reverse primer targeting exon W2 (Figs. 2b, amplicon 4; 3c). The overwhelming majority of RNA contained intact exon–intron junctions. The non-spliced form was also more abundant than the spliced form in EBV-free BJAB cells transduced to overexpress W1-sis1-W2. Likewise, amplification from the 3′ end of sisRNA-2 across W1-sisRNA-1-W2 to the 5′ end of the next sisRNA-2 intron (Fig. 2b, amplicon 5) demonstrated that the predominant form of RNA from the IR1 W-repeats retains intact exon–intron boundaries (Fig. 3d).
There are several possible interpretations of these findings. Firstly, W-repeat intronic RNAs could be embedded within non-polyadenylated [poly(A)] lncRNAs arising from this region. Secondly, splicing may have resulted in excision of an entire stable sis2-exonW1-sis1-exonW2 region (~ 5.9 kb) or an exonW1-sis1-exonW2-sis2 (~ 3 kb) region; the stability of such products could possibly be conferred through circularization and back splicing . Both interpretations are made plausible by early work that found prominent Northern blot bands of poly(A)− RNA up to 13 kb in size. . Thirdly, splicing may be inhibited leading to the accumulation of intron-containing pre-mRNA from the EBV W-repeat region, which has been observed previously . HNRNPL and HNRNPA1 are known spicing repressors , which supports this hypothesis. Likewise, lncRNAs spanning the W-repeat intronic regions may be polyadenylated. Previous analyses of poly(A)+ samples [33, 36], showed Northern blot bands overlapping poly(A)− RNAs as well as larger ones (up to ~ 22 kb). These possibilities are not mutually exclusive and the exact architecture of the ncRNAs from this region is likely to be complex: consisting of multiple sized species of spliced/unspliced RNAs with and without poly(A) tails. Supporting this idea, recent analyses of RNA-Seq data from EBV-infected cells found extensive read coverage of the ebv-sisRNA regions in both ribosome-depleted and poly(A)-selected samples .
We have shown that the ebv-sisRNAs exist within a complex milieu of non-coding transcripts. In addition to independent intronic units, sisRNAs are embedded within various sized ncRNAs that can arise from alternative splicing or initiation and that may be polyadenylated. Regardless of their transcript structure, sisRNAs are capable of interacting with a number of human regulatory proteins. We present a list of 81 proteins from in silico predictions and validate six highly interesting ones that associate with both ebv-sisRNA-1 and ebv-sisRNA-2 (FUS, HNRNPA1, HNRNPC, HNRNPL, HuR, LIN28). Furthermore, we identify additional interactions with HNRNPD/AUF1 and NONO, both of which interact with other EBV ncRNAs. Taken together, our results begin to unravel the RNP composition of this intriguing class of viral ncRNAs and suggest multiple leads for follow-up analyses into the functional consequences of deduced interactions: e.g., in EBV infection and the emergence of disease.
Additional work is needed to confirm other predicted interactions as well as to determine if validated binding proteins are interacting directly or indirectly. Most importantly, the exciting implications of the identified/validated RBPs have yet to be explored. We hope this report will facilitate additional research required to understand the functions of these interesting EBV ncRNAs.
Epstein–Barr virus nuclear antigen leader protein
false discovery rate
human herpesvirus 4
internal repeat 1
long non-coding RNA
quantitative polymerase chain reaction
RNA binding protein
reverse-transcriptase polymerase chain reaction
stable intronic sequence RNA
Tay MLI, Pek JW. Maternally inherited stable intronic sequence rna triggers a self-reinforcing feedback loop during development. Curr Biol. 2017;27:1062–7.
Marchese FP, Raimondi I, Huarte M. The multidimensional mechanisms of long noncoding rna function. Genome Biol. 2017;18:206.
Bhan A, Soleimani M, Mandal SS. Long noncoding rna and cancer: a new paradigm. Cancer Res. 2017;77:3965–81.
Osman I, Tay MLI, Pek JW. Stable intronic sequence rnas (sisrnas): a new layer of gene regulation. Cell Mol Life Sci. 2016;73:3507–19.
Gerstberger S, Hafner M, Ascano M, Tuschl T. Evolutionary conservation and expression of human rna-binding proteins and their role in human genetic disease. In: Yeo GW, editor. Systems biology of RNA binding proteins. New York: Springer; 2014. p. 1–55.
Pereira B, Billaud M, Almeida R. Rna-binding proteins in cancer: old players and new actors. Trends Cancer. 2017;3:506–28.
Moss WN, Steitz JA. Genome-wide analyses of epstein-barr virus reveal conserved rna structures and a novel stable intronic sequence rna. BMC Genomics. 2013;14:543.
Kelly GL, Milner AE, Tierney RJ, Croom-Carter DSG, Altmann M, Hammerschmidt W, et al. Epstein-barr virus nuclear antigen 2 (ebna2) gene deletion is consistently linked with ebna3a, -3b, and -3c expression in burkitt’s lymphoma cells and with increased resistance to apoptosis. J Virol. 2005;79:10709–17.
Mannick JB, Cohen JI, Birkenbach M, Marchini A, Kieff E. The epstein-barr virus nuclear protein encoded by the leader of the ebna rnas is important in b-lymphocyte transformation. J Virol. 1991;65:6826–37.
Thompson MP, Kurzrock R. Epstein-barr virus and cancer. Clin Cancer Res. 2004;10:803–21.
de Oliveira ED, Müller-Coan BG, Pagano JS. Viral carcinogenesis beyond malignant transformation: Ebv in the progression of human cancers. Trends Microbiol. 2016;24:649–64.
Szymula A, Palermo RD, Groves IJ, Ba abdullah M, Holder BS, White RE. Epstein-barr virus nuclear antigen ebna-lp is essential for transforming naive b cells, and facilitates recruitment of transcription factors to the viral genome. bioRxiv. 2017;176099. https://doi.org/10.1101/176099.
Hui J, Hung LH, Heiner M, Schreiner S, Neumüller N, Reither G, et al. Intronic ca-repeat and ca-rich elements: a new class of regulators of mammalian alternative splicing. EMBO J. 2005;24:1988–98.
Shankarling G, Cole BS, Mallory MJ, Lynch KW. Transcriptome-wide rna interaction profiling reveals physical and functional targets of hnrnp l in human t cells. Mol Cell Biol. 2014;34:71–83.
Yoon J-H, De S, Srikantan S, Abdelmohsen K, Grammatikakis I, Kim J, et al. Par-clip analysis uncovers auf1 impact on target rna fate and genome integrity. Nat Commun. 2014;5:5248.
Park HG, Yoon JY, Choi M. Heterogeneous nuclear ribonucleoprotein d/auf1 interacts with heterogeneous nuclear ribonucleoprotein l. J Biosci. 2007;32:1263–72.
Brewer G. An a + u-rich element rna-binding factor regulates c-myc mrna stability in vitro. Mol Cell Biol. 1991;11:2460–6.
Gouble A, Grazide S, Meggetto F, Mercier P, Delsol G, Morello D. A new player in oncogenesis: Auf1/hnrnpd overexpression leads to tumorigenesis in transgenic mice. Cancer Res. 2002;62:1489–95.
Liao B, Hu Y, Brewer G. Competitive binding of auf1 and tiar to myc mrna controls its translation. Nat Struct Mol Biol. 2007;14:511.
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:2073–82.
Goehe RW, Shultz JC, Murudkar C, Usanovic S, Lamour NF, Massey DH, et al. Hnrnp l regulates the tumorigenic capacity of lung cancer xenografts in mice via caspase-9 pre-mrna processing. J Clin Invest. 2010;120:3923–39.
Fei T, Chen Y, Xiao T, Li W, Cato L, Zhang P, et al. Genome-wide crispr screen identifies hnrnpl as a prostate cancer dependency regulating rna splicing. Proc Natl Acad Sci USA. 2017;114:E5207–15.
Niknafs YS, Han S, Ma T, Speers C, Zhang C, Wilder-Romans K, et al. The lncrna landscape of breast cancer reveals a role for dscam-as1 in breast cancer progression. Nat Commun. 2016;7:12791.
Rothrock CR, House AE, Lynch KW. Hnrnp l represses exon splicing via a regulated exonic splicing silencer. EMBO J. 2005;24:2792–802.
Paz I, Kosti I, Ares JM, Cline M, Mandel-Gutfreund Y. Rbpmap: a web server for mapping binding sites of rna-binding proteins. Nucleic Acids Res. 2014;42:W361–7.
Moss W. Analyses of non-coding rnas generated from the epstein–barr virus w repeat region. Proc IWBBIO. 2014;2014:238–52.
Cazalla D, Yario T, Steitz JA. Down-regulation of a host microrna by a Herpesvirus saimiri noncoding rna. Science. 2010;328:1563–6.
Cao S, Moss W, O’Grady T, Concha M, Strong MJ, Wang X, et al. 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. J Virol. 2015;89:7120–32.
Lee N, Yario TA, Gao JS, Steitz JA. Ebv noncoding rna eber2 interacts with host rna-binding proteins to regulate viral gene expression. Proc Natl Acad Sci USA. 2016;113:3221–6.
Lee N, Moss Walter N, Yario Therese A, Steitz Joan A. Ebv noncoding rna binds nascent rna to drive host pax5 to viral DNA. Cell. 2015;160:607–18.
David CJ, Manley JL. Alternative pre-mrna splicing regulation in cancer: pathways and programs unhinged. Genes Dev. 2010;24:2343–64.
Chen L-L. The biogenesis and emerging roles of circular rnas. Nat Rev Mol Cell Biol. 2016;17:205.
van Santen V, Cheung A, Hummel M, Kieff E. Rna encoded by the ir1-u2 region of epstein-barr virus DNA in latently infected, growth-transformed cells. J Virol. 1983;46:424–33.
Melčák I, Cermanová Š, Jirsová K, Koberna K, Malı́nský J, Raška I. Nuclear pre-mrna compartmentalization: trafficking of released transcripts to splicing factor reservoirs. Mol Biol Cell. 2000;11:497–510.
Fu X-D, Ares M Jr. Context-dependent control of alternative splicing by rna-binding proteins. Nat Rev Genet. 2014;15:689.
Rogers RP, Speck SH. Bidirectional transcription of the epstein-barr virus major internal repeat. J Virol. 1990;64:2426–9.
Cao S, Strong MJ, Wang X, Moss WN, Concha M, Lin Z, et al. High-throughput rna sequencing-based virome analysis of 50 lymphoma cell lines from the cancer cell line encyclopedia project. J Virol. 2015;89:713–29.
Snel B, Lehmann G, Bork P, Huynen MA. String: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000;28:3442–4.
VST—experimental design and execution, interpretation of results and primary authorship of the manuscript; DV executed experiments; WNM conceived of the work, oversaw its production, and contributed to the writing of the manuscript. All authors read and approved the final manuscript.
We would like to thank Ashwanth Lalithaa-Padmanabhan (Iowa State University) and Nuwanthika D. Kumarashinghe (Iowa State University) for technical contributions, as well as Jacob Schwartz (University of Arizona), Joan Steitz (Yale University), and Siegfried Janz (University of Iowa) for kindly providing antibodies and cells.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Consent for publication
Ethics approval and consent to participate
This work was supported by startup funds from the Iowa State University College of Agriculture and Life Sciences and the Roy J. Carver Charitable Trust, as well as Grant 4R00GM112877-02 from the NIH/NIGMS.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Methods and materials. Detailed description of methods and materials used to collect and analyze data.
Mass spectrometry results. List of proteins that bound biotinylated forms of ebv-sisRNA-1 after mass spectrometry from gel slices indicated in Fig. 1B. Hits resulting from either the top (60-65 kDa) or bottom (35-40 kDa) are indicated.
Fold enrichment individual experimental results from all RIPs.
RBPmap of EBV W1s1W2s2. Human or mouse proteins predicted (RBPmap) to bind W1s1W2s2 RNA (see Fig. 2A) from the IR1 region of EBV. Results from both EBV type 1 and type 2 sequences are included as well as a comparison.
STRING analyses of predicted RBPs. Results from STRING enrichment analyses of 81 RBPmap-identified proteins. Information includes annotations and enrichment data on biological processes, KEGG pathways, molecular functions, cellular components, InterPro, and Pfam.
Map of predicted binding sites for proteins assessed by RIP.
About this article
Cite this article
Tompkins, V.S., Valverde, D.P. & Moss, W.N. Human regulatory proteins associate with non-coding RNAs from the EBV IR1 region. BMC Res Notes 11, 139 (2018). https://doi.org/10.1186/s13104-018-3250-8