Profiling and Bioinformatics Analysis of Differentially Expressed circRNAs in Spinal Ligament Tissues of Patients with Ankylosing Spondylitis

Recent studies have reported that circular RNAs (circRNAs) play a crucial regulatory role in a variety of human diseases. However, the roles of circRNAs in ankylosing spondylitis (AS) remain unclear. In this study, we conducted circRNA expression profiling of the spinal ligament tissues of patients with AS by RNA sequencing (RNA-seq) and analyzed the potential functions of differentially expressed circRNA by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to investigate the potential mechanisms associated with AS. The results showed that a total of 1,172 circRNAs were detected in the spinal ligament tissue samples, of which 123 circRNAs were significantly differentially expressed by a fold change≥1.5 and p value < 0.05. Among these, 57 circRNAs were upregulated, and 66 were downregulated. GO and KEGG analyses demonstrated that the differentially expressed circRNAs were mainly involved in the regulation of biological processes of peptidyl-serine phosphorylation and human immune system that may be related to AS. In addition, the circRNA/miRNA interaction networks were established to predict the potential roles of differentially expressed circRNAs by bioinformatics analysis. Taken together, these results revealed the expression profiles of circRNAs and the potential functions of the differentially expressed circRNAs in the spinal ligament tissue of patients with AS, which may provide new clues for understanding the mechanisms associated with AS, and proceed to identify novel potential molecular targets for the diagnoses and treatment of AS.

Download Full-text

Differential Expression Profiles of Long Noncoding RNA and mRNA of Osteogenically Differentiated Mesenchymal Stem Cells in Ankylosing Spondylitis

The Journal of Rheumatology ◽

10.3899/jrheum.151181 ◽

2016 ◽

Vol 43 (8) ◽

pp. 1523-1531 ◽

Cited By ~ 24

Author(s):

Zhongyu Xie ◽

Jinteng Li ◽

Peng Wang ◽

Yuxi Li ◽

Xiaohua Wu ◽

...

Keyword(s):

Stem Cells ◽

Mesenchymal Stem Cells ◽

Ankylosing Spondylitis ◽

Osteogenic Differentiation ◽

Noncoding Rna ◽

Bioinformatics Analysis ◽

Expression Profiles ◽

Differentially Expressed ◽

Qrt Pcr ◽

Pcr Assays

Objective.We previously demonstrated that mesenchymal stem cells (MSC) from patients with ankylosing spondylitis (AS; ASMSC) have a greater osteogenic differentiation capacity than MSC from healthy donors (HDMSC) and that this difference underlies the pathogenesis of pathological osteogenesis in AS. Here we compared expression levels of long noncoding RNA (lncRNA) and mRNA between osteogenically differentiated ASMSC and HDMSC and explored the precise mechanism underlying abnormal osteogenic differentiation in ASMSC.Methods.HDMSC and ASMSC were induced with osteogenic differentiation medium for 10 days. Microarray analyses were then performed to identify lncRNA and mRNA differentially expressed between HDMSC and ASMSC, which were then subjected to bioinformatics analysis and confirmed by quantitative real-time PCR (qRT-PCR) assays. In addition, coding-non-coding gene co-expression (CNC) networks were constructed to examine the relationships between the lncRNA and mRNA expression patterns.Results.A total of 520 lncRNA and 665 mRNA were differentially expressed in osteogenically differentiated ASMSC compared with HDMSC. Bioinformatics analysis revealed 64 signaling pathways with significant differences, including transforming growth factor-β signaling. qRT-PCR assays confirmed the reliability of the microarray data. The CNC network indicated that 4 differentially expressed lncRNA, including lnc-ZNF354A-1, lnc-LIN54-1, lnc-FRG2C-3, and lnc-USP50-2 may be involved in the abnormal osteogenic differentiation of ASMSC.Conclusion.Our study characterized the differential lncRNA and mRNA expression profiles of osteogenically differentiated ASMSC and identified 4 lncRNA that may participate in the abnormal osteogenic differentiation of ASMSC. These results provide insight into the pathogenesis of pathological osteogenesis in AS.

Download Full-text

Construction and investigation of a circRNA-associated ceRNA regulatory network in Tetralogy of Fallot

BMC Cardiovascular Disorders ◽

10.1186/s12872-021-02217-w ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Haifei Yu ◽

Xinrui Wang ◽

Hua Cao

Keyword(s):

Tetralogy Of Fallot ◽

Heart Development ◽

Expression Profiles ◽

Screening Strategy ◽

Differentially Expressed ◽

Circular Rnas ◽

Tissue Samples ◽

Cerna Network ◽

Positive Correlations

Abstract Background As the most frequent type of cyanotic congenital heart disease (CHD), tetralogy of Fallot (TOF) has a relatively poor prognosis without corrective surgery. Circular RNAs (circRNAs) represent a novel class of endogenous noncoding RNAs that regulate target gene expression posttranscriptionally in heart development. Here, we investigated the potential role of the ceRNA network in the pathogenesis of TOF. Methods To identify circRNA expression profiles in TOF, microarrays were used to screen the differentially expressed circRNAs between 3 TOF and 3 control human myocardial tissue samples. Then, a dysregulated circRNA-associated ceRNA network was constructed using the established multistep screening strategy. Results In summary, a total of 276 differentially expressed circRNAs were identified, including 214 upregulated and 62 downregulated circRNAs in TOF samples. By constructing the circRNA-associated ceRNA network based on bioinformatics data, a total of 19 circRNAs, 9 miRNAs, and 34 mRNAs were further screened. Moreover, by enlarging the sample size, the qPCR results validated the positive correlations between hsa_circ_0007798 and HIF1A. Conclusions The findings in this study provide a comprehensive understanding of the ceRNA network involved in TOF biology, such as the hsa_circ_0007798/miR-199b-5p/HIF1A signalling axis, and may offer candidate diagnostic biomarkers or potential therapeutic targets for TOF. In addition, we propose that the ceRNA network regulates TOF progression.

Download Full-text

The circular RNA circSlc7a11 promotes bone cancer pain pathogenesis in rats by modulating LLC-WRC 256 cell proliferation and apoptosis

Molecular and Cellular Biochemistry ◽

10.1007/s11010-020-04020-1 ◽

2021 ◽

Author(s):

Han-Wen Chen ◽

Xiao-Xia Zhang ◽

Zhu-Ding Peng ◽

Zu-Min Xing ◽

Yi-Wen Zhang ◽

...

Keyword(s):

Cell Proliferation ◽

Cancer Pain ◽

Expression Profiles ◽

Bone Cancer ◽

Circular Rna ◽

Differentially Expressed ◽

Bone Cancer Pain ◽

Circular Rnas ◽

Altered Expression ◽

Proliferation And Apoptosis

AbstractTreatment of bone cancer pain (BCP) caused by bone metastasis in advanced cancers remains a challenge in clinical oncology, and the underlying mechanisms of BCP are poorly understood. This study aimed to investigate the pathogenic roles of circular RNAs (circRNAs) in regulating cancer cell proliferation and BCP development. Eight differentially expressed circRNAs in the rat spinal cord were validated by agarose gel electrophoresis and Sanger sequencing. Expression of circRNAs and mRNAs was detected by quantitative RT-PCR. MTS assay and flow cytometry were performed to analyze cell proliferation and apoptosis, respectively. Differentially expressed mRNA profiles were characterized by deep RNA sequencing, hierarchical clustering, and functional categorization. The interactions among circRNAs, microRNAs (miRNAs), and mRNAs were predicted using TargetScan. Additionally, western blot was performed to determine the protein levels of Pax8, Isg15, and Cxcl10. Multiple circRNAs were differentially expressed in the spinal cords of BCP model rats; of these, circSlc7a11 showed the greatest increase in expression. The overexpression of circSlc7a11 significantly promoted cell proliferation and repressed apoptosis of LLC-WRC 256 and UMR-106 cells, whereas circSlc7a11 silencing produced the opposite effects. Altered expression of circSlc7a11 also induced substantial changes in the mRNA expression profiles of LLC-WRC 256 cells; these changes were linked to multiple apoptotic processes and signaling pathways, such as the chemokine signaling pathway, and formed a complex circRNA/miRNA/mRNA network. Additionally, Pax8, Isg15, and Cxc110 protein level in LLC-WRC 256 cells was consistent with the mRNA results. The circRNA circSlc7a11 regulates rat BCP development by modulating LLC-WRC 256 cell proliferation and apoptosis through multiple-signaling mechanisms.

Download Full-text

Circular RNA hsa_circ_0000277 sequesters miR-4766-5p to upregulate LAMA1 and promote esophageal carcinoma progression

Cell Death and Disease ◽

10.1038/s41419-021-03911-5 ◽

2021 ◽

Vol 12 (7) ◽

Author(s):

Peng Li Zhou ◽

Zhengyang Wu ◽

Wenguang Zhang ◽

Miao Xu ◽

Jianzhuang Ren ◽

...

Keyword(s):

Molecular Mechanisms ◽

Bioinformatics Analysis ◽

Epithelial Mesenchymal Transition ◽

Expression Profiles ◽

Luciferase Reporter ◽

Migration And Invasion ◽

Circular Rnas ◽

Mesenchymal Transition ◽

Dismal Prognosis ◽

Advanced Tumor

AbstractGrowing evidence has indicated that circular RNAs (circRNAs) play a pivotal role as functional RNAs in diverse cancers. However, most circRNAs involved in esophageal squamous cell carcinoma (ESCC) remain undefined, and the underlying molecular mechanisms mediated by circRNAs are largely unclear. Here, we screened human circRNA expression profiles in ESCC tissues and found significantly increased expression of hsa_circ_0000277 (termed circPDE3B) in ESCC tissues and cell lines compared to the normal controls. Moreover, higher circPDE3B expression in patients with ESCC was correlated with advanced tumor-node-metastasis (TNM) stage and dismal prognosis. Functional experiments demonstrated that circPDE3B promoted the tumorigenesis and metastasis of ESCC cells in vitro and in vivo. Mechanistically, bioinformatics analysis, a dual-luciferase reporter assay, and anti-AGO2 RNA immunoprecipitation showed that circPDE3B could act as a competing endogenous RNA (ceRNA) by harboring miR-4766-5p to eliminate the inhibitory effect on the target gene laminin α1 (LAMA1). In addition, LAMA1 was significantly upregulated in ESCC tissues and was positively associated with the aggressive oncogenic phenotype. More importantly, rescue experiments revealed that the oncogenic role of circPDE3B in ESCC is partly dependent on the miR-4766-5p/LAMA1 axis. Furthermore, bioinformatics analysis combined with validation experiments showed that epithelial-mesenchymal transition (EMT) activation was involved in the oncogenic functions of the circPDE3B–miR-4766-5p/LAMA1 axis in ESCC. Taken together, we demonstrate for the first time that the circPDE3B/miR-4766-5p/LAMA1 axis functions as an oncogenic factor in promoting ESCC cell proliferation, migration, and invasion by inducing EMT, implying its potential prognostic and therapeutic significance in ESCC.

Download Full-text

Identification of Gastric Cancer-Related Circular RNA through Microarray Analysis and Bioinformatics Analysis

BioMed Research International ◽

10.1155/2018/2381680 ◽

2018 ◽

Vol 2018 ◽

pp. 1-9 ◽

Cited By ~ 12

Author(s):

Wei Gu ◽

Ying Sun ◽

Xiong Zheng ◽

Jin Ma ◽

Xiao-Ying Hu ◽

...

Keyword(s):

Gene Expression ◽

Gastric Cancer ◽

Microarray Analysis ◽

Cancer Progression ◽

Pathway Analysis ◽

Bioinformatics Analysis ◽

Malignant Tumors ◽

Expression Profiles ◽

Circular Rnas ◽

Target Mrnas

Gastric cancer is one of the common malignant tumors worldwide. Increasing studies have indicated that circular RNAs (circRNAs) play critical roles in the cancer progression and have shown great potential as useful markers and therapeutic targets. However, the precise mechanism and functions of most circRNAs are still unknown in gastric cancer. In the present study, we performed a microarray analysis to detect circRNA expression changes between tumor samples and adjacent nontumor samples. The miRNA expression profiles were obtained from the National Center of Biotechnology Information Gene Expression Omnibus (GEO). The differentially expressed circRNAs and miRNAs were identified through fold change filtering. The interactions between circRNAs and miRNAs were predicted by Arraystar’s home-made miRNA target prediction software. After circRNA-related miRNAs and dysregulated miRNAs were intersected, 23 miRNAs were selected. The target mRNAs of miRNAs were predicted by TarBase v7.0. Gene ontology (GO) enrichment analysis and pathway analysis were performed using standard enrichment computational methods for the target mRNAs. The results of pathway analysis showed that p53 signaling pathway and hippo signal pathway were significantly enriched and CCND2 was a cross-talk gene associated with them. Finally, a circRNA-miRNA-mRNA regulation network was constructed based on the gene expression profiles and bioinformatics analysis results to identify hub genes and hsa_circRNA_101504 played a central role in the network.

