Original Research ARTICLE
Microarray Expression Profiles of lncRNAs and mRNAs in Postoperative Cognitive Dysfunction
- 1Department of Anesthesiology, Affiliated Traditional Chinese Medicine Hospital, Southwest Medical University, Luzhou, China
- 2School of Pharmacy and Medical Sciences, Sansom Institute, Division of Health Sciences, University of South Australia, Adelaide, SA, Australia
- 3Department of Cardiothoracic Surgery, Affiliated Hospital, Southwest Medical University, Luzhou, China
Postoperative cognitive dysfunction (POCD) is serious disorder in the central nervous system common in aged patients after anesthesia. Although its clinical symptoms are well recognized, however, the molecular etiology of the POCD remains unrevealed. Similarly, neither gold standard molecular diagnosis nor effective treatment is available for POCD until the present. Therefore, we aimed to explore the molecular mechanism of this disorder through investigating lncRNAs and mRNAs associated with POCD human patients and investigate their underlying regulatory pathways. In this study, we recruited 200 patients requiring hip or knee replacement surgery. Their neurological functions were assessed at two time points, 1 day before the surgery and 30 days post-surgery. In parallel, serum samples were collected from the participants to analyze lncRNAs and mRNAs differential expression profile between POCD and non-POCD patients using microarray analysis. To further investigate the role differentially expressed mRNA and lncRNAs, Gene Ontology (GO), pathway analyses on mRNAs and lncRNA-mRNA interaction network were performed. As a result, 68 lncRNAs and 115 mRNAs were dysregulated in the POCD group compared to non-POCD group. Among them, the top 10 upregulated lncRNAs and 10 downregulated lncRNAs were listed for enrichment analysis. Interestingly, we found that these lncRNA and mRNA are involved in biological process, molecular function, and cellular component in addition to various signaling pathways, suggesting that the pathogenesis of POCD involves lncRNAs and mRNAs differential expression. Consequently, the genetic dysregulation between the non-POCD and POCD patients participates in the occurrence and development of POCD, and could be served as diagnostic biomarkers and drug targets for POCD treatment.
Postoperative cognitive dysfunction (POCD) is a severe central nervous system (CNS) complication in the cognitive function common in elderly patients after anesthesia and surgery. The condition involves symptoms rang from mild to serious cognitive deficits affecting life quality (Hartholt et al., 2012; Arora et al., 2014). POCD is clinically characterized by different levels of memory lost, emotional disorder such as anxiety, confusion, and depression, and in extreme cases, coma and death (Bekker and Weeks, 2003). While the neuropsychological test for POCD typically includes Mini-Mental State Examination (MMSE), Wechsler Adult Memory Scale (WMS), and Minnesota Multiphasic Personality Inventory (MMPI; Moller et al., 1998; Krenk et al., 2010; Monk and Price, 2011). Owing to differences in definitions and testing methods, the incidence of POCD has been reported globally with high prevalence. Generally, the incidence of POCD was estimated up to 26% at 7 days and 10% at 3 months in beyond than 65 years old patients (van Harten et al., 2012; Silbert et al., 2014). Disease pathogenicity is complex and multifactorial. It includes central cholinergic system degradation, central inflammation (Hu et al., 2014), tau phosphorylation (Planel et al., 2007), β-amyloid deposition (Xu et al., 2012), and neuronal apoptosis (Vacas et al., 2014). Although POCD occurs in high prevalence and has severe effect on the affected patients (Witlox et al., 2010). However, there is no standard treatment of the disease or measures that can be taken to prevent or reduce the incidence of POCD (Pappa et al., 2017).
Long non-coding RNAs (lncRNAs) are set of non-coding RNA with more than 200 nucleotides lengths (Marchese et al., 2017; Muret et al., 2017). Interestingly, increasing number of publications presented evidences of critical role in various biological processes, mainly on chromosome dosage compensation, genomic imprinting, target mimicry, and functional protein trafficking (Clark and Mattick, 2011; Nagano and Fraser, 2011; Wang and Chang, 2011). In addition, lncRNAs involve in various biological function, such as cell proliferation, differentiation, apoptosis, and regulating the expression of target genes during transcription, post-transcription, and epigenetic levels (Yoon et al., 2013). At the present, researchers mainly relay on computational and mathematical methods to predict lncRNA functional domains at structural levels. Among them, chromatin isolation by RNA purification (Ch IRP) and RNA-binding protein immunoprecipitation (RIP) have widely used to investigate the interaction between lncRNA with genes and proteins respectively (Yoon et al., 2013). Moreover, lncRNAs have been confirmed to associate with number of neurological disorders in addition to many other diseases like cancer and immune diseases (Ponting et al., 2009; Ulitsky and Bartel, 2013; Maass et al., 2014; Isin and Dalay, 2015). Nevertheless, the relation of lncRNAs and POCD is unclear and needs to be explored. Therefore, this study is to investigate the role of lncRNA and mRNA in the development of POCD.
We selected hip and knee replacement patients as clinical cases to detect the lncRNAs and mRNAs profiles in POCD patients non-POCD patients using microarray analysis. Furthermore, Gene Ontology (GO) and KEGG analyses were performed to study the potential functions of the dysregulated lncRNAs and mRNA. Finally, network interaction of lncRNAs–mRNAs in the POCD patients was constructed.
Materials and Methods
Ethical Considerations and Study Approval
The study was approved by the ethic committee of Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University (project number KY2018001) and registered with Chinese Ethics Committee of Registering Clinical Trials in 2018 (Unique identifier ChiCTR-INR-15007607). All participating patients were informed with the aim of the study and signed the related consent forms.
Selection of Patients
Two hundred patients admitted to Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University to carry out hip or knee replacement surgery have participated in this study. Patients selected based on firm criteria which were (1) age at >60 years old; (2) anesthetic procedure (American Society of Anesthesiologists Grade 1–3); and (3) mental health. It was decided to exclude patients beyond 85 years old, patients with nervous system diseases, preoperative stoke, and other central nerve system diseases cases (preoperative MMSE < 24 scores), and patients who suffered from multiple injuries and craniocerebral injury.
Anesthesia and Surgery
All patients were asked to fast from food 8 h before the operation and from water 3 h before anesthesia. The anesthetic procedure was combined of spinal and epidural anesthesia at the L3–L4 intervertebral disc in a side-lying position using an 18-gage needle and 25-gage spinal needle (Tuo Ren, Henan, China). Then, the 2 ml mixed solution consisting of 0.5% bupivacaine (Hua Lu, Shandong, China) and 3.33% glucose (Shuang He, Beijing, China) was slowly injected to subarachnoid space. The vital sign of patients was monitored during the surgery. All the participants were given the same analgesia doses 12 μg/kg fentanyl (Ren Fu, Yichang, China) and 10 mg tropisetron hydrochloride (Ren Fu, Yichang, China) for 2 days.
For cognitive function evaluation, a battery of neuropsychological tests were followed 1 day before and 1 month after the surgery by trained physician. The tests consisted of the MMSE, Symbol Digit Modalities Test (SDMT), Digit Span Test (DST), Trail Making Test A, and Clock Drawing Test (CDT). The cognitive functions were principally assessed including memory, learning concentration, and executive function. Briefly, MMSE was designed to evaluate and quantify the global cognitive state, assess memory, spatial–temporal orientation, and reasoning (Bertolucci et al., 1994). SDMT was used to test psychomotor speed and attentional control (Valentin et al., 2016). DST is a subtest of WAIS-RC, which instructed participants to memorize numbers, and repeat them reversely and was used for memory assessment (Gong, 1992). Trail Making Test A was adopted to assess executive function, selective and alternating attention (Valentin et al., 2016). While CDT was applied for visual space function, concentration, the planning, and execution function of the action assessment (Aprahamian et al., 2010).
A group 30 health volunteers were also assessed as normal reference of cognitive functions and were corresponding with their patient’s counterparts based on education level, age, and gender. The cognitive functions were evaluated twice in 7 days intervals. Patients were separated into POCD and non-POCD groups. The (SD) value, mean variation, and baseline score were calculated to statistically validate the data (Rasmussen et al., 2001). The learning effect and baseline score were deducted from the score of neuropsychological tests (Duan et al., 2018) and, divided by the SD of the baseline score of the healthy group. This result was defined as Z score, and the participants Z ≥ 1.96 were included in POCD (Silbert et al., 2014).
Three days after the surgery, 8 ml of peripheral blood samples were collected from the patients in pro-coagulant tube and centrifuged (at speed of 2500 rpm for 10 min). The supernatant was taken and stored in liquid nitrogen until the test. Random number table was applied to select three patients with early POCD (labeled A, B, C) and three patients with non-POCD was as control (labeled D, E, F). Total RNA was extracted using TRIzol reagent kit based on vendor instructions, the plasma. RNA integrity was assessed by Agilent Bioanalyzer 2100 (Agilent Technologies), and the purity was measured with NanoDrop2000 spectrophotometer (Thermo Scientific). RNA samples with OD260/280 > 1.8 and OD260/230 > 1.5 were taken for further testing.
RNA samples were sent to Shanghai OE Biotech Inc. (Shanghai, China) to carry out. The sample labeling, hybridization, and washing were performed according to service provider protocols. In brief, cDNAs was synthesized for from RNA samples then synthesized cRNAs. Second cycle cDNAs were synthesized from cRNAs and followed by fragmentation and biotin labeling, the second cycle cDNAs were hybridized onto the Affymetrix Human OE lncRNA microarray (Affymetrix, Santa Clara, CA, United States), comprising 66,741 human lncRNAs and 25,986 coding mRNA. The chips were then washed, stained, and scanned by Affymetrix Scanner 3000 (Affymetrix, Santa Clara, CA, United States).
Affymetrix Gene Chip Command Console (version 4.0, Affymetrix) software was used to extract raw data. Expression Console (version1.3.1, Affymetrix) software offered RMA normalization for both gene and exon level analysis. Then, Genesrping software (version 13.1; Agilent Technologies) was employed to carry out the basic analysis. Venn Diagram software was performed to generate high-resolution Venn plots.
To understand the function of lncRNAs associated with the PODC, we first established the role of differentially expressed mRNA through GO and KEGG analyses. The lncRNAs then were linked to their adjacent mRNAs to establish their potential function (Zhao et al., 2015). Pearson correlation and p-values were calculated and followed to measure pathway correlation to the conditions. Method of Storey (2002) was used to identify the false discovery rate. Venn diagram was applied to display the co-regulated mRNAs of lncRNAs. GO analysis was used for further insight on the biological function and molecular mechanism of the dysregulated genes. The workflow of the study is described in Figure 1A.
FIGURE 1. Samples self-test. (A) The workflow of the study. (B) Box–whisker plot shows the dispersion. (C) Scatter plot shows the total trend of dispersion. (D) PAC shows the rationality of POCD group (labeled A, B, C) and non-POCD group (labeled D, E, F). POCD, postoperative cognitive dysfunction. A, B, and C in the images represent three samples in POCD group, and D, E, and F represent three samples in non-POCD group.
In this study, SPSS software version 17.0 (SPSS, Inc., Chicago, IL, United States) and GraphPad Prism 5.0 (GraphPad Software, Inc., La Jolla, CA, United States) were used for data analysis and the experimental results were presented as mean values ± standard deviation (SD); P < 0.05 was viewed as statistically significant difference.
Characteristics and Neuropsychological Test of Participates
Of these 200 patients, 181 patients (91.5%) achieved the whole tests, while 19 patients (9.5%) did not complete the tests due personal reasons. The general characteristics of the patients including age, body mass index, gender, and education are presented in Table 1. No significance variations were noticed in the scores of the all participants at first day before surgery (P > 0.05). On the other hand, CDT scores were noticeably lower in some patients and they were taken as POCD group compared with non-POCD group after 1 month of the surgery (Table 2). Whereas, the general characteristics and intraoperative conditions were statistically insignificant. One month after hip or knee replacement surgery, 21 patients (11.6%) were diagnosed as cognitive dysfunction, performed remarkably worse verbal IQ, performance IQ, and working memory compared with other non-POCD patients (Table 3).
TABLE 1. General characteristics and cognitive functions of participants from study and control groups at baseline
TABLE 2. General characteristics, intraoperative indicators and cognitive functions of patients from non-POCD and POCD at baseline.
Samples Quality Detection
Chip data dispersion was illustrated in box–whisker plot, scatter plot, and PCA. From box–whisker plot, no statistical variance was seen neither in the lncRNAs distributions nor in mRNAs expression profile in the samples (Figure 1B). As shown in Figures 1C,D, the scatter plot and PCA were applied to measure dissimilarities in lncRNAs expression between POCD group and control group, the detection revealed that in scatter plot and PCA, the overall distribution of the two groups of data is concentrated (Figures 1C,D).
Differential Expression Profile of lncRNAs
Hierarchical clustering analysis data and heatmap (Figure 2A) have shown the distinct lncRNAs expression profile between the groups; 68 lncRNAs were observed to be differentially expressed in the POCD patients, including 29 and 39 upregulated and downregulated, respectively (see Supplementary Tables S1, S2). The top 10 upregulated lncRNAs were identified as lnc-FAM53B-2_1, ENST00000615844.1, lnc-MYL2-2_1, lnc-NKD2-1_1, lnc-EVC-3_1, lnc-SLC5A7-1_1, NR_024101, lnc-DCP1B-6_1, lnc-MNAT1-8_1, and lnc-TNFRSF19-4_1 (Figure 2B). Whereas the top 10 downregulated lncRNAs were NR_024187, lnc-PTPRU-7_1, lnc-LIPC-6_1, lnc-CTD-2144E22.5.1-26_1, lnc-OTX1-7_1, ENST00000618386.1, lnc-PRL-9_1, lnc-PAG1-9_1, lnc-ZNF729-4_1, and lnc-RAP2C-5_1 (Figure 2C). These lncRNAs were further studied by bioinformatics tools.
FIGURE 2. Differential expression of lncRNAs. (A) Heatmap of lncRNAs in POCD group vs. non-POCD group. (B) The top 10 upregulated lncRNAs. (C) The top 10 downregulated lncRNAs. POCD, postoperative cognitive dysfunction.
Bioinformatics Analysis of Upregulated mRNAs in POCD Group
Gene Ontology analysis for the main overexpressed mRNA has revealed that these transcriptomes are involved in biological processes through series of functions including deacetylation, cytoplasmic mRNA processing, stem cell population maintenance, and positive regulation of immunity against viruses by vesicle organization, regulation of mitotic cell cycle, peptidyl-threonine phosphorylation, endocytosis, histone deacetylation, and actin cytoskeleton reorganization (Figure 3A). At the same time, these mRNA sequences were predicted to target different cellular components including cytoplasmic stress granule, transcriptional repressor complex, clathrin-coated vesicle, NatA complex, cullin-RING ubiquitin ligase complex, cytosolic large ribosomal subunit, Cul5-RING ubiquitin ligase complex, Cdc73/Paf1 complex, clathrin complex, condensed chromosome, and outer kinetochore (Figure 3B). On the molecular level, the other upregulated RNAs were involved in modulating the cells by being nucleosomal DNA binding, RNA binding factors, protein deacetylase activity, poly(A), RNA polymerase II repressing factor, serine/threonine kinase, signaling receptors, interleukin-12 receptor binding factor, and 1-phosphatidylinositol-3-phosphate 4-kinase activity (Figure 3C). On other dimension, KEGG data of the upregulated mRNAs were predicted to take part in transcriptional misregulation in cancer, Epstein–Barr virus infection, pantothenate and CoA biosynthesis, endocytosis, spliceosome, RNA polymerase, protein processing in endoplasmic reticulum, alcoholism, N-glycan biosynthesis, and Huntington’s disease (Figure 3D).
FIGURE 3. Analysis of bioinformatics for the upregulated genes. (A) The map of the top 10 biological process. (B) The map of the top 10 cell components. (C) The map of the top 10 molecular functions. (D) The top 10 significant pathways of upregulated genes in POCD group.
Bioinformatics Analysis of Downregulated mRNAs in POCD Group
On the basis of the enrichment scores, the top 10 biological processes targeted by downregulated mRNAs were identified as cytoplasmic mRNA processing body assembly, snRNA transcription, polymerase II promoter, cell division, regulation of mitotic cell cycle, outflow tract morphogenesis, stem cell population maintenance, regulation of translation, regulation of fertilization, meiotic chromosome condensation, and regulation of myotube differentiation (Figure 4A). Whereas the highest 10 cellular components related to the downregulated mRNAs were cytoplasmic stress granule, Cajal body, condensin complex, cullin-RING ubiquitin ligase complex, cytosolic large ribosomal subunit, Cul5-RING ubiquitin ligase complex, Cdc73/Paf1 complex, nucleoplasm, outer dense fiber, and condensed chromosome outer kinetochore (Figure 4B). Similarly, the top 10 functions on molecular mechanism include poly(A) RNA binding, nucleosomal DNA binding, signaling protein receptor, serine/threonine kinase, phosphatase regulator, SUMO-specific protease activity, protein serine/threonine kinase activity, ATP binding, RNA polymerase II intronic transcription regulator, neuropeptide Y receptor activity, and structural constituent of ribosome (Figure 4C). Laterally, the top 10 affected pathways of mRNA sequences downregulated were Epstein–Barr virus infection, pantothenate and CoA biosynthesis, spliceosome, RNA polymerase, Protein processing in endoplasmic reticulum, transcriptional misregulation in cancer, alcoholism, N-Glycan biosynthesis, carcinogenesis viral RNA, and adipocytes lipolysis regulation (Figure 4D).
FIGURE 4. Analysis of bioinformatics for the downregulated genes. (A) The map of the top 10 biological process. (B) The map of the top 10 cell components. (C) The map of the top 10 molecular functions. (D) The top 10 significant pathways of downregulated genes in POCD group.
Co-expressed Analysis of lncRNA and mRNA in POCD Group
Based on the information of trans (trans-prediction) and Venn analyses, the lncRNA–mRNA interaction network was constructed. The data showed that the 10 upregulated lncRNAs were correlated with amounts of mRNAs through trans analysis. The lnc-FAM53B-2_1 in particular was the highest upregulated lncRNA. Interestingly, it interacted with the majority of differentially dysregulated mRNA sequences. It was also noticed that SMARCE1, ATXN2, SMC4, and IL23R were remarkably unregulated in POCD patients (Figure 5A). Venn diagram showed that the highest five upregulated lncRNAs were lnc-FAM53B-2_1, ENST00000615844.1, lnc-MYL2-2_1, lnc-NKD2-1_1, and lnc-EVC-3_1. While the main overexpressed mRNAs were AC009487.4, FANK1, OLFML2B, OXSR1, RNU6-18P, RP11-483E7.1, SMARCE1, and ZFR (Figure 5B).
FIGURE 5. Analysis of upregulated lncRNAs-mRNA co-expression. (A) The top 10 upregulated of regulatory network. Blue nodes represent lncRNAs and yellow nodes represent target genes. The pink lines represent positive correlations and the green lines represent negative correlations. (B) Top five upregulated lncRNAs co-regulated mRNA.
In addition, the top 10 downregulated lncRNAs were also associated with many mRNA sequences. As the most downregulated lncRNA in POCD, NR_024187 was correlated with that of NEK3, XKRX, VNN1, TBXAS1, and HDAC5 (Figure 6A). Venn diagram also showed that the main five downregulated lncRNAs, were NR_024187, lnc-PTPRU-7_1, lnc-LIPC-6_1, lnc-CTD-2144E22.5.1-26_1, and lnc-OTX1-7_1 and they were slightly correlated with RP5-834N19.1 and PARL mRNA (Figure 6B).
FIGURE 6. Analysis of downregulated lncRNAs-mRNA co-expression. (A) The top 10 downregulated of regulatory network. Blue nodes represent lncRNAs and yellow nodes represent target genes. The pink lines represent positive correlations and the green lines represent negative correlations. (B) Top five downregulated lncRNAs co-regulated mRNA.
Differential Expression of mRNAs
Microarray data analysis also showed the overexpressed genes. Total of 115 mRNA sequences were dysregulated between the POCD and control groups (FC > 2, P < 0.05). Among them, 39 mRNAs were upregulated and 76 were downregulated (P < 0.05). Variations in mRNA expression levels of the POCD and non-POCD group are illustrated in Figure 7.
FIGURE 7. Differential expression of mRNAs. Heatmaps of mRNAs in POCD patients vs. non-POCD patients. (A) The top 10 of upregulated mRNAs. (B) The top 10 of downregulated mRNAs.
Postoperative cognitive dysfunction is a postoperative complication with high incidence and poor prognosis (van Dijk et al., 2007). At present, the mechanism of POCD is focused on inflammation, but there is little report on lncRNA level while the specific molecular mechanism is not completely understood (Terrando et al., 2011; Zhang X. et al., 2016).
As a class of important ncRNA, lncRNA plays vital role in central system (Kong et al., 2017; Liu et al., 2018). In recent years, several reports suggested that lncRNA is involved in neural physiology and diseases like Alzheimer’s disease (AD) and Parkinson’s disease (PD; Wang et al., 2018). Side by side with the identification of numerous biomarkers linked with the POCD, for instance, Aβ and tau (Wu et al., 2018), pro-inflammatory protein S100A8 (Lu et al., 2015), and preoperative levels of immunoglobulin M (Mathew et al., 2003). Along with the previously reported, NF-kappa B, IL-1b, and TNF-α (Renard et al., 1997; Chandel et al., 2000). Meanwhile, the growing number of evidences have demonstrated the crucial role of lncRNAs in neurotransmission, memory structure, and synaptic plasticity in the brain (Wang et al., 2018), such as the newly discovered lncRNA, MALAT1 (Shang et al., 2018), NEAT1 (Yan et al., 2018), and H19 (Zhao et al., 2017). In a similar trend, lncRNAs may present biomarker for molecular diagnosis of diseases. Recent evidence has shown that lncRNAs exhibit high percentage of the total transcriptome (Guennewig and Cooper, 2014; Briggs et al., 2015; Chen et al., 2016) and strongly related to neural differentiation and synaptic plasticity (Briggs et al., 2015; Gudenas and Wang, 2015). In the meantime, many researchers have reported considerable progress toward thorough understanding of the relationship between lncRNAs and brain physiology and diseases (Gudenas and Wang, 2015; Chen et al., 2016).
In the present study, various neuropsychological tests were adopted to evaluate cognitive functions 1 day before and 30 days after knee or hip replacement surgery to distinguish between POCD and non-POCD patients. Then, lncRNA sequences were studied from all the patients using developed microarray chips. Microarray analyses revealed that there were total of 29 upregulated lncRNAs and 39 mRNAs, 39 downregulated lncRNAs, and 76 mRNAs associated with POCD. Likewise, the GO and KEGG were performed to predict the roles of the dysregulated lncRNAs, and the co-expression mRNA. Furthermore, network interaction of lncRNAs and mRNAs in the POCD patients was constructed to explore the relationship between lncRNAs and mRNAs with the disease. As POCD symptoms are memory decline, poor attention, learning disability, and deficit function execution. Therefore, this experiment uses a series of neuropsychology tests to assess these five categories of cognition functions to distinguish whether cognitive impairment existed before surgery it has occurred after the surgery. To increase the accuracy, all cognitive functions were measured by more than one test to diagnose the POCD patients.
Wei et al. (2017) were the first to report lncRNA linkage to POCD in animal model. Whereas, our study has discovered various lncRNAs and mRNAs associated with POCD human patients and have expanded the understanding of the role of non-coding RNA in POCD. Microarray was used to investigate differentially expressed lncRNAs and mRNAs which then become the main data for the subsequent bioinformatics analyses of the study. This working strategy has been applied and validate in studying many diseases including cancer (Han et al., 2018; Liang et al., 2018; Lou et al., 2018), cerebral ischemia-induced lung injury (Hu et al., 2017), and traumatic brain injury (Lyu et al., 2017). In our study, microarray was used to study the changes in expression of lncRNAs and mRNAs in POCD group. The results presented 68 differentially expressed lncRNAs and 115 mRNAs in POCD group compared to non-POCD group. lnc-FAM53B-2_1 and NR_024187 were the most obvious expressions of upregulated and downregulated lncRNA, respectively. This has confirmed the previous animal model study in showing the correlation of POCD with changing levels of lncRNAs and mRNAs expression. Microarray analysis combining with GO and KEGG analyses provide reliable platforms for subsequent study for further understanding for the occurrence and progression of POCD on molecular levels.
Bioinformatics analyses were aimed to indicate the potential functions of the dysregulated lncRNAs that have been identified in our study. GO analysis has illuminated the functions of genes through biological process, cellular component, and molecular function. While KEGG is normally used to study the possible roles of expressed RNA with correlation to cell signaling pathways (Hu et al., 2017).
In our analysis of biological process, it was noticed that the differentially expressed lncRNAs, mainly participated in process of protein and histone deacetylation, cytoplasmic mRNA processing body assembly, regulation of cell cycle, and stem cell population. Previous study has also pointed out the role of histone acetylation/deacetylation in neurodegenerative conditions (Verma et al., 2017). In concur with other recent works which presented evidences of histone deacetylase-2 as player in stress-induced cognitive impairment via histone deacetylation and PI3K/AKT pathway modulation (Verma et al., 2017; Wu et al., 2017), our results are consent with these published studies. However, the recognized biological processes leading to the pathogenesis of POCD, such as inflammation, stress, and apoptosis (Tian et al., 2015; Johansson et al., 2017; Zhang et al., 2018), are less prominent in our results. Sirtuin1 (SIRT1) is a deacetylase protein that reported to mediate hippocampal neuronal apoptosis in mice and causing cognitive deficit (Min et al., 2018). Our results have presented histone deacetylase’s critical role in the pathogenesis of POCD. In addition, molecular function analysis has shown that protein deacetylase activity, poly (A) RNA binding, interleukin-12 receptor binding, and serine/threonine kinase were also involved in the POCD. Histone deacetylase is translated from HDAC gene, which is regulated by lncRNAs that was identified as differentially transcripted RNA in our study, indicating that lncRNAs exert their functions through regulating histone deacetylase activity in POCD pathogenesis. Yi et al. (2015) reported the close relatedness of inflammatory response with the changes in cognitive function in cardiopulmonary bypass patients (Yi et al., 2015). In a similar manner, interleukin-12 and its receptor were demonstrated to influence cognitive function (Johansson et al., 2017; Vonder Haar et al., 2017; Wu et al., 2017). Likewise, cytokines and cytokines receptor participate apoptosis and immune regulation through Jak-STAT signaling pathway by specifically activate STAT4 (Darnell et al., 1994; Kang and Kang, 2008). In regard to the observation of the role serine/threonine kinase in the POCD, similar observation was reported diabetic related cognitive impairment. As the expression of phosphorylated serine/threonine kinase (p-AKT) increased under in those patients, reflecting the associated of serine/threonine kinase in cognitive impair function (Zhou et al., 2017).
Furthermore, pathway analysis revealed that there were 10 pathways implicated in POCD patients, Epstein–Barr virus infection, and pantothenate and CoA biosynthesis were the highly enriched pathways in POCD patients. Epstein–Barr virus, a human gamma-herpesvirus, potentially facilitates the damage of autoimmunity and CNS tissue, causing psychiatric symptoms and cognitive dysfunction (Wynne et al., 2010). There are evidences to demonstrate critical function of pantothenate kinase 2 (pank2) and CoA homeostasis in neuronal development and functioning in zebra fish (Zizioli et al., 2016). Our data have revealed the elevated level of enrichment of pantothenate and CoA biosynthesis in POCD group, suggesting the linkage of this pathway to POCD pathogenesis. The other moderately enriched-pathways like Jak-STAT signaling pathway, MAPK signaling pathway, and metabolic pathways in our study were also recorded to be involved in occurrence and development of POCD in rats’ model (Wei et al., 2017). Inflammatory reaction can trigger several metabolic pathways, which are closely related to the changes in cognitive function in patients undergoing cardiopulmonary bypass (Yi et al., 2015). It has been confirmed the phosphoinositide 3-kinase (PI3K)/AKT/glycogen synthase kinase-3β (GSK-3β) pathway are induced due to isoflurane anesthesia induced-inflammation, stress, and apoptosis, thus resulting to cognitive impairment (Zhang Y. et al., 2016). These studies combined with ours can provide a comprehensive insight of the mechanism of POCD.
Co-expression of lncRNA and mRNA
Until the present, the functions of lncRNAs are not entirely clear. The co-expression network of aberrantly expressed mRNAs and lncRNAs in POCD was constructed to further understand of the role of lncRNAs via calculation and inference.
In the present study, we have selected the top 10 up and downregulated lncRNAs to construct the lncRNA-mRNA co-expression network of interactions to predict the genes targeted by lncRNA. We found that the increased lncRNA, lnc-FAM53B-2_1, correlated with that of SMARCE1, ATXN2, SMC4, and IL23R, while main reduced lncRNA was NR_024187 and has influence on NEK3, XKRX, VNN1, TBXAS1, and HDAC5. Previous study showed SMARCE1 to be associated with the occurrence and development of non-syndromic intellectual disability (Santen et al., 2013). Increasing number of evidences found that ATXN2 gene plays key role in neural diseases like dementia and motor neuron disease which could affect cognitive function (van Blitterswijk et al., 2014). HDAC5 has been confirmed to be involved in cognitive function including learning and memory process (Taniguchi et al., 2017). Furthermore, we discovered that the top five upregulated lncRNAs commonly regulating eight mRNAs, two among them were the mainly affected mRNA sequences. These implied that POCD is a multigene induced disease and lncRNAs may participant in POCD via regulating mRNA expression. These genes are reported in POCD for the first time assuming they will serve as the base for further research in the future.
However, it will be worthwhile to evaluated neurological function in long-term manner in POCD patients. In addition, more hands on researches should be conducted to investigate the exact mechanisms and function of genes that have been reported in this study to examine the specific interaction mechanisms between lncRNAs and POCD.
We have recorded that lncRNAs and mRNAs in serum were differentially expressed in patients with POCD compared with patients without POCD after hip or knee replacement surgery. These vital findings may contribute to well understand the mechanism of POCD and provide new biomarkers for diagnosis and monitoring the progression of POCD and potential therapeutic target for POCD treatment.
F-XY and YZ conceived and designed the experiments. YZ, Y-XL, Q-XX, and QL performed the experiments. Y-XL, RD, ID, and JB analyzed the data. YZ, MA-H, and Y-XL wrote the manuscript. MA-H and ID participated in improving the paper and making the language more authentic during the revision.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We greatly appreciate professor Ting-Hua Wang for his faithful contribution in modifying the paper.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2018.00694/full#supplementary-material
Aprahamian, I., Martinelli, J. E., Neri, A. L., and Yassuda, M. S. (2010). The accuracy of the Clock Drawing Test compared to that of standard screening tests for Alzheimer’s disease: results from a study of Brazilian elderly with heterogeneous educational backgrounds. Int. Psychogeriatr. 22, 64–71. doi: 10.1017/S1041610209991141
Briggs, J. A., Wolvetang, E. J., Mattick, J. S., Rinn, J. L., and Barry, G. (2015). Mechanisms of long non-coding RNAs in mammalian nervous system development, plasticity, disease, and evolution. Neuron 88, 861–877. doi: 10.1016/j.neuron.2015.09.045
Chandel, N. S., Trzyna, W. C., Mcclintock, D. S., and Schumacker, P. T. (2000). Role of oxidants in NF-kappa B activation and TNF-alpha gene transcription induced by hypoxia and endotoxin. J. Immunol. 165, 1013–1021.
Chen, R., Liu, L., Xiao, M., Wang, F., and Lin, X. (2016). Microarray expression profile analysis of long noncoding RNAs in premature brain injury: a novel point of view. Neuroscience 319, 123–133. doi: 10.1016/j.neuroscience.2016.01.033
Duan, X., Zhu, T., Chen, C., Zhang, G., Zhang, J., Wang, L., et al. (2018). Serum glial cell line-derived neurotrophic factor levels and postoperative cognitive dysfunction after surgery for rheumatic heart disease. J. Thorac. Cardiovasc. Surg. 155, 958.e1–965.e1. doi: 10.1016/j.jtcvs.2017.07.073
Gudenas, B. L., and Wang, L. (2015). Gene coexpression networks in human brain developmental transcriptomes implicate the association of long noncoding RNAs with intellectual disability. Bioinform. Biol. Insights 9, 21–27. doi: 10.4137/BBI.S29435
Han, C. L., Ge, M., Liu, Y. P., Zhao, X. M., Wang, K. L., Chen, N., et al. (2018). Long non-coding RNA H19 contributes to apoptosis of hippocampal neurons by inhibiting let-7b in a rat model of temporal lobe epilepsy. Cell Death Dis. 9:617. doi: 10.1038/s41419-018-0496-y
Hu, N., Guo, D., Wang, H., Xie, K., Wang, C., Li, Y., et al. (2014). Involvement of the blood-brain barrier opening in cognitive decline in aged rats following orthopedic surgery and high concentration of sevoflurane inhalation. Brain Res. 1551, 13–24. doi: 10.1016/j.brainres.2014.01.015
Hu, Y., Xiong, L. L., Zhang, P., and Wang, T. H. (2017). Microarray expression profiles of genes in lung tissues of rats subjected to focal cerebral ischemia-induced lung injury following bone marrow-derived mesenchymal stem cell transplantation. Int. J. Mol. Med. 39, 57–70. doi: 10.3892/ijmm.2016.2819
Johansson, P., Almqvist, E. G., Wallin, A., Johansson, J. O., Andreasson, U., Blennow, K., et al. (2017). Reduced cerebrospinal fluid concentration of interleukin-12/23 subunit p40 in patients with cognitive impairment. PLoS One 12:e0176760. doi: 10.1371/journal.pone.0176760
Kang, M. K., and Kang, S. K. (2008). Interleukin-6 induces proliferation in adult spinal cord-derived neural progenitors via the JAK2/STAT3 pathway with EGF-induced MAPK phosphorylation. Cell Prolif. 41, 377–392. doi: 10.1111/j.1365-2184.2008.00537.x
Kong, Z. H., Chen, X., Hua, H. P., Liang, L., and Liu, L. J. (2017). The oral pretreatment of glycyrrhizin prevents surgery-induced cognitive impairment in aged mice by reducing neuroinflammation and Alzheimer’s-Related pathology via HMGB1 inhibition. J. Mol. Neurosci. 63, 385–395. doi: 10.1007/s12031-017-0989-7
Krenk, L., Rasmussen, L. S., and Kehlet, H. (2010). New insights into the pathophysiology of postoperative cognitive dysfunction. Acta Anaesthesiol. Scand. 54, 951–956. doi: 10.1111/j.1399-6576.2010.02268.x
Liang, Y., Song, X., Li, Y., Sang, Y., Zhang, N., Zhang, H., et al. (2018). A novel long non-coding RNA-PRLB acts as a tumor promoter through regulating miR-4766-5p/SIRT1 axis in breast cancer. Cell Death Dis. 9:563. doi: 10.1038/s41419-018-0582-1
Liu, X., Yu, Y., and Zhu, S. (2018). Inflammatory markers in postoperative delirium (POD) and cognitive dysfunction (POCD): a meta-analysis of observational studies. PLoS One 13:e0195659. doi: 10.1371/journal.pone.0195659
Lou, Y., Jiang, H., Cui, Z., Wang, X., Wang, L., and Han, Y. (2018). Gene microarray analysis of lncRNA and mRNA expression profiles in patients with highgrade ovarian serous cancer. Int. J. Mol. Med. 42, 91–104. doi: 10.3892/ijmm.2018.3588
Lu, S. M., Yu, C. J., Liu, Y. H., Dong, H. Q., Zhang, X., Zhang, S. S., et al. (2015). S100A8 contributes to postoperative cognitive dysfunction in mice undergoing tibial fracture surgery by activating the TLR4/MyD88 pathway. Brain Behav. Immun. 44, 221–234. doi: 10.1016/j.bbi.2014.10.011
Lyu, Q., Zhang, Z. B., Fu, S. J., Xiong, L. L., Liu, J., and Wang, T. H. (2017). Microarray expression profile of lncRNAs and mRNAs in rats with traumatic brain injury after A2B5 + Cell transplantation. Cell Transplant. 26, 1622–1635. doi: 10.1177/0963689717723014
Mathew, J. P., Grocott, H. P., Phillips-Bute, B., Stafford-Smith, M., Laskowitz, D. T., Rossignol, D., et al. (2003). Lower endotoxin immunity predicts increased cognitive dysfunction in elderly patients after cardiac surgery. Stroke 34, 508–513.
Min, S. W., Sohn, P. D., Li, Y., Devidze, N., Johnson, J. R., Krogan, N. J., et al. (2018). SIRT1 Deacetylates Tau and reduces pathogenic Tau spread in a mouse model of tauopathy. J. Neurosci. 38, 3680–3688. doi: 10.1523/JNEUROSCI.2369-17.2018
Moller, J. T., Cluitmans, P., Rasmussen, L. S., Houx, P., Rasmussen, H., Canet, J., et al. (1998). Long-term postoperative cognitive dysfunction in the elderly ISPOCD1 study. ISPOCD investigators. International study of post-operative cognitive dysfunction. Lancet 351, 857–861.
Muret, K., Klopp, C., Wucher, V., Esquerre, D., Legeai, F., Lecerf, F., et al. (2017). Long noncoding RNA repertoire in chicken liver and adipose tissue. Genet. Select. Evol. 49:6. doi: 10.1186/s12711-016-0275-0
Planel, E., Richter, K. E., Nolan, C. E., Finley, J. E., Liu, L., Wen, Y., et al. (2007). Anesthesia leads to tau hyperphosphorylation through inhibition of phosphatase activity by hypothermia. J. Neurosci. 27, 3090–3097.
Renard, P., Zachary, M. D., Bougelet, C., Mirault, M. E., Haegeman, G., Remacle, J., et al. (1997). Effects of antioxidant enzyme modulations on interleukin-1-induced nuclear factor kappa B activation. Biochem. Pharmacol. 53, 149–160.
Santen, G. W., Aten, E., Vulto-Van Silfhout, A. T., Pottinger, C., Van Bon, B. W., Van Minderhout, I. J., et al. (2013). Coffin-Siris syndrome and the BAF complex: genotype-phenotype study in 63 patients. Hum. Mutat. 34, 1519–1528. doi: 10.1002/humu.22394
Shang, J. L., Cheng, Q., Duan, S. J., Li, L., and Jia, L. Y. (2018). Cognitive improvement following ischemia/reperfusion injury induced by voluntary runningwheel exercise is associated with LncMALAT1mediated apoptosis inhibition. Int. J. Mol. Med. 41, 2715–2723. doi: 10.3892/ijmm.2018.3484
Silbert, B. S., Evered, L. A., and Scott, D. A. (2014). Incidence of postoperative cognitive dysfunction after general or spinal anaesthesia for extracorporeal shock wave lithotripsy. Br. J. Anaesth. 113, 784–791. doi: 10.1093/bja/aeu163
Taniguchi, M., Carreira, M. B., Cooper, Y. A., Bobadilla, A. C., Heinsbroek, J. A., Koike, N., et al. (2017). HDAC5 and Its Target Gene, Npas4, function in the nucleus accumbens to regulate cocaine-conditioned behaviors. Neuron 96, 130.e6–144.e6. doi: 10.1016/j.neuron.2017.09.015
Terrando, N., Eriksson, L. I., Ryu, J. K., Yang, T., Monaco, C., Feldmann, M., et al. (2011). Resolving postoperative neuroinflammation and cognitive decline. Ann. Neurol. 70, 986–995. doi: 10.1002/ana.22664
Tian, A., Ma, H., Zhang, R., Tan, W., Wang, X., Wu, B., et al. (2015). Interleukin17A promotes postoperative cognitive dysfunction by triggering β-amyloid accumulation via the transforming growth factor-β (TGFβ)/smad signaling pathway. PLoS One 10:e0141596. doi: 10.1371/journal.pone.0141596
Vacas, S., Degos, V., Tracey, K. J., and Maze, M. (2014). High-mobility group box 1 protein initiates postoperative cognitive decline by engaging bone marrow-derived macrophages. Anesthesiology 120, 1160–1167. doi: 10.1097/ALN.0000000000000045
Valentin, L. S., Pereira, V. F., Pietrobon, R. S., Schmidt, A. P., Oses, J. P., Portela, L. V., et al. (2016). Effects of single low dose of dexamethasone before noncardiac and nonneurologic surgery and general anesthesia on postoperative cognitive Dysfunction-A Phase III double blind. randomized clinical trial. PLoS One 11:e0152308. doi: 10.1371/journal.pone.0152308
van Blitterswijk, M., Mullen, B., Heckman, M. G., Baker, M. C., Dejesus-Hernandez, M., Brown, P. H., et al. (2014). Ataxin-2 as potential disease modifier in C9ORF72 expansion carriers. Neurobiol. Aging 35, 2421.e13–2427.e17. doi: 10.1016/j.neurobiolaging.2014.04.016
van Harten, A. E., Scheeren, T. W., and Absalom, A. R. (2012). A review of postoperative cognitive dysfunction and neuroinflammation associated with cardiac surgery and anaesthesia. Anaesthesia 67, 280–293.
Verma, T., Mallik, S. B., Ramalingayya, G. V., Nayak, P. G., Kishore, A., Pai, K. S. R., et al. (2017). Sodium valproate enhances doxorubicin-induced cognitive dysfunction in Wistar rats. Biomed. Pharmacother. 96, 736–741. doi: 10.1016/j.biopha.2017.09.150
Vonder Haar, C., Martens, K. M., Riparip, L. K., Rosi, S., Wellington, C. L., and Winstanley, C. A. (2017). Frontal traumatic brain injury increases impulsive decision making in rats: a potential role for the inflammatory cytokine Interleukin-12. J. Neurotrauma 34, 2790–2800. doi: 10.1089/neu.2016.4813
Wang, D. Q., Fu, P., Yao, C., Zhu, L. S., Hou, T. Y., Chen, J. G., et al. (2018). Long Non-coding RNAs, novel culprits, or bodyguards in neurodegenerative diseases. Mol. Ther. Nucleic Acids 10, 269–276. doi: 10.1016/j.omtn.2017.12.011
Wei, C., Luo, T., Zou, S., Zhou, X., Shen, W., Ji, X., et al. (2017). Differentially expressed lncRNAs and miRNAs with associated ceRNA networks in aged mice with postoperative cognitive dysfunction. Oncotarget 8, 55901–55914. doi: 10.18632/oncotarget.18362
Witlox, J., Eurelings, L. S., De Jonghe, J. F., Kalisvaart, K. J., Eikelenboom, P., and Van Gool, W. A. (2010). Delirium in elderly patients and the risk of postdischarge mortality, institutionalization, and dementia: a meta-analysis. JAMA 304, 443–451. doi: 10.1001/jama.2010.1013
Wu, J., Liu, C., Zhang, L., Qu, C. H., Sui, X. L., Zhu, H., et al. (2017). Histone deacetylase-2 is involved in stress-induced cognitive impairment via histone deacetylation and PI3K/AKT signaling pathway modification. Mol. Med. Rep. 16, 1846–1854. doi: 10.3892/mmr.2017.6840
Wu, Z., Zhang, M., Zhang, Z., Dong, W., Wang, Q., and Ren, J. (2018). Ratio of beta-amyloid protein (Abeta) and Tau predicts the postoperative cognitive dysfunction on patients undergoing total hip/knee replacement surgery. Exp. Ther. Med. 15, 878–884. doi: 10.3892/etm.2017.5480
Wynne, A., Kanwar, R. K., Khanna, R., and Kanwar, J. R. (2010). Recent advances on the possible neuroprotective activities of Epstein-Barr virus oncogene BARF1 protein in chronic inflammatory disorders of central nervous system. Curr. Neuropharmacol. 8, 268–275. doi: 10.2174/157015910792246191
Xu, J., Zhang, R., Zuo, P., Yang, N., Ji, C., Liu, W., et al. (2012). Aggravation effect of isoflurane on Abeta(25-35)-induced apoptosis and tau hyperphosphorylation in PC12 cells. Cell Mol. Neurobiol. 32, 1343–1351. doi: 10.1007/s10571-012-9860-0
Yan, W., Chen, Z. Y., Chen, J. Q., and Chen, H. M. (2018). LncRNA NEAT1 promotes autophagy in MPTP-induced Parkinson’s disease through stabilizing PINK1 protein. Biochem. Biophys. Res. Commun. 496, 1019–1024. doi: 10.1016/j.bbrc.2017.12.149
Yi, S. Q., Yang, M., and Duan, K. M. (2015). Immune-Mediated metabolic kynurenine pathways are involved in the postoperative cognitive dysfunction after cardiopulmonary bypass. Thorac. Cardiovasc. Surg. 63, 618–623. doi: 10.1055/s-0034-1393704
Zhang, Q., Li, Y., Bao, Y., Yin, C., Xin, X., Guo, Y., et al. (2018). Pretreatment with nimodipine reduces incidence of POCD by decreasing calcineurin mediated hippocampal neuroapoptosis in aged rats. BMC Anesthesiol. 18:42. doi: 10.1186/s12871-018-0501-0
Zhang, X., Dong, H., Li, N., Zhang, S., Sun, J., Zhang, S., et al. (2016). Activated brain mast cells contribute to postoperative cognitive dysfunction by evoking microglia activation and neuronal apoptosis. J. Neuroinflammation 13:127. doi: 10.1186/s12974-016-0592-9
Zhang, Y., Zhang, Z., Wang, H., Cai, N., Zhou, S., Zhao, Y., et al. (2016). Neuroprotective effect of ginsenoside Rg1 prevents cognitive impairment induced by isoflurane anesthesia in aged rats via antioxidant, anti-inflammatory and anti-apoptotic effects mediated by the PI3K/AKT/GSK-3beta pathway. Mol. Med. Rep. 14, 2778–2784. doi: 10.3892/mmr.2016.5556
Zhao, F., Qu, Y., Liu, J., Liu, H., Zhang, L., Feng, Y., et al. (2015). Microarray profiling and co-expression network analysis of LncRNAs and mRNAs in neonatal rats following hypoxic-ischemic brain damage. Sci. Rep. 5:13850. doi: 10.1038/srep13850
Zhao, Y. H., Ji, T. F., Luo, Q., and Yu, J. L. (2017). Long non-coding RNA H19 induces hippocampal neuronal apoptosis via Wnt signaling in a streptozotocin-induced rat model of diabetes mellitus. Oncotarget 8, 64827–64839. doi: 10.18632/oncotarget.17472
Zhou, X., Wang, S., Ding, X., Qin, L., Mao, Y., Chen, L., et al. (2017). Zeaxanthin improves diabetes-induced cognitive deficit in rats through activiting PI3K/AKT signaling pathway. Brain Res. Bull. 132, 190–198. doi: 10.1016/j.brainresbull.2017.06.001
Zizioli, D., Tiso, N., Guglielmi, A., Saraceno, C., Busolin, G., Giuliani, R., et al. (2016). Knock-down of pantothenate kinase 2 severely affects the development of the nervous and vascular system in zebrafish, providing new insights into PKAN disease. Neurobiol. Dis. 85, 35–48. doi: 10.1016/j.nbd.2015.10.010
Keywords: postoperative cognitive dysfunction, lncRNA, mRNA, Gene Ontology, genetic dysregulation, microarray analysis, bioinformatics analysis
Citation: Zhang Y, Liu Y-X, Xiao Q-X, Liu Q, Deng R, Bian J, Deng IB, Al-Hawwas M and Yu F-X (2018) Microarray Expression Profiles of lncRNAs and mRNAs in Postoperative Cognitive Dysfunction. Front. Neurosci. 12:694. doi: 10.3389/fnins.2018.00694
Received: 28 June 2018; Accepted: 18 September 2018;
Published: 08 October 2018.
Edited by:Laura Ferraiuolo, The University of Sheffield, United Kingdom
Reviewed by:Varun Kesherwani, University of Nebraska Medical Center, United States
Ru-Ping Dai, Central South University, China
Copyright © 2018 Zhang, Liu, Xiao, Liu, Deng, Bian, Deng, Al-Hawwas and Yu. 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: Feng-Xu Yu, firstname.lastname@example.org
†These authors have contributed equally to this work