key: cord-0017200-vde4rub0
authors: Avilala, Janardhan; Becnel, David; Abdelghani, Ramsy; Nanbo, Asuka; Kahn, Jacob; Li, Li; Lin, Zhen
title: Role of Virally Encoded Circular RNAs in the Pathogenicity of Human Oncogenic Viruses
date: 2021-04-20
journal: Front Microbiol
DOI: 10.3389/fmicb.2021.657036
sha: e4ecbb0afa593ae1518c691bbcc15adf61706dae
doc_id: 17200
cord_uid: vde4rub0

Human oncogenic viruses are a group of important pathogens that etiologically contribute to at least 12% of total cancer cases in the world. As an emerging class of non-linear regulatory RNA molecules, circular RNAs (circRNAs) have gained increasing attention as a crucial player in the regulation of signaling pathways involved in viral infection and oncogenesis. With the assistance of current circRNA enrichment and detection technologies, numerous novel virally-encoded circRNAs (vcircRNAs) have been identified in the human oncogenic viruses, initiating an exciting new era of vcircRNA research. In this review, we discuss the current understanding of the roles of vcircRNAs in the respective viral infection cycles and in virus-associated pathogenesis.

Circular RNAs (circRNAs) are a class of non-linear macromolecules that are single-stranded covalently closed cyclized RNAs with a great diversity in length, varied from less than 100 to thousands of nucleotides (Kolakofsky, 1976; Sanger et al., 1976; Hsu and Coca-Prados, 1979; Arnberg et al., 1980; Salzman et al., 2012 Salzman et al., , 2013 Jeck et al., 2013; Memczak et al., 2013; Jeck and Sharpless, 2014; Rybak-Wolf et al., 2015; Kristensen et al., 2018) . The first evidence of circRNA was reported in 1976 when Sanger et al. conducted the electron microscopy analysis of the denatured RNA molecules from the viroid of the pale fruit disease of cucumber (CPFV). The authors found that these viroid RNA existed as single-stranded circles (Sanger et al., 1976) . Three years later, endogenous circRNAs were further detected in the cytoplasmic RNA extracted from human HeLa cells, monkey CV-1 cells, and Chinese hamster ovary cells as well as total RNA isolated from Physarum polycephalum (a plasmodial slime mold) (Hsu and Coca-Prados, 1979) . However, at that time, the synthesis of circRNAs was considered to be a rare event and those circRNAs were treated as "splicing noise" and byproducts likely resulting from splicing errors (Cocquerelle et al., 1993) . Hence, the potential impact of this discovery did not receive broad recognition and only a handful of follow-up studies were conducted in the next three decades.

The dawn of the circRNA era occurred approximately 30 years later when high-throughput next-generation sequencing technology and associated informatics approaches have been utilized in the circRNA research field. Since the early 2010s, thousands of endogenous circRNAs have been identified in all domains of life including human, animals, plants, yeasts, bacteria, and viruses (Danan et al., 2012; Salzman et al., 2012 Salzman et al., , 2013 Jeck et al., 2013; Memczak et al., 2013; Jeck and Sharpless, 2014; Rybak-Wolf et al., 2015; Kristensen et al., 2018; Tagawa et al., 2018; Toptan et al., 2018; Ungerleider et al., 2018 Ungerleider et al., , 2019a Zhao et al., 2019) . To date, the list of identified circRNAs is still rapidly growing and circRNAs are now appreciated as an important evolutionarily conserved component of living organisms.

Circular RNAs are usually expressed in a tissue-and/or developmental stage-specific manner (Memczak et al., 2013) . The mature circRNAs can be detected in cytoplasm, nucleus, as well as extracellular vesicles (e.g., exosomes) . In general, circRNAs show ∼100-1,000 times lower abundance than their associated isogenic linear RNAs and are generated by more than 20% of expressed genes . However, in many cases, circRNAs can be the main gene products and exhibit more than 10-fold higher expression levels compared to related linear RNAs (Jeck et al., 2013) . Lacking free termini, circRNAs are resistant to the hydrolysis by numerous cellular exonucleases such as RNase R, and thus have significantly longer half-lives compared to their linear RNA counterparts (Jeck et al., 2013) .

Based on the parental gene component that's carried by the circRNAs, circRNAs can be classified to four major types (Figure 1) (Jeck et al., 2013; Zhang et al., 2013; Li Z. et al., 2015; Noto et al., 2017) : (1) the most common type is the exonic circRNAs (ecircRNAs) which contain only one or multiple exons; (2) the exon-intron circRNAs (EIciRNAs) carry both exons and introns; (3) the circular intronic RNAs (ciRNAs) carry introns only; and (4) the tRNA intronic circRNAs (tricRNAs) are formed by circularization of the excised tRNA introns. The detailed mechanisms of generating different types of circRNAs are still less understood. The current theory holds that the formation of ecircRNAs and EIciRNAs requires a special type of splicing known as backsplicing, which allows a downstream 5 splice donor to react with an upstream 3 splice acceptor. Thus, a 3 -5 phosphodiester bond can be established to circularize the RNA molecule (Jeck et al., 2013; Li Z. et al., 2015) . Unlike exon-derived circRNAs, ciRNAs are generated by lariat introns removed from pre-mRNAs by canonical splicing. A consensus 7-nucleotide GU-rich sequence near the 5 splice donor and an 11-nucleotide C-rich sequence near the branchpoint site can help prevent the lariat from debranching and exonuclease attacking. Stable ciRNAs can then be formed after trimming the 3 tail downstream from the branchpoint (Zhang et al., 2013) . Lastly, tricRNAs are formed by directly linking the free ends of introns excised from the pre-tRNAs (Noto et al., 2017) .

In general, the exon-derived circRNAs (ecircRNAs and EIciRNAs) are produced from pre-mRNA transcribed by RNA polymerase II. The circularization signals are located within the introns adjacent to the circularizable exons (Jeck et al., 2013; Zhang et al., 2014) . A growing number of RNAbinding proteins (RBPs) as well as spliceosome components are reported to be involved in the circularization process, including adenosine deaminases acting on RNA (ADAR), the splicing factor muscleblind (MBL), quaking (QKI), DExH-Box helicase 9 (DHX9), epithelial splicing regulatory protein 1 (ESRP1), FUS, serine/arginine-rich proteins, nuclear factors NF90/NF110, small ribonucleoprotein particle U1 subunit 70K (snRNP-U1-70K), small ribonucleoprotein particle U1 subunit C (snRNP-U1-C), pre-mRNA processing 8 (Prp8), cell division cycle 40 (CDC40) (Chen and Yang, 2015; Conn et al., 2015; Ivanov et al., 2015; Kramer et al., 2015; Rybak-Wolf et al., 2015; Aktas et al., 2017; Errichelli et al., 2017; Yu et al., 2017) .

In recent years, the way in which circRNAs execute their function has been the subject of active study. Although the biological significance of the majority of the identified circRNAs are still unclear, it's believed that circRNAs predominantly function as regulatory non-coding RNAs (ncRNAs) (Guo et al., 2014) . circRNAs can function in trans or in cis. For instance, there is increasing evidence showing that some circRNAs (such as CDR1as/ciRS-7) can serve as microRNA (miRNA) sponges to sequester endogenous miRNAs and modulate miRNA function in vivo (Hansen et al., 2013; Piwecka et al., 2017; Kleaveland et al., 2018) . Meanwhile, Ashwal-Fluss et al. (2014) have shown that some circRNAs can be co-transcribed with their associated linear mRNAs in a competitive manner. Thus, these circRNAs can function as an RNA trap for mRNA synthesis by competing with linear splicing in cis. Since circRNAs can bind to numerous RBPs, it's postulated that these stable RNA molecules can serve as decoys, scaffolds, or transporters for many RBPs and regulate their functions (Du et al., 2016; Holdt et al., 2016; Abdelmohsen et al., 2017; Pamudurti et al., 2018) . The discovery of circRNAs in the extracellular compartments such as exosomes indicates that circRNAs might also function as secretory intercellular signaling molecules to affect distant cells . Interestingly, several recent studies have shown that certain circRNAs (e.g., circZNF609) can be translated in a non-canonical cap-independent manner Pamudurti et al., 2017; Yang et al., 2017) . Although the biological impact of circRNA translation and the resulting products are still largely unclear, it's postulated that the truncated protein products/peptides may compete with the full-length protein products encoded by their associated linear mRNAs and affect relevant signaling pathways . Together, this evidence suggests that circRNAs may function as important gene regulators and modulate signaling transduction at different levels.