Download Full-text

Identification of differentially expressed microRNAs in outgrowth embryos compared with blastocysts and non-outgrowth embryos in mice

Reproduction Fertility and Development ◽

10.1071/rd18161 ◽

2019 ◽

Vol 31 (4) ◽

pp. 645 ◽

Cited By ~ 6

Author(s):

Jihyun Kim ◽

Jaewang Lee ◽

Jin Hyun Jun

Keyword(s):

Bioinformatics Analysis ◽

Preimplantation Embryo ◽

Expression Profiles ◽

In Silico Analysis ◽

Differentially Expressed ◽

Recurrent Implantation Failure ◽

Preimplantation Embryo Development ◽

Mirna Expression Profiles ◽

Differentially Expressed Mirnas ◽

Silico Analysis

Recurrent implantation failure (RIF) is one of the main causes for the repeated failure of IVF, and the major reason for RIF is thought to be a miscommunication between the embryo and uterus. However, the exact mechanism underlying embryo–uterus cross-talk is not fully understood. The aim of the present study was to identify differentially expressed microRNAs (miRNAs) among blastocysts, non-outgrowth and outgrowth embryos in mice using microarray analysis. A bioinformatics analysis was performed to predict the potential mechanisms of implantation. The miRNA expression profiles differed significantly between non-outgrowth and outgrowth embryos. In all, 3163 miRNAs were detected in blastocysts and outgrowth embryos. Of these, 10 miRNA candidates (let-7b, miR-23a, miR-27a, miR-92a, miR-183, miR-200c, miR-291a, miR-425, miR-429 and miR-652) were identified as significant differentially expressed miRNAs of outgrowth embryos by in silico analysis. The expression of the miRNA candidates was markedly changed during preimplantation embryo development. In particular, let-7b-5p, miR-200c-3p and miR-23a-3p were significantly upregulated in outgrowth embryos compared with non-outgrowth blastocysts. Overall, differentially expressed miRNAs in outgrowth embryos compared with blastocysts and non-outgrowth embryos could be involved in embryo attachment, and interaction between the embryo proper and maternal endometrium during the implantation process.

Download Full-text

Global microarray profiling identifiedhsa_circ_0064428as a potential immune-associated prognosis biomarker for hepatocellular carcinoma

Journal of Medical Genetics ◽

10.1136/jmedgenet-2018-105440 ◽

2018 ◽

Vol 56 (1) ◽

pp. 32-38 ◽

Cited By ~ 16

Author(s):

Qiaoyou Weng ◽

Minjiang Chen ◽

Maoquan Li ◽

Yong-Fa Zheng ◽

Guoliang Shao ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Tumour Size ◽

Differentially Expressed ◽

Circular Rnas ◽

Tissue Samples ◽

Qrt Pcr ◽

Immunohistochemistry Staining ◽

Microarray Profiling ◽

Spss Software ◽

Infiltrating Lymphocytes

BackgroundIncreasing evidence has shown that circular RNAs (circRNAs) are involved tumourigenesis and metastasis of hepatocellular carcinoma (HCC); however, progression about its function in HCC is relatively slow. Here, we aimed to investigate whether plasma circRNAs could reflect the tumour-infiltrating lymphocytes (TILs) in HCC tumour tissues and serve as prognosis biomarker for HCC.MethodsTissue samples of patients with HCC were subjected to immunohistochemistry staining against CD8 to examine the TILs. Then, we investigated the expression profile of circRNAs by microarray between plasma of patients with HCC with high TILs and low TILs, and the differentially expressed circRNAs were validated with qRT-PCR. Statistical analysis was performed with SPSS software and GraphPad Prism.ResultsWe have demonstrated that patients with HCC with high TILs exhibit a significant better overall survival, suggesting clinical outcome could be predicted by TILs. Global circRNA microarray between plasma of patients with HCC with high TILs and low TILs successfully identified six differentially expressed novel circRNAs. Among them, the expression ofhsa_circ_0064428was significantly reduced in patients with HCC with high TILs but increased in patients with low TILs. Moreover,hsa_circ_0064428was negatively correlated with patient’s survival, tumour size and metastasis.ConclusionThese findings together imply thathsa_circ_0064428could be considered as a potential HCC prognosis biomarker. Future in-depth research is required to further illustrate the involvement ofhsa_circ_0064428in HCC tumourigenesis and metastasis.

