key: cord-0682143-n3dayryc
authors: Barbagallo, Davide; Palermo, Concetta Ilenia; Barbagallo, Cristina; Battaglia, Rosalia; Caponnetto, Angela; Spina, Vittoria; Ragusa, Marco; Di Pietro, Cinzia; Scalia, Guido; Purrello, Michele
title: Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells
date: 2022-01-17
journal: Cell Mol Life Sci
DOI: 10.1007/s00018-021-04119-8
sha: 0c65321aadddc6fcaf58c5bd4f987f6acc482639
doc_id: 682143
cord_uid: n3dayryc

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

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is an enveloped virus classified as a new member of the Coronaviridae family, Betacoronavirus genus, whose genome consists of a single stranded positive ( +) RNA molecule about 30 kilobases long [1, 2] . SARS-CoV-2 is etiologically involved in the life-threatening coronavirus disease 2019 and is responsible for a pandemic outbreak in March, 2020 [3] . From December 2019 to September 2021, 3,508 SARS-CoV-2 genomes were sampled and 20 clades were identified around the world, according to the GISAID database (https:// nexts train. org/ ncov/ gisaid/ global). Emergence of new SARS-CoV-2 genotypes alerts the scientific community to the possibility that some variants of concern may bypass the immune barrier given by the vaccines currently used and highlights the importance of the identification of therapeutic targets helpful for the management of this clinically very relevant disease [4] .

Circular RNAs (circRNAs) are a recently discovered class of RNAs, mainly synthesized through backsplicing and characterized by the covalent bond between their 5' and 3' termini [5, 6] . CircRNAs follow tissue-and developmental-specific expression patterns and are mainly localized in the cell cytoplasm [7, 8] . The best characterized functions of circRNAs consist in sponging microRNAs (miRNAs) and RNA-binding proteins (RBPs) [5, 9] : in the first case, circRNAs may be typically involved in competitive endogenous RNA (ceRNA) networks [10] [11] [12] [13] ; in the second case, circRNAs may regulate biological processes within eukaryotic cells, such as assembly of preinitiation complex (PIC) at the beginning of transcription or also splicing [14, 15] . Moreover, circRNAs may function either as a template for the synthesis of generally short peptides, thanks to the presence of Internal Ribosomes Entry Sites (IRESs) within their sequences, or as a scaffold for the regulation of host gene transcription [16] [17] [18] . CircRNAs have been found aberrantly expressed in many cancers and degenerative diseases [19, 20] and are associated with several biological processes, both in physiological and pathological conditions [21] [22] [23] [24] . Due to their intrinsic resistance to the activity of exoribonucleases and their presence in several human body fluids as well as within extracellular vesicles, circRNAs have been suggested as good candidate diagnostic and prognostic biomarkers for several diseases [25] [26] [27] [28] .

Recent evidence has shown the etiological involvement of circRNAs in viral infections. Most specifically, cross-talk between host cell circRNA biogenesis and RBPs linked to immune response (e.g.: immune factors NF90/NF110) has been described [29] . Influenza virus-infected A549 cells showed the induction of a circRNA that acts as a sponge for miRNAs regulating the expression of interferon beta (IFN-β) enhanceosome [30] . At the same time, circRNAs have been demonstrated to be synthesized from the genome of several DNA viruses (e.g.: Herpesviruses), contributing to the infection's pathogenesis [31] [32] [33] .

Recently, thanks to RNA-seq data analysis from cells infected with Middle East Respiratory Syndrome Coronavirus (MERS-CoV), Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and SARS-CoV-2 RNA (+) betacoronaviruses, several circRNAs of viral origin have been identified and characterized [34, 35] . At the same time, analysis of human circRNAome revealed several differentially expressed (DE) circRNAs in lung epithelial cells as well as in the peripheral blood of individuals infected with SARS-CoV-2 [36, 37] . Based on gene expression datasets, perturbation of ceRNA networks (circRNAs-miRNAs-mRNAs) in host cells following SARS-CoV-2 infection has also been predicted [38] . To improve our knowledge of the molecular dynamics of SARS-CoV-2 infection and prospectively identify new candidate therapeutic targets, in this study we focused on circRNAs synthesized from the viral genome, suggesting their involvement in COVID-19 pathogenesis.

Three ml of universal transport medium (UTM™) (COPAN Italia SpA, Brescia, Italy) from nasopharyngeal swabs of individuals suspected to be infected by SARS-CoV-2 were used for diagnostic purposes. A residual 1 ml of UTM™ was centrifuged at 350×g for 5 min at 4 °C to pellet cell debris. Supernatants were discarded and cell pellets stored at − 80 °C until further processing. Nucleic acids were isolated directly from UTM™ through a ThermoFisher Flex apparatus by the Mag-MAXTM Viral Pathogen Kit (ThermoFisher Scientific, Monza, Italy) for diagnostic purposes. Diagnosis was performed by a multiplex real-time PCR, Allplex™ SARS-CoV-2 Master assay (Seegene Inc., Arrow Diagnostics, Genoa, Italy), following the instructions of the manufacturer. The method amplifies SARS-CoV-2 E, N, RdRp, and S genes, according to World Health Organization (WHO)'s guidelines. Amplification cycles were performed on a Bio-Rad CFX96 real-time PCR instrument (Bio-Rad, Segrate-Milan, Italy). Only samples showing cycle threshold (Ct) values ≤ 35 for all the transcripts assayed were considered positive. A total of 15 positive and 6 and negative samples were assayed in this retrospective study. Data on biological specimens anonymously collected in this study were processed in accordance with the ethical principles reported in the Declaration of Helsinki.

RNA was extracted from cells previously collected by nasopharyngeal swab using TRIzol ® (ThermoFisher Scientific, ThermoFisher Scientific, Waltham, MA, USA), according to the manufacturer's instruction [39] . RNA was quantified by a GenQuant pro spectrophotometer (Biochrom, Cambridge, UK). CircRNAs and mRNAs were amplified using a Power SYBR ® Green RNA-to-CT™ 1-Step Kit (ThermoFisher Scientific). MiRNAs were reverse transcribed into cDNA by TaqMan™ MicroRNA Reverse Transcription Kit (Ther-moFisher Scientific) and amplified through TaqMan™ Universal Master Mix II (ThermoFisher Scientific). PCRs were run on a 7900HT real-time PCR instrument (ThermoFisher Scientific). PCR products of the amplified circ_3205 were then sent to BMR Genomics, Padua, Italy (www. bmr-genom ics. it) for purification through ExoSap (Applied Biosystems™, ThermoFischer Scientific). Sanger sequencing was performed with BigDyeTM Terminator v3.1 Cycle Sequencing Kit (ThermoFisher Scientific) on a ABI 3730xl DNA Analyzer (Applied Biosystems™, ThermoFisher Scientific). Sequences and IDs of specific primer pairs and TaqMan probes used in this manuscript are listed in Supplemental Table 1 .

SARS-CoV-2 candidate circRNA sequences were retrieved from VirusCircBase (v. 1.1) [40] , based on RNA-Seq experiments performed on Calu3 cells infected with the virus for 12 h and 24 h (GSE148729). Candidate circRNAs were filtered through the following criteria: (i) backsplice junction reads counted by at least two circRNA predictive algorithms among CIRI2; circRNA_finder; find_circ [6, 41, 42] , for the same circRNA; (ii) abundancy of circRNA (at least two backsplice junction reads for each predictive algorithm).

Interactions between candidate circRNA and miRNAs were predicted by TarpMiR [43] , Analysis of Common Targets for circular RNAs (ACT) [44] , and STarMir [45] . FASTA sequences of human miRNAs (from miRBase 22 release [46] ) and candidate circRNA were given as input to each of the three predictive algorithms. TarpMiR was set choosing human model and a probability cutoff of 0.5. Human V-CLIP data were used to train STarMir predictions [47] . Only miRNAs predicted to interact with SARS-CoV-2 candidate circRNA by all the three tools were considered as potentially implied in the ceRNA network.

MiRNAs predicted to interact with candidate circRNA were given as input to DIANA miRPath v3.0 [48] and Biological Process (BP) GO's subcategory was analyzed. MicroT-CDS (MicroT and False Discovery Rate (FDR)-corrected p value thresholds = 0.8 and 0.05, respectively) was selected as the algorithm for the prediction of miRNA-mRNA interactions. BP GO analysis filtered miRNA targets based on: (i) their potential involvement in blood clotting and immune response pathways, known to be related to COVID-19; (ii) dysregulated expression in lung or nasopharyngeal cells from SARS-CoV-2 infected individuals.

First and second neighbor interactants of the candidate miRNA targets were retrieved from the HUman Reference protein Interactome (HuRI) database [49] . The list of interactions was given as input to Cytoscape (version 3.8.2) [50] and the network generated was analyzed through g:GOSt, within the g:Profiler web interface [51] , and the Cytoscape plugin cytoHubba [52] to study gene functional enrichment and topological features, respectively. Topological analysis focused on the centrality parameters: betweenness; bottleneck; closeness; eccentricity; radiality and stress; for each of them, a corresponding subnetwork has been generated and analyzed for GO BP enrichment. 

Real-time PCR data were analyzed through the 2 −DDCt method [53] and two-sided Student's t test was applied to identify DE transcripts. Spearman's correlation test was used to identify positive or negative correlations among transcripts. Modified Fisher's exact test followed by false discovery rate methodology was used to calculate p values for GO analysis. p values ≤ 0.05 were considered statistically significant.

Based on VirusCircBase, a total of 3473 circRNAs were detected by at least one tool among CIRI2; circRNA_finder; find_circ (Supplemental Table 2 ). SARS-CoV-2 circRNAs 363, 368, 2667, 2670, 2685, 2795, 3058, 3205 were the top eight most expressed circRNAs in Calu-3 infected with SARS-CoV-2 at least in one time point, according at least to two out of three algorithms among CIRI2; circRNA_finder; find_circ (Table 1) .

The expression of SARS-CoV-2 circRNAs 363, 368, 2667, 2670, 2685, 2795, 3058, and 3205 was assayed in a discovery cohort made of a pool of three positive samples and three UCs. Based on real-time PCR and gel electrophoresis data, we focused on circ_3205, a 286 nt-long circRNA whose sequence is embedded within the open reading frame (ORF) coding for the nucleocapside (N) protein of SARS-CoV-2 (nucleotide position 28,609-28,898, Wuhan-Hu-1 reference genome, NCBI Reference Sequence: NC_045512.2), clearly present and highly expressed only in positive samples ( Fig. 1A and  1B) . Sanger sequencing of the PCR amplification product of circ_3205 confirmed the presence of the 3'-5' junction, specific of this circRNA (Fig. 1C) . Gene expression assay in the validation cohort confirmed that circ_3205 was expressed only in positive samples and its amount positively correlated with that of Spike (S) mRNA and SARS-CoV-2 viral load (r values = 0.80952 and 0.84867, p values (two sided) = 0.015 and 0.016, respectively, Spearman's correlation test) (Fig. 2) . The expression of S mRNA positively correlated with SARS-CoV-2 viral load, too (r value = 0.92582, p values (two sided) = 0.003) (Fig. 2 ).

Hsa miRNAs 298, 3940-5p, 4640-5p, 6081 and 6133 were predicted to interact with circ_3205 by all the three predictive algorithms queried (TarpMiR, ACT and STarMir) (Fig. 3A) . BP GO analysis highlighted a potential involvement of the mRNA targets regulated by hsa miRNAs 298, 3940-5p, 4640-5p, 6081 and 6133 in SARS-CoV-2 infection-related processes (blood coagulation and immune response) (Fig. 3B ).

Based on (i) literature data, (ii) our GO analysis, and (iii) the high probability of interaction with SARS-CoV-2_circ_3205, we focused on hsa-mir-298 (Fig. 4A) .

Has-miR-298 was predicted to target 18 and 30 mRNAs Table 3 ). Ten predicted mRNA targets (CPB2; DCTN1; FN1; KCNMB4; MERTK; MYO1C; PDE1A; PIK3C3; PRKCE; SERPINB2) were selected as candidates because of their known upregulation in biological specimens from SARS-CoV-2-infected individuals, as compared with UCs (Supplemental Table 4 ). 

PPI network generated by HuRI consisted of 482 nodes and 1452 edges. GO analysis of the whole network revealed an over-representation of biological functions linked to blood coagulation, immune response, and inflammation (Supplemental Fig.1 ). The analysis of centralities of the network revealed a total of 27 most central proteins: among them, EGFR, HSP90AB1, YWHAZ occurred in five out of six subnetworks made of the most central nodes (Fig. 5) . The generated subnetworks revealed an enrichment in BPs related to SARS-CoV-2 infection program (Fig. 6 ).

The capability of viruses to synthesize circRNAs upon infection has been ascertained, especially in DNA viruses [54] [55] [56] [57] [58] [59] [60] . As an example, Epstein Barr virus (EBV), a double stranded DNA virus belonging to the Herpesviridae family, is known to produce about 30 circR-NAs during different phases of its infection [61] and some of them (e.g.: circBART2.2) contribute to virusinduced carcinogenesis through the immune escape of Fig. 3 A Venn diagrams showing the number of miRNAs predicted to interact with circ_3205 by TarpMir, ACT and STarMiR predictive tools. The total number of miRNAs predicted by each tool is reported in brackets, while the number of miRNAs, either specifically predicted by only one or by more than one predictive tool, is reported (without brackets) within the Venn diagrams. The name of the five candidate miRNAs predicted to interact with circ_3205 by all the three predictive tools is reported in extenso. B Horizontal bar chart representing the enrichment of biological processes carried out by the predicted targets of hsa miRNAs 298, 3940-5p, 4640-5p, 6081, and 6133. Biological processes are reported in y-axis; statistical significance is reported as -LOG (p value) (x-axis) ◂ Fig. 4 A Graphical representation of the interaction predicted by STarMiR between hsa-miR-298 and circ_3205. MiRNA seed region is showed in red. The numbers below vertical bars indicate the nucleotide position of circ_3205. B Heatmap showing the expression of seven predicted hsa-miR-298 targets whose expression was revealed in positive (P) and negative (N) samples by real-time PCR. Expression is reported as DCt: the lower is this value, the higher is the expression of the target and viceversa. * = statistically significant DE transcript (p value < 0.05, Student's t test) nasopharyngeal carcinoma cells [62] . Kaposi's sarcomaassociated herpesvirus (KSHV) generates circRNAs found to be inserted into virions and implied in several steps of the infection [63] .

In this study, we first assayed the expression of eight circRNAs predicted to be synthesized by SARS-CoV-2 upon host cell infection. Based on our experimental results, we then focused on circRNA 3205. Based on data stored in VirusCircBase, circ_3205 was predicted to be synthesized from the negative RNA strand of SARS-CoV-2, specifically from a sequence embedded in the ORF coding for the N protein of the virus. Although the mechanism of circRNA biogenesis from RNA (+) viruses is under investigation, some hypotheses may be proposed: in a recent study, it was found that SARS-CoV-2 genome may be in part reverse transcribed and integrated as DNA into the host genome through a LINE1-mediated mechanism, leading to the production of chimeric viral-host cellular transcripts [64, 65] . Based on this study, it is conceivable that viral circRNAs could be generated through backsplicing from these chimeric viral-host transcripts. Nevertheless, the hypothesis of an integration of SARS-CoV-2 genome (or a part) into the host genome is debated [66] . An alternative path of cir-cRNA biogenesis from RNA (+) viral genomes may consist in splicing-independent mechanisms occurring in the cytoplasm of host cells: this mechanism was previously described for IRE1alpha-mediated XBP1 mRNA splicing in mammalian cells, for IL1b transcripts in platelets and for the recently suggested miR-7-mediated circularization of the CDR1AS transcript [67] [68] [69] . Expression of viral circRNAs may perturb the ceRNA networks originally present within host cells or may create new ones [70] . Specifically, our data suggest that once synthesized within the host cell, circ_3205 may function as sponge for hsa-miR-298, allowing for the upregulation of targets involved in the progression through the infection (Fig. 7) . MiR-298, together with miR-296, belongs to a genomic locus that is imprinted both in mice and humans [71, 72] . Interestingly, hsa-miR-298 has been predicted to bind the 5'-UTR of the SARS-CoV-2 genome, potentially altering its secondary structure and negatively impacting on its capability to be translated after infection of the host cell [73] . Chopra N. et al. identified hsa-miR-298 as a potential therapeutic agent for Alzheimer's disease, because of its capability to negatively regulate the expression of human amyloid-β precursor protein (APP), β-site APP-converting enzyme 1 (BACE1) and specific tau protein isoforms [74] . Hsa-miR-298 has also been defined as oncomiRNA in several cancers, thanks to its ability to downregulate proapoptotic proteins such as BAX and PTEN [75, 76] . These literature data support our hypothesis that the circ_3205 sponge effect against hsa-miR-298 may contribute to the progression of the infection, by stabilizing the SARS-CoV-2 genome and triggering biological processes such as inflammation and apoptosis. Our data also shed light on two predicted targets of hsa-miR-298 (KCNMB4 and PRKCE), which we found to be upregulated in positive samples and whose expression positively correlated with that of circ_3205. KCNMB4 encodes a β4 subunit of a voltagedependent K + channel, belonging to the Ca2 + -activated Slo subfamily (BK) [77, 78] . Upregulation of KCNMB4 correlates with the increased intracellular concentration of Ca 2+ observed during SARS-CoV-2 infection [79] . Even though the function exerted by BK channels during the SARS-CoV-2 infection needs further investigation, the role of K + concentration and K + channels in facilitating the entry of some viruses into the host cells has been convincingly ascertained [80, 81] . Furthermore, abnormalities in electrolyte serum concentrations (especially sodium, potassium, calcium and chloride) have been found to be related with the prognosis of COVID-19 patients and with the possibility to develop blood clots [82] [83] [84] . PRKCE encodes a Ca 2+ -independent protein kinase, belonging to the subfamily of nonconventional protein kinase C (PKCs); it has been described as involved in SARS-CoV infection, through the calcium-independent PI3K/PKCε/JNK/CREB pathway; this, in turn, induces COX-2 expression upon the interaction between viral S protein and cellular receptors [85] . COX-2 has been found upregulated in several cell types also after SARS-CoV-2 infection [86] . PRKCE expression is further induced by Interferon-α (IFN-α), one of the first cytokines synthesized in infected cells on innate immune response [87] . Due to an abnormal recruitment of proinflammatory cells, IFN-α signaling over a prolonged period of time is known to cause an uncontrolled inflammatory response and potential organ failure in tissues infected by SARS-CoV-2 as well as other respiratory viruses [88, 89] . The study of the PPI network, generated starting from KCNMB4 and PRKCE, revealed Heat Shock Protein 90, Alpha family, class B member 1 (HSP90AB1) as one of the most central nodes, according to five out of six parameters of network centrality. HSP90AB1 is a first neighbor interactant of PRKCE and belongs to the Heat Shock Protein 90 (HSP90) family. Some members of HSP90 family have been recently suggested to foster MERS-CoV, SARS-CoV and SARS-CoV-2 replication and proinflammatory Fig. 5 Analysis of centralities of KCNMB4 and PRKCE's PPI network. Subnetworks show the top ten most central nodes calculated by A betweenness, B bottleneck, C closeness, D eccentricity, E radiality, F stress methods. For each subnetwork the scalebar reporting the scores of centrality is shown ◂ cytokine expression [90, 91] . Subnetworks generated by the study of centralities further revealed an enrichment in BPs strictly related to the local and systemic effects of SARS-CoV-2 infection, such as remodeling of protein trafficking within infected host cell [92] , mucus hypersecretion [93] , ErbB protein family and growth factor receptor signaling [94, 95] , and unfolded protein response [96] . Collectively, these findings corroborate our hypothesis of a functional involvement of KCNMB4 and PRKCE in SARS-CoV-2 infection.

Based on the integration of our experimental data and predictive analysis, we propose SARS-Cov-2_circ_3205/ hsa-miR-298/KCNMB4 and SARS-Cov-2_circ_3205/ hsa-miR-298/PRKCE molecular axes as involved in the progression of SARS-Cov-2 infection and, more in detail, in the related processes of blood clotting and immune response, respectively (Fig. 7) . 

A Novel coronavirus from patients with pneumonia in China

A pneumonia outbreak associated with a new coronavirus of probable bat origin

A new coronavirus associated with human respiratory disease in China

The novel coronavirus originating in Wuhan, China: challenges for global health governance

Natural RNA circles function as efficient microRNA sponges

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

Cell-type specific features of circular RNA expression

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

The GAU GAA Motif Is Responsible for the Binding between circS-MARCA5 and SRSF1 and Related Downstream Effects on Glioblastoma Multiforme Cell Migration and Angiogenic Potential

A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?

The multilayered complexity of ceRNA crosstalk and competition

CircRNA_000203 enhances the expression of fibrosis-associated genes by derepressing targets of miR-26b-5p, Col1a2 and CTGF, in cardiac fibroblasts

LncRNA UCA1, upregulated in CRC biopsies and downregulated in serum exosomes, controls mRNA expression by RNA-RNA interactions

The RNA binding protein quaking regulates formation of circRNAs

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

A novel protein encoded by the circular form of the SHPRH gene suppresses glioma tumorigenesis

Circular RNA circ0005276 promotes the proliferation and migration of prostate cancer cells by interacting with FUS to transcriptionally activate XIAP

Therapeutic targeting of circ-CUX1/ EWSR1/MAZ axis inhibits glycolysis and neuroblastoma progression

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

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations

circMine: a comprehensive database to integrate, analyze and visualize human disease-related circRNA transcriptome

CircNAPEPLD is expressed in human and murine spermatozoa and physically interacts with oocyte miRNAs

CircAtlas: an integrated resource of one million highly accurate circular RNAs from 1070 vertebrate transcriptomes

Circbank: a comprehensive database for circRNA with standard nomenclature

Serum extracellular vesicle-derived circHIPK3 and circS-MARCA5 are two novel diagnostic biomarkers for glioblastoma multiforme

Plasmatic exosome-derived circRNAs panel act as fingerprint for glioblastoma

Serum hsa_circ_101178 as a Potential Biomarker for Early Prediction of Osteoarthritis

Plasma-derived exosomal circular RNA hsa_ circ_0005540 as a novel diagnostic biomarker for coronary artery disease

Coordinated circRNA biogenesis and function with NF90/NF110 in viral infection

A novel intronic circular RNA antagonizes influenza virus by absorbing a microRNA that degrades CREBBP and accelerating IFN-β production

Identification of virus-encoded circular RNA

Characterizing expression and regulation of gamma-herpesviral circular RNAs

Role of virally encoded circular RNAs in the pathogenicity of human oncogenic viruses

Identification and characterization of circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2

Circular RNA profiling reveals abundant and diverse circRNAs of SARS-CoV-2, SARS-CoV and MERS-CoV origin

Differential host circRNA expression profiles in human lung epithelial cells infected with SARS-CoV-2

A study of differential circRNA and lncRNA expressions in COVID-19-infected peripheral blood

Circular RNA-Micro-RNA-MRNA interaction predictions in SARS-CoV-2 infection

MicroRNAs are stored in human MII oocyte and their expression profile changes in reproductive aging

VirusCircBase: a database of virus circular RNAs

Circular RNA identification based on multiple seed matching

Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation

TarPmiR: a new approach for micro-RNA target site prediction

Analysis of common targets for circular RNAs

STarMir Tools for Prediction of microRNA Binding Sites

miRBase: annotating high confidence microRNAs using deep sequencing data

A quantitative analysis of CLIP methods for identifying binding sites of RNA-binding proteins

DIANA-miRPath v3.0: deciphering microRNA function with experimental support

A reference map of the human binary protein interactome

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update)

cytoHubba: identifying hub objects and sub-networks from complex interactome

Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method

Viruses join the circular RNA world

Circular DNA tumor viruses make circular RNAs

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

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

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

RNA-Seq Revealed a circular RNA-microRNA-mRNA regulatory network in hantaan virus infection

Characterizing viral circRNAs and their application in identifying circRNAs in viruses

The Epstein Barr virus circRNAome

Epstein-barr virus-encoded circular RNA CircBART2.2 promotes immune escape of nasopharyngeal carcinoma by regulating PD-L1

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

Reverse-transcribed SARS-CoV-2 RNA can integrate into the genome of cultured human cells and can be expressed in patient-derived tissues

Exogenous coronavirus interacts with endogenous retrotransposon in human cells

No evidence of SARS-CoV-2 reverse transcription and integration as the origin of chimeric transcripts in patient tissues

Escaping the nuclear confines: signal-dependent pre-mRNA splicing in anucleate platelets

A new molecular mechanism of RNA circularization and the microRNA sponge formation

HumanViCe: host ceRNA network in virus infected cells in human

MicroRNAs 296 and 298 are imprinted and part of the GNAS/ Gnas cluster and miR-296 targets IKBKE and Tmed9

miR-296-3p, miR-298-5p and their downstream networks are causally involved in the higher resistance of mammalian pancreatic α cells to cytokine-induced apoptosis as compared to β cells

Computational analysis of targeting SARS-CoV-2, viral entry proteins ACE2 and TMPRSS2, and interferon genes by host microRNAs

Micro-RNA-298 reduces levels of human amyloid-β precursor protein (APP), β-site APP-converting enzyme 1 (BACE1) and specific tau protein moieties

Bufalin promotes apoptosis of gastric cancer by down-regulation of miR-298 targeting bax

2020) miR-298 plays a pivotal role in colon cancer invasiveness by targeting PTEN