Cancer is developed through a multi-step and multi-factorial process. As an emerging gene regulator, circRNAs are expected to be involved in the regulation of various carcinogenic pathways. To date, there is emerging evidence showing that disruption of circRNA expression correlates with a broad spectrum of human cancers, including hepatocellular carcinoma (HCC; Qin et al., 2016; Shang et al., 2016; Yu et al., 2016; Xu et al., 2017; Yao et al., 2017) , colorectal cancer Wang et al., 2015; Guo et al., 2016; Xie et al., 2016; Hsiao et al., 2017; Weng et al., 2017; Zhu et al., 2017) , gastric cancer (Li P. et al., 2015 (Li P. et al., , 2017 Li W.H. et al., 2017; Sui et al., 2017; Zhao et al., 2018) , gliomas (Barbagallo et al., 2016; Zheng et al., 2017) , lung cancer Zhang et al., 2020) , breast cancer (Liu Z. et al., 2019) , pancreatic ductal adenocarcinoma , laryngeal FIGURE 1 | Biogenesis and potential functions of circRNAs. Based on the parental gene component that's carried by the circRNAs, circRNAs can be classified to four major types. (1) The most common type is the exonic circRNAs (ecircRNAs) which contain only one or multiple exons; (2) The exon-intron circRNAs (EIciRNAs) carry both exons and introns; (3) The circular intronic RNAs (ciRNAs) carry introns only; and (4) The tRNA intronic circRNAs (tricRNAs) are formed by circularization of the excised tRNA introns. The formation of ecircRNAs and EIciRNAs requires a special type of splicing known as backsplicing, which allows a downstream 5 splice donor (SD) to react with an upstream 3 splice acceptor (SA). Thus, a 3 -5 phosphodiester bond can be established to circularize the RNA molecule. ciRNAs are generated by lariat introns removed from pre-mRNAs by canonical splicing. A consensus 7-nucleotide GU-rich sequence near the 5 splice donor (shown as red dot in the figure) and an 11-nucleotide C-rich sequence near the branchpoint site (shown as green dot in the figure) can help prevent the lariat from debranching and exonuclease attacking. Stable ciRNAs can then be formed after trimming the 3 tail downstream from the branchpoint. Lastly, tricRNAs are formed by directly linking the free ends of introns excised from the pre-tRNAs. circRNAs can function as molecular sponges to sequester miRNA and protein. In general, circRNAs mainly function as non-coding regulatory transcripts. However, some circRNAs can be translated through either IRES-mediated or m 6 A-mediated non-canonical cap-independent translation initiation.

cancer (Xuan et al., 2016) , esophageal squamous cell carcinoma (Xia et al., 2016) , basal cell carcinoma (Sand et al., 2016b) , cutaneous squamous cell carcinoma (Sand et al., 2016a) , ovarian carcinoma (Ahmed et al., 2016) , cervical cancer (Wu et al., 2020) , and nasopharyngeal carcinoma (NPC; Yang et al., 2020) . During the oncogenesis process, circRNAs can serve as either oncogenes or tumor suppressors by inhibiting critical protein-coding and/or non-coding genes involved in oncogenesis. Currently, the study of detailed mechanisms through which circRNAs facilitate cancer development is still in its infancy. Nevertheless, it is apparent that the discovery of circRNAs has created a new dimension for us to elucidate the mechanism of carcinogenesis.

The family of human oncogenic viruses currently comprises seven members: Epstein-Barr virus (EBV), Kaposi's sarcomaassociated herpesvirus (KSHV), human papillomavirus (HPV), Merkel cell polyomavirus (MCV), hepatitis B virus (HBV), hepatitis C virus (HCV), and human T-lymphotropic virus (also known as human T-cell leukemia virus type 1). Together, they account for 12-16% of total cancer cases in the world (White et al., 2014; Zapatka et al., 2020) . Notably, only a small portion of infected persons develop cancers, which usually happens many years after initial infection (Mesri et al., 2014) . It indicates that these oncogenic viruses are not fully oncogenic -they are necessary but not sufficient to cause cancers (Mesri et al., 2014) . Indeed, for all these oncogenic viruses, cancer formation is believed to be a biological accident, but not an adaptation during their co-evolution with the hosts (Moore and Chang, 2017) . After all, cancers are considered as evolutionarily deadend events which put both hosts and viruses in danger (Moore and Chang, 2017) . Viral transmission is the central point for the viral persistence and evolution. However, oncogenic viruses carried by these cancers are usually not permissive to produce sufficient infectious virions in order to support human-tohuman transmission (Moore and Chang, 2017) . The transmission of the oncogenic virus generally occurs among asymptomatic individuals. The dramatic increase of viral cancers in the immunocompromised individuals indicates that host immune surveillance plays a critical role in preventing viral cancer formation in the infected individuals (Moore and Chang, 2017) .

So, why do these oncogenic viruses cause cancers? There is a continual battle between hosts and oncogenic viruses. Hosts have evolved various innate and adaptive immune strategies to eliminate viral infections. Meanwhile, oncogenic viruses have co-evolved strong immunoevasion and replication programs to counteract hosts' defense systems (Mesri et al., 2014) . During the process of immunoevasion and replication, some of the viral strategies such as proliferation and antiapoptosis accidentally break host homeostasis and trigger cancer development. However, since the cancer development is a multifactorial and multistep process, additional co-factors such as chronic inflammation, environmental mutagens and immunosuppression are usually required for the formation of viral cancers (Mesri et al., 2014) .

To date, increasing evidence has indicated that circRNAs play important roles in the regulation of RNA biogenesis and cell signaling. We are not surprised to see oncogenic viruses once again exploit these pathways to suit their needs of survival. Next, we have discussed a number of recently identified vcircRNAs and the possible relevance of these vcircRNA in virus infection and the development of virus-associated human malignancies.

Epstein-Barr virus is also known as EBV or human herpesvirus 4. It is the first human tumor virus identified in 1964 and accounts for ∼200,000 new cancer cases annually (Rickinson and Kieff, 2007; Young et al., 2016; Shannon-Lowe and Rickinson, 2019) . EBV is etiologically related to a number of human malignant diseases including Burkitt's lymphoma (BL), Hodgkin's disease (HD), extranodal nasal-type natural killer/T cell lymphoma (NKTL), posttransplant lymphoproliferative disease (PTLD) and lymphoma in AIDS and immunocompromised individuals, NPC, gastric carcinoma (GC), as well as lung cancers (LC) (de Sanjose et al., 2007; Rickinson and Kieff, 2007; Strong et al., 2013a,b; Young et al., 2016; Kheir et al., 2019; Shannon-Lowe and Rickinson, 2019) .

Epstein-Barr virus belongs to the human γ-herpesvirus subfamily, in the same family with KSHV. There are two distinct stages in EBV's life cycle: latency and lytic cycle (also known as lytic reactivation). After primary infection, EBV co-exists with its host permanently and exhibits a life-long latency. During the latency, the viral genome remains as an episomal chromosome and only a limited repertoire of viral genes are expressed. In response to certain physiological stimuli, the latent genome can be reactivated to express more than 100 viral genes and produce new infectious virions.

Infectious virions can be transmitted through various types of body fluids, but they are mainly spread through saliva exchange. There's no anti-EBV vaccine currently available and it's estimated that ∼90% of the world population permanently carries EBV (Rickinson and Kieff, 2007) . To date, the detailed mechanism of EBV primary infection is still unclear. The current theory holds that EBV initially infects replication-permissive oropharyngeal epithelial cells in the host. The propagated infectious virions are then transmitted from the initial epithelial cells to the infiltrating B-cells in which the full set of viral latency transcripts are expressed (known as type III latency) ( Table 1) . The type III latency can subsequently activate B cells and cause a transient expansion of the EBV-infected B-cell pool. Majority of such B-cells will be eliminated by host immune surveillance. However, a sub-set of infected circulating B-cells can persist for the lifetime of the host. This is attributable to expression of no (type 0 latency) or only one viral protein (type I latency) ( Table 1) . It is believed that sporadic reactivation occurs in epithelial tissues as well as B cells differentiating into plasma cells (Laichalk and Thorley-Lawson, 2005; Sun and Thorley-Lawson, 2007) , and the progeny infectious virions are secreted to the saliva and then transmitted to new hosts.

Previously, EBV has been shown to utilize non-coding linear RNAs such as miRNAs, long non-coding RNAs (lncRNAs), and small non-coding EBV-encoded RNAs (EBERs) to facilitate its life cycle and oncogenesis (Pfeffer et al., 2004; Cai et al., 2006; Hutzinger et al., 2009; Concha et al., 2012; Moss and Steitz, 2013; O'Grady et al., 2014; Cao et al., 2015a,b; O'Grady et al., 2016) . Recently, increasing evidence has shown that EBV also encodes its own set of circRNAs to promote its survival. Since both the virus survival and tumor survival often share the same needs, the expression of EBV-encoded circRNAs can potentially contribute to the development of virus-associated cancers through the circRNA-dependent alteration of cancer pathways.

Human-to-human transmission is pivot to viral evolution, which requires establishment of long-term infection. Thus, powerful immunoevasion strategies are always needed for the viruses to accomplish their goal of persistence. Since circRNAs are largely non-immunogenic, vcircRNAs are viewed as ideal signal molecules to fulfill many virus needs by regulating host pathways without eliciting the adaptive/innate anti-viral immune responses. This may partially explain why more abundant vcircRNAs are observed in herpesviruses-infected cells. Furthermore, the circRNA biogenesis cascade is likely preferentially suited to DNA viruses. Hence, on one hand, the nuclear replication of these viruses allows them to easily access to many RBPs/splicesome components required for circRNA synthesis. On the other hand, the inherent bicistronic or polycistronic nature of viral transcripts also increases their chance of being backspliced by the splicesomes.

