1,738
Views
9
CrossRef citations to date
0
Altmetric
Research paper

Identification of circRNA_001846 as putative non–small cell lung cancer biomarker

, , , & ORCID Icon
Pages 8690-8697 | Received 11 Sep 2021, Accepted 06 Oct 2021, Published online: 25 Oct 2021

ABSTRACT

CircRNAs play diverse roles in the regulation of oncogenic processes, yet the roles of these circRNAs in non–small cell lung cancer (NSCLC) remain to be fully clarified. Herein, patterns of circRNA expression in NSCLC tissues and paracancerous tissues were assessed, and the relationships between these circRNAs and NSCLC patient clinical findings were assessed. NSCLC tissues were evaluated using a circRNA microarray approach, with Quantitative real‑time polymerase chain reaction (qPCR) qPCR being used to validate candidate circRNA expression levels in NSCLC patients peripheral blood samples. GEO2R was used to analyze the array data. SPSS23.0, GraphPad Prism, and Sigmaplot were utilized for statistical analyses. Overall, 114 circRNAs that were differentially expressed in NSCLC tissue relative to paracancerous tissue levels were identified, of which 77 and 37 were downregulated and upregulated, respectively. CircRNA_001846 were then chosen based on its differential expression score. Consistent with array findings, serum samples from NSCLC patients exhibited circRNA_001846 upregulation relative to those from healthy individuals. The circRNA_001846 was associated with tumor differentiation, lymph node metastasis, and node metastasis stage. Receiver operating characteristic (ROC) curves analyses revealed that levels of circRNA_001846 in the serum were capable of differentiating between individuals diagnosed with NSCLC and controls at a cut off of 3.9496, yielding respective sensitivity and specificity values of 78.2% and 81.1%, with an area under the curve (AUC) value of 0.872. When combined with carcinoembryonic antigen (CEA), this circRNA achieved an improved AUC value of 0.925. CircRNA_001846 may represent a promising diagnostic biomarker for NSCLC.

Introduction

Non–small cell lung cancer (NSCLC) is among the most common cancers in the world, causing serious harm and mortality in affected patients [Citation1,Citation2]. Approaches to diagnosing and treating this form of malignancy have markedly improved in recent years, such as fungal-derived materials [Citation3], yet patient survival rates remain poor. There is thus a clear need to define a reliable diagnostic biomarker for NSCLC and to explore the value of such biomarkers as targets for therapeutic intervention.

As a recently discovered form of non-coding RNA in the form of a covalently closed loop [Citation4,Citation5], circular RNAs (circRNAs) lack 5ʹ or 3ʹ ends, yet nonetheless harbor sites for microRNA (miRNA) binding [Citation6]. Most circRNAs are derived from exons of protein-coding genes and are produced via a back-splicing process. Owing to their highly stable structural makeup, circRNAs can be readily detected in patient serum samples and in harvested exosomes, making them promising diagnostic targets. Indeed, prior studies have highlighted the value of circRNA-based diagnostic analyses in breast [Citation7,Citation8], colorectal [Citation9,Citation10], liver [Citation11,Citation12], gastric [Citation13,Citation14], and non-small cell lung cancer [Citation15,Citation16]. Functionally, these circRNAs can regulate diverse oncogenic activities such as the migration, proliferation, and survival of tumor cells.

The purpose of this study was to investigate circRNA_001846 as non–small cell lung cancer biomarker. We explored patterns of circRNA_001846 expression in NSCLC patient tumors and paracancerous tissues, and we then evaluated the potential diagnostic utility of an identified NSCLC-related circRNA_001846. Overall, these analyses highlighted the promising value of circRNA_001846 in the diagnostic assessment of potential NSCLC patients.

Materials and methods

Sample collection

The the Ethics Committee of First Affiliated Hospital of Chengdu Medical College approved this study, which was conducted in a manner consistent with the declaration of Helsinki. Samples of serum were collected from 206 NSCLC patients and 206 healthy controls following the provision of informed written consent. Blood was collected into a 2 mL sterile RNase-free tube, and RNA was isolated within 30 minutes. Clinicopathological and molecular characteristics were recorded for each patient as appropriate.

RNA isolation

Trizol (Invitrogen) was used to extract RNA from human samples, after which amirVana Isolation Kit (Ambion, TX, USA) was used based on provided directions. A NanoDrop2000 instrument (Thermo Fisher Scientific) was used to assess RNA quantity and purity, ensuring that OD260/280 ratios were between 1.9 and 2.0. The integrity of isolated RNA was assessed via 1% denatured gel electrophoresis.

Microarray analysis

A previously published Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) dataset (GSE158695) was used to assess circRNA expression patterns in NSCLC, with GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) being used for analysis.

Quantitative real‑time polymerase chain reaction (qPCR)

A Superscript Reverse Transcription System (Invitrogen) was used to prepare cDNA, after which TB Green qPCR Mastermix (TaKaRa, Japan) and a LightCycler® 480 real-time PCR Platform (Roche) were used to conduct qPCR analyses, which were performed in a total 20 μL reaction volume containing 0.8 μL (10 μM) of forward/reverse primers, 10 μL of TB Green qPCR Mastermix, 2 μL of cDNA, and 6.4 μL double-distilled water. Normalization of gene expression was performed using GAPDH. Primers were prepared by Sangon Biotech (Shanghai) Co. Ltd, circRNA_001846-forward: 5ʹ-CGGCCCTAACAGGGCTCTC-3ʹ, circRNA_001846-reverse: 5ʹ-CCTCTGGCCCTAGTCTCAGAC-3ʹ. The ΔΔCT method was used to assess relative gene expression.

Enzyme linked immunosorbent assay (ELISA)

Serum CEA levels were measured with third-generation ELISA kits (Can Ag, Canada) based on provided directions. Samples were analyzed in triplicate.

Statistical analysis

SPSS23.0 (IBM, USA) and GraphPad Prism 7.0 (GraphPad Software, USA) were used to analyze all data, which are given as means ± SD, and were compared via independent samples t-tests or one-way ANOVAs. Receiver operating characteristic (ROC) curves were constructed, and an area under the curve (AUC) value of 0.5 was considered to be indicative of a lack of diagnostic utility. Kaplan-Meier curves were used to compare patient survival outcomes together with log-rank tests, while independent predictors of patient survival were established through univariate and multivariate Cox proportional hazards methods. P < 0.05 was the threshold of significance.

Results

This study attempted to investigate circRNA_001846 as non–small cell lung cancer biomarker. High levels of serum circRNA_001846 expression were observed in NSCLC. The circRNA_001846 was associated with tumor differentiation, lymph node metastasis, and node metastasis stage. ROC curves analyses revealed that serum circRNA_001846 represents an effective diagnostic biomarker for NSCLC, and the combined evaluation of CEA and serum circRNA_001846 improve an effective diagnostic for NSCLC.

NSCLC patient circRNA profiling

By analyzing a previously published GEO dataset (GSE158695), we identified 114 circRNAs that were differentially expressed in NSCLC patient tumor samples ()). These circRNAs were able to effectively differentiate between NSCLC samples and paracancerous tissues ()), and of these, circRNA_001846 was the most highly upregulated in NSCLC relative to normal tissue levels ()). It was thus selected as a target for subsequent study.

Figure 1. NSCLC patient circRNA profiling results. (a) Differentially expressed circRNAs. (b) Mean-difference plot for the GSE158695 dataset. (c) Differentially expressed genes in the GSE158695 dataset were presented using a Volcano plot. (d) Genes upregulated in the GSE158695 dataset were subjected to clustering analysis. (e) NSCLC samples exhibited circRNA_001846 is upregulation as compared to paracancerous tissue samples

