Decoding competitive endogenous RNA regulatory network in postoperative cognitive dysfunction

Abstract

Postoperative cognitive dysfunction (POCD) is a common postoperative neurological complication in elderly patients. Circular RNAs (circRNAs) are abundant in the mammalian brain and can probably regulate cognitive function. However, the competitive endogenous RNA (ceRNA) regulatory network in POCD remains illiterate. Transcriptomic signatures in the hippocampus of POCD mice derived from the Gene Expression Omnibus (GEO) dataset GSE190880, GSE95070, and GSE115440 were used to identify the circRNA, miRNA, and mRNA expression profiles of POCD mice compared with controls, respectively. A set of differentially expressed RNAs, including 119 circRNAs, 33 miRNAs, and 49 mRNAs were identified. Transcript validation showed the enhanced expression of circ_0001634, circ_0001345, and circ_0001493. A ceRNA regulatory network composed of three circRNAs, three miRNAs, and six mRNAs was established. The hub mRNAs in the ceRNA network were further found to be involved in the hormone catabolic process and regulation of canonical Wnt signaling pathway, revealing their crucial role in POCD. Finally, three miRNAs and four mRNAs were verified by qRT-PCR. These results based on bioinformatics and PCR array suggest that circ_0001634/miR-490-5p/Rbm47, circ_0001634/miR-490-5p/Sostdc1, circ_0001634/miR-7001-5p/Sostdc1, circ_0001345/miR-7001-5p/Sostdc1, and circ_0001493/miR-7001-5p/Sostdc1 may be novel diagnostic biomarkers and therapeutic targets for POCD.

Introduction

Postoperative cognitive dysfunction (POCD) is a neurocognitive disorder such as acute or persistent impairments in attention, learning, memory, and information processing that occurs predominantly in geriatric patients who undergo anesthesia and major surgery (Suwanabol et al., 2022). The long-term impact of POCD is associated with a worse overall quality of life, increased mortality, and a heavy burden on society and families (Boone et al., 2020; Tang et al., 2020). Recent evidence suggests that various elements, such as aberrant expression of apolipoprotein E4 genotype, blood-brain barrier (BBB) compromise, surgery-induced neuroinflammation, microglial activation, mitophagy impairment, and iron accumulation (McDonagh et al., 2010; Chen et al., 2020; Danielson et al., 2020; Wu et al., 2020; Yang et al., 2020), are intimately involved in all stages of postoperative cognitive decline. Despite major advances in comprehending pathogenesis and prevention methods, the overall curative effect remains inadequate. Consequently, it is imperative to screen promising therapeutic targets for POCD.

Circular RNA (circRNA), unlike conventional linear RNA, is a subtype of endogenous non-coding RNAs with covalently closed-loop structures that confer high resistance to degradation (Chen, 2020). CircRNAs have been proposed to modulate gene transcription, to act as microRNA (miRNA) or RNA-binding protein (RBP) decoys, and to function as protein scaffolds (Kristensen et al., 2019). However, the overwhelming majority of circRNAs are thought to act as competing endogenous RNAs (ceRNAs), which specify that circRNAs can compete with mRNAs for binding to the shared miRNAs and thereby indirectly modify gene expression at the post-transcriptional level (Chen, 2020). The abundance of circRNAs in the mammalian brain, particularly synaptoneurosomes, is highlighted by RNA sequencing data (Rybak-Wolf et al., 2015). Mounting evidence also suggests circRNAs play a vital role in cognitive impairments such as Alzheimer impisease (AD) (Dube et al., 2019). Downregulation of circCwc27 improved cognitive capacity in the AD mouse (Song et al., 2022). Knockdown of circTshz2-2 alleviated obesity-induced spatial memory decline via modulating BDNF/TrkB signaling pathway (Yoon et al., 2021). CircPtk2 contributed to sepsis-provoked cognitive dysfunction by serving as a sponge of miR-181c-5p to facilitate HMGB1 expression (Li et al., 2021). Niu et al. (2021) proposed that aerobic exercise mitigated vascular cognitive impairment by activating circRIMS2/miR-186/BDNF axis. Additionally, plasma circRNA-089763 is positively correlated with the occurrence of POCD (Wang et al., 2019; Zhou et al., 2020). A microarray analysis screened 210 differentially expressed circRNAs in POCD patients’ serum, such as circCPNE1, circUBE3B, and circITSN1 (Gao et al., 2020). Another microarray profiling highlighted three circRNAs (circ_22058, circ_44122, and circ_22673) as key elements in POCD (Wu et al., 2021). Circ_009789- and circ_004229-associated ceRNA networks were identified to elucidate the mechanism underlying susceptibility to POCD in aged mice (Zhang et al., 2022). CircShank3 may be involved in dexmedetomidine-mediated protection against POCD via targeting the p53 and NF-κB signaling pathways (Cao et al., 2020). Nevertheless, the role of circRNAs in POCD remains elusive.

MiRNAs are small (18–22 nucleotide long) non-coding RNAs that trigger post-transcriptional repression of gene expression by directly binding to the 3′-untranslated region (UTR) of mRNAs, profoundly governing neurodevelopment and neurodegeneration (Singh and Yadav, 2020). Multiple aberrantly expressed miRNAs, such as miR-124, miR-146a, and miR-381, were implicated in the development of POCD (Chen et al., 2019a,b; Wang et al., 2021). However, the underlying mechanisms of miRNAs in the neuropathogenesis of POCD are still unclear.