To date, vcircRNAs have been discovered in the transcriptomes of a number of human and non-human viruses including EBV, KSHV, HPV, HBV, MCV, rhesus macaque lymphocryptovirus (rLCV), murid herpesvirus 68 (MHV68), and rat polyomavirus 2 (RatPyV2) (Sekiba et al., 2018; Toptan et al., 2018; Ungerleider et al., 2018 Ungerleider et al., , 2019a Zhao et al., 2019; Abere et al., 2020a,b) . As shown in Table 2 , the herpesvirus family holds the record for the highest number of vcircRNAs. Intriguingly, EBV is the first oncogenic human virus reported to express its own circRNAs and is also one of the top vcircRNA encoders.

Two research teams almost simultaneously utilized the RNase R-sequencing approach to assess vcircRNA expression in a variety of EBV infected cell models and patient tumor samples. In RNase R-sequencing, ribodepleted RNAs are treated with RNase R to remove the linear RNAs before library preparation. Thus, the sensitivity of circRNA detection is improved by enriching the circRNA species. Together, more than 30 EBV-encoded circRNAs (i.e., virally-encoded exonic circRNA and Exon-Intron circRNA) were identified, based on the unique backsplice junctions (Toptan et al., 2018; Ungerleider et al., 2018) . Many of these vcircRNAs are expressed at comparable or higher levels to most cellular circRNAs in EBV-associated cancers including gastric cancer, NPC, BL, PTLD, and AIDS-associated lymphoma, indicating a potential functional significance of these vcircRNAs in EBV's pathogenesis (Toptan et al., 2018; Ungerleider et al., 2018; Huang et al., 2019; Liu Q. et al., 2019; Qiao et al., 2019; Gong et al., 2020) .

These identified vcircRNAs can be further grouped into 4 clusters. Below, we will discuss each cluster based on their parental gene loci.

A total of sixteen vcircRNAs were derived from the RPMS1 gene locus within the BamHI A rightward transcript (BART) region (Toptan et al., 2018; Ungerleider et al., 2018) . The viral BART gene was initially discovered in EBV positive NPC cells (Hitt et al., 1989) . Later on, BART transcripts were observed in both healthy viral carriers and patient with EBV-associated malignancies such as BL, HD, and gastric cancer (Deacon et al., 1993; Chiang et al., 1996; Sugiura et al., 1996; Tao et al., 1998; Webster-Cyriaque and Raab-Traub, 1998) . In General, BART region is highly conserved among EBV strains except the B95-8 laboratory EBV strain has a 12-kb deletion of BART locus (Raab-Traub et al., 1980) . BART transcripts are a set of highly diverse transcripts including long noncoding RNAs such as RPMS1 and more than 40 EBV viral miRNA spliced from the RPMS1 intronic sequences (Edwards et al., 2008; Yamamoto and Iwatsuki, 2012; Cao et al., 2015b) . Both RPMS1 and viral miRNAs are involved in regulating cancer-related signaling pathways.

Four vcircRNAs derived from this locus were observed in all the examined latency types and virally infected cancer biopsies, but not in the B95-8 EBV infected cells due to a BART/RPMS1 locus deletion. Among them, circRPMS1_E4_E2_E3a_I3a_E3b (circBART_1.1) and circRPMS1_E4_E2_E3a_E3b (circBART_1.2) were formed by backsplicing of RPMS1 exon 4 to exon 2. Whereas, circRPMS1_E4_E3a_I3a_E3b (circBART_2.1), circRPMS1_E4_E3a_E3b (circBART_2.2) were formed by backsplicing of RPMS1 exon 4 to exon 3a. Interestingly, the intron-retained vcircRNAs circBART_1.1 and circBART_2.1 show nuclear distribution. In contrast, the exonic circBART_1.2 and circBART_2.2 are localized in both cytoplasm and nucleus.

Three vcircRNAs are derived from the EBNA latency locus (Toptan et al., 2018; Ungerleider et al., 2018) . Among them, circEBNA_W1_C1 and circEBNA_W2_C1 are formed by backsplicing of Bam HI W repeat W1 and W2 exons to the C1 exon. The expression of EBNA C1 exon is driven by the Bam HI C promoter (Cp) during the type III latency. Thus, these 2 vcircRNAs were mainly detected in type III latency cells with active Cp activity, but not in the Type I or II latency cells with epigenetically-silenced Cp. Meanwhile, circEBNA_U is generated by backsplicing of the Bam HI U fragment exon. Since the expression of EBNA U exon can be driven by all three EBNA promoters Cp, Wp, and Qp during both the latency and lytic cycles, the backsplicing of exon U was detected in all the examined type I and type III latency cells as well as cells harboring lytic reactivation.