Figure 1. NSCLC patient circRNA profiling results. (a) Differentially expressed circRNAs. (b) Mean-difference plot for the GSE158695 dataset. (c) Differentially expressed genes in the GSE158695 dataset were presented using a Volcano plot. (d) Genes upregulated in the GSE158695 dataset were subjected to clustering analysis. (e) NSCLC samples exhibited circRNA_001846 is upregulation as compared to paracancerous tissue samples

Baseline patient characteristics

To evaluate circRNA_001846 in greater detail, we utilized samples from NSCLC patients and healthy controls (n = 206 each, ). Of these NSCLC patients, 111 and 95 were respectively diagnosed with adenocarcinomas and squamous cell carcinoma. Moreover, 107 patients exhibited lymph node metastases, and 105 and 101 patients, respectively, exhibited good to moderate differentiation and poor differentiation. Based on the TNM staging system, 109 patients had stage I/II disease, while 97 had stage III/IV disease. Patient clinical characteristics are compiled in .

Table 1. NSCLC and control patient clinicopathological characteristics

Table 2. Association between levels of circRNA_001846 expression and clinicopathological characteristics

NSCLC patients exhibit increased serum circRNA_001846 expression

A qPCR approach revealed that NSCLC patients exhibited higher serum circRNA_001846 levels as compared to healthy study participants ().

Figure 2. Serum circRNA_001846 in NSCLC patients. Levels of serum circRNA_001846 were higher in patients with NSCLC relative to healthy individuals. ***p < 0.001

Figure 2. Serum circRNA_001846 in NSCLC patients. Levels of serum circRNA_001846 were higher in patients with NSCLC relative to healthy individuals. ***p < 0.001

Serum circRNA_001846 levels are related to NSCLC patient clinicopathological characteristics

Next, NSCLC patients in our study cohort were stratified into those expressing low and high levels of serum circRNA_001846 (n = 103 each) based upon the median expression of this circRNA among these patients. Subsequent analyses of the clinical characteristics of patients in these two groups revealed high circRNA_001846 levels to be associated with lymph node metastasis (P = 0.002), poor tumor differentiation (P = 0.001), and advanced tumor node metastasis (TNM) stage (P = 0.001) (). These levels were unrelated to patient age (P = 0.158), sex (P = 0.128), histological type (P = 0.097), or smoking status (P = 0.209). Moreover, higher serum levels of this circRNA were linked to significantly reduced progression-free survival (PFS, )) and overall survival (OS, )) among NSCLC patients as compared to those of patients with low serum circRNA_001846 levels. Consistent with these results, serum circRNA_001846 levels were correlated with PFS (P = 0.002, ) and OS (P = 0.005, ) in univariate analyses, as and were also associated with significant reductions in patient PFS (P = 0.001, ) and OS (P = 0.008, ) in multivariate analyses.

Table 3. Univariate and multivariate Cox analysis of factors related to NSCLC patient PFS

Table 4. Univariate and multivariate Cox analysis of factors related to NSCLC patient OS

Figure 3. The relationship between levels of circRNA_001846 in NSCLC patient serum and associated clinicopathological characteristics. PFS (a) and OS (b) were compared between patients expressing low and high levels of circRNA_001846 via log-rank tests

Figure 3. The relationship between levels of circRNA_001846 in NSCLC patient serum and associated clinicopathological characteristics. PFS (a) and OS (b) were compared between patients expressing low and high levels of circRNA_001846 via log-rank tests

Assessment of the diagnostic utility of serum CEA in NSCLC

CEA is perhaps the best-studied serum biomarker used to detect NSCLC at an early stage. The mean serum CEA levels among NSCLC patients and controls in this study were 14.82 ± 17.24 ng/mL and 3.37 ± 2.74 ng/mL (P < 0.0001) ()). ROC curves revealed that these serum CEA levels were able to effectively differentiate between NSCLC patients and controls with an area under the curve (AUC) of 0.736 (95% confidence interval [CI]: 0.685–0.787) at a cutoff level of 5.731 ng/mL ()), yielding respective diagnostic sensitivity and specificity values of 62.6% and 85.9%.

