Next Article in Journal
Comprehensive Molecular Analysis Identified an SRSF Family-Based Score for Prognosis and Therapy Efficiency Prediction in Hepatocellular Carcinoma
Next Article in Special Issue
Plasma-Based microRNA Expression Analysis in Advanced Stage NSCLC Patients Treated with Nivolumab
Previous Article in Journal
The Impact of Race–Ethnicity and Diagnosis of Alzheimer’s Disease and Related Dementias on Mammography Use
Previous Article in Special Issue
A miRNA-Based Prognostic Model to Trace Thyroid Cancer Recurrence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review and In Silico Analysis of Tissue and Exosomal Circular RNAs: Opportunities and Challenges in Thyroid Cancer

1
Division of Endocrine and Oncologic Surgery, Department of Surgery, School of Medicine, Tulane University, New Orleans, LA 70112, USA
2
Genetics Unit, Department of Histology and Cell Biology, Suez Canal University, Ismailia 41522, Egypt
3
Department of Medical Biochemistry and Molecular Biology, Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
4
Department of Biochemistry, Faculty of Medicine, Northern Border University, Arar 91431, Saudi Arabia
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(19), 4728; https://doi.org/10.3390/cancers14194728
Submission received: 16 July 2022 / Revised: 22 September 2022 / Accepted: 26 September 2022 / Published: 28 September 2022

Abstract

:

Simple Summary

Thyroid cancer is the most common endocrine neoplasm. Recently, knowledge of the molecular genetic changes of thyroid cancer has dramatically improved. Understanding the roles of these molecular changes in thyroid cancer tumorigenesis and progression is essential in developing a successful treatment strategy and improving disease outcomes. As a family of non-coding RNAs, circular RNAs (circRNAs) have been involved in several aspects of the physiological and pathological processes of the cells. The roles of circRNAs in cancer development and progress are evident. In the current review, we aimed to explore the clinical potential of circRNAs as potential diagnostic, prognostic, and therapeutic targets in thyroid cancer. Furthermore, screening the genome-wide circRNAs and performing functional enrichment analyses for all associated dysregulated circRNAs in thyroid cancer have been done. Given the unique advantages circRNAs have, such as superior stability, higher abundance, and presence in different body fluids, this family of non-coding RNAs could be promising diagnostic and prognostic biomarkers and potential therapeutic targets for thyroid cancer.

Abstract

Thyroid cancer (TC) is the most common endocrine tumor. The genetic and epigenetic molecular alterations of TC have become more evident in recent years. However, a deeper understanding of the roles these molecular changes play in TC tumorigenesis and progression is essential in developing a successful treatment strategy and improving patients’ prognoses. Circular RNAs (circRNAs), a family of non-coding RNAs, have been implicated in several aspects of carcinogenesis in multiple cancers, including TC. In the current review, we aimed to explore the clinical potential of circRNAs as putative diagnostic, prognostic, and therapeutic targets in TC. The current analyses, including genome-wide circRNA screening and functional enrichment for all deregulated circRNA expression signatures, show that circRNAs display atypical contributions, such as sponging for microRNAs, regulating transcription and translation processes, and decoying for proteins. Given their exceptional clinical advantages, such as higher stability, wider abundance, and occurrence in several body fluids, circRNAs are promising prognostic and theranostic biomarkers for TC.

1. Introduction

Thyroid cancer (TC) is the most prevalent endocrine-related neoplasm, and its incidence continues to increase [1]. Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer, accounting for nearly 85% of cases [2]. Although the five-year overall survival rate of patients with PTC is nearly 95%, the cancer can metastasize into lymph nodes and distant organs, resulting in poor prognosis and recurrence in nearly 30% of patients [3]. It is critical that we unravel the mechanisms underlying cancer initiation and progression in order to design novel diagnostics and therapeutics.
Non-coding RNAs (ncRNAs) are non-protein coding transcripts that can be classified into two main types: small ncRNAs (18–200 nucleotides) and long ncRNAs (lncRNAs, >200 nucleotides) [4,5]. The microRNA family (miRNAs; 22–25 nucleotides) belongs to the former class, binding to miRNA response elements (MREs) in the 3′-untranslated region (3′-UTR) of the target messenger RNA (mRNA) to induce mRNA degradation and/or suppress its translation into protein [6,7].
Circular RNAs (circRNAs) have emerged as a large family of endogenous ncRNA transcripts widely expressed in different cellular contexts and can act as post-transcriptional regulators. They are generated by head-to-tail splicing that has neither 5′–3′ polarities nor a polyadenylated tail, thus exhibiting a covalently closed loop structure [8,9]. Most circRNAs are derived from the gene’s exon region, but a small percentage are created by intron cleavage [10]. They were first identified in humans nearly 40 years ago using electron microscopy [8]. Initially, circRNAs were misinterpreted as splicing errors and transcriptional artifacts [11,12,13]. Given the recent advances in high-throughput sequencing and bioinformatics analyses, thousands of circRNAs have been recognized as abundant stable ncRNAs in various species and cell lines [14].
Several studies have emphasized the essential roles of circRNAs in both normal cellular processes and pathological conditions, as depicted in diverse types of cancers [15,16,17,18,19]. Various circRNAs have been recognized as cancer-specific non-coding RNAs in tumorigenesis [20], such as circMTO1 in hepatocellular carcinoma, CCDC66 in colon cancer, hsa_circ_0067934 in esophageal squamous carcinoma, circAGFG1 in breast cancer, and cTFRC in bladder carcinoma [14,21,22,23,24,25,26]. The aberrant expression of circRNAs may impact their regulatory roles in gene expression, which contributes to carcinogenesis by disrupting several signaling pathways, including the Wnt/β-catenin pathway [27], phosphatidylinositol-3-kinase (PI3K)/AKT signaling transduction pathway [28], vascular endothelial growth factor A (VEGFA)/VEGF receptor-2 pathway, and Ras/extracellular signal-regulated kinase (ERK) cascade signaling pathway [29]. These pathways control multiple cellular processes that drive cancer development, including cell growth, cell cycle, migration, angiogenesis, survival, and apoptosis, among others.
A growing body of literature supports the association of circRNAs with PTC carcinogenesis [27,30,31,32]; however, the exact role of circRNA in thyroid tumorigenesis and progression has not been fully clarified. Overall, the current review aims to explore the clinical utility of circRNAs as putative prognostic and theranostic targets in TC by revising related information from the published literature on several data sources and running in silico analysis to provide a suggested list of differentially expressed circRNAs in thyroid cancer.
Throughout this review, a brief presentation of the circRNA nomenclature, biogenesis, biological activities, and the online related databases/resources has been carried out, followed by a review of the diagnostic and prognostic utility of circRNAs in patients with TC, how the circRNAs mediate tumor progression in vivo and in vitro, the tumor-cell-derived exosomal circRNAs, their implication in treatment resistance, functional enrichment for their deregulation, screening the genome-wide circular RNAs and identifying their microarray transcriptomic signature in PTC, and lastly, a discussion of the current challenges in clinical practice and future perspectives have been provided.

1.1. CircRNA Nomenclature

There are no specified HUGO standards for gene nomenclature for circRNAs, but recent guidelines recommend using “circ(gene symbol)-n”, where the gene symbol represents the unspliced “host” gene and n is a five-digit number listing the order of discovery for circRNAs [33]. For example, the 15 circRNAs that are derived from the EGFR gene are named “hsa_circEGFR_001” to “hsa_circEGFR_015” [20] with “hsa (homo sapiens)” indicating human-derived transcripts.

1.2. CircRNA Biogenesis, Export, and Turnover

CircRNA and linear RNA are usually are co-expressed [20]; however, nascent circRNAs are generally processed at a later stage than linear ones [34]. The size of a circRNA molecule can range from <100 to >4000 nucleotides [35]. They can be formed from spliced introns or ≥1 exons [36]. Based on their origin from different regions of the genome, circRNAs can be classified into three types: (1) exonic circRNA (EcircRNA), representing the majority of circRNAs (more than 80%), is derived from direct back-splicing or exon skipping; (2) circular intronic circRNAs (ciRNA) that consist of lariat introns that are generated due to the presence of specific elements near the branchpoint site and the 5’ splice site, these types of introns are much more likely to contain complementary inverted Alu repeats, which is the most common transposon in the genome [35]; and (3) exon–intron circRNA (EIciRNA) that contains both exons and introns [15,35,37] and is produced from exons, introns, or both [36,38] (Figure 1). Two other types of circRNAs, intergenic and antisense, are uncommon and remain poorly understood [39].
Several types of circRNAs are created from precursor mRNAs (pre-mRNAs) by distinct processes such as exon skipping, intron pairing, and RNA-binding proteins [15,39,40]. Given the endless nature of circRNAs, this adds exceptional stability and resistance to enzymatic digestion by exonucleases [20,35] and allows some circRNAs to accumulate to levels that may surpass the linear forms [41]. CircRNA half-lives can reach up to 48 h, relative to about 10 h for linear RNAs [35]. Rates of biogenesis, cellular export, and turnover, together regulate the steady-state intracellular levels of circRNAs [42]. EcircRNAs are predominantly found in the cytoplasm, whereas ciRNAs and EIciRNAs are more abundant in the nucleus [38]. How circRNAs are exported from the nucleus is not entirely clear, but the N6-methyladenosine (m6A) modification is often found in circRNAs and has been shown to affect their nuclear export [43]. CircRNAs were also found in secreted extracellular vesicles and were more abundant in exosomes than cancer cell lines [44,45,46].
CircRNAs often have tissue-restricted and cancer-specific expression profiles, which are regulated by specific trans-acting factors and cis-acting elements. They can accumulate to levels exceeding messenger RNAs and show cell/development-specific expression. In cancer, multiple factors can affect the balance between canonical RNA splicing for intron removal between adjacent exons, and back-splicing, leading to deregulation of circRNA expression [42]. The elucidation of deregulated circRNAs and identification of their biological roles are ongoing processes in cancer research.

1.3. Biological Activity of circRNAs

There is growing evidence that circRNAs can influence a range of biological pathways through diverse mechanisms of action, including serving as microRNA sponges [47,48,49], protein decoys [16,50], regulators of transcriptional activity [51] that can mediate the parental expression genes [52], regulators of alternative splicing [16], and modulators of mRNA stability [53]. In addition, circRNAs can act as templates for protein translation and sequester proteins from their indigenous subcellular localization [16,50,51,54] (Figure 2). CircRNAs are also known to function as mediators of RNA–protein interactions [55] and scaffolds for protein complex assembly [9].
Nuclear circRNAs have been implicated in regulation of transcription via protein recruitment or as structural scaffolds. Cytoplasmic circRNAs function in several ways, including as miRNA sponges and/or RNA-binding proteins and as scaffolds for proteins or enzyme–substrate complex assembly, facilitating posttranslational modifications [20]. CircRNAs have also recently been identified as essential components of the microvesicle cargo and the exosome, where they can be enriched related to their cell production level [56].

1.3.1. Epigenetic Regulators

A handful of circRNAs can indirectly affect transcription by modulating histone modification via miRNA binding. For example, circ_100284 is involved in malignant transformation via microRNA-217 regulation of “enhancer of zeste homolog 2 (EZH2)” [57], whereas circ_0071589 can intensify the expression of EZH2 methyltransferase through sponging miR-600 [52].
Furthermore, circRNAs can mediate epigenetic modification through DNA methylation. For example, circ-FECR1 can directly bind to the promoter of DNMT1 methyltransferase, leading to transcription inhibition [58].

1.3.2. Transcription Regulators

Several circRNAs remain at their nuclear sites of synthesis and seem to function as gene expression regulators for their parent coding genes [15]. Exonic circRNAs, which make up the majority of the explored circRNAs, are mainly located in the cytoplasm; therefore, they are unlikely to play a role in regulating the transcription process [59]. On the other hand, two EIciRNAs, circ-EIF3J and circ-PAIP2, were found to bind with U1 snRNP, forming a complex with RNA Pol II in the promoter sequence of parental genes, boosting their transcription through a positive feedback loop [52]. In addition, other circular introns have the potential to bind to the elongating RNA Pol II and enhance the transcription of their host gene [60]. Studies have also suggested that circRNAs can act as splicing regulatory factor processors, since the back-splicing mechanism during which they are formed competes for precursor mRNA with the canonical mechanism of linear RNAs [16,61].

1.3.3. Regulation of Transcription Factors

Emerging evidence shows that circRNAs regulate the activity of the transcriptional activators and/or repressors activity that control the expression of key players in tumor biology. The circular form of TP53 (circ-TP53) was found to inhibit the linear p53 transcription factor [30]. The basic helix-loop-helix transcription factors, TCF3 and TCF4, were downregulated by circ_102171 in TC [27]. The same circular RNA can suppress the DNA-binding transcription factor LEF1 [27], whereas circ-ABCB10 represses the expression level of the C2H2 zinc finger transcription factor KLF6 [47]. On the other hand, circ-FAT1(e2) causes overexpression of ZEB1, a C2H2 zinc finger transcription factor, via its sponging of miR-873 [62]. Circ_0001018 and circ_0005273 are involved in mediating the expression of the HMG box transcription factors SOX2 and SOX4 [63,64]. Similarly, circ-NRIP1 can directly interact with miR-195-5p to regulate the levels of the DNA-binding transcription factor STAT1 [65], and circ-PRKCI is involved in regulating E2F3 by directly binding to miR-335 [66]. Moreover, circ_0039411 promotes the expression of the MTA1 transcription factor via sponging miR-1179 and miR-1205 [49]. Given their role as transcriptional regulators, circRNAs are related to the expression of target miRNAs and proteins associated with TC cell proliferation, invasion, and metastasis, contributing to TC development (Table 1).
Table 1. Transcriptional regulation by circRNAs in thyroid cancer.
Table 1. Transcriptional regulation by circRNAs in thyroid cancer.
Transcription FactorsType of TFcircRNAParental GenemiRNA SpongeReference
Upregulated TF
↑MTA1Homeodomain TFcirc_0039411 ↓miR-1179, miR-1205[49]
↑ZEB1C2H2 zinc finger TFcirc-FAT1(e2)FATI↓miR-873[62]
↑SOX2HMG box TFcirc_0005273 ↓miR-1138[64]
↑SOX4HMG box TFcirc_0001018 ↓ miR-338-3p[63]
↑STAT1DNA-binding TFcirc-NRIP1NRIP1↓miR-195-5p[65]
↑E2F3General TFcirc-PRKCIPRKCI↓miR-335[66]
Downregulated TF
↓TCF3, TCF4, LEF1Basic helix-loop-helix TFcirc_102171 [27]
↓TP53P53-like TFcirc-TP53TP53↓miR-1233-3p[30]
↓KLF6C2H2 zinc finger TFcirc-ABCB10ABCB10 [47]

1.3.4. Acting as Molecular Sponges

