3.3. Transcriptomic Analyses
In an attempt to identify the molecular mechanisms that underlie the delayed behavioral effects of LPS and their modification by DEX pretreatment, transcriptomic analyses of hippocampus were performed 3 months after the drugs administration. To examine the effect of LPS, LPS-administered rats were compared with SAL group. Effect of DEX pretreatment was examined by comparing DEX+LPS and DEX+ SAL as well as DEX+LPS and LPS groups.
According to the criteria for screening of differentially expressed genes (DEGs): |log2 FC (fold change)| ≥ 1 and an adjusted P value (padj) < 0.05, 33 DEGs, 29 up-regulated and 4 down-regulated (
Alox5, LOC500035, Tbx19, H1f4), were revealed in LPS vs. SAL (
Supplementary Materials,
Table S1). Pretreatment of LPS rats with DEX increased the total number of DEGs that were 161 (128 up-regulated and 33 down-regulated) in DEX+LPS vs. DEX+SAL (
Supplementary Materials,
Table S2), whereas no DEGs were detected in DEX+ SAL vs. SAL. Comparison of DEX+LPS and LPS groups revealed 112 DEGs (110 up-regulated and 2 down-regulated). All DEGs of three compared pairs of groups, LPS and SAL, DEX+LPS and DEX+SAL, DEX+LPS and LPS, at the significant level padj < 0.05, but │log2 FC (fold change)│ ≥ 0, presented in volcano plots (
Figure 4). According to these criteria, only 3 DEGs,
Ccdc88a,
Ptprb and
RT1-N2, were found in DEX+ SAL vs. SAL (volcano plot not shown).
Using the Venn diagrams (
http://bioin forma tics.psb.ugent.be/webtools/Venn/), 22 common genes were revealed between DEGs among 33 DEGs in LPS vs. SAL and 161 DEGs in DEX+LPS vs. DEX+SAL, as well as 84 common genes among 161 DEGs in DEX+LPS vs. DEX+SAL and 112 DEGs in DEX+LPS vs. LPS (
Figure 5a;
Supplementary Materials,
Table S3).
Figure 5b shows expression profiles of 33 DEGs for LPS vs. SAL. Data were expressed as log2 transformation of fold changes. Asterisks indicate 22 DEGs shared between LPS vs. SAL and DEX+LPS vs. DEX+SAL.
Gene ontology (GO) analysis showed that DEGs after LPS administration (LPS vs. SAL) were predominantly associated with immune and inflammatory responses (
Supplementary Materials,
Table S4; p-value < 0.05). The 13 most enriched GO terms for biological processes (FDR (false discovery rate) < 0.05) were “antigen processing and presentation of exogenous peptide antigen via MHC class II” (GO:0019886), “antigen processing and presentation of peptide or polysaccharide antigen via MHC class II” (GO:0002504), “response to interferon-gamma“ (GO:0034341), “peptide antigen assembly with MHC class II protein complex” (GO:0002503), “positive regulation of immune response” (GO:0050778), “positive regulation of T cell activation” (GO:0050870), “adaptive immune response” (GO:0002250), “inflammatory response” (GO:0006954), “antigen processing and presentation” (GO:0019882), “immune response” (GO:0006955), “positive regulation of neutrophil chemotaxis” (GO:0090023), “immunoglobulin mediated immune response” (GO:0016064), “humoral immune response” (GO:0006959). Genes that were linked to these biological process terms included the major histocompatibility complex (MHC) class II members (
Cd74,
RT1-Ba, RT1-Bb,
RT1-Da,
RT1-Db1,
RT1-Db2) and their expression transactivator (
Ciita), adhesion receptor (
Adgre1), chemoattractants for B lymphocyte (
Cxcl13), complement C3 (
C3), neutrophil cytosolic factor 1 (
Ncf1), proinflammatory cytokines (
Il1b, Tnfsf13b), regulators of cellular responses to cytokines (
Csf2rb) and inflammatory response to IFNG/IFN-gamma (
Mefv), catalyzer of the first step in leukotriene biosynthesis (
Alox5), Ras-related C3 botulinum toxin substrate 2 (
Rac2), and ubiquitin-like protein modifier (
Ubd).
Contrary to our expectation, pretreatment with DEX did not significantly affect the responses of most of these immune/inflammatory-related genes to LPS.
Figure 6 shows similar changes in expression of several representative genes such as MHCII genes and
Ciita (a),
Il1b and
Tnfsf13b (b) after LPS alone and with DEX pretreatment. The RNA-sec results were confirmed by quantitative PCR for the
IL1b chosen for verification (c). Significant correlation (r = +0.93) was found between sec and PCR
IL1b expression values.
Due to this, GO analysis (DAVID) of DEGs that were revealed between DEX+LPS and DEX+SAL also showed that the top 8 biological processes (FDR < 0.05) mainly related to immune response (
Supplementary Materials,
Table S5): “antigen processing and presentation of exogenous peptide antigen via MHC class II” (GO:0019886), “immunoglobulin production involved in immunoglobulin mediated immune response” (GO:0002381), “peptide antigen assembly with MHC class II protein complex” (GO:0002503), “antigen processing and presentation of peptide or polysaccharide antigen via MHC class II” (GO:0002504), “positive regulation of T cell activation” (GO:0050870), “antigen processing and presentation” (GO:0019882), “response to interferon-gamma”(GO:0034341), “adaptive immune response” (GO:0002250).
To explore the candidate genes and pathways that might play important roles in preventing LPS-induced memory impairment by DEX, STRING (the Search Tool for the Retrieval of Interacting Genes/Proteins,
https://string-db.org/) analyses was performed for common and unique DEGs between LPS vs. SAL and DEX+LPS vs. DEX+SAL.
Among the 33 DEGs from LPS vs. SAL, 3 genes (Cd4, Cd74, Il1b) ranked with the interaction degree ≥ 10 connections/interactions, were served as the top hub genes (Table). Among DEGs from comparison of DEX+LPS and DEX+SAL, the number of top hub genes was 6, the list of which in addition to similar with LPS vs. SAL-related DEGs (Cd4, Cd74, Il1b), included Ciita, Creb1 and RT1-Bb (Table). The high level of connections (n = 9) was also observed for Grin2a gene. Of these hub genes, Cd4 exhibited the highest degree.
Of the 33 DEGs revealed between LPS and SAL, 9 did not fall into the Protein-protein interaction (PPI) network (the number of connections/interactions = 0), which, when made from only other DEGs, was constructed of 24 nodes and 62 edges. K-means clustering of these DEGs outlined two clusters depicted in
Figure 7a, the biological functions of each of that were characterized by the Gene Ontology (GO). The cluster 1 (green) included 17 DEGs together with
Cd4 and
Cd74 hub genes (
Cd22,
Cd4,
Cd74,
Ciita,
Csf2rb,
Cxcl13,
Emr1,
Itgal,
Mx1,
Plac8,
RT1-Ba,
RT1-Bb,
RT1-Da,
RT1-Db1,
RT1-Db2,
Rac2,
Tnfsf13b). These DEGs were mainly enriched in immune response and immune system processes. The second cluster consisted of 7 DEGs together with
Il1b hub gene (
Alox5,
C3,
F10,
Ifnlr1,
Il1b,
Mefv,
Ncf1), which were involved in inflammatory response and response to stress.
From 22 DEGs shared between LPS vs. SAL and DEX+LPS vs. DEX+SAL, 16 DEGs together with
Cd4,
Cd74 and
Il1b hub genes (
Alox5,
C3,
Cd22,
Cd4,
Cd74,
Ciita,
Cxcl13,
Emr1,
Il1b,
Itgal,
Ncf1,
RT1-Ba,
RT1-Bb,
RT1-Da,
RT1-Db1,
Tnfsf13b) that felled into the PPI network, were also predominantly related to immune and inflammatory responses (
Figure 7b).
In addition to similar with LPS vs. SAL GO biological processes associated with immune-inflammatory activation that were revealed by DAVID analyses, DEGs from DEX+LPS vs. DEX+SAL also enriched (p-value < 0.05) in “regulation of ion transmembrane transport” (GO:0034765), “cellular response to growth factor stimulus”(GO:0071363), “regulation of membrane potential” (GO:0042391), “regulation of long-term neuronal synaptic plasticity” (GO:0048169), “potassium ion transmembrane transport” (GO:0071805), “intracellular signal transduction” (GO:0035556), “calcium ion transport” (GO:0006816) (
Supplementary Materials,
Table S5). Genes associating with these terms changed their expression specifically in LPS-administered rats pretreated with DEX, but not in rats received LPS or DEX alone. The list of these genes included calcium channels (
Cacna1e,
Cacna2d1), Cyclic AMP-responsive element-binding protein 1 (
Creb1), Gamma-aminobutyric acid receptor subunits (
Gabrb2,
Gabrr2), Glutamate ionotropic receptor nmda type subunit 2a (
Grin2a) and Metabotropic glutamate receptor 5 (
Grm5), Potassium voltage-gated channels (
Kcnh5,
Kcnh7,
Kcnq3), Mitogen activated protein kinase kinase kinase 2 (
Map3k2).
Figure 8 illustrates the specific increasing effects of DEX on representative genes in LPS-treated rats.
Of 84 specific for DEX+LPS vs. DEX+SAL DEGs, 39 DEGs were felled into the string PPI network. Using k-means clustering, 21 DEGs, including
Grin2a hub gene, were grouped (
Figure 9). These genes were associated with such GO biological processes, as “Regulation of transmembrane transport” (GO:0034762), “Cognition” ( GO:0050890), “Learning” (GO:0007612), “Neurogenesis” (GO:0022008), “Nervous system development” (GO:0007399).