Skip to main content
  • Primary research
  • Open access
  • Published:

CircADARB1 serves as a new biomarker in natural killer T-cell lymphoma and a potential regulator of p-Stat3

Abstract

Background

Natural killer/T-cell lymphoma (NKTCL) is a rare and aggressive subtype of Non-Hodgkin’s Lymphoma. CircRNA has shown great potential to become a biomarker in plasma. In this study, we aimed to determine circRNA for its diagnostic and prognostic value and biological function in NKTCL.

Method

The circRNA microarray of plasma from NKTCL patients and healthy donors were conducted. The relative expressions of target circRNA were verified by qRT-PCR. We conducted function experiments in vitro and in vivo. Bioinformatics predicted the target miRNA of the target circRNA and the binding site was detected by the dual luciferase report assay. Downstream target protein was predicted and detected by western blot in vitro and immunohistochemistry in vivo.

Result

By analyzing the plasma circRNA microarrays in NKTCL, 6137 circRNAs were up-regulated and 6190 circRNAs were down-regulated. The relative expressions of circADARB1 were significantly higher in NKTCL patients. The knockdown of circADARB1 inhibited proliferation of NKTCL cells in vitro and in vivo. CircADARB1 could bind to miR-214-3p in the downstream and regulate the expression of p-Stat3. In nude mice tumor tissue, p-Stat3 was under-expressed in the circADARB1 knockdown group.

Conclusion

CircADARB1 was highly expressed in NKTCL plasma and circADARB1 was a potential biomarker to assist diagnosis and predict the response in NKTCL. CircADARB1 bound up to miR-214-3p and regulated p-Stat3.

Introduction

85–90% of Non-Hodgkin’s Lymphoma (NHL) originates from B lymphocytes, the remaining non-Hodgkin Lymphoma originates from T lymphocytes and Natural Killer (NK) lymphocytes [1]. Natural killer/T-cell lymphoma (NKTCL) is a rare and aggressive subtype of NHL originating from NK lymphocyte and cytotoxic T lymphocyte [2] and associated with Epstein-Barr virus infection [3]. NKTCL is endemic in East Asia and Latin America while rare in Europe and South America [4], usually associated with a poor outcome [5]. Due to its atypical early clinical manifestations, including general symptoms such as fever, night sweats, and fatigue, as well as related symptoms involving organs. The current diagnosis mainly depends on immunohistochemical staining, and the immunohistochemical characteristics of NKTCL cells include CD2 and CD56 positive, CD3 negative; cytoplasmic CD3ε and cytotoxic molecules (perforin, granzyme, TIA1) positive [6]. The prognostic index of natural killer cell lymphoma (PINK) and PINK-E [7] includes age greater than 60 years, stage 3 or 4 stage, involving distant lymph nodes and non-nasal diseases, and the addition of EBV DNA testing .

However, it is difficult to distinguish NKTCL since the early symptoms are not typical. Tumor tissue necrosis is common, which makes sampling difficult and misdiagnosed. In China, patients with pathological test of NKTCL more than once account for 37.1% [8]. In Western countries, the median time from symptom onset to diagnosis is 5 months, the longest is 36 months, and patients with more than once pathological examination account for 48% of patients with NKTCL [9]. Therefore, it is imperative to find new markers to assist in tumor diagnosis.

In addition, efficacy prediction and therapeutic targeting are extremely important in molecular tumor research. Wang et al. found that about 95% of ENKTL cases were CD38 positive and that patients with strong CD38 expression had a low CR rate after chemotherapy [10]. A subsequent clinical trial showed that daratumumab, a monoclonal antibody against CD38, achieved an ORR of only 25% in patients with R/R NKTCL, and none of the patients achieved CR [11]. Although there have been many molecular studies on NKTCL, there is still a lack of validated efficacy predictive biomarkers and therapeutic targets.

Circular (circ)RNA specializes in forming covalently closed continuous loops with back-splicing events [12]. The majority of circRNAs are generated from coding exons (usually between one and five exons) [13, 14]. CircRNAs were initially considered as RNA by-products [15]. Recent studies have revealed biological functions as miRNA sponges [14], regulating transcription or splicing [16, 17], and interacting with RNA binding proteins (RBPs) [18].

CircRNA has shown great potential to become a biomarker because of its conserved, stable, and stage- and tissues-specific characteristic [19, 20]. In addition, circRNAs exist widely in human body, not only in tissues, but also in the human peripheral blood [21, 22] and saliva [23] where can be collected easily [24]. In previous researches, several circRNAs in B-cell lymphoma [25, 26] and T-cell lymphoma [27] have been emphasized. However, clinical references for circRNAs profile and biological functions in NKTCL are not clear. In this study, we focused on circRNA for its diagnostic and prognostic value, and biological function in NKTCL.

Materials and methods

Patients and samples