RNA transcripts such as mRNAs, lncRNAs, and pseudogenes engage in mutual regulations with each other via competitively shared miRNAs [14]. The latter can downregulate gene expression, destabilize target RNAs, or inhibit translation [67]. In this way, depending on the target genes, miRNAs can act as oncogenes or tumor-suppressor genes. It has been speculated that all classes of RNAs may compete for miRNAs. This competing endogenous RNA (ceRNA) activity forms wide-scale trans-regulatory crosstalk across the transcriptome and plays essential roles in cancer [68].
It has been shown that circRNAs are an essential subtype of ceRNAs. They can act as miRNA molecular sponges via binding to miRNA response elements (MREs), attenuating the effect of miRNAs on target gene expression [50]. CircRNAs can also contain multiple MREs that regulate gene expression at the transcriptional and/or post-transcriptional level [50]. Because circular RNA sponges are characterized by high expression levels, stability, and many miRNA binding sites in their structure, they are likely to be more effective sponges than linear transcripts [69].
The “circRNA-miRNA-target gene” axis could play an essential role in regulation of several biological features of cancers. It has been found that numerous ceRNA networks are dysregulated during TC tumorigenesis, invasion/metastasis, epithelial–mesenchymal transition (EMT), and therapeutic resistance [70] (Figure 3). Many studies have suggested that circRNAs function as miRNA sponges, with subsequent loss of miRNA suppressive function on target genes. For example, circ-ITCH inhibits TC progression by sponging miR-22-3p and upregulating CBL level [71]; circ_0025033 accelerates TC proliferation through sponging miR-1233, and miR-1304 [48]; circ_0005273 acts as an oncogene through targeting miR-1138 [64]; circITGA7 acts as miR-198 ceRNA to regulate FGFR1 expression [72]; and overexpression of circular RNA circ-NEK6 abrogates the tumor suppressive effect of miR-370-3p on TC via the WNT signaling pathway [73]. Since deregulation of the ceRNAs networks is implicated in cancer initiation and progression, it can provide an opportunity for targeted therapeutics.

1.3.5. Acting as Regulators of Translation into Protein

Due to their lack of polarity, IRES, and 5′ cap, circRNAs were first thought to be incapable of forming proteins; however, evidence shows that some circRNAs can act as templates for protein translation due to harboring initiating sequences [61,74,75]. For example, “circ-SHPRH” can encode the SHPRH-146 amino acid protein, which acts as a tumor suppressor [76].

1.3.6. Acting as Protein Decoys

Beyond regulating microRNAs, circular RNAs may bind and sequester RNA-binding proteins (RBP) or even base pair with RNAs besides microRNAs, forming large RNA–protein complexes [77]. For example, ciRS-7 is a novel circRNA that regulates miRNA, mRNA, and RBP functions to influence specific physiological processes [77]. As an oncogene, ciRS-7 suppresses miR-7 and promotes tumor growth in several cancer types. Furthermore, circRNA ciRS-7 has been implicated in promoting EMT in PTC cells [78], and via modulating the miR-7/EGFR axis, ciRS-7 enhances PTC cell proliferation and migration cells [78]. Circ-PABPN1 binds to ELAV-like protein 1, one of the RBPs, preventing its binding to its cognate PABPN1 mRNA with subsequent PABPN1 translation suppression [79]. Circ_102171 activates the Wnt signaling pathway by binding to CTNNBIPI [27].

1.4. Online Databases and Resources

Multiple websites and freely accessible circRNA-related online databases have recently been established. For example, “CircAtlas v2.0, CIRCpedia v2, and TCSD” provide the expression levels of circRNA across different human tissues/cell lines. Also, “CircRiC” and “MiOncoCirc v2.0” explore circRNAs in cancer cell lines and human cancer samples, respectively. “CircAtlas v2.0, circBase, and circRNADisease”, in addition, present circRNA–miRNA interactions. “CircInteractome, CircRiC, CSCD, TSCD, CircFunBase, and Circ2Disease” predict “circRNA–RBP” interactions. Meanwhile, “CIRCpedia, circFunbase, Circad, NPInter (v4), and CircAtlas (v2)” databases cover circRNA-related information across different species. The complete list of these databases [80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107] with their description and weblinks were detailed in our previous work [108].

2. Evidence Acquisition

This study provides a review of 64 related published original articles in PubMed and Google Scholar databases through 23 June 2022, and by applying keywords such as “(circular OR circRNA) AND (thyroid) AND (cancer OR carcinoma OR tumor)” to explore the role of circular RNAs in thyroid tumorigenesis and progression in humans, and navigate the in vitro and in vivo studies to assess their diagnostic/prognostic performance as tissue and liquid biopsy biomarkers. References of all relevant review articles have also been cross-checked against the search.

3. Results

This section of the current review includes six subsections; the diagnostic and prognostic utility of circRNAs in TC patients, circRNAs mediation for tumor progression in vivo and in vitro, tumor-cell-derived exosomal circRNAs, circRNAs and treatment resistance, functional enrichment analysis for the deregulated markers, and genome-wide circular RNA screening.

3.1. Diagnostic and Prognostic Utility of circRNAs in TC Patients

In the last decade, circRNAs have received much attention as key players in cancer development and progression, supporting their utility as biomarkers for cancer diagnosis or monitoring [42]. Because of their closed-loop structure, circRNAs have a high level of stability and resistance to degradation [20,35] and are thus likely to be promising biomarkers for PTC. Many circRNAs show increased expression levels in TC tissues and cell lines compared with adjacent non-cancer tissues, normal tissues, and normal cell lines and promote tumor development.
The “circRNA pleckstrin and Sec7 domain containing 3 (circ_PSD3)” were aberrantly upregulated in PTC tumor tissues compared with adjacent normal tissues [109]. Overexpression of hsa_circ_0011290 in PTC compared with matched para-cancerous thyroid tissues was associated with advanced stage and poor prognosis of PTC [32]. Yao et al. demonstrated that overexpression of circ_0058124 in PTC tissues was associated with malignant features and poor outcomes in PTC patients in terms of advanced tumor-node-metastasis (TNM) stage, tumor size, LN metastasis, extrathyroidal extension, and distant metastasis [110]. Upregulation of circ_0005273 was reported in PTC tissues and was related to poor survival [64]. Overexpression of circ_BACH2 was reported in PTC tissues compared with adjacent tissues and was closely linked to larger tumor size, advanced TNM stage, lymph node metastasis, and shorter survival times [3]. Upregulated circ_0008274 was associated with poor PTC prognosis. Overexpression was associated with advanced TNM stage, tumor infiltration, and lymph node metastasis [111,112].
Similarly, circ_0067934 overexpression in TC tissues was associated with larger tumor size, higher AJCC stage, lymph node metastasis, and shorter survival times [113]. Wang et al. reported a lower level of circRNA itchy E3 ubiquitin-protein ligase (circ-ITCH) expressed in PTC tissues than in normal adjacent tissues [71]. In addition, circ_0137287 was downregulated in PTC tissues and related to aggressive clinical pathologic features in PTC patients, such as larger tumor size, advanced T stage, lymph node metastasis, and extrathyroidal extension [114] (Table 2).
Table 2. Role of circRNAs as diagnostic and prognostic biomarkers in PTC patients.
Table 2. Role of circRNAs as diagnostic and prognostic biomarkers in PTC patients.
circRNAsTNMETEClinical StagePoor SurvivalRef.
Upregulated biomarkers
circ_PSD3 [109]
circ-UMAD1 [115]
circ_0058124 [110]
circ-NRIP1 [65]
circ_0005273 [64]
circ_0067934 [113]
circ_BACH2 [3]
circ_0008274 [111,112]
circ_0004458 [116]
circ_0062389 [117]
circ_0011290 [32]
circ-NRIP1 [65]
circ_0001666 [118]
Downregulated biomarkers
circ_ITCH [71]
circ_0137287 [114]

3.2. Circular RNAs Mediate Tumor Progression In Vivo and In Vitro

Emerging evidence suggests that some circRNAs play a significant role in tumor progression via multiple mechanisms, such as cell proliferation, differentiation, autophagy [119,120], EMT [56,78,121,122], immune modulation [51], apoptosis [114], vascularization [123], invasion, and metastasis [120,124]. The circRNAs–miRNA–mRNA axis can deregulate various cancer-related signaling pathways [29,125] (Figure 4).

3.2.1. Tumor-Suppression-Related circRNAs

Some circRNAs were found to exert an anti-tumor function in PTC. In a microarray experiment, downregulated circRNA_100395 was observed in six PTC tumors compared with matched contralateral normal tissues. Circ_100395 showed interactive potential with two cancer-related miRNAs (miR-141-3p and miR-200a-3p), suggesting a role for the “hsa_circRNA_100395/miR-141-3p/miR-200a-3p” axis in PTC pathogenesis [130]. Moreover, downregulation of circ-ITCH was reported in “K1, IHH-4, and TCP-1” PTC cell lines [71]. Circ-ITCH promoted apoptosis in vitro and inhibited cancer cell proliferation in vivo, based on gain-of-functional assays. Circ-ITCH upregulates the expression of “casitas B-lineage lymphoma (CBL)”, an E3 ligase of nuclear β-catenin, by sponging miR-22-3p. Elevated levels of CBL suppressed the Wnt/β-catenin pathway activity with subsequent inhibition of PTC progression [71]. Another study analyzed microarray data from 18 thyroid samples and identified downregulated circRNA_100395 as a promising biomarker for papillary thyroid carcinoma through modulation of “miR-141-3p” and “miR-200a-3p” [130].

3.2.2. Tumorigenesis-Related circRNAs

