- Open Access
CircEAF2 counteracts Epstein-Barr virus-positive diffuse large B-cell lymphoma progression via miR-BART19-3p/APC/β-catenin axis
Molecular Cancer volume 20, Article number: 153 (2021)
Epstein-Barr virus (EBV) represents an important pathogenic factor of lymphoma and is significantly associated with poor clinical outcome of diffuse large B-cell lymphoma (DLBCL). Circular RNAs (circRNAs) play an essential role in lymphoma progression. However, the underlying mechanism of circRNA on DLBCL progression related to EBV remains largely unknown.
CircRNA was screened by high-throughput sequencing in tumor samples of 12 patients with DLBCL according to EBV infection status. Expression of circEAF2, as well as the relationship with clinical characteristics and prognosis, were further analyzed in tumor samples of 100 DLBCL patients using quantitative real-time PCR. Gain- and loss-of-function experiments were conducted to investigate the biological functions of circEAF2 both in vitro and in vivo. The underlying mechanism of circRNA on DLBCL progression were further determined by RNA sequencing, RNA pull down assay, dual-luciferase reporter assay, rescue experiments and western blotting.
We identified a novel circRNA circEAF2, which was downregulated in EBV + DLBCL and negatively correlated with EBV infection and DLBCL progression. In EBV-positive B lymphoma cells, circEAF2 overexpression induced lymphoma cell apoptosis and sensitized lymphoma cells to epirubicin. As mechanism of action, circEAF2 specifically targeted EBV-encoded miR-BART19-3p, upregulated APC, and suppressed downstream β-catenin expression, resulting in inactivation of Wnt signaling pathway and inhibition of EBV + DLBCL cell proliferation. In EBV-positive B-lymphoma murine models, xenografted tumors with circEAF2 overexpression presented decreased Ki-67 positivity, increased cell apoptosis and retarded tumor growth.
CircEAF2 counteracted EBV + DLBCL progression via miR-BART19-3p/APC/β-catenin axis, referring circEAF2 as a potential prognostic biomarker. Therapeutic targeting EBV-encoded miRNA may be a promising strategy in treating EBV-associated lymphoid malignancies.
Epstein-Barr virus (EBV) is an important human oncogenic virus and significantly related to pathogenesis and progression of diffuse large B-cell lymphoma (DLBCL) . Of note, EBV is the first virus reported to encode microRNAs (miRNAs) , which are derived from two viral genome segments, BamHI fragment H rightward open reading frame 1 (BHRF1) and Bam HI-A region rightward transcript (BART) . Recent study has shown that these EBV-encoded miRNAs contribute to B-cell transformation and lymphomagenesis . Functionally, miR-BHRF1-2 represses the transcriptional activity of PRDM1, a master regulator of B-cell terminal differentiation, to promote EBV lymphomagenesis . MiR-BARTs also play a crucial role in B-cell lymphoma progression. MiR-BART3, miR-BART9, and miR-BART17-5p act as the post-transcriptional regulators of BCL6 expression in DLBCL . Given the abovementioned roles of EBV-encoded miRNAs, targeting EBV as a therapeutic approach appears to be promising in DLBCL.
Circular RNAs (circRNAs) are formed as closed loops by precursor RNA splicing. Owing to their circular structure to avoid degradation by RNases, circRNAs are highly stable in cells . With the development of high-throughput sequencing and in-depth mechanistic researches, circRNAs have been increasingly recognized as modulators of pathological processes in lymphomas . For instance, circLAMP1 enhances T-cell lymphoblastic lymphoma cell proliferation and inhibits cell apoptosis by regulating miR-615-5p/DDR2 axis , while circAPC sponges miR-888, initiates APC translation, and counteracts DLBCL progression . CircRNAs may thus function as oncogenes or tumor suppressors, respectively, by directly sponging miRNAs or modulating their target genes. However, it is still unknown how circRNA interacts with EBV-encoded miRNAs for the pathogenesis and progression of DLBCL.
In the present study, we performed circRNA high-throughput sequencing in DLBCL according to EBV infection status and identified that a novel circRNA circEAF2 was negatively correlated with EBV infection and DLBCL progression. Meanwhile, the biological functions of circEAF2 on tumor progression were evaluated both in vitro and in vivo. Mechanistic analysis showed that circEAF2 specifically targeted EBV-encoded miR-BART19-3p, upregulated APC, suppressed downstream β-catenin expression, and counteracted EBV + DLBCL progression.
Patients and methods
A total of 100 patients with newly diagnosed DLBCL were enrolled in this study. The histological diagnosis was established according to World Health Organization (WHO) classification. All patients were treated by R-CHOP (rituximab combined with cyclophosphamide, doxorubicin, vincristine and prednisone)-based immunochemotherapy. The study was approved by the Shanghai Rui Jin Hospital Review Board with informed consent obtained from all subjects in accordance with the Declaration of Helsinki.
Human B-lymphoma cell lines (Farage, SU-DHL-4, OCI-LY-10, OCI-LY-7, Toledo), and EBV-producing marmoset B-cell line B95-8 were purchased from ATCC (American Type Culture Collection. Manassas, VA, USA). Cells were cultured in RPMI-1640 (Gibco, Shanghai, China) and humidified atmosphere of 95% air and 5% CO2 at 37 °C. HEK-293 T was cultured in DMEM (Gibco, Shanghai, China). EBV-infected human B-lymphoma cell models were established using the method previously described . Briefly, EBV virion was prepared from EBV-producing B95-8 cells. SU-DHL-4 and OCI-LY-10 cells were exposed to the supernatant of B95-8 cells at 37 °C for acute EBV infection (within 4 days) or long-term EBV infection (for 3 weeks). SU-DHL-4 and OCI-LY-10 cells cultured in EBV-free growth medium were used as negative controls.
EBV DNA quantification and EBER in-situ hybridization
DNA was extracted from frozen tumor samples and cell lines using the QIAamp DNA Mini Kit (Qiagen, Valencia, CA, USA). Quantification of EBV-specific sequences was performed by real-time quantitative PCR with 7500HT Fast Real-time PCR system (Applied Biosystems, Foster City, CA, USA) using EBV PCR Fluorescence Quantitative Diagnostic Kit (DaAn Gene Co, Sun Yat-sen University, Guangzhou, China). The copy number of EBV DNA in each sample was calculated from a standard curve with a cut-off value of 700 copies/ml in tumor samples . The primer sequences were obtained from DaAn Gene Co. and covered by patent. EBER (EBV-encoded RNA) in-situ hybridization was performed on 3 μm paraffin sections with 20% cut-off value of EBER-positive cells .
CircRNA-Seq high-throughput sequencing
CircRNA-Seq high-throughput sequencing and subsequent bioinformatics analysis were performed by Cloud-Seq Biotech (Shanghai, China). Total RNA was used to prepare the CircRNA sequencing library by circRNA Enrichment Kit (Cloud-seq Biotech, Shanghai, China). RNA libraries were constructed by using pretreated RNAs with TruSeq Stranded Total RNA Library Prep Kit (Illumina, San Diego, CA, USA) and quality controlled by the BioAnalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). Paired-end reads were harvested from Illumina HiSeq 4000 sequencer and quality controlled by Q30. The reads were aligned to the reference genome/transcriptome with STAR software and circRNAs were detected with DCC software and annotated by circBase database .
RNA extraction, reverse transcription, and quantitative real-time PCR
Total RNA was isolated from tumor samples and cell lines using TRIzol reagent (Invitrogen, Life Technologies, Carlsbad, CA, USA) and genomic DNA (gDNA) were isolated using FastPure DNA Isolation (Vazyme, Nanjing, China). For circRNA and mRNA, reverse transcriptions were performed using the PrimeScript RT Master Mix (Takara, Shiga, Japan) with random primers. For miRNA, reverse transcriptions were performed using the PrimeScript RT Reagent Kit (Takara) with specific stem-loop primers and cDNA amplification performed using TB Green Premix Ex Taq II (Takara). GAPDH and U6 were used as internal control and each sample was repeated in triplicate. The relative fold-change in expression with respect to a control sample was calculated by the 2-ΔΔCt method. Divergent primers were used to detect backsplice junction of circRNA and convergent primers were used to detect linear mRNA. The primers for miR-BART19-3p and U6 small nuclear RNA were obtained from RiboBio Company (Guangzhou, China). The primers were listed in Table S1.
CircRNA plasmid construction and transfection
To construct circEAF2-overexpressing plasmids, circEAF2 cDNA, as well as a specific flanking sequence of upstream and downstream (Listed in Fig. S1A), were synthesized and cloned into the pCDH-CMV-EF1а-GFP-T2A-puro vector (Lingke Biotech, Shanghai, China). Plasmids were transfected into HEK-293 T cells to package lentivirus using Lipofectamine 3000 reagent (Invitrogen). Farage and EBV-infected SU-DHL-4 cells were infected with the packaged lentivirus and then selected with 2 μg/ml puromycin (Sigma, St. Louis, MO, USA) for 14 days until circEAF2 was stably overexpressed. Circularization of the circEAF2 lentiviral construct was further validated in Fig. S1. MiRNA mimics, miRNA mutants, as well as siRNA targeting circEAF2 and control siRNA, were purchased from GenePharma (Shanghai, China) and transfected with Lipofectamine 3000 reagent (Invitrogen). The oligos were shown in Table S2.
RNase R and actinomycin D treatment
RNase R treatment was used to identify and confirm the circRNA. Total RNA (2 μg) was incubated for 30 min at 37 °C with 3 U/μg of RNase R (Epicenter, Madison, WI, USA). Transcription was prevented by the addition of 2 μg/ml actinomycin D (Glpbio, Montclair, CA, USA) for indicated time. After treatment with actinomycin D and RNase R, RNA expression levels of circEAF2 and EAF2 mRNA were detected by qRT-PCR.
Nucleic acid electrophoresis
The qRT-PCR products of cDNA and gDNA were investigated using 2% agarose gel electrophoresis with TAE running buffer. DNA was separated by electrophoresis at 100 V for 30 min. The DNA marker was DL2000 DNA Marker (2000-100 bp) (Takara). The bands were examined by UV irradiation.
Luciferase report assay
The sequences of miR-binding site of circEAF2 or 3′-UTR of APC, as well as their corresponding mutants, were designed, synthesized, and inserted into luciferase reporter vector pmirGLO (Lingke Biotech). HEK-293 T cells were seeded into 24-well plates and transfected with luciferase reporter vector and miRNA mimics using the Lipofectamine 3000 reagent. After 48 h incubation, the Firefly and Renilla luciferase activities were quantified with a dual-luciferase reporter assay (Promega, Madison, WI, USA). The groups transfected wild-type circEAF2 luciferase reporter without miRNA mimics and scramble mimics were used as control group (Ctrl) and negative control group (NC), respectively.
RNA pull-down assay
The pull-down assay was performed using Pierce Magnetic RNA-Protein Pull-Down Kit as previously described . In brief, Farage cells that stably expressed circ-EAF2 were harvested, lysed, and sonicated. To pull down miRNA by circRNA, biotinylated circ-EAF2 probe was synthesized by Cloud-Seq Biotech, and oligo probe was referred as a control. After 2 h incubation with magnetic beads, cell lysates were incubated with the probes overnight. After incubation, the bound RNAs were washed and purified for further analysis. To pull down circRNA by miRNA, biotinylated miRNA mimics or their mutants were synthesized by Cloud-Seq Biotech and incubated with magnetic beads. Cell lysates were incubated with the probes overnight, and bound RNAs were detected for the abundance of circEAF2 by qRT-PCR. The sequences of probes were shown in Table S3.
Cell proliferation and in vitro chemosensitivity assay
B lymphoma cells were seeded into 96-well plates at an initial density of 1 × 104 cells per well and incubated with the indicated concentrations of epirubicin at 37 °C. After incubation, cell counting Kit-8 (CCK-8; Dojindo, Japan) was added and the absorbance of each well was measured at 450 nm by spectrophotometry. Cell proliferation curves were plotted using fold change of the absorbance at each time point and the 50% inhibitory concentration (IC50) values of epirubicin were calculated using Graphpad.
Cell apoptosis assay
Apoptosis assay was performed using Annexin V-PE Apoptosis Detection Kit (Invitrogen) and TUNEL Cell Apoptosis Detection Kit (Beyotime, Shanghai, China) according to the manufacturer’s protocol. In brief, either transfected Farage or EBV-infected SU-DHL-4 cells (4 × 105 per well) were treated with binding buffer with 5 μl Annexin V-PE and analyzed for apoptosis rates by flow cytometry (BD Biosciences, San Jose, CA, USA). Apoptotic cells were stained with 50 μL of TUNEL for 30 min at room temperature and captured under a fluorescence microscope (Leica, Wetzlar, Germany).
RNA sequencing was downloaded from the National Omics Data Encyclopedia (NODE, http://www.biosino.org/node) under accession number OEP0001143 . The raw RNA sequencing reads were aligned to human reference genome hg19 using Hisat2 (v2.0.5) and STAR (v2.5.2b). Transcript counts table files were generated by the HTSeq. Package ‘limma’ (v3.38.3) was used to normalize the raw reads and obtain differentially expressed genes (DEGs). Normalized reads were analyzed by the package ‘clusterProfiler’ (v3.10.1) and gene set enrichment analysis (GSEA) was performed. Metascape, a gene-list analysis tool, was used to verify the biological process enrichment of DEGs associated with EBV infection and low circEAF2 expression . CatRAPID, an algorithm to estimate the binding propensity of protein-RNA pairs, was used to predict the binding site between targeted proteins and parental EAF2 transcript .
Cytoplasmic and nuclear proteins were extracted from cells using BeyoECL Plus Nuclear and Cytoplasmic Protein Extraction Kit (Beyotime) supplemented with protease inhibitors and phosphatase inhibitors. Western blotting was performed as previously described . Primary antibodies specific to GAPDH (1:1000 dilution, Cell Signaling Technology, Danvers, Massachusetts, USA), H3 (1:1000 dilution, Cell Signaling Technology), and APC (1:1000 dilution, Abcam, Cambridge, UK) were used. The blots were then incubated with goat anti-rabbit or anti-mouse secondary antibody (1:3000 dilution, Abcam) and visualized by Immobilon Western Chemiluminescent HRP Substrate (Millipore, St. Louis, MO, USA).
In vivo tumorigenesis assay
Female NOD-SCID mice (4 weeks old) were randomly divided into two groups, with ten mice in each group. Farage cells stably transfected with circEAF2 overexpression plasmids or control vectors (3 × 106) were injected subcutaneously into the right flank of mice randomly. Treatments were started after the tumor volumes reached about 0.5 × 0.5 cm. The untreated group received an equal volume of PBS, whereas the epirubicin group received epirubicin (1.50 mg/kg) twice weekly for 2 weeks [19, 20]. Tumor volumes were calculated as 0.5 × a (length) × b (width)2.
After 2 weeks, the tumors were dissected and fixed in 10% formalin. Tumor tissues were cut into 4 μm sections. The slides were treated for antigen retrieval and incubated with the Ki-67 antibody (1:100 dilution, Dako, Glostrup, Denmark) overnight at 4 °C, following by HRP-conjugated secondary antibody (1:3000 dilution, Dako) for 2 h at room temperature.
Difference of circEAF2 expression among groups was calculated using Mann-Whitney U test. Survival curves were plotted using Kaplan-Meier method and compared by log-rank test. The Cox proportional hazards regression model was constructed for univariate and multivariate analysis. In vitro experimental results were expressed as mean ± S.D. of data obtained from three separate experiments and determined by t-test to compare variance. Statistical significance was defined as P < 0.05.
Identification and characterization of circEAF2 in EBV + DLBCL
We first performed circRNA high-throughput sequencing in frozen tissue samples of 12 DLBCL patients, including 8 EBV + DLBCL and 4 EBV-DLBCL, as previously defined . In total, 107 circRNAs were found significantly downregulated in EBV + DLBCL (P < 0.05). To identify potentially functional circRNAs, more stringent filtering criteria (FDR < 0.1, P < 0.01, Fold change<− 2) were used and 8 candidate circRNAs were obtained (Fig. 1A). Among these circRNAs (Figs. S2 and S3), circEAF2 expression was the most significantly and consistently decreased in EBV-positive B lymphoma cell lines (P = 0.0001) and EBV + DLBCL patients (P = 0.0015, Fig. 1B). Particularly, when OCI-LY-10 and SU-DHL-4 cells were acutely infected with EBV, circEAF2 expression was significantly downregulated, as compared to negative control cells. During long-term EBV infection, circEAF2 expression remained low expression level, suggesting that circEAF2 downregulation was significantly associated with EBV infection in DLBCL (Fig. 1C).
CircEAF2 is produced by reverse splicing from exons 2, 3, 4, and 5 of the EAF2 gene locus (CircBase ID: hsa_circ_0066971, referred to as circEAF2, http://www.circbase.org/). To further determine the characteristics of circEAF2, we validated the putative circEAF2 junction using circRNA-specific divergent primers and Sanger sequencing (Fig. 1D). CircEAF2 was more resistant to RNase R than the liner form of EAF2 (mRNA EAF2). And the half-life of circEAF2 was longer than that of mRNA EAF2 (Fig. 1E). Moreover, the convergent primers could amplify mRNA EAF2 using both cDNA and gDNA as templates, while the divergent primers could only amplify circEAF2 using cDNA templates, but not gDNA templates (Fig. 1F). These results indicated that circEAF2 was in a circular form and functionally stable in cells.
CircEAF2 expression is associated with EBV infection and tumor progression in DLBCL
To determine whether circEAF2 was involved in EBV infection of DLBCL, the expression of circEAF2 was further assessed in tumor samples of 100 DLBCL patients. According to EBV infection status (Fig. 2A), the relative expression of circEAF2 was significantly reduced in patients with tumor EBER positivity (P = 0.0003), tumor EBV DNA positivity (P < 0.0001), and serum EBV DNA positivity (P = 0.0017). Moreover, circEAF2 expression showed negatively correlation with EBV DNA copy number in tumor and in serum (Fig. 2B). In addition, tumor EBV DNA was significantly associated with serum EBV DNA, poor progression-free survival (PFS) and overall survival (OS) of DLBCL patients (Fig. S4).
To evaluate the clinical relevance of circEAF2 expression in DLBCL, the median value of relative circEAF2 expression was defined as the cutoff value for dividing DLBCL patients into high- and low-expression groups (Table 1). Indeed, circEAF2 expression was negatively correlated with poor performance status (ECOG≥2, P = 0.0095), elevated serum LDH (P = 0.0152), high-risk international prognostic index (IPI 3-5, P = 0.0029) and low rate of complete remission (CR, P = 0.0206). To evaluate the prognostic value of circEAF2 expression in DLBCL, the effects of circEAF2 expression on PFS and OS of the patients were analyzed. Kaplan-Meier curves indicated that patients with low circEAF2 expression had a remarkably shorter 2-year PFS (39.2% vs 63.3%, P < 0.0001) and 2-year OS (52.9% vs 65.3%, P = 0.0110) than those with high circEAF2 expression (Fig. 2C). Moreover, univariate analysis showed that age ≥ 60 yr, ECOG≥2, advanced disease stage, multiple extranodal involvements, elevated serum LDH, and IPI score > 2 were significant predictors of inferior PFS and OS, while high circEAF2 expression was significant predictive on favorable PFS and OS (Fig. 2D). In multivariate analysis, high circEAF2 expression remained an independent favorable predictor of PFS (P = 0.0031, Fig. 2D). Together, circEAF2 indicated EBV infection but could counteract DLBCL progression.
CircEAF2 counteracts tumor progression of EBV + DLBCL in vitro
To examine the biological functions of circEAF2 on EBV + DLBCL, EBV-positive B lymphoma cells (Farage cells) were stably transfected with lentiviruses using circEAF2 plasmid and showed significantly higher expression of circEAF2 than those transfected with vector (P < 0.0001). EBV-negative B lymphoma cells (SU-DHL-4 cells) were transfected with circEAF2 siRNA and showed significantly lower expression of circEAF2 than those transfected with control siRNA (P < 0.0001, Fig. 3A). After 48-h, ectopic circEAF2 expression resulted in a significant decrease in Farage cell proliferation, as compared to vector-transfected cells (Fig. 3B). Flow cytometry (Fig. 3C) and TUNEL assay (Fig. 3D) revealed that circEAF2 overexpression induced more apoptosis in Farage cells. Moreover, increased expression of circEAF2 sensitized Farage cells to epirubicin, with a significant decrease in the IC50 value, as compared to vector-transfected cells (42.4 ± 3.0 nM vs 88.2 ± 16.7 nM, P = 0.0095; Fig. 3E). In contrary, molecular silencing of circEAF2 had no significant effect on SU-DHL-4 cell proliferation, apoptosis and sensitivity to epirubicin (Fig. 3B-E). These data showed that circEAF2 could exert a protective role in EBV + DLBCL cells through promoting cell apoptosis and enhancing chemosensitivity.
CircEAF2 functions as a sponge of miR-BART-19 in EBV + DLBCL
As circRNAs mostly function as miRNA sponges to regulate downstream genes, we next performed biotinylated-circEAF2 probe pull-down assay and miRNA screening to identify circEAF2-interacting miRNAs. As shown in Fig. 4A, the specific enrichment of human miRNAs (miR-4286, miR-660, miR-301a, miR-579, miR-1307, miR-212, miR-942, miR-93, miR-550, miR-3173) and EBV-encoded miRNAs (miR-BART19-3p, miR-BART4, miR-BART7) were detected in the biotinylated-circEAF2 probe, as compared to the oligo probe in Farage cells. Among these miRNAs, only EBV-encoded miR-BART19-3p mimics could significantly reduce the circEAF2 luciferase reporter activity of wild-type circEAF2 in a dose-dependent manner, but not mutant circEAF2 (Fig. 4B and C). The sponge effect of circEAF2 on miR-BART19-3p was further confirmed by pulling down biotinylated-miR-BART19-3p (Fig. 4D). CircEAF2 was abundantly pulled down by biotinylated-miR-BART19-3p probe rather than oligo probe in Farage cells.
To further investigate the functional interaction between circEAF2 and miR-BART19-3p, rescue experiments were conducted by co-transfecting miR-BART19-3p mimics (miR-BART19-WT) with circEAF2 plasmid (circEAF2-WT) in Farage cells (P = 0.0155) or circEAF2 siRNA (si-circEAF2) in SU-DHL-4 cells (P < 0.0001; Fig. 4E). As expected, miR-BART19-3p mimics (miR-BART19-WT) could partially attenuate inhibition of Farage cell proliferation (Fig. 4F), induction of apoptosis (Fig. 4G) and chemosensitivity to epirubicin (Fig. 4H) mediated by circEAF2 overexpression. In contrary, only after circEAF2 silencing in SU-DHL-4 cells, miR-BART19-WT significantly promoted SU-DHL-4 cell proliferation (Fig. 4F), inhibited apoptosis (Fig. 4G) and chemosensitivity to epirubicin (Fig. 4H). These data demonstrated that circEAF2 could function as a sponge for EBV-encoded miR-BART19-3p to modulate EBV + DLBCL progression.
CirEAF2 regulates the Wnt signaling pathway through miR-BART19-3p/ APC/β-catenin axis
To investigate the underlying mechanism of circEAF2 on EBV + DLBCL, we analyzed RNA sequencing data to identify DEGs between low circEAF2 expression group and high circEAF2 expression group. GSEA revealed Wnt signaling pathway as the most significantly altered pathway (P = 0.0020, NES = 1.85, Fig. 5A). Previous reports have shown that miR-BART19-3p promotes EBV-negative cell proliferation and regulates the expression of APC, a key regulator of the Wnt signaling pathway [21,22,23]. We assessed the expression of circEAF2, miR-BART19-3p and APC in tumor samples of DLBCL by qRT-PCR (Fig. 5B) and found significant negative correlation between circEAF2 and miR-BART19-3p (P = 0.0493), positive correlation between circEAF2 and APC (P < 0.0001), and negative correlation between miR-BART19-3p and APC (P = 0.0443). Therefore, we proposed that APC could be a target of circEAF2/miR-BART19-3p regulatory network in EBV + DLBCL. The putative binding sites of APC with miR-BART19-3p were further tested in a dual-luciferase assay. The results showed that APC-3’UTR luciferase reporter activity was significantly inhibited by miR-BART19-3p mimics in a dose-dependent manner, but not by circEAF2. However, circEAF2 rescued the inhibition of miR-BATR19 on APC-3’UTR luciferase reporter activity in a dose-dependent manner (Fig. 5C). In Farage cells, APC-3’UTR luciferase reporter activity was significantly decreased by miR-BART19 mimics and partially rescued by co-transfection with circEAF2 (Fig. 5D). These findings confirmed that circEAF2 sponges miR-BART-19 to regulate APC expression.
As confirmed in qRT-PCR (Fig. 5E) and western blotting (Fig. 5F), the expression of APC was upregulated in circEAF2-overexpressed Farage cells at both mRNA and protein levels, whereas was partially rescued by co-transfection with miR-BART19-3p. Accordingly, nuclear expression of β-catenin was decreased after circEAF2 overexpression and partially restored by miR-BART19-3p mimics (Fig. 5F). Therefore, the mechanistic analysis showed that circEAF2 modulated miR-BART19-3p/APC/β-catenin axis to counteract EBV + DLBCL progression (Fig. 5G).
CircEAF2 counteracts tumor progression of EBV + DLBCL in vivo
In order to search for in vivo evidence of circEAF2 on EBV + DLBCL progression, murine xenograft models were established with subcutaneous injection of Farage cells stably transfected with circEAF2 or vector as control. Consistent with in vitro data, reduced volume of xenograft tumors was observed in the circEAF2-overexpressing group, as compared to that of the vector group. Moreover, the circEAF2-overexpressing group treated with epirubicin showed greater reductions in tumor growth (Fig. 6A). Accordingly, immunohistochemistry study and TUNEL assay showed that the circEAF2-overexpressing tumors, treated with or without epirubicin, exhibited decreased positivity of Ki-67 and increased rate of apoptotic cells, as compared to the vector group (Fig. 6B).
CircEAF2 may be regulated by EBV through RNA splicing
To elucidate the underlying mechanism of circEAF2 downregulation in EBV + DLBCL, RNA sequencing data were further analyzed for DEGs, according to EBV infection status and relative circEAF2 expression (Fig. S5A). In both EBV + DLBCL patients and low circEAF2 expression group, upregulated genes were identified and significantly enriched in mRNA metabolic process and RNA splicing (Fig. S5B). Considering circRNAs are generated by the backsplicing of precursor mRNA (pre-mRNA) transcripts, catRAPID was used to estimate the binding propensity of proteins with pre-mRNA EAF2. Among the upregulated genes, UPF1, RBM15B and DHX38 were predicted to have intense interaction with intronic complementary elements of pre-mRNA EAF2, suggesting the disruption of circEAF2 backsplicing formation (Fig. S5C). Moreover, relative EAF2 mRNA expression by RNA sequencing showed no difference, according to EBV infection status and relative circEAF2 expression in DLBCL (Fig. S5D). Neither upon acute nor upon long-term EBV infection, the EAF2 expression had no obvious change in SU-DHL-4 cells and OCI-LY-10 cells (Fig. S5E), suggesting that EBV infection may reduce expression of circEAF2 by RNA splicing, but not expression of linear EAF2. Together, these results supported the notion that EBV infection in DLBCL may alter RNA splicing and result in disruption of circEAF2 formation.
CircRNAs have important functions in oncogenesis and tumor progression [24, 25]. In DLBCL, circAPC suppresses tumor cell proliferation by sponging miR-888, upregulates APC expression, and eventually inactivates Wnt/β-catenin pathway. Clinically, high circAPC expression was associated with favorable prognosis . More recent study has reported that EBV-encoded circRNAs are implicated in EBV biology and cancers . However, the function and regulatory mechanism of circRNAs in EBV + DLBCL remain largely unknown.
In the present study, circEAF2 was identified as a novel circRNA significantly downregulated in EBV + DLBCL, negatively correlated with EBV infection, and indicated good clinical outcome of the patients. In vitro and in vivo studies further revealed that circEAF2 counteracted EBV + DLBCL progression through inducing lymphoma cell apoptosis and sensitizing lymphoma cell to epirubicin. Recently, several studies implicated the role of circRNAs in tumor sensitivity to chemotherapy . Circ-LARP4 is downregulated and correlated with prolonged survival in osteosarcoma, due to its ability to increase tumor cell sensitivity to doxorubicin by sponging miR-424 . Similar observations were reported in breast cancer, where circKDM4C suppresses tumor progression and enhances doxorubicin chemosensitivity by regulating miR-548p/PBLD axis . Epirubicin is a major anthracycline cytotoxic agent for DLBCL . EBV + DLBCL tends to have poor response and inferior survival upon treatment with standard anthracycline-based chemotherapy . Here we demonstrated that circEAF2 overexpression was associated with favorable prognosis of the patients and good response to anthracycline-based immunochemotherapy in DLBCL. Functionally, circEAF2 overexpression sensitized EBV-positive lymphoma cells to epirubicin both in vitro and in vivo. Together, circEAF2 exhibits a tumor suppressive role and could be considered as a favorable prognostic biomarker of EBV + DLBCL.
EBV-encoded miRNAs are critically involved in EBV-associated oncogenesis and function as oncogenic drivers , modulating host mRNAs involved in transcription regulation, apoptosis, and cell cycle control . This is, to our knowledge, the first study to explore the potential interaction between circRNA and EBV-encoded miRNA in DLBCL. Indeed, we screened oncogenic miRNAs and confirmed that circEAF2 was able to bind to EBV-encoded miR-BART19-3p. High expression of miR-BART19-3p was previously found in Hodgkin lymphoma  and NK/T-cell lymphoma . MiR-BART19-3p is implicated in the activation of Wnt signaling pathway by targeting tumor suppressor APC, enhancing proliferation and suppressing apoptosis of EBV-infected tumor cells . APC can downregulate Wnt signaling pathway by reducing the accumulation of β-catenin in the nucleus, thereby inhibit tumor growth and chemoresistance in multiple cancer types [34, 35]. Here in EBV-positive B-lymphoma cells, circEAF2 specifically targeted miR-BART19-3p, upregulated APC, and suppressed downstream β-catenin, indicative a viral miRNA-mediated mechanism underlying EBV + DLBCL progression. In EBV-associated nasopharyngeal carcinoma, circRPMS1 mediates oncogenesis through sponging multiple miRNAs . Our results elucidated not only the suppressive role of circEAF2 in EBV + DLBCL progression, but also the interaction between circEAF2 and miR-BART19-3p to become a potential therapeutic target for EBV + DLBCL treatment.
Previous research indicated that circRNA formation is regulated by RNA splicing regulatory proteins, such as eIF4A3 and ADAR1 [37, 38]. In this study, we found that RNA splicing-related genes were significantly enriched in EBV + DLBCL patients with low circEAF2 expression group. Among them, UPF1, RBM15B and DHX38 were predicted to bind to intronic complementary elements of pre-mRNA EAF2. UPF1, a nonsense-mediated RNA decay factor, is critically involved in host cellular response to EBV infection [24, 25]. RBM15B, initially identified as a binding partner of the EBV mRNA export factor, mediates RNA N6-methyladenosine methylation (m6A) and regulates RNA alternative splicing . DHX38, a pre-mRNA splicing-related DEAH box RNA helicase, interacts with satellite I noncoding RNA . Based on these observations, one possibility is that EBV infection in DLBCL may alter RNA splicing and result in disruption of circEAF2 formation. Upon both acute and long-term EBV infection, circEAF2 expression was significantly downregulated, while linear EAF2 expression had no obvious change. These findings supported that EBV infection actively represses circEAF2 formation by RNA splicing. Another possibility is that miR-BART-19-3p indirectly represses circEAF2 through negative selection as circEAF2 low-expressing cells outcompete circEAF2 high-expressing cells. As previously reported, EBV-encoded BART miRNAs conferred a selective growth advantage to EBV-positive tumor cells in vivo . Thus it is also likely that such oncogenic properties reported for BART-miRNAs could indirectly represses circEAF2 through negative selection. However, future study is needed to further illustrate the regulatory mechanisms on circEAF2 in EBV + DLBCL.
Understanding the precise function of circEAF2 in the progression of EBV + DLBCL may increase our knowledge of EBV lymphomagenesis and provide clinical rationale of developing new therapeutic strategy in treating EBV-associated lymphoid malignancies.
Availability of data and materials
Further information and requests for reagents and resources should be directed to and will be made available by the lead contact Wei-li Zhao (firstname.lastname@example.org) upon reasonable request.
Diffuse large B-cell lymphoma
BamHI fragment H rightward open reading frame 1
Bam HI-A region rightward transcript
World Health Organization
American type culture collection
Rituximab combined with cyclophosphamide, doxorubicin, vincristine and prednisone
Complement deoxyrubonucleic acid
International prognostic index
Eastern Cooperative Oncology Group
Germinal center B-cell
Polymerase chain reaction
Quantitative real-time PCR
Cell counting Kit-8
Terminal dexynucleotidyl transferase (TdT)-mediated dUTP nick end labeling
Half maximal inhibitory concentration
Differentially expressed genes
Gene set enrichment analysis
Kyoto Encyclopedia of Genes and Genomes
Beltran BE, Castro D, Paredes S, Miranda RN, Castillo JJ. EBV-positive diffuse large B-cell lymphoma, not otherwise specified: 2020 update on diagnosis, risk-stratification and management. Am J Hematol. 2020;95:435–45.
Pfeffer S, Zavolan M, Grässer FA, Chien M, Russo JJ, Ju J, et al. Identification of virus-encoded microRNAs. Science. 2004;304:734–6.
Grundhoff A, Sullivan CS, Ganem D. A combined computational and microarray-based approach identifies novel microRNAs encoded by human gamma-herpesviruses. RNA. 2006;12:733–50.
Caetano BFR, Jorge BAS, Müller-Coan BG, Elgui de Oliveira D. Epstein-Barr virus microRNAs in the pathogenesis of human cancers. Cancer Lett. 2021;499:14–23.
Ma J, Nie K, Redmond D, Liu Y, Elemento O, Knowles DM, et al. EBV-miR-BHRF1-2 targets PRDM1/Blimp1: potential role in EBV lymphomagenesis. Leukemia. 2016;30:594–604.
Martín-Pérez D, Vargiu P, Montes-Moreno S, León EA, Rodríguez-Pinilla SM, Lisio LD, et al. Epstein-Barr virus microRNAs repress BCL6 expression in diffuse large B-cell lymphoma. Leukemia. 2012;26:180–3.
Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20:675–91.
Perez de Acha O, Rossi M, Gorospe M. Circular RNAs in blood malignancies. Front Mol Biosci. 2020;7:109.
Deng L, Liu G, Zheng C, Zhang L, Kang Y, Yang F. Circ-LAMP1 promotes T-cell lymphoblastic lymphoma progression via acting as a ceRNA for miR-615-5p to regulate DDR2 expression. Gene. 2019;701:146–51.
Hu Y, Zhao Y, Shi C, Ren P, Wei B, Guo Y, et al. A circular RNA from APC inhibits the proliferation of diffuse large B-cell lymphoma by inactivating Wnt/β-catenin signaling via interacting with TET1 and miR-888. Aging (Albany NY). 2019;11:8068–84.
Zheng Z, Sun R, Zhao HJ, Fu D, Zhong HJ, Weng XQ, et al. MiR155 sensitized B-lymphoma cells to anti-PD-L1 antibody via PD-1/PD-L1-mediated lymphoma cell interaction with CD8+T cells. Mol Cancer. 2019;18:54.
Okamoto A, Yanada M, Inaguma Y, Tokuda M, Morishima S, Kanie T, et al. The prognostic significance of EBV DNA load and EBER status in diagnostic specimens from diffuse large B-cell lymphoma patients. Hematol Oncol. 2017;35:87–93.
Park S, Lee J, Ko YH, Han A, Jun HJ, Lee SC, et al. The impact of Epstein-Barr virus status on clinical outcome in diffuse large B-cell lymphoma. Blood. 2007;110:972–8.
Glažar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. RNA. 2014;20:1666–70.
Zhang X, Liu X, Jing Z, Bi J, Li Z, Liu X, et al. The circINTS4/miR-146b/CARMA3 axis promotes tumorigenesis in bladder cancer. Cancer Gene Ther. 2020;27:189–202.
Shen R, Xu PP, Wang N, Yi HM, Dong L, Fu D, et al. Influence of oncogenic mutations and tumor microenvironment alterations on extranodal invasion in diffuse large B-cell lymphoma. Clin Transl Med. 2020;10:e221.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Armaos A, Colantoni A, Proietti G, Rupert J, Tartaglia GG. catRAPID omics v2.0: going deeper and wider in the prediction of protein-RNA interactions. Nucleic Acids Res. 2021;49:W72–w79.
Ganzina F. 4′-epi-doxorubicin, a new analogue of doxorubicin: a preliminary overview of preclinical and clinical data. Cancer Treat Rev. 1983;10:1–22.
Yazdian-Robati R, Ramezani M, Jalalian SH, Abnous K, Taghdisi SM. Targeted delivery of epirubicin to cancer cells by polyvalent aptamer system in vitro and in vivo. Pharm Res. 2016;33:2289–97.
Zhang Q, Luo D, Xie Z, He H, Duan Z. The oncogenic role of miR-BART19-3p in Epstein-Barr virus-associated diseases. Biomed Res Int. 2020;2020:5217039.
Zhao Z, Liu W, Liu J, Wang J, Luo B. The effect of EBV on WIF1, NLK, and APC gene methylation and expression in gastric carcinoma and nasopharyngeal cancer. J Med Virol. 2017;89:1844–51.
Clevers H, Nusse R. Wnt/β-catenin signaling and disease. Cell. 2012;149:1192–205.
Arnaiz E, Sole C, Manterola L, Iparraguirre L, Otaegui D, Lawrie CH. CircRNAs and cancer: biomarkers and master regulators. Semin Cancer Biol. 2019;58:90–9.
Wang Y, Mo Y, Gong Z, Yang X, Yang M, Zhang S, et al. Circular RNAs in human cancer. Mol Cancer. 2017;16:25.
Ungerleider N, Concha M, Lin Z, Roberts C, Wang X, Cao S, et al. The Epstein Barr virus circRNAome. PLoS Pathog. 2018;14:e1007206.
Cui C, Yang J, Li X, Liu D, Fu L, Wang X. Functions and mechanisms of circular RNAs in cancer radiotherapy and chemotherapy resistance. Mol Cancer. 2020;19:58.
Hu Y, Gu J, Shen H, Shao T, Li S, Wang W, et al. Circular RNA LARP4 correlates with decreased Enneking stage, better histological response, and prolonged survival profiles, and it elevates chemosensitivity to cisplatin and doxorubicin via sponging microRNA-424 in osteosarcoma. J Clin Lab Anal. 2020;34:e23045.
Liang Y, Song X, Li Y, Su P, Han D, Ma T, et al. circKDM4C suppresses tumor progression and attenuates doxorubicin resistance by regulating miR-548p/PBLD axis in breast cancer. Oncogene. 2019;38:6850–66.
Stathis A, Zucca E. Use of epirubicin to treat diffuse large B-cell lymphoma. Lancet Haematol. 2019;6:e292–3.
Riley KJ, Rabinowitz GS, Yario TA, Luna JM, Darnell RB, Steitz JA. EBV and human microRNAs co-target oncogenic and apoptotic viral and human genes during latency. EMBO J. 2012;31:2207–21.
Sakamoto K, Sekizuka T, Uehara T, Hishima T, Mine S, Fukumoto H, et al. Next-generation sequencing of miRNAs in clinical samples of Epstein-Barr virus-associated B-cell lymphomas. Cancer Med. 2017;6:605–18.
Motsch N, Alles J, Imig J, Zhu J, Barth S, Reineke T, et al. MicroRNA profiling of Epstein-Barr virus-associated NK/T-cell lymphomas by deep sequencing. PLoS One. 2012;7:e42193.
Stefanski CD, Prosperi JR. Wnt-independent and Wnt-dependent effects of APC loss on the chemotherapeutic response. Int J Mol Sci. 2020;21(21):7844.
Saito-Diaz K, Benchabane H, Tiwari A, Tian A, Li B, Thompson JJ, et al. APC inhibits ligand-independent Wnt signaling by the Clathrin endocytic pathway. Dev Cell. 2018;44:566–581.e568.
Liu Q, Shuai M, Xia Y. Knockdown of EBV-encoded circRNA circRPMS1 suppresses nasopharyngeal carcinoma cell proliferation and metastasis through sponging multiple miRNAs. Cancer Manag Res. 2019;11:8023–31.
Zheng X, Huang M, Xing L, Yang R, Wang X, Jiang R, et al. The circRNA circSEPT9 mediated by E2F1 and EIF4A3 facilitates the carcinogenesis and development of triple-negative breast cancer. Mol Cancer. 2020;19:73.
Ivanov A, Memczak S, Wyler E, Torti F, Porath HT, Orejuela MR, et al. Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. Cell Rep. 2015;10:170–7.
Qiu J, Smith P, Leahy L, Thorley-Lawson DA. The Epstein-Barr virus encoded BART miRNAs potentiate tumor growth in vivo. PLoS Pathog. 2015;11:e1004561.
This study was supported, in part, by research funding from the National Natural Science Foundation of China (82130004 and 81830007 to W.-L.Z., U2005204 and 81870135 to J.-D.H., and 81600155 to Z.-X.Y.), Clinical Research Plan of Shanghai Hospital Development Center (SHDC2020CR1032B to W.-L.Z.), Shanghai Rising-Star Program (19QA145600 to Z.-X.Y.), Municipal Human Resources Development Program for Outstanding Young Talents in Medical and Health Sciences in Shanghai (2017YQ075 to Z.-X.Y.), Program for Innovative Research Team in Science and Technology in Fujian Medical University (2017XQ2030 to C.-X.Z. and 2019QH2012 to J.-J.W), Joint Funds for the Innovation of Science and Technology, Fujian province (2016Y9032 to J.-D.H.).
Ethics approval and consent to participate
The study was approved by the Shanghai Rui Jin Hospital Review Board with informed consent obtained from all subjects in accordance with the Declaration of Helsinki.
Consent for publication
All authors agree to submit for consideration for publication in the journal. All authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The primer sequences used for real-time quantitative PCR.
Construction and verification for circEAF2 plasmid.
Sequence information of miRNA mimics and circ siRNA.
Sequence information of biotinylated probes.
Verification of EBV-related circRNAs in B-lymphoma cell lines.
Verification of EBV-related circRNAs in tissues samples of 54 DLBCL patients.
Tumor EBV DNA is associated with serum EBV DNA an d clinical outcomes in DLBCL patients.
CircEAF2 may be regulated by EBV infection through RNA splicing.
About this article
Cite this article
Zhao, Cx., Yan, Zx., Wen, Jj. et al. CircEAF2 counteracts Epstein-Barr virus-positive diffuse large B-cell lymphoma progression via miR-BART19-3p/APC/β-catenin axis. Mol Cancer 20, 153 (2021). https://doi.org/10.1186/s12943-021-01458-9
- Epstein-Barr virus
- Diffuse large B-cell lymphoma
- Wnt signaling pathway