key: cord-0854509-mnnwaaoq
authors: Li, Jing; Teng, Peiying; Yang, Fan; Ou, Xia; Zhang, Jihong; Chen, Wei
title: Bioinformatics and Screening of a Circular RNA-microRNA-mRNA Regulatory Network Induced by Coxsackievirus Group B5 in Human Rhabdomyosarcoma Cells
date: 2022-04-22
journal: Int J Mol Sci
DOI: 10.3390/ijms23094628
sha: 7bd5f77b326ae81593758959e4f8f47e6b674fca
doc_id: 854509
cord_uid: mnnwaaoq

Hand, foot and mouth disease (HFMD) caused by Coxsackievirus Group B5 (CVB5) is one of the most common herpetic diseases in human infants and children. The pathogenesis of CVB5 remains unknown. Circular RNAs (CircRNAs), as novel noncoding RNAs, have been shown to play a key role in many pathogenic processes in different species; however, their functions during the process of CVB5 infection remain unclear. In the present study, we investigated the expression profiles of circRNAs using RNA sequencing technology in CVB5-infected and mock-infected human rhabdomyosarcoma cells (CVB5 virus that had been isolated from clinical specimens). In addition, several differentially expressed circRNAs were validated by RT-qPCR. Moreover, the innate immune responses related to circRNA-miRNA-mRNA interaction networks were constructed and verified. A total of 5461 circRNAs were identified at different genomic locations in CVB5 infections and controls, of which 235 were differentially expressed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis demonstrated that the differentially expressed circRNAs were principally involved in specific signaling pathways related to ErbB, TNF, and innate immunity. We further predicted that novel_circ_0002006 might act as a molecular sponge for miR-152-3p through the IFN-I pathway to inhibit CVB5 replication, and that novel_circ_0001066 might act as a molecular sponge for miR-29b-3p via the NF-κB pathway and for the inhibition of CVB5 replication. These findings will help to elucidate the biological functions of circRNAs in the progression of CVB5-related HFMD and identify prospective biomarkers and therapeutic targets for this disease.

Coxsackievirus Group B5 (CVB5) is a non-enveloped virus with linear single-stranded RNA that belongs to a member of the genus Enterovirus from the Picornaviridae family. The CVB5 genome is approximately 7400 bp long, with only one open reading frame (ORF), which ultimately constitutes of four structural proteins (VP1, VP2, VP3, and VP4) and seven non-structural proteins (2A, 2B, 2C, 3A, 3B, 3C, and 3D) [1] . CVB5 is one of the common pathogens causes of hand, foot and mouth disease (HFMD) [2, 3] . HFMD is generally considered to be a common global exanthemata and febrile illness that is mostly seen in children under 5 years old. This disease is characterized by a maculopapular rash or blisters on the hands, feet, groin, and buttocks and is associated with painful ulcerative lesions of the mouth [4] . Moreover, CVB5 is a neurotropic virus, tends to be more severe, and is more likely to produce severe complications in the central nervous system (CNS), thus causing infections such as aseptic meningitis or, acute brainstem encephalitis. Children who recover from brainstem encephalitis can be left with significant neurological aftereffects [5, 6] . Thus, exploring the mechanisms of the CVB5 infection might provide new insights into the pathology of HFMD and novel targets for the treatment of meningitis. Circular RNAs (circRNAs) represent an emerging type of RNA with single-stranded RNA and is expressed in almost all species. Unlike linear RNAs, circRNAs are characterized by back-splicing and the terminated 5 caps and 3 tails are not found within circRNAs [7] . MicroRNAs (miRNAs) are sponged by circRNAs as competitive endogenous RNAs (ceR-NAs), thus preventing miRNAs from binding to their target mRNAs and resulting in the increased expression of targets, such as circRNAs, which could suppress the production of viral proteins and promote viral clearance by binding miR-122 during a hepatitis C virus (HCV) infection [8] . This regulatory mechanism is the most common mechanism for circRNAs. In addition, RNA-sequencing is revealing an increasing number of circRNA expression profiles after virus infection; evidence suggests that some circRNAs can regulate host-virus interaction through an innate immune response. For example, Qu et al. found that circRNA-AIVR mediated the expression of the CREB Binding Protein (CREBBP) by binding to miRNA, which promoted the production of IFN-β and antagonized the replication of the influenza virus [9] . CircRNA regulated immune responses in viral infections have gradually been demonstrated, but without a focus on the functions and mechanisms involved in the regulation of pathogenesis for CVB5 infection.

In this study, we performed a circRNA expression analysis of CVB5-infected human rhabdomyosarcoma (RD) cells. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were adopted to predict potential physiological activities and related signal pathways. We also performed a circRNA-miRNA-mRNA interaction analysis. To better understand host-CVB5 interactions, we determined that novel_circ_0002006 and novel_circ_0001066 exerted effects on viral replication, and we investigated their effects on target miRNA expression; thus, we report a comprehensive catalog of circRNAs and provide a transcriptome to identify novel molecular targets and pathways for the treatment of HFMD.

To identify the circRNA expression profiles associated with CVB5 infection, we obtained sequencing data through the Illumina sequencing platform. As shown in Table S1 , a total of 5461 circRNAs were identified in all samples, and information relating to sequence characteristics was obtained. Of these, 3331 circRNAs were observed to be commonly expressed in both treatment groups; however, 493 and 1637 circRNAs were identified as being exclusively expressed in CVB5-infected cells and mock-infected cells, respectively ( Figure 1A ). The identified circRNAs were evenly distributed across different chromosomes, except for chromosome Y ( Figure 1B ). The circRNAs were mainly derived from exons, with a small proportion formed by introns and intergenic regions ( Figure 1C ). The circRNA sequences ranged from 26 nucleotides (nt) to 1323 nt in length, and the majority of circRNAs were 100 nt to 500 nt in length ( Figure 1D ). 

The circRNA expression levels in CVB5-infected and mock-infected RD cells were compared. There were 235 differentially expressed circRNAs (63 were upregulated and 172 were downregulated) for which the expression level changed more than twofold after CVB5 infection, relative to mock-infected cells (Figure 2A ). Inter-sample correlation analysis demonstrated that there was an obvious distinction in circRNA expression patterns between the two groups. Only two circRNAs were simultaneously detected in the cells for different treatments. Some changes to circRNA expression involved a change from absent to present ( Figure 2B ). 

The circRNA expression levels in CVB5-infected and mock-infected RD cells were compared. There were 235 differentially expressed circRNAs (63 were upregulated and 172 were downregulated) for which the expression level changed more than twofold after CVB5 infection, relative to mock-infected cells ( Figure 2A ). Inter-sample correlation analysis demonstrated that there was an obvious distinction in circRNA expression patterns between the two groups. Only two circRNAs were simultaneously detected in the cells for different treatments. Some changes to circRNA expression involved a change from absent to present ( Figure 2B ).

To reveal the biological functions of the differentially expressed circRNAs, GO assignments and KEGG pathway analysis were used to analyze the functions of the host genes of differentially expressed circRNAs. The most significantly enriched GO terms for a cellular component (CC) were 'intracellular' and 'intracellular part'. The biological process categories were mainly involved in the 'cellular macromolecule metabolic process' and 'macromolecule metabolic process' ( Figure 2C , Table S2 ). KEGG pathway analysis indicated that differentially expressed circRNAs were mainly enriched in ubiquitin mediated proteolysis, tight junctions, the ErbB signaling pathway, and the TNF signaling pathway ( Figure 2D , Table S3 ). To reveal the biological functions of the differentially expressed circRNAs, GO assignments and KEGG pathway analysis were used to analyze the functions of the host genes of differentially expressed circRNAs. The most significantly enriched GO terms for a cellular component (CC) were 'intracellular' and 'intracellular part'. The biological process categories were mainly involved in the 'cellular macromolecule metabolic process' and 'macromolecule metabolic process' ( Figure 2C , Table S2 ). KEGG pathway analysis indicated that differentially expressed circRNAs were mainly enriched in ubiquitin mediated proteolysis, tight junctions, the ErbB signaling pathway, and the TNF signaling pathway ( Figure 2D , Table S3 ).