A total of eight vcircRNAs were derived from the LMP2 locus (Toptan et al., 2018; Ungerleider et al., 2019b) . LMP2 is mainly a Type II/III latency protein that is expressed in a wide array of EBV associated B-and epithelial cell tumors where it likely contributes to the tumor phenotype (Yuan et al., 2006) . Notably, the LMP2 protein is also expressed during the lytic stage, indicating a potential lytic function of these LMP2 products Yuan et al., 2006; Concha et al., 2012) . Furthermore, previous studies have shown that LMP2 locus undergoes extensive alternative splicing to generate various linear transcript isoforms during reactivation (Concha et al., 2012; O'Grady et al., 2014 O'Grady et al., , 2016 . Among circLMP2s, high levels of circLMP2_E8_E2 were detected in reactivated Akata BL cells. This vcircRNA was generated by backsplicing of exon 8 to exon 2 of the LMP2 gene. Surprisingly, although Type III latency cells also constitutively express LMP2, circLMP2s were barely detected. It's postulated that certain tissue specific splicing factors are likely involved in the circLMP2 biogenesis.

circBHLF1 is derived from the intra-exonic backsplicing of the BHLF1 gene (Toptan et al., 2018; Ungerleider et al., 2018) . Increasing evidence has indicated that BHLF1 can be expressed in both latency and lytic cycles and may work as a non-coding gene for various EBV strains (Lin et al., 2010 (Lin et al., , 2013 Yetming et al., 2020) . The expression level of circBHLF1 is largely consistent with the isogenic linear transcript and it can be detected in most of the latency and reactivation conditions examined.

Notably, although most of the vcircRNAs were encoded by the latency genes, many of these vcircRNAs can be detected in both latency and reactivation conditions. Nevertheless, higher levels of vcircRNAs were detected in the lytic reactivation. Liang et al. (2017) have shown that cell stress and associated depletion of spliceosomal components can increase production of circRNA by re-directing nascent RNAs into circRNA biogenesis pathways. Since reactivation elicits strong stress to the host transcription machinery, it's possible that EBV might utilize similar strategy to promote vcircRNA production. However, the level of the backsplicing of cellular genes is not induced during reactivation (Ungerleider et al., 2018) . It suggests that the expression bias is likely due to virus-related mechanisms but not simply resulting from attenuated linear canonical splicing induced by the reactivation-related stress responses. Among the candidate viral factors regulating the biogenesis/transportation of these vcircRNAs, BMLF1 likely attracts the most attentions. BMLF is a viral lytic gene and it is actively involved in splicing, polyadenylation, RNA transport, and translation by interacting with various host factors including splicing regulators (Semmes et al., 1998; Hiriart et al., 2005; Sergeant et al., 2008; Juillard et al., 2009 Juillard et al., , 2012 Mure et al., 2018) . However, whether BMLF1 can indeed regulate the biogenesis of these vcircRNAs remains to be explored in the future.

Epstein-Barr virus circRNAs likely perform a wide spectrum of nuclear and cytoplasmic functions to facilitate viral pathogenesis. Given the universal expression of circRPMS1 in EBV-associated tumor tissues/cells as well as strong evolutionary conservation between EBV circRPMS1 and its rLCV homologue, it's reasonable to hypothesize that circRPMS1 may play an important role to facilitate EBV pathogenesis and associated tumor development (Ungerleider et al., 2019a) .

Indeed, Liu Q. et al. (2019) found higher levels of circRPMS1 in the metastatic NPC tissues collected from patients with a shorter overall survival. Further, the authors provided both in vitro and in vivo evidence showing that circRPMS1 likely functions as key oncogene to facilitate NPC development. By sequestering/sponging a number of host miRNAs including miR-203, miR-31 and miR-451, circRPMS1 can promote proliferation, apoptosis-resistance, and epithelial-mesenchymaltransition (EMT) of the tumor cells and thus enhance the oncogenic processes (Liu Q. et al., 2019) .

In another work, Huang et al. utilized computational approaches to identify a number of human and EBV miRNAs that may be directly sequestered by circRPMS1 . Some potential miRNA candidates have been shown to negatively regulate tumor development (e.g., hsa-miR-28-5p) or modulate EBV lytic cycles (e.g., ebv-miR-BART20-5p). Moreover, enhanced expression of circRPMS1 in a gastric cancer cell model significantly downregulated host miRNA candidates and increased cell motilities. This evidence supports an oncogenic role of circRPMS1.

Recently, Qiao et al. (2019) did a more comprehensive bioinformatics analysis of miRNAs candidates sponged by the EBV circRNAs including circRPMS1. They predicted that these targeted miRNAs are likely involved in herpesvirus infection, cell cycle regulation and various oncogenic pathways.

Meanwhile, since the backsplicing sites for circRPMS1_E4_E3a_I3a_E3b (circBART_2.1) and circRPMS1_E4_E3a_E3b (circBART_2.2) are in close proximity to the intron encoding EBV miRNAs, it's postulated that the backsplicing is likely to regulate these miRNA biosynthesis. Although circRPMS1 may still function as protein-coding RNAs in certain circumstances (e.g., in the in vivo environment), they are not associated with translated polysome fractions in the examined cells (Toptan et al., 2018) .

Given that both circEBNA_W1_C1 and circEBNA_W2_C1 are enriched in the nucleus, they are less likely to function as proteincoding transcripts. Further, the restricted expression of these vcircRNAs in Type III latency suggests that they may exhibit some latency-specific regulatory functions. For instance, they may be involved in the regulation of Cp activity and/or Cp-associated transcript diversity.

The evolutionary conservation of circEBNA_U between EBV and rLCV as well as similar expression of homologous rLCV_circEBNA_U in monkey lymphoma samples highly suggest a functional importance of this vcircRNA in viral pathogenesis (Ungerleider et al., 2019a) . Interestingly, circEBNA_U contains an open reading frame (ORF) encoding a short 6-amino acid peptide. Moreover, ribosome profiling analysis showed possible translation initiation from this ORF (Bencun et al., 2018) . With a potential IRES (internal ribosome entry site) translation initiation site, circEBNA_U may indeed encode some short peptides to regulate viral pathogenesis.

Generally, LMP2 gene encodes two major isoforms, LMP2A and LMP2B. LMP2A is encoded by exons 1 to 9, and functions as a transmembrane signaling molecule (Longnecker and Kieff, 1990) . In contrast, LMP2B is encoded by exons 2 to 9 (Longnecker and Kieff, 1990) . Without the first exon encoding a cytoplasmic signaling domain, LMP2B mainly functions as a dominant-negative transmembrane inhibitor of LMP2A (Rovedo and Longnecker, 2007) . Structure analysis shows that the circLMP2_E8_E2 likely carries all the coding sequence of the LMP2B. Since circLMP2_E8_E2 shows a significant cytoplasmic distribution, it's likely that this vcircRNA could function as a novel isoform of LMP2B to regulate LMP2A's function. In addition, the active expression of circLMP2_E8_E2 during the lytic cycle also indicates a potential role for this vcircRNA in reactivation.

Having analyzed the highly malignant SNU-4th cells isolated from the gastric cancer xenografts, Gong et al. provided evidence demonstrating that the EBV-encoded circLMP2A_E5_E3 likely functions as a miRNA sponge to induce and maintain cancer stemness by suppressing the anti-tumor effects of the miR-3908/TRIM59/p53 signaling, which is likely responsible for the metastasis and poor prognosis of EBV-associated gastric cancer .

Cell fraction analysis shows that the circBHLF1 is mainly localized in the nucleus, indicating a nuclear function of this vcircRNA. Interestingly, previous studies have shown that BHLF1 linear transcript is associated with the lytic replication complex and promotes lytic viral DNA replication (Rennekamp and Lieberman, 2011) . Thus, circBHLF1 may similarly interact with the proximal oriLyt DNA sequences and regulate viral DNA replication. The expression of circBHLF1 during the latency indicates that this vcircRNA may exhibit certain latency-related function as well.

circBHLF1 also carries an ORF encoding a putative 200 amino acid peptide. Evidence of translation initiation at this genomic location was reported previously (Bencun et al., 2018) . However, it still remains unclear if the low abundance of circBHLF in the cytoplasm can produce this putative peptide, or if the product is functional.

Kaposi's sarcoma-associated herpesvirus is also known as human herpesvirus 8 and is causatively associated with a number of human malignant diseases including Kaposi's sarcoma (KS), primary effusion lymphoma (PEL), and multicentric Castleman's disease (MCD) (Chang et al., 1994; Cesarman et al., 1995 Cesarman et al., , 1996 Soulier et al., 1995) . As a human γ-herpesvirus, KSHV also has two distinct stages, latency and lytic phase in its life cycle. During the latency, relatively restricted sets of viral genes are expressed in order to minimize the mounted host adaptive immunity (Staskus et al., 1997; Katano et al., 2000; Parravicini et al., 2000; Hosseinipour et al., 2014) . In the lytic cycle, more than 90 viral genes are expressed to facilitate the production of new infectious virions.

Like EBV, KSHV encodes its own non-coding linear RNAs to facilitate their needs for long-term persistence and transmission. At least 25 mature miRNAs have been identified in KSHV. This cluster of miRNAs is derived from a viral latency-associated genomic region and plays critical roles to support viral latency and associated oncogenesis (Cai et al., 2005; Pfeffer et al., 2005; Samols et al., 2005; Hansen et al., 2010; Lei et al., 2010a,b; Qin et al., 2010; Umbach and Cullen, 2010; Gottwein et al., 2011) .

KSHV also encodes a number of lncRNAs including PAN (polyadenylated nuclear) RNA, K7.3, T3.0, T1.2, and ALT (antisense-to-latency transcript) (Staskus et al., 1997; Sarid et al., 1998; Wang et al., 2004; Taylor et al., 2005; Chandriani et al., 2010; Majerciak et al., 2013; Arias et al., 2014; Rossetto and Pari, 2014; Schifano et al., 2017) . Most of these viral lncRNAs show increased expression during the lytic cycle and PAN RNA accounts for more than 80% of the viral transcripts during viral reactivation (Rossetto and Pari, 2014) . Further, PAN RNA is also the most thoroughly studied KSHV lncRNA and it's involved in transcriptional and epigenetic regulation of host and viral gene expression Pari, 2011, 2012; Rossetto et al., 2013) . For instance, PAN facilitates late lytic gene expression in the nucleus and PAN may also modulate the host immune response to facilitate viral infection (Borah et al., 2011; Pari, 2011, 2014; Rossetto et al., 2013; Campbell et al., 2014a,b) . Intriguingly, cytoplasmic PAN transcripts can be actively loaded into KSHV virions and likely function as immediate early RNA in the newly infected cells (Bechtel et al., 2005) . K7.3 RNA is derived from the opposite strand to PAN. K7.3 RNA is also increased during lytic cycle and it's predicted to exhibit some non-coding function (Dresang et al., 2011; Schifano et al., 2017) .

Recently, several groups have reported that KSHV encoded a repertoire of circRNAs in virally-infected cell lines and clinical specimens including primary KSHV-associated tumor tissues and blood samples (Tagawa et al., 2018; Toptan et al., 2018; Ungerleider et al., 2019a; Abere et al., 2020a) . These vcircRNAs are mainly derived from 2 genomic loci. In addition, a number of low-abundant vcircRNAs (such as kcirc3, kcirc38, kcirc54, kcirc55, kcirc57, and kcirc97) derived from several other genomic regions were also detected in certain cell model systems (Tagawa et al., 2018) .

More than 100 low-abundance vcircRNAs are derived from the PAN/K7.3 loci (Tagawa et al., 2018; Toptan et al., 2018; Ungerleider et al., 2019a; Abere et al., 2020a) . These vcircRNAs can be further categorized into two distinct clusters: circPAN and circK7.3 based on their parental sense (PAN) as well as antisense (K7.3) transcripts, respectively. Although the level of individual vcircRNA within each cluster is low, the accumulated abundance of the circPAN and circK7.3 is very high and they become a major component of the KSHV circRNAome. Both circPAN and circK7.3 were detected in the analyzed PEL cell lines regardless of KSHV life cycle, as well as KSHV-associated tumor tissues (e.g., KS and MCD) and blood samples. Further, both circPAN and circK7.3 are distributed in both cytoplasm and nucleus.

vIRF Locus circ-vIRF4s are derived from the vIRF4 (viral interferon regulatory factor 4)/K10 gene locus. The parental gene vIRF4 together with three other vIRFs namely vIRF1, vIRF2, and vIRF3 locate in a 10 kb genomic segment. The vIRF4 is homologous to cellular IRF4 and can inhibit the cellular IRF4 and it's also involved in cell cycle control and oncogenic process (Lee et al., 2009; Jacobs and Damania, 2011) .

circ-vIRF4s are generated by backsplicing of exon 2 to exon 1. Two isoforms are observed. circ-vIRF4 is a 530 bp exonic circRNA, while circ-vIRF4.IR is a 632 bp Exon-Intron circRNA. Differential subcellular distributions were observed for the two circ-IRF4 isoforms. The intron retention form (circ-vIRF4.IR) is predominantly enriched in the nucleus. The exononly form (circ-vIRF4) is mainly located in the cytoplasm. RNA decay analysis showed that the half-life of circ-vIRF4 (∼5 h) is significantly longer than the one of linear vIRF4 mRNA (∼2 h). circ-vIRF4s can be detected in PEL cell lines, KSHV-associated tumor tissues (e.g., KS and MCD) and blood samples.

The overlapping nature of circPAN and circK7.3 makes the characterization work more challenging. None of those backsplicing junctions are utilized to splice the known linear RNA transcripts from the PAN/K7.3 loci. To date, the mechanism for such hypervariable and bidirectional synthesis of circPAN/K7.3 remains unclear. Notably, circPAN and circK7.3 can be detected in both latency and lytic reactivation. After lytic induction, the levels of circPAN and circK7.3 are increased in parallel with the linear PAN and K7.3 transcripts. Thus, the synthesis of circPAN/K7.3 is likely at least partially regulated by the transcription initiation of their parental genes.

circ-vIRF4s are constantly expressed in both latency and lytic cycle. Although the parental gene vIRF4 is induced after lytic induction, the levels of circ-vIRF4s are barely increased or even decreased in some cell models. In certain cell models (e.g., BCP-1 and BBG1), higher levels of circ-IRF4s were observed than the isogenic linear transcript. This indicates that the observed latent circ-vIRF4s are regulated by an independent mechanism compared to the biosynthesis of linear vIRF4 transcript. For instance, the expression of circ-vIRF4s may be regulated by an unrecognized latency promoter.

To date, the roles of these vcircRNAs in KSHV's pathogenesis are still largely unclear. However, given the important roles of the parental genes of these vcircRNAs (as discussed earlier) and their (constitutive) expression profiles, it's postulated that these non-polysome-associated vcircRNAs are most likely to function as non-coding regulatory RNAs.

Intriguingly, circPAN, circK7.3, and cytoplasmic circ-vIRF4 (but not the nucleus circ-vIRF4.IR) can be actively loaded into the KSHV infectious virions. It's hypothesized that these vcircRNAs, just like viral tegument proteins carried by the progeny virions (Sathish et al., 2012) , may function as immediate-early gene products and host-immune modulators to help establish the infection in the newly infected cells.

Recently, Tagawa et al. showed that enhanced expression of some individual vcircRNAs such as kcirc54, kcirc55, and kcirc97 suppressed the transcription of viral lytic gene RTA and also induced an altered growth behaviors. This suggests a potential role of these vcircRNAs in KSHV pathogenesis (Tagawa et al., 2018) .

Human papillomaviruses are small non-enveloped, doublestranded DNA viruses (Brentjens et al., 2002; Roden and Stern, 2018) . To date, more than 200 types of HPV have been identified and they show epitheliotropic and mainly infect mucosal and cutaneous epithelia (Roden and Stern, 2018) . Based on the clinical consequence of the viral infection, HPV can be classified as low-risk HPVs (LR-HPVs) that generally cause benign lesions such as warts, as well as high-risk HPVs (HR-HPVs) that can transform the infected cells and promote cancer development (Roden and Stern, 2018) . Among HR-HPVs, HPV16, and HPV18 are the most prevalent high-risk types and are etiologically associated with a large number of human cancers that occur in the cervical, anogenital, and head and neck regions. The recent development and implementation of anti-HPV vaccines successfully help to reduce new infection caused by several HR-HPV and LR-HPV types, but unfortunately do not have any effect on the existing persistent infection (Roden and Stern, 2018) . Since the development of HPV-associated cancers is a chronic process, it is predicted that the HPV-associated malignancies will continuously cause significant medical problems in the next 50 years at least.

Human papillomavirus has a small ∼8 kb DNA genome which carries both the early (E) gene loci and late (L) gene loci (Brentjens et al., 2002) . The E loci are regulated by the early-promoter and they encode proteins (E1-E7) required for viral replication. Meanwhile, the late-promoter-regulated L loci encode the structural proteins (L1-L2) for progeny virion assembly (Brentjens et al., 2002) . HR-HPV E6 and E7 are the two most characterized viral oncogenes. They can inhibit major tumor suppressor proteins such as p53 and Rb and promote cancer development (Brentjens et al., 2002) . So far, HPV has not been reported to encode any linear non-coding transcripts such as lncRNA and miRNA.

Recently, a number of HPV-encoded circRNAs have been detected in HR-HPV-infected cell model systems and clinical cancer samples Zhao et al., 2019) . These vcircRNAs are generated by backsplicing of various splice sites located in both the E and L loci. Among these vcircRNAs, the most abundant form is known as circE7 which accounts for ∼1-3% of total E7 transcript and it is first identified in the HPV16 infected cells and cancer tissues. Analogous circE7 has been detected in other types of HR-HPVs including HPV18, 31, 33, 35, 45, and 58 . This 472 bp vcircRNA is generated by joining a downstream splice donor in E1 gene to an upstream splice acceptor in E6 gene. circE7 also carries the entire ORF of E7.

Current evidence shows that the biogenesis of circE7 is likely partially regulated by its m 6 A (N 6 -methyladenosine) sites . Indeed, mutation of circE7's m 6 A sites leads to a decreased level of circE7 transcripts, indicating that these sites may be targeted by splicing factors required for circRNA biogenesis.

The identification of multiple candidate miRNA binding sites on circE7 suggests that this vcircRNA may function as a miRNA sponge. Further, knockdown of HPV16 circE7 leads to an increased level of the linear E6/E7 transcripts, suggesting that circE7 may regulate the biogenesis of its isogenic transcripts by competing for their shared splice sites.

Besides its potential non-coding function, circE7 is likely to function as a protein-coding RNA, since it is enriched in the cytoplasm and also associated with translated polysome fractions. circE7 has been shown to produce a full length E7 protein in various cell models. The translation is apparently initiated through an m 6 A (N 6 -methyladenosine)-mediated capindependent mechanism. HPV16 circ-E7-mediated expression of E7 protein has been shown to regulate cell proliferation, anchorage-independent cell growth, as well as oncogenesis in both in vitro cell models and in vivo tumor xenografts. Further, the expression of circ-E7 in virally infected non-cancer cells indicates that this vcircRNA may also play a role in viral infection cycle.

Merkel cell polyomavirus is causally associated with a majority of the Merkel cell carcinoma, an aggressive skin cancer mainly affecting elderly and immunocompromised individuals Shuda et al., 2008) . MCV is a nonenveloped virus with a circular double-stranded DNA genome and belongs to the human polyomavirus family Shuda et al., 2008) . In the infected cells, the MCV genome can remain either as a free episomal DNA undergoing active replication, or as a replication-deficient form after integrating into the host genome Shuda et al., 2008) . It is believed that genome integration plays a critical role in MCVmediated oncogenesis Shuda et al., 2008) . Its relatively small DNA genome can be simply divided into two gene loci, the early region (ER) and late region (LR). Generally, the ER-encoded gene products such as large T antigen (LT) are involved in the regulation of viral replication and oncogenesis. For instance, the LT can directly target tumor suppressor protein Rb. Meanwhile, it can bind to the viral origin of replication and initiate viral genome replication with its intrinsic helicase/ATPase activity Shuda et al., 2008) . In contrast, the LR loci encode viral structural proteins such as VP1 and VP2 required for new virion production (Tolstov et al., 2009) .