Other upregulated circRNAs are known to promote tumor growth (Table 3). CircTP53 is upregulated in TC cell lines, and this overexpression is associated with increased cancer cell proliferation and viability of TPC-1 cells. CircTP53 sponge miR-1233-3p thus increased “mouse double minute 2 (MDM2)” expression and downregulates the protein level of p53. Knockdown of circTP53 inhibited the expression of MDM2 and increased the protein level of p53 [30]. Jin and colleagues reported that circ_0004458 was overexpressed in PTC tissues and was associated with tumor size, invasion, lymphatic infiltration, and distal metastasis. Furthermore, expression levels were upregulated in different TC cells (BCPAP, TPC-1, K1, and IHH4) compared with NORI cells. In contrast, circ_0004458 silenced suppressed cell growth and enhanced apoptosis and cell cycle arrest in PTC cell lines in vitro, and hampered tumor growth in nude mice by inhibiting miR-885-5p and activating “Rac family small GTPase 1 (RAC1)” [116].
Han et al. reported circ_ABCB10 upregulation in PTC tissues relative to adjacent tissues. Circ_ABCB10 overexpression induced the proliferation/invasion of TC cells by targeting KLF6. The silencing of this circRNA promotes apoptosis both in vitro and in vivo and, therefore, can be considered for future targeted therapy in PTC [47]. Similarly, circ_0025033 was elevated in PTC. Pan et al. confirmed the role of circ_0025033 overexpression in promoting cell growth, invasion, and metastasis. Circ_0025033 sponges miR-1231 and miR-1304 to stimulate clone formation and inhibit apoptosis [48]. In another study, circ_0039411 was elevated in PTC and cell lines (K1, TPC-1, FTC-133, SW579). At the same time, circ_0039411 knockdown inhibited cell proliferation/motility and accelerated apoptosis in PTC cells by regulating miR-1179 and miR-1205, which regulate the expression of the lipid “ATP-binding cassette transporter A9 (ABCA9)” and “metastasis-associated 1 (MTA1)” [49].
Yao and colleagues identified overexpression of circ_0058124 in PTC tissues and cell lines (K1, TPC-1). Circ_0058124 can modulate miR-218-5p and NUMB (i.e., target gene) expression, which in turn regulates the “NOTCH3/GATAD2A” signaling axis both in vitro and in vivo, suggesting that circ_0058124 may be considered as a promising target for therapeutic intervening in the progression of PTC [110]. The circRNA BTB domain and CNC homolog 2 (circ-BACH2) were overexpressed in PTC cell lines relative to normal thyroid follicular epithelial cells (Nthyori 3–1). Circ-BACH2 was predominately localized in the cytoplasm and may sponge miR-139-5p with subsequent “LIM domain only 4 (LMO4)” upregulation and PTC development [3]. In another study by Liu and colleagues, the circ_0058124 was upregulated in PTC tissues and cells (IHH-4 and TPC-1), acting as an oncogene by mediating the “miR-370-3p/LMO4” axis. Its knockdown inhibited cell viability, colony formation, migration/invasion and blocked tumor growth in vivo, and promoted apoptosis in vitro [131].
Table 3. Upregulated circRNAs promote thyroid tumor progression.
Table 3. Upregulated circRNAs promote thyroid tumor progression.
CircRNAParental GeneModelSourcemiRNATarget GeneMechanismRef.
circ-TP53TP53In vitro, in vivoTPC-1↓miR-1233-3pMDM2, ↓TP53Promotes cancer cell viability and proliferation[30]
circ_0004458 In vitro, in vivo, human (N = 48)BCPAP, TPC-1, K1, IHH4↓miR-885-5pRAC1Promotes cell growth and inhibits apoptosis and cell cycle arrest[116]
circ-ABCB10ABCB10In vitro, in vivo, human (N = 40) KLF6Enhances proliferation and invasion[47]
circ_0025033 In vitro ↓miR-1233, miR-1304 Promotes cell growth, invasion, and metastasis, triggers clone formation and prevents apoptosis[48]
circ_0039411 In vitro, human (N = 46)K1, TPC-1, FTC-133, SW579↓miR-1179, miR-1205MTA1, ABCA9Promotes cell proliferation/motility and inhibits apoptosis[49]
circ-BACH2BACH2In vitro ↓miR-139-5pLMO4Initiates tumor[3]
circ_0058124 In vivo, in vitro, human (N = 92)TPC-1, K1↓miR-218-5pNUMB, NOTCH3, GATAD2AActs as oncogene[110]
circ_0058124 In vitro, in vivoTPC-1, IHH-4↓miR-370-3pLMO4Promotes progression, colony formation, migration, invasion, inhibits apoptosis and tumor growth[131]
circ_0058124 In vitro, in vivo, human (N = 51)BCPAP, TPC-1, IHH-4, HTH83↓miR-940 Acts as oncogene[28]
circ-LDLRLDLRIn vitro, in vivo, human (N = 60) ↓miR-195-5pLIPHPromotes malignant behavior, colony formation, proliferation, migration, and invasion, and inhibits apoptosis[132]
circ_ZFRZFRIn vitro, human (N = 41)k1, TPC-1, SW579, FTC133↓miR-1261C8orf4Triggers tumor growth[31]
circ_102171 In vitro, in vivo, human (N = 47)TPC-1, NPA87, KAT-5 CTNNB1P1, TCF3, TCF4, LEF1Promotes progression, cell proliferation, migration, activation of the beta-catenin pathway, inhibits apoptosis[27]
circ_0011290 In vitroNIM-1, HTH83, TPC-1, K1, BCPAP↓miR-1252FSTL1Promotes cell viability, proliferation, apoptosis, and glucose metabolisms[32]
circ_0058129 In vitro, in vivo, human (N = 70)TPC-1, SNU-790↓miR-873-5pFSTL1Fosters cell proliferation, invasion, and migration[133]
circ_0005273 In vitroKTC-1, IHH-4, BCPAP, TPC-1↓miR-1138SOX2Promotes cell viability, proliferation, invasion, and migration[64]
circ_0005273 In vitro, human ↓miR-1138 Acts as oncogene[113]
circ-NEK6NEK6 ↓miR-370-3pFZD8Promotes progression, invasion, and metastasis, activates Wnt signaling pathway[73]
circ-FAT1(e2)FATI ↓miR-873ZEB1Promotes PTC cell growth, migration, and invasion[62]
circ-PSD3PSD3In vitro, in vivo ↓miR-637HEMGNActivates PI3K/Akt signaling, promote cell cycle progression, proliferation and metastasis and impedes apoptosis[109]
circ_0008274 In vitro, in vivo ↓miR-154-3pSLC7A11Regulates apoptosis, migration, and adhesion[111]
circ_0008274 In vitro, human AMPK, mTORPromotes proliferation/invasion through the “AMPK/mTOR” signaling pathway[112]
circ_0067934 In vitro, human (N = 57) PI3K, AKTPromotes EMT, PI3K/AKT signaling pathway, cell proliferation, invasion, migration, and inhibits apoptosis.[121]
circ_0067934 In vitro, in vivo ↓miR-1304CXCR1, PRKCIIncreases cell proliferation, migration, and invasion[113]
circ_0067934
circ-DOCK1DOCK1In vitroFTC-133, TPC-1↓miR-124 Induces tumorigenesis[134]
circ-PRMT5PRMT5In vitroK1, TPC-1, IHH4, BCPAP cells↓miR-30c Induces tumorigenesis[135]
circ-FOXM1FOXM1In vitro, in vivo, human (N = 78)K1, IHH-4, BCPAP, TCP-1↓miR-1179 Increases tumor growth[136]
circ_0011385 In vitro, in vivoBCPAP↓miR-361-3p Increases cell proliferation, invasion, and migration, inhibits apoptosis and regulates metastasis-related proteins[137]
circ-RAPGEF5RAPGEF5In vitro, in vivoBCPAP, KTC-1, K1, HEK293↓miR-198FGFR1Knockdown inhibits cell proliferation, migration, and invasion[138]
circ_0103552 In vitro, human (N = 56)TPC-1, SW579, 8505C Regulator for invasion and migration[139]
circ-UMAD1UMAD1Human (N = 50) ↓miR-873Gal3Poor prognosis, lymph node metastasis[115]
circ_IPCEF1 Human (N = 50) ↓miR-3619-5pCASR, CDC25B, NFκB1 and SHOC2Of 158 deregulated circRNAs, it showed the best predictive power[140]
circ_0001666 In vitro, in vivo, human ↓miR-330-5p, miR-193a-5p, miR-326ETV4Enhances G1 phase cell cycle progression[118]
circ-UBAP2UBAP2In vitro ↓miR-370-3p Regulates proliferation, apoptosis, and invasion[141]
circ-PVT1PVT1In vitro ↓ miR-195VEGFAActivates the Wnt/β-catenin signaling pathway[142]
circ-PUM1PUM1In vitro ↓miR-21MAPK1Promotes cell growth, metastasis, and glycolysis[143]
circ-PRKCIPRKCI ↓miR-335E2F3, MMP2, MMP9, snailPromotes cell progression and glycolysis[66]
circ_0079558 In vitro, human ↓miR-26b-5pMET/AKTFacilitates the proliferation and motility[144]
circ-RPS28RPS28In vitro ↓miR-345-5pFZD8, MMP2, MMP9, BAX1Regulates cell growth and motility[145]
circ-NRIP1NRIP1In vitro, in vivo, human ↓miR-195-5pp38, JAK2, STAT1Promotes cell proliferation/invasion and inhibits apoptosis[65]
circ_0001018 In vitro ↓ miR-338-3pSOX4Facilitates cell survival, invasion, G1/S cell cycle progression, and represses cell apoptosis[63]
circ_0067934 In vitro, in vivo, human ↓ miR-1301-3pHMGB1Increases growth, colony formation, migration, invasion, EMT[146]
circ_0062389 In vitro, human ↓miR-1179HMGB1Promotes proliferation, migration, and EMT[117]
circ-HIPK3 In vitro ↓miR-338-3pRAB23Promotes invasiveness, migration, invasion, and proliferation[147]
circ_0059354 In vitro, in vivo ↓miR-766-3pARFGEF1Aggravates cell proliferation, migration, invasion, and angiogenesis[123]
circ-NEURL4NEURL4 ↓ miR-1278LATS1Enhances cell proliferation and invasion[148]
circ-TIAM1TIAM1In vitro, in vivo, human (N = 60)K1, TPC-1, IHH-4, B-CPAP↓miR-646hnRNPA1Promotes cell migration, adhesion, growth, and polarity[149]
circ-ITGA7ITGA7In vitroTPC1 and CAL-62↓miR-198FGFR1Promotes migration and invasion[72]
ciRS-7 ↓miR-7EGFRPromotes tumor growth, EMT[78]
circ_0001681 ↓miR-942-5pTWIST1 [150]
Circ-LDLR is overexpressed in PTC tissues/cells. This overexpression regulated Lipase H (LIPH) expression by sponging miR-195-5p. Knockdown of this circRNA suppressed PTC cells’ malignant behaviors and promoted PTC cell colony apoptosis in vivo [132]. Circ_ZFR was also implicated in thyroid cancer. It can sponge and silence miR-1261, allowing the activation of the transcriptional and immune response regulator C8orf4 and triggering cancer cell growth [31]. The role of circ_102171 in PTC tissues and cell lines (TPC-1, NPA87, KAT-5) has also been examined. Higher expression was shown in PTC samples. Circ_102171 accelerated PTC progression by modulating “CTNNB1P1-dependent beta-catenin” pathway activation. It can target CTNNBIP1 by suppressing the “β-catenin/TCF3/TCF4/LEF1” pathway to initiate the “Wnt/β-catenin” pathway and progression of PTC. CircRNA_102171 downregulation decreased cell proliferation/migration and triggered apoptosis in vivo and in vitro [27].
Overexpression of circ_0011290 was reported in PTC cell lines (NIM-1, HTH83, TPC-1, K1, BCPAP). In vitro, circ_0011290 could positively modulate the protease inhibitor follistatin-like 1 (FSTL1) expression through sponging miR-1252. Suppression of miR-1252 reversed the phenotypic impacts of circ_0011290 silencing, including cell viability/proliferation and glucose metabolism [32]. Similarly, in PTC tissues and cells, circ_0058129 and FSTL1 abundances were increased, whereas the miR-873-5p content was decreased relative to control groups. Circ_0058129 shortage suppressed PTC cell proliferation/invasion and migration. Moreover, miR-873-5p repressed PTC cell malignancy by binding to FSTL1. Circ_0058129 targeted miR-873-5p to regulate FSTL1 [133]. Circ_0005273 was overexpressed in PTC cell lines (KTC-1, IHH-4, BCPAP, TPC-1). It can bind to miR-1138 and relieve suppression of “SRY (sex-determining region Y)-box 2 (SOX2)”. Knockdown of circ_0005273 hindered cell viability, proliferation, invasion, and migration of PTC. Silencing of circ_0005273 could suppress PTC tumor growth in vivo, making this circRNA a potential new therapeutic target for PTC [64,113].
By targeting miR-370-3p, Circ-NEK6 promotes cancer cell invasion, metastasis, and disease progression via transmembrane signaling receptor “frizzled family receptor 8 (FZD8)” upregulation and “Wnt signaling” pathway activation in TC [73]. CircFAT1(e2), derived from exon 2 of the FAT1 gene, was upregulated in PTC cell lines, and its knockdown inhibited PTC cell growth/invasion and migration. CircFAT1(e2) can sponge miR-873, promoting ZEB1 expression, which modulates cancer progression, invasion, and metastasis [62]. Circ_PSD3 was upregulated in PTC cell lines (TPC-1, CGTH-W3, IHH-4, SW579). Circ_PSD3 enhances PTC progression by regulating the “miR-637/HEMGN” axis and activating “PI3K/Akt” signaling. Circ_PSD3 and HEMGN facilitate cell cycle progression, proliferation, and metastasis and impedes PTC cell apoptosis. Circ_PSD3 silencing hampered tumor growth in vivo and facilitated apoptosis in vitro [109].
Zhou et al. suggested circ_0008274 to be a putative tumor-promoting circRNA in PTC. It was found to regulate PTC cell migration/adhesion and apoptosis in vitro. Its knockdown inhibited PTC cell migration/adhesion and enhanced apoptosis in vitro, and restrained tumor growth in vivo. At the cellular level, circ_0008274 modulated “solute carrier family 7-member A11 (SLC7A11)” expression by acting as a sponge of miR-154-3p, releasing the inhibiting effect on the transmembrane transporter SLC7A11. Hence, regulating the “miR-154-3p/SLC7A11” axis might provide a promising therapeutic target for PTC treatment [111]. Similar results were obtained by Zhou et al.; circ_0008274 was also significantly higher in PTC tissues and cell lines and was associated with poor prognosis. Circ_0008274 enhances PTC cell proliferation/invasion by regulating the “AMPK/mTOR” signaling pathway [112].
Circ_0067934 was highly expressed in TC cell lines. It acts as an oncogene, leading to the development of TC by promoting the EMT and “phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT)” signaling pathways [113]. Knockdown of circ_0067934 inhibited cell proliferation/migration and invasion and enhanced apoptosis [121]. By analogy to cell lines, circ_0067934 RNA was also highly expressed in TC compared with adjacent non-cancer tissues. Its overexpression was associated with shorter patient survival. Silencing of this circRNA inhibited cell proliferation and triggered apoptosis, as evidenced by cell counting/Transwell migration assays. Western blot results suggested this knockdown also prevented EMT and PI3K/AKT-pathway-related protein synthesis [113,121].
Circ_0011385 is overexpressed in TC. Its silencing inhibited TC cell proliferation/migration and invasion and supported cell cycle arrest/apoptosis. It has been reported to act as a sponge for miR-361-3p with negative regulation in vitro. It was also implicated in regulating several apoptosis/metastasis-related proteins in TC [137], supporting the value of the “circ_0011385/miR-361-3p” axis as a novel therapeutic opportunity for TC.
Circ-RAPGEF5 is a transcript of a five-exon RAPGEF5 gene expressed in PTC cells with uncertain oncogenic properties. Liu et al. evaluated the potential participation of this circRNA in TC pathogenesis and explored its upregulation in PTC tissues and cultures. They claimed that RAPGEF5 suppression could attenuate the transmembrane signaling receptor “fibroblast growth factor receptor 1 (FGFR1)” and thus suppress tumor growth by binding to miR-198, which sponged the 3′-UTR of FGFR1 via circRAPGEF5 targets in vitro. Xenograft experiments suggest that silencing circRAPGEF5 by prompting miR-198 expression can inhibit the oncogenic potential of FGFR1, making the “circRAPGEF5/miR-198/FGFR1” axis a promising target for PTC treatment [138]. Further studies showed that circ-0103552 was involved in the invasion/migration of TC via miR-127 sponging [139].
The concomitant presence of circulating circRNA-UMAD1 and Galectin-3 is considered a “co-biomarker” for predicting TC lymph node metastasis [115]. CircRNA profiling by microarray revealed a potential role of circ_IPCEF1 in PTC and showed interactions with four cancer-related genes, “CASR, CDC25B, NFκB1, and SHOC2” through sponging miR-3619-5p [140]. Additionally, circ_0001666 was upregulated in both PTC tissues and cell lines, and its expression was associated with lymph node metastasis. Silencing circ_0001666 inhibited cell proliferation, increased cell-cycle-associated proteins, arrested cell cycle progression at the G1 phase, and enhanced the expression of proapoptotic proteins. The oncogenic effect of circ_0001666 on tumor growth was also verified in xenograft experiments [118].
In another study, circUBAP2 was upregulated in TC tissues/cells, whereas its target miR-370-3p was downregulated. Silencing of this circRNA suppressed cell proliferation/invasion and induced apoptosis in TC cells [141]. Circular RNA Pvt1 oncogene (CircPVT1) supports the PTC progression by triggering the “Wnt/β-catenin” signaling pathway and modulating the ratio of miR-195/VEGFA expression [142]. Silencing of circPUM1 impairs cell growth, tissue glycolysis, and metastasis of PTC by sponging of miR-21-5p ensuing MAPK1 upregulation [143].
Knockdown of circRNA protein kinase C iota (circ-PRKCI) inhibits cell progression/glycolysis of PTC via the “circ-PRKCI/miR-335/E2F3” ceRNA axis [66]. Xenograft experiments indicated that the knockdown of this circRNA could hinder PTC cell proliferation and tumor growth in vivo [66]. Circ_0079558 stimulates PTC progression and facilitated proliferation/motility by sponging miR-26b-5p with activation of MET/AKT signaling [144]. CircRPS28 (hsa_circ_0049055) is an emerging player in regulating cell growth/motility of PTC by sponging miR-345-5p to regulate FZD8 [145]. Silencing this type of circRNA retained PTC cell viability, wound-healing rate, and numbers of colony formation and migration/invasion cells associated with apoptosis promotion. These effects are in line with decreased “B-cell lymphoma (Bcl)-2” levels and increased “Bcl-2-associated X protein (Bax), matrix metalloproteinase-2 (MMP2), and MMP9” levels as explored by Western blotting [145].
CircRNA NRIP1 was upregulated in PTC tissues/cells, and its high levels were associated with advanced PTC stage. Silencing of this circRNA suppressed PTC cell proliferation/invasion while accelerating apoptosis. Overexpression of miR-195-5p inhibited proliferation/invasion capabilities, triggering apoptosis of PTC cell lines, and retained the growth of tumor xenografts. These functions were reversed following circRNA NRIP1 upregulation in PTC cells/tumor xenografts. Protein levels of p38/JAK2/STAT1 were markedly downregulated in miR-195-5p overexpressed PTC cells/tumor xenografts, whereas circRNA NRIP1 overexpression negated these impacts [65].
Circ_0001018 promotes PTC by promoting cell survival/invasion, G1/S cell cycle progression, and inhibiting cell apoptosis via crosstalk with miR-338-3p and SOX4 [63]. Circ_0067934 supports the PTC cells’ progression via the “miR-1301-3p/high mobility group box 1 (HMGB1)” axis. Silencing of circ_0067934 inhibited growth, colony formation, migration/invasion, EMT, and growth of tumor xenograft and triggered apoptosis of PTC cells [146]. Circ_0062389 regulates PTC progression via the “miR-1179/HMGB1” axis and is significantly overexpressed in PTC tissues/cell lines. Upregulation of this circRNA was associated with large tumor size and positive lymph node metastasis. Silencing circ_0062389 could impede the PTC cells’ proliferation/migration and EMT [117].
In TC, another study suggested the potential role of circ_0001681 in the EMT process. This cytoplasmic circRNA was upregulated in TC cells and inversely related to the miR-942-5p level [150]. By targeting miR-942-5p and upregulating the TWIST1 expression, silencing of circ_0001681 significantly suppressed TC progression [150]. These findings indicate that circ_0001681 is a promising treatment candidate for TC.
The oncogenic circ_0059354 has been implicated in PTC progression by upregulating ARFGEF1 through miR-766-3p sponging. Circ_0059354 silencing facilitated apoptosis and suppressed cell proliferation/migration, invasion, and angiogenesis in PTC cells [123]. Inhibition of miR-766-3p reversed circ_0059354-knockdown-mediated effects, whereas its overexpression restrained the PTC cells’ malignant behaviors. At the same time, downregulation of circ_0059354 blocked tumor growth in vivo [123].
Circ-NRIP1 promotes the progress of PTC via sponging mir-195-5p and modulating the “P38 MAPK and JAK/STAT” pathways [65]. This circRNA was overexpressed in PTC tissues/cells, and its levels correlated with advanced stages of PTC. Silencing of this circRNA suppresses PTC cell proliferation/invasion while accelerating apoptosis. Upregulation of miR-195-5p repressed proliferation/invasion capabilities, accelerated apoptosis of PTC cell lines, and restrained the growth of tumor xenografts. In addition, protein levels of “p-p38, p-JAK2, and p-STAT1” were downregulated in miR-195-5p-overexpressed PTC cells and tumor xenografts, whereas CircRNA NRIP1 upregulation reversed these effects [65]. Furthermore, circHIPK3 sponges miR-338-3p and enhance RAB23 expression in TC cells, promoting cancer invasiveness. Silencing of circHIPK3 significantly reduced proliferation/invasion and migration of TC cells [147].
Circ_0079558 and MET were reported to be upregulated in PTC tissues/cell lines. This circRNA promoted PTC cell proliferation/migration and decreased apoptosis through the “miR-26b-5p/MET/AKT” axis [144]. Circ-TIAM1 induced cell migration of PTC via the “miR-646/heterogeneous ribonucleoprotein A1 (HNRNPA1)” axis [149]. Meanwhile, circNEURL4 downregulation in the PTC samples was directly associated with an unfavorable prognosis. Through binding to miR-1278, circNEURL4 might suppress cell proliferation/invasion of PTC in vivo and in vitro by indirectly increasing the expression of LATS1 [148]. Thus, this circRNA could have a clinical diagnostic and/or prognostic utility for PTC patients and a potential targeted therapeutic for PTC via the “miR-1278/LATS1” axis.