Cloning and functional characterization of novel large conductance calcium-activated potassium channel beta subunits, hKC-NMB3 and hKCNMB4

Mechanism of beta4 subunit modulation of BK channels

SARS-CoV-2, hypoxia, and calcium signaling: the consequences and therapeutic options

Potassium is a trigger for conformational change in the fusion spike of an enveloped RNA virus

Baseline electrolyte abnormalities would be related to poor prognosis in hospitalized coronavirus disease 2019 patients

Renal histopathological analysis of 26 postmortem findings of patients with COVID-19 in China

Epidemiology, prognosis and management of potassium disorders in Covid-19

Spike protein of SARS-CoV stimulates cyclooxygenase-2 expression via both calcium-dependent and calcium-independent protein kinase C pathways

Cyclooxgenase-2 is induced by SARS-CoV-2 infection but does not affect viral entry or replication

Interferon alpha induces protein kinase C-epsilon (PKC-epsilon) gene expression and a 4.7-kb PKC-epsilon-related transcript

Mouse model of SARS-CoV-2 reveals inflammatory role of type I interferon signaling

Human coronavirus dependency on host heat shock protein 90 reveals an antiviral target

Transcriptomic profiling of SARS-CoV-2 infected human cell lines identifies HSP90 as target for COVID-19 therapy

SARS-CoV-2 disrupts splicing, translation, and protein trafficking to suppress host defenses

Cytokine storm and mucus hypersecretion in COVID-19: review of mechanisms

Discovery of pan-ErbB inhibitors protecting from SARS-CoV-2 replication, inflammation, and lung injury by a drug repurposing screen

Growth factor receptor signaling inhibition prevents SARS-CoV-2 replication

Modulation of the unfolded protein response by the severe acute respiratory syndrome coronavirus spike protein

Acknowledgements We thank Dr Elisabetta Tendi for her kind support in the first phases of this work. We also thank the Scientific Bureau of the University of Catania for language support. Dr Cristina Barbagallo is actually receiving a protempore contribution (assegno di ricerca)

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

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00018-021-04119-8.Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells