- Open Access
Microarray analysis and functional prediction of differentially expressed circular RNAs in acquired middle ear cholesteatoma
BioMedical Engineering OnLine volume 20, Article number: 129 (2021)
Middle ear cholesteatoma is characterized by hyper-proliferation of keratinocytes. Circular RNA (circRNA) plays an essential role in the pathogenesis of many proliferative diseases. However, the role of circRNA in the etiopathogenesis of middle ear cholesteatoma is rarely investigated so far. We aimed to investigate the differential expression profiling of circRNAs between acquired middle ear cholesteatoma and normal skin, and to identify potential circRNAs contributing to the etiopathogenesis of middle ear cholesteatoma. Microarray analysis and functional prediction were performed to investigate the circRNA expression profiling between middle ear cholesteatoma and normal skin. Validation of differentially expressed circRNAs was conducted by qRT-PCR. Prediction of m6A modification was also carried out.
Microarray analysis displayed that totally 93 up-regulated and 85 down-regulated circRNAs were identified in middle ear cholesteatoma. Through validation, expressions of hsa_circRNA_104327 and hsa_circRNA_404655 were significantly higher, while hsa_circRNA_000319 was significantly down-regulated in cholesteatoma. GO classification, KEGG pathway, and ceRNA network analyses suggested that these differentially expressed circRNAs might play important roles in the etiopathogenesis of middle ear cholesteatoma. Prediction of m6A modification exhibited that hsa_circRNA_000319 possessed 4 m6A sites with very high confidence, and hsa_circRNA_404655 had 3 m6A sites with high confidence.
Our study revealed that these differentially expressed circRNAs might contribute to the etiopathogenesis of middle ear cholesteatoma. Further researches should be conducted to investigate the exact mechanism of these differentially expressed circRNAs in the etiopathogenesis of middle ear cholesteatoma. Targeting on these circRNAs may provide a new strategy for middle ear cholesteatoma therapy in the future.
Middle ear cholesteatoma is a locally destructive squamous epithelial lesion with clinical features such as hyper-proliferation, migration, aggressiveness, and recurrence. Although middle ear cholesteatoma is pathologically benign, it progressively develops and leads to erosion of adjacent bony structures, resulting in vestibular dysfunction, labyrinthine fistulae, facial paralysis, hearing loss and serious intracranial complications [1, 2]. The incidence of middle ear cholesteatoma has been reported as high as 9.2 per 100,000 populations in Europe, while even higher in Asians . Although a wide spectrum of hypotheses have been put forward, the exact etiopathogenesis of middle ear cholesteatoma is not clear yet. Since there is no effective drug for treating middle ear cholesteatoma, surgical resection is still the only option presently. However, surgical resection is often accompanied by dysfunction of middle and inner ear, which may cause further hearing loss . Due to its recurrence, a few patients have to receive repeated surgeries. Meanwhile, for patients with poor general conditions, such as severe respiratory diseases, cardiovascular and cerebrovascular diseases, surgical resection cannot be performed.
Circular RNA (circRNA) is a special endogenous non-coding RNA with a unique circular structure, derived from exons or introns . MicroRNA (miRNA) regulates the translation of message RNA (mRNA) at the post-transcriptional level via targeting complementary mRNAs and affecting several biological processes, such as cell differentiation, proliferation, and apoptosis . Through binding to their targeted messenger RNA (mRNA) at the 3’-untranslated region, miRNA initiates degradation of mRNAs . CircRNA is rich in miRNA response element (MRE), and functions as efficient miRNA sponges. Acting as a competing endogenous RNA (ceRNA), circRNA is able to suppress the activity of miRNA, release the inhibition of miRNA on its targeted mRNA, and promote the expression of target gene, through binding to miRNA . CircRNA plays an important role in cell proliferation, differentiation, apoptosis, migration and body metabolism .
CircRNA participates in the pathogenesis of many proliferative diseases, such as hyper-proliferation of keratinocytes in psoriasis , hypoxia-induced proliferation of vascular endothelial cells , skeletal myoblast proliferation , and so on. Currently, hyper-proliferation of keratinocytes has been discovered in middle ear cholesteatoma, demonstrated by significantly high expression of numerous epithelial cell proliferation markers such as proliferating cell nuclear antigen (PCNA), Ki-67, cytokeratin (CK) 13, and CK16 in middle ear cholesteatoma epithelium [13,14,15]. However, the role of circRNA in the etiopathogenesis of middle ear cholesteatoma is rarely investigated so far.
In our research, we investigated the differential expression profiling of circRNAs between acquired middle ear cholesteatoma and normal skin by microarray technology. Expression of several selected circRNAs was validated through quantitative real-time polymerase chain reaction (qRT-PCR), purposing to recognize meaningful circRNAs contributing to the etiopathogenesis of middle ear cholesteatoma.
Differential expressions of circRNAs in cholesteatoma and normal skin
In our study, a total of 13,562 expressed circRNAs were detected through microarray analysis. Compared with normal skin, totally 93 up-regulated (hsa_circRNA_104327, hsa_circRNA_404655, etc.) and 85 down-regulated (hsa_circRNA_000319, hsa_circRNA_048764, etc.) circRNAs were identified in middle ear cholesteatoma tissues with a value of fold change (FC) ≥ 2 (P < 0.05). The differentially expressed circRNAs were displayed in a volcano plot (Fig. 1A), and then subjected to a scatter plot (Fig. 1B) and hierarchical clustering analysis (Fig. 1C), to display the variation and differentiation between middle ear cholesteatoma and normal skin. Moreover, among these 178 significantly differentially expressed circRNAs, a total of 145 were exons. Discoveries in our study indicated that these differentially expressed circRNAs might contribute to the etiopathogenesis of middle ear cholesteatoma.
Validation of differentially expressed circRNAs
To validate the circRNAs microarray result, a total of 8 differentially expressed circRNAs were selected for qRT-PCR analysis. In detail, 5 up-regulated circRNAs (hsa_circRNA_103670, hsa_circRNA_048764, hsa_circRNA_404864, hsa_circRNA_104327, hsa_circRNA_404655), and 3 down-regulated circRNAs (hsa_circRNA_101965, hsa_circRNA_000319, hsa_circRNA_100927) were validated. The validation result was in accordance with the microarray discoveries partly: compared with normal skin, expressions of hsa_circRNA_104327 and hsa_circRNA_404655 were significantly higher in middle ear cholesteatoma (P < 0.05), while the expression of hsa_circRNA_000319 significantly decreased (P < 0.05) (Fig. 2).
Gene ontology (GO) and Kyoto Encyclopedia Genes and Genomes (KEGG)
To preliminarily investigate the function of differentially expressed circRNAs in middle ear cholesteatoma, GO classification and KEGG pathway analysis were carried out on presumptive target genes of validated circRNAs (hsa_circRNA_104327, hsa_circRNA_404655, and hsa_circRNA_000319). Three domains including biological process (BP), cellular component (CC) and molecular function (MF) in GO classification were investigated. In this research, a total of 205 BP, 28 CC, and 28 MF domains enriched significantly (P < 0.05). The top-10 enriched BP, CC, and MF domains in GO classification are displayed in Fig. 3. Meanwhile, as for KEGG pathway analysis, a total of 11 signaling pathways enriched significantly (Fig. 4).
ceRNA network analysis
Based on the ceRNA hypothesis, circRNAs can bind to miRNAs and suppress the activity of miRNAs through MRE. In our study, 3 validated differentially expressed circRNAs with statistical significance including hsa_circRNA_000319, hsa_circRNA_104327, and hsa_circRNA_404655 were selected to perform a circRNA–miRNA–mRNA ceRNA network analysis.
The ceRNA network was composed of 3 circRNAs, 309 miRNAs, and 131 mRNAs (Additional file 1: Fig. S1). We discovered that hsa_circRNA_000319 interacted with hsa-miRNA-4436b-5p, hsa_circRNA_104327 interacted with miRNA-152-5p, and hsa_circRNA_404655 interacted with miRNA-3664-3p, with a complementarity based on 7mer-m8, 6mer, offset 6mer, or 8mer sites (Fig. 5). Moreover, in our previous study , compared with normal skin, the expression of miRNA-4436b-5p was demonstrated to be significantly up-regulated in middle ear cholesteatoma (FC > 2, P < 0.05), while miRNA-152-5p and miRNA-3664-3p were obviously down-regulated (FC < 0.5, P < 0.05), which were synchronous with expressions of hsa_circRNA_000319, hsa_circRNA_104327, and hsa_circRNA_404655 in middle ear cholesteatoma. Additionally, several targeted mRNAs for miRNA-3664-3p were discovered, including IQ motif-containing GTPase activating protein 3 (IQGAP3), receptor tyrosine kinase-like orphan receptor 2 (ROR2), ephrin type-B receptor 1 (EPHB1), nucleophosmin 1 (NPM1), chromosome condensation 2 (RCC2), nuclear receptor coactivator 5 (NCOA5), and so on.
Prediction of N6-methyladenosine (m6A) on differentially expressed circRNAs
m6A is a post-transcriptional methylation modification that widely presents at the adenosine bases of RNA transcripts. This modification has been suggested to regulate the function of circRNAs. In our study, we discovered that hsa_circRNA_000319 possessed 4 m6A sites with very high confidence in total, and hsa_circRNA_404655 had 3 m6A sites with high confidence, while no m6A site with very high or high confidence was found for hsa_circRNA_104327. Details are displayed in Table 1 and Fig. 6.
CircRNA is a special endogenous non-coding RNA with a unique ring structure, which regulates gene expression at the post-transcriptional level and plays an important role in a series of life processes, such as cell proliferation, differentiation, apoptosis, migration and metabolism [5, 9]. However, there are few studies on circRNAs in middle ear cholesteatoma, and the role of circRNAs in the etiopathogenesis of middle ear cholesteatoma is rarely investigated so far. We totally identified 93 significantly up-regulated and 85 down-regulated circRNAs in cholesteatoma compared with normal skin. Moreover, qRT-PCR was performed to verify the microarray discoveries by detecting the expression of 8 selected circRNAs. Expressions of hsa_circRNA_104327, hsa_circRNA_404655 and hsa_circRNA_000319 were in line with the microarray results, suggesting that these three circRNAs might play a pivotal role in the etiopathogenesis of middle ear cholesteatoma. The difference between microarray results and qRT-PCR validation might be caused by the false-positive results from microarray analysis. Furthermore, a validation with a large sample size will be conducted in our further research.
CircRNA is closely related to the pathogenesis of a variety of proliferative diseases, such as psoriasis , hypoxia-induced vascular endothelial cell proliferation  and so on. In the process of these proliferative diseases, circRNA regulates the expression of related target genes as miRNA "sponge" and participates in the etiopathogenesis of these diseases. Wang et al.  found that circTitin activated insulin-like growth factor 2 (IGF2)/phosphatidylinositol 3-kinase (PI3K)/Akt signal pathway through adsorption of miRNA-432, which led to the proliferation and differentiation of skeletal myoblasts. In addition, circRNA is also closely related to tumorigenesis. Studies discovered an abnormal expression of circRNAs in many kinds of tumor tissues, including breast cancer , liver cancer  and esophageal cancer . Additionally, circRNAs were demonstrated to participate in tumor cell hyper-proliferation, invasion and metastasis. For example, the expression of circSMARCA5 was significantly down-regulated in hepatocellular carcinoma, and was closely related to the hyper-proliferation, invasion and metastasis of hepatocellular carcinoma cells. Transfection of circSMARCA5 into the hepatocellular carcinoma cells up-regulated the expression of tissue inhibitor of matrix metalloproteinase 3 (TIMP3) through adsorption of miRNA-17-3p and miRNA-181b-5p, and ultimately inhibited the proliferation, invasion and metastasis of hepatocellular carcinoma cells . A study found that hsa_circRNA_102747 was significantly down-regulated in cholesteatoma and acted as the endogenous sponge of miRNA-21-3p , a miRNA promoting excessive proliferation of cholesteatoma keratinocytes . Meanwhile, hsa_circRNA_101458 was also obviously down-regulated in cholesteatoma and was able to interact with miRNA-let-7a-3p , a miRNA which was confirmed to suppress proliferation of keratinocytes in cholesteatoma . Consequently, circRNAs might play an important role in the etiopathogenesis of middle ear cholesteatoma.
So far, the exact function of circRNAs in middle ear cholesteatoma still remains unclear. According to this research and our previous research , expression of hsa_circRNA_404655 significantly increased in middle ear cholesteatoma, while miRNA-3664-3p obviously decreased and interacted with hsa_circRNA_404655. The decreased expression of miRNA-3664-3p might lead to an increased expression of growth differentiation factor 15 (GDF15), which consequently increased the expression of p53 and promoted tumor progression in spinal cord glioblastoma . Although cholesteatoma is pathologically benign, it is locally invasive and can lead to bone resorption. We speculate that hsa_circRNA_404655, as a sponge of miRNA-3664-3p, may participate in the proliferation and local invasion of cholesteatoma. Moreover, we also discovered several targeted mRNAs for hsa_circRNA_404655 and miRNA-3664-3p, including IQGAP3, ROR2, EPHB1, NPM1, RCC2, NCOA5, and so on. Although these targeted mRNAs have not been investigated in middle ear cholesteatoma, their functions were explored in other proliferative and invasive diseases. For instance, IQGAP3, EPHB1, ROR2 and RCC2 were demonstrated to obviously promote cell migration, invasion, and epithelial-to-mesenchymal transition in hepatocellular carcinoma, medulloblastoma and breast carcinoma, respectively [24,25,26,27]. Moreover, NPM1 contributed to cell proliferation through phosphatase and tensin homologue deleted on chromosome ten (PTEN) inactivation , while NCOA5 promoted proliferation, migration and invasion of colorectal cancer cells via activation of PI3K/Akt pathway . Similarly, epithelial-to-mesenchymal transition, PTEN inactivation and PI3K/Akt pathway activation were also noticed in middle ear cholesteatoma [30, 31]. These results provide new insights into the pathogenesis of cholesteatoma and suggest that circRNAs may be a potential target for the treatment of cholesteatoma. The exact mechanism needs to be further investigated.
Through KEGG pathway analysis, autophagy and hypoxia inducible factor-1 (HIF-1) signaling pathway were identified. Li et al.  proposed that enhanced autophagy might contribute to the pathogenesis of middle ear cholesteatoma through PI3K/Akt pathway activation. On the contrary, Ho et al.  put forward that autophagy was significantly suppressed in cholesteatoma. Moreover, activated HIF-1 pathway was detected in middle ear cholesteatoma, indicating that cholesteatoma might be hypoxic . Furthermore, KEGG pathway analysis exhibited several potential signaling pathways, which broadened our horizons on the etiopathogenesis of middle ear cholesteatoma.
More and more evidence shows that m6A modification can regulate the function of circRNAs. For instance, m6A modification was able to initiate circRNAs translation through the cooperation of eukaryotic translation initiation factor 4 gamma 2 (eIF4G2) and m6A reader YTH domain family 3 (YTHDF3) . Park et al. . demonstrated that m6A modification could reduce the stability of circRNAs. In addition, m6A modification can also regulate cytoplasmic output. In previous research, YTH domain-containing protein 1 (YTHDC1) mediated the nuclear output of circNSUN2 in a m6A methylation-dependent manner, thereby increasing cytoplasmic expression of circNSUN2 and promoting metastasis of colorectal cancer . Moreover, it has been reported that foreign circRNAs were able to induce T cell activation, antibody production, and anti-tumor immunity, and m6A modification on circRNAs inhibited immune gene activation . In our study, we predicted that hsa_circRNA_000319 possessed 4 m6A sites with very high confidence in total, and hsa_circRNA_404655 had 3 m6A sites with high confidence. Experiments to verify m6A modification on differentially expressed circRNAs will be performed in our further study.
In summary, 93 up-regulated and 85 down-regulated circRNAs were identified in middle ear cholesteatoma tissues in the current study. Compared with normal skin, the qRT-PCR validation for differentially expressed circRNAs demonstrated that hsa_circRNA_104327 and hsa_circRNA_404655 were expressed significantly higher while hsa_circRNA_000319 expression significantly decreased in cholesteatoma tissue, suggesting that circRNAs may contribute to the pathogenesis of middle ear cholesteatoma. Further researches should be conducted to investigate the exact mechanism of these differentially expressed circRNAs in the pathogenesis of middle ear cholesteatoma. Targeting on theses circRNAs may provide a new strategy for middle ear cholesteatoma therapy in the future.
Materials and methods
Patients and samples
A total of 16 acquired middle ear cholesteatoma tissues were obtained from patients undergoing cholesteatoma surgery between August 2019 and December 2019. Meanwhile, 14 normal skin (from incision ethier endaural or postaturicular) samples were collected to be recognized as controls. All samples were stored in liquid nitrogen at once for microarray analysis and qRT-PCR validation.
Total RNAs were extracted from 6 middle ear cholesteatoma and 4 normal skin tissues using TRIzol reagent (Invitrogen; USA) according to the manufacturer's instructions. Then, RNasey mini kit (QIAGEN, German) was used to purify RNAs. Meanwhile, RNA concentrations were determined through NanoDrop ND-1000 instrument (Thermo, USA), and RNA integrity was assessed through gel electrophoresis.
CircRNA labeling and microarray hybridization
The sample preparation and circRNA microarray hybridization were performed based on the Arraystar’s standard protocols. Briefly, total RNAs were digested with Rnase R (Epicentre, Inc.) to remove linear RNAs and enrich circRNAs. Then, the enriched circRNAs were amplified and transcribed into fluorescent cRNA utilizing a random priming method (Arraystar Super RNA Labeling Kit; Arraystar). The labeled cRNAs were hybridized onto the Arraystar Human circRNA Array V2 (8 × 15K, Arraystar), and incubated for 17 h at 65 °C in an Agilent Hybridization Oven. After having washed the slides, the arrays were scanned by the Agilent Scanner G2505C.
Image acquisition and microarray analysis of circRNAs
Agilent Feature Extraction software (version 220.127.116.11) was used to analyze acquired array images. Quantile normalization and subsequent data processing were performed using the R software limma package. Differentially expressed circRNAs with statistical significance between middle ear cholesteatoma and normal skin were identified using FC cutoff or volcano plot filtering, respectively. The value of FC ≥ 2 and P < 0.05 was regarded as significantly differentially expressed. Hierarchical Clustering was performed to show the distinguishable circRNAs expression pattern between middle ear cholesteatoma and normal skin samples. Finally, the circRNA–miRNA interaction was predicted with Arraystar's home-made miRNA target prediction software based on TargetScan (http://www.targetscan.org) and miRanda (http://www.microrna.org).
Validation for differentially expressed circRNAs by qRT-PCR
Total cellular RNA was, respectively, purified from 10 cholesteatoma and 10 normal skin tissues using Trizol reagent (Invitrogen, USA), precipitated with alcohol, and then stored at − 70 °C until use. Extracted RNA was quantified by measuring optical density at 260 nm, and then 0.4 µg RNA was reversely transcribed to cDNA in a 20-µl reaction system according to the manufacturer’s protocol. Then, quantitative PCR reactions were conducted using QuantStudioTM 5 System (Applied Biosystems, USA) with 2 µl cDNA sample adding into a 8-µl reaction system. Amplification conditions were controlled as follows: 95 °C for 10 min for 1 cycle, followed by 40 cycles of 95 °C for 10 s and 60 °C for 60 s. The specific amplification primer pairs are shown in Table 2.
ceRNA network analysis
Significantly differentially expressed circRNAs including hsa_circRNA_000319, hsa_circRNA_104327, and hsa_circRNA_404655 were selected to perform ceRNA network analysis. The cytoscape software V3.5.0 (San Diego, CA, USA) was adopted to construct circRNA-associated ceRNA network.
GO and KEGG pathway analysis
GO classification (http://www.geneontology.org) and KEGG pathway analysis (http://www.genome.jp/kegg/) were applied to explore functions of target genes in the pathogenesis of middle ear cholesteatoma. The GO classification covers three domains including BP, CC, and MF, and describes the gene functions. Pathway analysis was conducted to predict molecular interactions and reaction networks through mapping genes to KEGG pathways. A P < 0.05 was considered to be statistically significant.
Prediction of m6A modification on differentially expressed circRNAs
Sequence-based RNA adenosine methylation site predictor (SRAMP) website (http://www.cuilab.cn/sramp) was searched in order to explore the m6A modification on differentially expressed circRNAs. Significantly differentially expressed circRNAs including hsa_circRNA_000319, hsa_circRNA_104327, and hsa_circRNA_404655 were selected to conduct the prediction of m6A modification sites on the RNA sequences of these 3 circRNAs.
Availability of data and materials
All data generated and analyzed during this study are included in this published article, detailed data have been uploaded in https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE171006.
MiRNA response element
Competing endogenous RNA
Proliferating cell nuclear antigen
Quantitative real-time polymerase chain reaction
Kyoto Encyclopedia Genes and Genomes
IQ motif-containing GTPase activating protein 3
Receptor tyrosine kinase-like orphan receptor 2
Ephrin type-B receptor 1
Chromosome condensation 2
Nuclear receptor coactivator 5
Insulin-like growth factor 2
Tissue inhibitor of matrix metalloproteinase 3
Growth differentiation factor 15
Phosphatase and tensin homologue deleted on chromosome ten
Hypoxia inducible factor-1
Eukaryotic translation initiation factor 4 gamma 2
YTH domain family 3
YTH domain-containing protein 1
Retinoic acid-inducible gene I
Castle JT. Cholesteatoma pearls: practical points and update. Head Neck Pathol. 2018;12(3):419–29.
Xie S, Wang X, Ren J, Liu W. The role of bone resorption in the etiopathogenesis of acquired middle ear cholesteatoma. Eur Arch Otorhinolaryngol. 2017;274(5):2071–8.
Bhutta MF, Williamson IG, Sudhoff HH. Cholesteatoma. BMJ. 2011;342:d1088.
Zhang ZG, Chen SJ, Sun W, Zheng YQ. Classification and surgical management of localized attic cholesteatoma: single-institution experience and follow-up. ORL. 2010;72(2):96–100.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32(5):453–61.
Galvão-Lima LJ, Morais AHF, Valentim RAM, Barreto EJSS. miRNAs as biomarkers for early cancer detection and their application in the development of new diagnostic tools. Biomed Eng Online. 2021;20(1):21.
Dai R, Ahmed SA. MicroRNA, a new paradigm for understanding immunoregulation, inflammation, and autoimmune diseases. Transl Res. 2011;157(4):163–79.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Hansen TB, Kjems J, Damgaard CK. Circular RNA and miR-7 in cancer. Cancer Res. 2013;73(18):5609–12.
Moldovan LI, Hansen TB, Venø MT, Okholm TLH, Andersen TL, Hager H, et al. High-throughput RNA sequencing from paired lesional- and non-lesional skin reveals major alterations in the psoriasis circRNAome. BMC Med Genomics. 2019;12(1):174.
Dang RY, Liu FL, Li Y. Circular RNA hsa_circ_0010729 regulates vascular endothelial cell proliferation and apoptosis by targeting the miR-186/HIF-1α axis. Biochem Biophys Res Commun. 2017;490(2):104–10.
Wang X, Cao X, Dong D, Shen X, Cheng J, Jiang R, et al. Circular RNA TTN Acts As a miR-432 sponge to facilitate proliferation and differentiation of myoblasts via the IGF2/PI3K/AKT signaling pathway. Mol Ther Nucleic Acids. 2019;18:966–80.
Hamed MA, Nakata S, Shiogama K, Suzuki K, Sayed RH, Nishimura Y, et al. Cytokeratin 13, cytokeratin 17, and Ki-67 expression in human acquired cholesteatoma and their correlation with its destructive capacity. Clin Exp Otorhinolaryngol. 2017;10(3):213–20.
Liu W, Ren H, Ren J, Yin T, Hu B, Xie S, et al. The role of EGFR/PI3K/Akt/cyclinD1 signaling pathway in acquired middle ear cholesteatoma. Mediators Inflamm. 2013;2013:651207.
Liu W, Yin T, Ren J, Li L, Xiao Z, Chen X, et al. Activation of the EGFR/Akt/NF-κB/cyclinD1 survival signaling pathway in human cholesteatoma epithelium. Eur Arch Otorhinolaryngol. 2014;271(2):265–73.
Xie S, Liu X, Pan Z, Chen X, Peng A, Yin T, et al. Microarray analysis of differentially-expressed MicroRNAs in acquired middle ear cholesteatoma. Int J Med Sci. 2018;15(13):1547–54.
Pan G, Mao A, Liu J, Lu J, Ding J, Liu W. Circular RNA hsa_circ_0061825 (circ-TFF1) contributes to breast cancer progression through targeting miR-326/TFF1 signalling. Cell Prolif. 2020;53(2):e12720.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, et al. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68(6):1214–27.
Li F, Zhang L, Li W, Deng J, Zheng J, An M, et al. Circular RNA ITCH has inhibitory effect on ESCC by suppressing the Wnt/β-catenin pathway. Oncotarget. 2015;6(8):6001–13.
Gao J, Tang Q, Xue R, Zhu X, Wang S, Zhang Y, et al. Comprehensive circular RNA expression profiling with associated ceRNA network reveals their therapeutic potential in cholesteatoma. Oncol Rep. 2020;43(4):1234–44.
Chen X, Li X, Qin Z. MicroRNA-21 promotes the proliferation and invasion of cholesteatoma keratinocytes. Acta Otolaryngol. 2016;136(12):1261–6.
Zhang W, Chen X, Qin Z. MicroRNA let-7a suppresses the growth and invasion of cholesteatoma keratinocytes. Mol Med Rep. 2015;11(3):2097–103.
An T, Fan T, Zhang XQ, Liu YF, Huang J, Liang C, et al. Comparison of Alterations in miRNA Expression in Matched Tissue and Blood Samples during Spinal Cord Glioma Progression. Sci Rep. 2019;9(1):9169.
Shi Y, Qin N, Zhou Q, Chen Y, Huang S, Chen B, et al. Role of IQGAP3 in metastasis and epithelial-mesenchymal transition in human hepatocellular carcinoma. J Transl Med. 2017;15(1):176.
Xu J, Shi J, Tang W, Jiang P, Guo M, Zhang B, et al. ROR2 promotes the epithelial-mesenchymal transition by regulating MAPK/p38 signaling pathway in breast cancer. J Cell Biochem. 2020;121(10):4142–53.
Chen Z, Wu W, Huang Y, Xie L, Li Y, Chen H, et al. RCC2 promotes breast cancer progression through regulation of Wnt signaling and inducing EMT. J Cancer. 2019;10(27):6837–47.
Bhatia S, Baig NA, Timofeeva O, Pasquale EB, Hirsch K, MacDonald TJ, et al. Knockdown of EphB1 receptor decreases medulloblastoma cell growth and migration and increases cellular radiosensitization. Oncotarget. 2015;6(11):8929–46.
Boudra R, Lagrafeuille R, Lours-Calet C, de Joussineau C, Loubeau-Legros G, Chaveroux C, et al. mTOR transcriptionally and post-transcriptionally regulates Npm1 gene expression to contribute to enhanced proliferation in cells with Pten inactivation. Cell Cycle. 2016;15(10):1352–62.
Sun K, Wang S, He J, Xie Y, He Y, Wang Z, et al. NCOA5 promotes proliferation, migration and invasion of colorectal cancer cells via activation of PI3K/AKT pathway. Oncotarget. 2017;8(64):107932–46.
Takahashi M, Yamamoto-Fukuda T, Akiyama N, Motegi M, Yamamoto K, Tanaka Y, et al. Partial epithelial-mesenchymal transition was observed under p63 expression in acquired middle ear cholesteatoma and congenital cholesteatoma. Otol Neurotol. 2019;40(8):e803–11.
Liu W, Yin T, Ren J, Wang Y, Shen X, Li L, et al. Expression of PTEN, P-Akt and NF-κB in middle ear cholesteatoma epithelium and its significance. Chinese J Otol. 2011;9(2):206–11.
Li Q, Ao Y, Yu Q, Zhou S. Role of autophagy in acquired cholesteatoma. Otol Neurotol. 2019;40(10):e993–8.
Ho KY, Huang CJ, Hung CC, Wu YR, Chiu CC, Chien CY, et al. Autophagy is deficient and may be negatively regulated by SERPINB3 in middle ear cholesteatoma. Otol Neurotol. 2020;41(7):e881–8.
Adunka O, Gstoettner W, Knecht R, Kierner AC. Expression of hypoxia inducible factor 1 alpha and Von Hippel Lindau protein in human middle ear cholesteatoma. Laryngoscope. 2003;113(7):1210–5.
Yang Y, Fan X, Mao M, Song X, Wu P, Zhang Y, et al. Extensive translation of circular RNAs driven by N6-methyladenosine. Cell Res. 2017;27(5):626–41.
Park OH, Ha H, Lee Y, Boo SH, Kwon DH, Song HK, et al. Endoribonucleolytic cleavage of m6A-containing RNAs by RNase P/MRP complex. Mol Cell. 2019;74(3):494-507.e8.
Chen RX, Chen X, Xia LP, Zhang JX, Pan ZZ, Ma XD, et al. N6-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019;10(1):4695.
Chen YG, Chen R, Ahmad S, Verma R, Kasturi SP, Amaya L, et al. N6-methyladenosine modification controls circular RNA immunity. Mol Cell. 2019;76(1):96-109.e9.
This work was supported by the National Natural Science Foundation of China [Grant Number 82071036, Wei Liu, Grant Number 82000973, Shumin Xie], Natural Science Foundation of Hunan Province [Grant Number 2019JJ50967, Shumin Xie, 2017JJ3441, Wei Liu].
Ethics approval and consent to participate
This study procedure was approved by the Ethics Committee of The Second Xiangya Hospital of Central South University, and written informed consent was acquired from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xie, S., Jin, L., Yin, T. et al. Microarray analysis and functional prediction of differentially expressed circular RNAs in acquired middle ear cholesteatoma. BioMed Eng OnLine 20, 129 (2021). https://doi.org/10.1186/s12938-021-00960-x
- Microarray analysis