Identification of a lncRNA/circRNA-miRNA-mRNA ceRNA Network in Alzheimer’s Disease

Background : Alzheimer’s disease (AD) occurs in the elderly and pre-elderly, characterized by decline of memory, cognitive dysfunction, impairment of learning capacity, and motor dysfunction. Recently a competitive endogenous RNA (ceRNA) network has been found to be related to AD progression, but there is still little understanding of the ceRNA regulatory network in AD. This study aims to explore the important regulatory mechanisms of ceRNA regulatory networks containing long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in AD. Methods : Data from the gene expression omnibus (GEO) database were used for the analysis. To study enrichment function for the upregulated and downregulated mRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the Metascape database, respectively. Based on the STRING database and Cytoscape software 3.9.1, a protein-protein interaction (PPI) network was constructed. The hub genes in this network were identified utilizing the CytoHubba plugin in Cytoscape. The TargetScan, miRWalk, and miRDB were selected to calculate the regulatory interaction between miRNAs and the hub genes. LncRNAs were predicted using RNA22. Ad-ditionally, circRNA prediction was executed using the circBank database. Results : 711 downregulated and 670 upregulated overlapping mRNAs were identified between AD and control samples. 32 downregulated and 340 upregulated miRNAs were obtained from AD samples compared with control samples. 78 upregulated and 205 downregulated circRNAs were screened. 275 upregulated lncRNAs and 209 downregulated lncRNAs were found between AD samples and control samples. The PPI network constructed consists of 1016 nodes and 13,946 edges. Ten hub genes were selected to identify target miRNAs and ceRNAs. On the basis of the ceRNA hypothesis, a circRNA/lncRNA-miRNA-mRNA network was established. It included five lncRNAs (TRHDE-AS1, SNHG10, OIP5-AS, LINC00926 and LINC00662), 26 circRNAs, five miRNAs (hsa-miR-3158-3p, hsa-miR-4435, hsa-let-7d-3p, hsa-miR-330-5p and hsa-miR-3605-3p), and ten mRNAs ( RPL11 , RPL34 , RPL21 , RPL22 , RPL6 , RPL32 , RPL24 , RPL35 , RPL31 , and RPL35A ). RPL35 and RPL35A were found to be significantly associated with AD pathology in tau and A β line AD models by the AlzData database. The study discovered the significance of several lncRNA-miRNA-mRNA axes and circRNA-miRNA-mRNA axes that included RPL35A and RPL35 . Conclu-sions : ceRNAs were found to be important regulators in the development of AD and provide potential biological therapy targets for AD management.