Figure 4. CEA expression and diagnostic utility in NSCLC. (a) Levels of CEA in serum samples were higher in individuals with NSCLC as compared to controls. (b) Serum CEA levels were capable of differentiating between individuals with NSCLC and controls in an ROC curve analysis

Figure 4. CEA expression and diagnostic utility in NSCLC. (a) Levels of CEA in serum samples were higher in individuals with NSCLC as compared to controls. (b) Serum CEA levels were capable of differentiating between individuals with NSCLC and controls in an ROC curve analysis

Serum circRNA_001846 represents an effective diagnostic biomarker for NSCLC

A similar analysis of the diagnostic value of serum circRNA_001846 levels in the evaluation of NSCLC patients was next performed, revealing that these levels could readily differentiate between NSCLC patients and controls with an AUC of 0.872 (95% CI: 0.839–0.906) at an optimal cutoff value of 3.9496 ()), yielding respective sensitivity and specificity values of 78.2% and 81.1%. As such, these results suggested that circRNA_001846 levels may be more sensitive than those of CEA when evaluating patients for NSCLC.

Figure 5. Serum circRNA_001846 diagnostic utility in NSCLC. (a) Serum circRNA_001846 levels enabled effective differentiation between individuals with NSCLC and controls in an ROC curve analysis. (b) An ROC curve was generated for a combination of both CEA and circRNA_001846 as an approach to differentiating NSCLC patients from healthy controls

Figure 5. Serum circRNA_001846 diagnostic utility in NSCLC. (a) Serum circRNA_001846 levels enabled effective differentiation between individuals with NSCLC and controls in an ROC curve analysis. (b) An ROC curve was generated for a combination of both CEA and circRNA_001846 as an approach to differentiating NSCLC patients from healthy controls

To expand upon these findings, we assessed the combined diagnostic value of serum CEA and circRNA_001846 in NSCLC patients. This combined analysis yielded an AUC of 0.925 (95% CI: 0.899–0.950), with respective sensitivity and specificity values of 88.8% and 85.9% at an optimal cutoff of 4.9467, with predictive probabilities above this cutoff level being considered positive for NSCLC.

Discussion

NSCLC is among the most prevalent and deadly tumors in the world, yet patients are often only diagnosed when their disease has already reached an advanced stage at which time they exhibit relatively poor 5-year survival outcomes [Citation17]. The early diagnosis and treatment of these patients is thus vital to ensure their improved treatment outcomes. Currently, approaches used to screen for NSCLC tumors include computed tomography, magnetic resonance imaging, positron emission tomography, and sputum cytology. However, these approaches can be invasive, expensive, and can expose patients to radioactivity [Citation18–20]. In addition, circulating tumor biomarkers may manifest in the periphery at an earlier time point than imaging can detect lung tumors, offering an opportunity to detect otherwise asymptomatic early-stage disease [Citation21]. Evaluating peripheral blood biomarkers thus represents an attractive, noninvasive, inexpensive means of diagnosing NSCLC.

The identification of circRNAs as stable and highly abundant non-coding RNAs with direct oncogenic relevance has been a focus of intensive research interest in recent years as a means of gauging cancer patient prognosis and guiding diagnostic efforts [Citation22–24]. These circRNAs remain stable in patient plasma and serum, and are resistant to the effects of pH shifts or freeze-thawing [Citation25,Citation26]. As such, they are ideal targets for the blood testing-based evaluation of patients in a minimally invasive manner for a range of cancers. For example, Pan et al. found that levels of has_circ_0004771 in patient serum were predictive of colorectal cancer status [Citation27], while serum has_Circ_0141633 has shown promise as a gastric cancer biomarker [Citation28], and Sun et al. found circRNA_001846 to be a promising gastric carcinoma-related biomarker [Citation29]. Functionally in cervical cancer cells, the upregulation of circRNA_001846 led to shifts in miR-1296 activity and consequent cyclin-dependent kinases 2 (CDK2) upregulation, ultimately driving proliferative activity [Citation30]. No prior studies have reported on the diagnostic utility of circRNA_001846 in NSCLC.

Here, we explored the potential diagnostic utility of circulating circRNA_001846 in the context of NSCLC. We found that the upregulation of this circRNA in patient serum was associated with lymph node metastasis, tumor differentiation, and TNM staging, suggesting a potential functional role for circRNA_001846 as a regulator of NSCLC onset or progression. Moreover, higher serum circRNA_001846 levels were linked to worse NSCLC patient OS and PFS, suggesting that it may represent a valuable prognostic biomarker capable of guiding patient management efforts. The diagnostic value of serum circRNA_001846 was further evaluated using ROC curves, revealing that it was able to differentiate between NSCLC patients and healthy controls with excellent sensitivity and specificity that were superior to those of CEA. Moreover, the combined evaluation of CEA levels and serum circRNA_001846 expression achieved 90% specificity in the diagnostic evaluation of NSCLC patients at a cutoff value of 0.678. Such combined biomarker-based analytical approaches may thus represent a potent means of reliably identifying NSCLC patients in the clinical setting.

Conclusion

In summary, circRNA_001846 upregulation was evident in individuals with NSCLC, and such upregulation may be linked to disease incidence. The serum levels of this circRNA may offer value as a diagnostic biomarker for early-stage NSCLC. By combining this biomarker with other standard clinical markers, it may be possible to achieve superior diagnostic reliability. Despite the limitations of this analysis, our results nonetheless highlight directions for future research evaluating the clinical utility of circRNA_001846 when screening for patients with early-stage NSCLC.

Highlights

  1. NSCLC patients exhibit increased serum circRNA_001846

  2. Serum circRNA_001846 are related to NSCLC patient clinicopathological characteristics

  3. Serum circRNA_001846 represents an effective diagnostic biomarker for NSCLC

The combined evaluation of CEA and serum circRNA_001846 improve an effective diagnostic for NSCLC

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was supported by the Study of molecular mechanism of the protective effects of resveratrol on the myocardial injury of pulmonary embolism rats [CYFY2017GLPHX07].