In this study, we extracted the differentially expressed circRNAs, miRNAs, and mRNAs in POCD from Gene Expression Omnibus (GEO) databases. According to the flowchart diagram (Figure 1), circRNA-miRNA pairs and miRNA-mRNA pairs were successively identified for constructing a ceRNA regulatory network. Furthermore, functional enrichment analysis and transcript validation were implemented to identify the hub genes and interpret potential regulatory mechanisms in the development of POCD. The present findings may strengthen our understanding of POCD.

FIGURE 1

Materials and methods

Expression microarray data

The circRNA dataset (GSE190880) was retrieved from the GPL21826 platform and contained 3 pairs of hippocampus from POCD and control mice (Ran et al., 2022). The miRNA dataset (GSE95070) was retrieved from the GPL19117 platform and contained 5 pairs of hippocampus from POCD and control mice (Wei et al., 2017). Meanwhile, we also downloaded the mRNA expression profiles from GSE115440 (GPL11533) (Yang et al., 2019), including 3 pairs of hippocampus tissues from POCD and control mice, to further construct the ceRNA network about POCD. The profiles of these three microarray datasets (GSE190880, GSE95070, and GSE115440) regarding circRNAs, miRNAs, and mRNAs were listed in Table 1. In these three microarray datasets, aseptic tibia fracture in C57BL/6 male mice was applied to mimic POCD.

Identification of the differentially expressed circular RNAs, microRNAs, and mRNAs

Differentially expressed (DE) analysis between POCD and control hippocampus samples were performed using GEO2R (Barrett et al., 2013), an online program based on R language that can analyze any GEO series to compare two groups of data. The cut-off standard of DEcircRNAs was set to p < 0.05 and |logFC| > 1.5, the cut-off standard of DEmiRNAs was set to p < 0.05 and |logFC| > 0.5, and the cut-off standard of DEmRNAs was set to p < 0.05 and |logFC| > 1. Principal component analysis (PCA) and heatmap (linkage method: row clustering; distance measure: euclidean) were visualized using the Xiantao search tool.¹

Establishment of circular RNAs– microRNA–mRNA network

The circRNA identification (ID) of DEcircRNAs were first converted into the circRNA ID in the circBase (Glažar et al., 2014),² an online database that integrates thousands of public circRNAs across five species and their genomic profiles. Only genomic sequences in FASTA format of annotated circRNAs were further analyzed using the circMIR1.0 software, a forecasting instrument of circRNA-miRNA interactions based on miRanda and RNAhybrid that can visualize the binding sites of miRNA adsorption. Overlapping miRNAs between the target miRNAs and DEmiRNAs screened in the GSE95070 were considered candidate miRNAs. Moreover, miRDB (Chen and Wang, 2020)³ and Targetscan (McGeary et al., 2019) (version 8.0)⁴ databases were used to predict the target mRNAs of the candidate miRNAs. The Venn diagram drawn in the Xiantao search tool was utilized to overlap the miRNA-forecasted mRNAs and DEmRNAs in the GSE115440. Finally, the ceRNA network was constructed using Cytoscape (Otasek et al., 2019) (version 3.9.1), by integrating circRNA-miRNA pairs and miRNA-mRNA pairs.

Functional classifications and pathway enrichment analysis

The Gene Ontology (GO) project is an ontological annotation resource that describes gene product function based on multiple databases (Ashburner et al., 2000). Kyoto Encyclopedia of Genes and Genomes (KEGG)⁵ is an integrated tool that provides a new perspective on molecular-level functions, diseases and drugs (Kanehisa et al., 2016). The DAVID Knowledgebase⁶ is a bioinformatics data resource based on a single-linkage method for the agglomeration of millions of genes that can be entered into DAVID gene clusters (Huang da et al., 2009a,b). GO annotation and KEGG pathway enrichment analysis were implemented using the Xiantao search tool. A p-value of <0.05 was set as the cut-off point.

Animals

Twelve-month-old C57BL/6 male mice weighing 25–35 g were purchased from Hubei Provincial Center for Disease Control and Prevention. Mice with ad libitum access to food and water were housed in a colony room with a temperature of 22–25°C, a humidity of 50%, and a 12-h light/dark cycle. The animals were acclimated to the environment for 7 days before experimental manipulation. All experimental procedures were approved by the Laboratory Animal Ethics Committee in Renmin Hospital of Wuhan University (No. WDRM-20210709).

Animal model of postoperative cognitive dysfunction

