The Epstein Barr virus circRNAome

Characterization of the SARS-CoV-2 Genome 3'-Untranslated Region Interactions with Host MicroRNAs

The 2019 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has marked the spread of a novel human coronavirus. SARS-CoV-2 has exhibited increased disease severity and immune evasion across its variants, and the molecular mechanisms behind these phenomena remain largely unknown. Conserved elements of the viral genome, such as secondary structures within the 3'-untranslated region (UTR), could prove crucial in furthering our understanding of the host-virus interface. Analysis of the SARS-CoV-2 viral genome 3'-UTR revealed the potential for host microRNA (miR) binding sites, allowing for sequence-specific interactions. In this study, we demonstrate that the SARS-CoV-2 genome 3'-UTR binds the host cellular miRs miR-34a-5p, miR-34b-5p, and miR-760-3p in vitro. Native gel electrophoresis and steady-state fluorescence spectroscopy were utilized to biophysically characterize the binding of these miRs to their predicted sites within the SARS-CoV-2 genome 3'-UTR. Additionally, we investigated 2'-fluoro-d-arabinonucleic acid (FANA) analogs as competitive binding inhibitors for these interactions. These miRs modulate the translation of granulin (GRN), interleukin-6 (IL-6), and the IL-6 receptor (IL-6R), all of which are key modulators and activators of JAK/STAT3 signaling and are implicated in regulation of the immune response. Thus, we propose that hijacking of these miRs by SARS-CoV-2 could identify a mechanism of host immune modulation by the virus. The mechanisms detailed in this study have the potential to drive the development of antiviral treatments for SARS-CoV-2, through direct targeting of the virus-host interface.

HIV-associated cancers and lymphoproliferative disorders caused by Kaposi sarcoma herpesvirus and Epstein-Barr virus

SUMMARYWithin weeks of the first report of acquired immunodeficiency syndrome (AIDS) in 1981, it was observed that these patients often had Kaposi sarcoma (KS), a hitherto rarely seen skin tumor in the USA. It soon became apparent that AIDS was also associated with an increased incidence of high-grade lymphomas caused by Epstein-Barr virus (EBV). The association of AIDS with KS remained a mystery for more than a decade until Kaposi sarcoma-associated herpesvirus (KSHV) was discovered and found to be the cause of KS. KSHV was subsequently found to cause several other diseases associated with AIDS and human immunodeficiency virus (HIV) infection. People living with HIV/AIDS continue to have an increased incidence of certain cancers, and many of these cancers are caused by EBV and/or KSHV. In this review, we discuss the epidemiology, virology, pathogenesis, clinical manifestations, and treatment of cancers caused by EBV and KSHV in persons living with HIV.

Identification and characterization of Varicella Zoster Virus circular RNA in lytic infection

This study investigates the role of circular RNAs (circRNAs) in the context of Varicella-Zoster Virus (VZV) lytic infection. We employ two sequencing technologies, short-read sequencing and long-read sequencing, following RNase R treatment on VZV-infected neuroblastoma cells to identify and characterize both cellular and viral circRNAs. Our large scanning analysis identifies and subsequent experiments confirm 200 VZV circRNAs. Moreover, we discover numerous VZV latency-associated transcripts (VLTs)-like circRNAs (circVLTslytic), which contain multiple exons and different isoforms within the same back-splicing breakpoint. To understand the functional significance of these circVLTslytic, we utilize the Bacteria Artificial Chromosome system to disrupt the expression of viral circRNAs in genomic DNA location. We reveal that the sequence flanking circVLTs' 5' splice donor plays a pivotal role as a cis-acting element in the formation of circVLTslytic. The circVLTslytic is dispensable for VZV replication, but the mutation downstream of circVLTslytic exon 5 leads to increased acyclovir sensitivity in VZV infection models. This suggests that circVLTslytic may have a role in modulating the sensitivity to antiviral treatment. The findings shed new insight into the regulation of cellular and viral transcription during VZV lytic infection, emphasizing the intricate interplay between circRNAs and viral processes.

Interferon induced circRNAs escape herpesvirus host shutoff and suppress lytic infection

To globally profile circRNAs, we employ RNA-Sequencing paired with chimeric junction analysis for alpha-, beta-, and gamma-herpesvirus infection. We find circRNAs are, as a population, resistant to host shutoff. We validate this observation using ectopic expression assays of human and murine herpesvirus endoribonucleases. During lytic infection, four circRNAs are commonly induced across all subfamilies of human herpesviruses, suggesting a shared mechanism of regulation. We test one such mechanism, namely how interferon-stimulation influences circRNA expression. 67 circRNAs are upregulated by either interferon-β or -γ treatment, with half of these also upregulated during lytic infection. Using gain and loss of function studies we find an interferon-stimulated circRNA, circRELL1, inhibits lytic Herpes Simplex Virus-1 infection. We previously reported circRELL1 inhibits lytic Kaposi sarcoma-associated herpesvirus infection, suggesting a pan-herpesvirus antiviral activity. We propose a two-pronged model in which interferon-stimulated genes may encode both mRNA and circRNA with antiviral activity. This is critical in cases of host shutoff, such as alpha- and gamma-herpesvirus infection, where the mRNA products are degraded but circRNAs escape.

Hijacking and rewiring of host CircRNA/miRNA/mRNA competitive endogenous RNA (ceRNA) regulatory networks by oncoviruses during development of viral cancers

A significant portion of human cancers are caused by oncoviruses (12%-25%). Oncoviruses employ various strategies to promote their replication and induce tumourigenesis in host cells, one of which involves modifying the gene expression patterns of the host cells, leading to the rewiring of genes and resulting in significant changes in cellular processes and signalling pathways. In recent studies, a specific mode of gene regulation known as circular RNA (circRNA)-mediated competing endogenous RNA (ceRNA) networks has emerged as a key player in this context. CircRNAs, a class of non-coding RNA molecules, can interact with other RNA molecules, such as mRNAs and microRNAs (miRNAs), through a process known as ceRNA crosstalk. This interaction occurs when circRNAs, acting as sponges, sequester miRNAs, thereby preventing them from binding to their target mRNAs and modulating their expression. By rewiring the host cell genome, oncoviruses have the ability to manipulate the expression and activity of circRNAs, thereby influencing the ceRNA networks that can profoundly impact cellular processes such as cell proliferation, differentiation, apoptosis, and immune responses. This review focuses on a comprehensive evaluation of the latest findings on the involvement of virus-induced reprogramming of host circRNA-mediated ceRNA networks in the development and pathophysiology of human viral cancers, including cervical cancer, gastric cancer, nasopharyngeal carcinoma, Kaposi's sarcoma, hepatocellular carcinoma, and diffuse large B cell lymphoma. Understanding these mechanisms can improve our knowledge of how oncoviruses contribute to human tumourigenesis and identify potential targets for developing optimised therapies and diagnostic tools for viral cancers.

Virus-Encoded Circular RNAs: Role and Significance in Viral Infections

Circular RNAs (circRNAs) have been the focus of intense scientific research to understand their biogenesis, mechanisms of action and regulatory functions. CircRNAs are single stranded, covalently closed RNA molecules lacking the 5'-terminal cap and the 3'-terminal polyadenine chain, characteristics that make them very stable and resistant. Synthesised by both cells and viruses, in the past circRNAs were considered to have no precise function. Today, increasing evidence shows that circRNAs are ubiquitous, some of them are tissue- and cell-specific, and critical in multiple regulatory processes (i.e., infections, inflammation, oncogenesis, gene expression). Moreover, circRNAs are emerging as important biomarkers of viral infection and disease progression. In this review, we provided an updated overview of current understanding of virus-encoded and cellular-encoded circRNAs and their involvement in cellular pathways during viral infection.

Critical involvement of circular RNAs in virus-associated cancers