Besides viral protein products, MCV also encodes two miRNAs, namely MCV-miR-M1-5p and MCV-miR-M1-3p, which are derived from the pri-miRNA transcript synthesized from the negative strand of the ER (Seo et al., 2009) . These 22 nucleotide viral miRNAs have been shown to inhibit MCV lytic replication by knocking-down the LT transcripts in order to facilitate maintenance of the episomal genome (Seo et al., 2009; Theiss et al., 2015) .

Recently, MCV has been shown to encode four circRNAs (Abere et al., 2020b) . Among them, three vcircRNA are derived from the ER loci and one is encoded by the LR loci. circMCV-T (circMCV_1622_861) is the most abundant one and accounts for more than 70% of the total vcircRNA population. circMCV-T is a 762 bp vcircRNA carrying the entire exon II of the T antigen. It's generated by backsplicing of the downstream 5 splice donor (nt 1622) to the upstream 3 splice acceptor (nt 861). These splice sites are also required for the splicing of all the linear T antigen isoforms including LT, small T antigen (sT), 57,000-molecular-weight T antigen (57 kT), and alternative LT open reading frame (ALTO) transcripts. circMCV_2955_861 and circMCV_3337_861 use the same 3 splice acceptor sites but different 5 splice donor sites at nt 2955 and 3337, respectively. circMCV_3913_3271 shows lowest abundance and is derived from the LR loci. To date, among these vcircRNAs, only the most abundant circMCV-T was further characterized.

The circMCV-T transcripts are detected mainly in the cells harboring actively replicating episomal viral genomes, but rarely in the cells carrying more dormant integrated viral genome. This suggests that the biogenesis of circMCV-T is at least partially regulated by the transcription activity of its parental ER loci.