3.3. Tumor-Cell-Derived Exosomal circRNAs

Exosomes are created in multivesicular endosomes and can be excreted from several cellular types. They are implicated in intercellular communication by transferring intracellular cargoes, such as proteins/nucleic acids [151]. RNA-seq analyses indicated that circRNAs were enriched in the exosomal compartment compared with the original cells [44,45], and a further study has identified that circRNA levels in exosomes could mirror their levels in cells/tissues [57]. Also, circRNAs sorting to exosomes may be controlled by changes in the levels of associated miRNA in producer cells, and the biological activities could be transferred to recipient cells. The abundance of tumor-derived serum “exosomal circRNAs (exo-circRNA)” in xenografted mice was correlated with tumor mass [44]. Exo-circRNAs could discriminate patients with colon cancer from healthy controls, suggesting the putative biological function of exosomal circRNAs [92]. Emerging evidence shows the excellent promise of exosomal circRNAs in tumor immunity regulation, patient outcome prediction, and drug efficacy evaluation [152].
Exosomal circ_007293 was reported to promote proliferation/invasion, migration, and EMT of PTC cells through regulation of the “miR-653-5p/paired box 6 (PAX6)” axis [153], suggesting that tumor-derived exosomal circRNAs may function as potential biomarkers for cancer progression. Additionally, “circRNA-UMAD1 and Gal3” were recognized as having more substantial co-biomarker potential with significant upregulation in the sera of patients with LNM, relative to those with primary tumors, as verified by the serum RNA expression levels of 50 patients with PTC associated with or without LNM by “quantitative real-time PCR” [115]. Further studies using liquid biopsy samples are recommended to facilitate dynamic monitoring of cancer and postoperative surveillance.

3.4. Circular RNAs and Treatment Resistance

Since the sensitivity of tumors to chemotherapy can impact the survival/prognosis of patients, deciphering the mechanisms underlying drug resistance is critical. Dysregulation of circRNAs in cancer was found to induce chemoresistance. CircHIPK3 overexpression was found to promote gemcitabine (GEM) resistance in pancreatic cancer cells by targeting “RASSF1” through miR-330-5p and was proposed to be a novel biomarker in GEM-resistant PC [154]. Similarly, “circRNA Cdr1as” induced apoptosis and increased the cisplatin chemosensitivity of bladder cancer cells both in vitro and in vivo by upregulating “APAF1” expression through miR-1270 inhibition [155]. Sang et al. indicated that circRNA_0025202 modulated tamoxifen sensitivity and tumor progress through the “miR-182-5p/FOXO3a” signaling pathway in breast cancer [156].
In thyroid carcinoma, Liu et al. [119] recognized that “circRNA EIF6 (hsa_circ_0060060)” sponges miR-144-3p to promote resistance to cisplatin in PTC and ATC cell lines by regulation of autophagy. CircEIF6 targets hsa-miR-144-3p, which modulates the expression of transforming growth factor α (TGF-α) and other autophagy-related proteins. In contrast, circEIF6 silencing augmented cisplatin sensitivity in vivo by mediating autophagy [119]. The above studies indicate that circRNAs are promising therapeutic targets for reversing drug resistance in tumors (Figure 5).

3.5. Functional Enrichment Analysis for the Deregulated Markers

Prior publications have shown the role of deregulated circular RNAs in various hallmarks of cancer (Figure 6). Circ_0001018 activates cell proliferation by regulating the miR-338-3p/SOX4 axis [63]. Circ-ITCH and similar circRNAs promote cancer cells, evading antigrowth signals by enhancing other oncogenes, such as CBL [71]. CircEIF6/miR-144-3p promotes cancer cells by evading cell death via regulating cellular apoptosis or autophagy [119]. Circ-PSD3 limits the replicative potential of cancer cells and impedes apoptosis by regulating HEMGN [109]. Circ_0059354/miR-766-3p sustains angiogenesis through regulating ARFGEF1 [123]. CiRS-7/miR-7/EGFR [78], circ_0067934/PI3K/AKT [121], circ_0067934/miR-1301-3p/HMGB1 [146], and circ_0062389/miR-1179/HMGB1 [117] regulate the EMT process and thus cancer tissue invasion/metastasis.
Collectively, 34 microRNAs and 51 target genes are known to be deregulated by circRNAs in thyroid cancer. Pathway enrichment analysis of 51 target genes in the reviewed articles demonstrated their role in inflammation mediated by chemokine and cytokine signaling (MAPK1, JAK2, TCF3, STAT1, RAC1, CXCR1, and p38), the Wnt signaling pathway (TCF3, PRKC1, TP53, FSTL1, FZD8, and LEF1), angiogenesis (MAPK1, FGFR1, VEGFA, p38, STAT1, and PRKCI), CCKR signaling map (MAPK1, JAK1, TCF4, MP9, p38, and RAC1), and the PDGF signaling pathway (MAPK1, JAK2, mTOR < STAT1, and p38) (Figure 7).

3.6. Genome-Wide Circular RNA Screening

3.6.1. Data Source and Processing

Microarray transcriptomic signature of circRNAs in human PTC was obtained from the “Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/)”. A systematic search revealed a single dataset for circRNA analysis in TC (GSE93522). It contains 18 samples, of which 6 PTCs were compared with their matching contralateral normal thyroid tissues. The platform used was “GPL19978 Agilent-069978 Arraystar Human CircRNA microarray V1”.
The raw data were analyzed via GEO2R software (http://www.ncbi.nlm.nih.gov/geo/geo2r/) using R version 3.2.3 software packages (Biobase 2.30.0, GEOquery 2.40.0, limma 3.26.8), and the results were confirmed using NetworkAnalyst (www.networkanalyst.ca). Differentially expressed circRNA (DEC) analysis was performed using a p-value threshold of <0.05 and log fold change (FC) >1. Annotation files from the platform were used to identify the corresponding circRNA names of the probes, and genes from which circRNAs are derived were identified using circBase (www.circBase.org). A volcano plot was used to visualize the significantly differential circRNAs between PTC and controls. A heatmap was constructed to display the differential circRNA expression patterns among the samples under R version 4.0.3 using the matrix data. DEC was plotted on ideograms using Phenograms (https://ritchielab.org/software/phenogram) (all online tools were last accessed on 23 June 2022).

3.6.2. Identification of Differentially Expressed circRNAs (DEC)

Characteristics of human thyroid specimens are shown in Figure 8A–C. In the current analysis, we identified 137 significant DEC in PTC compared with controls (Figure 8D). The expression pattern of these DEC enables clear demarcation between cancer and normal thyroid tissues (Figure 8E–F). The top 10 downregulated circRNAs were circ_0000104, circ_0012171, circ_0014161, circ_0000031, circ_0011969, circ_0010832, circ_0010777, circ_0012226, circ_0010395, circ_0011213. The top upregulated genes were circ_0013809, circ_0010146, circ_0011287, circ_0012451, circ_0014700, circ_0040024, circ_0011748, circ_0014566, circ_0014940, circ_0014387 (Figure 8G). Their spliced mature sequence length ranged from 71 to 141,103 base pairs. A total of 74 are transcribed from the antisense strand, whereas 63 DEC are generated from the sense strand. Of 2,895 studied circRNAs, 92.2% (N = 2668) were located along chromosome 1, and 2.38% of genes existed in chromosome 16 (N = 69) (Figure 8H). These included 115 upregulated genes, of which 105 circRNAs lie on chromosome 1. In contrast, the 22 downregulated circRNAs are in chromosome 1. Significant circRNAs were clustered at the short arm of chromosome 1 (1p13.3, 1p32.2), 1q21.3, and 1p34.2 to 1p36.13), long arm of chromosome 1 (1q22), and chromosome 16q22.1. (Figure 8I). The list of DECs and their annotation is shown in Supplementary Table S1.

4. Discussion of the Current Challenges in Clinical Practice and Future Perspectives

The first wave of human disease circRNA research showed great potential in human diseases. Moving forward, a better functional and clinical understanding of circRNA requires the implementation of best practice principles in future research; however, some experimental challenges should be considered. For purification, standard RNA preparation kits can isolate circRNAs [20], but circRNA will be largely missing in poly(A)-enriched samples as they lack poly(A) tails.
On the other hand, the removal of linear RNA challenges the assessment of whether the linear host genome and the circRNA expression are independent of one another [20]. It may be beneficial to begin with total RNA but allow enrichment for circRNA using an RNAase R treatment with Li+ (rather than K+) [157], poly(A) counterselection [158], and rRNA removal to increase the number of circRNA reads [20]. Assessment of RNA quantity, quality, and purity requires different profiling methods for circRNA [159].
When quantifying circRNA in RNA-seq data, the overlay is preferred between more than one bioinformatic algorithm (such as find_circ, CIRI2, CIRCexplorer2, and circRNA_finder) to minimize the potential for false-positive circRNAs [160]. Additional validation measures to examine the circular structure should also be undertaken. When detecting circRNAs using a microscope, fixed cells should be treated with a protease prior to probe hybridization to remove their association with ribosome binding proteins (RBP), which may inhibit probe binding and result in the absence of signals. Positive and negative controls should be included to check probe specificity [20].
Since circRNA expression is often specific to cell type, developing techniques for single-cell profiling or spatial transcriptomics for circRNAs offers resolution at a single-cell level, which is especially important in complex tumor tissues. Validating and predicting gene and drug targets will provide comprehensive mechanistic approaches. Although detection methods are still in their infancy, primarily confined to microarray analysis, early generation sequencing analyses suggest the possibility of circRNA detection transcriptome wide. Currently, circRNA validation is based mainly on circular structure confirmation and detection of its back splice junction. Moving forward, more accurate circRNA detection and validation methods are an area of interest in this field.
One of the limitations we had in the current analysis was the lack of circRNA detection in patients with other histopathological variants such as anaplastic, follicular, or medullary thyroid cancers. Moreover, studies tracing dynamic longitudinal monitoring of circRNAs during surveillance, postoperatively, or following radiotherapy were unavailable. The discovery of other circRNAs and elucidating the factors mentioned above will significantly improve the future of gene-targeted therapy.

5. Conclusions

Cancer remains a significant challenge for healthcare providers. Advances in its diagnosis and treatment are critical moving forward. This review aimed to elucidate TC development, progression, and treatment concerning circRNAs. The literature has demonstrated the benefits of using circRNAs as stable and tissue-specific biomarkers. Although the field is still in its early stages, we believe that circRNAs could serve essential roles in the diagnosis, prognosis, and treatment of cancers.
This work summarized the role of circRNAs in thyroid cancer and their ability to act as miRNA sponges, regulate transcription and translation, and serve as protein decoys. Consequently, subsequent determination of circRNA deregulation mechanisms could optimize their utility in the diagnosis, prognosis, and treatment of thyroid cancer. Importantly, circRNAs can be detected in patient plasma and potentially be used to track cancer progression. The advantages of circRNA, including its higher stability, wider abundance, and occurrence in several body fluids, potentiate its use as promising prognostic and theranostic biomarkers in TC.
It has been established that circRNAs interact with miRNAs to influence downstream targets playing vital roles in TC, such as metastasis, angiogenesis, and therapeutic resistance. Future work to elucidate these regulatory networks, as well as the significance of their manipulation and exploitation, will provide a comprehensive “Rosetta Stone” of this RNA dialect.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14194728/s1, Table S1: Fold change and annotation of differentially expressed circRNAs in PTC compared with normal thyroid tissues.

Author Contributions

All authors contributed significantly to the creation of this manuscript, each having fulfilled criteria established by the ICMJE. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a research grant (THYROIDGRANT2021-0000000232) from “Bite Me Cancer” and facilitated by the “American Thyroid Association” (to E.T.) and “Tulane University Bridge Fund” (to E.K.).

Acknowledgments