The plasma samples of the control group involved in this experiment were obtained from healthy donors, and the plasma samples of the experimental group were obtained from patients with pathological diagnosis of NK/T-cell lymphoma. Inclusion criteria for the experimental group: (1) those who attended the First Affiliated Hospital of Zhengzhou University between 2012 and 2019, (2) patients with confirmed NK/T-cell lymphoma diagnosed or consulted by the Department of Pathology of the First Affiliated Hospital of Zhengzhou University, (3) availability of all clinical data, and (4) plasma samples collected at the first admission and signed an informed consent form. Exclusion criteria: (1) patients were previously treated with antineoplastic therapy, (2) patients with an unclear pathological diagnosis or a diagnosis of mixed malignant lymphoma. Clinical data such as age, gender, B symptoms, staging [28], PINK-E, ECOG score, response to therapy and follow-up information were collected.

The initial study involved six patients with a recent diagnosis of NKTCL and six healthy donors as negative controls. Validation experiments involved 50 patients diagnosed with NKTCL and 50 healthy donors as negative controls. Plasma samples from patients were collected before chemotherapy and preserved at − 80 °C until analyzed.

Total RNA extraction

Total RNA was extracted from plasma samples using TRIzol LS reagent (Invitrogen, Carlsbad, CA, USA) and from cells using TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. The purity and concentration of RNA samples were determined with a Nanodrop 1000 spectrophotometry (Thermo Scientific, Wilmington, DE, USA).

Circular RNA preparation and microarray

Total RNAs were digested with Rnase R (Epicentre, Inc.) to remove linear RNAs and enrich circular RNAs. Then, the enriched circular RNAs were amplified and transcribed into fluorescent cRNA utilizing a random priming method (Arraystar Super RNA Labeling Kit; Arraystar). The labeled cRNAs were hybridized onto the Arraystar Human circRNA Array V2 (8 × 15 K, Arraystar). After having washed the slides, the arrays were scanned by the Agilent Scanner G2505C.

Quantitative reverse transcription polymerase chain reaction (qRT-PCR)

cDNA was synthesized by reverse transcription using the PrimeScript RT Master Mix with random primers, and qRT-PCR was performed according to the manufacturer’s protocols (TaKaRa Bio, Shiga, Japan). Levels of circRNA and mRNA were normalized to GAPDH. Relative expression level was calculated with the 2−ΔΔCt method.

Cell culture

HEK293T and YT cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) and Roswell Park Memorial Institute (RPMI) 1640, respectively (ThermoFisher, MA, USA) containing 10% fetal bovine serum (Clark Bioscience, VA, USA). SNK-6 cell line was cultured in RPMI 1640 containing 5% advanced cell culture supplement (AventaCell BioMedical Co., GA, USA) and 1000 U/ml recombinant human IL-2. In addition, 100 U/ml penicillin and streptomycin (ThemoFisher Scientific, MA, USA) were added in cell culture medium at 37 °C with 5% CO2.

CircRNA plasmid construction and stable transfection

circADARB1 cDNA was synthesized and cloned into the GV493 vector (Genechem, Shanghai, China) to construct three short hairpin RNA plasmids. The plasmid constructions were confirmed by sequencing, then transfected into HEK293T cells for lentiviral packaging using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s protocol. After the lentivirus harvested, NKTCL cell lines YT and SNK-6 were transfected and selected for seven days with 3 μg/mL puromycin.

Cell viability assay

YT and SNK-6 cells were cultured in 96-well plates coated with 20 μl of Cell Counting Kit-8 (CCK-8, Dojindo, Tokyo, Japan) per well and incubated at 37 °C supplied with 5% CO2 for 2 h. Next, cell viability was presented by the wavelength of OD450 nm with a Multiskan FC microplate reader (Thermo Scientific, Waltham, MA, USA).

Cell apoptosis analysis

Cell apoptosis assay was performed using the Annexin V-APC Apoptosis Detection Kit (BD Biosciences, USA) following the manufacturer’s instructions. Cells were stained with Annexin V-APC and 7AAD, followed by flow cytometry in a FACS Canto II (BD Biosciences).

Bioinformatics analysis

Targetscan (http://www.targetscan.org) [29], ENCORI (formerly known as starBase) [30] and miRanda [31] for bioinformatics analysis and prediction.

Dual luciferase reporter assay

The different fragment sequences were synthesized and then inserted into the psiCHECK2 vector (Hanbio, China). Luciferase activity was assessed using the Luciferase Assay Reagent II (LAR II) (Luciferase Assay Reagent, Progema) according to the manufacturer’s instructions.

Western blotting

Total protein was extracted from cells using RIPA with Protease and Phosphatase Inhibitor cocktail (CWBIO, China) and EDTA cocktail. The protein concentrations were detected by BCA Protein Assay Kit (CWBIO, China) according to the manufacturer’s instructions. Electrophoresis was performed with sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE), and protein was transferred onto polyvinylidene fluoride (PVDF) membranes (Amersham Biosciences, Piscataway, NJ, USA). Primary antibodies used were phospho-Stat3 (1:1000), Bax (1:1000), GAPDH (1:2000) (Cell Signaling Technology, Boston, MA, USA) and the secondary antibodies (ProteinTech, Chicago, IL, USA). Detection was carried out with the eECL Western Blot Kit (CWBIO, China). The band images were digitally captured and quantified with a ChemiDoc. XRS + system (Bio-Rad Laboratories, Hercules, CA, USA).

In vivo tumor xenograft

Ten female BALB/c nude mice (6-week-old) were maintained under specific pathogen-free conditions from the Nanjing Medical Experimental Animal Care Commission. Transfected YT cells were harvested into two groups. 1 × 107 cells were subcutaneously injected into a single flank of each mouse. The mice were CO2 euthanasia at a rate of 30% chamber volume/min using an euthanasia system according to AVMA Guidelines for the Euthanasia of Animals and all procedures complied with the animal care guidelines from the First Affiliated Hospital of Zhengzhou University.

Immunohistochemistry

Immunohistochemistry was performed as previously described [29]. Sections were incubated with antibodies: phospho-Stat3 (Ser473) (diluted 1:100; Cell Signaling Technology) and Ki-67 (diluted 1:50). Expression levels were scored using the mean optical density (AOD) by Image-Pro Plus 6.0.

TUNEL assay

We used TUNEL staining using a TUNEL Apoptosis Assay Kit (keyGen BioTECH, Nanjing, China) following the manufacturer’s instruction.

Results

CircRNA profiling of NKTCL patients’ plasma

We performed microarrays to characterize the profiles of circRNAs in the plasma from 6 NKTCL patients and six healthy donors (Additional file 2). The results showed that 6137 up-regulated circRNAs and 6190 down-regulated circRNAs with differential expression (FC value > 1.0). Scatter plot, Volcano map and cluster analysis plot were drawn according to the results (see Fig. 1A–C). According to the conditions of FC value > 2.0 and P < 0.05, the differentially expressed circRNA was further screened, of which 170 were up-regulated and 188 were down-regulated. The top five upregulated circRNAs are hsa_circ_0087631, hsa_circ_0045932, hsa_circ_0052372, hsa_circ_0005037, hsa_circ_0007059, of which circADARB1 (ID hsa_circ_0005037 in circbase) was stably expressed and the difference was significant in small sample qPCR experiment (P = 0.029).

Fig. 1
figure 1

Differential expression of circRNA in NKTCL patients and healthy controls. A Scatter plot; The CircRNAs above the top green line and below the bottom green line indicated more than 2.0-fold change of circRNAs between the two groups. B Volcano plot; The vertical lines correspond to a 2.0-fold increase and decrease, respectively, and the horizontal line represents a p-value of 0.05. The red dots represent circRNAs with statistically significant differential expression. C Hierarchical clustering; Hierarchical clustering was performed based on all targets value. The result of hierarchical clustering shows a distinguishable circRNA expression profiling among samples

CircADARB1 is highly expressed in the plasma of NKTCL patients

The expressions of circADARB1 in plasma of 50 NKTCL patients and 50 healthy controls were measured by qRT-PCR. The results suggested that circADARB1 was significantly up-regulated in the plasma of NKTCL patients (P < 0.001) (Fig. 2A). The receiver operating characteristic curve (ROC) was drawn, and the area under the curve (Area Under Curve, AUC) was calculated. The results showed (Fig. 2B) that the AUC was 0.742 ± 0.050 (P < 0.001), the 95% confidence interval was (0.644, 0.840), the sensitivity was 0.860, and the specificity was 0.600. The results showed that the relative expression of circADARB1 in plasma may assist the diagnosis of NKTCL.

Fig. 2
figure 2

Analysis of relative expression of circADARB1 and clinical data of NKTCL patients. A Comparison of the relative expression of circADARB1 in the plasma of NKTCL patients and healthy controls; B ROC curve of circADARB1; C Comparison of the relative expression of circADARB1 in NKTCL patients with different efficacy; D three-year survival curve. *P < 0.05, ***P < 0.001

We next evaluated the association between circADARB1 and clinical pathological parameters. The results showed that the relative expression of circADARB1 was no statistically difference grouped by age, gender, B symptoms, stage, PINK-E prognostic index and ECOG score (Table 1). But the expression of circADARB1 was associated with efficacy (classified into CR + PR and SD + PD) (P = 0.041, Fig. 2C), and the relative expression of circADARB1 was higher in patients with SD and PD. The results suggest that circADARB1 may be a potential biomarker for predicting efficacy. Using the median expression level of circADARB1 as cutoff value, 50 NKTCL patients were divided into two groups: high expression and low expression of circADARB1. Survival analysis was conducted based on follow-up data, and a 3-year survival curve was drawn (Fig. 2D). The overall survival rates of the high- and low-expression group were 60.1% and 82.9%, respectively, but the difference was not statistically significant (P = 0.153).