A circMCV-T homologue, namely circRatPyV2_4112_4468 was recently detected in the RatPyV2 infected cells (Abere et al., 2020b) . The evolution conservation of circMCV-T further indicates a functional significance of this RNA molecule. Increasing evidence has shown that the MCV-encoded transcripts including circMCV-T, miR-M1, and linear T antigen transcripts form a dynamic regulatory circuit to control MCV's life cycle. Due to a 100% sequence complementarity between miR-M1 and linear T antigen mRNAs, the miRNA-mRNA interaction can knock down linear T antigen transcripts and inhibit viral DNA replication. The presence of the new player circMCV-T further increases the complexity of the network. Due to a 100% sequence complementarity between miR-M1 and circMCV-T, the reciprocal interaction between circMCV-T and miR-M1 can lead to two scenarios. In the first scenario, miR-M1 can decrease the stability of circMCV-T by linearizing it and thus reduce the genome replication. In the second, circMCV-T can increase the level of linear T antigen transcripts and facilitate viral replication by sequestering miR-M1. Hence, the overall effects will be determined by the relatively amount of miR-M1 and circMCV-T: a condition of high circMCV-T and low miR-M1 will favor active MCV replication. Conversely, a low circMCV-T and high miR-M1 condition will suppress active viral replication and promote a long-term persistent infection.

Besides the non-coding function, the circMCV-T transcript may potentially function as a protein-coding RNA. The in situ staining data show that circMCV-T is predominantly located in the cytoplasm. It also carries the entire exon 2 of T-antigen as well as the start codon for the ALTO ORF. However, some data show that circMCV-T transcripts are not associated with any translated polysomes in certain cell models.

Hepatitis B virus is a small oncogenic human virus closely associated with liver cirrhosis and HCC (Torresi et al., 2019) . HBV has a 3.2 kb partially double-stranded relaxed circular DNA (rcDNA) genome. In the infected cells, a covalently closed circular DNA (cccDNA) genome is formed and remains as a persistent episomal genome. The cccDNA is used as templates to generate key transcripts including the pregenomic RNA (pgRNA) transcript for viral replication.

An HBV-encoded circRNA has been recently identified in both HBV-infected cell models and HCC tissues (Sekiba et al., 2018; Zhu et al., 2020) . Zhu et al. (2020) has shown that this 2.5 kb vcircRNA, namely HBV_circ_1 is derived by the intronless pgRNA. Surprisingly, instead of utilizing the backsplicing, it is believed that this vcircRNA is generated by a mechanism involving the homologous recombination of the inverted repeat sequences at both 3 and 5 ends of the pgRNA (Zhu et al., 2020) . In addition, HBV_circ_1 is predominantly located in the cytoplasm (Zhu et al., 2020) . Sekiba et al. (2018) showed that knockdown of an RBP DXH9 increases the levels of HBV_circ_1. It is proposed that the expression of HBV_circ_1 can be negatively regulated by the direct interaction between DHX9 and the inverted repeat regions of the pgRNA. To date, the biological function of HBV_circ_1 remains unclear.

To date, detection of viral genomic material has been widely implemented in the clinical practice. For instance, circulating EBV DNA has been used as a good biomarker for the diagnosis and prognosis of EBV-associated NPC (Lo et al., 1999; Leung et al., 2006; Chan et al., 2017) . Recently, circRNAs have also received increasing attention for its potential clinical applications. These highly stable vcircRNAs that exist in various types of body fluids and tissues may serve as good diagnostic and prognostic markers for many virally associated diseases.

Current evidence showed that vcircRNAs might be superior to the current RNA diagnostic marker. For EBV, EBERs are used as a gold standard for diagnosis of EBV infection (Ambinder and Mann, 1994) . However, data showed that circRPMS1 were detected in certain PTLD tissues carrying reasonable number of EBV genomes, where EBER cannot be detected, indicating potential false-negative results coming from the EBER diagnostic test (Toptan et al., 2018) . For KSHV, circPAN/K7.3 and circ-vIRF4 showed a higher detection rate than the linear viral transcript LANA (latency-associated nuclear antigen) in the analyzed primary KS tumor and associated blood samples (Abere et al., 2020a) . Lastly, in the case of HPV, PCR-based detection of HR-HPV has been used for cervical cancer screening. The current data show that circE7 can be detected in fixed archive tissue samples as well . It is warranted to investigate if circE7 may serve as a promising biomarker for HPV infection and HPV-related cancers.

Due to the recent development of circRNA enrichment and high-throughput detection technologies, we are now in a rapidly growing era of vcircRNA research. Numerous novel vcircRNAs have been identified in human oncogenic viruses, and they are no longer viewed as insignificant splicing noise. Instead, vcircRNAs are now seen for what they are: an important new class of regulatory molecules in viral infection cycle and associated malignant diseases. The addition of this new member has created an entirely new layer of complexity of the viral pathogenic mechanism.

Nevertheless, the vcircRNA research is in its nascent stage where many puzzles remain to be solved. Currently, the underlying mechanism for the biogenesis of vcircRNAs is still largely unclear. It's postulated that, in addition to the general splicing machinery, both tissue-specific host factors and viral factors are likely to be involved in the formation of vcircRNAs. For instance, EBV circLMP2_E8_E2 can be detected in reactivated BL cells, but not in the type III latency cell lines that constitutively express parental linear LMP2A transcripts. It indicates that some tissue-specific factors may regulate the backsplicing by targeting the cis regulatory elements in the primary LMP2 transcript (Ungerleider et al., 2018) . Meanwhile, EBV viral lytic gene BMLF1 is known to play critical roles in the maturation of viral linear transcripts and likely contribute to the backsplicing process. Future research is warranted to elucidate the detailed mechanism of vcircRNA synthesis.

Notably, within the family of human oncogenic viruses, it is still unclear if HCV or human T-lymphotropic virus (human T-cell leukemia virus type 1 or HTLV-1) encodes any vcircRNAs. Future research is surely warranted to investigate if any vcircRNAs are present in the HCV or HTLV-1positive biospecimen. The RNase R-sequencing will be an ideal approach to evaluate the vcircRNAome in a high-throughput and unbiased manner.

Among all the human oncogenic viruses investigated so far, herpesviruses are undoubtedly the largest producer of vcircRNAs and account for >85% of the identified vcircRNAs. This expression bias may be partially attributable to herpesviruses' unique bi-phasic life cycles: they can exist in both latent and lytic phases. Particularly during the latency, to escape the host immune surveillance, herpesviruses cannot express antigenic viral proteins, while they still need to regulate themselves and host environment. Hence, expression of noncoding RNAs which lack immunogenicity is likely an ideal strategy for the viruses to regulate themselves and the host environment when an immunologically transparent state is required. Increasing evidence has shown that vcircRNAs are less immunogenic, compared to linear transcripts. Chen et al. reported that unmodified exogenous circRNAs (without any nucleoside modifications) are recognized as "nonself " RNAs by the cellular sensors such as RIG-I and thus trigger the innate immune signaling. However, the circRNAs spliced by endogenous human spliceosomes are viewed as "self " RNAs and lose the immunogenicity, due to the N 6 -methyladenosine (m 6 A) modification during the biogenesis (Chen Y.G. et al., 2017; Chen et al., 2019) . Based on the high prevalence of m 6 A modification of viral transcripts and evidence that vcircRNAs are largely generated by the same host splicesomes used to synthesize "self " circRNAs, it is reasonably to believe that vcircRNAs are non-immunogenic (Wu F. et al., 2019; Lu et al., 2020) . Notably, Wesselhoeft et al. recently reported a discordant observation that immunogenicity of circRNA is independent of m 6 A modification. They found that unmodified exogenous circRNAs are also non-immunogenic and can avoid antiviral defense induction upon cellular entry by dodging immune sensors such as RIG-I and Toll-like receptor (Wesselhoeft et al., 2019) . Although this new discovery still supports the notion that vcircRNAs are non-immunogenic, more efforts are warranted to investigate the discrepant data on the mechanism of immunoevasion of foreign circRNAs.

To date, only a handful of identified vcircRNAs are functionally investigated. A majority of these characterized vcircRNAs are enriched in the cytoplasm and are proposed to function as either miRNA sponges or protein-coding transcripts to regulate viral pathogenesis. These conclusions are largely based on the data from the in vitro cell models. However, the functional impacts of the large number of nuclear vcircRNAs have not yet been reported. These nuclear vcircRNAs such as EBV circRPMS1_E4_E2_E3a_I3a_E3b (circBART_1.1) and circRPMS1_E4_E3a_I3a_E3b (circBART_2.1) generally carry retained-introns from their parental genes. Due to their high abundance and nuclear distribution, they may play important regulatory roles such as modulating the levels of their isogenic transcripts by competing the splice sites and limited splicing factors. Further investigation is surely warranted.

The detailed transcript structures of these vcircRNA candidates are largely unclear. The size and sequence of most of these vcircRNAs are predicted simply based on their backsplicing sites. In-depth sequencing-based assembly and validation of the vcircRNAs is surely warranted, which is essential for the subsequent characterization work to understand the biogenesis and function of each vcircRNA isoform. It will be interesting to see if some new computation pipelines designed for reconstruction of full-length circRNA can be successfully implemented in vcircRNA research (Wu J. et al., 2019; Zheng et al., 2019) .

To date, all the reported vcircRNAomes of these human oncogenic viruses are generated by RNA-seq analyses of bulk cell cultures and/or clinical tissues. Due to the inevitable heterogeneity of tissue/cell culture, the vcircRNAomes are likely different at the single-cell level. Thus, the expression and function of these vcircRNAs might vary among single cells. Implementation of single-cell sequencing technologies which are compatible with circRNAs might help us better decipher the biogenesis mechanism and function of vcircRNAs (Fan et al., 2015) .