Virus-related cancer is cancer where viral infection leads to the malignant transformation of the host's infected cells. Seven viruses (e.g., human papillomavirus (HPV), Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV), Hepatitis B virus (HBV), Hepatitis C virus (HCV), Human T-lymphotropic virus (HTLV), and Merkel cell polyomavirus (MCV)) that infect humans have been identified as an oncogene and have been associated with several human malignancies. Recently, growing attention has been attracted to exploring the pathogenesis of virus-related cancers. One of the most mysterious molecules involved in carcinogenesis and progression of virus-related cancers is circular RNAs (circRNA). These emerging non-coding RNAs (ncRNAs), due to the absence of 5' and 3' ends, have high stability than linear RNAs and are found in some species across the eukaryotic organisms. Compelling evidence has revealed that viruses also encode a repertoire of circRNAs, as well as dysregulation of these viral circRNAs play a critical role in the pathogenesis and progression of different types of virus-related cancers. Therefore, understanding the exact role and function of the virally encoded circRNAs with virus-associated cancers will open a new road for increasing our knowledge about the RNA world. Hence, in this review, we will focus on emerging roles of virus-encoded circRNAs in multiple cancers, including cervical cancer, gastric cancer, Merkel cell carcinoma, nasopharyngeal carcinoma, Kaposi cancer, and liver cancer.

Interferon induced circRNAs escape herpesvirus host shutoff and suppress lytic infection

A first line of defense during infection is expression of interferon (IFN)-stimulated gene products which suppress viral lytic infection. To combat this, herpesviruses express endoribonucleases to deplete host RNAs. Here we demonstrate that IFN-induced circular RNAs (circRNAs) can escape viral-mediated degradation. We performed comparative circRNA expression profiling for representative alpha- (Herpes simplex virus-1, HSV-1), beta- (human cytomegalovirus, HCMV), and gamma-herpesviruses (Kaposi sarcoma herpesvirus, KSHV; murine gamma-herpesvirus 68, MHV68). Strikingly, we found that circRNAs are, as a population, resistant to host shutoff. This observation was confirmed by ectopic expression assays of human and murine herpesvirus endoribonucleases. During primary lytic infection, ten circRNAs were commonly regulated across all subfamilies of human herpesviruses, suggesting a common mechanism of regulation. We tested one such mechanism, namely how interferon-stimulation influences circRNA expression. 67 circRNAs were upregulated by either IFN-β or -γ treatment, with half of these also upregulated during lytic infection. Using gain and loss of function studies we found an interferon-stimulated circRNA, circRELL1, inhibited lytic HSV-1 infection. We have previously reported circRELL1 inhibits lytic KSHV infection, suggesting a pan-herpesvirus antiviral activity. We propose a two-pronged model in which interferon-stimulated genes may encode both mRNA and circRNA with antiviral activity. This is critical in cases of host shutoff, such as alpha- and gamma-herpesvirus infection, where the mRNA products are degraded but circRNAs escape.

Insights into the Transcriptome of Human Cytomegalovirus: A Comprehensive Review

Human cytomegalovirus (HCMV) is a widespread pathogen that poses significant risks to immunocompromised individuals. Its genome spans over 230 kbp and potentially encodes over 200 open-reading frames. The HCMV transcriptome consists of various types of RNAs, including messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs), with emerging insights into their biological functions. HCMV mRNAs are involved in crucial viral processes, such as viral replication, transcription, and translation regulation, as well as immune modulation and other effects on host cells. Additionally, four lncRNAs (RNA1.2, RNA2.7, RNA4.9, and RNA5.0) have been identified in HCMV, which play important roles in lytic replication like bypassing acute antiviral responses, promoting cell movement and viral spread, and maintaining HCMV latency. CircRNAs have gained attention for their important and diverse biological functions, including association with different diseases, acting as microRNA sponges, regulating parental gene expression, and serving as translation templates. Remarkably, HCMV encodes miRNAs which play critical roles in silencing human genes and other functions. This review gives an overview of human cytomegalovirus and current research on the HCMV transcriptome during lytic and latent infection.

An Epstein-Barr virus protein interaction map reveals NLRP3 inflammasome evasion via MAVS UFMylation

Epstein-Barr virus (EBV) causes infectious mononucleosis, triggers multiple sclerosis, and is associated with 200,000 cancers/year. EBV colonizes the human B cell compartment and periodically reactivates, inducing expression of 80 viral proteins. However, much remains unknown about how EBV remodels host cells and dismantles key antiviral responses. We therefore created a map of EBV-host and EBV-EBV interactions in B cells undergoing EBV replication, uncovering conserved herpesvirus versus EBV-specific host cell targets. The EBV-encoded G-protein-coupled receptor BILF1 associated with MAVS and the UFM1 E3 ligase UFL1. Although UFMylation of 14-3-3 proteins drives RIG-I/MAVS signaling, BILF1-directed MAVS UFMylation instead triggered MAVS packaging into mitochondrial-derived vesicles and lysosomal proteolysis. In the absence of BILF1, EBV replication activated the NLRP3 inflammasome, which impaired viral replication and triggered pyroptosis. Our results provide a viral protein interaction network resource, reveal a UFM1-dependent pathway for selective degradation of mitochondrial cargo, and highlight BILF1 as a novel therapeutic target.

Functions of Circular RNA in Human Diseases and Illnesses

Circular RNAs (circRNAs) represent single-stranded RNA species that contain covalently closed 3' and 5' ends that provide them more stability than linear RNA, which has free ends. Emerging evidence indicates that circRNAs perform essential functions in many DNA viruses, including coronaviruses, Epstein-Barr viruses, cytomegalovirus, and Kaposi sarcoma viruses. Recent studies have confirmed that circRNAs are present in viruses, including DNA and RNA viruses, and play various important functions such as evading host immune response, disease pathogenesis, protein translation, miRNA sponges, regulating cell proliferation, and virus replication. Studies have confirmed that circRNAs can be biological signatures or pathological markers for autoimmune diseases, neurological diseases, and cancers. However, our understanding of circRNAs in DNA and RNA viruses is still limited, and functional evaluation of viral and host circRNAs is essential to completely understand their biological functions. In the present review, we describe the metabolism and cellular roles of circRNA, including its roles in various diseases and viral and cellular circRNA functions. Circular RNAs are found to interact with RNA, proteins, and DNA, and thus can modulate cellular processes, including translation, transcription, splicing, and other functions. Circular RNAs interfere with various signaling pathways and take part in vital functions in various biological, physiological, cellular, and pathophysiological processes. We also summarize recent evidence demonstrating cellular and viral circRNA's roles in DNA and RNA viruses in this growing field of research.

A review on circular RNAs and bacterial infections

Bacterial infections and related diseases have been a major burden on social public health and economic stability around the world. However, the effective diagnostic methods and therapeutic approaches to treat bacterial infections are still limited. As a group of non-coding RNA, circular RNAs (circRNAs) that were expressed specifically in host cells and played a key regulatory role have the potential to be of diagnostic and therapeutic value. In this review, we systematically summarize the role of circRNAs in common bacterial infections and their potential roles as diagnostic markers and therapeutic targets.

[The circular life of viruses]

Circular RNAs (circRNA), as ancient as the first viruses, take an important part in the host-pathogen relationship. After the first description of dysregulated cellular circRNAs upon viral infection, numerous circRNAs of viral origin were identified and characterized. They are impacting both viral and cellular cycles and are associated with virus-induced oncogenesis, immune system regulation and cell differentiation. While the naïve reader might get swamped by discovering this new field of RNA biology, it seems that these RNA rings are actually full of surprises and wonders at both a functional and a biogenesis level.

Virological Markers in Epstein–Barr Virus-Associated Diseases

Epstein–Barr virus (EBV) is an oncogenic virus infecting more than 95% of the world’s population. After primary infection—responsible for infectious mononucleosis in young adults—the virus persists lifelong in the infected host, especially in memory B cells. Viral persistence is usually without clinical consequences, although it can lead to EBV-associated cancers such as lymphoma or carcinoma. Recent reports also suggest a link between EBV infection and multiple sclerosis. In the absence of vaccines, research efforts have focused on virological markers applicable in clinical practice for the management of patients with EBV-associated diseases. Nasopharyngeal carcinoma is an EBV-associated malignancy for which serological and molecular markers are widely used in clinical practice. Measuring blood EBV DNA load is additionally, useful for preventing lymphoproliferative disorders in transplant patients, with this marker also being explored in various other EBV-associated lymphomas. New technologies based on next-generation sequencing offer the opportunity to explore other biomarkers such as the EBV DNA methylome, strain diversity, or viral miRNA. Here, we review the clinical utility of different virological markers in EBV-associated diseases. Indeed, evaluating existing or new markers in EBV-associated malignancies or immune-mediated inflammatory diseases triggered by EBV infection continues to be a challenge.

