Abstract
Circular RNAs are a class of RNAs with a covalently closed configuration, and several members of them have been reported to be capable of regulating various biological processes and predicting the outcome of disease. Among them, circular RNA circ-ITCH has been identified to be aberrantly expressed and associated with disease progression in diverse cancers. However, the correlation of circ-ITCH expression with clinicopathological features, as well as the prognosis of cancers, remains inconclusive. Therefore, a meta-analysis was performed to investigate the clinical significance of circ-ITCH in cancers by systematically summarizing all eligible literatures. Up to August 31, 2020, relevant articles were searched in PubMed, Web of Science, Cochrane library, Embase, CNKI, and Wanfang databases. Pooled hazard ratios (HRs) and odds ratios (ORs) with corresponding 95% confidence intervals (CIs) were calculated. A total of 1604 patients from 14 studies were included in this meta-analysis. The results showed that cancer patients with low circ-ITCH expression were more susceptible to develop lymph node metastasis (, 95% CI: 1.67-3.02, ), larger tumor size (, 95% CI: 2.01-4.52, ), advanced TNM stage (, 95% CI: 1.92-4.14, ), and poor overall survival (OS) (, 95% CI: 2.07–2.90, , univariate analysis; , 95% CI: 1.82-3.96, , multivariate analysis). Thus, low circ-ITCH expression was significantly associated with aggressive clinicopathological features and unfavorable outcome in various cancers. Therefore, circ-ITCH may serve as a molecular therapy target and a prognostic marker in human cancers.
1. Introduction
Circular RNAs (circRNAs) are a class of newly discovered RNAs with a covalently closed configuration that exist in various organisms [1]. circRNAs are generated from precursor mRNAs (pre-mRNAs) via the back-splicing of exons, introns, or both, to form a covalently closed continuous loop with no 5 caps and 3 poly (A) tails [2]. Initially, circRNAs were considered byproducts of splicing errors. With the advances in the field of high-throughput sequencing, an increasing number of circRNAs and their functions have been characterized [3]. Functionally, circRNAs exert their regulation role through multiple mechanisms. Some circRNAs can work as competing endogenous RNA (ceRNA) to sequester miRNAs and reduce their activity, which in turn positively regulates the expression of miRNA-related target genes [4]. Besides, several circRNAs have been revealed to bind to target proteins by acting as protein decoys [5]. Also, circRNAs can regulate the transcription of targeted genes through interacting with transcription factors [6]. Importantly, recent reports have proved that a subset of circRNAs can be translated [7]. In particular, the essential role of circRNAs in cancer development has been illustrated. Several circRNAs have been reported to be aberrantly expressed and exert oncogenic or tumor-suppressor function in cancers, possibly by acting as sponges for miRNAs. For instance, one of the earliest and best-characterized circRNAs, ciRS-7, contains more than 70 binding sites for miR-7 [2]. Therefore, ciRS-7 was proved to promote cancer progression by acting as miR-7 sponges to upregulate the direct target genes of miR-7 and activate cancer-related pathways. ciRS-7 was found to facilitate the more aggressive phenotype of gastric cancer via suppressing the miR-7-mediated PTEN/PI3K/AKT signaling pathway [8]; in colorectal carcinoma, ciRS-7 was reported to activate the EGFR/RAF1/MAPK pathway via antagonizing miR-7 activity [9]; the ciRS-7/miR-7/NF-κB axis was demonstrated to play a crucial role in accelerating development of esophageal squamous cell carcinoma and lung cancer [10, 11]. circRNAs are also implicated in angiogenesis of cancers. circRNA-MYLK was found to activate vascular endothelial growth factor A (VEGFA) through acting as a sponge of miR-29a, thereby promoting tumor growth, metastasis, and angiogenesis of bladder cancer [12]. Furthermore, circRNAs were implicated in cancer therapy as well. Zhu et al. reported that knockdown of circPVT1, which was found upregulated in osteosarcoma (OS) tissues and chemoresistant cell lines, impairs the resistance to doxorubicin and cisplatin of OS cells a of classical drug resistance-related gene ABCB1 [13]. Furthermore, significant association between the expression of circRNAs and the progression of cancers has been found, indicating the potential of circRNAs to serve as a biomarker for predicting the outcome of cancers [14–17].
Recently, circular RNA Itchy E3 ubiquitin protein ligase (circ-ITCH), a novel circular RNA originated from exons of gene itchy E3 ubiquitin protein ligase (ITCH), located on chromosome 20q11.22, was reported to be lower expressed in several cancers [18]. So far, circ-ITCH has been proved to be implicated in prostate cancer [19, 20], ovarian cancer [21–24], bladder cancer [25], breast cancer [26], lung cancer [27], oral squamous cell carcinoma [28], gastric cancer [29], hepatocellular carcinoma [30], glioma [31], and multiple myeloma [32]. Accumulating evidence has implied that circ-ITCH exert tumor-suppressor function in these cancers by acting as a sponge for oncogenic microRNAs. Although the correlation between circ-ITCH expression and cancer progression has been investigated by these studies above, most individual studies have been limited by inconsistent conclusions or small sample sizes. Thus, we performed this quantitative meta-analysis by systematically evaluating the relationship between circ-ITCH expression and the clinicopathological parameters as well as prognosis of cancers with all eligible articles.
2. Materials and Methods
2.1. Publication Search
Our literature search was performed following the preferred reporting items for systemic reviews and meta-analyses (PRISMA) statement criteria [33]. A comprehensive electronic search was performed in PubMed, Web of Science, Cochrane library, Embase, CNKI, and Wanfang databases updated to August 31, 2020. The keywords during the literature search were “circular RNA ITCH” or “circ-ITCH” or “cir-ITCH” or “circular RNA Itchy E3 ubiquitin protein ligase”. The search strategy in PubMed was “circular RNA ITCH [All Fields] OR circ-ITCH [All Fields] OR cir-ITCH [All Fields] OR circular RNA Itchy E3 ubiquitin protein ligase [All Fields].” The search strategy in Web of Science was “TS = (circular RNA ITCH OR circ-ITCH OR cir-ITCH OR circular RNA Itchy E3 ubiquitin protein ligase)”. In addition, the citation lists of retrieved articles were screened manually for potential eligible studies.
2.2. Inclusion and Exclusion Criteria
Studies were considered eligible if they fulfilled the inclusion criteria as follows: (1) articles investigated the correlation of circ-ITCH expression with cancer progression and/or clinicopathological factors, (2) the expression of circ-ITCH in cancerous tissues was measured, (3) cancer patients were divided into high/low groups according to the circ-ITCH expression, and (4) related clinicopathological parameters and/or prognostic results were described. Exclusion criteria of the present meta-analysis are the following: (1) duplicated publications; (2) reviews, letters, comments, and conference articles; (3) articles irrelevant to the present study; or (4) studies without available data.
2.3. Data Extract
Two investigators (Xiao-Dong Sun and Chen Huan) performed the data extraction from the eligible studies independently. Discrepancies were resolved by discussion with a third investigator (Da-Wei Sun) to reach a consensus. The following items were collected from each included study: first author, year of publication, origin of patients, cancer type, number of patients, detecting method of circ-ITCH expression, cutoff value for grouping, number of patients with lymph node metastasis (LNM), larger tumor size and advanced TNM stage in each group, follow-up period, survival analytical method (multivariate or univariate), and hazard ratio (HR) with 95% confidence interval (CI) for overall survival (OS). When the prognosis was plotted as a Kaplan-Meier curve, the software Engauge Digitizer version 4.1 (http://digieizer.sourceforge.net/) was applied to digitize the data, and HR with 95% CI was calculated as described [34].
2.4. Quality Assessment
The methodological quality of the included studies was assessed with Newcastle-Ottawa Scale (NOS) criteria, which is scored based on subject selection, comparability of subject, and clinical outcome [35]. The final scores of NOS ranged from 0 to 9, and studies with were considered to be of high quality.
2.5. Data Analysis
Meta-analyses were conducted using Stata SE12.0 (Stata Corporation, College Station, Texas). Pooled indicated that low circ-ITCH correlated with poor prognosis, and indicated low circ-ITCH correlated with poor clinicopathological outcomes. The heterogeneity among the included studies was assessed through -based test and statistics. When the value > 50% and/or , indicating that the heterogeneity was significant, a random-effects model was applied; otherwise, the fixed-effects model was adopted. To assess the publication bias of included studies, Begg’s and Egger’s tests were conducted to estimate the potential publication bias quantificationally, where a value < 0.05 illustrates significant publication bias [36, 37]. In addition, funnel plots were used to present the distribution of included studies’ results when more than 8 studies were included, which represented the publication bias based on visual inspection. Sensitivity analysis was also performed to evaluate the effect of each individual study on the overall effect of meta-analysis results. All tests were two-sided; values < 0.05 were considered statistically significant.
3. Results
3.1. Literature Information
The flow diagram for literature screening and selection was shown in Figure 1. A total of 110 records were retrieved by searching the databases, and 96 articles were excluded according to the inclusion and exclusion criteria. Finally, 14 articles comprising 1604 patients were identified as eligible and included in the present meta-analysis.

3.2. Study Characteristics
The main characteristics of eligible studies are summarized in Table 1. These 14 enrolled articles were published between 2017 and 2020 with sample sizes ranging from 20 to 288. Most of the populations were from China and divided into a high or low group based on the median value of circ-ITCH expression. The expression of circ-ITCH was detected with the method of quantitative reverse transcription-polymerase chain reaction (qRT-PCR) in all 14 populations. According to the NOS criteria, all of the included studies got , indicating their high methodological quality.
3.3. Association between circ-ITCH and Clinicopathological Parameters
As shown in Figure 2, pooled meta-analysis was performed to estimate the relationship between circ-ITCH expression and clinicopathological features of cancers. Since there was no significant heterogeneity among these studies, a fixed-effects model was exploited. The pooled OR with 95% CI indicated that cancer patients with low circ-ITCH expression were more susceptible to develop LNM (, 95% CI: 1.67-3.02, ) and advanced TNM stage (, 95% CI: 1.92-4.14, ), as well as larger tumor size (, 95% CI: 2.01-4.52, ), suggesting that a low circ-ITCH level may serve as an indicator of aggressive clinicopathological features for cancer patients.

(a)

(b)

(c)
3.4. Association between circ-ITCH and OS
On one hand, 13 studies comprising a total number of 1574 patients investigated the association between circ-ITCH expression and OS through univariate analysis. The fixed-effects model was used to assess the pooled HR and its 95% CI since no heterogeneity was found among these studies (, ). We found that low circ-ITCH expression was significantly associated with poor OS (, 95% CI: 2.07–2.90, ) (Figure 3(a)). Besides, subgroup meta-analysis was also conducted. The results showed that low circ-ITCH expression was a significant prognostic indicator of poor OS for patients with different types of cancers: prostate cancer (, 95% CI: 1.89-4.39, ), ovarian cancer (, 95% CI: 1.82-4.02, ), and other types of cancers (, 95% CI: 1.86-2.82, ). Meanwhile, the significant association between low ITCH expression and unfavorable OS was stable despite the variation of sample size and different cutoff values (Table 2).

(a)

(b)
On the other hand, 4 studies with a total number of 623 patients investigated the association between circ-ITCH expression and OS through multivariate analysis. Since there was no heterogeneity among these studies (, ), the fix-effects model was used to assess the pooled HR and its 95% CI. We found that low circ-ITCH expression was also significantly associated with poor OS (, 95% CI: 1.82-3.96, ) (Figure 3(b)).
3.5. Sensitive Analysis
To assess the robustness of our results, sensitivity analysis was conducted by omitting each individual included study. As illustrated in Figure 4, removing any of the enrolled studies did not change the overall meta-analysis effect of circ-ITCH on the pooled ORs and HRs, indicating that our findings were relatively stable.

(a)

(b)

(c)

(d)

(e)
3.6. Publication Bias
In this meta-analysis, both Begg’s and Egger’s value tests were used to assess the potential publication bias. No publication bias was found in most analyses, including the studies with LNM (, 0.761), TNM stage (, 0.890), tumor size (, 0.727), and OS (, 0.101, multivariate analysis). Publication bias was found in the studies with OS (, 0.011, univariate analysis). Besides, the funnel plots of OS from univariate analysis (Figure 5) were largely symmetrical. Therefore, we speculate that most of our meta-analysis results are reliable.

4. Discussion
Recently, as increasing studies have demonstrated the participation of circRNAs in carcinogenesis, the potential of circRNA to predict cancer progression has been suggested due to the correlation between their expression and clinicopathological characteristics as well as the outcome of cancers. circ-ITCH was a newly identified circRNA; many studies have investigated the association between circ-ITCH expression and prognosis in cancers. However, the sample sizes of most studies are small. Besides, there is no consensus about the prognostic value of circ-ITCH expression in cancers. Here, we performed this meta-analysis to investigate the clinical and prognostic value of circ-ITCH in cancers.
We included 14 studies with a total of 1604 cancer patients in this meta-analysis. The pooled ORs with their 95% CIs showed that low circ-ITCH expression was significantly associated with larger tumor size, increased LNM, and advanced TNM stage, indicating that low circ-ITCH expression was an indicator of aggressive clinicopathological parameters. Moreover, the pooled HRs with their 95% CIs showed that low circ-ITCH expression was also significantly correlated with poor OS, implying that low circ-ITCH expression may serve as an indicator of unfavorable prognosis of cancers. Since circRNAs lack free 3 or 5 tails, they are more resistant to exonuclease RNase R-induced degradation and confer longer half-lives than that of linear mRNAs [2, 38, 39]. Most circRNAs exhibit a half-life longer than 48 h, while an average value of mRNAs is 10 h [40]. Meanwhile, circRNAs are illustrated to be more abundant than their linear isoforms in human cells in a cell type-specific manner [41]; lots of studies have demonstrated that the presence and abundance of circRNAs in different types of cancer cells are distinctive as well [42]. Thus, these remarkable characteristics of circRNAs above could make them serve as detectable biomarkers in diseases [43, 44]. Taken together, circ-ITCH could serve as a biomarker for predicting the progression and outcome of cancers.
Functionally, accumulating evidence has implied a tumor-suppressor role circ-ITCH in diverse cancers. For instance, by acting as a sponge for oncogenic miR-214 and miR-17, circ-ITCH significantly enhances expression of its ITCH linear isoform via competitive interacting with microRNAs, thereby inactivating Wnt/beta-catenin signaling in various cancers [26, 31, 45–47]. Meanwhile, circ-ITCH was also reported to act as ceRNAs of other microRNAs, like microRNA-93-5p and miR-145, to execute its tumor suppressive activity in cervical cancer and ovarian carcinoma, respectively [21, 48]. Additionally, circ-ITCH was further found to inhibit tumorigenesis through other mechanisms rather than acting as a microRNA sponge. In melanoma, circ-ITCH suppresses cancer cell proliferation via impairing glucose uptake of cancer cells [49]. Meanwhile, it has been demonstrated that ectopic expression of circ-ITCH was capable of inhibiting cancer growth in vivo [19, 25], hinting that circ-ITCH might be a potential approach for cancer treatment. The development of RNA-related therapeutics, especially the investigation of strategies to manipulate circRNA levels, may facilitate the circRNA-based therapeutic strategies in the near future [50].
To the best of our knowledge, this is the first meta-analysis to investigate the clinical significance of circ-ITCH in cancer patients. Nevertheless, some limitations of this meta-analysis should be declared. The primary concern is that most of the included studies were conducted based on the population from China; thus, the results should be substantiated by additional studies in the worldwide population. Secondly, publication bias was observed in the studies with OS via univariate analysis. For some studies, the OS results were not available, which may contribute to the publication bias. Based on these limitations above, prospective and well-designed clinical studies with large scale, as well as studies based on other populations beyond China, are still warranted to investigate the role of circ-ITCH in cancer patients.
5. Conclusion
In summary, our meta-analysis revealed that low circ-ITCH expression was significantly associated with larger tumor sizes, advanced TNM stage, increased LNM, and poor survival rate in cancers. Therefore, circ-ITCH may serve as a prognostic biomarker and a promising molecular therapy target in cancers.
Data Availability
The data supporting this meta-analysis are from previously reported studies that have been cited. The processed data are available from the corresponding author upon request.
Disclosure
This manuscript has been presented as a preprint (https://www.researchsquare.com/article/rs-29660/v1).
Conflicts of Interest
All the authors declared no conflicts of interests in this work.
Authors’ Contributions
Xiao-Dong Sun and Chen Huan have contributed equally to this work.
Acknowledgments
This work was supported by research funding from the Department of Finance of Jilin Province (No. 2018SCZWSZX-044).