Wensi Wu1†Yongpai Peng2†Jiaxin Zhou1Xiaojun Zhang1Lin Cao1
Wei-jye Lin3,4
Yanan Lu1Jing Wen1
Zhi Wang1*- 1Department of Anesthesiology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 2Department of Gynecological Oncology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 3Guangdong Provincial Key Laboratory of Malignant Tumor Epigenetics and Gene Regulation, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 4Medical Research Center, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
Postoperative neurocognitive disorder (PND) is one of the most common postoperative neurological complications in aged patients, characterized by mental disorder, anxiety, personality changes, and impaired memory. At present, the molecular mechanism of PND remains largely unclear, and the ideal biomarker for clinical diagnosis and prognosis are lacking. Circular RNA (circRNA) and microRNA (miRNA), as unique non-coding RNAs, affecting the regulation of miRNAs on genes and further intervening in the progression of diseases through the sponge action between the two. Besides, it could be served as novel biomarkers in various diseases. In order to detect the differential expression profiles of genes caused by PND, a total of 26 18-month-old male C57BL/6 mice were randomly assigned to control group and PND group. Behavioral tests showed that mice in the PND group had impaired cognitive function compared with the control group. Three mice in each group were randomly selected to harvest the brain for analysis the expressions of circRNAs, miRNAs, and mRNAs in the prefrontal cortex by next-generation sequencing (NGS) technology. Differentially expressed genes, including 1192 circRNAs, 27 miRNAs, and 266 mRNAs were identified, and its accuracy was further confirmed by qRT-PCR. Bioinformatics analysis results suggested that neuroinflammation was the main pathological mechanism of PND. The construction of competitive endogenous RNA (ceRNA) networks and the identification of hub genes provided possible therapeutic targets for PND. Cinnarizine and Clemastine were predicted to have the potential therapeutic effects on PND. This is the first study to explore the differential expression profiles of genes and their regulation mechanisms in PND, our results provided new clues and targets for the treatment of this refractory disease.
Introduction
Postoperative neurocognitive disorder (PND) is mainly manifested as reduced learning and memory ability, impaired social ability, and cognitive ability after non-brain operation (Evered et al., 2018). PND is a cognitive disorder disease that has attracted much attention in clinical work. Advanced age is considered as an independent risk factor of PND (Luo et al., 2019). Among aged patients over 60-year old, the overall incidence of PND is 10–15%, which is significantly higher than other postoperative complications such as pulmonary infection, stroke, and heart failure (Urits et al., 2019). PND increases the risk of Alzheimer’s disease and might lead to dementia (Evered et al., 2016), which seriously affects the quality of life of patients and imposes a heavy burden on society and patients’ families. Many studies have been performed to explore possible molecular mechanisms of PND. Neuroinflammation, neuronal apoptosis, amyloid-Aβ deposition, and hyperphosphorylation of tau proteins have been reported as main mediators in PND (Xie et al., 2013; Chen et al., 2019; Li M. et al., 2019). At present, the pathogenesis of PND remains elusive, and there is no effective intervention method.
The next-generation sequencing (NGS) technology, also known as high-throughput sequencing technology, can perform large-scale genomics researches. This application can reveal the potential pathological mechanism and provide new ideas for the treatment of disease (Jakel and Williams, 2020). Circular RNA (circRNA) and microRNA (miRNA) belong to non-coding RNAs. In terms of function, recent studies have shown that circRNA molecules are rich in miRNA binding sites and served as miRNA sponges, thereby releasing the inhibitory effect of miRNA on its target genes and increasing the expression level of target genes. This is the competitive endogenous RNA (ceRNA) mechanism. The ceRNA networks are widely involved in the process of various cognitive disorder diseases such as Alzheimer’s disease (Zhang et al., 2017; Wang et al., 2018b; Ma et al., 2019). Previous studies showed that different expressions of circRNAs were detected in the hippocampus of rats with PND (Cao et al., 2020). The differentially expressed circRNAs also had been determined in the serum of aged PND patients. Therefore, circRNAs may be related to the pathogenesis of PND (Gao et al., 2020). Meanwhile, the up-regulated expression of miRNA-572 in hippocampal neurons contributes to the recovery of cognitive function in rats with PND (Yu et al., 2015), and the miRNA-190a also plays a key role in the pathogenesis of PND (Liu et al., 2019). However, the differential expression profiles and possible ceRNA networks caused by PND remains unclear.
Therefore, in this study, the differential expression profiles in prefrontal cortex of aged mice with PND were analyzed by NGS. The ceRNA networks and hub genes were identified, and relevant functional enrichment analyses were performed by bioinformatics. In addition, some small molecular drugs with potential therapeutic effects on PND were predicted. These results might help provide new biomarkers and therapeutic targets for PND.
Materials and Methods
Animals
All the animal experiments in this study were approved by the Institutional Animal Care and Use Committee (Approval No.: SYSU-IACUC-2020-000326) and the Laboratory Animal Ethics Committee of Sun Yat-sen University.
The Sun Yat-sen University (Guangzhou, China) provided 18-month-old male C57BL/6 mice with a weight between 45 and 50 g for the experiment. The mice were randomly divided into two groups: isoflurane plus exploratory laparotomy (ISO + EL) group and control (CON) group. Each group had 13 mice. In a colony room with a temperature of 19–22°C, a humidity of 40–60%, and a light/dark cycle of 12 h (lighting time is 07:00–19:00), 5 mice were housed in each cage with access to food and water at will. The experiments began after all the animals had adapted to the environment for 2 weeks.
Operation and Anesthesia
Isoflurane anesthesia plus exploratory laparotomy has been proved to be an effective method to construct PND model (Qiu et al., 2020). Before exploratory laparotomy, the mice were anesthetized by exposure to an oxygen chamber prefilled with 1.5% isoflurane for 30 min. A median incision approximately 2 cm in the abdomen was made to enter the abdominal cavity and explore the abdominal organs such as the liver, spleen, and intestine. Sterile 5-0 operation sutures were used to suture the peritoneum and skin. Operation was also performed under isoflurane anesthesia and lasted 30 min. At the end of operation and every day within 3 days after operation, 2.5% lidocaine cream was applied to the incision to alleviate the postoperative pain, and povidone iodine solution was applied to prevent infection. Mice in the control group were not treated with anesthesia or operation.
Open Field
A black opaque plastic chamber (60 cm × 60 cm × 50 cm, ZH-ZFT, Anhui Zhenghua Biological Instrument equipment Co., Ltd., Anhui, China) was used as the open field arena. The open field test was performed 24 h after operation to evaluate the locomotor activity of the mice. Each mouse was placed in the center of the field and allowed to explore freely for 5 min while a video tracking system (Smart v3.0.06, Panlab Harvard Apparatus, Barcelona, ES) automatically recorded its movements and analyzed the total distance to assess the locomotor activity. Meanwhile, the time spent in the central area of mice were analyzed to evaluate whether postoperative anxiety occurred. During each test interval, the field was cleaned with 75% ethanol to eliminate feces and odors.
Fear Conditioning
One day before the operation, each mouse was placed into the conditioning chamber (Freeze Monitor, San Diego Instruments, San Diego, CA, United States) and allowed to explore the room freely for 180 s. Then they were given a 30 s tone (70 db), followed by a 2 s foot shock (0.7 mA), and the next tone-shock stimulation cycle was entered at an interval of 60 s. A total of 3 cycles were performed. On postoperative days 1, 2 h after the open field test, each mouse was placed into the conditioning chamber without any tone or electrical stimulation for 360 s and the environment was identical with that before operation, the time of freezing behavior was recorded to evaluate the context-related memory. Two-hours later, they were placed into a new environment completely different from that before the operation and explored room freely for 180 s. Then they were given the same tone stimulation (without electrical stimulation) as before the operation, the time of freezing behavior was recorded to evaluate tone-related memory. Freezing behavior means there is no visible movement other than breathing. During each test interval, the conditioning chamber was cleaned with 75% ethanol to eliminate feces and odors.
Barnes Maze
On postoperative days 2–6, we performed Barnes maze to evaluate the spatial learning and memory of mice. Each mouse was placed in the center of a circular platform (Anhui Zhenghua Biological Instrument equipment Co., Ltd., Anhui, China) with a diameter of 92 centimeters, which had 20 equally spaced holes. Among all the holes, only one was linked to a small dark recessed chamber. The mice were expected to find the hole and enter the dark chamber under the bright light (200 W). The mice were trained for 4 days with 3 trials per day, with each trial lasting 3 min, and an interval of 15 min between each trial. If the mice could not find the correct hole and enter the dark chamber within 3 min, they were guided to the correct location. During each test interval, the platform and dark chamber were cleaned with 75% ethanol to eliminate feces and odors. The last day, the escape latency was recorded and measured by a video tracking system (Smart v3.0.06, Panlab Harvard Apparatus, Barcelona, ES), which was used to evaluate the spatial learning and memory retention of mice.
Harvesting of Tissue
Half an hour after the Barnes maze test, the mice were deeply anesthetized with isoflurane, and perfused transcardially with normal saline. We dissected the brain in the air to obtain the prefrontal cortex, and then frozen tissues in liquid nitrogen. Three cortical tissues were randomly selected from each group for subsequent sequencing analysis. The others were preserved for qRT-PCR.
RNA Extraction and Purification for Next-Generation Sequencing
Total RNAs were extracted by miRNeasy Mini Kit (217004, Qiagen, GmBH, Germany) and checked for a RIN number to inspect RNA integrity by an Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, United States). Qualified total RNAs were further purified by RNAClean XP Kit (A63987, Beckman Coulter, Inc., Kraemer Boulevard Brea, CA, United States) and RNase-Free DNase Set (79254, Qiagen, GmBH, Germany). The purity of total RNAs were assessed by NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States) and Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, United States). The integrity of RNAs were determined by standard denaturing agarose gel electrophoresis.
Library Preparation and Next-Generation Sequencing of circRNAs
Ribosomal RNAs (rRNAs) were removed from total RNAs, and then RNase R (RNR07250, Epicenter, United States) was used to digest linear RNAs. Subsequently, purified RNAs fragments were subjected to the first strand and second strand cDNA synthesis, and then adapter-ligated and enriched were performed by VAHTS Total RNA-seq (H/M/R) Library Prep Kit for Illumina (NR603-02, Vazyme, CHN). The quality and concentration of library were measured by the Qubit dsDNA HS Assay Kit (Q32854, Invitrogen, Carlsbad, CA, United States). Illumina HiSeq X ten platform was used for circRNAs sequencing to obtain the raw reads. Raw reads were filtered with Seqtk1 and obtained the clean reads. The results of clean reads were mapped to the reference genome by BWA-MEM2. CircRNAs were detected by CIRI (Beijing Institutes of Life Science, Beijing, China) (Gao et al., 2015), based on the circRNAs location information, the same circRNAs in each sample were merged, renumbered, and compared with the circBase database3. At present, the complete sequence cannot be obtained for most circRNAs (You et al., 2015), so we used SRPBM to normalize the clean reads (Li Y. et al., 2015). We used edgeR package of R software (version 3.6.04) to analyze the different circRNAs between samples and get the p-value (Robinson et al., 2010). At the same time, we calculated the differentially expressed multiples based on the SRPBM value, which was Fold-change. The differential expressions of circRNAs were identified using the Limma package of R Software. The differentially expressed circRNAs were screened with | Fold-change| > 2 and p-value < 0.05 as the cut-off criteria. Volcano map and heatmap with differential expression of circRNAs were presented by ggpubr, ggthemes, and pheatmap packages of R software.
Library Preparation and Next-Generation Sequencing of miRNAs
The adapter-ligated and enriched of purified total RNAs were performed by the NEBNext Multiplex Small RNA Library Prep Set for Illumina (Set 2) (E7580S, NEB, United States). After the amplification, the quality, and concentration of library were measured by the Qubit dsDNA HS Assay Kit (Q32854, Invitrogen, Carlsbad, CA, United States). Illumina HiSeq X ten platform was used for miRNAs sequencing to obtain the raw reads. Raw reads were filtered with fastx5 and obtained the clean reads. Clean reads of 18–40 nt length were compared with the reference genome by Bowtie (Langmead et al., 2009). And miRNAs were identified based on the location information of known miRNAs in miRBase (Kozomara and Griffiths-Jones, 2014). In order to make the expression levels of miRNAs between different miRNAs and different samples comparable, we normalized the number of clean reads mapped to each miRNA by the trimmed mean of M values (TMM) (Robinson and Oshlack, 2010), and then converted it into the transcripts per million (TPM, the formula is: the number of reads on a miRNA × 106/total reads) for standardization of miRNAs expression (Mortazavi et al., 2008). We used edgeR package of R software to analyze the different miRNAs between samples and get the p-value. At the same time, we calculated the differentially expressed multiples based on the TPM value, which is Fold-change. The differential expressions of miRNAs were identified using the Limma package of R Software. The differentially expressed miRNAs were screened with | Fold-change| > 2 and p-value < 0.05 as the cut-off criteria. Volcano map and heatmap with differential expression of miRNAs were presented by ggpubr, ggthemes, and pheatmap packages of R software.
Library Preparation and Next-Generation Sequencing of mRNAs
Similar to the circRNAs sequencing process, rRNAs were removed from total RNAs, purified RNAs fragments were subjected to the first strand and second strand cDNA synthesis, and then adapter-ligated and enriched were performed by VAHTS Total RNA-seq (H/M/R) Library Prep Kit for Illumina (NR603-02, Vazyme, CHN). The quality and concentration of library were measured by the Qubit dsDNA HS Assay Kit (Q32854, Invitrogen, Carlsbad, CA, United States). Illumina HiSeq X ten platform was used for mRNAs sequencing to obtain the raw reads. Raw reads were filtered with Seqtk and clean reads was obtained. Hisat2 spliced mapping algorithm was applied to perform genome mapping of pre-processed reads (Kim et al., 2015). In order to make data comparable, processed reads were transformed into fragments per kilobase of exon model per million mapped reads (FPKM, the formula is: total exon fragments/mapped reads/exon length) for the standardization of mRNAs expression. We used StringTie to count the number of fragments of each gene after Hisat2 comparison (Pertea et al., 2015, 2016), then used TMM value to make normalize, and finally used perl (version 5.32.16) to calculate the number of fragments for each gene FPKM value. EdgeR package of R software was also used to analyze the different mRNAs between samples, and the differentially expressed multiples were calculated based on the FPKM value. The differential expressions of mRNAs were identified using the Limma package of R Software. The differentially expressed mRNAs were screened with | Fold-change| > 2 and p-value < 0.05 as the cut-off criteria. Volcano map and heatmap with differential expression of mRNAs were presented by ggpubr, ggthemes and pheatmap packages of R software.
Functional Enrichment Analysis
Gene ontology (GO) enrichment analysis is a commonly used method of gene annotation, including biological pathways (BP), molecular functions (MF), and cell components (CC). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was used to identify the key biological pathways involved in candidate genes. The clusterProfiler package of R software was used to perform GO and KEGG pathway enrichment analysis was based on the differentially expressed mRNAs with adjusted p-value < 0.05 as the cut-off criterion. The p-value was adjusted by false discovery rate (FDR) method.
Identification of Protein-Protein Interaction Networks and Hub Genes
The Search Tool for the Retrieval of Interacting Genes (STRING, version 11.07) database was used to predict protein-protein interactions (PPI) between candidate genes. In order to explore the potential correlation between these differentially expressed mRNAs, we used STRING to obtain the data of PPI. The protein-protein minimum required interaction score was set to 0.7, indicating high confidence. Furthermore, hub genes were identified based on the count of interaction pairs between proteins. Cytoscape (version 3.8.0, San Diego, CA, United States) was used to visualize the PPI network data and identify the hub genes.
Construction of ceRNA Networks
The multiMiR package of R software collects nearly 50 million human and mouse miRNA-mRNA interaction records from 14 different databases (Ru et al., 2014). We used multiMiR to retrieve miRNA-mRNA interactions. Based on 108 CLIP-Seq (PAR-CLIP, HITS-CLIP, iCLIP, and CLASH) data sets, starBase database (v2.08) identified approximately 10,000 ceRNA pairs from 37 independent studies. The starBase database was applied to retrieve miRNA-circRNA interactions (Li J. H. et al., 2014). The obtained genes were intersected with the sequencing results via an online Venn tool9. Based on the inter-relationships of circRNAs-miRNAs-mRNAs, ceRNA networks were then constructed by Cytoscape to visualize the genes regulation networks.
cDNA Synthesis and Quantitative Real-Time PCR (qRT-PCR)
To verify the accuracy and reliability of the sequencing results, 5 circRNAs (circRalgapa1, cird3, circPpme1, circKansl1l, and circRnf150), 2 miRNAs (miR-200b-3p and miR-141-3p), and 5 mRNAs (Slco1b2, Ino80d, Chrnb1, Sncg, and Agt) were randomly selected for qRT-PCR, which were selected from the ceRNA network, the most obvious differential expression genes and the hub genes. Total RNA was extracted from the prefrontal cortex of the PND group and CON group (n = 10, each group) using RNA Quick Purification kit (RN001, Esscience, CHN). For circRNAs and mRNAs, the reverse transcription to cDNA was achieved with Hifair® III 1st Strand cDNA Synthesis SuperMix for qPCR (11141ES60, Yeasen, CHN). For miRNAs, the reverse transcription was achieved with Bulge-Loop miRNA qRT-PCR Starter Kit (C10211-1, RiboBio, and CHN). And qRT-PCR was conducted with Roche LightCycler 480 II Real-Time PCR System (Roche, United States). The relative expression levels were calculated by the 2–ΔΔCT method. GAPDH was used as an internal reference gene for circRNAs and mRNAs expression, and U6 was used as an internal reference gene for miRNAs expression. Sequence-specific primers for U6 and miRNAs were synthesized by RiboBio (Guangzhou, China). The sequences of primers of mRNAs, circRNAs, and GAPDH were presented in Table 1.
Table 1. The primer sequences used for qRT-PCR.
Prediction of Potential Therapeutic Small Molecular Drugs
Connectivity Map (CMap10) was used to measure connectivity between disease gene expression characteristics and compound induced gene expression profiles, which has become a valuable tool in understanding drug mechanism of action and discovering new indications for drugs (Cheng et al., 2014). These mRNAs significantly up-regulated or down-regulated were uploaded onto the CMap database to predict the small molecular drugs with potential therapeutic effects on PND.
Statistical Analysis
All results were expressed as mean ± SD. Statistical analyses were performed using the GraphPad Prism 8.0 (San Diego, CA, United States) and R software. Student’s t-test was used to determine the statistical significance between two groups. A p-value < 0.05 was considered statistically significant.
Results
Isoflurane Plus Exploratory Laparotomy Caused PND in Aged Mice
Our previous study has confirmed that isoflurane plus exploratory laparotomy (ISO + EL) successfully induced PND model in aged mice (Wang et al., 2018a). In the present study, there was no statistical difference in the total distance between the control (CON) group and the ISO + EL group at the same time window in the open field test. The result indicating that anesthesia and operation factors have no effect on the locomotor activity of the aged mice (Figure 1A). At the same time, the time spent in the central area had no statistically different, indicating that anesthesia and operation factors did not induce postoperative anxiety behaviors of aged mice, which was consistent with the results of previous studies (Qiu et al., 2020; Figure 1B). Fear conditioning test was divided into two parts: the context-related test and the tone-related test. In the context test, we found that aged mice in the ISO + EL group performed less freezing behavior in the new context than mice in the CON group, which indicated anesthesia and operation factors impaired the context-related memory of aged mice (Figure 1C). However, in the tone test, we found that the performance of mice was similar between in the ISO + EL group and in the CON group (Figure 1D). During the training phase of the Barnes maze, compared with day one, the time required for all mice to find the target hole significantly decreased on day four (Figure 1E). This suggested that all the mice got performance development during the training process. However, we found that mice in the ISO + EL group took longer time to identify the target hole than mice in the CON group during test phase (Figure 1F). This suggested that anesthesia and operation factors impaired the spatial learning and memory in older mice. In general, ISO + EL impaired the spatial learning and memory and the context-related memory in aged mice, which met the criteria of PND, our results were similar to previous studies (Lin and Zuo, 2011; Li J. et al., 2019; Lin et al., 2020). Therefore, the ISO + EL group was PND group.
Figure 1. Anesthesia and operation-induced PND in aged mice. (A) Total distance of aged mice in open field test. (B) Time spent in the central area of aged mice in open field test. (C) Performance of aged mice in context test. (D) Performance of aged mice in tone test. (E) Performance of aged mice in Barnes maze training phase. (F) Performance of aged mice in Barnes maze test phase. The data are presented as mean ± S.D (n = 10). **P < 0.01 compared with the CON group; ****P < 0.001 compared with the CON group, #P < 0.05 compared with the corresponding data in the first trial on day 1.
Identification of Differentially Expressed circRNAs
In this study, NGS technology was used to identify the differentially expressed circRNAs in the prefrontal cortex of aged mice with PND. A total of 71481 circRNAs were identified. Compared with the CON group, we had screened 1192 differentially expressed circRNAs in the PND group, which included 1025 up-regulated circRNAs and 167 down-regulated circRNAs. The hierarchical clustering heatmaps revealed distinguishable circRNAs expression profiles (Figure 2A), and the volcano plots showed a significant difference in the expression levels of circRNAs between the PND group and CON group (Figure 2B). In the prefrontal cortex of aged mice with PND, most of the differentially expressed circRNAs were originated from exons (Figure 2C).
Figure 2. Altered expression profile of circRNAs in the prefrontal cortex of aged mice. (A) Hierarchical clustering analysis of all differentially expressed circRNAs in the prefrontal cortex of CON group and PND group: the red represents a higher fold change, while the blue represents a lower fold change. (B) The volcano plots of all the detected circRNAs in both the CON and PND groups: The vertical line represents to | log2 (fold change) | > 1 up and down, respectively, and the horizontal line corresponds to a p-value of 0.05 (–log10 scaled). The red points represent the significantly up-regulated circRNAs and the blue points represent the significantly down-regulated circRNAs. The top ten differentially expressed circRNAs were also indicated in the figure. (C) The gene sources of differentially levels of circRNAs. The red represents the significantly up-regulated circRNAs and the blue represents the significantly down-regulated circRNAs.
Identification of Differentially Expressed miRNAs
Compared with the CON group, a total of 1189 miRNAs were identified in the prefrontal cortex of aged mice with PND. We had screened 27 differentially expressed miRNAs in the PND group, which included 12 up-regulated miRNAs and 15 down-regulated miRNAs. The hierarchical clustering heatmaps revealed distinguishable miRNAs expression profiles (Figure 3A), and the volcano plots showed a significant difference in the expression levels of miRNAs between the PND group and CON group (Figure 3B).
Figure 3. Altered expression profile of miRNAs in the prefrontal cortex of aged mice. (A) Hierarchical clustering analysis of all differentially expressed miRNAs in the prefrontal cortex of CON group and PND group: the red represents a higher fold change, while the blue represents a lower fold change. (B) The volcano plots of all the detected miRNAs in both the CON and PND groups: The vertical line represents to | log2 (fold change) | > 1 up and down, respectively, and the horizontal line corresponds to a p-value of 0.05 (–log10 scaled). The red points represent the significantly up-regulated miRNAs and the blue points represent the significantly down-regulated miRNAs. The top ten differentially expressed miRNAs were also indicated in this figure.
Identification of Differentially Expressed mRNAs
Compared with the CON group, a total of 21006 mRNAs were identified in the prefrontal cortex of aged mice with PND. We had screened 266 differentially expressed mRNAs in the PND group, which included 160 up-regulated mRNAs and 106 down-regulated mRNAs. The hierarchical clustering heatmaps revealed distinguishable mRNAs expression profiles (Figure 4A), and the volcano plots showed mRNAs levels with significant difference between the PND group and CON group (Figure 4B).
Figure 4. Altered expression profile of mRNAs in the prefrontal cortex of aged mice. (A) Hierarchical clustering analysis of all differentially expressed mRNAs in the prefrontal cortex of CON group and PND group: the red represents a higher fold change, while the blue represents a lower fold change. (B) The volcano plots of all the detected mRNAs in both the CON and PND groups: The vertical line represents to | log2 (fold change) | > 1 up and down, respectively, and the horizontal line corresponds to a p-value of 0.05 (–log10 scaled). The red points represent the significantly up-regulated mRNAs and the blue points represent the significantly down-regulated mRNAs. The top ten differentially expressed mRNAs were also indicated in this figure.
Bioinformatics Analysis of Differentially Expressed mRNAs
Although differentially expressed circRNAs and miRNAs have been identified in prefrontal cortex, proteins transcribed from mRNAs actually play key biological functions in vivo and in vitro. Therefore, we performed GO enrichment analysis and KEGG pathway analysis to predict the possible biological functions of the differentially expressed mRNAs. The top GO terms in the biological processes (BP) were included as below: I-kappaB kinase/NF-kappaB signaling, proteasome-mediated ubiquitin-dependent protein catabolic process, lipid localization, feeding behavior, and lipid transport. The top GO terms in the cellular components (CC) were listed as below: receptor complex, integrin complex, and protein complex involved in cell adhesion, interstitial matrix, and plasma membrane signaling receptor complex. The top GO terms in the molecular function (MF) were include as below: pattern recognition receptor activity, postsynaptic neurotransmitter receptor activity, extracellular ligand-gated ion channel activity, ketosteroid monooxygenase activity, and hormone activity (Figure 5A). Through KEGG pathway analysis, the top 10 enriched pathways found as below: neuroactive ligand-receptor interaction, insulin resistance, pathways of neurodegeneration-multiple diseases, vascular smooth muscle contraction, regulation of actin cytoskeleton, hippo signaling pathway, herpes simplex virus 1 infection, inflammatory mediator regulation of TRP channels, adipocytokine signaling pathway, and proteoglycans in cancer (Figure 5B).
Figure 5. Bioinformatics analysis of differentially expressed mRNAs in the prefrontal cortex of aged mice. (A) Gene ontology (GO) analysis of the differentially expressed mRNAs in the prefrontal cortex of aged mice. (B) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the differentially expressed mRNAs in the prefrontal cortex of aged mice.
Construction of PPI Network
The 266 differentially expressed mRNAs were upload onto the STRING database to explore the inter-relationships among the different proteins, a total of 198 nodes and 36 edges were obtained. We visualized these data by Cytoscape. And a total of 42 mRNAs were mapped into the network, which included 23 up-regulated mRNAs and 19 down-regulated mRNAs (Figure 6A). And then, we identified the top 6 hub genes (Agt, Chrm5, Grp, Hcrtr1, Trh, and Top2a) by the count of interaction pairs between different proteins, indicating that these genes may have important roles in PND (Figure 6B).
Figure 6. Construction of PPI network and identification of hub genes. (A) Visualization of PPI network. The red points represent 23 up-regulated genes and the blue points represent 19 down-regulated genes. (B) Identification of the top 6 hub genes. The red bars represent up-regulated genes and the blue bars represent down-regulated genes.
Construction of ceRNA Networks
Most circRNAs can act as miRNAs sponges to eliminate the inhibition of miRNA on their downstream mRNAs expressions. And this is the main mechanism of ceRNA networks. Based on differentially expressed miRNAs, 917 downstream mRNAs were obtained by the multiMiR package of R software. We intersected these mRNAs with the differentially expressed mRNAs which obtained by NGS. And 5 mRNAs were identified as the downstream mRNAs of miRNAs (Figure 7A). Their corresponding relationships with miRNAs were shown in Table 2. By the starBase database, 40 circRNAs which were obtained could competitively bind with the differentially expressed miRNAs. Meanwhile, 13 circRNAs were identified as the ceRNAs of miRNAs (Figure 7B). The possible corresponding relationships between miRNAs and circRNAs were shown in Table 3. Based on the relationships between miRNAs-mRNAs and miRNAs-circRNAs obtained by databases, combined with the expression trend of genes, a ceRNA network were visualized by Cytoscape (Figure 7C). In this network, we found that the circTcf4, circPlcb1, circMemo1, circCtnnd2, circRalgapa1/miR-200b-3p/Slco1b2 axes, and the circPsd3, circRabgap1, circTulp4, circZbtb20/miR-141-3p/Ino80d axes conform to the general gene expression law of ceRNA network, which might be involved in the mechanisms of PND regulation.
Figure 7. Construction of ceRNA networks. (A) A total of 5 mRNAs were identified with miRNA targeting binding sites. The blue circle represents the results of bioinformatics databases and the red circle represents the NGS results. (B) A total of 13 circRNAs were identified with miRNA targeting binding sites. The blue circle represents the results of bioinformatics databases and the red circle represents the NGS results. (C) Visualization of ceRNA networks. The red marks represent the up-regulated genes and the blue marks represent the down-regulated genes.
Table 2. Correspondence relationships between miRNAs and mRNAs.
Table 3. Correspondence relationships between miRNAs and circRNAs.
Validation of the Differentially Expressed Genes
To validate the accuracy and reliability of the NGS results, we randomly selected 5 circRNAs, 2 miRNAs, and 5 mRNAs for qRT-PCR. Compared with CON group, the expressions of circRalgapa1, circPsd3, circPpme1, circKansl1l, cinrcRnf150, Slco1b2, Ino80d, chrnb1, Sncg, and Agt in PND group were significantly up-regulated (Figure 8). The expressions of miR-200b-3p and miR-141-3p were significantly down-regulated (Figure 9). The results of qRT-PCR largely verified the accuracy of NGS results and provided potential therapeutic targets for the treatment of PND.
Figure 8. The validation results of circRalgapa1, circPsd3, circPpme1, circKansl1l, circRnf150, Slco1b2, Ino80d, Chrnb1, Sncg, and Agt by qRT-PCR in CON and PND groups (n = 10, each group). *P < 0.05 compared with the CON group, **P < 0.01 compared with the CON group, ****P < 0.005 compared with the CON group.
Figure 9. The validation results of miR-200b-3p and miR-141-3p by qRT-PCR in CON and PND groups (n = 10, each group). ****P < 0.001 and **P < 0.01 compared with the CON group.
Prediction of the Potential Therapeutic Drugs
The differentially expressed mRNAs in the prefrontal cortex of aged mice with PND were upload onto the CMap database to predict the small molecular drugs that have potential therapeutic effects on PND. Negative drug’s mean value indicating the drug has a potential to improve the cognitive function of PND and reverse the expression of differentially expressed genes included in the ceRNA network. The top ten small molecular drugs were shown in Table 4. We found that Cinnarizine and Clemastine had the highest negative mean value. These results indicated that these two compounds might have potential therapeutic effects on PND, and the chemical structure of these two drugs were shown in Figures 10A,B.
Table 4. The top 10 small molecular drugs with high negative correlations with PND.
Figure 10. The chemical structure of potential therapeutic drugs. (A) The chemical structure of Cinnarizine. (B) The chemical structure of Clemastine.
Discussion
Among all the risk factors closely associated with PND, age is the only one factor associated with long-term PND (more than 3 months after operation) (Le et al., 2014). Meanwhile, advanced age is considered as an independent risk factor of PND (Luo et al., 2019), which is not conducive to neurogenesis, neuron survival, and restoration of neuronal function following injury (Koks et al., 2016). Based on these theory, 18-month-old male C57BL/6 mice were used as experimental subjects in this study. The PND animal model was established by exploratory laparotomy plus isoflurane anesthesia, which has been generally proved to be effective in previous studies (Wang et al., 2018a; Qiu et al., 2020). The hippocampus and the cortex, especially the prefrontal cortex, play important roles in regulating cognitive function (Watanabe and Funahashi, 2014; Mansouri et al., 2017; Martorell et al., 2019). In behavioral tests, we found that the spatial learning memory and context-related memory were impaired in aged mice after anesthesia and operation without any post-operative anxious behavior. The phenomenon is regulated by both the hippocampus and the prefrontal cortex (Raybuck and Lattal, 2011), which strongly supports our view that anesthetic and operating factors contribute to the development of PND in aged mice.
As a serious perioperative neurological complication, PND impair patients’ postoperative memory, attention, language comprehension, and social skills, increasing the incidence of dementia and affecting the prognosis of patients (Evered et al., 2016, 2018). Although PND has a profound impact on patients, it is difficult to determine the specific pathogenesis. Many pathological processes have been involved in PND, such as, neuroinflammation (Luo et al., 2019), mitochondrial dysfunction, oxidative stress (Netto et al., 2018), the blood-brain barrier damage (Zhu et al., 2018), neurotrophic support impairment (Fan et al., 2016), and synaptic damage (Xiao et al., 2018). Due to its complex pathogenesis, there is a lack of effective treatments. The application of NGS technology has greatly improved our cognition to identify and diagnose the causes of disease (Rego and Snyder, 2019). Therefore, the NGS technology help us to obtain new clues and treatment targets for PND.
Currently, by RNA sequencing in the hippocampus of aged rats with PND, some studies have identified the differentially expressed circRNAs. And these circRNAs are related to the P53 signaling pathway and the NF-KB pathway (Cao et al., 2020). Like PND, Alzheimer’s disease is also associated with progressive impairment of cognitive function. Differential expression of circRNAs has been found in the hippocampus of Alzheimer’s disease patients and mouse models (Lu et al., 2019; Lo et al., 2020; Ma et al., 2020). These results indicated that the differentially expressed circRNAs were related to cognitive function changes. Moreover, the role of miRNAs in PND has been analyzed in many studies. The miR-572 improves early postoperative cognitive dysfunction by down-regulating NCAM1 (Yu et al., 2015). Meanwhile, the miR-190a improves postoperative cognitive dysfunction in mice by regulating the expression of Tiam1 (Liu et al., 2019). The miRNA-146a protects mice from cognitive impairment caused by operation by inhibiting the neuroinflammation (Chen L. et al., 2019). These results indicate that genes including circRNAs and miRNAs are closely related to PND.
Many circRNAs contain miRNAs binding sites and have the characteristics of miRNA sponge. The inhibition of miRNAs on downstream gene expression is eliminated by binding to target miRNAs. This is the ceRNA mechanism (Thomson and Dinger, 2016), which plays a key role in gene transcription. In this study, 1192 circRNAs and 27 miRNAs with differential expression were identified. Based on our sequencing results, 14 ceRNA pairs were obtained by bioinformatics databases. In addition, the results suggested that not all circRNAs and miRNAs had ceRNA relationships, and the function of most circRNAs need to be further explored. The miRNAs are short non-coding RNAs that regulate gene expression through mRNA degradation or translation inhibition (Zhang et al., 2020). Based on 266 differentially expressed mRNAs, 5 miRNA-mRNA axes were further identified and circRNA-miRNA-mRNA networks were constructed. Genes in these networks, together with other differentially expressed genes, may play key roles in the regulation of the pathogenesis of PND in aged mice.
From our sequencing results, it is not difficult to observe that among all differentially expressed genes, only a small part of genes can form a ceRNA network. In other words, the upstream of differentially expressed mRNAs, such as miRNAs and circRNAs, may not be differentially expressed. This also reflected that the regulation of genes expression in the brain is an extremely complex process.
Gene ontology enrichment and KEGG pathway analysis were performed to explore the roles of differentially expressed mRNAs in the pathogenesis of PND. In the terms of BP, we found that I-kappaB kinase/NF-kappaB signaling was significantly enriched. NF-kappaB signaling is a classical well-known neuroinflammatory signaling pathway, which induces the occurrence and development of PND (Zhu et al., 2017; Cao et al., 2018). The over-activation of NF-kappaB signaling leads to blood-brain barrier damage and the brain injury aggravation (Li Z. et al., 2016). The inhibition of NF-kappaB signaling reduces the neuroinflammation and oxidative stress, which could improve the cognitive impairment associated with PND (Yang and Yuan, 2018). In the terms of MF, we found that postsynaptic neurotransmitter receptor activity was significantly enriched. Postsynaptic neurotransmitter receptor transport is the basic mechanism for the dynamic regulation of synaptic strength, which is crucial for the formation and regulation of neuronal synaptic function (Papadopoulos et al., 2017). The dysregulation of the formation and regulation of neuronal synapse function is closely related to the impaired cognitive function of PND patients (Gao et al., 2021). In the results of KEGG pathway analysis, insulin resistance, pathways of neurodegeneration-multiple diseases, hippo signaling pathway and inflammatory mediator regulation of TRP channels were significantly enriched. Insulin resistance accelerated the development of PND (Tang et al., 2017; He et al., 2019), intranasal insulin treatment prevents anesthesia-induced cognitive impairments (Li X. et al., 2019). Neurodegeneration is one of the main features of PND (Pan et al., 2016). Hippo signaling pathway is involved in neuroinflammation, neuronal cell differentiation, and neuronal death (Cheng et al., 2020), and its intervention might contribute to the treatment of cognitive dysfunction (Yu et al., 2019). TRP channels are expressed in various cell types in the brain and playing a pathological role in a variety of neuroinflammatory diseases (Hu et al., 2020). The inhibitor of TRP channels significantly improve cognitive dysfunction in different diseases (Hazalin et al., 2020; Hu et al., 2020; Thapak et al., 2020). These results reveal the mechanism by which differentially expressed genes in the prefrontal cortex may regulate PND.
Based on the differentially expressed mRNAs, we constructed a PPI network to identify the hub genes of PND. The results indicated that Agt had the largest count of interaction pairs between proteins. The full name of Agt is Angiotensinogen, which is a basic component of the renin-angiotensin system. Its overexpression reflected the activation of renin-angiotensin system, which is closely related to the induction of PND (Li Z. et al., 2014). The expression of Agt was significantly up-regulated in the cerebrospinal fluid of patients with Alzheimer’s disease, which was associated with damage to the blood-brain barrier and impairment of cognitive function (Mateos et al., 2011). The activation of renin-angiotensin system can lead to cognitive decline and brain damage caused by chronic cerebral ischemia (Dong et al., 2011). The identification of hub genes provided potential therapeutic targets for PND.
Meanwhile, Cinnarizine and Clemastine were predicted to have potential therapeutic effects on PND. Cinnarizine is an antihistamine and T-type calcium channel inhibitor that promotes cerebral blood flow. Clomastine is a histamine H1 receptor antagonist. The application of antihistamines could ameliorated cognitive function in mice with Alzheimer’s disease (Steele et al., 2009; Wang et al., 2011). There are many similarities between PND and Alzheimer’s disease in view of pathogenesis, such as oxidative stress, neurodegeneration, and neuroinflammation. It has been reported that Cinnarizine can improve haloperidol-induced memory impairment and oxidative stress (Abdel-Salam et al., 2012). Clemastine shows more potential in cognitive enhancement than Cinnarizine. It could reduce the IL-1β mediated neuroinflammation (Liu et al., 2016; Xie et al., 2020), save mice from age-related memory deficits (Cree et al., 2018; Wang et al., 2020), and reverse spatial learning and memory impairment induced by isoflurane anesthesia (Li Q. et al., 2019). Based on the abovementioned studies, these two drugs, especially Clemastine, might have a potentiality to ameliorate cognitive impairment caused by anesthesia and surgery. However, there is no relevant research has been published so far. This also points out the direction for our follow-up research.
There are also some limitations in this study. First, we cannot identify the source of differentially expressed circRNAs, miRNAs, and mRNAs. Because the prefrontal cortex contains astrocytes, microglia, neurons, and other types of cells. Our sequencing results provided an overall expression trend of the differentially expressed genes. Second, NGS technology has been performed this study in aged mice, the differential expression of genes has been partly verified by qRT-PCR. It remains unclear whether the genes are differentially expressed in homo sapiens brains. Finally, we used bioinformatics analysis to predict the possible mechanisms of differentially expressed genes, and many of the mechanisms that may be related to PND need to be further verified by experiments.
Conclusion
In summary, this study is the first one to explore the differential expression profiles of genes caused by PND. We analyzed the differentially expressed circRNAs, miRNAs, and mRNAs in the prefrontal cortex of aged mice with PND by next-generation sequencing technology and performed some validation work. Meanwhile, through bioinformatics analysis, we predicted the biological functions of differentially expressed genes, constructed a ceRNA network, provided potential targets and small molecule drugs for the treatment of PND.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE174413.
Ethics Statement
The animal study was reviewed and approved by the Institutional Animal Care and Use Committee (Approval No. SYSU-IACUC-2020-000326) and the Laboratory Animal Ethics Committee of Sun Yat-sen University.
Author Contributions
ZW conceived and designed the study and interpreted experiments. WW and XZ performed the experiments and prepared the initial draft of the manuscript. YP, JZ, LC, W-jL, YL, and JW supervised the project. All authors read and approved the final submission.
Funding
This study was supported by grants from the National Natural Science Foundation of China (81571196) and grants from Guangdong Science and Technology Department (2019A1515012147).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as apotential conflict of interest.
Acknowledgments
WW wished to thank, in particular, the patience, care, and support from XZ over the passed years.
Footnotes
- ^ https://github.com/lh3/seqtk
- ^ https://github.com/lh3/bwa
- ^ http://circrna.org/
- ^ https://www.r-project.org/
- ^ http://hannonlab.cshl.edu/fastx_toolkit/index.html
- ^ https://www.perl.org/
- ^ https://www.string-db.org/
- ^ http://starbase.sysu.edu.cn/
- ^ http://bioinformatics.psb.ugent.be/webtools/Venn/
- ^ https://portals.broadinstitute.org/cmap/index.jsp
References
Abdel-Salam, O. M., El-Sayed El-Shamarka, M., Salem, N. A., El-Mosallamy, A. E., and Sleem, A. A. (2012). Amelioration of the haloperidol-induced memory impairment and brain oxidative stress by cinnarizine. EXCLI J. 11, 517–530.
Cao, C., Deng, F., and Hu, Y. (2020). Dexmedetomidine alleviates postoperative cognitive dysfunction through circular RNA in aged rats. 3 Biotech 10:176.
Cao, Y., Li, Z., Ma, L., Ni, C., Li, L., Yang, N., et al. (2018). Isofluraneinduced postoperative cognitive dysfunction is mediated by hypoxiainducible factor1alphadependent neuroinflammation in aged rats. Mol. Med. Rep. 17, 7730–7736.
Chen, C., Gao, R., Li, M., Wang, Q., Chen, H., Zhang, S., et al. (2019). Extracellular RNAs-TLR3 signaling contributes to cognitive decline in a mouse model of postoperative cognitive dysfunction. Brain Behav. Immun. 80, 439–451. doi: 10.1016/j.bbi.2019.04.024
Chen, L., Dong, R., Lu, Y., Zhou, Y., Li, K., Zhang, Z., et al. (2019). MicroRNA-146a protects against cognitive decline induced by surgical trauma by suppressing hippocampal neuroinflammation in mice. Brain Behav. Immun. 78, 188–201. doi: 10.1016/j.bbi.2019.01.020
Cheng, J., Wang, S., Dong, Y., and Yuan, Z. (2020). The role and regulatory mechanism of hippo signaling components in the neuronal system. Front. Immunol. 11:281. doi: 10.3389/fimmu.2020.00281
Cheng, J., Yang, L., Kumar, V., and Agarwal, P. (2014). Systematic evaluation of connectivity map for disease indications. Genome Med. 6:540.
Cree, B. A. C., Niu, J., Hoi, K. K., Zhao, C., Caganap, S. D., Henry, R. G., et al. (2018). Clemastine rescues myelination defects and promotes functional recovery in hypoxic brain injury. Brain 141, 85–98. doi: 10.1093/brain/awx312
Dong, Y. F., Kataoka, K., Toyama, K., Sueta, D., Koibuchi, N., Yamamoto, E., et al. (2011). Attenuation of brain damage and cognitive impairment by direct renin inhibition in mice with chronic cerebral hypoperfusion. Hypertension 58, 635–642. doi: 10.1161/hypertensionaha.111.173534
Evered, L., Silbert, B., Knopman, D. S., Scott, D. A., DeKosky, S. T., Rasmussen, L. S., et al. (2018). Recommendations for the nomenclature of cognitive change associated with Anaesthesia and surgery-2018. Anesthesiology 129, 872–879. doi: 10.1097/aln.0000000000002334
Evered, L. A., Silbert, B. S., Scott, D. A., Maruff, P., and Ames, D. (2016). Prevalence of dementia 7.5 years after coronary artery bypass graft surgery. Anesthesiology 125, 62–71. doi: 10.1097/aln.0000000000001143
Fan, D., Li, J., Zheng, B., Hua, L., and Zuo, Z. (2016). Enriched environment attenuates surgery-induced impairment of learning, memory, and neurogenesis possibly by preserving BDNF expression. Mol. Neurobiol. 53, 344–354. doi: 10.1007/s12035-014-9013-1
Gao, R., Chen, C., Zhao, Q., Li, M., Wang, Q., Zhou, L., et al. (2020). Identification of the potential key circular RNAs in elderly patients with postoperative cognitive dysfunction. Front. Aging Neurosci. 12:165. doi: 10.3389/fnagi.2020.00165
Gao, S., Zhang, S., Zhou, H., Tao, X., Ni, Y., Pei, D., et al. (2021). Role of mTOR-regulated autophagy in synaptic plasticity related proteins downregulation and the reference memory deficits induced by anesthesia/surgery in aged mice. Front. Aging Neurosci. 13:628541. doi: 10.3389/fnagi.2021.628541
Gao, Y., Wang, J., and Zhao, F. (2015). CIRI: an efficient and unbiased algorithm for de novo circular RNA identification. Genome Biol. 16:4.
Hazalin, N., Liao, P., and Hassan, Z. (2020). TRPM4 inhibition improves spatial memory impairment and hippocampal long-term potentiation deficit in chronic cerebral hypoperfused rats. Behav. Brain Res. 393:112781. doi: 10.1016/j.bbr.2020.112781
He, X., Long, G., Quan, C., Zhang, B., Chen, J., and Ouyang, W. (2019). Insulin resistance predicts postoperative cognitive dysfunction in elderly gastrointestinal patients. Front. Aging Neurosci. 11:197. doi: 10.3389/fnagi.2019.00197
Hu, H., Zhu, T., Gong, L., Zhao, Y., Shao, Y., Li, S., et al. (2020). Transient receptor potential melastatin 2 contributes to neuroinflammation and negatively regulates cognitive outcomes in a pilocarpine-induced mouse model of epilepsy. Int. Immunopharmacol. 87:106824. doi: 10.1016/j.intimp.2020.106824
Jakel, S., and Williams, A. (2020). What have advances in transcriptomic technologies taught us about human white matter pathologies? Front. Cell. Neurosci. 14:238. doi: 10.3389/fncel.2020.00238
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Koks, S., Dogan, S., Tuna, B. G., Gonzalez-Navarro, H., Potter, P., and Vandenbroucke, R. E. (2016). Mouse models of ageing and their relevance to disease. Mech. Ageing Dev. 160, 41–53. doi: 10.1016/j.mad.2016.10.001
Kozomara, A., and Griffiths-Jones, S. (2014). miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 42, D68–D73.
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25.
Le, Y., Liu, S., Peng, M., Tan, C., Liao, Q., Duan, K., et al. (2014). Aging differentially affects the loss of neuronal dendritic spine, neuroinflammation and memory impairment at rats after surgery. PLoS One 9:e106837. doi: 10.1371/journal.pone.0106837
Li, J., Zhu, X., Yang, S., Xu, H., Guo, M., Yao, Y., et al. (2019). Lidocaine attenuates cognitive impairment after isoflurane anesthesia by reducing mitochondrial damage. Neurochem. Res. 44, 1703–1714. doi: 10.1007/s11064-019-02799-0
Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 42, D92–D97.
Li, M., Chen, C., Zhang, W., Gao, R., Wang, Q., Chen, H., et al. (2019). Identification of the potential key long non-coding RNAs in aged mice with postoperative cognitive dysfunction. Front. Aging Neurosci. 11:181. doi: 10.3389/fnagi.2019.00181
Li, Q., Mathena, R. P., Xu, J., Eregha, O. N., Wen, J., and Mintz, C. D. (2019). Early postnatal exposure to isoflurane disrupts oligodendrocyte development and myelin formation in the Mouse Hippocampus. Anesthesiology 131, 1077–1091. doi: 10.1097/aln.0000000000002904
Li, X., Run, X., Wei, Z., Zeng, K., Liang, Z., Huang, F., et al. (2019). Intranasal insulin prevents anesthesia-induced cognitive impairments in aged mice. Curr. Alzheimer Res. 16, 8–18. doi: 10.2174/1567205015666181031145045
Li, Y., Zheng, Q., Bao, C., Li, S., Guo, W., Zhao, J., et al. (2015). Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. 25, 981–984. doi: 10.1038/cr.2015.82
Li, Z., Cao, Y., Li, L., Liang, Y., Tian, X., Mo, N., et al. (2014). Prophylactic angiotensin type 1 receptor antagonism confers neuroprotection in an aged rat model of postoperative cognitive dysfunction. Biochem. Biophys. Res. Commun. 449, 74–80. doi: 10.1016/j.bbrc.2014.04.153
Li, Z., Mo, N., Li, L., Cao, Y., Wang, W., Liang, Y., et al. (2016). Surgery-induced hippocampal angiotensin II elevation causes blood-brain barrier disruption via MMP/TIMP in aged rats. Front. Cell. Neurosci. 10:105. doi: 10.3389/fncel.2016.00105
Lin, D., and Zuo, Z. (2011). Isoflurane induces hippocampal cell injury and cognitive impairments in adult rats. Neuropharmacology 61, 1354–1359. doi: 10.1016/j.neuropharm.2011.08.011
Lin, F., Zheng, Y., Pan, L., and Zuo, Z. (2020). Attenuation of noisy environment-induced neuroinflammation and dysfunction of learning and memory by minocycline during perioperative period in mice. Brain Res. Bull. 159, 16–24. doi: 10.1016/j.brainresbull.2020.03.004
Liu, J., Dupree, J. L., Gacias, M., Frawley, R., Sikder, T., Naik, P., et al. (2016). Clemastine enhances myelination in the prefrontal cortex and rescues behavioral changes in socially isolated mice. J. Neurosci. 36, 957–962. doi: 10.1523/jneurosci.3608-15.2016
Liu, Q., Hou, A., Zhang, Y., Guo, Y., Li, J., Yao, Y., et al. (2019). MiR-190a potentially ameliorates postoperative cognitive dysfunction by regulating Tiam1. BMC Genomics 20:670. doi: 10.1186/s12864-019-6035-0
Lo, I., Hill, J., Vilhjalmsson, B. J., and Kjems, J. (2020). Linking the association between circRNAs and Alzheimer’s disease progression by multi-tissue circular RNA characterization. RNA Biol. 17, 1789–1797. doi: 10.1080/15476286.2020.1783487
Lu, Y., Tan, L., and Wang, X. (2019). Circular HDAC9/microRNA-138/Sirtuin-1 pathway mediates synaptic and amyloid precursor protein processing deficits in Alzheimer’s disease. Neurosci. Bull. 35, 877–888. doi: 10.1007/s12264-019-00361-0
Luo, A., Yan, J., Tang, X., Zhao, Y., Zhou, B., and Li, S. (2019). Postoperative cognitive dysfunction in the aged: the collision of neuroinflammaging with perioperative neuroinflammation. Inflammopharmacology 27, 27–37. doi: 10.1007/s10787-018-00559-0
Ma, N., Pan, J., Ye, X., Yu, B., Zhang, W., and Wan, J. (2019). Whole-transcriptome analysis of APP/PS1 mouse brain and identification of circRNA-miRNA-mRNA networks to investigate AD pathogenesis. Mol. Ther. Nucleic Acids 18, 1049–1062. doi: 10.1016/j.omtn.2019.10.030
Ma, N., Tie, C., Yu, B., Zhang, W., and Wan, J. (2020). Circular RNAs regulate its parental genes transcription in the AD mouse model using two methods of library construction. FASEB J. 34, 10342–10356. doi: 10.1096/fj.201903157r
Mansouri, F. A., Koechlin, E., Rosa, M. G. P., and Buckley, M. J. (2017). Managing competing goals - a key role for the frontopolar cortex. Nat. Rev. Neurosci. 18, 645–657. doi: 10.1038/nrn.2017.111
Martorell, A. J., Paulson, A. L., Suk, H. J., Abdurrob, F., Drummond, G. T., Guan, W., et al. (2019). Multi-sensory gamma stimulation ameliorates Alzheimer’s-associated pathology and improves cognition. Cell 177, 256–271.e22.
Mateos, L., Ismail, M. A., Gil-Bea, F. J., Leoni, V., Winblad, B., Bjorkhem, I., et al. (2011). Upregulation of brain renin angiotensin system by 27-hydroxycholesterol in Alzheimer’s disease. J. Alzheimers Dis. 24, 669–679. doi: 10.3233/jad-2011-101512
Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., and Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621–628. doi: 10.1038/nmeth.1226
Netto, M. B., de Oliveira Junior, A. N., Goldim, M., Mathias, K., Fileti, M. E., da Rosa, N., et al. (2018). Oxidative stress and mitochondrial dysfunction contributes to postoperative cognitive dysfunction in elderly rats. Brain Behav. Immun. 73, 661–669. doi: 10.1016/j.bbi.2018.07.016
Pan, K., Li, X., Chen, Y., Zhu, D., Li, Y., Tao, G., et al. (2016). Deferoxamine pre-treatment protects against postoperative cognitive dysfunction of aged rats by depressing microglial activation via ameliorating iron accumulation in hippocampus. Neuropharmacology 111, 180–194. doi: 10.1016/j.neuropharm.2016.09.004
Papadopoulos, T., Rhee, H. J., Subramanian, D., Paraskevopoulou, F., Mueller, R., Schultz, C., et al. (2017). Endosomal phosphatidylinositol 3-phosphate promotes gephyrin clustering and GABAergic neurotransmission at inhibitory postsynapses. J. Biol. Chem. 292, 1160–1177. doi: 10.1074/jbc.m116.771592
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi: 10.1038/nprot.2016.095
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Qiu, L. L., Pan, W., Luo, D., Zhang, G. F., Zhou, Z. Q., Sun, X. Y., et al. (2020). Dysregulation of BDNF/TrkB signaling mediated by NMDAR/Ca(2+)/calpain might contribute to postoperative cognitive dysfunction in aging mice. J. Neuroinflammation 17:23.
Raybuck, J. D., and Lattal, K. M. (2011). Double dissociation of amygdala and hippocampal contributions to trace and delay fear conditioning. PLoS One 6:e15982. doi: 10.1371/journal.pone.0015982
Rego, S. M., and Snyder, M. P. (2019). High throughput sequencing and assessing disease risk. Cold Spring Harb. Perspect. Med. 9:a026849. doi: 10.1101/cshperspect.a026849
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Robinson, M. D., and Oshlack, A. (2010). A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11:R25.
Ru, Y., Kechris, K. J., Tabakoff, B., Hoffman, P., Radcliffe, R. A., Bowler, R., et al. (2014). The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res. 42:e133. doi: 10.1093/nar/gku631
Steele, J. W., Kim, S. H., Cirrito, J. R., Verges, D. K., Restivo, J. L., Westaway, D., et al. (2009). Acute dosing of latrepirdine (Dimebon), a possible Alzheimer therapeutic, elevates extracellular amyloid-beta levels in vitro and in vivo. Mol. Neurodegener. 4:51. doi: 10.1186/1750-1326-4-51
Tang, N., Jiang, R., Wang, X., Wen, J., Liu, L., Wu, J., et al. (2017). Insulin resistance plays a potential role in postoperative cognitive dysfunction in patients following cardiac valve surgery. Brain Res. 1657, 377–382. doi: 10.1016/j.brainres.2016.12.027
Thapak, P., Bishnoi, M., and Sharma, S. S. (2020). Amelioration of diabetes-induced cognitive impairment by Transient Receptor Potential Vanilloid 2 (TRPV2) channel inhibitor: Behavioral and mechanistic study. Neurochem. Int. 139:104783. doi: 10.1016/j.neuint.2020.104783
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA sponges: evidence and controversy. Nat. Rev. Genet. 17, 272–283. doi: 10.1038/nrg.2016.20
Urits, I., Orhurhu, V., Jones, M., Hoyt, D., Seats, A., and Viswanath, O. (2019). Current perspectives on postoperative cognitive dysfunction in the ageing population. Turk. J. Anaesthesiol. Reanim. 47, 439–447.
Wang, F., Ren, S. Y., Chen, J. F., Liu, K., Li, R. X., Li, Z. F., et al. (2020). Myelin degeneration and diminished myelin renewal contribute to age-related deficits in memory. Nat. Neurosci. 23, 481–486. doi: 10.1038/s41593-020-0588-8
Wang, J., Ferruzzi, M. G., Varghese, M., Qian, X., Cheng, A., Xie, M., et al. (2011). Preclinical study of dimebon on beta-amyloid-mediated neuropathology in Alzheimer’s disease. Mol. Neurodegener. 6:7. doi: 10.1186/1750-1326-6-7
Wang, Z., Meng, S., Cao, L., Chen, Y., Zuo, Z., and Peng, S. (2018a). Critical role of NLRP3-caspase-1 pathway in age-dependent isoflurane-induced microglial inflammatory response and cognitive impairment. J. Neuroinflammation 15:109.
Wang, Z., Xu, P., Chen, B., Zhang, Z., Zhang, C., Zhan, Q., et al. (2018b). Identifying circRNA-associated-ceRNA networks in the hippocampus of Abeta1-42-induced Alzheimer’s disease-like rats using microarray analysis. Aging 10, 775–788. doi: 10.18632/aging.101427
Watanabe, K., and Funahashi, S. (2014). Neural mechanisms of dual-task interference and cognitive capacity limitation in the prefrontal cortex. Nat. Neurosci. 17, 601–611. doi: 10.1038/nn.3667
Xiao, J. Y., Xiong, B. R., Zhang, W., Zhou, W. C., Yang, H., Gao, F., et al. (2018). PGE2-EP3 signaling exacerbates hippocampus-dependent cognitive impairment after laparotomy by reducing expression levels of hippocampal synaptic plasticity-related proteins in aged mice. CNS Neurosci. Ther. 24, 917–929. doi: 10.1111/cns.12832
Xie, D., Ge, X., Ma, Y., Tang, J., Wang, Y., Zhu, Y., et al. (2020). Clemastine improves hypomyelination in rats with hypoxic-ischemic brain injury by reducing microglia-derived IL-1beta via P38 signaling pathway. J. Neuroinflammation 17:57.
Xie, Z., McAuliffe, S., Swain, C. A., Ward, S. A., Crosby, C. A., Zheng, H., et al. (2013). Cerebrospinal fluid abeta to tau ratio and postoperative cognitive change. Ann. Surg. 258, 364–369. doi: 10.1097/sla.0b013e318298b077
Yang, Z. Y., and Yuan, C. X. (2018). IL-17A promotes the neuroinflammation and cognitive function in sevoflurane anesthetized aged rats via activation of NF-kappaB signaling pathway. BMC Anesthesiol. 18:147. doi: 10.1186/s12871-018-0607-4
You, X., Vlatkovic, I., Babic, A., Will, T., Epstein, I., Tushev, G., et al. (2015). Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat. Neurosci. 18, 603–610. doi: 10.1038/nn.3975
Yu, F., Han, W., Zhan, G., Li, S., Jiang, X., Xiang, S., et al. (2019). Differential levels of hippo signaling in selected brain and peripheral tissues in streptozotocin-induced cognitive dysfunction in mice. Neuroscience 421, 48–58. doi: 10.1016/j.neuroscience.2019.09.018
Yu, X., Liu, S., Li, J., Fan, X., Chen, Y., Bi, X., et al. (2015). MicroRNA-572 improves early post-operative cognitive dysfunction by down-regulating neural cell adhesion molecule 1. PLoS One. 10:e0118511. doi: 10.1371/journal.pone.0118511
Zhang, B., Tian, L., Xie, J., Chen, G., and Wang, F. (2020). Targeting miRNAs by natural products: a new way for cancer therapy. Biomed. Pharmacother. 130:110546. doi: 10.1016/j.biopha.2020.110546
Zhang, S., Zhu, D., Li, H., Li, H., Feng, C., and Zhang, W. (2017). Characterization of circRNA-Associated-ceRNA networks in a senescence-accelerated mouse prone 8 brain. Mol. Ther. 25, 2053–2061. doi: 10.1016/j.ymthe.2017.06.009
Zhu, H., Liu, W., and Fang, H. (2018). Inflammation caused by peripheral immune cells across into injured mouse blood brain barrier can worsen postoperative cognitive dysfunction induced by isoflurane. BMC Cell Biol. 19:23. doi: 10.1186/s12860-018-0172-1
Zhu, Y., Wang, Y., Yao, R., Hao, T., Cao, J., Huang, H., et al. (2017). Enhanced neuroinflammation mediated by DNA methylation of the glucocorticoid receptor triggers cognitive dysfunction after sevoflurane anesthesia in adult rats subjected to maternal separation during the neonatal period. J. Neuroinflammation 14:6.
Keywords: postoperative neurocognitive disorder (PND), bioinformatics, competitive endogenous RNA network, function enrichment analysis, therapeutic target
Citation: Wu W, Peng Y, Zhou J, Zhang X, Cao L, Lin W-j, Lu Y, Wen J and Wang Z (2021) Identification of the Potential Gene Regulatory Networks and Therapeutics in Aged Mice With Postoperative Neurocognitive Disorder. Front. Neurosci. 15:689188. doi: 10.3389/fnins.2021.689188
Received: 31 March 2021; Accepted: 28 May 2021;
Published: 24 June 2021.
Edited by:
Mingshu Mo, The First Affiliated Hospital of Guangzhou Medical University, ChinaCopyright © 2021 Wu, Peng, Zhou, Zhang, Cao, Lin, Lu, Wen and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhi Wang, wangzhi@mail.sysu.edu.cn
†These authors have contributed equally to this work