The animals were randomly assigned to two groups: the POCD group (n = 6) and control group (n = 6). Intramedullary fixation for open tibial fracture under isoflurane anesthesia has been widely used as a POCD model (Yang et al., 2019). Briefly, mice were exposed to 3% isoflurane for the anesthesia induction, followed by 1.5% isoflurane in 100% oxygen for the maintenance. After shaving and disinfecting the animal, a skin incision was made on the lateral tibia, and a 0.3-mm intramedullary fixation pin was inserted into the medullary cavity. Next, an osteotomy was performed at the middle and distal third of the bone, and the wound was closed with 5-0 Vicryl suture (Ethicon, Somerville, NJ, USA). The entire procedure from the induction to the end of surgery lasted 30 min. A warm pad was utilized to keep the body temperature around 37°C throughout the surgery. Then, 2% lidocaine solution and 1% tetracaine hydrochloride mucilage was applied locally for postoperative analgesia twice daily until 3 days after surgery. The control group received 100% oxygen without surgery or anesthesia.

Open field test

The open field test (OFT) was employed to assess the locomotor activity of the mice 3 days after the surgery. The apparatus was an opaque plastic cube box with a side length of 45 cm. Each mouse was gently placed in the central area and left to move freely for 5 min. The total distance and time spent in the center area were recorded using SuperMaze software (XinRuan Information Technology, Shanghai, China).

Morris water maze

Spatial learning and memory was assessed using Morris water maze test 4 days after the surgery. The apparatus was a circular plastic pool of 120-cm diameter and 50-cm height surrounded by four curtains. The pool was equipped with a 6-cm-diameter hidden platform placed in the center of the fourth quadrant. Prior to the experiment, the pool was filled with opaque water at around 25°C to a depth of about 30 cm, approximately 1 cm above the platform. Each mouse received four daily trials with a 30-min intertrial interval for four consecutive days to find the platform, and was randomly placed at one of the four quadrants. Mice that failed to find the platform within 90 s were manually directed to the platform. The escape latency (time taken to reach the platform) in the four quadrants on the same day was averaged. The platform was removed on the fifth day to allow for probe trial. During the 90-s session, the total time spent in target quadrant and the frequency of crossing over the target quadrant were recorded. The swimming trajectory of the animals were recorded automatically via a video tracking system (XinRuan Information Technology, Shanghai, China).

RNA extraction and quantitative real-time polymerase chain reaction

Firstly, five circRNAs (circ_0001634, circ_0001345, circ_0001493, circ_0000487, and circ_0001468) whose host genes may be related to cognitive function were selected for quantitative real-time polymerase chain reaction (qRT-PCR). To verify the reliability and accuracy of the predicted ceRNA network, three miRNAs (miR-6912-5p, miR-490-5p, and miR-7001-5p), and four mRNAs (Rbm47, Sostdc1, Cdh3, and Sfrp5) were selected for qRT-PCR. Total RNA was extracted from the hippocampus tissue of the control and POCD group using TRIpure reagent (ELK, Wuhan, China). For circRNAs and mRNAs, cDNA was generated from template RNA using EntiLink™ Reverse Transcriptase Kit (ELK, Wuhan, China), and analyzed with EnTurbo™ SYBR Green PCR SuperMix (ELK, Wuhan, China) using StepOne™ Real-Time PCR system (Life Technologies, Carlsbad, CA, USA). The reaction procedure was as follows: 95°C for 1 min, followed by 40 cycles of 95°C for 15 s, 58°C for 20 s and 72°C for 45 s. For miRNAs, cDNA was generated using EntiLink™ 1st Strand cDNA Synthesis Kit (ELK, Wuhan, China), and analyzed with EnTurbo™ SYBR Green PCR SuperMix (ELK, Wuhan, China) using QuantStudio 6 Flex real-time PCR system (Life Technologies, Carlsbad, CA, USA). The reaction procedure was as follows: 95°C for 30 s, followed by 95°C for 10 s, 58°C for 30 s, and 72°C for 30 s. The relative expression levels of these RNAs were calculated via the 2^–ΔΔCt method. GAPDH was used as an internal control gene for circRNA and mRNA expression, while U6 was used as an internal control gene for miRNA expression. The primers for these RNAs are shown as Table 2.

Statistical analysis

Quantitative data were presented as the mean ± standard deviation (SD) or median (range). GraphPad Prism version 9.0 (GraphPad Software Inc., San Diego, California, USA) was used for statistical analyses. The statistical significance between the two groups was determined using Student’s t-test or non-parametric test. A p-value of <0.05 was considered statistically significant.

Results

Isoflurane plus orthopedic surgery caused postoperative cognitive dysfunction in mice

The open feld test showed no statistical difference in the total distance and duration of the central area between the two groups, indicating that anesthesia/surgery did not affect the locomotor activity of the mice (Figures 2A,B). To address whether isoflurane/surgery impaired cognitive function in mice, we performed Morris water maze to assess learning and memory following orthopedic surgery. During hidden platform test, the mice undergoing anesthesia/surgery showed significantly longer escape latency on day 2, 3, and 4 than the untreated control (Figure 2C). During the probe test, time spent in target quadrant and platform crossing times in the mice undergoing anesthesia/surgery were less than those in the control group (Figures 2D–F). These results suggested that isoflurane plus orthopedic surgery caused POCD.

FIGURE 2

Data preprocessing

The accuracy and reliability of microarray data were evaluated using PCA. By comparing the distribution patterns of all identified circRNAs (Figure 3A) and miRNAs (Figure 4A), we could completely separate POCD samples from control hippocampus.

FIGURE 3

FIGURE 4

Circular RNA expression profiles in postoperative cognitive dysfunction