Table 1 Relationship between the relative expression of circADARB1 in plasma of NKTCL patients and the characteristics of clinical data

Knockdown of circADARB1 inhibits the proliferation and promotes apoptosis of NKTCL cell lines in vitro

According to the results of CCK-8, it was found that the proliferation of YT cells after the knockdown of circADARB1 was significantly lower than that of the control group (Fig. 3A, P = 0.001 on day 7), and similar results were found in SNK-6 (Fig. 3A, P = 0.043 on day 5, P = 0.047 on day 7). The results showed that after knockdown of circADARB1, the proliferation of NKTCL cells decreased.

Fig. 3
figure 3

Functional experiments after transfection of NKTCL cells. A CCK-8 experiment in YT and SNK-6 (from 1 to 7 days) cells; B Apoptosis experiment results of NC-circADARB1 and sh-circADARB1 transfected YT and SNK-6 cells. C Western Blot detection of Bax in NC-circADARB1 and sh-circADARB1 transfected YT and SNK-6 cells. Note: sh-circADARB1 stands for the knockdown of circADARB1; NC-circADARB1 stands for transfected with control vectors. *P < 0.05, **P < 0.01

The cell apoptosis analysis showed that the apoptosis rate of YT cells [(14.63 ± 0.3033) %] with knockdown of circADRB1 was significantly higher than that of control group [(7.803 ± 0.8872) %] (Fig. 3B, P = 0.002). The results were similar in SNK-6 (Fig. 3B, P = 0.010). The western blot results showed an increase in Bax protein expression after knockdown of circADARB1 (Fig. 3C). The results indicated that knockdown of circADARB1 promoted apoptosis of NKTCL cells.

Knockdown of circADARB1 inhibits the proliferation of NKTCL cells in vivo

The transfected YT cells were inoculated into the right armpit of nude mice to establish the NKTCL tumor-bearing mouse models. The tumor volumes in the knockdown of circADARB1 group were significantly smaller (Fig. 4A), and weight of tumors was lighter in comparison with control group (P = 0.019, Fig. 4B). The results show that the knockdown of circADARB1 can inhibit the proliferation of NKTCL cells in vivo. The relative expression of circADARB1 from knockdown group was significantly lower than that of the control group (P < 0.001, Additional file 1: Fig. S1A).

Fig. 4
figure 4

Knockdown of circADARB1 inhibits the proliferation of NKTCL cells in vivo; A Volume of subcutaneous tumor formation in nude mice between NC-circADARB1 and sh-circADARB1 group; B Comparison of the quality of subcutaneous tumor formation between NC-circADARB1 and sh-circADARB1 group; C Ki-67 and TUNEL assay of nude mice tumor formation experiment between NC-circADARB1 and sh-circADARB1 group. sh-circADARB1 stands for the knockdown of circADARB1; NC-circADARB1 stands for transfected with control vectors.*P < 0.05

The results of TUNEL assay (Fig. 4C) showed that the tissues from knockdown group was positive while the control group was weakly positive (Additional file 1: Fig. S1B, P = 0.013). The Ki-67 in the knockdown group was weakly positive while the Ki-67 in the control group was positive (Fig. 4C, Additional file 1: Fig. S1C, P = 0.032).

circADARB1 serves as a sponge for miR-214 and downstream signaling STAT3

Through TargetScan and miRanda bioinformatics analysis, it is predicted that circADARB1 can combine with miR-214-3p (Fig. 5A) to act as a competitive endogenous RNA and play the role of “miRNA sponge”.

Fig. 5
figure 5

Downstream of circADARB1. A Analysis of circADARB1 and miR-214-3p binding sites; B Dual luciferase reporter design for miR-214-3p and circADARB1 target site; C Relative luciferase activity reported by circADARB1 and miR-214-3p dual luciferase reporter assay; D Western Blot detection of p-Stat3; E p-Stat3 immunohistochemistry in animal experiments. Note: sh-circADARB1 stands for the knockdown of circADARB1; NC-circADARB1 stands for transfected with control vectors. *P < 0.05, **P < 0.01

To further explore the relationship, wild-type and mutant plasmids of circADARB1 (Fig. 5B) were designed based on the co-binding sites of circADARB1 and miR-214-3p, and a dual luciferase reporter assay was conducted. The results showed that after co-transfection of wild-type circADARB1 and miR-214-3p mimic, the relative luciferase activity was significantly reduced (P = 0.010) comparing with co-transfection of mutant circADARB1 and miR-214-3p mimic group. There was no significant difference between the mutant and control mimic groups (Fig. 5C). The experimental results indicated that circADARB1 interacts with miR-214-3p.