To validate the differential expression of the RNA-seq results, we randomly selected six circRNAs (three upregulated and three downregulated circRNAs) for qPCR validation according to the expression heat map ( Figure 3A , Table S4 ). As shown in Figure 3B , these results indicated that the novel_circ_0010492, hsa_circ_0075346, novel_circ_0001066 were 

To validate the differential expression of the RNA-seq results, we randomly selected six circRNAs (three upregulated and three downregulated circRNAs) for qPCR validation according to the expression heat map ( Figure 3A , Table S4 ). As shown in Figure 3B , these results indicated that the novel_circ_0010492, hsa_circ_0075346, novel_circ_0001066 were upregulated, whereas hsa_circ_0008378, hsa_circ_0004390, and novel_circ_0002006 were downregulated. All results were consistent with the results arising from RNA-seq ( Figure 3C ). upregulated, whereas hsa_circ_0008378, hsa_circ_0004390, and novel_circ_0002006 were downregulated. All results were consistent with the results arising from RNA-seq ( Figure 3C ). 

To investigate the potential for circRNAs to act as ceRNAs in a CVB5 infection, we selected five enriched circRNAs (hsa_circ_0008378, hsa_circ_0031831, novel_circ_0002006, hsa_circ_0075346, and novel_circ_0001066) that were involved in the innate immune pathway based on GO and KEGG analyses and established two circRNA-miRNA-mRNA interaction networks. Four circRNAs (hsa_circ_0008378, hsa_circ_0031831, novel_circ_0002006, hsa_circ_0075346) in the IFN signaling pathway were found to target 42 miRNAs and 129 target genes ( Figure 4A , Table S5 ). Meanwhile, novel_circ_0001066 was found to be involved in the NF-κB signaling pathway-related network and was associated with 19 miRNAs and 78 target genes ( Figure 4B , Table S5 ). These data show that the circRNA-miRNA-mRNA axes may form a complex interaction network to play a role in CVB5 infection. 

To investigate the potential for circRNAs to act as ceRNAs in a CVB5 infection, we selected five enriched circRNAs (hsa_circ_0008378, hsa_circ_0031831, novel_circ_0002006, hsa_circ_0075346, and novel_circ_0001066) that were involved in the innate immune pathway based on GO and KEGG analyses and established two circRNA-miRNA-mRNA interaction networks. Four circRNAs (hsa_circ_0008378, hsa_circ_0031831, novel_circ_0002006, hsa_circ_0075346) in the IFN signaling pathway were found to target 42 miRNAs and 129 target genes ( Figure 4A , Table S5 ). Meanwhile, novel_circ_0001066 was found to be involved in the NF-κB signaling pathway-related network and was associated with 19 miRNAs and 78 target genes ( Figure 4B , Table S5 ). These data show that the circRNA-miRNA-mRNA axes may form a complex interaction network to play a role in CVB5 infection. 

We transfected pcDNA3.1_novel_circ_0002006 and pcDNA3.1_novel_circ_0001066 into RD cells, respectively. The results showed that novel_circ_0002006 significantly inhibited CVB5 replication, and upregulated PPRs (RIG-I and MDA5), IFN-β, and some critical interferon-stimulating factors (ISGs) (OASL, MXA, ISG15, IFIT1, IFIT2, and IFITM3) ( Figure 5A ,B). Moreover, novel_circ_0002006 also significantly inhibited CVB5 replication and upregulated inflammatory factors (TNF-α and IL-6) ( Figure 5C ,D). These results suggest that the two circRNAs exhibit antiviral effects through different innate immune pathways. 