Cytoplasmic poly(A)-binding protein 1 (PABPC1) is a prognostic biomarker to predict survival in nasopharyngeal carcinoma regardless of chemoradiotherapy

BACKGROUND

Nasopharyngeal carcinoma (NPC), especially the nonkeratinizing type, is a malignant tumor primarily occurring in southern China and Southeast Asia. Chemotherapy (CT) and combined radiotherapy (RT) is used to treat NPC. However, the mortality rate is high in recurrent and metastatic NPC. We developed a molecular marker, analyzed its correlation with clinical characteristics, and assessed the prognostic value among NPC patients with or without chemoradiotherapy.

METHODS

A total of 157 NPC patients were included in this study, with 120 undergoing treatment and 37 without treatment. EBER1/2 expression was investigated using in situ hybridization (ISH). Expression of PABPC1, Ki-67, and p53 was detected with immunohistochemistry. The correlations of EBER1/2 and the expression of the three proteins having clinical features and prognosis were evaluated.

RESULTS

The expression of PABPC1 was associated with age, recurrence, and treatment but not with gender, TNM classification, or the expression of Ki-67, p53, or EBER. High expression of PABPC1 was associated with poor overall survival (OS) and disease-free survival (DFS) and was an independent predictor depending on multivariate analysis. Comparatively, no significant correlation was observed between the expression of p53, Ki-67, and EBER and survival. In this study, 120 patients received treatments and revealed significantly better OS and DFS than the untreated 37 patients. PABPC1 high expression was an independent predictor of shorter OS in the treated (HR = 4.012 (1.238-13.522), 95% CI, p = 0.021) and the untreated groups (HR = 5.473 (1.051-28.508), 95% CI, p = 0.044). However, it was not an independent predictor of shorter DFS in either the treated or the untreated groups. No significant survival difference was observed between patients with docetaxel-based induction chemotherapy (IC) + concurrent chemoradiotherapy (CCRT) and those with paclitaxel-based IC + CCRT. However, when combined with treatment and PABPC1 expression, patients with paclitaxel-added chemoradiotherapy plus PABPC1 low expression had significantly better OS than those who underwent chemoradiotherapy (p = 0.036).

CONCLUSIONS

High expression of PABPC1 is associated with poorer OS and DFS among NPC patients. Patients with PABPC1 having low expression revealed good survival irrespective of the treatment received, indicating that PABPC1 could be a potential biomarker for triaging NPC patients.

BmNPV circular RNA-encoded peptide VSP39 promotes viral replication

CircRNAs are covalently closed single-stranded circular RNA molecules, which are not easily degraded by endonucleases and play vital roles in many biological processes. Currently, most studies on circRNAs focus on endogenous circRNAs in cells, and there are few studies on virus-encoded circRNAs. In this study, a viral circRNA (circRNA-000010) derived from the region (-/bp: 114514-115,319) of the complementary strand of Bombyx mori Nucleopolyhedrovirus (BmNPV) genome was identified with the circRNA-sequencing. The authenticity of viral circRNA-000010 was further confirmed by reverse transcription PCR, reverse transcription-rolling circle amplification (TCA), in situ hybridization, immunofluorescent staining, and Northern blotting. The results of overexpression and knockdown experiments showed that circRNA-000010 promoted viral replication. Furthermore, a viral small peptide VSP39 with 39 amino acid residues translated by circRNA-000010 but not its linear molecule was confirmed. Finally, VSP39 was found to promote viral replication. Our findings indicated that a viral circRNA encoded by BmNPV promoted viral replication. These findings will provide new clues for further understanding coding information of the BmNPV genome and open a new insight for investigating host-virus interactions.

Regulatory circular RNAs in viral diseases: applications in diagnosis and therapy

Circular RNA (circRNA) forms closed loops via back-splicing in precursor mRNA, resisting exonuclease degradation. In higher eukaryotes, protein-coding genes create circRNAs through exon back-splicing. Unlike mRNAs, circRNAs possess unique production and structural traits, bestowing distinct cellular functions and biomedical potential. In this review, we explore the pivotal roles of viral circRNAs and associated RNA in various biological processes. Analysing the interactions between viral circRNA and host cellular machinery yields fresh insights into antiviral immunity, catalysing the development of potential therapeutics. Furthermore, circRNAs serve as enduring biomarkers in viral diseases due to their stable translation within specific tissues. Additionally, a deeper understanding of translational circRNA could expedite the establishment of circRNA-based expression platforms, meeting the rising demand for broad-spectrum viral vaccines. We also highlight the applications of circular RNA in biomarker studies as well as circRNA-based therapeutics. Prospectively, we expect a technological revolution in combating viral infections using circRNA.

Viral miRNA regulation of host gene expression

Viruses have evolved a multitude of mechanisms to combat barriers to productive infection in the host cell. Virally-encoded miRNAs are one such means to regulate host gene expression in ways that benefit the virus lifecycle. miRNAs are small non-coding RNAs that regulate protein expression but do not trigger the adaptive immune response, making them powerful tools encoded by viruses to regulate cellular processes. Diverse viruses encode for miRNAs but little sequence homology exists between miRNAs of different viral species. Despite this, common cellular pathways are targeted for regulation, including apoptosis, immune evasion, cell growth and differentiation. Herein we will highlight the viruses that encode miRNAs and provide mechanistic insight into how viral miRNAs aid in lytic and latent infection by targeting common cellular processes. We also highlight how viral miRNAs can mimic host cell miRNAs as well as how viral miRNAs have evolved to regulate host miRNA expression. These studies dispel the myth that viral miRNAs are subtle regulators of gene expression, and highlight the critical importance of viral miRNAs to the virus lifecycle.

HSV-1 and influenza infection induce linear and circular splicing of the long NEAT1 isoform

The herpes simplex virus 1 (HSV-1) virion host shut-off (vhs) protein cleaves both cellular and viral mRNAs by a translation-initiation-dependent mechanism, which should spare circular RNAs (circRNAs). Here, we show that vhs-mediated degradation of linear mRNAs leads to an enrichment of circRNAs relative to linear mRNAs during HSV-1 infection. This was also observed in influenza A virus (IAV) infection, likely due to degradation of linear host mRNAs mediated by the IAV PA-X protein and cap-snatching RNA-dependent RNA polymerase. For most circRNAs, enrichment was not due to increased circRNA synthesis but due to a general loss of linear RNAs. In contrast, biogenesis of a circRNA originating from the long isoform (NEAT1_2) of the nuclear paraspeckle assembly transcript 1 (NEAT1) was induced both in HSV-1 infection-in a vhs-independent manner-and in IAV infection. This was associated with induction of novel linear splicing of NEAT1_2 both within and downstream of the circRNA. NEAT1_2 forms a scaffold for paraspeckles, nuclear bodies located in the interchromatin space, must likely remain unspliced for paraspeckle assembly and is up-regulated in HSV-1 and IAV infection. We show that NEAT1_2 splicing and up-regulation can be induced by ectopic co-expression of the HSV-1 immediate-early proteins ICP22 and ICP27, potentially linking increased expression and splicing of NEAT1_2. To identify other conditions with NEAT1_2 splicing, we performed a large-scale screen of published RNA-seq data. This uncovered both induction of NEAT1_2 splicing and poly(A) read-through similar to HSV-1 and IAV infection in cancer cells upon inhibition or knockdown of CDK7 or the MED1 subunit of the Mediator complex phosphorylated by CDK7. In summary, our study reveals induction of novel circular and linear NEAT1_2 splicing isoforms as a common characteristic of HSV-1 and IAV infection and highlights a potential role of CDK7 in HSV-1 or IAV infection.

Biogenesis and Regulatory Roles of Circular RNAs