References

  • Sung H, Ferlay J, Siegel RL, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–249.
  • Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
  • How CW, Ong YS, Low SS, et al. How far have we explored fungi to fight cancer? Semin Cancer Biol. 2021. DOI:10.1016/j.semcancer.2021.03.009
  • Ebbesen KK, Hansen TB, Kjems J. Insights into circular RNA biology. RNA Biol. 2017;14(8):1035–1045.
  • Zhou M, Xiao MS, Li Z, et al. New progresses of circular RNA biology: from nuclear export to degradation. RNA Biol. 2020; 18(10):1–9.
  • Zhang Y, Yang L, Chen LL. Life without A tail: new formats of long noncoding RNAs. Int J Biochem Cell Biol. 2014;54:338–349.
  • Wang WB, Ren P, Ren FH, et al. Circ_0000526 blocks the progression of breast cancer by sponging miR-492. Cancer Biother Radiopharm. 2021;36(6):467–476.
  • Wu D, Jia H, Zhang Z, et al. Circ-PRMT5 promotes breast cancer by the miR-509-3p/TCF7L2 axis activating the PI3K/AKT pathway. J Gene Med. 2021;23(2):e3300.
  • Huang B, Cui D, Ren Y, et al. Circ_0006174 promotes colorectal cancer progression by sponging microRNA-142-3p and regulating X-linked inhibitor of apoptosis expression. Int J Biol Markers. 2021;17246008211034178.
  • Li PF, Zhang ZX, Yuan X, et al. Hsa_circ_0001696 modulates cell proliferation and migration in colorectal cancer. Oncol Lett. 2021;21(2):154.
  • Ji C, Hong X, Lan B, et al. Circ_0091581 promotes the progression of hepatocellular carcinoma through targeting miR-591/FOSL2 axis. Dig Dis Sci. 2021;66(9):3074–3085.
  • Yan F, Fan B, Wang J, et al. Circ_0008305-mediated miR-660/BAG5 axis contributes to hepatocellular carcinoma tumorigenesis. Cancer Med. 2021;10(3):833–842.
  • Li H, Shan C, Wang J, et al. CircRNA Hsa_circ_0001017 inhibited gastric cancer progression via acting as a sponge of miR-197. Dig Dis Sci. 2021;66(7):3074–3085.
  • Wang D, Jiang X, Liu Y, et al. Circular RNA circ_HN1 facilitates gastric cancer progression through modulation of the miR-302b-3p/ROCK2 axis. Mol Cell Biochem. 2021;476(1):199–212.
  • Zhu Z, Wu Q, Zhang M, et al. Hsa_circ_0016760 exacerbates the malignant development of nonsmall cell lung cancer by sponging miR1455p/FGF5. Oncol Rep. 2021;45(2):501–512.
  • Ma Q, Huai B, Liu Y, et al. Circular RNA circ_0020123 promotes non-small cell lung cancer progression through miR-384/TRIM44 axis. Cancer Manag Res. 2021;13:75–87.
  • Bade BC, Dela CCS. Lung cancer 2020: epidemiology, etiology, and prevention. Clin Chest Med. 2020;41(1):1–24.
  • Liu S, Lin Z, Rao W, et al. Upregulated expression of serum exosomal hsa_circ_0026611 is associated with lymph node metastasis and poor prognosis of esophageal squamous cell carcinoma. J Cancer. 2021;12(3):918–926.
  • Houston T. Screening for lung cancer. Med Clin North Am. 2020;104(6):1037–1050.
  • Wadowska K, Bil-Lula I, Trembecki L, et al. Genetic markers in lung cancer diagnosis: a review. Int J Mol Sci. 2020; 21(13):4569.
  • Chu GCW, Lazare K, Sullivan F. Serum and blood based biomarkers for lung cancer screening: a systematic review. BMC Cancer. 2018;18(1):181.
  • Meng S, Zhou H, Feng Z, et al. CircRNA: functions and properties of a novel potential biomarker for cancer. Mol Cancer. 2017;16(1):94.
  • Miao L, Yin RX, Zhang QH, et al. A novel circRNA-miRNA-mRNA network identifies circ-YOD1 as a biomarker for coronary artery disease. Sci Rep. 2019;9(1):18314.
  • Wang Y, Mo Y, Gong Z, et al. Circular RNAs in human cancer. Mol Cancer. 2017;16(1):25.
  • Vea A, Llorente-Cortes V, de Gonzalo-Calvo D. Circular RNAs in blood. Adv Exp Med Biol. 2018;1087:119–130.
  • Guo S, Hu C, Zhai X, et al. Circular RNA 0006602 in plasma exosomes: a new potential diagnostic biomarker for hepatocellular carcinoma. Am J Transl Res. 2021;13:6001–6015.
  • 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.
  • Yin G, Zheng S, Shan Y, et al. Circulating circRNA Has_Circ_0141633 serves as a potential biomarker for gastric cancer. Clin Lab. 2020; 66(9).
  • Sun H, Tang W, Rong D, et al. Hsa_circ_0000520, a potential new circular RNA biomarker, is involved in gastric carcinoma. Cancer Biomark. 2018;21(2):299–306.
  • Zheng Q, Zhang J, Zhang T, et al. Hsa_circ_0000520 overexpression increases CDK2 expression via miR-1296 to facilitate cervical cancer cell proliferation. J Transl Med. 2021;19(1):314.