We transfected pcDNA3.1_novel_circ_0002006 and pcDNA3.1_novel_circ_0001066 into RD cells, respectively. The results showed that novel_circ_0002006 significantly inhibited CVB5 replication, and upregulated PPRs (RIG-I and MDA5), IFN-β, and some critical interferon-stimulating factors (ISGs) (OASL, MXA, ISG15, IFIT1, IFIT2, and IFITM3) ( Figure 5A,B) . Moreover, novel_circ_0002006 also significantly inhibited CVB5 replication and upregulated inflammatory factors (TNF-α and IL-6) ( Figure 5C,D) . These results suggest that the two circRNAs exhibit antiviral effects through different innate immune pathways. experiments (B,D) . Ns, not significant, * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001, two-tailed unpaired Student's t test.

Using RT-qPCR to detect target miRNAs, according to novel_circ_0002006 and novel_ circ_0001066 networks, found that the expression levels of miR-4526 and miR-3123 were upregulated while those of miR-152-3p were inhibited in terms of novel_circ_0002006 overexpression ( Figure 6A ). The overexpression of novel_circ_0001066 inhibited miR-29b-3p and miR-567 but promoted miR-15b-3p ( Figure 6B ). These results indicated that complicated circRNA-miRNA-mRNA networks were active during CVB5 infection and suggested that these interactive networks regulated viral replication. 

CircRNAs have received widespread attention in the field of RNA due to their highly stable closed-loop structure and highly conserved nature [10] . HFMD caused by viral infection is one of the most common diseases in human infants and children [4] . Remarkably, circRNAs are closely associated with viral infections [11] [12] [13] . CircSIAE targets TAOK2 via miR-331-3p sponge adsorption and inhibits Coxsackie virus B3 proliferation [14] . In our current study, we systematically analyzed the expression profiles of circRNAs infected with CVB5 in RD cells. In addition, we constructed circRNA-miRNA-mRNA networks focused on innate immune response. Moreover, the effect of novel_circ_0002006 and novel_circ_0001066 on CVB5 replication were verified, and the expression of target proteins downstream of the enrichment pathway (IFN and NF-κB signaling pathway) were tested. Furthermore, we verified the possibility that circRNA acted through sponges.

We used RNA-seq to reveal the expression profiles of circRNAs during CVB5 infection. 5461 circRNAs were identified from different genomic locations. Our results further suggested that most circRNAs originate from exons (Table S4) , a finding that is consistent with the spectrum of circRNAs that have been found to be associated with CVA16 and EV71 infection [15, 16] ; therefore, we speculated that this may be related to the virus rather than cell specificity. Furthermore, the length distribution of circRNAs was consistent for the majority of circRNA expression profiles, except for some circRNAs with more than 1000 nucleotides, although this difference may depend on the specificity of the virus or host tissue. Compared with controls, we identified 235 circRNAs that were differentially expressed after CVB5 infection (Table S4 ). KEGG analysis of differentially expressed circRNA host genes identified 100 pathways, for which the p-value was < 0.05. We found that host genes corresponding to differentially expressed circRNAs were mainly involved in ubiquitin mediated proteolysis, the ErbB signaling pathway, aminoacyl-tRNA biosynthesis, and the TNF signaling pathway. These results suggested that circRNAs may be involved in the regulation of these cellular processes. This was unlike the circRNA enrich- 

CircRNAs have received widespread attention in the field of RNA due to their highly stable closed-loop structure and highly conserved nature [10] . HFMD caused by viral infection is one of the most common diseases in human infants and children [4] . Remarkably, circRNAs are closely associated with viral infections [11] [12] [13] . CircSIAE targets TAOK2 via miR-331-3p sponge adsorption and inhibits Coxsackie virus B3 proliferation [14] . In our current study, we systematically analyzed the expression profiles of circRNAs infected with CVB5 in RD cells. In addition, we constructed circRNA-miRNA-mRNA networks focused on innate immune response. Moreover, the effect of novel_circ_0002006 and novel_circ_0001066 on CVB5 replication were verified, and the expression of target proteins downstream of the enrichment pathway (IFN and NF-κB signaling pathway) were tested. Furthermore, we verified the possibility that circRNA acted through sponges.

