Skip to main content

Genome-wide analysis of circular RNAs and validation of hsa_circ_0086354 as a promising biomarker for early diagnosis of cerebral palsy

Abstract

Background

Cerebral palsy (CP) is a spectrum of non-progressive motor disorders caused by brain injury during fetal or postnatal periods. Current diagnosis of CP mainly relies on neuroimaging and motor assessment. Here, we aimed to explore novel biomarkers for early diagnosis of CP.

Methods

Blood plasma from five children with CP and their healthy twin brothers/sisters was analyzed by gene microarray to screen out differentially expressed RNAs. Selected differentially expressed circular RNAs (circRNAs) were further validated using quantitative real-time PCR. Receiver operating characteristic (ROC) curve analysis was used to assess the specificity and sensitivity of hsa_circ_0086354 in discriminating children with CP and healthy controls.

Results

43 up-regulated circRNAs and 2 down-regulated circRNAs were obtained by difference analysis (fold change > 2, p < 0.05), among which five circRNAs related to neuron differentiation and neurogenesis were chosen for further validation. Additional 30 pairs of children with CP and healthy controls were recruited and five selected circRNAs were further detected, showing that hsa_circ_0086354 was significantly down-regulated in CP plasma compared with control, which was highly in accord with microarray analysis. ROC curve analysis showed that the area under curve (AUC) to discriminate children with CP and healthy controls using hsa_circ_0086354 was 0.967, the sensitivity was 0.833 and the specificity was 0.966. Moreover, hsa_circ_0086354 was predicted as a competitive endogenous RNA for miR-181a, and hsa_circ_0086354 expression was negatively correlated to miR-181a expression in children with CP.

Conclusion

Hsa_circ_0086354 was significantly down-regulated in blood plasma of children with CP, which may be a novel competent biomarker for early diagnosis of CP.

Peer Review reports

Background

Since W.J. Little first described in the 1840s, the concept of cerebral palsy (CP) has been revised for several times and is now defined as a non-progressive motor disorder induced by brain injury during prenatal (80%), perinatal (10%) or postnatal (10%) [1, 2]. The incidence of CP is 1.25 per 1000 neonates in China and 2–3 per 1000 neonates worldwide [3, 4]. The brain injury in children with CP results in activity limitation in most cases, accompanying with impaired communication and cognition [5, 6]. To date, CP has no cure and would cost millions of healthcare expenditure, making CP as a severe public health problem that brings enormous burden for patient families [2, 7]. Preterm birth and asphyxia result from dystocia are the most common risk factors for CP [8, 9]. Administration of magnesium sulfate for women at risks of premature delivery and cooling therapy for infants at high risks of CP are considered to be effective preventive methods [10,11,12]. Unfortunately, existing diagnosis by comprehensive analysis of neonatal encephalopathy history, neuroimaging and neurodevelopmental assessment is limited and needs further researches [4, 13]. Therefore, a better understanding of CP aetiology and pursuit of more accurate early diagnostic methods are of great importance.

Noncoding RNAs represents more than 98% of all human transcripts, among which circular RNAs (circRNAs) are a special subtype without 5′ cap or 3′ poly-A tail [14, 15]. circRNAs become a new research hotspot in the past decade owing to their diverse physiological functions: circRNAs sponge microRNAs according to the “competing endogenous RNA” (ceRNA) [16]; they also act as protein scaffolds or templates for protein translation [17, 18]. Besides, increasing evidence indicates that circRNAs are implicated in the regulation of various human diseases including cardiovascular diseases, cancers and neurological diseases [19,20,21]. circRNAs may also serve as potent biomarkers for the early detection of specific diseases attributing to its stability and easy accessibility [22, 23].

With the rapid development of next-generation sequencing, over 1000 circRNAs in human serum exosomes were identified [23, 24]. In present study, we screened out differential expressed circRNAs between children with CP and their healthy controls using microarray technology, to select novel biomarkers for early diagnosis and intervention of CP as well as provide a better understanding of CP etiology.

Methods

Sample preparation