A sincere thanks to Loula Burton from Tulane’s Research Proposal Development Office for her diligent editing and proofreading of this paper. Many thanks are also extended to Peter Issa (A medical student at LSU) for language revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cabanillas, M.E.; McFadden, D.G.; Durante, C. Thyroid cancer. Lancet 2016, 388, 2783–2795. [Google Scholar] [CrossRef]
  2. Kitahara, C.M.; Sosa, J.A. The changing incidence of thyroid cancer. Nat. Rev. Endocrinol. 2016, 12, 646–653. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, X.; Zhao, Z.; Dong, J.; Lv, Q.; Yun, B.; Liu, J.; Shen, Y.; Kang, J.; Li, J. Circular RNA circBACH2 plays a role in papillary thyroid carcinoma by sponging miR-139-5p and regulating LMO4 expression. Cell Death Dis. 2019, 10, 1–12. [Google Scholar] [CrossRef]
  4. Toraih, E.A.; Fawzy, M.S.; Mohammed, E.A.; Hussein, M.H.; El-Labban, M.M. MicroRNA-196a2 Biomarker and Targetome Network Analysis in Solid Tumors. Mol. Diagn. Ther. 2016, 20, 559–577. [Google Scholar] [CrossRef] [PubMed]
  5. Toraih, E.A.; Aly, N.M.; Abdallah, H.Y.; Al-Qahtani, S.A.; Shaalan, A.A.; Hussein, M.H.; Fawzy, M.S. MicroRNA-target cross-talks: Key players in glioblastoma multiforme. Tumour Biol. 2017, 39, 1010428317726842. [Google Scholar] [CrossRef]
  6. Toraih, E.A.; Mohammed, E.A.; Farrag, S.; Ramsis, N.; Hosny, S. Pilot Study of Serum MicroRNA-21 as a Diagnostic and Prognostic Biomarker in Egyptian Breast Cancer Patients. Mol. Diagn. Ther. 2015, 19, 179–190. [Google Scholar] [CrossRef] [PubMed]
  7. Toraih, E.A.; Alghamdi, S.A.; El-Wazir, A.; Hosny, M.; Hussein, M.H.; Khashana, M.; Fawzy, M.S. Dual biomarkers long non-coding RNA GAS5 and microRNA-34a co-expression signature in common solid tumors. PLoS ONE 2018, 13, e0198231. [Google Scholar] [CrossRef]
  8. Hsu, M.-T.; Coca-Prados, M. Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells. Nature 1979, 280, 339–340. [Google Scholar] [CrossRef]
  9. Du, W.W.; Fang, L.; Yang, W.; Wu, N.; Awan, F.M.; Yang, Z.; Yang, B.B. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ. 2016, 24, 357–370. [Google Scholar] [CrossRef]
  10. Werfel, S.; Nothjunge, S.; Schwarzmayr, T.; Strom, T.-M.; Meitinger, T.; Engelhardt, S. Characterization of circular RNAs in human, mouse and rat hearts. J. Mol. Cell. Cardiol. 2016, 98, 103–107. [Google Scholar] [CrossRef] [Green Version]
  11. Lasda, E.; Parker, R. Circular RNAs: Diversity of form and function. RNA 2014, 20, 1829–1842. [Google Scholar] [CrossRef]
  12. Capel, B.; Swain, A.; Nicolis, S.; Hacker, A.; Walter, M.; Koopman, P.; Goodfellow, P.; Lovell-Badge, R. Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 1993, 73, 1019–1030. [Google Scholar] [CrossRef]
  13. Cocquerelle, C.; Mascrez, B.; Hétuin, D.; Bailleul, B. Mis-splicing yields circular RNA molecules. FASEB J. 1993, 7, 155–160. [Google Scholar] [CrossRef]
  14. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Zhang, X.-O.; Chen, T.; Xiang, J.-F.; Yin, Q.-F.; Xing, Y.-H.; Zhu, S.; Yang, L.; Chen, L.-L. Circular Intronic Long Noncoding RNAs. Mol. Cell 2013, 51, 792–806. [Google Scholar] [CrossRef]
  16. Ashwal-Fluss, R.; Meyer, M.; Pamudurti, N.R.; Ivanov, A.; Bartok, O.; Hanan, M.; Evantal, N.; Memczak, S.; Rajewsky, N.; Kadener, S. circRNA biogenesis competes with pre-mRNA splicing. Mol. Cell 2014, 56, 55–66. [Google Scholar] [CrossRef]
  17. Kristensen, L.S.; Hansen, T.B.; Venø, M.T.; Kjems, J. Circular RNAs in cancer: Opportunities and challenges in the field. Oncogene 2018, 37, 555–565. [Google Scholar] [CrossRef]
  18. Liu, J.; Zhang, X.; Yan, M.; Li, H. Emerging Role of Circular RNAs in Cancer. Front. Oncol. 2020, 10, 663. [Google Scholar] [CrossRef]
  19. Brown, J.R.; Chinnaiyan, A.M. The Potential of Circular RNAs as Cancer Biomarkers. Cancer Epidemiol. Biomark. Prev. 2020, 29, 2541–2555. [Google Scholar] [CrossRef]
  20. Nielsen, A.F.; Bindereif, A.; Bozzoni, I.; Hanan, M.; Hansen, T.B.; Irimia, M.; Kadener, S.; Kristensen, L.S.; Legnini, I.; Morlando, M.; et al. Best practice standards for circular RNA research. Nat. Methods 2022, 1–13. [Google Scholar] [CrossRef]
  21. Hu, W.; Bi, Z.; Chen, Z.-L.; Liu, C.; Li, L.-L.; Zhang, F.; Zhou, Q.; Song, Y.-Y.; Zhan, B.-T.; Zhang, Q.; et al. Emerging landscape of circular RNAs in lung cancer. Cancer Lett. 2018, 427, 18–27. [Google Scholar] [CrossRef]
  22. Han, D.; Li, J.; Wang, H.; Su, X.; Hou, J.; Gu, Y.; Qian, C.; Lin, Y.; Liu, X.; Huang, M.; et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology 2017, 66, 1151–1164. [Google Scholar] [CrossRef]
  23. Hsiao, K.-Y.; Lin, Y.-C.; Gupta, S.K.; Chang, N.; Yen, L.; Sun, H.S.; Tsai, S.-J. Noncoding Effects of Circular RNA CCDC66 Promote Colon Cancer Growth and Metastasis. Cancer Res. 2017, 77, 2339–2350. [Google Scholar] [CrossRef]
  24. Su, H.; Tao, T.; Yang, Z.; Kang, X.; Zhang, X.; Kang, D.; Wu, S.; Li, C. Circular RNA cTFRC acts as the sponge of MicroRNA-107 to promote bladder carcinoma progression. Mol. Cancer 2019, 18, 1–15. [Google Scholar] [CrossRef]
  25. Xia, W.; Qiu, M.; Chen, R.; Wang, S.; Leng, X.; Wang, J.; Xu, Y.; Hu, J.; Dong, G.; Xu, L.; et al. Circular RNA has_circ_0067934 is upregulated in esophageal squamous cell carcinoma and promoted proliferation. Sci. Rep. 2016, 6, 35576. [Google Scholar] [CrossRef]
  26. Yang, R.; Xing, L.; Zheng, X.; Sun, Y.; Wang, X.; Chen, J. RETRACTED ARTICLE: The circRNA circAGFG1 acts as a sponge of miR-195-5p to promote triple-negative breast cancer progression through regulating CCNE1 expression. Mol. Cancer 2019, 18, 1–19. [Google Scholar] [CrossRef]
  27. Bi, W.; Huang, J.; Nie, C.; Liu, B.; He, G.; Han, J.; Pang, R.; Ding, Z.; Xu, J.; Zhang, J. CircRNA circRNA_102171 promotes papillary thyroid cancer progression through modulating CTNNBIP1-dependent activation of β-catenin pathway. J. Exp. Clin. Cancer Res. 2018, 37, 1–9. [Google Scholar] [CrossRef]
  28. Borran, S.; Ahmadi, G.; Rezaei, S.; Anari, M.M.; Modabberi, M.; Azarash, Z.; Razaviyan, J.; Derakhshan, M.; Akhbari, M.; Mirzaei, H. Circular RNAs: New players in thyroid cancer. Pathol.-Res. Pract. 2020, 216, 153217. [Google Scholar] [CrossRef]
  29. Zhong, Z.; Huang, M.; Lv, M.; He, Y.; Duan, C.; Zhang, L.; Chen, J. Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway. Cancer Lett. 2017, 403, 305–317. [Google Scholar] [CrossRef]
  30. Ma, W.; Zhao, P.; Zang, L.; Zhang, K.; Liao, H.; Hu, Z. CircTP53 promotes the proliferation of thyroid cancer via targeting miR-1233-3p/MDM2 axis. J. Endocrinol. Investig. 2020, 44, 353–362. [Google Scholar] [CrossRef]
  31. Wei, H.; Pan, L.; Tao, D.; Li, R. Circular RNA circZFR contributes to papillary thyroid cancer cell proliferation and invasion by sponging miR-1261 and facilitating C8orf4 expression. Biochem. Biophys. Res. Commun. 2018, 503, 56–61. [Google Scholar] [CrossRef] [PubMed]
  32. Hu, Z.; Zhao, P.; Zhang, K.; Zang, L.; Liao, H.; Ma, W. Hsa_circ_0011290 regulates proliferation, apoptosis and glycolytic phenotype in papillary thyroid cancer via miR-1252/ FSTL1 signal pathway. Arch. Biochem. Biophys. 2020, 685, 108353. [Google Scholar] [CrossRef] [PubMed]
  33. Seal, R.L.; Chen, L.; Griffiths-Jones, S.; Lowe, T.M.; Mathews, M.B.; O’Reilly, D.; Pierce, A.J.; Stadler, P.F.; Ulitsky, I.; Wolin, S.L.; et al. A guide to naming human non-coding RNA genes. EMBO J. 2020, 39, e103777. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, Y.; Xue, W.; Li, X.; Zhang, J.; Chen, S.; Zhang, J.-L.; Yang, L.; Chen, L.-L. The Biogenesis of Nascent Circular RNAs. Cell Rep. 2016, 15, 611–624. [Google Scholar] [CrossRef] [PubMed]
  35. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef] [PubMed]
  36. Kristensen, L.S.; Andersen, M.S.; Stagsted, L.V.W.; Ebbesen, K.K.; Hansen, T.B.; Kjems, J. The biogenesis, biology and characterization of circular RNAs. Nat. Rev. Genet. 2019, 20, 675–691. [Google Scholar] [CrossRef]
  37. Salzman, J.; Chen, R.E.; Olsen, M.N.; Wang, P.L.; Brown, P.O. Cell-type specific features of circular RNA expression. PLoS Genet. 2013, 9, e1003777. [Google Scholar] [CrossRef]
  38. Huang, A.; Zheng, H.; Wu, Z.; Chen, M.; Huang, Y. Circular RNA-protein interactions: Functions, mechanisms, and identification. Theranostics 2020, 10, 3503–3517. [Google Scholar] [CrossRef]
  39. Wu, Q.; Li, P.; Wu, M.; Liu, Q. Deregulation of Circular RNAs in Cancer from the Perspectives of Aberrant Biogenesis, Transport and Removal. Front. Genet. 2019, 10, 16. [Google Scholar] [CrossRef]
  40. Chen, L.-L.; Yang, L. Regulation of circRNA biogenesis. RNA Biol. 2015, 12, 381–388. [Google Scholar] [CrossRef]
  41. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef]
  42. Kristensen, L.S.; Jakobsen, T.; Hager, H.; Kjems, J. The emerging roles of circRNAs in cancer and oncology. Nat. Rev. Clin. Oncol. 2021, 19, 188–206. [Google Scholar] [CrossRef]
  43. Zhou, C.; Molinie, B.; Daneshvar, K.; Pondick, J.V.; Wang, J.; Van Wittenberghe, N.; Xing, Y.; Giallourakis, C.C.; Mullen, A.C. Genome-Wide Maps of m6A circRNAs Identify Widespread and Cell-Type-Specific Methylation Patterns That Are Distinct from mRNAs. Cell Rep. 2017, 20, 2262–2276. [Google Scholar] [CrossRef]
  44. Li, Y.; Zheng, Q.; Bao, C.; Li, S.; Guo, W.; Zhao, J.; Chen, D.; Gu, J.; He, X.; Huang, S. Circular RNA is enriched and stable in exosomes: A promising biomarker for cancer diagnosis. Cell Res. 2015, 25, 981–984. [Google Scholar] [CrossRef]
  45. Dou, Y.; Cha, D.J.; Franklin, J.L.; Higginbotham, J.N.; Jeppesen, D.K.; Weaver, A.M.; Prasad, N.; Levy, S.; Coffey, R.J.; Patton, J.G.; et al. Circular RNAs are down-regulated in KRAS mutant colon cancer cells and can be transferred to exosomes. Sci. Rep. 2016, 6, 37982. [Google Scholar] [CrossRef]
  46. Wang, S.; Dong, Y.; Gong, A.; Kong, H.; Gao, J.; Hao, X.; Liu, Y.; Wang, Z.; Fan, Y.; Liu, C.; et al. Exosomal circRNAs as novel cancer biomarkers: Challenges and opportunities. Int. J. Biol. Sci. 2021, 17, 562–573. [Google Scholar] [CrossRef]
  47. Han, X.-T.; Jiang, J.-Q.; Li, M.-Z.; Cong, Q.-M. Circular RNA circ-ABCB10 promotes the proliferation and invasion of thyroid cancer by targeting KLF6. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 1271–1277. [Google Scholar]
  48. Pan, Y.; Xu, T.; Liu, Y.; Li, W.; Zhang, W. Upregulated circular RNA circ_0025033 promotes papillary thyroid cancer cell proliferation and invasion via sponging miR-1231 and miR-1304. Biochem. Biophys. Res. Commun. 2019, 510, 334–338. [Google Scholar] [CrossRef]
  49. Yang, Y.; Ding, L.; Li, Y.; Xuan, C. Hsa_circ_0039411 promotes tumorigenesis and progression of papillary thyroid cancer by miR-1179/ABCA9 and miR-1205/MTA1 signaling pathways. J. Cell. Physiol. 2020, 235, 1321–1329. [Google Scholar] [CrossRef]
  50. Zheng, Q.; Bao, C.; Guo, W.; Li, S.; Chen, J.; Chen, B.; Luo, Y.; Lyu, D.; Li, Y.; Shi, G.; et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat. Commun. 2016, 7, 11215. [Google Scholar] [CrossRef]
  51. Liu, C.-X.; Li, X.; Nan, F.; Jiang, S.; Gao, X.; Guo, S.-K.; Xue, W.; Cui, Y.; Dong, K.; Ding, H.; et al. Structure and Degradation of Circular RNAs Regulate PKR Activation in Innate Immunity. Cell 2019, 177, 865–880.e21. [Google Scholar] [CrossRef] [PubMed]
  52. Li, Z.; Huang, C.; Bao, C.; Chen, L.; Lin, M.; Wang, X.; Zhong, G.; Yu, B.; Hu, W.; Dai, L.; et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat. Struct. Mol. Biol. 2015, 22, 256–264, Erratum in Nat. Struct. Mol. Biol. 2017, 24, 194. [Google Scholar] [CrossRef] [PubMed]
  53. Meng, S.; Zhou, H.; Feng, Z.; Xu, Z.; Tang, Y.; Li, P.; Wu, M. CircRNA: Functions and properties of a novel potential biomarker for cancer. Mol. Cancer 2017, 16, 1–8. [Google Scholar] [CrossRef] [PubMed]
  54. Armakola, M.; Higgins, M.J.; Figley, M.D.; Barmada, S.J.; Scarborough, E.A.; Diaz, Z.; Fang, X.; Shorter, J.; Krogan, N.J.; Finkbeiner, S.; et al. Inhibition of RNA lariat debranching enzyme suppresses TDP-43 toxicity in ALS disease models. Nat. Genet. 2012, 44, 1302–1309. [Google Scholar] [CrossRef]
  55. Du, W.W.; Yang, W.; Liu, E.; Yang, Z.; Dhaliwal, P.; Yang, B.B. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016, 44, 2846–2858. [Google Scholar] [CrossRef]
  56. Amicone, L.; Marchetti, A.; Cicchini, C. Exosome-Associated circRNAs as Key Regulators of EMT in Cancer. Cells 2022, 11, 1716. [Google Scholar] [CrossRef]
  57. Dai, X.; Chen, C.; Yang, Q.; Xue, J.; Chen, X.; Sun, B.; Luo, F.; Liu, X.; Xiao, T.; Xu, H.; et al. Exosomal circRNA_100284 from arsenite-transformed cells, via microRNA-217 regulation of EZH2, is involved in the malignant transformation of human hepatic cells by accelerating the cell cycle and promoting cell proliferation. Cell Death Dis. 2018, 9, 1–14. [Google Scholar] [CrossRef]
  58. Chen, N.; Zhao, G.; Yan, X.; Lv, Z.; Yin, H.; Zhang, S.; Song, W.; Li, X.; Li, L.; Du, Z.; et al. A novel FLI1 exonic circular RNA promotes metastasis in breast cancer by coordinately regulating TET1 and DNMT1. Genome Biol. 2018, 19, 1–14. [Google Scholar] [CrossRef]
  59. Zhao, X.; Cai, Y.; Xu, J. Circular RNAs: Biogenesis, Mechanism, and Function in Human Cancers. Int. J. Mol. Sci. 2019, 20, 3926. [Google Scholar] [CrossRef]
  60. Bolisetty, M.T.; Graveley, B.R. Circuitous Route to Transcription Regulation. Mol. Cell 2013, 51, 705–706. [Google Scholar] [CrossRef]
  61. Gao, X.; Ma, X.-K.; Li, X.; Li, G.-W.; Liu, C.-X.; Zhang, J.; Wang, Y.; Wei, J.; Chen, J.; Chen, L.-L.; et al. Knockout of circRNAs by base editing back-splice sites of circularized exons. Genome Biol. 2022, 23, 1–22. [Google Scholar] [CrossRef]
  62. Liu, J.; Li, H.; Wei, C.; Ding, J.; Lu, J.; Pan, G.; Mao, A. circFAT1(e2) Promotes Papillary Thyroid Cancer Proliferation, Migration, and Invasion via the miRNA-873/ZEB1 Axis. Comput. Math. Methods Med. 2020, 2020, 1–9. [Google Scholar] [CrossRef]
  63. Luo, Q.; Guo, F.; Fu, Q.; Sui, G. hsa_circ_0001018 promotes papillary thyroid cancer by facilitating cell survival, invasion, G1/S cell cycle progression, and repressing cell apoptosis via crosstalk with miR-338-3p and SOX4. Mol. Ther.-Nucleic Acids 2021, 24, 591–609. [Google Scholar] [CrossRef]
  64. Zhang, W.; Zhang, H.; Zhao, X. circ_0005273 promotes thyroid carcinoma progression by SOX2 expression. Endocr.-Relat. Cancer 2020, 27, 11–21. [Google Scholar] [CrossRef]
  65. Li, C.; Zhu, L.; Fu, L.; Han, M.; Li, Y.; Meng, Z.; Qiu, X. CircRNA NRIP1 promotes papillary thyroid carcinoma progression by sponging mir-195-5p and modulating the P38 MAPK and JAK/STAT pathways. Diagn. Pathol. 2021, 16, 1–11. [Google Scholar] [CrossRef]
  66. Liu, Y.; Chen, G.; Wang, B.; Wu, H.; Zhang, Y.; Ye, H. Silencing circRNA protein kinase C iota (circ-PRKCI) suppresses cell progression and glycolysis of human papillary thyroid cancer through circ-PRKCI/miR-335/E2F3 ceRNA axis. Endocr. J. 2021, 68, 713–727. [Google Scholar] [CrossRef]
  67. Fabian, M.R.; Sonenberg, N.; Filipowicz, W. Regulation of mRNA Translation and Stability by microRNAs. Annu. Rev. Biochem. 2010, 79, 351–379. [Google Scholar] [CrossRef]
  68. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef]
  69. Shifman, B.M.; Platonova, N.M.; Vasilyev, E.V.; Abdulkhabirova, F.M.; Kachko, V.A. Circular RNAs and thyroid cancer: Closed molecules, open possibilities. Crit. Rev. Oncol. 2022, 173, 103662. [Google Scholar] [CrossRef]
  70. Bai, Q.; Pan, Z.; Nabi, G.; Rashid, F.; Liu, Y.; Khan, S. Emerging role of competing endogenous RNA and associated noncoding RNAs in thyroid cancer. Am. J. Cancer Res. 2022, 12, 961–973. [Google Scholar]
  71. Wang, M.; Chen, B.; Ru, Z.; Cong, L. CircRNA circ-ITCH suppresses papillary thyroid cancer progression through miR-22-3p/CBL/β-catenin pathway. Biochem. Biophys. Res. Commun. 2018, 504, 283–288. [Google Scholar] [CrossRef]
  72. Li, S.; Yang, J.; Liu, X.; Guo, R.; Zhang, R. circITGA7 Functions as an Oncogene by Sponging miR-198 and Upregulating FGFR1 Expression in Thyroid Cancer. BioMed Res. Int. 2020, 2020, 1–8. [Google Scholar] [CrossRef]
  73. Chen, F.; Feng, Z.; Zhu, J.; Liu, P.; Yang, C.; Huang, R.; Deng, Z. Emerging roles of circRNA_NEK6 targeting miR-370-3p in the proliferation and invasion of thyroid cancer via Wnt signaling pathway. Cancer Biol. Ther. 2018, 19, 1139–1152. [Google Scholar] [CrossRef]
  74. Abe, N.; Matsumoto, K.; Nishihara, M.; Nakano, Y.; Shibata, A.; Maruyama, H.; Shuto, S.; Matsuda, A.; Yoshida, M.; Ito, Y.; et al. Rolling Circle Translation of Circular RNA in Living Human Cells. Sci. Rep. 2015, 5, 16435. [Google Scholar] [CrossRef]
  75. Chen, C.-Y.; Sarnow, P. Initiation of Protein Synthesis by the Eukaryotic Translational Apparatus on Circular RNAs. Science 1995, 268, 415–417. [Google Scholar] [CrossRef]
  76. Zhang, M.; Huang, N.; Yang, X.; Luo, J.; Yan, S.; Xiao, F.; Chen, W.; Gao, X.; Zhao, K.; Zhou, H.; et al. A novel protein encoded by the circular form of the SHPRH gene suppresses glioma tumorigenesis. Oncogene 2018, 37, 1805–1814. [Google Scholar] [CrossRef]
  77. Peng, L.; Yuan, X.Q.; Li, G.C. The emerging landscape of circular RNA ciRS-7 in cancer (Review). Oncol. Rep. 2015, 33, 2669–2674. [Google Scholar] [CrossRef]
  78. Han, J.Y.; Guo, S.; Wei, N.; Xue, R.; Li, W.; Dong, G.; Li, J.; Tian, X.; Chen, C.; Qiu, S.; et al. ciRS-7 Promotes the Proliferation and Migration of Papillary Thyroid Cancer by Negatively Regulating the miR-7/Epidermal Growth Factor Receptor Axis. Biomed. Res. Int. 2020, 2020, 9875636. [Google Scholar] [CrossRef]
  79. Abdelmohsen, K.; Panda, A.C.; Munk, R.; Grammatikakis, I.; Dudekula, D.B.; De, S.; Kim, J.; Noh, J.H.; Kim, K.M.; Martindale, J.L.; et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017, 14, 361–369. [Google Scholar] [CrossRef]
  80. Ghosal, S.; Das, S.; Sen, R.; Basak, P.; Chakrabarti, J. Circ2Traits: A comprehensive database for circular RNA potentially associated with disease and traits. Front. Genet. 2013, 4, 283. [Google Scholar] [CrossRef]
  81. Dudekula, D.B.; Panda, A.C.; Grammatikakis, I.; De, S.; Abdelmohsen, K.; Gorospe, M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016, 13, 34–42. [Google Scholar] [CrossRef] [PubMed]
  82. Liu, Y.-C.; Li, J.-R.; Sun, C.-H.; Andrews, E.; Chao, R.-F.; Lin, F.-M.; Weng, S.-L.; Hsu, S.-D.; Huang, C.-C.; Cheng, C.; et al. CircNet: A database of circular RNAs derived from transcriptome sequencing data. Nucleic Acids Res. 2015, 44, D209–D215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Chen, X.; Han, P.; Zhou, T.; Guo, X.; Song, X.; Li, Y. circRNADb: A comprehensive database for human circular RNAs with protein-coding annotations. Sci. Rep. 2016, 6, 34985. [Google Scholar] [CrossRef] [PubMed]
  84. Xia, S.; Feng, J.; Chen, K.; Ma, Y.; Gong, J.; Cai, F.; Jin, Y.; Gao, Y.; Xia, L.; Chang, H.; et al. CSCD: A database for cancer-specific circular RNAs. Nucleic Acids Res. 2017, 46, D925–D929. [Google Scholar] [CrossRef]
  85. Xia, S.; Feng, J.; Lei, L.; Hu, J.; Xia, L.; Wang, J.; Xiang, Y.; Liu, L.; Zhong, S.; Han, L.; et al. Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes. Brief. Bioinform. 2016, 18, 984–992. [Google Scholar] [CrossRef]
  86. Dong, R.; Ma, X.K.; Li, G.W.; Yang, L. CIRCpedia v2: An Updated Database for Comprehensive Circular RNA Annotation and Expression Comparison. Genom. Proteom. Bioinform. 2018, 16, 226–233. [Google Scholar] [CrossRef]
  87. Yao, D.; Zhang, L.; Zheng, M.; Sun, X.; Lu, Y.; Liu, P. Circ2Disease: A manually curated database of experimentally validated circRNAs in human disease. Sci. Rep. 2018, 8, 1–6. [Google Scholar] [CrossRef]
  88. Zhao, Z.; Wang, K.; Wu, F.; Wang, W.; Zhang, K.; Hu, H.; Liu, Y.; Jiang, T. circRNA disease: A manually curated database of experimentally supported circRNA-disease associations. Cell Death Dis. 2018, 9, 1–2. [Google Scholar] [CrossRef]
  89. Meng, X.; Hu, D.; Zhang, P.; Chen, Q.; Chen, M. CircFunBase: A database for functional circular RNAs. Database 2019, 2019, baz003. [Google Scholar] [CrossRef]
  90. 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, D92–D97. [Google Scholar] [CrossRef]
  91. Gokool, A.; Anwar, F.; Voineagu, I. The Landscape of Circular RNA Expression in the Human Brain. Biol. Psychiatry 2020, 87, 294–304. [Google Scholar] [CrossRef]
  92. Liu, M.; Wang, Q.; Shen, J.; Yang, B.B.; Ding, X. Circbank: A comprehensive database for circRNA with standard nomenclature. RNA Biol. 2019, 16, 899–905. [Google Scholar] [CrossRef]
  93. Meng, X.; Chen, Q.; Zhang, P.; Chen, M. CircPro: An integrated tool for the identification of circRNAs with protein-coding potential. Bioinformatics 2017, 33, 3314–3316. [Google Scholar] [CrossRef]
  94. Cardenas, J.; Balaji, U.; Gu, J. Cerina: Systematic circRNA functional annotation based on integrative analysis of ceRNA interactions. Sci. Rep. 2020, 10, 1–14. [Google Scholar] [CrossRef]
  95. Wu, W.; Ji, P.; Zhao, F. CircAtlas: An integrated resource of one million highly accurate circular RNAs from 1070 vertebrate transcriptomes. Genome Biol. 2020, 21, 1–14. [Google Scholar] [CrossRef]
  96. Zheng, L.-L.; Ling-Ling, Z.; Wu, J.; Sun, W.-J.; Liu, S.; Wang, Z.-L.; Zhou, H.; Yang, J.-H.; Qu, L.-H. deepBase v2.0: Identification, expression, evolution and function of small RNAs, LncRNAs and circular RNAs from deep-sequencing data. Nucleic Acids Res. 2015, 44, D196–D202. [Google Scholar] [CrossRef]
  97. Rophina, M.; Sharma, D.; Poojary, M.; Scaria, V. Circad: A comprehensive manually curated resource of circular RNA associated with diseases. Database 2020, 2020, baaa019. [Google Scholar] [CrossRef]
  98. Lan, W.; Zhu, M.; Chen, Q.; Chen, B.; Liu, J.; Li, M.; Chen, Y.-P.P. CircR2Cancer: A manually curated database of associations between circRNAs and cancers. Database 2020, 2020, baaa085. [Google Scholar] [CrossRef]
  99. Li, S.; Li, Y.; Chen, B.; Zhao, J.; Yu, S.; Tang, Y.; Zheng, Q.; Li, Y.; Wang, P.; He, X.; et al. exoRBase: A database of circRNA, lncRNA and mRNA in human blood exosomes. Nucleic Acids Res. 2017, 46, D106–D112. [Google Scholar] [CrossRef]
  100. Vo, J.N.; Cieslik, M.; Zhang, Y.; Shukla, S.; Xiao, L.; Zhang, Y.; Wu, Y.-M.; Dhanasekaran, S.M.; Engelke, C.G.; Cao, X.; et al. The Landscape of Circular RNA in Cancer. Cell 2019, 176, 869–881.e13. [Google Scholar] [CrossRef]
  101. Xin, R.; Gao, Y.; Wang, R.; Kadash-Edmondson, K.E.; Liu, B.; Wang, Y.; Lin, L.; Xing, Y. isoCirc catalogs full-length circular RNA isoforms in human transcriptomes. Nat. Commun. 2021, 12, 1–11. [Google Scholar] [CrossRef]
  102. Wu, J.; Li, Y.; Wang, C.; Cui, Y.; Xu, T.; Wang, C.; Wang, X.; Sha, J.; Jiang, B.; Wang, K.; et al. CircAST: Full-length Assembly and Quantification of Alternatively Spliced Isoforms in Circular RNAs. Genom. Proteom. Bioinform. 2019, 17, 522–534. [Google Scholar] [CrossRef]
  103. Ruan, H.; Xiang, Y.; Ko, J.; Li, S.; Jing, Y.; Zhu, X.; Han, L. Comprehensive characterization of circular RNAs in ~1000 human cancer cell lines. Genome Med. 2019. 11, 55. [CrossRef]
  104. Tang, Z.; Li, X.; Zhao, J.; Qian, F.; Feng, C.; Li, Y.; Zhang, J.; Jiang, Y.; Yang, Y.; Wang, Q.; et al. TRCirc: A resource for transcriptional regulation information of circRNAs. Brief. Bioinform. 2018, 20, 2327–2333. [Google Scholar] [CrossRef]
  105. Wang, P.; Li, X.; Gao, Y.; Guo, Q.; Wang, Y.; Fang, Y.; Ma, X.; Zhi, H.; Zhou, D.; Shen, W.; et al. LncACTdb 2.0: An updated database of experimentally supported ceRNA interactions curated from low- and high-throughput experiments. Nucleic Acids Res. 2018, 47, D121–D127. [Google Scholar] [CrossRef]
  106. Zhang, W.; Yao, G.; Wang, J.; Yang, M.; Wang, J.; Zhang, H.; Li, W. ncRPheno: A comprehensive database platform for identification and validation of disease related noncoding RNAs. RNA Biol. 2020, 17, 943–955. [Google Scholar] [CrossRef]
  107. Teng, X.; Chen, X.; Xue, H.; Tang, Y.; Zhang, P.; Kang, Q.; Hao, Y.; Chen, R.; Zhao, Y.; He, S. NPInter v4.0: An integrated database of ncRNA interactions. Nucleic Acids Res. 2019, 48, D160–D165. [Google Scholar] [CrossRef]
  108. Toraih, E.A.; Sedhom, J.A.; Haidari, M.; Fawzy, M.S. Applications of noncoding RNAs in renal cancer patients. In Clinical Applications of Non-Coding RNAs in Cancer; Elsevier: Amsterdam, The Netherlands, 2022; pp. 211–284. [Google Scholar] [CrossRef]
  109. Li, Z.; Huang, X.; Liu, A.; Xu, J.; Lai, J.; Guan, H.; Ma, J. Circ_PSD3 promotes the progression of papillary thyroid carcinoma via the miR-637/HEMGN axis. Life Sci. 2021, 264, 118622. [Google Scholar] [CrossRef]
  110. Yao, Y.; Chen, X.; Yang, H.; Chen, W.; Qian, Y.; Yan, Z.; Liao, T.; Yao, W.; Wu, W.; Yu, T.; et al. Hsa_circ_0058124 promotes papillary thyroid cancer tumorigenesis and invasiveness through the NOTCH3/GATAD2A axis. J. Exp. Clin. Cancer Res. 2019, 38, 318. [Google Scholar] [CrossRef]
  111. Zhou, G.-K.; Zhang, G.-Y.; Yuan, Z.-N.; Pei, R.; Liu, D.-M. Has_circ_0008274 promotes cell proliferation and invasion involving AMPK/mTOR signaling pathway in papillary thyroid carcinoma. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 8772–8780. [Google Scholar]
  112. Ma, J.; Kan, Z. Circular RNA circ_0008274 enhances the malignant progression of papillary thyroid carcinoma via modulating solute carrier family 7 member 11 by sponging miR-154-3p. Endocr. J. 2021, 68, 543–552. [Google Scholar] [CrossRef] [PubMed]
  113. Zhang, H.; Ma, X.-P.; Li, X.; Deng, F.-S. Circular RNA circ_0067934 exhaustion expedites cell apoptosis and represses cell proliferation, migration and invasion in thyroid cancer via sponging miR-1304 and regulating CXCR1 expression. Eur. Rev. Med. Pharmacol. Sci. 2019, 23, 10851–10866. [Google Scholar] [PubMed]
  114. Lan, X.; Cao, J.; Xu, J.; Chen, C.; Zheng, C.; Wang, J.; Zhu, X.; Zhu, X.; Ge, M. Decreased expression of hsa_circ_0137287 predicts aggressive clinicopathologic characteristics in papillary thyroid carcinoma. J. Clin. Lab. Anal. 2018, 32, e22573. [Google Scholar] [CrossRef] [PubMed]
  115. Yu, W.; Ma, B.; Zhao, W.; Liu, J.; Yu, H.; Tian, Z.; Fan, Z.; Han, H. The combination of circRNA-UMAD1 and Galectin-3 in peripheral circulation is a co-biomarker for predicting lymph node metastasis of thyroid carcinoma. Am. J. Transl. Res. 2020, 12, 5399–5415. [Google Scholar] [PubMed]
  116. Jin, X.; Wang, Z.; Pang, W.; Zhou, J.; Liang, Y.; Yang, J.; Yang, L.; Zhang, Q. Upregulated hsa_circ_0004458 Contributes to Progression of Papillary Thyroid Carcinoma by Inhibition of miR-885-5p and Activation of RAC1. Med. Sci. Monit. 2018, 24, 5488–5500. [Google Scholar] [CrossRef] [PubMed]
  117. Wang, Y.; Zong, H.; Zhou, H. Circular RNA circ_0062389 modulates papillary thyroid carcinoma progression via the miR-1179/high mobility group box 1 axis. Bioengineered 2021, 12, 1484–1494. [Google Scholar] [CrossRef] [PubMed]
  118. Qi, Y.; He, J.; Zhang, Y.; Wang, L.; Yu, Y.; Yao, B.; Tian, Z. Circular RNA hsa_circ_0001666 sponges miR-330-5p, miR-193a-5p and miR-326, and promotes papillary thyroid carcinoma progression via upregulation of ETV4. Oncol. Rep. 2021, 45, 1–12. [Google Scholar] [CrossRef]
  119. Liu, F.; Zhang, J.; Qin, L.; Yang, Z.; Xiong, J.; Zhang, Y.; Li, R.; Li, S.; Wang, H.; Yu, B.; et al. Circular RNA EIF6 (Hsa_circ_0060060) sponges miR-144-3p to promote the cisplatin-resistance of human thyroid carcinoma cells by autophagy regulation. Aging 2018, 10, 3806–3820. [Google Scholar] [CrossRef]
  120. Rong, D.; Sun, H.; Li, Z.; Liu, S.; Dong, C.; Fu, K.; Tang, W.; Cao, H. An emerging function of circRNA-miRNAs-mRNA axis in human diseases. Oncotarget 2017, 8, 73271–73281. [Google Scholar] [CrossRef]
  121. Wang, H.; Yan, X.; Zhang, H.; Zhan, X. CircRNA circ_0067934 Overexpression Correlates with Poor Prognosis and Promotes Thyroid Carcinoma Progression. Med. Sci. Monit. 2019, 25, 1342–1349. [Google Scholar] [CrossRef]
  122. Liao, J.-Y.; Wu, J.; Wang, Y.-J.; He, J.-H.; Deng, W.-X.; Hu, K.; Zhang, Y.-C.; Zhang, Y.; Yan, H.; Wang, D.-L.; et al. Deep sequencing reveals a global reprogramming of lncRNA transcriptome during EMT. Biochim. Biophys. Acta Mol. Cell Res. 2017, 1864, 1703–1713. [Google Scholar] [CrossRef]
  123. Li, Z.; Xu, J.; Guan, H.; Lai, J.; Yang, X.; Ma, J. Circ_0059354 aggravates the progression of papillary thyroid carcinoma by elevating ARFGEF1 through sponging miR-766-3p. J. Endocrinol. Investig. 2021, 45, 825–836. [Google Scholar] [CrossRef]
  124. Cai, Z.; Li, H. Circular RNAs and Bladder Cancer. OncoTargets Ther. 2020, 13, 9573–9586. [Google Scholar] [CrossRef]
  125. Zeng, Y.; Xu, Y.; Shu, R.; Sun, L.; Tian, Y.; Shi, C.; Zheng, Z.; Wang, K.; Luo, H. Altered expression profiles of circular RNA in colorectal cancer tissues from patients with lung metastasis. Int. J. Mol. Med. 2017, 40, 1818–1828. [Google Scholar] [CrossRef]
  126. Schweppe, R.E.; Klopper, J.P.; Korch, C.; Pugazhenthi, U.; Benezra, M.; Knauf, J.A.; Fagin, J.A.; Marlow, L.A.; Copland, J.A.; Smallridge, R.C.; et al. Deoxyribonucleic Acid Profiling Analysis of 40 Human Thyroid Cancer Cell Lines Reveals Cross-Contamination Resulting in Cell Line Redundancy and Misidentification. J. Clin. Endocrinol. Metab. 2008, 93, 4331–4341. [Google Scholar] [CrossRef]
  127. Fan, Y.X.; Shi, H.Y.; Hu, Y.L.; Jin, X.L. Circ_0000144 facilitates the progression of thyroid cancer via the miR-217/AKT3 pathway. J. Gene Med. 2020, 22, e3269. [Google Scholar] [CrossRef]
  128. Landa, I.; Pozdeyev, N.; Korch, C.; Marlow, L.A.; Smallridge, R.C.; Copland, J.A.; Henderson, Y.C.; Lai, S.Y.; Clayman, G.L.; Onoda, N.; et al. Comprehensive Genetic Characterization of Human Thyroid Cancer Cell Lines: A Validated Panel for Preclinical Studies. Clin. Cancer Res. 2019, 25, 3141–3151. [Google Scholar] [CrossRef]
  129. Fogh, J. Human Tumor Lines for Cancer Research. Cancer Investig. 1986, 4, 157–184. [Google Scholar] [CrossRef]
  130. Peng, N.; Shi, L.; Zhang, Q.; Hu, Y.; Wang, N.; Ye, H. Microarray profiling of circular RNAs in human papillary thyroid carcinoma. PLoS ONE 2017, 12, e0170287. [Google Scholar] [CrossRef]
  131. Liu, L.; Yan, C.; Tao, S.; Wang, H. Circ_0058124 Aggravates the Progression of Papillary Thyroid Carcinoma by Activating LMO4 Expression via Targeting miR-370-3p. Cancer Manag. Res. 2020, 12, 9459–9470. [Google Scholar] [CrossRef]
  132. Gui, X.; Li, Y.; Zhang, X.; Su, K.; Cao, W. Circ_LDLR promoted the development of papillary thyroid carcinoma via regulating miR-195-5p/LIPH axis. Cancer Cell Int. 2020, 20, 1–14. [Google Scholar] [CrossRef]
  133. Tan, X.; Zhao, J.; Lou, J.; Zheng, W.; Wang, P. Hsa_circ_0058129 regulates papillary thyroid cancer development via miR-873-5p/follistatin-like 1 axis. J. Clin. Lab. Anal. 2022, 36, e24401. [Google Scholar] [CrossRef]
  134. Cui, W.; Xue, J. Circular RNA DOCK1 downregulates microRNA -124 to induce the growth of human thyroid cancer cell lines. BioFactors 2020, 46, 591–599. [Google Scholar] [CrossRef]
  135. Xue, C.; Cheng, Y.; Wu, J.; Ke, K.; Miao, C.; Chen, E.; Zhang, L. Circular RNA CircPRMT5 Accelerates Proliferation and Invasion of Papillary Thyroid Cancer Through Regulation of miR-30c/E2F3 Axis. Cancer Manag. Res. 2020, 12, 3285–3291. [Google Scholar] [CrossRef]
  136. Ye, M.; Hou, H.; Shen, M.; Dong, S.; Zhang, T. Circular RNA circFOXM1 Plays a Role in Papillary Thyroid Carcinoma by Sponging miR-1179 and Regulating HMGB1 Expression. Mol. Ther.-Nucleic Acids 2019, 19, 741–750. [Google Scholar] [CrossRef]
  137. Xia, F.; Chen, Y.; Jiang, B.; Bai, N.; Li, X. Hsa_circ_0011385 accelerates the progression of thyroid cancer by targeting miR-361-3p. Cancer Cell Int. 2020, 20, 1–14. [Google Scholar] [CrossRef]
  138. Liu, W.; Zhao, J.; Jin, M.; Zhou, M. circRAPGEF5 Contributes to Papillary Thyroid Proliferation and Metastatis by Regulation miR-198/FGFR1. Mol. Ther. Nucleic Acids 2019, 14, 609–616. [Google Scholar] [CrossRef]
  139. Zheng, F.-B.; Chen, D.; Ding, Y.-Y.; Wang, S.-R.; Shi, D.-D.; Zhu, Z.-P. Circular RNA circ_0103552 promotes the invasion and migration of thyroid carcinoma cells by sponging miR-127. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 2572–2578. [Google Scholar]
  140. Guo, M.; Sun, Y.; Ding, J.; Li, Y.; Yang, S.; Zhao, Y.; Jin, X.; Li, S.-S. Circular RNA profiling reveals a potential role of hsa_circ_IPCEF1 in papillary thyroid carcinoma. Mol. Med. Rep. 2021, 24, 1–15. [Google Scholar] [CrossRef]
  141. Xiong, H.; Yu, J.; Jia, G.; Su, Y.; Zhang, J.; Xu, Q.; Sun, X. Emerging roles of circUBAP2 targeting miR-370-3p in proliferation, apoptosis, and invasion of papillary thyroid cancer cells. Hum. Cell 2021, 34, 1866–1877. [Google Scholar] [CrossRef] [PubMed]
  142. Zeng, L.; Yuan, S.; Zhou, P.; Gong, J.; Kong, X.; Wu, M. Circular RNA Pvt1 oncogene (CircPVT1) promotes the progression of papillary thyroid carcinoma by activating the Wnt/β-catenin signaling pathway and modulating the ratio of microRNA-195 (miR-195) to vascular endothelial growth factor A (VEGFA) expression. Bioengineered 2021, 12, 11795–11810. [Google Scholar] [CrossRef] [PubMed]
  143. Li, Y.; Qin, J.; He, Z.; Cui, G.; Zhang, K.; Wu, B. Knockdown of circPUM1 impedes cell growth, metastasis and glycolysis of papillary thyroid cancer via enhancing MAPK1 expression by serving as the sponge of miR-21-5p. Genes Genom. 2021, 43, 141–150. [Google Scholar] [CrossRef] [PubMed]
  144. Zheng, H.; Fu, Q.; Ma, K.; Shi, S.; Fu, Y. Circ_0079558 promotes papillary thyroid cancer progression by binding to miR-26b-5p to activate MET/AKT signaling. Endocr. J. 2021, 68, 1247–1266. [Google Scholar] [CrossRef]
  145. Mao, Y.; Huo, Y.; Li, J.; Zhao, Y.; Wang, Y.; Sun, L.; Kang, Z. circRPS28 (hsa_circ_0049055) is a novel contributor for papillary thyroid carcinoma by regulating cell growth and motility via functioning as ceRNA for miR-345-5p to regulate frizzled family receptor 8 (FZD8). Endocr. J. 2021, 68, 1267–1281. [Google Scholar] [CrossRef]
  146. Dong, L.-P.; Chen, L.-Y.; Bai, B.; Qi, X.-F.; Liu, J.-N.; Qin, S. circ_0067934 promotes the progression of papillary thyroid carcinoma cells through miR-1301-3p/HMGB1 axis. Neoplasma 2022, 69, 1–15. [Google Scholar] [CrossRef]
  147. Shu, T.; Yang, L.; Sun, L.; Lu, J.; Zhan, X. CircHIPK3 Promotes Thyroid Cancer Tumorigenesis and Invasion through the Mirna-338-3p/RAB23 Axis. Med. Princ. Pr. 2020, 30, 412–420. [Google Scholar] [CrossRef]
  148. Ding, W.; Shi, Y.; Zhang, H. Circular RNA circNEURL4 inhibits cell proliferation and invasion of papillary thyroid carcinoma by sponging miR-1278 and regulating LATS1 expression. Am. J. Transl. Res. 2021, 13, 5911–5927. [Google Scholar]
  149. Zhang, D.; Tao, L.; Xu, N.; Lu, X.; Wang, J.; He, G.; Chu, J. CircRNA circTIAM1 promotes papillary thyroid cancer progression through the miR-646/HNRNPA1 signaling pathway. Cell Death Discov. 2022, 8, 21. [Google Scholar] [CrossRef]
  150. Wang, T.; Huang, Y. Downregulation of hsa_circ_0001681 suppresses epithelial-mesenchymal transition in thyroid carcinoma via targeting to miR-942-5p/TWIST1 signaling pathway. J. Bioenerg. Biomembr. 2021, 53, 609–620. [Google Scholar] [CrossRef]
  151. Keller, S.; Sanderson, M.P.; Stoeck, A.; Altevogt, P. Exosomes: From biogenesis and secretion to biological function. Immunol. Lett. 2006, 107, 102–108. [Google Scholar] [CrossRef]
  152. Tuo, B.; Chen, Z.; Dang, Q.; Chen, C.; Zhang, H.; Hu, S.; Sun, Z. Roles of exosomal circRNAs in tumour immunity and cancer progression. Cell Death Dis. 2022, 13, 1–9. [Google Scholar] [CrossRef]
  153. Lin, Q.; Qi, Q.; Hou, S.; Chen, Z.; Jiang, N.; Zhang, L.; Lin, C. Exosomal circular RNA hsa_circ_007293 promotes proliferation, migration, invasion, and epithelial-mesenchymal transition of papillary thyroid carcinoma cells through regulation of the microRNA-653-5p/paired box 6 axis. Bioengineered 2021, 12, 10136–10149. [Google Scholar] [CrossRef]
  154. Liu, Y.; Xia, L.; Dong, L.; Wang, J.; Xiao, Q.; Yu, X.; Zhu, H. CircHIPK3 Promotes Gemcitabine (GEM) Resistance in Pancreatic Cancer Cells by Sponging miR-330-5p and Targets RASSF1. Cancer Manag. Res. 2020, 12, 921–929. [Google Scholar] [CrossRef]
  155. Yuan, W.; Zhou, R.; Wang, J.; Han, J.; Yang, X.; Yu, H.; Lu, H.; Zhang, X.; Li, P.; Tao, J.; et al. Circular RNA Cdr1as sensitizes bladder cancer to cisplatin by upregulating APAF1 expression through miR-1270 inhibition. Mol. Oncol. 2019, 13, 1559–1576. [Google Scholar] [CrossRef] [Green Version]
  156. Sang, Y.; Chen, B.; Song, X.; Li, Y.; Liang, Y.; Han, D.; Zhang, N.; Zhang, H.; Liu, Y.; Chen, T.; et al. circRNA_0025202 Regulates Tamoxifen Sensitivity and Tumor Progression via Regulating the miR-182-5p/FOXO3a Axis in Breast Cancer. Mol. Ther. 2019, 27, 1638–1652. [Google Scholar] [CrossRef]
  157. Suzuki, H.; Zuo, Y.; Wang, J.; Zhang, M.Q.; Malhotra, A.; Mayeda, A. Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing. Nucleic Acids Res. 2006, 34, e63. [Google Scholar] [CrossRef]
  158. Panda, A.; De, S.; Grammatikakis, I.; Munk, R.; Yang, X.; Piao, Y.; Dudekula, D.B.; Abdelmohsen, K.; Gorospe, M. High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs. Nucleic Acids Res. 2017, 45, e116. [Google Scholar] [CrossRef]
  159. Conn, V.; Conn, S.J. SplintQuant: A method for accurately quantifying circular RNA transcript abundance without reverse transcription bias. RNA 2019, 25, 1202–1210. [Google Scholar] [CrossRef]
  160. Hansen, T.B. Improved circRNA Identification by Combining Prediction Algorithms. Front. Cell Dev. Biol. 2018, 6, 20. [Google Scholar] [CrossRef]