We used RNA-seq to reveal the expression profiles of circRNAs during CVB5 infection. 5461 circRNAs were identified from different genomic locations. Our results further suggested that most circRNAs originate from exons (Table S4) , a finding that is consistent with the spectrum of circRNAs that have been found to be associated with CVA16 and EV71 infection [15, 16] ; therefore, we speculated that this may be related to the virus rather than cell specificity. Furthermore, the length distribution of circRNAs was consistent for the majority of circRNA expression profiles, except for some circRNAs with more than 1000 nucleotides, although this difference may depend on the specificity of the virus or host tissue. Compared with controls, we identified 235 circRNAs that were differentially expressed after CVB5 infection (Table S4 ). KEGG analysis of differentially expressed circRNA host genes identified 100 pathways, for which the p-value was < 0.05. We found that host genes corresponding to differentially expressed circRNAs were mainly involved in ubiquitin mediated proteolysis, the ErbB signaling pathway, aminoacyl-tRNA biosynthesis, and the TNF signaling pathway. These results suggested that circRNAs may be involved in the regulation of these cellular processes. This was unlike the circRNA enrichment pathway associated with EV71 and CVA16 infections in SH-SY5Y cells, such as the Wnt signaling pathway, angiogenesis, the p53 pathway, and the PDGF signaling pathway [15, 16] . This may be related to the cellular specificity of the host's circRNA. RD cells are often used to isolate the HFMD virus from clinical specimens, and are also highly susceptible to CVB5 infection. Since the expression of circRNAs underwent alterations after CVB5 infection, it appeared that circRNAs were associated with CVB5 specific replication. Furthermore, the enriched pathways were also demonstrated to be closely associated with the viral infection process. For instance, ubiquitin-mediated protein hydrolysis has been shown to be involved in a number of viral processes, including viral replication, maturation, and release [17, 18] . Our current findings may provide new avenues to explore the pathogenesis of the CVB5 infection.

CircRNAs can act as "miRNA sponges", thereby inhibiting the binding of miRNAs to target genes and thus regulating the expression of mRNAs in certain physiological and pathogenic activities; therefore, the elucidation of circRNA-miRNA-mRNA networks after CVB5 infection in host cells is important for understanding the pathogenesis of CVB5. The host innate immune system is the first line of defense against viral invasion; however, viruses can successfully survive via immune escape. The interaction of the IFN signaling pathway with viral infections has been extensively documented, for example, in the Hepatitis B virus (HBV), Influenza A viruses (IAV), and the Japanese encephalitis virus (JEV) [19] [20] [21] . Thus, we focused on circRNAs that were enriched in innate immune system processes, and constructed circRNA-associated circRNA-miRNA-mRNA networks by probing potential miRNA binding sites. Hsa_circ_0008378, hsa_circ_0031831, novel_circ_0002006, and hsa_circ_0075346 were involved in the IFN signaling pathwayrelated network in which 129 mRNAs competed for 42 miRNAs. The novel_circ_0001066 was involved in the NF-κB signaling pathway-related network, in which 78 mRNAs competed for 19 miRNAs. Of these, we focused on miR-152-3p, which plays a role in multiple biological processes. It has been shown that the downregulation of miR-152 can induce impaired hepatic gluconeogenesis by targeting PTEN, which is involved in miR-152-3pmediated hepatocyte gluconeogenesis via the regulation of the AKT/GSK pathway [22] . In addition, miR-152-3p was aberrantly expressed in hepatitis B virus-associated hepatocellular carcinoma [23] . Furthermore, miR-152-3p plays an important role in regulating the cell cycle, controlling tumor growth, and reducing apoptosis [23] [24] [25] .