Our understanding of the biogenesis mechanism and function of these vcircRNAs are largely based on the data derived from in vitro cell cultures. These data need to be further validated and examined in more physiological relevant in vivo model systems such as animal models and clinical biospecimens.

The altered expression of vcircRNAs in infected host cells and diseased tissues as well as the long half-life of vcircRNAs make them potential biomarkers for clinical diagnosis of viral infection and prognosis of disease outcomes. Some of the oncogenic vcircRNAs such as circLMP2A might serve as potential targets for precision therapy. Given the successful implementation of mRNA-based vaccines to prevent SARS-CoV-2 infection, future highly stable circRNA-based vaccines might become a new tool to better prevent human oncogenic virus infection and related diseases.

The discoveries of viral miRNAs and lncRNAs and their roles have dramatically increased our knowledge of how oncogenic viruses interact with hosts and facilitate cancer development. Just like in the early years of viral miRNA and lncRNA research, the vcircRNA research has now started gaining increased attention and becoming an important topic in the viral oncology field. With further improvement of the circRNA analytical techniques, we expect that more vcircRNAs will be identified and successfully characterized in the coming years. These developments will help us better understand the roles of human oncogenic viruses in their associated cancers and improve patient prognosis in the fight against these deadly diseases.

ZL, JA, DB, and RA contributed to the conceptualization, drafting, and revision of the manuscript. AN, JK, and LL contributed to the conceptualization and manuscript revision. All authors contributed to the article and approved the submitted version.

Research reported in this publication was supported by the National Cancer Institute of the National Institutes of Health under Award Number R01CA261258 to ZL. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. This work was also supported by a National Institutes of Health COBRE grant P20GM121288, a U.S.-Japan Cooperative Medical Sciences Program Collaborative Award from the National Institute of Allergy and Infectious Diseases and CRDF Global (grant number DAA3-19-65602-1), a Tulane school of medicine faculty research pilot grant, and a Carol Lavin Bernick faculty grant to ZL as well as the Japan Society for the Promotion of Science (grant number JP17K4567), the Takeda Science Foundation, and the Shiseido Female Researcher Science Grant to AN.

Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1

Kaposi's sarcoma-associated herpesvirus-encoded circrnas are expressed in infected tumor tissues and are incorporated into virions

Merkel cell polyomavirus encodes circular RNAs (circRNAs) enabling a dynamic circRNA/microRNA/mRNA regulatory network

Altered expression pattern of circular RNAs in primary and metastatic sites of epithelial ovarian carcinoma

DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome

Epstein-barr-encoded RNA in situ hybridization: diagnostic applications

KSHV 2.0: a comprehensive annotation of the Kaposi's sarcomaassociated herpesvirus genome using next-generation sequencing reveals novel genomic and functional features

Some yeast mitochondrial RNAs are circular

circRNA biogenesis competes with pre-mRNA splicing

Dysregulated miR-671-5p / CDR1-AS / CDR1 / VSNL1 axis is involved in glioblastoma multiforme

RNAs in the virion of Kaposi's sarcoma-associated herpesvirus

Translational profiling of B cells infected with the Epstein-Barr virus reveals 5' leader ribosome recruitment through upstream open reading frames

A viral nuclear noncoding RNA binds re-localized poly(A) binding protein and is required for late KSHV gene expression

Human papillomavirus: a review

Kaposi's sarcoma-associated herpesvirus expresses an array of viral microRNAs in latently infected cells

virus microRNAs are evolutionarily conserved and differentially expressed

A lytic viral long noncoding RNA modulates the function of a latent protein

Long non-coding RNA and epigenetic gene regulation of KSHV

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

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

Kaposi's sarcoma-associated herpesvirus-like DNA sequences in AIDS-related body-cavity-based lymphomas

Kaposi's sarcoma-associated herpesvirus in non-AIDS related lymphomas occurring in body cavities

Assessment of circularized E7 RNA, GLUT1, and PD-L1 in anal squamous cell carcinoma

Analysis of Plasma Epstein-Barr virus DNA to screen for nasopharyngeal cancer

The lytic transcriptome of Kaposi's sarcoma-associated herpesvirus reveals extensive transcription of noncoding regions, including regions antisense to important genes

Identification of herpesvirus-like DNA sequences in AIDS-associated Kaposi's sarcoma

Circular RNA profile identifies circPVT1 as a proliferative factor and prognostic marker in gastric cancer

Using circular RNA hsa_circ_0000190 as a new biomarker in the diagnosis of gastric cancer

Sensing self and foreign circular RNAs by intron identity

Regulation of circRNA biogenesis

N6-methyladenosine modification controls circular RNA immunity

Nasal NK-and T-cell lymphomas share the same type of Epstein-Barr virus latency as nasopharyngeal carcinoma and Hodgkin's disease

Mis-splicing yields circular RNA molecules

Identification of new viral genes and transcript isoforms during Epstein-Barr virus reactivation using RNA-Seq

The RNA binding protein quaking regulates formation of circRNAs

Transcriptome-wide discovery of circular RNAs in Archaea

Epstein-Barr virus infection and risk of lymphoma: immunoblot analysis of antibody responses against EBV-related proteins in a large series of lymphoma subjects and matched controls

Epstein-Barr virus and Hodgkin's disease: transcriptional analysis of virus latency in the malignant cells

Coupled transcriptome and proteome analysis of human lymphotropic tumor viruses: insights on the detection and discovery of viral genes

Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2

Epstein-Barr virus BART microRNAs are produced from a large intron prior to splicing

FUS affects circular RNA expression in murine embryonic stem cell-derived motor neurons

Singlecell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos

Clonal integration of a polyomavirus in human Merkel cell carcinoma

Epstein-Barr virus-derived circular RNA LMP2A induces stemness in EBV-associated gastric cancer

Viral microRNA targetome of KSHV-infected primary effusion lymphoma cell lines

Comprehensive profile of differentially expressed circular RNAs reveals that hsa_circ_0000069 is upregulated and promotes cell proliferation, migration, and invasion in colorectal cancer

Expanded identification and characterization of mammalian circular RNAs

KSHV-encoded miRNAs target MAF to induce endothelial cell reprogramming

Natural RNA circles function as efficient microRNA sponges

Interaction of the Epstein-Barr virus mRNA export factor EB2 with human Spen proteins SHARP, OTT1, and a novel member of the family, OTT3, links Spen proteins with splicing regulation and mRNA export

EBV gene expression in an NPC-related tumour

Circular non-coding RNA ANRIL modulates ribosomal RNA maturation and atherosclerosis in humans

Viral profiling identifies multiple subtypes of Kaposi's sarcoma

Noncoding effects of circular RNA CCDC66 promote colon cancer growth and metastasis

Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells

cir-ITCH plays an inhibitory role in colorectal cancer by regulating the Wnt/beta-catenin pathway

Identification of virus-encoded circular RNA

The emerging role of circular RNAs in transcriptome regulation

Expression and processing of a small nucleolar RNA from the Epstein-Barr virus genome

Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals

The viral interferon regulatory factors of KSHV: immunosuppressors or oncogenes?

Detecting and characterizing circular RNAs

Circular RNAs are abundant, conserved, and associated with ALU repeats

Epstein-Barr virus protein EB2 stimulates cytoplasmic mRNA accumulation by counteracting the deleterious effects of SRp20 on viral mRNAs

Epstein-Barr virus protein EB2 contains an N-terminal transferable nuclear export signal that promotes nucleocytoplasmic export by directly binding TAP/NXF1

Expression and localization of human herpesvirus 8-encoded proteins in primary effusion lymphoma, Kaposi's sarcoma, and multicentric Castleman's disease

Detection of epstein-barr virus infection in non-small cell lung cancer

A network of noncoding regulatory RNAs acts in the mammalian brain

Isolation and characterization of sendai virus DI-RNAs

Combinatorial control of Drosophila circular RNA expression by intronic repeats, hnRNPs, and SR proteins

Circular RNAs are abundantly expressed and upregulated during human epidermal stem cell differentiation

Terminal differentiation into plasma cells initiates the replicative cycle of Epstein-Barr virus in vivo

Viral interferon regulatory factors

Circ-ZNF609 is a circular RNA that can be translated and functions in myogenesis

Regulation of herpesvirus lifecycle by viral microRNAs

Regulation of NF-kappaB inhibitor IkappaBalpha and viral replication by a KSHV microRNA

Plasma Epstein-Barr viral deoxyribonucleic acid quantitation complements tumornode-metastasis staging prognostication in nasopharyngeal carcinoma

Circular RNA expression profile of pancreatic ductal adenocarcinoma revealed by microarray

Circular RNA 0000096 affects cell growth and migration in gastric cancer

Decreased expression of Hsa_circ_00001649 in gastric cancer and its clinical significance

Using circular RNA as a novel type of biomarker in the screening of gastric cancer

Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis

Exon-intron circular RNAs regulate transcription in the nucleus