Covalently closed, single-stranded circular RNAs can be produced from viral RNA genomes as well as from the processing of cellular housekeeping noncoding RNAs and precursor messenger RNAs. Recent transcriptomic studies have surprisingly uncovered that many protein-coding genes can be subjected to backsplicing, leading to widespread expression of a specific type of circular RNAs (circRNAs) in eukaryotic cells. Here, we discuss experimental strategies used to discover and characterize diverse circRNAs at both the genome and individual gene scales. We further highlight the current understanding of how circRNAs are generated and how the mature transcripts function. Some circRNAs act as noncoding RNAs to impact gene regulation by serving as decoys or competitors for microRNAs and proteins. Others form extensive networks of ribonucleoprotein complexes or encode functional peptides that are translated in response to certain cellular stresses. Overall, circRNAs have emerged as an important class of RNAmolecules in gene expression regulation that impact many physiological processes, including early development, immune responses, neurogenesis, and tumorigenesis.

Circ-Udg Derived from Cyprinid Herpesvirus 2 Promotes Viral Replication

Cyprinid herpesvirus 2 (CyHV-2) has caused great losses to the gibel carp (Carassius auratus gibelio) industry. Previous studies showed that certain DNA viruses can encode circular RNAs (circRNAs) to regulate virus infection, which provides new clues for the treatment of viral disease. Whether CyHV-2 can encode circRNAs is still unknown. Here, 10 CyHV-2-derived circRNAs were identified, and the function of circ-udg, a circRNA derived from the CyHV-2 uracil DNA glycosylase (udg) gene, was studied. Although the expression level of circ-udg was lower than that of the parental gene, udg, its expression level was elevated in tandem with the proliferation of CyHV-2 and udg. In vitro experiments confirmed that circ-udg could promote the proliferation of CyHV-2. Moreover, circ-udg could encode a truncated UDG protein consisting of 147-amino-acid residues (termed circ-udg-P147). Both UDG and circ-udg-P147 were found to promote CyHV-2 proliferation, but the promoting effect of circ-udg on CyHV-2 proliferation was attenuated after circ-udg lost the ability to encode circ-udg-P147. Also, circ-udg-P147 could not change the transcription level of the udg gene. Interestingly, the UDG protein level was increased by circ-udg-P147. These results deepen the understanding of the genetic information carried by the genome of CyHV-2 and provide a new target for the treatment of gibel carp bleeding disease caused by CyHV-2. IMPORTANCE The outbreak of C. auratus gibelio gill hemorrhagic disease caused by CyHV-2 brought great losses to the gibel carp industry. Therefore, exploring the interaction between CyHV-2 and host and the molecular mechanism of viral infection is of great significance in preventing and treating the gibel carp gill hemorrhagic disease. Although some progress has been made in the study of CyHV-2, the mechanism of interaction between CyHV-2 and crucian carp is still unclear. In this study, we found that CyHV-2 can encode circRNA to regulate virus replication. Our study provides novel information on CyHV-2 functional genomics, a reference for research into the circRNA of other viruses, and theoretical guidance for preventing and treating gibel carp bleeding disease.

SARS-CoV-2, SARS-CoV, and MERS-CoV encode circular RNAs of spliceosome-independent origin

Circular RNAs (circRNAs) are a newly recognized component of the transcriptome with critical roles in autoimmune diseases and viral pathogenesis. To address the importance of circRNA in RNA viral transcriptome, we systematically identified and characterized circRNAs encoded by the RNA genomes of betacoronaviruses using both bioinformatical and experimental approaches. We predicted 351, 224, and 2764 circRNAs derived from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), SARS-CoV, and Middle East respiratory syndrome coronavirus, respectively. We experimentally identified 75 potential SARS-CoV-2 circRNAs from RNA samples extracted from SARS-CoV-2-infected Vero E6 cells. A systematic comparison of viral and host circRNA features, including abundance, strand preference, length distribution, circular exon numbers, and breakpoint sequences, demonstrated that coronavirus-derived circRNAs had a spliceosome-independent origin. We further showed that back-splice junctions (BSJs) captured by inverse reverse-transcription polymerase chain reaction have different level of resistance to RNase R. Through northern blotting with a BSJ-spanning probe targeting N gene, we identified three RNase R-resistant bands that represent SARS-CoV-2 circRNAs that are detected cytoplasmic by single-molecule and amplified fluorescence in situ hybridization assays. Lastly, analyses of 169 sequenced BSJs showed that both back-splice and forward-splice junctions were flanked by homologous and reverse complementary sequences, including but not limited to the canonical transcriptional regulatory sequences. Our findings highlight circRNAs as an important component of the coronavirus transcriptome, offer important evaluation of bioinformatic tools in the analysis of circRNAs from an RNA genome, and shed light on the mechanism of discontinuous RNA synthesis.

Circular RNAs Represent a Novel Class of Human Cytomegalovirus Transcripts

Human cytomegalovirus (HCMV) infects a large portion of the human population globally. Several HCMV-derived noncoding RNAs are involved in the regulation of viral gene expression and the virus life cycle. Here, we reported that circRNAs are a new class of HCMV transcripts. We bioinformatically predict 704 candidate circRNAs encoded by the TB40/E strain and 230 encoded by the HAN strain. We also systematically compare circRNA features, including the breakpoint sequence consensus, strand preference, length distribution, and exon numbers between host genome-encoded circRNAs and viral circRNAs, and showed that the unique characteristics of viral circRNAs are correlated with their genome types. Furthermore, we experimentally confirmed 324 back-splice junctions (BSJs) from three HCMV strains, Towne, TB40/E, and Toledo, and identified 4 representative HCMV circRNAs by RNase R treatment. Interestingly, we also showed that HCMV contains alternative back-splicing circRNAs. We developed a new amplified FISH method that allowed us to visualize circRNAs and quantify the number of circRNA molecules in the infected cells. The competitive endogenous RNA network analysis suggests that HCMV circRNAs play important roles in viral DNA synthesis via circRNA-miRNA-mRNA networks. Our findings highlight that circRNAs are an important component of the HCMV transcriptome that may contribute to viral replication and pathogenesis. IMPORTANCE HCMV infects 40% to 100% of the human population globally and may be a life-threatening pathogen in immunocompromised individuals. CircRNA is a family of unique RNA that is the most newly found and remains unknown in many aspects. Our current studies computationally identified HCMV-encoded circRNAs and confirmed the existence of the HCMV circRNAs in the infected cells. We systematically compared the features between host and different viral circRNAs and found that the unique characteristics of circRNAs were correlated with their genome types. We also first reported that HCMV contained alternative back-splicing circRNAs. More importantly, we developed a new amplified FISH method which allowed us for the first time not only to visualize circRNAs but also to quantify the number of circRNA molecules in the infected cells. This work describes a novel component of HCMV transcriptome bringing a new understanding of HCMV biology and disease.

The Cellular and Viral circRNAs Induced by Fowl Adenovirus Serotype 4 Infection

Circular RNAs (circRNAs) are a new class of noncoding RNAs that play vital roles in many biological processes. Virus infection induces modifications in cellular circRNA transcriptomes and expresses viral circRNAs. The outbreaks of Hydropericardium-hepatitis syndrome (HHS) caused by fowl adenovirus serotype 4 (FAdV-4) have resulted in huge economic losses to the poultry industry worldwide. To investigate the expression of circRNAs during FAdV-4 infection, we performed transcriptome analysis of FAdV-4-infected leghorn male hepatoma (LMH) cells. In total, 19,154 cellular circRNAs and 135 differentially expressed (DE) cellular circRNAs were identified. The characteristics of the DE cellular circRNAs were analyzed and most of them were related to multiple biological processes according to GO and KEGG enrichment analysis. The accuracy of 10 cellular circRNAs were verified by semiquantitative RT-PCR and sequencing. The change trend was consistent with the RNA sequencing results. Moreover, 2014 viral circRNAs were identified and 10 circRNAs were verified by the same methods. Our analysis showed that seven circRNAs with the same 3′ terminal and variable 5′ terminal regions were located at pTP protein and DNA pol protein of FAdV-4, which may be generated via alternative splicing events. Moreover, the expression level of viral circRNAs was closely related to the replication efficiency of the virus and partial of the viral circRNAs promoted the replication of FAdV-4. Competing endogenous RNA analysis further showed that the effects of cellular and viral circRNAs on host or viral genes may act via miRNAs. Collectively, our findings first indicate that FAdV-4 infection induced the differential expression of cellular circRNAs and FAdV-4 also expressed viral circRNAs, some of which affected FAdV-4 replication. These findings will provide new clues for further understanding FAdV-4 and provide a basis for investigating host-virus interactions.