We analyzed the targets of miR-214-3p using bioinformatics tools ENCORI and Targetscan and predicted that miR-214-3p could bind to STAT3 in combination with the analysis of important NKTCL signaling pathways in previous studies.

Then we detected p-Stat3 expression by Western Blot (Fig. 5D). The results showed that after knockdown of circADARB1, the expression of p-Stat3 in both YT and SNK-6 was reduced. Immunohistochemical detection of p-Stat3 was performed in nude mice in vivo model tumor tissue. It was found that p-Stat3 was strongly positive in control group and was reduced in knockdown of circADARB1 group (Fig. 5E).

Discussion

NKTCL is a rare subtype of NHL with the characteristics of aggressive invasiveness and rapid progress. The occurrence of NKTCL has obvious geographical characteristics, and it is more common in East Asia [1]. Its incidence in China is about 11% of lymphomas [32]. NKTCL has a poor prognosis, with a median survival of 59 months and a median progression-free survival of 20 months [33]. CircRNA is a member of non-coding RNA. Due to its circular structure, it has unique biological properties that make it suitable as a biomarker. Therefore, this study aims to explore the possibility of circular RNA in NKTCL as a new molecular biomarker.

There are a lot of abnormally expressed circRNAs in the plasma or serum of tumor patients. Li et al. [34] found that many circRNAs are presented in serum exosomes and thus have the potential to be circulating biomarkers in cancer diagnosis. Hang et al. [35] found circFARSA may promote lung cancer development by sponging miR-330-5p/miR-326 and attenuating their repression of FASN. They also found higher abundance of circFARSA in plasma than in plasma exosomes. They therefore speculated that there may be other ways to release circFARSA into the circulation [35]. Fan et al. [36] reported that hsa_circ_0001946 is highly expressed in patients’ preoperative plasma and correlates with recurrence rates in ESCC patients. Xu et al. [37] figured out that hsa_circ_0005037 (circADARB1) is downregulated in colorectal cancer tissues. In a previous study by our group, we found that circCDYL was highly expressed in the plasma of mantle cell lymphoma patients compared to healthy donors [38].

In this study, the expression levels of circADARB1 were significantly higher in plasma of NKTCL patients and correlated with the response to therapy, i.e., the relative expression of circADARB1 was higher in plasma of patients with SD or PD. The results of our study suggest that circADARB1 may be an auxiliary diagnostic and predictive marker for efficacy.

To explore the biological function of circADARB1 in NKTCL, we performed experiments in vitro and in vivo. We found that the knockdown of circADARB1 inhibits the proliferation of NKTCL cells and promotes apoptosis. The result was similar in subcutaneous tumor xenograft in nude mice, which indicates that circADARB1 has a role in promoting tumor growth in NKTCL, and circADARB1 has the potential as a therapeutic target for NKTCL.

In this study, circADARB1 was found to bind to miR-214-3p by bioinformatic analysis, as evidenced by the use of a dual luciferase reporter assay. Previous experiments have shown that miR-214 is significantly increased in activated T cells and can reduce the expression of Pten [39]. miR-214 can bind to the 3′UTR of Twist mRNA, thereby inhibiting the expression of epithelial-mesenchymal transition related gene TWIST and regulating the metastasis of intrahepatic cholangiocarcinoma [40]. The inhibition of NF-κB increase the expression of miR-214 in hepatocellular carcinoma, suggesting that NF-κB maybe a negative regulator of miR-214 expression [41]. In different tumors, miR-214 exhibits dual effects of promotion or suppression of malignant tumor, indicating that it has a complex biological mechanism during tumorigenesis [42]. miR-214 is significantly increased in patients with Sezary Syndrome of cutaneous T-cell lymphoma [43], and miR-214 inhibits apoptosis of cutaneous T-cell lymphoma cell line HUT78 [44]. The level of miR-214 is significantly higher than that of healthy donors, and miR-214 is expected to be a therapeutic target for cutaneous T-cell lymphoma [45]. Currently, there are no reports on miR-214 in NKTCL yet.