The output of protein-coding genes shifts to circular RNAs when the Pre-mRNA processing machinery is limiting

Whole-genome sequencing of the Akata and Mutu Epstein-Barr virus strains

Quantitative and qualitative RNA-Seq-based evaluation of Epstein-Barr virus transcription in type I latency Burkitt's lymphoma cells

Knockdown of EBV-encoded circRNA circRPMS1 suppresses nasopharyngeal carcinoma cell proliferation and metastasis through sponging multiple miRNAs

Circular RNA hsa_circ_001783 regulates breast cancer progression via sponging miR-200c-3p

Quantitative analysis of cell-free Epstein-Barr virus DNA in plasma of patients with nasopharyngeal carcinoma

A second Epstein-Barr virus membrane protein (LMP2) is expressed in latent infection and colocalizes with LMP1

Genome-wide transcription program and expression of the Rta responsive gene of Epstein-Barr virus

N(6)-methyladenosine modification enables viral RNA to escape recognition by RNA sensor RIG-I

A viral genome landscape of RNA polyadenylation from KSHV latent to lytic infection

Circular RNAs are a large class of animal RNAs with regulatory potency

Human viral oncogenesis: a cancer hallmarks analysis

Common commensal cancer viruses

Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA

Epstein-Barr virus protein EB2 stimulates translation initiation of mRNAs through direct interactions with both Poly(A)-binding protein and eukaryotic initiation factor 4G

Engineering and expressing circular RNAs via tRNA splicing

Global bidirectional transcription of the Epstein-Barr virus genome during reactivation

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

Translation of CircRNAs

An in vivo knockdown strategy reveals multiple functions for circMbl. bioRxiv

Differential viral protein expression in Kaposi's sarcoma-associated herpesvirus-infected diseases: Kaposi's sarcoma, primary effusion lymphoma, and multicentric Castleman's disease

Identification of microRNAs of the herpesvirus family

Identification of virus-encoded microRNAs

Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function

Epstein-Barr virus circRNAome as host miRNA sponge regulates virus infection, cell cycle, and oncogenesis

Hsa_circ_0001649: a circular RNA and potential novel biomarker for hepatocellular carcinoma

Upregulation of xCT by KSHV-encoded microRNAs facilitates KSHV dissemination and persistence in an environment of oxidative stress

DNA of Epstein-Barr virus VIII: B95-8, the previous prototype, is an unusual deletion derivative

Initiation of Epstein-Barr virus lytic replication requires transcription and the formation of a stable RNA-DNA hybrid molecule at OriLyt

Epstein-Barr virus

Opportunities and challenges for human papillomavirus vaccination in cancer

KSHV PAN RNA associates with demethylases UTX and JMJD3 to activate lytic replication through a physical interaction with the virus genome

Kaposi's sarcoma-associated herpesvirus noncoding polyadenylated nuclear RNA interacts with virus-and host cellencoded proteins and suppresses expression of genes involved in immune modulation

PAN's labyrinth: molecular biology of Kaposi's sarcoma-associated herpesvirus (KSHV) PAN RNA, a multifunctional long noncoding RNA

Regulation of viral and cellular gene expression by Kaposi's sarcomaassociated herpesvirus polyadenylated nuclear RNA

Epstein-barr virus latent membrane protein 2B (LMP2B) modulates LMP2A activity

Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed

Cell-type specific features of circular RNA expression

Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types

Cloning and identification of a microRNA cluster within the latency-associated region of Kaposi's sarcoma-associated herpesvirus

Circular RNA expression in cutaneous squamous cell carcinoma

Circular RNA expression in basal cell carcinoma

Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures

Transcription mapping of the Kaposi's sarcoma-associated herpesvirus (human herpesvirus 8) genome in a body cavity-based lymphoma cell line

Tegument proteins of Kaposi's sarcoma-associated herpesvirus and related gamma-herpesviruses

Expression of the antisense-to-latency transcript long noncoding RNA in Kaposi's sarcomaassociated herpesvirus

DHX9 regulates production of hepatitis B virus-derived circular RNA and viral protein levels

Mta has properties of an RNA export protein and increases cytoplasmic accumulation of Epstein-Barr virus replication gene mRNA

Merkel cell polyomavirus encodes a microRNA with the ability to autoregulate viral gene expression

The Epstein-Barr virus (EBV) protein EB is an mRNA export factor essential for virus production

Comprehensive circular RNA profiling reveals that hsa_circ_0005075, a new circular RNA biomarker. is involved in hepatocellular crcinoma development

The global landscape of EBVassociated tumors

T antigen mutations are a human tumor-specific signature for Merkel cell polyomavirus

Kaposi's sarcoma-associated herpesvirus-like DNA sequences in multicentric Castleman's disease

Kaposi's sarcoma-associated herpesvirus gene expression in endothelial (spindle) tumor cells

Epstein-Barr virus and human herpesvirus 6 detection in a non-Hodgkin's diffuse large B-cell lymphoma cohort by using RNA sequencing

Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy

Transcriptional analysis of Epstein-Barr virus gene expression in EBVpositive gastric carcinoma: unique viral latency in the tumour cells

Circular RNA and gene expression profiles in gastric cancer based on microarray chip technology

Plasma cell-specific transcription factor XBP-1s binds to and transactivates the Epstein-Barr virus BZLF1 promoter

Discovery of Kaposi's sarcoma herpesvirus-encoded circular RNAs and a human antiviral circular RNA

virus (EBV) in endemic Burkitt's lymphoma: molecular analysis of primary tumor tissue

Transcriptional analysis of latent and inducible Kaposi's sarcoma-associated herpesvirus transcripts in the K4 to K7 region

A comprehensive analysis of replicating merkel cell polyomavirus genomes delineates the viral transcription program and suggests a role for mcv-miR-M1 in episomal persistence

Human Merkel cell polyomavirus infection II. MCV is a common human infection that can be detected by conformational capsid epitope immunoassays

Circular DNA tumor viruses make circular RNAs

HBV-related hepatocarcinogenesis: the role of signalling pathways and innovative ex vivo research models

In-depth analysis of Kaposi's sarcomaassociated herpesvirus microRNA expression provides insights into the mammalian microRNA-processing machinery

The Epstein Barr virus circRNAome

Comparative analysis of gammaherpesvirus circular RNA repertoires: conserved and unique viral circular RNAs

Circular RNA-ITCH Suppresses Lung Cancer Proliferation via Inhibiting the Wnt/beta-Catenin Pathway

Decreased expression of hsa_circ_001988 in colorectal cancer and its clinical significances

Kaposi's sarcoma-associated herpesvirus ori-Lyt-dependent DNA replication: cis-acting requirements for replication and ori-Lyt-associated RNA transcription

Transcription of Epstein-Barr virus latent cycle genes in oral hairy leukoplakia

Circular RNA ciRS-7-a promising prognostic biomarker and a potential therapeutic target in colorectal cancer

RNA circularization diminishes immunogenicity and can extend translation duration in vivo

Viruses and human cancers: a long road of discovery of molecular paradigms

Association of N6-methyladenosine with viruses and related diseases

CircAST: full-length assembly and quantification of alternatively spliced isoforms in circular RNAs

Circular RNA HIPK3 plays a carcinogenic role in cervical cancer progression via regulating miR-485-3p/FGF2 axis

Circular RNA has_circ_0067934 is upregulated in esophageal squamous cell carcinoma and promoted proliferation

Emerging roles of circRNA_001569 targeting miR-145 in the proliferation and invasion of colorectal cancer

The circular RNA ciRS-7 (Cdr1as) acts as a risk factor of hepatic microvascular invasion in hepatocellular carcinoma

Circular RNA: a novel biomarker for progressive laryngeal cancer

Diversity of Epstein-Barr virus BamHI-A rightward transcripts and their expression patterns in lytic and latent infections

Circular RNA expression profiles in nasopharyngeal carcinoma by sequence analysis

Extensive translation of circular RNAs driven by N(6)-methyladenosine

ZKSCAN1 gene and its related circular RNA (circZKSCAN1) both inhibit hepatocellular carcinoma cell growth, migration, and invasion but through different signaling pathways

The BHLF1 locus of epstein-barr virus contributes to viral latency and B-cell immortalization

Epstein-Barr virus: more than 50 years old and still providing surprises

The circular RNA circBIRC6 participates in the molecular circuitry controlling human pluripotency

The circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 expression

Virus and cell RNAs expressed during Epstein-Barr virus replication

The landscape of viral associations in human cancers

RNA-Seq profiling of circular RNAs in human small cell lung cancer

Complementary sequence-mediated exon circularization

Circular intronic long noncoding RNAs

Transforming activity of an oncoprotein-encoding circular RNA from human papillomavirus

Clinical values of circular RNA 0000181 in the screening of gastric cancer

TTBK2 circular RNA promotes glioma malignancy by regulating miR-217/HNF1beta/Derlin-1 pathway

Reconstruction of full-length circular RNAs enables isoform-level quantification

HBV pgRNA can generate a circRNA with two junction sites. bioRxiv

Circular BANP, an upregulated circular RNA that modulates cell proliferation in colorectal cancer

The authors declare that the research was conducted in the