In the GSE190880 dataset, a total of 119 DEcircRNAs with thresholds of |logFC| > 1.5 and p < 0.05 were detected in the hippocampus of POCD mice, of which 103 were up-regulated and 16 were down-regulated. Additionally, 76.47% of DEcircRNAs originated from exonic regions, while 8.40% originated from intronic regions. Sense-overlapping and intergenic circRNAs accounted for 13.45 and 1.68%, respectively (Figure 3B). Chromosomal distribution demonstrated that DEcircRNAs are highly abundant in chr1-19 and chrX, but absent in chr20 and chrY (Figure 3C).

GO and KEGG pathway analyses were used to categorize and annotate the host genes of DEcircRNAs in order to further characterize these in POCD. The enriched GO terms were mainly associated with locomotory behavior (gene ratio = 8/105, p = 2.06E-05), vesicle-mediated transport in synapse (gene ratio = 8/105, p = 1.64E-05), regulation of synaptic transmission, glutamatergic (gene ratio = 5/105, p = 2.63E-05), intrinsic component of synaptic membrane (gene ratio = 9/106, p = 1.69E-06), integral component of synaptic membrane (gene ratio = 7/106, p = 7.67E-05), Schaffer collateral-CA1 synapse (gene ratio = 6/106, p = 1.27E-05), guanyl-nucleotide exchange factor activity (gene ratio = 7/106, p = 3.92E-05), calcium ion transmembrane transporter activity (gene ratio = 6/106, p = 4.49E-05), and cation:cation antiporter activity (gene ratio = 3/106, p = 4.49E-05) (Figure 3D and Supplementary Table 1). Moreover, KEGG pathway enrichment analysis concluded that DEcircRNAs were principally enriched in retrograde endocannabinoid signaling (gene ratio = 5/43, p = 0.0007), dilated cardiomyopathy (gene ratio = 4/43, p = 0.0010), and glutamatergic synapse (gene ratio = 4/43, p = 0.0021) (Figure 3D and Supplementary Table 1). Ulteriorly, the DEcircRNAs were blasted by circBase and 16 annotated circRNAs were identified, of which 14 were up-regulated and 2 were down-regulated, as shown by the hierarchical clustering heatmap (Figure 3E and Table 3).

Validation of differentially expressed circRNAs

To validate the expression of DEcircRNAs, five circRNAs whose host genes are possibly related to cognitive function were selected for verification by qRT-PCR. The expression of mmu_circ_0001634 (circFam53b), mmu_circ_0001345 (circFbxl5), and circ_0001493 (circIqsec1) were found to be upregulated in the hippocampus of POCD mice, which was consistent with the sequencing results. There was no significant difference in the expression of mmu_circ_0000487 (circXrcc4) between the two groups. Contrary to the sequencing result, isoflurane/surgery obviously diminished mmu_circ_0001468 (circHipk2) expression (Figure 5). To unveil insights into the potential functional mechanisms of circFam53b, circFbxl5, and circIqsec1, we predicted 987 miRNAs interacting with circRNAs and the putative binding domains with the aid of circMIR1.0 software (Figure 4B).

FIGURE 5

MicroRNA expression profiles in postoperative cognitive dysfunction

A total of 33 significant DEmiRNAs with cut-off criteria of |logFC| > 0.5 and p < 0.05 were screened in the hippocampus of POCD mice. We identified 12 overlapping DEmiRNAs by integrating GSE95070 data and predicted results, of which five were upregulated and seven were downregulated (Figures 4B,C and Supplementary Table 2). According to the ceRNA hypothesis that circRNAs could compete with mRNAs for the same miRNAs, we predicted the mRNAs downstream of 12 miRNAs. The results reflected that 2,488 target genes of the 12 DEmiRNAs were acquired by the databases of miRDB and TargetScan. Furthermore, a total of 49 significant DEmRNAs were extracted in GSE115440. The Venn diagram revealed that seven mRNAs were shared by the predicted genes of the above-mentioned circRNA-targeted miRNAs and the DEmRNAs from GSE115440 (Figure 4D and Supplementary Table 3).

Construction of the competitive endogenous RNA regulatory network and experimental verification

A ceRNA regulatory network based on 3 circRNAs, 12 miRNAs, and 7 mRNAs was created to further investigate the mechanisms by which circRNAs and miRNAs affect the occurrence and development of POCD. Eventually, 10 circRNA-miRNA pairs and 7 miRNA-mRNA pairs were identified, which were composed of three circRNAs, three miRNAs, and six mRNAs (Figure 6A). The enriched GO terms were chiefly associated with the hormone catabolic process, regulation of canonical Wnt signaling pathway, sperm midpiece, catenin complex, dipeptidyl-peptidase activity, and chloride ion binding. Moreover, the enriched KEGG pathway were renin-angiotensin system and renin secretion (Figure 6B and Supplementary Table 4). To substantiate the potential interaction in the ceRNA, these three miRNAs and Wnt-related genes were selected for verification by qRT-PCR. The enhanced expression of miR-6912-5p, Rbm47 and Sostdc1 were observed in the POCD group, while the expression of miR-490-5p, miR-7001-5p, Cdh3, and Sfrp5 were significantly decreased in the POCD group (Figures 6C,D).

FIGURE 6

Discussion

In our study, we outlined comprehensive transcriptome profiles of POCD. A total of 119 circRNAs, 33 miRNAs, and 49 mRNAs were identified as differentially expressed. The circRNA-miRNA-mRNA triple regulatory network consisted of three circRNAs, three miRNAs, and six mRNAs. Biological process enrichment analysis revealed that the bulk of DEmRNAs in the ceRNA network were involved in the “hormone catabolic process” and “regulation of canonical Wnt signaling pathway.” Transcript validation suggested that circ_0001634/miR-490-5p/Rbm47, circ_0001634/miR-490-5p/Sostdc1, circ_0001634/miR-7001-5p/Sostdc1, circ_0001345/miR-7001-5p/Sostdc1, and circ_0001493/miR-7001-5p/Sostdc1 axis were likely to participate in the development of POCD.

The vast majority of circRNAs are derived from a single or multiple exons of known coding genes and thereby located in the cytoplasm (Kristensen et al., 2019). Our analysis showed that up to 76.47% DEcircRNAs were from exons, and 14 out of 16 annotated circRNAs were of exon origin, which implied that circRNAs probably function as competing endogenous RNAs in the pathogenesis of POCD because the ceRNA network merely exists in the cytoplasm (Chen, 2020). Enrichment analysis showed the corresponding parental or host genes of DEcircRNAs were implicated in synaptic plasticity and mainly located in the Schaffer collateral-CA1 synapse, suggesting circRNAs are extensively involved in the regulation of cognitive function.

Circ_0001493, circ_0000487, circ_0001468, and circ_0001345 are spliced from IQ motif and sec7 domain-containing protein 1 (Iqsec1), x-ray repair cross complementing 4 (Xrcc4), homeodomain-interacting protein kinase 2 (Hipk2), and F-box/LRR repeat protein 5 (Fbxl5), respectively, which may be linked to neurocognitive disorders (Zhang et al., 2013; Gerez et al., 2019; Liang et al., 2020; Briševac et al., 2021). Moreover, circ_0001634 originates through back-splicing events from Fam53b, which has been reported to regulate Wnt signal transduction by altering β-catenin nuclear localization (Kizil et al., 2014). Wnt signaling has been implicated in the modulation of synaptogenesis, long-term potentiation (LTP), and dendrite arborization (Narvaes and Furini, 2022). Intranuclear accumulation of β-catenin marks activation of canonical Wnt signaling that sequentially leads to enhance the transcription of Wnt target genes via the interaction between β-catenin and T-cell factor/lymphoid enhancer factor (TCF/LEF) (Narvaes and Furini, 2022). Hu et al. (2016) found that prolonged exposure to 3.6% sevoflurane could disrupt BBB components via suppressing Wnt/β-catenin/Annexin A1 pathway in brain microvascular endothelial cells, indicating that promoting β-catenin synthesis can alleviate POCD. Hence, these five circRNAs were selected for verification. Quantitative RT-PCR array revealed that the augmented levels of circ_0001634, mmu_circ_0001345, and circ_0001493 were consistent with the sequencing results.

Based on ceRNA hypothesis that a protein-coding RNA and a non-coding RNA compete for the same miRNA through the shared miRNA response elements (MREs), miRNA serves as a bridge for circRNA-induced translation and/or stabilization of the mRNAs (Chen, 2020). In our study, miR-6912-5p, miR-490-5p, and miR-7001-5p were screened as hub nodes to compete with three circRNAs for governing the expression of six mRNAs. GO classification and enrichment analysis uncovered that cadherin 3 (Cdh3), secreted frizzled related protein 5 (Sfrp5), and sclerostin domain containing 1 (Sostdc1) in the ceRNA network participate in the regulation of the canonical Wnt signaling pathway. PCR analysis showed that anesthesia/surgery increased the expression of miR-6912-5p, RNA binding motif protein 47 (Rbm47), and Sostdc1 in the hippocampus, while decreased the expression of miR-490-5p, miR-7001-5p, Cdh3, and Sfrp5. Due to a positive expression correlation between genes and their corresponding circRNAs, circ_0001634/miR-490-5p/Rbm47, circ_0001634/miR-490-5p/Sostdc1, circ_0001634/miR-7001-5p /Sostdc1, circ_0001345/miR-7001-5p/Sostdc1, and circ_000149 3/miR-7001-5p/Sostdc1 axis might be existent in the pathogenesis of POCD. Pei et al. (2022) found that CircFAM53B, highly homologous to mmu_circ_0001634, impeded glioma cell apoptosis through sponging miR-532-3p, suggesting a potential role in neurological disorders. Interestingly, circFbxl5 (homologous with circ_0001345) could function as a miR-146a sponge in mouse cardiomyocyte (Li et al., 2022b), and miR-146a has been proven to ameliorate surgery-induced cognitive decline (Chen et al., 2019a). Accordingly, it is speculated that depletion of circ_0001345 may be an emerging therapeutic perspective on POCD. Most notably, miR-490-5p has been shown to ameliorate stroke-induced neurological dysfunction by repressing cyclin-dependent kinases 1 (CDK1) (Skovira et al., 2016; Ding et al., 2021), revealing a cognition-protective role for miR-490-5p. Rbm47, a central mediator of mRNA alternative splicing and stability, has been elaborated to suppress Wnt activity in cancer cells via maintaining AXIN1 or DKK1 mRNA stability (Vanharanta et al., 2014; Shen et al., 2020). Moreover, Sostdc1 has been documented to negatively modulate Wnt signaling (Li et al., 2022a). These imply that inhibition of Rbm47 or Sostdc1 may mitigate postoperative neurocognitive disorder.