In addition, there is accumulating evidence that inflammation is a major driving force in virus infection. For example, circ_ 0000220 and miR-326-3p have been shown to regulate the production of inflammatory cytokines during infection [26] . In the circRNA-miRNA-mRNA triple network diagram, novel_circ_0001066 was involved in the NF-κB signaling pathway-related network in which 78 mRNAs competed for 19 miRNAs. With regard to miR-29b-3p, IFN-I signaling-related miRNAs (ISR-miRNAs), have improved significantly during acute severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection [27] ; however, the role of miR-29b-3p in inflammation is yet to be investigated; therefore, the altered expression profiles of circRNAs in CVB5-infected cells may exert their biological functionality by interacting with specific microRNAs, for instance, miR-152-3p with miR-29b-3p; however, the expression of miR-152-3p and miR-29b-3p can also affect viral replication. There are currently neither characteristic nor functionality reports about novel_circ_0002006 and novel_circ_0001066.

By performing an in-depth study of circRNA function, it is possible to gain a better understanding of the mechanism of disease occurrence and improve the prevention and treatment of viral diseases. Based on these complex networks, we validated the role of circRNAs on CVB5 infection and the expression of related pathway target genes. We found that novel_circ_0002006 dramatically inhibited the expression of CVB5 VP1 via the activation of ISGs by using the IFN signaling pathway. Novel_circ_0001066 acted via the NF-κB signaling pathway to influence the inhibitory effect of inflammatory factors on CVB5 infection. These results suggested that novel_circ_0002006 and novel_circ_0001066 are representative of potential inhibitors of CVB5 infection. Furthermore, circRNAs have regulatory effects on NF-κB-related and IFN-related IRF3 with ISGs.

In summary, our data demonstrated expression profiles of circRNAs during CVB5 infection RD. Differentially expressed circRNAs may be involved in the mechanisms underlying innate immunity and could act as microRNA sponges to regulate the expression of mRNAs. These results may facilitate future studies investigating the molecular functions of circRNAs in viral pathogenesis.

Human striated muscle sarcoma cells (RD, American Type Culture Collection, CCL-136), the green monkey kidney cell line (Vero, American Type Culture Collection, CCL-81), and human embryonic kidney 293T (HEK-293T, American Type Culture Collection, CRL-3216) cells, were cultured in Dulbecco's modified Eagle medium (DMEM, HyClone, Logan, UT, USA), supplemented with 10% heat-inactivated fetal bovine serum (FBS) (AusGeneX, Brisbane, QLD, Australia), and incubated in CO 2 incubator at 37 • C.

The CVB5 strain (GeneBank: Mh201081.1) was isolated in the Kunming, Yunnan province in 2014, and preserved in our laboratory. CVB5 virus titer detection was performed in Vero cells. To investigate the effects of circRNAs on CVB5, RD (7 × 10 5 /well) was placed in a 6 well culture and infected with 1 multiplicity of infection (MOI) of CVB5 diluted in DMEM. After 2 h, the DMEM containing viral particles was exchanged with a 2% FBS medium for an additional 24 h. Control cells were treated equally but with the absence of the virus.

Total RNAs were isolated from infected, as well as uninfected RD cells, using RNAiso Plus (Takara, Tokyo, Japan), in accordance with the manufacturer's instructions. The purity, concentration, and integrity of extracted total RNAs were evaluated using a NanoPhotomete (IMPLEN, Westlake Village, CA, USA), Qubit ® 2.0 (Life Technologies, Carlsbad, CA, USA) and Agilent 2100 (Agilent Technologies, Santa Clara, CA, USA). A total amount of 3 µg of RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using a NEBNext ® UltraTM RNA Library Prep Kit for Illumina ® (NEB, Ipswich, MA, USA) following manufacturer's recommendations. In brief, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. First strand cDNA was then synthesized using random hexamer primers and M-MuLV Reverse Transcriptase (RNase H). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. Following the adenylation of 3 ends of DNA fragments, NEBNext Adaptor with a hairpin loop structure were ligated to prepare for hybridization. To select cDNA fragments that were preferentially 150-200 bp in length, the library fragments were purified with a AMPure XP system (Beckman Coulter, Beverly, MA, USA). Then 3 µL of USER Enzyme (NEB, Ipswich, MA, USA) was incubated with size-selected and adaptor-ligated cDNA at 37 • C for 15 min, followed by 5 min at 95 • C before PCR. PCR products were purified (AMPure XP system) and library quality was assessed using an Agilent Bioanalyzer 2100 system.