Five children with CP and their healthy twins were selected in our study to minimize individual differences (Average age: 3.3 ± 1.5, average birth weight: 2.9 ± 0.4 kg). Detailed clinical information of participants was provided in Additional file 1: Table S1. Additional 30 pairs of children with CP and healthy controls (without any congenital or acquired disease) were recruited for subsequent validation of differential expressed circRNAs (Age: 4.2 ± 1.6; average birth weight: 3.1 ± 0.6 kg). The diagnostic criterion for CP were the combination of following clinical findings: motor dysfunction assessed using the Hammersmith Infant Neurological Examination (HINE); abnormal neuroimaging detected by magnetic resonance imaging (MRI); comprehensive assessment of clinical history and high risks for CP including prematurity and low birthweight [13, 25]. Inclusive criteria: undergo no drug therapy; with complete clinical information; Exclusive criteria: acute/chronic infectious diseases, connective tissue diseases or malignant tumor; recent use of immunosuppressant; injury of liver and kidney function. Whole blood sample (5 ml per patient) was collected in Na2EDTA tubes. Plasma was isolated by centrifugation, followed by total RNA extraction using TRIzol reagent (ThermoFisher Scientific, Waltham, MA, USA). All blood samples were collected with the consent of parents of children with CP. And all experiments performed in this study were in accord with the ethical guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Xi’an International Medical Center Hospital.

Microarray analysis

After RNA integrity assessment using Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, USA), total RNAs were reversely transcribed into cDNA, which was further used to generate biotinylated cRNAs. Then cRNAs were hybridized with Hybridization Slides (Agilent technologies, Santa Clara, CA, USA) in a Hybridization Oven at 65 °C for 17 h. Sides were scanned under a Microarray Scanner (Agilent technologies, Santa Clara, CA, USA) and raw data were obtained by the Feature Extraction software 10.7 (Agilent technologies, Santa Clara, CA, USA), followed by raw data normalization using Quantile algorithm. cirRNAs with a fold change > 2, p value < 0.05 were presented as heatmap plots using R package “pheatmap”. Then differential expressed circRNAs with flag-signal of “Absent” in CP group or healthy controls were removed. Gene Ontology (GO) enrichment analysis were performed use Fisher's exact test by a R package “clusterProfiler” of the target genes. For CP etiology and biomarker investigations, five circRNAs regarding neuron differentiation and neurogenesis were chosen for further quantitative real-time PCR verification. The TargetScan prediction tool was used to identify interactions between hsa_circ_0086354 and target miRNAs. miRNAs that had perfect nucleotide pairing with hsa_circ_0086354 were selected. Further Pearson correlation was carried out to analyze the correlation between hsa_circ_0086354 and miRNAs, only interactions with significant negative correlation was retained. The circRNAs-miRNAs network was visualized by Cytoscape software (version 3.7.0; http://www.cytoscape.org) [26].

Quantitative real-time PCR

Additional thirty pairs of children with CP and their healthy controls were recruited to verify the differential expressed circRNAs screened by the microarray. In brief, total RNAs of plasma were extracted using UNIQ-10 RNA extraction kit (Sangon Biotech, Shanghai, China) and reversely transcribed into cDNA using Maxima Reverse Transcriptase (ThermoFisher Scientific, Waltham, MA, USA). Then cDNAs were quantified using Fast qPCR Master Mix (High Rox) (Sangon Biotech, Shanghai, China) in an ABI Stepone plus PCR instrument. Similar methods were used to detect miR-181a level. 18S ribosomal RNA was used as internal control for hsa-circRNAs and RNU6B was used as an internal control for miR-181a. All data were analyzed using the 2 method. Specific primers used for circRNAs detection were listed in Table 1.

Table 1 Primers used for quantitative real-time PCR in this study

Statistical analysis

Data from quantitative real-time PCR was analyzed using the 2 method and mean values were compared using unpaired t-test (Graphpad Prism 8.0, USA). All experiments were repeated for at least three times and p values less than 0.05 were regarded as statistically significant.

Results

Differentially expressed circRNAs screened by microarray