For the downstream proteins, we predicted the target of miR-214-3p through a combination of bioinformatics and the results of important signaling pathways in the previous NKTCL research. We found the binding site between miR-214-3p and STAT3. p-Stat3 is the activated form of Stat3. In this study, knocking down circADARB1 would reduce the expression of p-Stat3. The mechanism by which p-Stat3 exerts its biological functions. Previous studies have found that in the analysis of transcriptome differences between NKTCL and normal NK cells, STAT3 is a highly differentially expressed transcription factor in NKTCL [46]. In NKTCL, p-Stat3 is highly expressed in patient samples and cell lines, indicating that p-Stat3 plays an important role in the occurrence and development of NKTCL [47]. The main reason for the abnormally high expression of p-Stat3 is the mutation of STAT3 [48] and the activation of JAK/STAT3 pathway [49]. The mutation rate of the SH2 domain of STAT3 in NKTCL was 5.9%, and cell lines with this mutation were found to up-regulate p-Stat3 [48]. Another study obtained results with higher mutation rates. The rate of missense mutations of STAT3 in NKTCL was 18.9%, and the silent mutations was 5.4%. Among all patients with STAT3 missense mutations, the expression of p-Stat3 shows higher [49]. The expression profile of JAK/STAT3 pathway in NKTCL patient tissues and normal NK cells was significantly different, and further affected multiple cell functions [50]. These studies have shown that the abnormally high expression of p-Stat3 plays an important role in the development of NKTCL. In peripheral T-cell lymphoma cell line Hut78, knocking down STAT3 reduced Bcl-XL expression and promoting apoptosis [51]. In ALK + anaplastic large cell lymphoma, a strong correlation was found between the positive expression of p-Stat3 and the expression of Survivin [52]. STAT3, as an oncogene of transcription factors, can regulate the expression of MYC, VEGFA, BCL2L1, BIRC5, HGF, MMP2, MMP9, and CDK5, suggesting its pathogenesis in NKTCL [50]. In this study, the results showed that circADARB1 regulated miR-214-3p and p-Stat3, thereby promoting the proliferation of NKTCL cells. It supplemented and expanded the previous researches on p-Stat3 in NKTCL.

This study is an attempt to find biomarkers in plasma circular RNA. It is proposed that circADARB1 may become a potential biomarker of NKTCL, which is a supplement to previous studies on NKTCL. However, this study still has certain limitations. This trial lacks the design of prospective clinical studies to verify the validity of circRNA as a biomarker. Since this study only involved Chinese NKTCL patients, there might be a population bias. In addition, the other targets or the biological mechanism of miR-214 is also not fully understood. Other biological functions of circADARB1 besides the function of “miRNA sponge” also need to be further explored.

In conclusion, circADARB1 is highly expressed in the plasma of NKTCL patients and is a potential biomarker for adjuvant diagnosis and prediction of efficacy. Knockdown of circADARB1 inhibits NKTCL cell proliferation in vitro and in vivo. circADARB1 binds up to miR-214-3p and regulates p-Stat3.

Availability of data and materials

We have attached the original data (Additional file 2).