The clustering of index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v3-cBot-HS (Illumia, San Diego, CA, USA) in accordance with the manufacturer's instructions. Following cluster generation, the library preparations were sequenced on an Illumina Hiseq platform and 150 bp paired-end reads were generated.

Human genome sequences and gene annotations were downloaded from the NCBI database (https://www.ncbi.nlm.nih.gov/projects/genome/guide/human/index.shtml accessed on 12 June 2019). The paired-end clean reads were aligned to the reference genome using find_circ and CIRI2, and then junction reads at the same location were then identified as circRNA [28] .

Expressed circRNAs were quantified by TPM (transcript-per-million). Differential expression analysis of two groups (three biological replicates per condition) was performed using the DESeq2 R package (1.10.1) [29] . CircRNAs with p-values < 0.05 and fold changes (FC) > 2 were defined as differentially expressed Gene Ontology (GO) enrichment analysis of differentially expressed circRNAs and statistical enrichment in KEGG pathways were implemented by the cluster Profiler R package 2110. CircRNAs with p-values < 0.05 were considered to be significantly enriched [30] .

Total RNAs were extracted with RNAiso, plus (Takara, Tokyo, Japan) mock-infected and CVB5-infected RD cells. Reverse transcription was carried out using Hifair ® II 1st Strand cDNA Synthesis SuperMix for qPCR (Yeasen Biotechnology Co., Ltd., Shanghai, China) in accordance with the manufacturer's protocol. Thereafter, qPCR was performed using Hieff ® qPCR SYBR Green Master Mix (Yeasen Biotechnology Co., Ltd., Shanghai, China) on a 7500 real-time PCR system programmed with the following cycling conditions: 95 • C for 5 min, followed by 40 cycles of 95 • C for 5 s, extension and annealing at 60 • C for 34 s. The 2 -∆∆Ct method was used for quantification. The expression levels of circRNAs were normalized to those of ACTB. Primer sequences were given in Table S6A . The expression levels of miRNAs were normalized to those of U6. Primer sequences were given in Table S6B .

Putative miRNA binding sites in six circRNAs that were enriched in the innate immune pathway, based on GO and KEGG analyses, were identified by employing the miRanda algorithm. CircRNA-miRNA-mRNA interaction networks were constructed using Cytoscape (v3.4.0) [31] .

The cDNA of circRNA0001066 and circRNA0002006 were cloned into the pcDNA3.1 vector to yield pcDNA3.1_circRNA0001066 and pcDNA3.1_ circRNA0002006. Transfection was carried out using the TransIntroTM EL Transfection Reagent (TRANSgene, Beijing, China) in accordance with the manufacturer's instructions.

Total protein was extracted by radioimmunoprecipitation assay (RIPA) lysis buffer (Solarbio, Beijing, China), and the protein concentration for each sample was calculated by a bicinchoninic acid (BCA) protein assay kit. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) for protein separation and polyvinylidene difluoride filter (PVDF) membranes were used to transfer separated proteins. PVDF membranes with protein bands were blocked with 5% (w/v) nonfat dry milk, and probed with anti-α-Tubulin antibody and anti-CVB5 VP1 antibody. Then, the membranes were incubated with diluted secondary antibodies for 1 h at room temperature. Eventually, the membranes were scanned by chemiluminescence image analysis system Tanon 5200 (Tanon, Shanghai, China).