In order to minimize individual differences, the blood samples from five pairs of twin children (one CP and one healthy) were collected in our study. Sino human ceRNA array V3.0 which includes 53,625 human circRNAs was used to screen out differentially expressed circRNAs between the twins. Volcano plot showed that 134 circRNAs were differentially expressed in children with CP compared to their healthy controls, among which 77 circRNAs were up-regulated and 57 were down-regulated (fold change > 2, p < 0.05) (Fig. 1A). After further flag-signal screening, 45 differentially expressed human circRNAs were obtained (fold change > 2, p < 0.05). As listed in Table 2, 43 circRNAs were up-regulated and 2 circRNAs were down-regulated. According to GO enrichment analysis, we selected five differentially expressed circRNAs that were mainly involved in neuron differentiation and neurogenesis for further quantitative real-time PCR verification (Additional file 2: Fig. S1). As clustering analysis of heatmap shows, hsa_circ_0042123 (host gene: peripheral myelin protein 22 (PMP22)), hsa_circ_0083264 (host gene: Rho guanine nucleotide exchange factor 10 (ARHGEF10)), hsa_circ_0035127 (host gene: myelin expression factor 2 (MYEF2)) and hsa_circ_0015069 (host gene: PBX homeobox 1 (PBX1)) were up-regulated in CP group while hsa_circ_0086354 (host gene: protein tyrosine phosphatase receptor type D (PTPRD)) was down-regulated versus the control (Fig. 1B).

Fig. 1
figure 1

Differentially expressed circRNAs between children with CP and their healthy twins. A Differentially expressed circRNAs between children with CP and their healthy controls was shown in Volcano plot. “Red” represents up-regulated circRNAs, “Blue” represents down-regulated circRNAs (fold change > 2, p < 0.05). B Heatmap clustering analysis was performed to display 5 selected circRNAs. Rows represent differential circRNAs and columns represent five pairs of CP samples and healthy controls. “Green” represents down-regulation and “Red” represents up-regulation of circRNAs in each sample. CP cerebral palsy

Table 2 Microarray analysis of differential expressed circRNAs in 5 CP children compared with their healthy twins

Hsa_circ_0086354 is a potential biomarker for early diagnosis of CP

Further quantitative real-time PCR validation showed that the fold changes of CP versus Control were as follow: hsa_circ_0042123 was − 2.067 (microarray: 2.003), hsa_circ_0083264 was − 1.031 (microarray: 2.039), hsa_circ_0035127 was − 1.408 (microarray: 2.144), hsa_circ_0015069 was 1.76 (microarray: 2.113) and hsa_circ_0086354 was − 6.15 (microarray: − 3.676) (Fig. 2A). The expression pattern of hsa_circ_0086354 validated by real-time PCR was highly in accord with that detected by microarray, showing that hsa_circ_0086354 was significantly down-regulated in CP group (Fig. 2B). Further receiver operating characteristic (ROC) curve analysis showed that the area under the curve (AUC) to discriminate CP and healthy controls using hsa_circ_0086354 level was 0.967, the sensitivity was 0.833 and the specificity was 0.966 (Fig. 2C), suggesting that hsa_circ_0086354 is a potential biomarker for CP diagnosis.

Fig. 2
figure 2

Hsa_circ_0086354 is a potential biomarker for early diagnosis of CP. Five circRNAs screened out by microarray was further validated by quantitative real-time PCR. A All expression levels in Control group were normalized to “1”. Fold changes of CP versus Control obtained from microarray and real-time PCR were shown. B Relative expression of hsa_circ_0086354 in 30 pairs of CP and Control samples detected by quantitative real-time PCR was shown. C ROC curve analysis was carried out to assess the value of hsa_circ_0086354 in discriminating children with CP and healthy controls. The AUC was 0.967, the sensitivity was 0.833 and the specificity was 0.966. “ns”: not significant; ***p < 0.001. CP cerebral palsy, ROC receiver operating characteristic, AUC area under the curve

miR-181a is a downstream target of hsa_circ_0086354 in CP