References

  1. Shankland KR, Armitage JO, Hancock BW. Non-Hodgkin lymphoma. Lancet. 2012;380(9844):848–57.

    Article  PubMed  Google Scholar 

  2. Yamaguchi M, Suzuki R, Oguchi M. Advances in the treatment of extranodal NK/T-cell lymphoma, nasal type. Blood. 2018;131(23):2528–40.

    Article  CAS  PubMed  Google Scholar 

  3. Matutes E. The 2017 WHO update on mature T- and natural killer (NK) cell neoplasms. Int J Lab Hematol. 2018;40 Suppl 1:97–103.

    Article  CAS  PubMed  Google Scholar 

  4. Kumai T, Kobayashi H, Harabuchi Y. Novel targets for natural killer/T-cell lymphoma immunotherapy. Immunotherapy. 2016;8(1):45–55.

    Article  CAS  PubMed  Google Scholar 

  5. Mei M, Wang Y, Zhang M. Causes of mortality in cases with extra nodal natural killer/T-cell lymphoma, nasal type: a cohort study. PLoS ONE. 2019;14(4): e0214860.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Somasundaram N, Lim JQ, Ong CK, Lim ST. Pathogenesis and biomarkers of natural killer T cell lymphoma (NKTL). J Hematol Oncol. 2019;12(1):28.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Kim SJ, Yoon DH, Jaccard A, et al. A prognostic index for natural killer cell lymphoma after non-anthracycline-based treatment: a multicentre, retrospective analysis. Lancet Oncol. 2016;17(3):389–400.

    Article  CAS  PubMed  Google Scholar 

  8. Liu PL, Cheng ZX, Lu MP, Zhang LQ, Chen HB. Misdiagnosis analysis: 120 patients with nasal extranodal NK/T cell lymphoma in head and neck. Lin Chung Er Bi Yan Hou Tou Jing Wai Ke Za Zhi. 2020;34(1):73–8.

    CAS  PubMed  Google Scholar 

  9. Haverkos BM, Pan Z, Gru AA, et al. Extranodal NK/T cell lymphoma, nasal type (ENKTL-NT): an update on epidemiology, clinical presentation, and natural history in North American and European cases. Curr Hematol Malig Rep. 2016;11(6):514–27.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wang L, Wang H, Li PF, et al. CD38 expression predicts poor prognosis and might be a potential therapy target in extranodal NK/T cell lymphoma, nasal type. Ann Hematol. 2015;94(8):1381–8.

    Article  CAS  PubMed  Google Scholar 

  11. Huang H, Zhu J, Yao M, et al. Daratumumab monotherapy for patients with relapsed or refractory natural killer/T-cell lymphoma, nasal type: an open-label, single-arm, multicenter, phase 2 study. J Hematol Oncol. 2021;14(1):25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Lasda E, Parker R. Circular RNAs: diversity of form and function. RNA. 2014;20(12):1829–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wilusz JE, Sharp PA. Molecular biology. A circuitous route to noncoding RNA. Science. 2013;340(6131):440–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.

    Article  CAS  PubMed  Google Scholar 

  15. Kolakofsky D. Isolation and characterization of Sendai virus DI-RNAs. Cell. 1976;8(4):547–55.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang Y, Zhang XO, Chen T, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51(6):792–806.

    Article  CAS  PubMed  Google Scholar 

  17. Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE. 2012;7(2): e30733.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ashwal-Fluss R, Meyer M, Pamudurti NR, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014;56(1):55–66.

    Article  CAS  PubMed  Google Scholar 

  19. Chen LL, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12(4):381–8.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Rybak-Wolf A, Stottmeister C, Glazar P, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58(5):870–85.

    Article  CAS  PubMed  Google Scholar 

  21. Haque S, Ames RM, Moore K, et al. circRNAs expressed in human peripheral blood are associated with human aging phenotypes, cellular senescence and mouse lifespan. GeroScience. 2020;42(1):183–99.

    Article  CAS  PubMed  Google Scholar 

  22. Pan B, Qin J, Liu X, et al. Identification of serum exosomal hsa-circ-0004771 as a novel diagnostic biomarker of colorectal cancer. Front Genet. 2019;10:1096.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ghods FJ. Circular RNA in saliva. Adv Exp Med Biol. 2018;1087:131–9.

    Article  CAS  Google Scholar 

  24. Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846–58.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Hu Y, Zhao Y, Shi C, et al. A circular RNA from inhibits the proliferation of diffuse large B-cell lymphoma by inactivating Wnt/β-catenin signaling via interacting with TET1 and miR-888. Aging. 2019;11(19):8068–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Dahl M, Daugaard I, Andersen MS, et al. Enzyme-free digital counting of endogenous circular RNA molecules in B-cell malignancies. Lab Invest. 2018;98(12):1657–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  CAS  PubMed  Google Scholar 

  28. Hong H, Li Y, Lim ST, et al. A proposal for a new staging system for extranodal natural killer T-cell lymphoma: a multicenter study from China and Asia Lymphoma Study Group. Leukemia. 2020;34(8):2243–8.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4: e05005.

    Article  PubMed Central  Google Scholar 

  30. Li J-H, Liu S, Zhou H, Qu L-H, Yang J-H. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.

    Article  CAS  PubMed  Google Scholar 

  31. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5(1):R1.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Sun J, Yang Q, Lu Z, et al. Distribution of lymphoid neoplasms in China: analysis of 4,638 cases according to the World Health Organization classification. Am J Clin Pathol. 2012;138(3):429–34.

    Article  PubMed  Google Scholar 

  33. Fox CP, Civallero M, Ko YH, et al. Survival outcomes of patients with extranodal natural-killer T-cell lymphoma: a prospective cohort study from the international T-cell Project. Lancet Haematol. 2020;7(4):e284–94.

    Article  PubMed  Google Scholar 

  34. Li Y, Zheng Q, Bao C, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. 2015;25(8):981–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Hang D, Zhou J, Qin N, et al. A novel plasma circular RNA circFARSA is a potential biomarker for non-small cell lung cancer. Cancer Med. 2018;7(6):2783–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Fan L, Cao Q, Liu J, Zhang J, Li B. Circular RNA profiling and its potential for esophageal squamous cell cancer diagnosis and prognosis. Mol Cancer. 2019;18(1):16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Xu D, Wu Y, Wang X, et al. Identification of functional circRNA/miRNA/mRNA regulatory network for exploring prospective therapy strategy of colorectal cancer. J Cell Biochem. 2020;121:4785–97.

    Article  CAS  Google Scholar 

  38. Mei M, Wang Y, Wang Q, Liu Y, Song W, Zhang M. CircCDYL serves as a new biomarker in mantle cell lymphoma and promotes cell proliferation. Cancer Manag Res. 2019;11:10215–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Jindra PT, Bagley J, Godwin JG, Iacomini J. Costimulation-dependent expression of microRNA-214 increases the ability of T cells to proliferate by targeting Pten. J Immunol. 2010;185(2):990–7.

    Article  CAS  PubMed  Google Scholar 

  40. Li B, Han Q, Zhu Y, Yu Y, Wang J, Jiang X. Down-regulation of miR-214 contributes to intrahepatic cholangiocarcinoma metastasis by targeting Twist. FEBS J. 2012;279(13):2393–8.

    Article  CAS  PubMed  Google Scholar 

  41. Duan Q, Wang X, Gong W, et al. ER stress negatively modulates the expression of the miR-199a/214 cluster to regulates tumor survival and progression in human hepatocellular cancer. PLoS ONE. 2012;7(2): e31518.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Sharma T, Hamilton R, Mandal CC. miR-214: a potential biomarker and therapeutic for different cancers. Future Oncol. 2015;11(2):349–63.

    Article  CAS  PubMed  Google Scholar 

  43. Ballabio E, Mitchell T, van Kester MS, et al. MicroRNA expression in Sezary syndrome: identification, function, and diagnostic potential. Blood. 2010;116(7):1105–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Narducci MG, Arcelli D, Picchio MC, et al. MicroRNA profiling reveals that miR-21, miR486 and miR-214 are upregulated and involved in cell survival in Sézary syndrome. Cell Death Dis. 2011;2(4):e151.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Kohnken R, McNeil B, Wen J, et al. Preclinical Targeting of MicroRNA-214 in Cutaneous T-Cell Lymphoma. J Invest Dermatol. 2019;139(9):1966.e3–74.e3.

    Article  CAS  Google Scholar 

  46. Ng S-B, Selvarajan V, Huang G, et al. Activated oncogenic pathways and therapeutic targets in extranodal nasal-type NK/T cell lymphoma revealed by gene expression profiling. J Pathol. 2011;223(4):496–510.

    Article  CAS  PubMed  Google Scholar 

  47. Coppo P, Gouilleux-Gruart V, Huang Y, et al. STAT3 transcription factor is constitutively activated and is oncogenic in nasal-type NK/T-cell lymphoma. Leukemia. 2009;23(9):1667–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Küçük C, Jiang B, Hu X, et al. Activating mutations of STAT5B and STAT3 in lymphomas derived from γδ-T or NK cells. Nat Commun. 2015;6:6025.

    Article  PubMed  CAS  Google Scholar 

  49. Liu J, Liang L, Li D, et al. JAK3/STAT3 oncogenic pathway and PRDM1 expression stratify clinicopathologic features of extranodal NK/T-cell lymphoma, nasal type. Oncol Rep. 2019;41(6):3219–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Huang Y, de Reynies A, de Leval L, et al. Gene expression profiling identifies emerging oncogenic pathways operating in extranodal NK/T-cell lymphoma, nasal type. Blood. 2010;115(6):1226–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Verma NK, Davies AM, Long A, Kelleher D, Volkov Y. STAT3 knockdown by siRNA induces apoptosis in human cutaneous T-cell lymphoma line Hut78 via downregulation of Bcl-xL. Cell Mol Biol Lett. 2010;15(2):342–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Nasr MR, Laver JH, Chang M, Hutchison RE. Expression of anaplastic lymphoma kinase, tyrosine-phosphorylated STAT3, and associated factors in pediatric anaplastic large cell lymphoma: a report from the children’s oncology group. Am J Clin Pathol. 2007;127(5):770–8.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