Nonetheless, some limitations of this study must be considered. First, these three datasets are retrieved from different platform and samples, which certainly reduces confidence in the regulator-target pairs outlined. Second, the ideal scenario would be the same logFC threshold and adjusted p-value. However, adjusted p-value is either absent (GSE190880) or greater than 0.9 (GSE95070 and GSE115440), which may be attributed to the relative small sample size and unknown distribution of their genomic data. There are many DEcircRNAs but few DEmiRNAs in the case of p < 0.05 and |logFC| o > 1.0. Hence, p-value and different logFC thresholds were employed for screening DERNAs, as previously described (Wang et al., 2022). Third, limited RNA ID information of other datasets available such as GSE165798 (Wu et al., 2021) hampers the integration of multiple datasets. Finally, the circRNA-miRNA and miRNA-mRNA interaction relationships in the ceRNA network were both based upon a prediction algorithm that required further experimental verification.

In summary, a ceRNA regulatory network in the hippocampus of POCD mice included three circRNAs, three miRNAs, and six mRNAs. Through comprehensive bioinformatic analysis, the present study broadens our horizon on the occurrence and progression of POCD.

References

• 1

  AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene ontology: tool for the unification of biology, the gene ontology consortium.Nat. Genet.2525–29. 10.1038/75556

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BarrettT.WilhiteS. E.LedouxP.EvangelistaC.KimI. F.TomashevskyM.et al (2013). NCBI GEO: archive for functional genomics data sets–update.Nucleic Acids Res.41D991–D995. 10.1093/nar/gks1193

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BooneM. D.SitesB.Von RecklinghausenF. M.MuellerA.TaenzerA. H.ShaefiS. (2020). Economic burden of postoperative neurocognitive disorders among US medicare patients.JAMA Network Open3:e208931. 10.1001/jamanetworkopen.2020.8931

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BriševacD.ScholzR.DuD.ElagabaniM. N.KöhrG.KornauH. C. (2021). The small GTPase Arf6 is dysregulated in a mouse model for fragile X syndrome.J. Neurochem.157666–683. 10.1111/jnc.15230

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CaoC.DengF.HuY. (2020). Dexmedetomidine alleviates postoperative cognitive dysfunction through circular RNA in aged rats.3 Biotech10:176. 10.1007/s13205-020-2163-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChenL.DongR.LuY.ZhouY.LiK.ZhangZ.et al (2019a). MicroRNA-146a protects against cognitive decline induced by surgical trauma by suppressing hippocampal neuroinflammation in mice.Brain Behav. Immun.78188–201. 10.1016/j.bbi.2019.01.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  ChenY.SunJ. X.ChenW. K.WuG. C.WangY. Q.ZhuK. Y.et al (2019b). miR-124/VAMP3 is a novel therapeutic target for mitigation of surgical trauma-induced microglial activation.Signal Transduct. Target Ther.4:27. 10.1038/s41392-019-0061-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  ChenL. L. (2020). The expanding regulatory mechanisms and cellular functions of circular RNAs.Nat. Rev. Mol. Cell Biol.21475–490. 10.1038/s41580-020-0243-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  ChenY.WangX. (2020). miRDB: an online database for prediction of functional microRNA targets.Nucleic Acids Res.48D127–D131. 10.1093/nar/gkz757

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ChenY.ZhangP.LinX.ZhangH.MiaoJ.ZhouY.et al (2020). Mitophagy impairment is involved in sevoflurane-induced cognitive dysfunction in aged rats.Aging1217235–17256. 10.18632/aging.103673

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DanielsonM.WiklundA.GranathF.BlennowK.MkrtchianS.NellgårdB.et al (2020). Neuroinflammatory markers associate with cognitive decline after major surgery: findings of an explorative study.Ann. Neurol.87370–382. 10.1002/ana.25678

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DingL.NingJ.WangQ.LuB.KeH. (2021). Sevoflurane improves nerve regeneration and repair of neurological deficit in brain damage rats via microRNA-490-5p/CDK1 axis.Life Sci.271:119111. 10.1016/j.lfs.2021.119111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DubeU.Del-AguilaJ. L.LiZ.BuddeJ. P.JiangS.HsuS.et al (2019). An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations.Nat. Neurosci.221903–1912. 10.1038/s41593-019-0501-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  GaoR.ChenC.ZhaoQ.LiM.WangQ.ZhouL.et al (2020). Identification of the Potential Key Circular RNAs in Elderly Patients With Postoperative Cognitive Dysfunction.Front. Aging Neurosci.12:165. 10.3389/fnagi.2020.00165

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  GerezJ. A.PrymaczokN. C.RockensteinE.HerrmannU. S.SchwarzP.AdameA.et al (2019). A cullin-RING ubiquitin ligase targets exogenous α-synuclein and inhibits Lewy body-like pathology.Sci. Transl. Med.11:eaau6722. 10.1126/scitranslmed.aau6722

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  GlažarP.PapavasileiouP.RajewskyN. (2014). circBase: a database for circular RNAs.RNA201666–1670. 10.1261/rna.043687.113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HuN.WangC.ZhengY.AoJ.ZhangC.XieK.et al (2016). The role of the Wnt/β-catenin-Annexin A1 pathway in the process of sevoflurane-induced cognitive dysfunction.J. Neurochem.137240–252. 10.1111/jnc.13569

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  Huang daW.ShermanB. T.LempickiR. A. (2009a). Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists.Nucleic Acids Res.371–13. 10.1093/nar/gkn923

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  Huang daW.ShermanB. T.LempickiR. A. (2009b). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat. Protoc.444–57. 10.1038/nprot.2008.211

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  KanehisaM.SatoY.MorishimaK. (2016). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences.J. Mol. Biol.428726–731. 10.1016/j.jmb.2015.11.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  KizilC.KüchlerB.YanJ. J.ÖzhanG.MoroE.ArgentonF.et al (2014). Simplet/Fam53b is required for Wnt signal transduction by regulating β-catenin nuclear localization.Development1413529–3539. 10.1242/dev.108415

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  KristensenL. S.AndersenM. S.StagstedL. V. W.EbbesenK. K.HansenT. B.KjemsJ. (2019). The biogenesis, biology and characterization of circular RNAs.Nat. Rev. Genet.20675–691. 10.1038/s41576-019-0158-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  LiD.YouJ.MaoC.ZhouE.HanZ.ZhangJ.et al (2022b). Circular RNA Fbxl5 Regulates Cardiomyocyte Apoptosis During Ischemia Reperfusion Injury via Sponging microRNA-146a.J. Inflamm. Res.152539–2550. 10.2147/jir.S360129

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  LiC.WangM.ShiY.XinH. (2022a). SOSTDC1 acts as a tumor inhibitor in acute myeloid leukemia by downregulating the Wnt/β-catenin pathway.Environ. Toxicol.371934–1943. 10.1002/tox.23540

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  LiM.HuJ.PengY.LiJ.RenR. (2021). CircPTK2-miR-181c-5p-HMGB1: a new regulatory pathway for microglia activation and hippocampal neuronal apoptosis induced by sepsis.Mol. Med.27:45. 10.1186/s10020-021-00305-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  LiangL.XieR.LuR.MaR.WangX.WangF.et al (2020). Involvement of homodomain interacting protein kinase 2-c-Jun N-terminal kinase/c-Jun cascade in the long-term synaptic toxicity and cognition impairment induced by neonatal Sevoflurane exposure.J. Neurochem.154372–388. 10.1111/jnc.14910

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  McDonaghD. L.MathewJ. P.WhiteW. D.Phillips-ButeB.LaskowitzD. T.PodgoreanuM. V.et al (2010). Cognitive function after major noncardiac surgery, apolipoprotein E4 genotype, and biomarkers of brain injury.Anesthesiology112852–859. 10.1097/ALN.0b013e3181d31fd7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  McGearyS. E.LinK. S.ShiC. Y.PhamT. M.BisariaN.KelleyG. M.et al (2019). The biochemical basis of microRNA targeting efficacy.Science366:eaav1741. 10.1126/science.aav1741

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  NarvaesR. F.FuriniC. R. G. (2022). Role of Wnt signaling in synaptic plasticity and memory.Neurobiol. Learn. Mem.187:107558. 10.1016/j.nlm.2021.107558

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  NiuY.WanC.ZhangJ.ZhangS.ZhaoZ.ZhuL.et al (2021). Aerobic exercise improves VCI through circRIMS2/miR-186/BDNF-mediated neuronal apoptosis.Mol. Med.27:4. 10.1186/s10020-020-00258-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  OtasekD.MorrisJ. H.BouçasJ.PicoA. R.DemchakB. (2019). Cytoscape automation: empowering workflow-based network analysis.Genome Biol.20:185. 10.1186/s13059-019-1758-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  PeiJ.DouH.DengX. (2022). CircFAM53B promotes the proliferation and metastasis of glioma through activating the c-MET/PI3K/AKT pathway via sponging miR-532-3p.Cell Cycle21462–476. 10.1080/15384101.2021.2014738

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  RanW.LiangN.YuanR.WangZ.GaoJ. (2022). Identification of potential key circRNAs in aged mice with postoperative delirium.Front. Mol. Neurosci.15:836534. 10.3389/fnmol.2022.836534

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  Rybak-WolfA.StottmeisterC.GlažarP.JensM.PinoN.GiustiS.et al (2015). Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed.Mol. Cell58870–885. 10.1016/j.molcel.2015.03.027

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  ShenD. J.JiangY. H.LiJ. Q.XuL. W.TaoK. Y. (2020). The RNA-binding protein RBM47 inhibits non-small cell lung carcinoma metastasis through modulation of AXIN1 mRNA stability and Wnt/β-catentin signaling.Surg. Oncol.3431–39. 10.1016/j.suronc.2020.02.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  SinghT.YadavS. (2020). Role of microRNAs in neurodegeneration induced by environmental neurotoxicants and aging.Ageing Res. Rev.60:101068. 10.1016/j.arr.2020.101068

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  SkoviraJ. W.WuJ.MatyasJ. J.KumarA.HanscomM.KabadiS. V.et al (2016). Cell cycle inhibition reduces inflammatory responses, neuronal loss, and cognitive deficits induced by hypobaria exposure following traumatic brain injury.J. Neuroinflammation13:299. 10.1186/s12974-016-0769-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  SongC.ZhangY.HuangW.ShiJ.HuangQ.JiangM.et al (2022). Circular RNA Cwc27 contributes to Alzheimer’s disease pathogenesis by repressing Pur-α activity.Cell Death Differ.29393–406. 10.1038/s41418-021-00865-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  SuwanabolP. A.LiY.AbrahamseP.De RooA. C.VuJ. V.SilveiraM. J.et al (2022). Functional and cognitive decline among older adults after high-risk surgery.Ann. Surg.275e132–e139. 10.1097/sla.0000000000003950

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  TangV. L.JingB.BoscardinJ.NgoS.SilvestriniM.FinlaysonE.et al (2020). Association of functional, cognitive, and psychological measures with 1-year mortality in patients undergoing major surgery.JAMA Surg.155412–418. 10.1001/jamasurg.2020.0091

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  VanharantaS.MarneyC. B.ShuW.ValienteM.ZouY.MeleA.et al (2014). Loss of the multifunctional RNA-binding protein RBM47 as a source of selectable metastatic traits in breast cancer.eLife3:e02734. 10.7554/eLife.02734

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  WangJ.NiuY.LuoL.LuZ.ChenQ.ZhangS.et al (2022). Decoding ceRNA regulatory network in the pulmonary artery of hypoxia-induced pulmonary hypertension (HPH) rat model.Cell Biosci.12:27. 10.1186/s13578-022-00762-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  WangM.SuP.LiuY.ZhangX.YanJ.AnX.et al (2019). Abnormal expression of circRNA_089763 in the plasma exosomes of patients with post-operative cognitive dysfunction after coronary artery bypass grafting.Mol. Med. Rep.202549–2562. 10.3892/mmr.2019.10521

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  WangY. L.ZhangY.CaiD. S. (2021). Dexmedetomidine ameliorates postoperative cognitive dysfunction via the MicroRNA-381-Mediated EGR1/p53 Axis.Mol. Neurobiol.585052–5066. 10.1007/s12035-021-02417-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  WeiC.LuoT.ZouS.ZhouX.ShenW.JiX.et al (2017). Differentially expressed lncRNAs and miRNAs with associated ceRNA networks in aged mice with postoperative cognitive dysfunction.Oncotarget855901–55914. 10.18632/oncotarget.18362

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  WuJ.YangJ. J.CaoY.LiH.ZhaoH.YangS.et al (2020). Iron overload contributes to general anaesthesia-induced neurotoxicity and cognitive deficits.J. Neuroinflammation17:110. 10.1186/s12974-020-01777-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  WuY. Q.LiuQ.WangH. B.ChenC.HuangH.SunY. M.et al (2021). Microarray analysis identifies key differentially expressed circular RNAs in aged mice with postoperative cognitive dysfunction.Front. Aging Neurosci.13:716383. 10.3389/fnagi.2021.716383

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  YangT.VelagapudiR.TerrandoN. (2020). Neuroinflammation after surgery: from mechanisms to therapeutic targets.Nat. Immunol.211319–1326. 10.1038/s41590-020-00812-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  YangT.XuG.NewtonP. T.ChaginA. S.MkrtchianS.CarlströmM.et al (2019). Maresin 1 attenuates neuroinflammation in a mouse model of perioperative neurocognitive disorders.Br. J. Anaesth.122350–360. 10.1016/j.bja.2018.10.062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  YoonG.LimY. H.JoD.RyuJ.SongJ.KimY. K. (2021). Obesity-linked circular RNA circTshz2-2 regulates the neuronal cell cycle and spatial memory in the brain.Mol. Psychiatry266350–6364. 10.1038/s41380-021-01303-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  ZhangL. Y.ChenL. S.SunR.JiS. J.DingY. Y.WuJ.et al (2013). Effects of expression level of DNA repair-related genes involved in the NHEJ pathway on radiation-induced cognitive impairment.J. Radiat. Res.54235–242. 10.1093/jrr/rrs095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  ZhangM. X.LinJ. R.YangS. T.ZouJ.XueY.FengC. Z.et al (2022). Characterization of circRNA-Associated-ceRNA networks involved in the pathogenesis of postoperative cognitive dysfunction in aging mice.Front. Aging Neurosci.14:727805. 10.3389/fnagi.2022.727805

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  ZhouH.LiF.YeW.WangM.ZhouX.FengJ.et al (2020). Correlation between plasma CircRNA-089763 and postoperative cognitive dysfunction in elderly patients undergoing non-cardiac surgery.Front. Behav. Neurosci.14:587715. 10.3389/fnbeh.2020.587715

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article