Hsa_circ_0086354 associated ceRNA network was obtained using Cytoscape analysis. miR-181a, miR-4741 and miR-4656 were down-stream target microRNAs of hsa_circ_0086354 (Fig. 3A). Further quantitative real-time PCR assay showed that miR-181a level was significantly up-regulated in children with CP (Fig. 3B). Besides, the miR-181a level was negatively correlated to hsa_circ_0086354 level in children with CP (Fig. 3C). These results implied that miR-181a is a downstream target of hsa_circ_0086354 in CP.

Fig. 3
figure 3

miR-181a is a downstream target of hsa_circ_0086354 in CP. A Cytoscape analysis was performed to show the hsa_circ_0086354 associated ceRNA network. B The expression of miR-181a was detected using quantitative real-time PCR. C The correlation between hsa_circ_0086354 expression and miR-181a expression was analyzed using Graphpad Prism 8.0. ceRNA competitive endogenous RNA, CP cerebral palsy

Discussion

Owing to its enigmatic etiology, the diagnosis of CP can barely rely on neuroimaging and assessment of motor dysfunction [27]. CirRNAs were first considered as byproducts of mis-splicing, yet increasing evidence indicated that circRNAs are implicated in various molecular processes as well as human diseases: circRNAs regulate gene expression via regulating gene transcription, gene splicing or sponging microRNAs; circRNAs are involved in the regulation of neuronal diseases, cardiovascular disease and cancer progression. Of note, ciRS-7 regulates α-synuclein expression through co-expressing and co-localizing with miR-7 to further regulate brain development [24]. Besides, majority of identified circRNAs are abundantly detected in brain tissues and neurons, which inspired us to explore specific biomarkers for CP diagnosis.

In the present study, blood samples from five children with CP and their twin brothers/sisters were collected to screen out differentially expressed circRNAs using microarray. Twin participants at identical preterm conditions can exclude additional risk factors of CP, which makes our results more reliable. Five circRNAs enriched in neuron differentiation and neurogenesis were selected from 45 differentially expressed circRNAs for further validation. Another 30 pairs of plasma samples from children with CP and healthy controls were collected, and the expression levels of five selected circRNAs were quantified. It was remarkable that the expression pattern of hsa_circ_0086354 measured by quantitative real-time PCR was highly in consistent with that detected by microarray. Yet the expression differences between children with CP and healthy controls of hsa_circ_0042123, hsa_circ_0083264, hsa_circ_0035127 and hsa_circ_0015069 were either not significant or contradictory with microarray analysis. Therefore, our findings suggest that hsa_circ_0086354 might serve as a promising biomarker for CP diagnosis.

circRNAs have been reported to serve as competent biomarkers for diagnosis of various diseases. For instance, plasma hsa_circRNA_002453 was a potential biomarker for severity of renal involvement and diagnosis of lupus nephritis with an AUC of 0.906 [28]. Hsa_circRNA_0000520 is remarkably down-regulated in gastric cancer and may serve as a potential biomarker for early diagnosis [29]. Hsa_circRNA_0001649 is a novel specific biomarker for colorectal cancer assessment [30]. circRNAs display high stability owing to their covalent loop structure, which helps them get rid of de-adenylation, de-capping and RNases degradation. The tissue-specific expression pattern of circRNAs enables them to serve as specific biomarkers for specific diseases [31, 32]. The application of circRNAs as biomarkers has always been a controversial topic, and the abundance of circRNAs is the major concern. Indeed, generally, the abundance of circRNAs is relatively low compared to their linear RNA product in body fluids. However, others demonstrated that some circRNAs are detected at comparable, even higher expression to their linear RNA [33, 34]. Besides, the rapid development of next-generation sequencing will provide substantial technical support for circRNA detection. Dong, R concluded that majority of annotated circRNAs are identified in brain tissues and neurons [35]. In the present study, hsa_circ_0086354 was significantly down-regulated in CP plasma with an AUC of 0.967, suggesting hsa_circ_0086354 may be a promising biomarker for the early diagnosis of CP. In addition, the host gene of hsa_circ_0086354 is PTPRD, which is highly expressed in brain tissues and regulated neurite growth and neurons axon guidance, indicating that PTPRD and hsa_circ_0086354 might involve in CP etiology [36, 37].

We further discovered that hsa_circ_0086354 acts as a ceRNA of miR-181a. miR-181a is up-regulated in patients with mild cognitive impairment which later progressed to Alzheimer’s disease [38]. miR-181a is also up-regulated in rats after ischemia/reperfusion induced cerebral injury [39]. On the contrary, miR-181a silencing exerts neuroprotective effects through suppressing neuronal apoptosis and neuronal loss both in a rat model and in epilepsy children [40, 41]. MiR-181a silencing also promotes neuronal growth via regulating the Smad signaling in Parkinson’s disease [42]. Besides, miR-181a contributes to neural stem cell differentiation and promotes generation of neurons [43, 44]. Here we found that miR-181a was significantly up-regulated in children with CP, and miR-181a level was negatively correlated to hsa_circ_0086354 level. All these findings imply that hsa_circ_0086354 might be involved in the regulation of neuronal survival and neuronal differentiation through targeting miR-181a.

Conclusion

Hsa_circ_0086354 is significantly down-regulated in children with CP in contrast with their healthy control with an AUC of 0.967, making it as a promising biomarker for the early diagnosis of CP. Hsa_circ_0086354 may also be involved in the etiology of CP through targeting miR-181a.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the Gene Expression Omnibus (GEO) database at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE183021 (accession number: GSE183021).

Abbreviations

CP:

Cerebral palsy

circRNAs:

Circular RNAs

ceRNA”:

Competing endogenous RNA”

ROC:

Receiver operating characteristic

AUC:

Under the curve

PTPRD:

Protein tyrosine phosphatase receptor type D