Download Full-text

Identification of differentially expressed circular RNAs in HeLa cells infected with Chlamydia trachomatis

Pathogens and Disease ◽

10.1093/femspd/ftz062 ◽

2019 ◽

Vol 77 (7) ◽

Author(s):

Yuanjun Liu ◽

Chunmin Hu ◽

Yina Sun ◽

Haoqing Wu ◽

Xiaojun Chen ◽

...

Keyword(s):

Chlamydia Trachomatis ◽

Signaling Pathways ◽

Hela Cells ◽

Expression Profiles ◽

Quantitative Polymerase Chain Reaction ◽

Host Cells ◽

Differentially Expressed ◽

Chlamydia Infection ◽

Circular Rnas ◽

Mrna Microarray

ABSTRACT Non-coding circular RNAs (circRNAs) have been shown to have important roles in many diseases; however, no study has indicated circRNAs are involved in Chlamydia trachomatis infection. In this study, we used circRNA microarray to measure the global circRNA expression profiles in HeLa cells with or without C. trachomatis serovar E (Ct.E) infection. CircRNA/miRNA/mRNA interactions were predicted and bioinformatics analyses were performed. The differentially expressed circRNAs were selected according to our criterion for validation by reverse-transcription and quantitative polymerase chain reaction (RT-qPCR). The mRNA microarray was used to detect the mRNA expression profiles after Ct.E infection. Among 853 differentially expressed circRNAs, 453 were upregulated and 400 were downregulated after Ct.E infection. Target miRNAs and miRNA-targeted mRNAs of these circRNAs were predicted. RT-qPCR analysis indicated hsa_circRNA_001226, hsa_circRNA_007046 and hsa_circRNA_400027 were elevated similar to those determined in the circRNA microarray analysis. The mRNA microarray results showed 915 genes were upregulated and 619 genes were downregulated after Ct.E infection. Thirty-four differentially expressed genes overlapped in the bioinformatics and mRNA microarray results. KEGG pathway analysis revealed several signaling pathways, including endocytosis, MAPK and PI3P-Akt signaling pathways, that were targeted by circRNAs may play important roles in Chlamydia infection. This study provides evidence that circRNAs in host cells are involved in the process of Chlamydia infection.

Download Full-text

MicroRNA expression signatures of atrial fibrillation: The critical systematic review and bioinformatics analysis

Experimental Biology and Medicine ◽

10.1177/1535370219890303 ◽

2019 ◽

Vol 245 (1) ◽

pp. 42-53 ◽

Cited By ~ 2

Author(s):

Nan-Nan Shen ◽

Chi Zhang ◽

Zheng Li ◽

Ling-Cong Kong ◽

Xin-Hua Wang ◽

...

Keyword(s):

Atrial Fibrillation ◽

Systematic Review ◽

Sensitivity Analysis ◽

Mirna Expression ◽

Bioinformatics Analysis ◽

Detection Method ◽

Expression Profiles ◽

Differentially Expressed ◽

Cochrane Library ◽

Differentially Expressed Mirnas

Association between microRNA (miRNA) expression signatures and atrial fibrillation has been evaluated with inconsistent findings in different studies. This study aims to identify miRNAs that actually play vital role in pathophysiological process of atrial fibrillation and explore miRNA-targeted genes and the involved pathways. Relevant studies were retrieved from the electronic databases of Embase, Medline, and Cochrane Library to determine the miRNA expression profiles between atrial fibrillation subjects and non-atrial fibrillation controls. Robustness of results was assessed using sensitivity analysis. Subgroup analyses were performed based on species, miRNA detection method, sample source, and ethnicity. Quality assessment of studies was independently conducted according to QUADAS-2. Bioinformatics analysis was applied to explore the potential genes and pathways associated with atrial fibrillation, which were targeted by differentially expressed miRNAs. Form of pooled results was shown as log10 odds ratios (logORs) with 95% confidence intervals (CI), and random-effects model was used. In total, 40 articles involving 283 differentially expressed miRNAs were reported. And 51 significantly dysregulated miRNAs were identified in consistent direction, with 22 upregulated and 29 downregulated. Among above-mentioned miRNAs, miR-223-3p (logOR 6.473; P < 0.001) was the most upregulated, while miR-1-5p (logOR 7.290; P < 0.001) was the most downregulated. Subgroup analysis confirmed 53 significantly dysregulated miRNAs (21 upregulated and 32 downregulated) in cardiac tissue, with miRNA-1-5p and miRNA-223-3p being the most upregulated and downregulated miRNAs, respectively. Additionally, miR-328 and miR-1-5p were highly blood-specific, and miR-133 was animal-specific. In the detection method sub-groups, miRNA-29b and miRNA-223-3p were differentially expressed consistently. Four miRNAs, including miRNA-223-3p, miRNA-21, miRNA-328, and miRNA-1-5p, were consistently dysregulated in both Asian and non-Asian. Results of sensitivity analysis showed that 47 out of 51 (92.16%) miRNAs were dysregulated consistently. Totally, 51 consistently dysregulated miRNAs associated with atrial fibrillation were confirmed in this study. Five important miRNAs, including miR-29b, miR-328, miR-1-5p, miR-21, and miR-223-3p may act as potential biomarkers for atrial fibrillation. Impact statement Atrial fibrillation (AF) is considered as the most common arrhythmia, and it subsequently causes serious complications including thrombosis and heart failure that increase the social burden. The definite mechanisms underlying AF pathogenesis remain complicated and unclear. Many studies attempted to discover the transcriptomic changes using microarray technologies, and the present studies for this hot topic have assessed individual miRNAs profiles for AF. However, results of different articles are controversial and not each reported miRNA is actually associated with the pathogenesis of AF. The present systematic review and meta-analysis identified that 51 consistently dysregulated miRNAs were associated with AF. Of these miRNAs, five miRNAs (miRNA-1-5p, miRNA-328, miRNA-29b, miRNA-21, and miRNA-223-3p) may act as novel biomarkers for AF. The findings could offer a better description of the biological characteristics of miRNAs, meanwhile might serve as new target for the intervention and monitoring AF in future studies.

Download Full-text

Profiling and bioinformatics analysis of differentially expressed circular RNAs in human intervertebral disc degeneration

Acta Biochimica et Biophysica Sinica ◽

10.1093/abbs/gmz036 ◽

2019 ◽

Vol 51 (6) ◽

pp. 571-579 ◽

Cited By ~ 7

Author(s):

Shunmin Wang ◽

Jingchuan Sun ◽

Haisong Yang ◽

Weiguo Zou ◽

Bing Zheng ◽

...

Keyword(s):

Intervertebral Disc ◽

Disc Degeneration ◽

Intervertebral Disc Degeneration ◽

Expression Profiles ◽

Differentially Expressed ◽

Circular Rnas ◽

Functional Changes ◽

Qrt Pcr ◽

Reverse Transcription Pcr ◽

Microarray Expression

AbstractThe functional changes of nucleus pulposus (NP) cells are considered to be the initiating factors of intervertebral disc degeneration (IDD), and the differentially expressed circRNAs in NP cells may play an important role in the process of IDD. To identify circular RNAs (circRNAs) associated with human IDD, we isolated the NP cells from human degenerated and non-degenerated intervertebral disc and identified NP cells by microscopy and cell proliferation. CircRNA microarray expression profiles were obtained from NP cells of degenerated and non-degenerated intervertebral disc and further validated by quantitative reverse transcription PCR (qRT-PCR). The expression data were analyzed by bioinformatics. Microarray analysis identified 7294 circRNAs differentially expressed in degenerated human IDD NP cells. Among them, 3724 circRNAs were up-regulated and 3570 circRNAs were down-regulated by more than 2 folds. After validating by qRT-PCR, we predicted the possible miRNAs of the top dysregulated circRNAs using TargetScan, and miRanda. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the most modulated circRNAs regulate the viability, degradation, apoptosis and oxidative stress in NP cells, and the possible mechanism underlying IDD was discussed. These results revealed that circRNAs may play a role in IDD and might be a promising candidate molecular target for gene therapy.

Download Full-text