Viral Circular RNAs and Their Possible Roles in Virus-Host Interaction

Circular RNAs (circRNAs) as novel regulatory molecules have been recognized in diverse species, including viruses. The virus-derived circRNAs play various roles in the host biological process and the life cycle of the viruses. This review summarized the circRNAs from the DNA and RNA viruses and discussed the biogenesis of viral and host circRNAs, the potential roles of viral circRNAs, and their future perspective. This review will elaborate on new insights gained on viruses encoded circRNAs during virus infection.

Identification and characterization of virus-encoded circular RNAs in host cells

Emerging evidence has identified viral circular RNAs (circRNAs) in human cells infected by viruses, interfering with the immune system and inducing diseases including human cancer. However, the biogenesis and regulatory mechanisms of virus-encoded circRNAs in host cells remain unknown. In this study, we used the circRNA detection tool CIRI2 to systematically determine the virus-encoded circRNAs in virus-infected cancer cell lines and cancer patients, by analysing RNA-Seq datasets derived from RNase R-treated samples. Based on the thousands of viral circRNAs we identified, the biological characteristics and potential roles of viral circRNAs in regulating host cell function were determined. In addition, we developed a Viral-circRNA Database (http://www.hywanglab.cn/vcRNAdb/), which is open to all users to search, browse and download information on circRNAs encoded by viruses upon infection.

Micropeptide vsp21 translated by Reovirus circular RNA 000048 attenuates viral replication

To date, some DNA viruses and single-stranded RNA viruses have been found to generate circRNAs. However, the reports on circRNAs produced by double-stranded RNA viruses are very limited. In this study, Bombyx mori cypovirus (BmCPV), a typical double-stranded RNA virus belonging to the Reoviridae, was demonstrated to generate viral circRNAs (vcircRNAs) and a vcircRNA_000048 whose sequence corresponds with the region 164-1245 nt on the BmCPV genomic dsRNA S5 segment (GQ294468.1) was validated by PCR, Sanger sequencing, reverse transcription-rolling circle amplification, and Northern blotting. Furthermore, we verified that vcircRNA_000048 translates a micropeptide vsp21 with 21 amino acid residues in an IRES-dependent manner, and vsp21 attenuates the viral replication. These findings provided a novel clue to understanding the regulation of viral multiplication and interaction of reovirus with the host.

Marek's Disease Virus Virulence Genes Encode Circular RNAs

Circular RNAs (circRNAs) are a recently rediscovered class of functional noncoding RNAs that are involved in gene regulation and cancer development. Next-generation sequencing approaches identified circRNA fragments and sequences underlying circularization events in virus-induced cancers. In the present study, we performed viral circRNA expression analysis and full-length sequencing in infections with Marek’s disease virus (MDV), which serves as a model for herpesvirus-induced tumorigenesis. We established inverse PCRs to identify and characterize circRNA expression from the repeat regions of the MDV genome during viral replication, latency, and reactivation. We identified a large variety of viral circRNAs through precise mapping of full-length circular transcripts and detected matching sequences with several viral genes. Hot spots of circRNA expression included the transcriptional unit of the major viral oncogene encoding the Meq protein and the latency-associated transcripts (LATs). Moreover, we performed genome-wide bioinformatic analyses to extract back-splice junctions from lymphoma-derived samples. Using this strategy, we found that circRNAs were abundantly expressed in vivo from the same key virulence genes. Strikingly, the observed back-splice junctions do not follow a unique canonical pattern, compatible with the U2-dependent splicing machinery. Numerous noncanonical junctions were observed in viral circRNA sequences characterized from in vitro and in vivo infections. Given the importance of the genes involved in the transcription of these circRNAs, our study contributes to our understanding and complexity of this deadly pathogen.

IMPORTANCE Circular RNAs (circRNAs) were rediscovered in recent years both in physiological and pathological contexts, such as in cancer. Viral circRNAs are encoded by at least two human herpesviruses, the Epstein Barr virus and the Kaposi’s Sarcoma-associated herpesvirus, both associated with the development of lymphoma. Marek’s disease virus (MDV) is a well-established animal model to study virus-induced lymphoma but circRNA expression has not been reported for MDV yet. Our study provided the first evidence of viral circRNAs that were expressed at key steps of the MDV lifecycle using genome-wide analyses of circRNAs. These circRNAs were primarily found in transcriptional units that corresponded to the major MDV virulence factors. In addition, we established a bioinformatics pipeline that offers a new tool to identify circular RNAs in other herpesviruses. This study on the circRNAs provided important insights into major MDV virulence genes and herpesviruses-mediated gene dysregulation.

ebv-circRPMS1 promotes the progression of EBV-associated gastric carcinoma via Sam68-dependent activation of METTL3

Epstein-Barr virus (EBV) is a tumor virus that is associated with a variety of neoplasms, including EBV-associated gastric carcinoma (EBVaGC). Recently, EBV was reported to generate various circular RNAs (circRNAs). CircRNAs are important regulators of tumorigenesis by modulating the malignant behaviors of tumor cells. However, to date, the functions of ebv-circRNAs in EBVaGC remain poorly understood. In the present study, we observed high ebv-circRPMS1 expression in EBVaGC and showed that ebv-circRPMS1 promoted the proliferation, migration, and invasion and inhibited the apoptosis of EBVaGC cells. In addition, METTL3 was upregulated in GC cells overexpressing ebv-circRPMS1. Mechanistically, ebv-circRPMS1 bound to Sam68 to facilitate its physical interaction with the METTL3 promotor, resulting in the transactivation of METTL3 and cancer progression. In clinical EBVaGC samples, ebv-circRPMS1 was associated with distant metastasis and a poor prognosis. Based on these findings, ebv-circRPMS1 contributed to EBVaGC progression by recruiting Sam68 to the METTL3 promoter to induce METTL3 expression. ebv-circRPMS1, Sam68, and METTL3 might serve as therapeutic targets for EBVaGC.

Genomic Bootstrap Barcodes and Their Application to Study the Evolution of Sarbecoviruses

Recombination creates mosaic genomes containing regions with mixed ancestry, and the accumulation of such events over time can complicate greatly many aspects of evolutionary inference. Here, we developed a sliding window bootstrap (SWB) method to generate genomic bootstrap (GB) barcodes to highlight the regions supporting phylogenetic relationships. The method was applied to an alignment of 56 sarbecoviruses, including SARS-CoV and SARS-CoV-2, responsible for the SARS epidemic and COVID-19 pandemic, respectively. The SWB analyses were also used to construct a consensus tree showing the most reliable relationships and better interpret hidden phylogenetic signals. Our results revealed that most relationships were supported by just a few genomic regions and confirmed that three divergent lineages could be found in bats from Yunnan: SCoVrC, which groups SARS-CoV related coronaviruses from China; SCoV2rC, which includes SARS-CoV-2 related coronaviruses from Southeast Asia and Yunnan; and YunSar, which contains a few highly divergent viruses recently described in Yunnan. The GB barcodes showed evidence for ancient recombination between SCoV2rC and YunSar genomes, as well as more recent recombination events between SCoVrC and SCoV2rC genomes. The recombination and phylogeographic patterns suggest a strong host-dependent selection of the viral RNA-dependent RNA polymerase. In addition, SARS-CoV-2 appears as a mosaic genome composed of regions sharing recent ancestry with three bat SCoV2rCs from Yunnan (RmYN02, RpYN06, and RaTG13) or related to more ancient ancestors in bats from Yunnan and Southeast Asia. Finally, our results suggest that viral circular RNAs may be key molecules for the mechanism of recombination.

Hypoxia-induced ebv-circLMP2A promotes angiogenesis in EBV-associated gastric carcinoma through the KHSRP/VHL/HIF1α/VEGFA pathway

EBV-encoded circular RNA LMP2A (ebv-circLMP2A) was found to be expressed in EBV-associated gastric carcinoma (EBVaGC) and associated with distant metastasis and poor prognosis. Angiogenesis is a key step in tumor invasion and metastasis and plays a crucial role in tumor progression. However, it is unclear whether and how ebv-circLMP2A is involved in angiogenesis. In this study, we showed that MVD, HIF1α, and VEGFA expression was increased in EBVaGC mouse xenografts with high expression of ebv-circLMP2A. The expression of ebv-circLMP2A was positively correlated with MVD, HIF1α, and VEGFA expression in clinical samples of EBVaGC. Knockdown of ebv-circLMP2A repressed tube formation and migration of HUVECs and decreased VEGFA and HIF1α expression in cancer cells under hypoxia, while ectopic expression of ebv-circLMP2A reversed these effects. Additionally, knockdown of HIF1α blocked the upregulation of ebv-circLMP2A by hypoxia, and ebv-circLMP2A interacted with KHSRP to enhance KHSRP-mediated decay of VHL mRNA, leading to the accumulation of HIF1α under hypoxia. There was a positive feedback loop between HIF1α and ebv-circLMP2A that promotes angiogenesis under hypoxia. ebv-circLMP2A was essential in regulating tumor angiogenesis in EBVaGC and might provide a valuable therapeutic target for EBVaGC.

Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a new member of the Betacoronaviridae family, responsible for the recent pandemic outbreak of COVID-19. To start exploring the molecular events that follow host cell infection, we queried VirusCircBase and identified a circular RNA (circRNA) predicted to be synthesized by SARS-CoV-2, circ_3205, which we used to probe: (i) a training cohort comprised of two pools of cells from three nasopharyngeal swabs of SARS-CoV-2 infected (positive) or uninfected (negative, UCs) individuals; (ii) a validation cohort made up of 12 positive and 3 negative samples. The expression of circRNAs, miRNAs and miRNA targets was assayed through real-time PCR. CircRNA-miRNA interactions were predicted by TarpMiR, Analysis of Common Targets for circular RNAs (ACT), and STarMir tools. Enrichment of the biological processes and the list of predicted miRNA targets were retrieved from DIANA miRPath v3.0. Our results showed that the predicted SARS-CoV-2 circ_3205 was expressed only in positive samples and its amount positively correlated with that of SARS-CoV-2 Spike (S) mRNA and the viral load (r values = 0.80952 and 0.84867, Spearman's correlation test, respectively). Human (hsa) miR-298 was predicted to interact with circ_3205 by all three predictive tools. KCNMB4 and PRKCE were predicted as hsa-miR-298 targets. Interestingly, the function of both is correlated with blood coagulation and immune response. KCNMB4 and PRKCE mRNAs were upregulated in positive samples as compared to UCs (6 and 8.1-fold, p values = 0.049 and 0.02, Student's t test, respectively) and their expression positively correlated with that of circ_3205 (r values = 0.6 and 0.25, Spearman's correlation test, respectively). We propose that our results convincingly suggest that circ_3205 is a circRNA synthesized by SARS-CoV-2 upon host cell infection and that it may behave as a competitive endogenous RNA (ceRNA), sponging hsa-miR-298 and contributing to the upregulation of KCNMB4 and PRKCE mRNAs.

Integrative profiling of Epstein-Barr virus transcriptome using a multiplatform approach

BACKGROUND

Epstein-Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.

METHODS

In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other's data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.

RESULTS

This study detected novel genes embedded into longer host genes containing 5'-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.

CONCLUSIONS

An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

Bombyx mori Akirin hijacks a viral peptide vSP27 encoded by BmCPV circRNA and activates the ROS-NF-κB pathway against viral infection

Bombyx mori cypovirus (BmCPV), a member of the family Reoviridae, is a model of Cypovirus, has a 10 segmented double-stranded RNA genome. However, so far, only one viral small peptide vSP27 with negative regulation on viral infection was identified; the mechanisms underlying host-BmCPV interaction are still unknown. Here, we identified that vSP27 was translated from a BmCPV derived circular RNA (circRNA-vSP27). Subsequently, results showed that vSP27 induced generation of ROS activated the NF-κB signaling pathway, induced the expression of antimicrobial peptides, and suppressed BmCPV infection. On the other hand, we identified a nuclear protein Akirin that could hijack vSP27, positively regulate the NF-κB pathway, and lead to inhibiting the viral infection. Altogether, our data suggested that BmCPV derived circRNA-vSP27 with small peptide translation activity may be employed by the host immunity in defense against the BmCPV infection.

The Cellular and Viral circRNAome Induced by Respiratory Syncytial Virus Infection

Circular RNAs (circRNAs) are a new class of noncoding RNAs that have gained increased attention. DNA virus infections have been reported to induce modifications in cellular circRNA transcriptomes and express viral circRNAs. However, the identification and expression of cellular and viral circRNAs are unknown in the context of respiratory syncytial virus (RSV), a human RNA virus with no effective treatments or vaccines. Here, we report a comprehensive identification of the cellular and viral circRNAs induced by RSV infection in A549 cells with high-throughput sequencing. In total, 53,719 cellular circRNAs and 2,280 differentially expressed cellular circRNAs were identified. Trend analysis further identified three significant expression pattern clusters, which were related to the antiviral immune response according to gene enrichment analysis. Subsequent results showed that not only RSV infection but also poly(I·C) treatment and another RNA virus infection induced the upregulation of the top 10 circRNAs from the focused cluster. The top 10 circRNAs generally inhibit RSV replication in turn. Moreover, 1,254 viral circRNAs were identified by the same circRNA sequencing. The induced expression of viral circRNAs by RSV infection was found not only in A549 cells but also in HEp-2 cells. Additionally, we profiled the general characteristics of both cellular and viral circRNAs such as back-splicing signals, etc. Collectively, RSV infection induced the differential expression of cellular circRNAs, some of which affected RSV infection, and RSV also expressed viral circRNAs. Our study reveals novel layers of host-RSV interactions and identifies cellular or viral circRNAs that may be novel therapeutic targets or biomarkers. IMPORTANCE Noncoding RNAs (ncRNAs) demonstrate substantial roles in cell-virus interactions. Circular RNAs (circRNAs) are a newly identified class of ncRNAs that have gained increased attention recently. DNA virus infections have been reported to induce modifications in cellular circRNA transcriptomes and express viral circRNAs. However, the identification and expression of cellular and viral circRNAs are unknown in the context of respiratory syncytial virus (RSV), a human RNA virus with no effective treatments or vaccines. Here, we report a comprehensive identification of the cellular and viral circRNAs induced by RSV infection by high-throughput sequencing. We revealed that RSV infection induces the differential expression of cellular circRNAs, some of which affected RSV infection, and that RSV also expresses viral circRNAs. Our study reveals novel layers of host-RSV interactions and identifies cellular or viral circRNAs that may be novel therapeutic targets or biomarkers.

CircEAF2 counteracts Epstein-Barr virus-positive diffuse large B-cell lymphoma progression via miR-BART19-3p/APC/β-catenin axis

Background

Epstein-Barr virus (EBV) represents an important pathogenic factor of lymphoma and is significantly associated with poor clinical outcome of diffuse large B-cell lymphoma (DLBCL). Circular RNAs (circRNAs) play an essential role in lymphoma progression. However, the underlying mechanism of circRNA on DLBCL progression related to EBV remains largely unknown.

Methods

CircRNA was screened by high-throughput sequencing in tumor samples of 12 patients with DLBCL according to EBV infection status. Expression of circEAF2, as well as the relationship with clinical characteristics and prognosis, were further analyzed in tumor samples of 100 DLBCL patients using quantitative real-time PCR. Gain- and loss-of-function experiments were conducted to investigate the biological functions of circEAF2 both in vitro and in vivo. The underlying mechanism of circRNA on DLBCL progression were further determined by RNA sequencing, RNA pull down assay, dual-luciferase reporter assay, rescue experiments and western blotting.

Results

We identified a novel circRNA circEAF2, which was downregulated in EBV + DLBCL and negatively correlated with EBV infection and DLBCL progression. In EBV-positive B lymphoma cells, circEAF2 overexpression induced lymphoma cell apoptosis and sensitized lymphoma cells to epirubicin. As mechanism of action, circEAF2 specifically targeted EBV-encoded miR-BART19-3p, upregulated APC, and suppressed downstream β-catenin expression, resulting in inactivation of Wnt signaling pathway and inhibition of EBV + DLBCL cell proliferation. In EBV-positive B-lymphoma murine models, xenografted tumors with circEAF2 overexpression presented decreased Ki-67 positivity, increased cell apoptosis and retarded tumor growth.

Conclusions

CircEAF2 counteracted EBV + DLBCL progression via miR-BART19-3p/APC/β-catenin axis, referring circEAF2 as a potential prognostic biomarker. Therapeutic targeting EBV-encoded miRNA may be a promising strategy in treating EBV-associated lymphoid malignancies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01458-9.

Quantifying full-length circular RNAs in cancer

Circular RNAs (circRNAs) are abundantly expressed in cancer. Their resistance to exonucleases enables them to have potentially stable interactions with different types of biomolecules. Alternative splicing can create different circRNA isoforms that have different sequences and unequal interaction potentials. The study of circRNA function thus requires knowledge of complete circRNA sequences. Here we describe psirc, a method that can identify full-length circRNA isoforms and quantify their expression levels from RNA sequencing data. We confirm the effectiveness and computational efficiency of psirc using both simulated and actual experimental data. Applying psirc on transcriptome profiles from nasopharyngeal carcinoma and normal nasopharynx samples, we discover and validate circRNA isoforms differentially expressed between the two groups. Compared with the assumed circular isoforms derived from linear transcript annotations, some of the alternatively spliced circular isoforms have 100 times higher expression and contain substantially fewer microRNA response elements, showing the importance of quantifying full-length circRNA isoforms.

A Previously Undiscovered Circular RNA, circTNFAIP3, and Its Role in Coronavirus Replication

Circular RNAs (circRNAs) are a newly discovered class of noncoding RNAs (ncRNAs) present in various tissues and cells. However, the functions of most circRNAs have not been verified experimentally. Here, using deltacoronavirus as a model, differentially expressed circRNAs in cells with or without deltacoronavirus infection were analyzed by RNA sequencing to characterize the cellular responses to RNA virus infection. More than 57,000 circRNA candidates were detected, and seven significantly dysregulated circRNAs were quantitated by real-time reverse transcription-PCR. We discovered a previously unidentified circRNA derived from the TNFAIP3 gene, named circTNFAIP3, which is distributed and expressed widely in various tissues. RNA viruses, including deltacoronaviruses, rather than DNA viruses tend to activate the expression of endogenous circTNFAIP3. Overexpression of circTNFAIP3 promoted deltacoronavirus replication by reducing the apoptosis, while silencing of circTNFAIP3 inhibited deltacoronavirus replication by enhancing the apoptosis. In summary, our work provides useful circRNA-related information to facilitate investigation of the underlying mechanism of deltacoronavirus infection and identifies a novel circTNFAIP3 that promotes deltacoronavirus replication via regulating apoptosis.

Comparative Transcriptome Provides a Systematic Perspective on Epstein-Barr Virus-Associated Gastric Carcinoma Cell Lines

Purpose

Epstein–Barr virus (EBV) is widely recognised to cause various tumours, and EBV-associated gastric carcinoma (EBVaGC) is a special type of GC. It has obviously different clinical features and pathological manifestations from EBV-negative gastric carcinoma, but its progression remains elusive. The underlying cancer progression of viral infection detected by genome-wide transcriptome analysis has been demonstrated in numerous diseases.

Methods

We performed comparative RNA sequencing to identify gene expression signatures between GC and EBVaGC cell lines. The differentially expressed (DE) genes were analysed using gene ontology and pathway enrichment.

Results

A total of 4438 DE mRNAs, 3650 DE long non-coding RNAs (lncRNAs), and 248 DE circular RNAs (circRNAs) were detected in GC cells after EBV infection, most of which were highly related to oncogenesis. Likewise, EBV-coding RNA and non-coding RNA were also well-supplemented in EBVaGC. According to bioinformatics, DE mRNAs may contribute to the completion of EBV-infected host cells and modulate mitosis. Binding to actin and participating in adherens junctions to promote contact between the virus and cells are a potential function of DE lncRNAs. The roles of DE circRNAs were enriched in DNA repair and protein modification, and a typical example of this is acting as an miRNA sponge. The establishment of a circRNA-miRNA-mRNA network helps to determine the key elements in the progression of EBVaGC.

Conclusion

This study is the first to systematically reveal the transcriptome landscape of EBVaGC, which will provide an essential resource for genomic, genetic, and molecular mechanisms in the future.

Epstein-Barr Virus-Encoded Circular RNA CircBART2.2 Promotes Immune Escape of Nasopharyngeal Carcinoma by Regulating PD-L1

Epstein-Barr virus (EBV) infection is an established cause of nasopharyngeal carcinoma (NPC) and is involved in a variety of malignant phenotypes, including tumor immune escape. EBV can encode a variety of circular RNAs (circRNA), however, little is known regarding the biological functions of these circRNAs in NPC. In this study, EBV-encoded circBART2.2 was found to be highly expressed in NPC where it upregulated PD-L1 expression and inhibited T-cell function in vitro and in vivo. circBART2.2 promoted transcription of PD-L1 by binding the helicase domain of RIG-I and activating transcription factors IRF3 and NF-κB, resulting in tumor immune escape. These results elucidate the biological function of circBART2.2, explain a novel mechanism of immune escape caused by EBV infection, and provide a new immunotherapy target for treating NPC. SIGNIFICANCE: This work demonstrates that circBART2.2 binding to RIG-I is essential for the regulation of PD-L1 and subsequent immune escape in nasopharyngeal carcinoma.

FLAME: long-read bioinformatics tool for comprehensive spliceome characterization

Comprehensive characterization of differentially spliced RNA transcripts with nanopore sequencing is limited by bioinformatics tools that are reliant on existing annotations. We have developed FLAME, a bioinformatics pipeline for alternative splicing analysis of gene-specific or transcriptome-wide long-read sequencing data. FLAME is a Python-based tool aimed at providing comprehensible quantification of full-length splice variants, reliable de novo recognition of splice sites and exons, and representation of consecutive exon connectivity in the form of a weighted adjacency matrix. Notably, this workflow circumvents issues related to inadequate reference annotations and allows for incorporation of short-read sequencing data to improve the confidence of nanopore sequencing reads. In this study, the Epstein-Barr virus long noncoding RNA RPMS1 was used to demonstrate the utility of the pipeline. RPMS1 is ubiquitously expressed in Epstein-Barr virus associated cancer and known to undergo ample differential splicing. To fully resolve the RPMS1 spliceome, we combined gene-specific nanopore sequencing reads from a primary gastric adenocarcinoma and a nasopharyngeal carcinoma cell line with matched publicly available short-read sequencing data sets. All previously reported splice variants, including putative ORFs, were detected using FLAME. In addition, 32 novel exons, including two intron retentions and a cassette exon, were discovered within the RPMS1 gene.

Characterizing viral circRNAs and their application in identifying circRNAs in viruses

Circular RNAs (circRNAs) are non-coding RNAs with a special circular structure produced formed by the reverse splicing mechanism, which play an important role in a variety of biological activities. Viruses can encode circRNA, and viral circRNAs have been found in multiple single-stranded and double-stranded viruses. However, the characteristics and functions of viral circRNAs remain unknown. Sequence alignment showed that viral circRNAs are less conserved than circRNAs in animal, indicating that the viral circRNAs may evolve rapidly. Through the analysis of the sequence characteristics of viral circRNAs and circRNAs in animal, it was found that viral circRNAs and animals circRNAs are similar in nucleic acid composition, but have obvious differences in secondary structure and autocorrelation characteristics. Based on these characteristics of viral circRNAs, machine learning algorithms were employed to construct a prediction model to identify viral circRNA. Additionally, analysis of the interaction between viral circRNA and miRNAs showed that viral circRNA is expected to interact with 518 human miRNAs, and preliminary analysis of the role of viral circRNA. And it has been also found that viral circRNAs may be involved in many KEGG pathways related to nervous system and cancer. We curated an online server, and the data and code are available: http://server.malab.cn/viral-CircRNA/.

Hepatocellular carcinoma progression mediated by hepatitis B virus-encoded circRNA HBV_circ_1 through interaction with CDK1

Hepatitis B virus (HBV) produces circular RNA (circRNA), whose functions have not yet been clearly elucidated. In this study, a novel circRNA HBV_circ_1 produced by HBV was identified in HBV-positive HepG2.2.15 cells and HBV-related hepatocellular carcinoma (HCC) tissue (HCCT). Microarray analysis of 68 HCCT samples showed that HBV_circ_1 abundance was significantly higher than that in paracancerous tissues. In addition, survival rate of HBV_circ_1-positive patients was significantly lower compared with HBV_circ_1-negative patients. Transient expression indicated that HBV_circ_1 enhanced cell proliferation, migration, and invasion and inhibited apoptosis in vitro. Furthermore, ectopical HBV_circ_1 expression increased tumor size in vivo. HBV_circ_1 was confirmed to interact with cyclin-dependent kinase 1 (CDK1) to regulate cell proliferation. These results suggest that HCC progression may be promoted by interaction of HBV_circ_1 with CDK1. Our data not only showed a novel clue to understand carcinogenesis and progress of HBV-related HCC but also provided a new target for the development of therapeutic drugs.

Role of Virus-Induced Host Cell Epigenetic Changes in Cancer

The tumor viruses human T-lymphotropic virus 1 (HTLV-1), hepatitis C virus (HCV), Merkel cell polyomavirus (MCPyV), high-risk human papillomaviruses (HR-HPVs), Epstein-Barr virus (EBV), Kaposi’s sarcoma-associated herpes virus (KSHV) and hepatitis B virus (HBV) account for approximately 15% of all human cancers. Although the oncoproteins of these tumor viruses display no sequence similarity to one another, they use the same mechanisms to convey cancer hallmarks on the infected cell. Perturbed gene expression is one of the underlying mechanisms to induce cancer hallmarks. Epigenetic processes, including DNA methylation, histone modification and chromatin remodeling, microRNA, long noncoding RNA, and circular RNA affect gene expression without introducing changes in the DNA sequence. Increasing evidence demonstrates that oncoviruses cause epigenetic modifications, which play a pivotal role in carcinogenesis. In this review, recent advances in the role of host cell epigenetic changes in virus-induced cancers are summarized.

MicroRNA and Other Non-Coding RNAs in Epstein-Barr Virus-Associated Cancers

EBV is a direct causative agent in around 1.5% of all cancers. The oncogenic properties of EBV are related to its ability to activate processes needed for cellular proliferation, survival, migration, and immune evasion. The EBV latency program is required for the immortalization of infected B cells and involves the expression of non-coding RNAs (ncRNAs), including viral microRNAs. These ncRNAs have different functions that contribute to virus persistence in the asymptomatic host and to the development of EBV-associated cancers. In this review, we discuss the function and potential clinical utility of EBV microRNAs and other ncRNAs in EBV-associated malignancies. This review is not intended to be comprehensive, but rather to provide examples of the importance of ncRNAs.

Herpesvirus Epigenetic Reprogramming and Oncogenesis

Among all of the known biological carcinogens, Epstein-Barr virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV) are two of the classical oncogenic herpesviruses known to induce the oncogenic phenotype. Many studies have revealed important functions related to epigenetic alterations of the EBV and KSHV genomes that mediate oncogenesis, but the detailed mechanisms are not fully understood. It is also challenging to fully describe the critical cellular events that drive oncogenesis as well as a comprehensive map of the molecular contributors. This review introduces the roles of epigenetic modifications of these viral genomes, including DNA methylation, histone modification, chromatin remodeling, and noncoding RNA expression, and elucidates potential strategies utilized for inducing oncogenesis by these human gammaherpesviruses.

Identification and characterization of a novel Epstein-Barr Virus-encoded circular RNA from LMP-2 Gene

Epstein-Barr virus (EBV) has been recently found to generate novel circular RNAs (circRNAs) through backsplicing. However, comprehensive catalogs of EBV circRNAs in other cell lines and their functional characterization are still lacking. In this study, we have identified a list of putative EBV circRNAs in GM12878, an EBV-transformed lymphoblastoid cell line, with a significant majority encoded from the EBV latent genes. A novel EBV circRNA derived from the exon 5 of LMP-2 gene which exhibited highest prevalence, was further validated using RNase R assay and Sanger sequencing. This circRNA, which we term circLMP-2_e5, can be universally detected in a panel of EBV-positive cell lines modelling different latency programs. It ranges from lower expression in nasopharyngeal carcinoma (NPC) cells to higher expression in B cells, and is localized to both the cytoplasm and the nucleus. We provide evidence that circLMP-2_e5 is expressed concomitantly with its cognate linear LMP-2 RNA upon EBV lytic reactivation, and may be produced as a result of exon skipping, with its circularization possibly occurring without the involvement of cis elements in the short flanking introns. Furthermore, we show that circLMP-2_e5 is not involved in regulating cell proliferation, host innate immune response, its linear parental transcripts, or EBV lytic reactivation. Taken together, our study expands the current repertoire of putative EBV circRNAs, broadens our understanding of the biology of EBV circRNAs, and lays the foundation for further investigation of their function in the EBV life cycle and disease development.

Characterizing Expression and Regulation of Gamma-Herpesviral Circular RNAs

Multiple herpesviruses have been recently found to encode viral circular RNAs. Like cellular circular RNAs, these RNAs lack poly-A tails and their 5' and 3' ends have been joined, which confers protection from RNA exonucleases. We examined the expression patterns of circular RNAs from Kaposi's sarcoma herpesvirus (KSHV) in various environments. We performed deep sequencing of circRNA-enriched total RNA from a KSHV-positive patient lymph node for comparison with previous circRNA-Seq results. We found that circvIRF4 is highly expressed in the KSHV-positive patient sample relative to both B cell lines and de novo infected primary vascular and lymphatic endothelial cells (LECs). Overall, this patient sample showed a viral circRNA expression pattern more similar to the pattern from B cell lines, but we also discovered new back-spliced junctions and additional viral circular RNAs in this patient sample. We validated some of these back-spliced junctions as circular RNAs with standard assays. Differential expression patterns of circular RNAs in different cell types led us to investigate what cellular factors might be influencing the ratio of viral linear mRNAs to circular RNAs. We found that repression of certain RNA-binding proteins shifted the balance between viral linear mRNAs and circular RNAs. Taken together, examining viral circular RNA expression patterns may become useful tools for discovering their functions, the regulators of their expression, and determining the stage and cell types of infection in humans.

Long non-coding RNAs in Epstein-Barr virus-related cancer

Epstein Barr-virus (EBV) is related to several cancers. Long non-coding RNAs (lncRNAs) act by regulating target genes and are involved in tumourigenesis. However, the role of lncRNAs in EBV-associated cancers is rarely reported. Understanding the role and mechanism of lncRNAs in EBV-associated cancers may contribute to diagnosis, prognosis and clinical therapy in the future. EBV encodes not only miRNAs, but also BART lncRNAs during latency and the BHLF1 lncRNA during both the latent and lytic phases. These lncRNAs can be targeted regulate inflammation, invasion, and migration and thus tumourigenesis. The products of EBV also directly and indirectly regulate host lncRNAs, including LINC00312, NORAD CYTOR, SHNG8, SHNG5, MINCR, lncRNA-BC200, LINC00672, MALATI1, LINC00982, LINC02067, IGFBP7-AS1, LOC100505716, LOC100128494, NAG7 and RP4-794H19.1, to facilitate tumourigenesis using different mechanisms. Additionally, lncRNAs have been previously validated to interact with microRNAs (miRNAs), and lncRNAs and miRNAs mutually suppress each other. The EBV-miR-BART6-3p/LOC553103/STMN1 axis inhibits EBV-associated tumour cell proliferation. Additionally, H. pylori-EBV co-infection promotes inflammatory lesions and results in EMT. HPV-EBV co-infection inhibits the transition from latency to lytic replication. KSHV-EBV co-infection aggravates tumourigenesis in huNSG mice. COVID-19-EBV co-infection may activate the immune system to destroy a tumour, although this situation is rare and the mechanism requires further confirmation. Hopefully, this information will shed some light on tumour therapy strategies tumourigenesis. Additionally, this strategy benefits for infected patients by preventing latency to lytic replication. Understanding the role and expression of lnRNAs in these two phases of EBV is critical to control the transition from latency to the lytic replication phase. This review presents differential expressed lncRNAs in EBV-associated cancers and provides resources to aid in developing superior strategies for clinical therapy.