Figure 1. The biogenesis of circRNAs. There are three major classes of circRNAs: Class (1) Exonic circular RNA (ecircRNA): Initially, intron No. 1 is removed, and the 5′ splice site of exon No. 2 is brought close to 3′ splice site of Exon No. 1, creating ecircRNA that contains multiple exons. Exons can also skip splicing; exon No. 1 can also link with exon No. 3. Class (2) Circular intronic RNA (ciRNA): Reverse complementary sequences of lariat intron spliced from pre-mRNA can pair to produce a closed-loop structure termed ciRNA. Class (3) Exon–intron circRNAs (EIciRNAs): Intron No. 4 can be retained with Exon No. 4 and Exon No. 5 to form an EIciRNA (Created by Biorender.com).
Figure 1. The biogenesis of circRNAs. There are three major classes of circRNAs: Class (1) Exonic circular RNA (ecircRNA): Initially, intron No. 1 is removed, and the 5′ splice site of exon No. 2 is brought close to 3′ splice site of Exon No. 1, creating ecircRNA that contains multiple exons. Exons can also skip splicing; exon No. 1 can also link with exon No. 3. Class (2) Circular intronic RNA (ciRNA): Reverse complementary sequences of lariat intron spliced from pre-mRNA can pair to produce a closed-loop structure termed ciRNA. Class (3) Exon–intron circRNAs (EIciRNAs): Intron No. 4 can be retained with Exon No. 4 and Exon No. 5 to form an EIciRNA (Created by Biorender.com).
Cancers 14 04728 g001
Figure 2. The putative functions of circRNAs. (1) The nuclear type can act as regulators of their host genes’ transcription, either by interacting with U1 small nuclear ribonucleoprotein (U1 snRNP) and enhancing RNA polymerase II (Pol II) function, or by recruiting methyl cytosine dioxygenase (TET1) to the promoter sequences [51]. (2) Epigenetic regulation, including DNA and histone modifications [51]. (3) In cytoplasm, circRNAs can act as microRNA sponges, leading to a disrupted competing endogenous RNA (ceRNA) network [16,50]. (4) CircRNAs can function as RNA binding protein (RBP) decoys or otherwise can act as modulators for specific RBPs’ half-life, promoting or reducing their degradation by proteasomes. Enabling the colocalization of enzymes to their substrates facilitates reaction kinetics. CircRNAs can also act as protein scaffolds. (5) CircRNAs with internal ribosome entry site (IRES) elements and the initial methionine code (AUG) for translation sites can be translated into proteins through a CAP-independent mechanism, facilitated by the presence of methyl adenosine (m6A) and other translation-specific factors. (6) CircRNAs have been implicated in several posttranslational modifications, including phosphorylation (P), ubiquitylation (Ub) and acetylation (Ac). (7) Finally, circRNAs can be secreted in extracellular vesicles regulating intercellular communication [56].
Figure 2. The putative functions of circRNAs. (1) The nuclear type can act as regulators of their host genes’ transcription, either by interacting with U1 small nuclear ribonucleoprotein (U1 snRNP) and enhancing RNA polymerase II (Pol II) function, or by recruiting methyl cytosine dioxygenase (TET1) to the promoter sequences [51]. (2) Epigenetic regulation, including DNA and histone modifications [51]. (3) In cytoplasm, circRNAs can act as microRNA sponges, leading to a disrupted competing endogenous RNA (ceRNA) network [16,50]. (4) CircRNAs can function as RNA binding protein (RBP) decoys or otherwise can act as modulators for specific RBPs’ half-life, promoting or reducing their degradation by proteasomes. Enabling the colocalization of enzymes to their substrates facilitates reaction kinetics. CircRNAs can also act as protein scaffolds. (5) CircRNAs with internal ribosome entry site (IRES) elements and the initial methionine code (AUG) for translation sites can be translated into proteins through a CAP-independent mechanism, facilitated by the presence of methyl adenosine (m6A) and other translation-specific factors. (6) CircRNAs have been implicated in several posttranslational modifications, including phosphorylation (P), ubiquitylation (Ub) and acetylation (Ac). (7) Finally, circRNAs can be secreted in extracellular vesicles regulating intercellular communication [56].
Cancers 14 04728 g002
Figure 3. CircRNAs serve as miRNA sponges.
Figure 3. CircRNAs serve as miRNA sponges.
Cancers 14 04728 g003
Figure 4. In vitro cell lines were used to study circRNA expression [126,127,128,129].
Figure 4. In vitro cell lines were used to study circRNA expression [126,127,128,129].
Cancers 14 04728 g004
Figure 5. CircRNAs and chemotherapy drugs in thyroid cancer. (A) The role of overexpressed “circEIF6”, its target miR-144-3p, and downstream “transforming growth factor-α (TGF-α)” in cisplatin resistance in papillary and anaplastic thyroid cancer cell lines. CircEIF6 sponges miR-144-3p to promote the cisplatin resistance of human TC cells by autophagy regulation. TGF-α activation accelerates the carcinogenesis process and cisplatin resistance. (B) Knockout of “circEIF6” reverses the processes mentioned above, sensitizing the papillary and anaplastic thyroid cancer cell lines to cisplatin (Created by Biorender.com).
Figure 5. CircRNAs and chemotherapy drugs in thyroid cancer. (A) The role of overexpressed “circEIF6”, its target miR-144-3p, and downstream “transforming growth factor-α (TGF-α)” in cisplatin resistance in papillary and anaplastic thyroid cancer cell lines. CircEIF6 sponges miR-144-3p to promote the cisplatin resistance of human TC cells by autophagy regulation. TGF-α activation accelerates the carcinogenesis process and cisplatin resistance. (B) Knockout of “circEIF6” reverses the processes mentioned above, sensitizing the papillary and anaplastic thyroid cancer cell lines to cisplatin (Created by Biorender.com).
Cancers 14 04728 g005
Figure 6. CircRNAs regulate various hallmarks of cancer.
Figure 6. CircRNAs regulate various hallmarks of cancer.
Cancers 14 04728 g006
Figure 7. Pathway enrichment analysis for overexpressed genes caused by upregulated circRNAs. Data source: “PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification System (http://pantherdb.org (accessed on 23 June 2022))”.
Figure 7. Pathway enrichment analysis for overexpressed genes caused by upregulated circRNAs. Data source: “PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification System (http://pantherdb.org (accessed on 23 June 2022))”.
Cancers 14 04728 g007
Figure 8. Microarray transcriptomic signature of circRNAs in human papillary thyroid carcinoma. (A) Characteristics of six PTC patients were retrieved from the Gene Expression Omnibus database (GSE93522). PTC was compared with their paired contralateral normal thyroid tissues (B) Boxplots showing the distribution of the values of the selected samples. Median-centered values are indicated after log transform and quantile normalization. (C) Principal component analysis for data exploration. The geometrical projection for samples is based on the similarities and differences in the expression of 2895 circRNAs. Samples are plotted across two coordinates; axis 1 explained 52.1% of the variance, whereas axis 2 demonstrated 17% of the data variability. (D) Volcano plot displays statistical significance (−log10 p-value) versus magnitude of change (log2 fold change) in PTC compared with controls. Each dot represents a circRNA. The plot shows 137 significantly differentially expressed circRNAs (DEC): 115 upregulated (red dots) and 22 downregulated (blue dots). Analysis was performed using GEO2R software, with the p-value threshold at <0.05 and log fold changes (FC) at >1. (E) Heatmap for the 12 thyroid tissue specimens using the 137 differentially expressed circRNAs. There was an upregulation pattern of genes in patient samples. Clear differentiation of four PTC tissues and incomplete separation of two others are shown. (F) Principal component analysis to visualize how samples are related to each other. By using the deregulated genes, a complete demarcation was observed between PTC and normal tissues. Axes 1 and 2 explain 66.6% and 9.7% of the variability, respectively. (G) The top-up and downregulated circRNAs in PTC tissues. Derived protein-coding genes from which circRNAs were formed are shown. Fold change values are depicted at the end of the bars. (H) Frequency of tested circRNAs in the microarray according to their chromosomal localization. Y axis is log10 transformed, and the X axis represents the human chromosomes with removed missing chromosomes (chromosomes 18, 21, and Y). (I) Chromosomal ideograms represent the differentially expressed circRNAs annotated with lines in color at specific base-pair locations: blue for upregulated DECs and green for downregulated DECs.
Figure 8. Microarray transcriptomic signature of circRNAs in human papillary thyroid carcinoma. (A) Characteristics of six PTC patients were retrieved from the Gene Expression Omnibus database (GSE93522). PTC was compared with their paired contralateral normal thyroid tissues (B) Boxplots showing the distribution of the values of the selected samples. Median-centered values are indicated after log transform and quantile normalization. (C) Principal component analysis for data exploration. The geometrical projection for samples is based on the similarities and differences in the expression of 2895 circRNAs. Samples are plotted across two coordinates; axis 1 explained 52.1% of the variance, whereas axis 2 demonstrated 17% of the data variability. (D) Volcano plot displays statistical significance (−log10 p-value) versus magnitude of change (log2 fold change) in PTC compared with controls. Each dot represents a circRNA. The plot shows 137 significantly differentially expressed circRNAs (DEC): 115 upregulated (red dots) and 22 downregulated (blue dots). Analysis was performed using GEO2R software, with the p-value threshold at <0.05 and log fold changes (FC) at >1. (E) Heatmap for the 12 thyroid tissue specimens using the 137 differentially expressed circRNAs. There was an upregulation pattern of genes in patient samples. Clear differentiation of four PTC tissues and incomplete separation of two others are shown. (F) Principal component analysis to visualize how samples are related to each other. By using the deregulated genes, a complete demarcation was observed between PTC and normal tissues. Axes 1 and 2 explain 66.6% and 9.7% of the variability, respectively. (G) The top-up and downregulated circRNAs in PTC tissues. Derived protein-coding genes from which circRNAs were formed are shown. Fold change values are depicted at the end of the bars. (H) Frequency of tested circRNAs in the microarray according to their chromosomal localization. Y axis is log10 transformed, and the X axis represents the human chromosomes with removed missing chromosomes (chromosomes 18, 21, and Y). (I) Chromosomal ideograms represent the differentially expressed circRNAs annotated with lines in color at specific base-pair locations: blue for upregulated DECs and green for downregulated DECs.
Cancers 14 04728 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Toraih, E.A.; Hussein, M.H.; Fawzy, M.S.; Kandil, E. A Review and In Silico Analysis of Tissue and Exosomal Circular RNAs: Opportunities and Challenges in Thyroid Cancer. Cancers 2022, 14, 4728. https://doi.org/10.3390/cancers14194728

AMA Style

Toraih EA, Hussein MH, Fawzy MS, Kandil E. A Review and In Silico Analysis of Tissue and Exosomal Circular RNAs: Opportunities and Challenges in Thyroid Cancer. Cancers. 2022; 14(19):4728. https://doi.org/10.3390/cancers14194728

Chicago/Turabian Style

Toraih, Eman A., Mohammad H. Hussein, Manal S. Fawzy, and Emad Kandil. 2022. "A Review and In Silico Analysis of Tissue and Exosomal Circular RNAs: Opportunities and Challenges in Thyroid Cancer" Cancers 14, no. 19: 4728. https://doi.org/10.3390/cancers14194728

APA Style

Toraih, E. A., Hussein, M. H., Fawzy, M. S., & Kandil, E. (2022). A Review and In Silico Analysis of Tissue and Exosomal Circular RNAs: Opportunities and Challenges in Thyroid Cancer. Cancers, 14(19), 4728. https://doi.org/10.3390/cancers14194728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop