## 3 Results and Discussion
The enriched GO terms were divided into biological process (BP), molecular function (MF), and cellular component (CC) ontologies (Tables S1-S4). The results of the GO analysis indicated that DEGs were mainly enriched in BP, including mRNA metabolic process, oxidative RNA demethylation, regulation of core promoter binding, and glandular epithelial cell development (Figure S1A). MF analysis revealed that the DEGs were significantly enriched endoribonuclease activity, producing 3\({}^{\prime}\)-phosphomoesters, guanylate kinase activity, cAMP binding, and microtubule plus-end binding (Figure S1B). For CC, the DEGs were enriched in a collage-containing extracellular matrix, a basement membrane and endosome lumen (Figure S1C). The results of the Reactome pathway analysis showed that DEGs were mainly enriched in pathways such as in assembly of collagen fibrils extracellular matrix organization, collagen formation, and cell-cell communication (Figure S1D). More precisely, a _Dst_ (dystonin) gene encodes a member of the plakin protein family of adhesion junction plaque proteins acting as an integrator of intermediate filaments and microtubule cytoskeleton networks. A _Dst_ regulates the organization and stability of the microtubule network of sensory neurons to allow axonal transport and mediates docking of the dynein/dynactin motor complex to vesicle cargos for retrograde axonal transport through its interaction with TMEM108 and DCTN1, while loss of function can cause hereditary sensory, autonomic neuropathy type 6, and epidermolysis bullosa simplex [27]. Biological process analysis indicated that a _Dst_ was significantly enriched in hemidesmosome assembly and retrograde axonal transport (Figure S2A, Table S6). Deficits in the latter are associated with the pathogenesis of multiple neurodegenerative diseases, including amyotrophic lateral sclerosis [28]. Pathway enrichment analysis showed that a _Dst_ was mainly enriched in the RHO family of GTPases (Figure S2D, Table S8) which are involved in the regulation of cell migration and cell adhesion and play important roles in neuronal development [29]. \begin{table}
\begin{tabular}{c c c c c c c} \hline
**Filter Data** & **Accuracy** & **Kappa** & **Precision** & **Recall** & **F1** & **
\begin{tabular}{c} **Elapsed Time** \\ **(m)** \\ \end{tabular} \\ \hline Initial Data & 0.77 & 0.55 & 0.73 & 0.86 & 0.79 & 4.25 \\ Variable Data & 0.69 & 0.39 & 0.72 & 0.63 & 0.67 & 1.03 \\ Importance Data & 0.81 & 0.63 & 0.83 & 0.78 & 0.81 & 1.3 \\ Statistical Data & 0.60 & 0.21 & 0.72 & 0.34 & 0.46 & 1.01 \\ EFSM & 0.82 & 0.65 & 0.83 & 0.81 & 0.82 & 1 \\ \hline \end{tabular}
\end{table} 
Table 1: Prediction accuracy with kNN classifier prediction accuracy based on 200 informative genes by each gene selection method and initial data. Figure 2: A representation of the receiver operating characteristics (ROC) indicating that the classifier can be organized and its performance better understood. (**A**) AUC-ROC, (**B**) AUC-PRG, (**C**) calibration curves. A _Rapgef4_ (Rap guanine nucleotide exchange factor 4) gene is involved in the regulation of neuronal action potential and the development of the nervous system and plays a role in postsynaptic density and glutamatergic synapse. For BF, _Rapgef4_ was enriched in the regulation of peptide hormone secretion, the regulation of insulin and protein secretion, respectively (Figure S3A, Table S9). Peptide hormones are also the hypothalamic peptide hormones such as CRH, GHIGH, and TRH and their function is implicated in the regulation of peptide-containing secretory neurons [30]. The results of MF analysis revealed that _Rapgef4_ was mainly enriched cyclic adenosine monophosphate (cAMP) binding, cyclic nucleotide binding, and guanyl-nucleotide exchange factor activity (Figure S3B, Table S10). It should be noticed that cAMP is involved in neuronal functionality and is known to activate and integrate a variety of downstream pathways.