3.2 LIPUS and CUS altered the hippocampus and prefrontal cortex transcriptome profile
Previously published work indicated that LIPUS effectively mitigated anhedonia in CUS rats, as evidenced by increased sucrose consumption in the sucrose preference test and reduced immobility time in FST. In this study, we tested the mRNA expression profiles in the PFCs and hippocampus, two brain regions that are closely associated with depression, in control, CUS, and LIPUS-treated rats. The number of differentially expressed genes (DEGs) among the groups is shown in Table 2.
Table 2
Statistics of the differentially expressed genes
Comparisons
|
Hippocampus
|
Prefrontal cortex
|
Up
|
Down
|
Up
|
Down
|
Ctrl vs. CUS
|
1805
|
1659
|
923
|
858
|
Ctrl vs. LIPUS
|
1583
|
1304
|
1466
|
1482
|
CUS vs. LIPUS
|
897
|
700
|
1029
|
1081
|
Ctrl, control group; CUS, chronic unpredictable stress group; LIPUS, low-intensity pulsed ultrasound group.
In the hippocampus, the CUS group exhibited 3,464 DEGs, including 1,805 up-regulated and 1,659 down-regulated genes (Fig. 2A), when compared to the control group. Following ultrasound stimulation, we observed significant differential expression of 1,597 genes in the LIPUS group relative to CUS rats, comprising 897 up-regulated and 700 down-regulated genes (Fig. 2B). In the PFC, we identified 1,781 DEGs in the comparison between the CUS and control groups (Fig. 2C). Similarly, in the comparison between the LIPUS and CUS groups, 2,110 DEGs were identified (Fig. 2D). Among the DEGs, 923 genes were up-regulated and 858 were down-regulated in the CUS group, while in the LIPUS group, 1,029 genes were up-regulated and 1,081 were down-regulated. Venn diagrams show the overlaps of regulated genes (Fig. 2E, F). Notably, in the hippocampus, ultrasound stimulation reversed the expression of 592 genes, with 359 being up-regulated and 233 being down-regulated in the LIPUS groups. Similarly, in the PFC, ultrasound stimulation reversed the expression of 254 genes in the CUS group, of which 118 were up-regulated and 136 were down-regulated (Table 3).
Table 3
Statistics of the gene changes induced by LIPUS, compared with CUS
Position
|
Down in CUS but up in LIPUS
|
Up in CUS but down in LIPUS
|
Hippocampus
|
359
|
233
|
Prefrontalcortex
|
136
|
118
|
Ctrl, control group; CUS, chronic unpredictable stress group; LIPUS, low-intensity pulsed ultrasound group.
3.3 Functional analysis of DEGs in LIPUS and CUS rats
To explore the functional changes induced by LIPUS, we performed pathway analysis using the identified DEGs. We utilized the Gene Ontology (GO) database, which classifies gene functions into three categories: biological process (BP), cellular component (CC), and molecular function (MF). Enriched GO terms with a p-adj < 0.05 were considered to be significantly enriched.
In the hippocampus, the down-regulated genes in the depressive model rats, compared to the control group, were primarily associated with catalytic activity, oxygen carrier activity, and cofactor binding in the MF category (Fig. 3A). Conversely, the up-regulated genes in the MF category were related to GTPase binding and protein binding. The differentially expressed mRNAs in the BP category were mainly involved in catabolic processes and regulation of neuron projection development. In terms of the CC category, the differentially expressed mRNAs were associated with nuclear speck and envelope (Fig. 3B). Following ultrasound stimulation, the LIPUS rats showed 897 up-regulated genes compared to the CUS group. In the BP category, the enriched genes were significantly associated with oxygen and gas transport. In the CC category, most of the genes were anchored components of the synaptic membrane and plasma membrane, while some were involved in oxygen carrier activity and oxygen binding in the MF category (Fig. 3D). Additionally, we analyzed the 700 down-regulated genes in the LIPUS group. According to the GO enrichment analysis, these genes were associated with various biological processes, including epithelial cell migration, ameboidal-type cell migration, tissue migration, and substrate adhesion-dependent cell spreading (Fig. 3C).
In the PFC, we analyzed the top ten enriched terms of down-regulated and up-regulated genes in the CUS and LIPUS groups (Fig. 3E, H). The down-regulated genes in the depressive model rats, compared to the control group, were primarily associated with organic anion transmembrane transporter activity and anion transmembrane transporter activity in the MF category. In terms of biological processes, these genes were involved in sterol metabolic process, sterol biosynthetic process, secondary alcohol biosynthetic process, and antigen processing and presentation in the BP category (Fig. 3E). Conversely, the up-regulated genes identified in the PFC were associated with protein folding, the establishment of protein localization to organelle, and response to the topologically incorrect protein in the BP category. In the MF category, these genes were involved in unfolded protein binding, heat shock protein binding, and in the CC category, they were associated with ribonucleoprotein granules (Fig. 3F).
After ultrasound stimulation, we observed 1,029 up-regulated genes and 1,081 down-regulated genes in the PFC. The down-regulated genes were associated with functions such as dendritic spine membrane, C2H2 zinc finger domain binding, metallopeptidase activity, cell adhesion molecular binding, and heart development, among others (Fig. 3G). On the other hand, the up-regulated genes were involved in defense response to other organisms, response to viruses, innate immune response, production of molecular mediators of the immune response, immunoglobulin production, double-stranded RNA binding, and others (Fig. 3H).
We utilized the Kyoto Encyclopedia Genes and Genomes (KEGG) database to analyze the enriched pathways with p-adj < 0.05. In the hippocampus, the CUS group exhibited several inhibited pathways compared to the control group, including amino sugar and nucleotide sugar metabolism, spliceosome, pentose phosphate pathway, and steroid biosynthesis (Fig. 4A). Conversely, several pathways were activated in the CUS group, including protein processing in endoplasmic reticulum, longevity regulating pathways (multiple species), the Forkhead box, class O signaling pathway, autophagy, longevity regulating pathway, and Rap1 signaling pathway (Fig. 4B). In the comparison between the LIPUS and CUS groups, ultrasound stimulation down-regulated several pathways in the hippocampus, including proteoglycans in cancer, transforming growth factor-𝛽 signaling pathway, mitogen-activated protein kinase (MAPK) signaling pathway, phosphoinositide 3-kinase-protein kinase B (PI3K-Akt) signaling pathway, focal adhesion, and amino sugar and nucleotide sugar metabolism (Fig. 4C). On the other hand, ultrasound stimulation activated pathways such as African trypanosomiasis, malaria, terpenoid backbone biosynthesis, leishmaniasis, circadian entrainment, and steroid biosynthesis (Fig. 4D).
In the PFC, the KEGG pathway analysis revealed that the CUS group exhibited inhibition of pathways such as steroid biosynthesis, Epstein-Barr virus infection, type1 diabetes mellitus, retinoic acid-inducible gene-1 (RIG-1)-like receptor signaling pathway, ATP binding cassette transporters, and Parkinson’s disease (Fig. 5A). The up-regulated genes in the CUS group were enriched in pathways including protein processing in the endoplasmic reticulum, antigen processing and presentation, mRNA surveillance pathway, RNA degradation, and circadian rhythm (Fig. 5B).
In comparison, ultrasound stimulation down-regulated pathways such as the MAPK signaling pathway, PI3K-Akt signaling pathway, and human papillomavirus infection in the PFC (Fig. 5C). Conversely, ultrasound stimulation activated pathways including the nucleotide oligomerization domain (NOD)-like receptor signaling pathway, Herpes simplex virus 1 infection, Epstein-Barr virus infection, Kaposi sarcoma-associated herpesvirus infection, Hepatitis C, and influenza A (Fig. 5D).
3.4 Identification of DEGs under LIPUS Treatment in the CUS model.
To confirm the results from RNA-sequencing, we performed qRT-PCR analysis on a selected set of ten RNAs. These genes were chosen based on the following criteria: 1. log2 fold change > 1 or log2 fold change < -1; 2. p-adj < 0.1; 3. reversal of gene expression by LIPUS treatment. The qRT-PCR results revealed that claudin2 (Cldn2), occludin (Ocln), growth hormone 1 (GH-1), glucokinase (GCK), and hypothetical protein LOC688459 exhibited high expression levels in the hippocampi of CUS rats, all of which were reversed after LIPUS treatment. Furthermore, LIPUS also reversed the reduced expression of hemoglobin subunit β (Hbb), ATP synthase peripheral stalk-membrane subunit b (ATP5BP), and hemoglobin subunit α (Hba) in CUS rats. No significant changes in the expression of phosphoglycerate kinase 1 (PGK1) were observed following CUS, but LIPUS down-regulated its expression (Fig. 6A).
In the PFC, the qPR-PCR analysis showed prominent reversals of gene expression by LIPUS on Hbb, 𝛽-globin, Hba, LOC688459, GH-1, PKG1, and GCK. Additionally, although the down-regulated expression of ATP5BP in the CUS group did not show a significant difference compared to the control group, LIPUS significantly improved its expression in the PFC. In contrast to the hippocampus, where Cldn2 and Ocln were down-regulated in the CUS group compared to the control group, the expression of these genes was further down-regulated by LIPUS in the PFC (Fig. 6B). Notably, LIPUS reversed the expression of Hbb, GH1, and GCK in both the hippocampus and PFC.