YT was kindly adopted from professor Wing C. Chan from city of Hope medical center. SNK-6 was kindly adopted from Dr. Norio Shimizu and Yu Zhang from Chiba University.

Funding

This study was funded by National natural Science Foundation of China (81970184).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: MM and MZ; Data curation: MM; Formal analysis: MM and QW; Funding acquisition: MZ; Investigation: MM and JL; Methodology: MM, YW, WS and ZL, Project administration: MZ and ZL. Resources: MZ and ZL; Software: MM and QW; Supervision: MZ; Validation: MZ; Visualization: MM and MZ; Roles/Writing—original draft: MM; Writing—review and editing: MZ, MM, ZL, JL, QW, YW and WS. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mingzhi Zhang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University (2020-KY-298) and met international standards for patient confidentiality. Written informed consent was obtained from all patients.

Consent for publication

All listed authors have participated in the study and approved the submitted manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Knocking down circADARB1 inhibits the proliferation of NKTCL in vivo; (A) Relative expression of circADARB1 in nude mouse tumor tissue. (B) Average optical density (AOD) of TUNEL assay between NC-circADARB1 and sh-circADARB1 group. (C) AOD of Ki-67 between NC-circADARB1 and sh-circADARB1 group. * P < 0.05, *** P < 0.001.

Additional file 2.

 Raw data for circRNA profiling.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mei, M., Wang, Y., Song, W. et al. CircADARB1 serves as a new biomarker in natural killer T-cell lymphoma and a potential regulator of p-Stat3. Cancer Cell Int 21, 594 (2021). https://doi.org/10.1186/s12935-021-02296-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-021-02296-x

Keywords