Alzheimer's disease (AD) is considered to be a complex heterogeneous disease.Its exact pathogenesis remains unknown, and has no effective treatment method [4].Thus, there is a need for further research on underlying disease processes to better understand AD and aid the development of feasible treatment options.Considerable evidence shows that ncRNAs, particularly lncRNAs, miRNAs, and circR-NAs, are a major class of regulatory molecules involved in the pathophysiology associated with AD [5].
LncRNAs are ncRNAs with a length of more than 200 nts and have been reported as important for cell regulation.Cytoplasmic lncRNA act as miRNA sponges or miRNA precursors to alter the expression and function of miRNA.They also recognize target messenger RNA (mRNA) and interact with cellular translation mechanisms.In SH-SY5Y cells, the upregulation of nuclear enriched abundant transcript 1 (lncRNA NEAT1) has been related to the downregulation of miR-27a-3p and upregulation of amyloid (Aβ), β-amyloid-precursor-protein-cleaver-enzyme 1 (BACE1) protein, amyloid precursor protein (APP), and tau protein [6].Studies suggest that LncRNA Ribonuclease P RNA component H1 (Rpph1) participate in the regulation of the onset and progression of AD by competing with miR-330-5p and miR326-3p [7].
MiRNAs are small molecules with transcripts of 21-23 nts.Changes in the levels of several miRNAs and their targeted mRNA have been reported to affect the development and progression of the disease such as APP production, Tau aggregation, Aβ production, protein misfolding, neuroinflammation, and oxidative damage [8].For example, compared with controls, the expression of miR-485-3p was significantly increased in the brain tissues of AD patients [9], and MiR-485-3p antisense oligonucleotide (ASO) reduced the expression level of Aβ via phagocytosis of Aβ mediated by CD36 [9].Additionally, it is reported that miR-146a overexpression reduces cognitive impairment, secretion of proinflammatory cytokines, Aβ plaque accumulation and protects the neurons of APP/PS1 transgenic mice from damage [10].
CircRNAs are non-coding RNA molecules that do not have a 5 ′ end cap and a 3 ′ end poly (A) tail that form a ring structure by a covalent bond.CircRNAs have been reported to play a significant role in various biological functions as a miRNA sponge and transcription regulator.For example, circRNA_0003611 is augmented in AD patients.It was reported that silenced circRNA_0003611 might reduce neuronal damage triggered by Aβ via adjusting the miR-383-5p/kinesin family member 1B (KIF1B) axis [11].CircRNA Cwc27 is a neuronal-rich circRNA substantially expressed in the brain and significantly increases in AD patients and mice [12].CircRNA Cwc27 deficiency significantly ameliorates AD-related pathological features and increases cognitive function [12].CircRNA Cwc27 was shown to directly bind to purine-rich element-binding protein A (Purα) and repress Pur-α recruitment to AD gene promoter clusters [13].
Competitive endogenous RNA (ceRNA) includes mRNA of coding-protein lncRNA and circRNA.Numerous studies have emphasized that these ceRNA molecules regulate miRNA target gene expression via competitively binding to the same miRNA by miRNA responsive elements (MREs).Recently, it was reported that the ceRNA network plays a critical role in many diseases, such as cardiovascular, neuroimmune and neurodegenerative disease [14][15][16].Thus, the ceRNA network offers new approaches to understanding AD.Though some ceRNAs have been found related to AD progression, there is currently little discussion of the role of the ceRNA regulatory network in AD.In this study, data were obtained from the gene expression omnibus (GEO) database.Differentially expressed circRNAs (DEcircRNAs), mRNAs (DEmRNAs), lncRNAs (DElncR-NAs), and miRNAs (DEmiRNAs) were screened from sequencing data.The targeted miRNAs and their related cir-cRNAs and lncRNAs were forecast, an integrated ceRNA network was constructed, and their potential value for AD diagnosis was explored.The key pathways and regulatory mechanisms of central genes were analyzed and discussed using bioinformatics tools.MiRNAs, targeted mRNAs, related lncRNAs, and circRNAs involved in the ceRNA network may be possible biomarkers and treatment targets for AD.The workflow is shown in Fig. 1.