Data are presented as means ± SD. CircRNA expressions were compared between CVB5 infected-cells and controls using Student's t-test. The expressions of mRNA and miRNA were compared between control cells and circRNA overexpression cells using Student's t-test as well. p-values < 0.05 were considered to be statistically significant.

In this study, we identified, for the first time, the expression profiles of circRNAs that are associated the interaction of CVB5 and host cells, and constructed circRNA-associated regulatory networks relating to a host's first line of defense against viral infection-innate immunity. We further investigated the role of circRNA on the replication of CVB5 and the possible mechanisms underlying its antiviral effect. These findings provided new insights for our understanding of virus-host interactions, and revealed that circRNAs may play an important role in CVB5 infection. 

Picornavirus and enterovirus diversity with associated human diseases

Ongoing change of severe hand, foot, and mouth disease pathogens in Yunnan

Molecular epidemiology of enterovirus from children with herpangina or hand, foot, and mouth disease in Hangzhou

Emerging epidemiology and the need for a vaccine strategy

Genome Sequences of Coxsackievirus B5 Isolates from Two Children with Meningitis in Australia

An outbreak of aseptic meningitis caused by a distinct lineage of coxsackievirus B5 in China

Circular RNAs: A novel type of non-coding RNA and their potential implications in antiviral immunity

Functional sequestration of microRNA-122 from Hepatitis C Virus by circular RNA sponges

A Novel Intronic Circular RNA Antagonizes Influenza Virus by Absorbing a microRNA That Degrades CREBBP and Accelerating IFN-β Production

The emerging role and significance of circular RNAs in viral infections and antiviral immune responses: Possible implication as theranostic agents

Identification of virus-encoded circular RNA

Circular RNA Expression Profiles and the Pro-tumorigenic Function of CircRNA_10156 in Hepatitis B Virus-Related Liver Cancer

The role of circular RNAs in viral infection and related diseases

The circRNA circSIAE Inhibits Replication of Coxsackie Virus B3 by Targeting miR-331-3p and Thousand and One Amino-Acid Kinase Front

circRNA expression patterns and circRNA-miRNA-mRNA networks during CV-A16 infection of SH-SY5Y cells

Comprehensive analysis of the circRNA expression profile and circRNA-miRNA-mRNA network in the pathogenesis of EV-A71 infection

E3 ligase Nedd4l promotes antiviral innate immunity by catalyzing K29-linked cysteine ubiquitination of TRAF

Regulation of the linear ubiquitination of STAT1 controls antiviral interferon signaling

HBcrAg Levels Are Associated With Virological Response to Treatment With Interferon in Patients With Hepatitis Delta

Porcine TRIM35 positively regulate TRAF3-mediated IFN-β production and inhibit Japanese encephalitis virus replication

MicroRNA 152 regulates hepatic glycogenesis by targeting PTEN

Identification of miRNAs associated with dendritic cell dysfunction during acute and chronic hepatitis B virus infection

Long non-coding RNA LINC00174 promotes glycolysis and tumor progression by regulating miR-152-3p/SLC2A1 axis in glioma

Mutant NPM1-regulated lncRNA HOTAIRM1 promotes leukemia cell autophagy and proliferation by targeting EGR1 and ULKJ

Genome-wide profiling of host-encoded circular RNAs highlights their potential role during the Japanese encephalitis virus-induced neuroinflammatory response

Expression of plasma IFN signalingrelated miRNAs during acute SARS-CoV-2 infection and its association with RBD-IgG antibody response

Circular RNA identification based on multiple seed matching

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

clusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters

Cytoscape: A software environment for integrated models of Biomolecular Interaction Networks

Institutional Review Board Statement: Not applicable.

Data Availability Statement: The genetic data presented in this paper are publicly available in the GenBank database under accession No. GEO180816.

All authors declare no conflict of interest.