References

  1. Little WJ. The classic: hospital for the cure of deformities: course of lectures on the deformities of the human frame. Clin Orthop Relat Res. 2012;470(5):1252–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Sewell MD, Eastwood DM, Wimalasundera N. Managing common symptoms of cerebral palsy in children. BMJ. 2014;349:g5474.

    Article  PubMed  Google Scholar 

  3. He P, Chen G, Wang Z, Guo C, Zheng X. Children with motor impairment related to cerebral palsy: prevalence, severity and concurrent impairments in China. J Paediatr Child Health. 2017;53(5):480–4.

    Article  PubMed  Google Scholar 

  4. Korzeniewski SJ, Slaughter J, Lenski M, Haak P, Paneth N. The complex aetiology of cerebral palsy. Nat Rev Neurol. 2018;14(9):528–43.

    Article  PubMed  Google Scholar 

  5. Sadowska M, Sarecka-Hujar B, Kopyta I. Cerebral palsy: current opinions on definition, epidemiology, risk factors, classification and treatment options. Neuropsychiatr Dis Treat. 2020;16:1505–18.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Pakula AT, Van Naarden Braun K, Yeargin-Allsopp M. Cerebral palsy: classification and epidemiology. Phys Med Rehabil Clin N Am. 2009;20(3):425–52.

    Article  PubMed  Google Scholar 

  7. Kancherla V, Amendah DD, Grosse SD, Yeargin-Allsopp M, Van Naarden Braun K. Medical expenditures attributable to cerebral palsy and intellectual disability among Medicaid-enrolled children. Res Dev Disabil. 2012;33(3):832–40.

    Article  PubMed  Google Scholar 

  8. Frank R, Garfinkle J, Oskoui M, Shevell MI. Clinical profile of children with cerebral palsy born term compared with late- and post-term: a retrospective cohort study. BJOG. 2017;124(11):1738–45.

    Article  CAS  PubMed  Google Scholar 

  9. Moster D, Wilcox AJ, Vollset SE, Markestad T, Lie RT. Cerebral palsy among term and postterm births. JAMA. 2010;304(9):976–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Nelson KB, Grether JK. Can magnesium sulfate reduce the risk of cerebral palsy in very low birthweight infants? Pediatrics. 1995;95(2):263–9.

    Article  CAS  PubMed  Google Scholar 

  11. Schendel DE, Berg CJ, Yeargin-Allsopp M, Boyle CA, Decoufle P. Prenatal magnesium sulfate exposure and the risk for cerebral palsy or mental retardation among very low-birth-weight children aged 3 to 5 years. JAMA. 1996;276(22):1805–10.

    Article  CAS  PubMed  Google Scholar 

  12. Rao R, Trivedi S, Vesoulis Z, Liao SM, Smyser CD, Mathur AM. Safety and short-term outcomes of therapeutic hypothermia in preterm neonates 34–35 weeks gestational age with hypoxic-ischemic encephalopathy. J Pediatr. 2017;183:37–42.

    Article  PubMed  Google Scholar 

  13. Novak I, Morgan C, Adde L, Blackman J, Boyd RN, Brunstrom-Hernandez J, et al. Early, accurate diagnosis and early intervention in cerebral palsy: advances in diagnosis and treatment. JAMA Pediatr. 2017;171(9):897–907.

    Article  PubMed  Google Scholar 

  14. Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32(5):453–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447(7146):799–816.

    Article  CAS  PubMed  Google Scholar 

  16. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014;56(1):55–66.

    Article  CAS  PubMed  Google Scholar 

  18. Pamudurti NR, Bartok O, Jens M, Ashwal-Fluss R, Stottmeister C, Ruhe L, et al. Translation of CircRNAs. Mol Cell. 2017;66(1):9-21 e7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Altesha MA, Ni T, Khan A, Liu K, Zheng X. Circular RNA in cardiovascular disease. J Cell Physiol. 2019;234(5):5588–600.

    Article  CAS  PubMed  Google Scholar 

  20. Li Y, Ge YZ, Xu L, Jia R. Circular RNA ITCH: a novel tumor suppressor in multiple cancers. Life Sci. 2020;254:117176.

    Article  CAS  PubMed  Google Scholar 

  21. Li TR, Jia YJ, Wang Q, Shao XQ, Lv RJ. Circular RNA: a new star in neurological diseases. Int J Neurosci. 2017;127(8):726–34.

    Article  CAS  PubMed  Google Scholar 

  22. Lei B, Tian Z, Fan W, Ni B. Circular RNA: a novel biomarker and therapeutic target for human cancers. Int J Med Sci. 2019;16(2):292–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Li Y, Zheng Q, Bao C, Li S, Guo W, Zhao J, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. 2015;25(8):981–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Hsiao KY, Sun HS, Tsai SJ. Circular RNA—new member of noncoding RNA with novel functions. Exp Biol Med (Maywood). 2017;242(11):1136–41.

    Article  CAS  Google Scholar 

  25. Patel DR, Neelakantan M, Pandher K, Merrick J. Cerebral palsy in children: a clinical overview. Transl Pediatr. 2020;9(Suppl 1):S125–35.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Li Y, Sun R, Li R, Chen Y, Du H. Prognostic nomogram based on circular RNA-associated competing endogenous RNA network for patients with lung adenocarcinoma. Oxid Med Cell Longev. 2021;2021:9978206.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ashwal S, Russman BS, Blasco PA, Miller G, Sandler A, Shevell M, et al. Practice parameter: diagnostic assessment of the child with cerebral palsy: report of the Quality Standards Subcommittee of the American Academy of Neurology and the Practice Committee of the Child Neurology Society. Neurology. 2004;62(6):851–63.

    Article  CAS  PubMed  Google Scholar 

  28. Ouyang Q, Huang Q, Jiang Z, Zhao J, Shi GP, Yang M. Using plasma circRNA_002453 as a novel biomarker in the diagnosis of lupus nephritis. Mol Immunol. 2018;101:531–8.

    Article  CAS  PubMed  Google Scholar 

  29. Sun H, Tang W, Rong D, Jin H, Fu K, Zhang W, et al. Hsa_circ_0000520, a potential new circular RNA biomarker, is involved in gastric carcinoma. Cancer Biomark. 2018;21(2):299–306.

    Article  CAS  PubMed  Google Scholar 

  30. Ji W, Qiu C, Wang M, Mao N, Wu S, Dai Y. Hsa_circ_0001649: a circular RNA and potential novel biomarker for colorectal cancer. Biochem Biophys Res Commun. 2018;497(1):122–6.

    Article  CAS  PubMed  Google Scholar 

  31. Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE. 2012;7(2):e30733.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Salzman J, Chen RE, Olsen MN, Wang PL, Brown PO. Cell-type specific features of circular RNA expression. PLoS Genet. 2013;9(9):e1003777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Salzman J. Circular RNA expression: its potential regulation and function. Trends Genet. 2016;32(5):309–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19(2):141–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Dong R, Ma XK, Chen LL, Yang L. Increased complexity of circRNA expression during species evolution. RNA Biol. 2017;14(8):1064–74.

    Article  PubMed  Google Scholar 

  36. Tomita H, Cornejo F, Aranda-Pino B, Woodard CL, Rioseco CC, Neel BG, et al. The protein tyrosine phosphatase receptor delta regulates developmental neurogenesis. Cell Rep. 2020;30(1):215–28.

    Article  CAS  PubMed  Google Scholar 

  37. Uhl GR, Martinez MJ. PTPRD: neurobiology, genetics, and initial pharmacology of a pleiotropic contributor to brain phenotypes. Ann N Y Acad Sci. 2019;1451(1):112–29.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ansari A, Maffioletti E, Milanesi E, Marizzoni M, Frisoni GB, Blin O, et al. miR-146a and miR-181a are involved in the progression of mild cognitive impairment to Alzheimer’s disease. Neurobiol Aging. 2019;82:102–9.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang Y, Shan Z, Zhao Y, Ai Y. Sevoflurane prevents miR-181a-induced cerebral ischemia/reperfusion injury. Chem Biol Interact. 2019;308:332–8.

    Article  CAS  PubMed  Google Scholar 

  40. Ren L, Zhu R, Li X. Silencing miR-181a produces neuroprotection against hippocampus neuron cell apoptosis post-status epilepticus in a rat model and in children with temporal lobe epilepsy. Genet Mol Res. 2016;15(1):1–11.

    Article  Google Scholar 

  41. Moon JM, Xu L, Giffard RG. Inhibition of microRNA-181 reduces forebrain ischemia-induced neuronal loss. J Cereb Blood Flow Metab. 2013;33(12):1976–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Hegarty SV, Sullivan AM, O’Keeffe GW. Inhibition of miR-181a promotes midbrain neuronal growth through a Smad1/5-dependent mechanism: implications for Parkinson’s disease. Neuronal Signal. 2018;2(1):NS20170181.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Stappert L, Borghese L, Roese-Koerner B, Weinhold S, Koch P, Terstegge S, et al. MicroRNA-based promotion of human neuronal differentiation and subtype specification. PLoS ONE. 2013;8(3):e59011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Hamada N, Fujita Y, Kojima T, Kitamoto A, Akao Y, Nozawa Y, et al. MicroRNA expression profiling of NGF-treated PC12 cells revealed a critical role for miR-221 in neuronal differentiation. Neurochem Int. 2012;60(8):743–50.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: SS and XB; methodology: YH; formal analysis: YH; Sample collection: YH, CW and YW; writing—original draft preparation: XB; writing—review and editing: YW, XG and SZ; All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Shiquan Sun.

Ethics declarations

Ethics approval and consent to participate

All experiments performed in this study were in accord with the ethical guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Xijing University. All blood samples were collected with the informed consent of legal guardians of children with CP.

Consent for publication

Written consent to publish this information was obtained from legal guardians of participants by the study coordinators. All legal guardians were informed about sample collection, data publication and personal information confidentiality.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Relative clinical information of children with cerebral palsy and their healthy controls.

Additional file 2: Fig. S1.

Top 30 of biological_process, cellular_component and molecular_function obtained using Gene Ontology enrichment. Plot size refers to gene number.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, Y., Bian, X., Wu, C. et al. Genome-wide analysis of circular RNAs and validation of hsa_circ_0086354 as a promising biomarker for early diagnosis of cerebral palsy. BMC Med Genomics 15, 13 (2022). https://doi.org/10.1186/s12920-022-01163-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-022-01163-6

Keywords