Microarray Datasets
MRNA, miRNA, circRNA, and lncRNA expression datasets were obtained from the GEO database (https:// www.ncbi.nlm.nih.gov/geo/).GSE63063 was selected for analysis of differentially expressed mRNAs from AD patients compared with controls, in which 329 samples (145 AD cases and 104 controls) based on GPL6947 platform (chip A) were collected and 382 samples (139 AD cases and 134 controls) based on the GPL10558 platform (chip B) were selected.GSE120584 was downloaded to derive DEmiRNAs between AD patients and controls, in which 1309 data sets were downloaded (1021 AD cases and 288 controls) for analysis.The circRNA expression data were derived from GSE161435 (3 AD cases and 3 control cases).The lncRNAs expression data were obtained from GSE182910 (3 AD cases and 3 control cases).
Details of the microarray data selected are described in Table 1.

Identification of DEmRNAs, DEmiRNAs, DecircRNAs, and DElncRNAs
The differential expression analysis of mRNAs, miR-NAs, circRNAs, and lncRNAs among different groups was performed using GEO2R.GEO2R completed the analysis using limma R software packages (version 3.56.0)(limma R software, University of Melbourne, Parkville, Victoria, Australia).p < 0.05, false discovery rate (FDR) <0.05 and an absolute value of log 2 -fold ≥0.01 were used as thresholds to screen DEmRNAs.mRNAs without symbol annotation were excluded.When different probes fitted the same mRNA, the mRNAs with the highest absolute value of log 2fold were selected.
The criteria for the selection of DEmiRNAs were FDR <0.05, p < 0.05 and the absolute value of log 2 -fold >0.2.The miRNA with the highest absolute value of log 2 -fold was selected when different probes fitted the same miRNA.
DEcircRNAs were identified with the criterion pvalue < 0.05 and log 2 -fold >0.4 were considered significant.
DElncRNAs with p < 0.05 and log 2 -fold >1.5 were considered to indicate significant differences in expression.

PPI Network Integration of Modules and Their Analysis
The interaction relationships among DEmRNAs were constructed by STRING version 11.5 (https://cn.string-db.org/) [17].High confidence (0.7) was selected as the cutoff criterion for the minimum required interaction score.The protein-protein interaction (PPI) network of DEmR-NAs was visualized by Cytoscape 3.9.1 (Cytoscape, Institute for Systems Biology, Seattle, WA, USA) [18].Key modules of the PPI network were screened using the Molecular Complex Detection (MCODE) plugin in Cytoscape 3.9.1 (Cytoscape, Institute for Systems Biology, Seattle, WA, USA).Cutoff = 2, node score cutoff = 0.2, k-core = 2, and maximum depth = 100 were regarded as selection criteria.

Identification and Verification of Hub Genes
Hub genes of the network were identified by applying CytoHubba which is a plugin for Cytoscape 3.9.1 (Cytoscape, Institute for Systems Biology) [19].The top ten nodes ranked by the maximum group centrality (MCC) algorithm were selected as hub genes.The AlzData database, a one-stop database, covers high-throughput omics data and high-confident functional data of AD [20].The hub genes identified in this study were inserted into the AlzData database (http://www.alzdata.org/)to obtain each hub generelated cell type and correlation with AD pathology (Aβ and tau).KEGG pathway enrichment of the hub genes was analysed via Metascape database (v3.5.20230501) (http: //metascape.org/gp/index.html#/main/step1).

CeRNA Network of DEmiRNA, Hub Genes, DElncRNA, and DEcircRNA
A ceRNA network was constructed based on the regulatory mechanism among DEcircRNAs, DElncRNAs, DEmiRNAs, and DEmRNAs.
DEmiRNA binding sites of the circRNA were predicted via the circBank database (http://www.circbank.cn/)and RNA22 database was selected to calculate DEmiRNAs binding sites of the lncRNAs.Finally, a ceRNA network of hub genes, DEmiRNAs, DElncRNAs, and DEcircRNAs was constructed.

Statistical Analysis
Statistical analyses were performed using GraphPad Prism 8 (GraphPad Software, Inc., San Diego, CA, USA).A p-value < 0.05 was considered statistically significant.The differential expression of mRNAs and miRNAs was analyzed using limma R packages (version 3.56.0)(limma R software, University of Melbourne, Parkville, Victoria, Australia).A p-value < 0.05 and FDR <0.05 were considered statistically significant.Differential analyses of lncRNAs and circRNAs expression were performed using limma R packages (version 3.56.0)(limma R software, University of Melbourne).A p-value < 0.05 was considered statistically significant.

Functional Enrichment Analysis of DEmRNAs
Functional enrichment analysis showed that the upregulated DEmRNAs were mainly enriched in the regulation of response to biotic stimulus, regulation of defense response and regulation of innate immune response in the biological process (BP) category; secretory granule mem-brane, tertiary granule, specific granule and specific granule membrane in the cellular component (CC) category; GT-Pase regulator activity, nucleoside-triphosphatase regulator activity, enzyme activator activity, and small GTPase activator activity in the molecular function (MF) category (Fig. 3A,C,E and Supplementary Table 5).For the downregulated DEmRNAs, functional enrichment analysis indicated that the enriched GO terms mainly included translation, peptide biosynthetic process, and peptide metabolic process in the BP category; ribosomal subunit, focal adhesion, cell-substrate junction, and mitochondrial proteincontaining complex in the CC category; structural constituent of ribosome, rRNA binding, oxidoreduction-driven active transmembrane transporter activity and nicotinamide adenine dinucleotide (NADH) dehydrogenase (ubiquinone) activity in the MF category (Fig. 3B,D,F and Supplementary Table 6).
As shown in Fig. 3G,H and Supplementary Tables 7,8, KEGG pathways for the upregulated DEmRNAs mainly included osteoclast differentiation, NOD-like receptor signaling pathway, lipid and atherosclerosis and epstein-Barr virus infection (Supplementary Table 7), while the downregulated DEmRNAs were mainly enriched in various metabolic pathways such as ribosome, oxidative phosphorylation and Alzheimer disease (Supplementary Table 8).

Construction of PPI Interaction Network, Module Analysis, and Identification of Hub Genes
A PPI network was constructed, which consisted of 1016 nodes and 13,946 edges.The MCODE plug-in was utilized for finding key subnetworks and genes based on the relationship between edges and nodes in the network for downstream analysis.Based on MCODE analysis, 44 modules were obtained (Supplementary Table 9), 6 modules of which were selected for Fig. 4A-F.The highest score of module 1 is 53.724, which was significantly enriched in ribosomes and protein export (Fig. 4G).Module 2 was mainly related to the ribosome (Fig. 4H).Similarly, module 3 was significantly enriched in oxidative phosphorylation, and so on (Fig. 4I).Module 4 was mainly assembled in ribosome biogenesis (Fig. 4J).Module 5 was primarily associated with the proteasome (Fig. 4K).Module 6 was mainly associated with the spliceosome (Fig. 4L).
Based on these results, it was shown that these hub genes may have a significant role in AD pathogenesis.

The DEmiRNAs-Hub Gene Network
MiRNAs are endogenous non-coding RNAs with regulatory functions.The size of miRNAs is about 21~23 nts.They inhibit target mRNA expression by recognizing target mRNA through base complementary pairing.Based on this, the DEmiRNAs-hub genes pairs with an up-down mode were selected for further study.
A total of 337 overlapped upregulated DEmiRNAs related to the above ten downregulated hub genes were collected using TargetScan, miRWalk, and miRDB databases.Based on the regulatory relationship, a miRNA-mRNA network was constructed, as given in Supplementary Table 10.
To further verify the 337 overlapped upregulated DEmiRNAs, a GSE46579 dataset for studying miRNA expression profiles in the blood of AD patients and controls was collected.In total, GSE46579 detected 180 differentially expressed miRNAs, including 90 downregulated and 90 upregulated miRNAs in AD samples compared with controls.Then, the common DEmiRNAs were screened between the 337 DEmiRNAs and the 180 differentially expressed miRNAs, which revealed 5 common upregulated DEmiRNAs.Information in detail is shown in Table 4.

CeRNA Network of DEmiRNAs, Hub Genes, DElncRNAs, and DEcircRNAs
In addition to mRNA, lncRNA and circRNA also have MREs.When lncRNA or circRNA and mRNA have the     same MREs, they will competitively bind to the same miR-NAs.In other words, lncRNA or circRNA expression level in cells directly influences the expression of miRNAs that bind to the corresponding mRNA.That is, through the MREs bridge, lncRNA or circRNA indirectly regulating the mRNA expression level, thus regulate cell function.
On the basis of the aforementioned 5 overlapped upregulated DEmiRNAs, 36 DElncRNA-DEmiRNA regulatory pairs were obtained using RNA22 and 53 DEcircRNA-DEmiRNA regulatory pairs were obtained using CircBank.
Based on the ceRNA hypothesis, the DEcircRNAs/DElncRNAs-DEmiRNAs-hub genes with a mode of down-up-down were chosen for further investigation.
Generally, the lncRNAs in cytoplasm competitively bound to miRNA to regulate mRNA expression level.The RNALocate v2.0 database was selected to perform the intracellular localization of these lncRNAs (Fig. 7).It should be noted that the lncRNAs not located in cytoplasm or unclear were excluded from the network.The network of DElncRNAs-DEmiRNAs-hub genes (including 5 lncR-NAs, 3 miRNAs and 10 mRNAs) is given in Table 6 and Fig. 8. AD-associated lncRNAs were then identified from the 5 DElncRNAs using the LncRNADisease v2.0 online software.Out of the 5 DElncRNAs, 3 have previously been shown to be associated with AD (SNHG10, OIP5-AS1 and LINC00662).

Discussion
Various ceRNAs are abnormally expressed in AD [23][24][25] and mounting evidence has shown the relationships among dysregulated lncRNAs, circRNAs, miRNAs, and mRNAs and the construction of AD-related ceRNA networks [26].The ceRNA network is beneficial for a deeper understanding of the AD mechanism.In this study, transcriptome data of AD patients from the GEO database was used that contained mRNAs, circRNA, lncRNAs, and miRNAs.A total of 1381 DEmRNAs, 283 DEcircRNA, 484 DElncRNAs, and 372 DEmiRNAs were identified.A circRNA/lncRNA-miRNA-mRNA ceRNA network with biologic functions in AD was constructed following bioinformatic analysis.It is noteworthy that the ceRNA hypothesis indicates that there is a negative correlation between the ceRNA and miRNA expression levels but a positive correlation with mRNA expression.Therefore, the predicted relationships were integrated with the corresponding expression data to obtain more reliable results.
Through functional enrichment analysis, the upregulated DEmRNAs were mainly associated with regulation of response to biotic stimulus, regulation of defense response and regulation of innate immune response.Additionally, the downregulated DEmRNAs were mainly involved in translation, peptide biosynthetic process and peptide metabolic process.Importantly, ten hub genes were identified (RPL11, RPL34, RPL21, RPL22, RPL6, RPL32, RPL24, RPL35, RPL31 and RPL35A) using the CytoHubba tool.These ten hub genes were found to be downregulated in AD samples.The results of enrichment analysis showed that hub genes were mainly involved in ribosomes.These hub genes are members of a ribosomal protein family, which is important for the end and correctness of translation in the ribosome [27].Ribosomal proteins were reported to be significant parts of ribosomes and participate in regulating bioprocesses [28].Recently, it was reported that ribosomal proteins are tightly related to the progress of AD [29][30][31].
Previous investigation has shown that the low expression level of RPL11 is related to ribosome pathway dysfunction during AD development, which is the same as the conclusion of this analysis.There was evidence to suggest that RPL11 was correlated with brain maturation [32].RPL11 contributed to the cortical neuron function damaged by ribosomal stress in neurodegenerative diseases [33].Thus, it was speculated here that RPL11 may influence AD progression via the ribosome pathway.RPL34, a member of the ribosomal protein L34E family, is a highly conserved protein and has been demonstrated to participate in the occurrence and development of human malignancies.No evidence showed a role for RPL34 in AD development.The RPL21 gene was identified to be involved in Parkinson's disease [34].RPL22, a ribosomal protein, is induced by Cre recombinase, which primarily participates in the formation of the 60S subunit [35].Recent research has demonstrated that RPL22 plays a significant role in normal T cell development to arrest new protein synthesis by regulating endoplasmic reticulum stress signal transduction [36] and deactivation or deficiency of the single allele of RPL22 has been found [37].However, its association with neurodegenerative diseases such as AD has not been evaluated.RPL6 is a gene encoding a 60S subribosome, which plays an important role in oxidative phosphorylation and neuronal signal transduction in neurodegenerative diseases such as Parkinson's disease [38].RPL6 has been reported to be downregulated in Parkinson's disease and induce 5 ′ -adenosine monophosphate (AMP)-activated protein kinase signaling activation [39].Similarly, down-regulation of RPL6 was observed in Parkinson's disease peripheral blood mononuclear cell samples [40].However, its function in AD is unclear.Knockdown of RPL32 significantly represses hepa-tocellular carcinoma proliferation [41].RPL24 is another translation factor related to tumorigenesis.RPL24 reduction was found to be associated with over-expression of p53, suggesting that RPL24 in tumorigenesis is p53 dependent [42].RPL35 encodes a protein of the 60s ribosomal subunit, which is part of the ribosomal protein L29 family.The L29 family proteins have been found to be usually related to the activation of oncogenes and tumor suppressor genes.However, little study of RPL35 function on neurodegenerative disease development has been found [43,44].RPL31, a 60S large ribosomal subunit protein, was demonstrated to be correlated with the cell cycle and regulate cancer development and progression via the p53 pathway [45].Prior literature has shown that RPL35A was the first to be involved in Diamond-Blackfan anemia, which was found to act as a biomarker in tumor angiogenesis [46].RPL35A has also been reported to be a significant down-regulator in an APP/PS1 mouse model of AD [47].However, little is known about the relationship between RPL35A and AD development.
Three upregulated miRNAs (hsa-miR-3158-3p, hsa-miR-4435, and hsa-miR-3605-3p) were screened from AD samples in this ceRNA network.Among them, previous research has reported the upregulation of hsa-miR-4435 in AD samples [48].This study was consistent with what is reported here.Additionally, the relationship of hsa-miR-3158-3p and hsa-miR-3605-3p with the pathogenesis of AD has not previously been reported.Similarly, five downregulated lncRNAs (TRHDE-AS1, SNHG10, OIP5-AS, LINC00926, and LINC00662) were found in AD samples.Three of them (SNHG10, OIP5-AS, and LINC00662) have been found to be correlated to the progression of AD by LncRNADisease v2.0.LINC00662 is a newly discovered long noncoding RNA.Previous study has shown that LINC00662 knockdown lowers the permeability of the blood brain barrier in the AD microenvironment due to enhanced tight junction-related protein expression levels [49].Moreover, the ceRNA network constructed here predicted 26 differentially expressed circRNAs (see Table 4) that have been identified to be downregulated in AD.Among them, circRNA_0051620 was found to facilitate the occurrence and development of gastric cancer via acting as a miR-338-3p sponge and inhibiting ADAM17 [50].However, there is little research on other circRNAs.
Notably, this study identified the clinical significance of several lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA axes that included RPL35A and RPL35.LncRNA SNHG10 and OIP5-AS1 might be helpful for the development of AD via acting as a hsa-miR-3158-3p sponge and inhibiting the target mRNA RPL35A expression level.LncRNA LINC00662 might help the development of AD via acting as a sponge of hsa-miR-3605-3p and inhibiting the target mRNA RPL35A expression level.LncRNA OIP5-AS1 and LINC00662 might be helpful in the development of AD via acting as a sponge of hsa-miR-4435 and regulating target mRNA RPL35 expression level.Six differentially expressed circRNAs (hsa_circRNA_0090682, hsa_circRNA_0100723, hsa_circRNA_0076336, hsa_circRNA_0100174, hsa_circRNA_0001959, and hsa_circRNA_0104352) might regulate the development of AD via acting as a sponge of hsa-miR-3158-3p and inhibiting the target mRNA RPL35A expression.11 circR-NAs with differential expression (hsa_circRNA_0101760, hsa_circRNA_0002172, hsa_circRNA_0101758, hsa_circRNA_0103108, hsa_circRNA_0020990, hsa_circRNA_0051620, hsa_circRNA_0101472, hsa_circRNA_0002640, hsa_circRNA_0003057, hsa_circRNA_0102589, and hsa_circRNA_0104911) might regulate the development of AD via acting as a sponge of hsa-miR-4435 and inhibiting the expression level of target mRNA RPL35A and RPL35.Two circRNAs with differential expression (hsa_circRNA_0073593 and hsa_circRNA_0082335) might regulate the development of AD by sponging hsa-let-7d-3p and inhibiting the target mRNA RPL35A and RPL35 expression level.Ten differentially expressed circRNAs (hsa_circRNA_0090682, hsa_circRNA_0100946, hsa_circRNA_0104911, hsa_circRNA_0001276, hsa_circRNA_0100021, hsa_circRNA_0001399, hsa_circRNA_0026074, hsa_circRNA_0066331, hsa_circRNA_0076336, and hsa_circRNA_0100723) were found to be correlated to AD development via acting as a sponge of hsa-miR-330-5p and inhibiting the expression level of the target mRNA RPL35A.Hsa_circRNA_0025332 was found to be related to AD progression via acting as a sponge of hsa-miR-3605-3p and inhibiting the target mRNA RPL35A expression.Additionally, this analysis using the AlzData database showed that the target mRNA RPL35 and RPL35A were significantly associated with the clinicopathological features of AD.This study could provide new insights for understanding the interaction of ceRNA regulation and the pathogenesis of AD.

Note:
eQTL: target gene expression is affected by variation of the AD gene (p < 1 −3 ).GWAS: p < 1 −3 .PPI: interaction between target genes and APP, PSEN1, PSEN2, APOE, or MAPT is significant (p < 0.05); Early_DEG: differential expression of target gene in AD mouse models before the onset of AD pathology; Pathology correlation (Aβ): correlation between expression of target gene and AD pathobiology in AD mouse models of Aβ line; Pathology correlation (tau): correlation between expression of target gene and AD pathobiology in AD mouse models of tau line.ns, Non-significant; * p < 0.05; *** p < 0.001; CFG: total CFG score of target gene.CFG score is increased based on the importance of the above evidence.

Fig. 4 .
Fig. 4. Construction of the PPI network and analysis of modules.(A-F) PPI network of for the selected 6 modules.Red: Upregulated DEmRNAs.Green: Downregulated DEmRNAs.(G-L) Enrichment of KEGG pathway of these modules.DEmRNAs, differentially expressed mRNAs; KEGG, Kyoto Encyclopedia of Genes and Genomes.

Fig. 5 .
Fig. 5. Verification of hub genes.(A) Regulatory network of the hub genes.(B) KEGG pathway enrichment of hub genes.(C-K) Hub gene expression in different parts of the brain.

Table 2 . The top ten genes listed as hub genes through the maximum group centrality (MCC) algorithm.
× 10 53 downregulated