Differentially expressed gene (DEG) based protein-protein interaction (PPI) network identifies a spectrum of gene interactome, transcriptome and correlated miRNA in nondisjunction Down syndrome

A comprehensive integrated gene network construction to explore the essential role of Notch 1 in lung adenocarcinoma (LUAD)

The heterogeneous biological landscape of non-small cell lung cancer (NSCLC) is largely attributed to the activation of Notch signalling pathway. Among the Notch family transmembrane proteins, neurogenic locus notch homolog protein1 (NOTCH1) is a putative oncogene in NSCLC which activates the pathway as negative prognostic factor. This study aims to explore integrated network approach in lung adenocarcinoma (LUAD) especially linked to the notch pathway and its receptors. Our gene set enrichment analysis reveals the key Notch pathway genes are predominantly down regulated in LUAD. There were 675 genes with a total of 6517 functional interactions and 6 densely connected clusters of 38 miRNAs, 84 transcription factors with 156 edges identified through network construction. Here we report five key genes namely NOTCH1, CDH1, ERBB2, GAPDH and COL1A1 significantly enriched in Notch pathway which are further validated through the KM plot, box plots, stage plots and TIMER analysis. In addition, the NOTCH1 receptor is strongly linked to the immune checkpoint inhibitor CD274 (PD-L1) and can be considered as prognostic marker and tumour suppressor gene in LUAD which surely provide the basis for early diagnosis and futuristic immunotherapeutic targets for LUAD.Communicated by Ramaswamy H. Sarma.

FYN: emerging biological roles and potential therapeutic targets in cancer

Src family protein kinases (SFKs) play a key role in cell adhesion, invasion, proliferation, survival, apoptosis, and angiogenesis during tumor development. In humans, SFKs consists of eight family members with similar structure and function. There is a high level of overexpression or hyperactivity of SFKs in tumor, and they play an important role in multiple signaling pathways involved in tumorigenesis. FYN is a member of the SFKs that regulate normal cellular processes. Additionally, FYN is highly expressed in many cancers and promotes cancer growth and metastasis through diverse biological functions such as cell growth, apoptosis, and motility migration, as well as the development of drug resistance in many tumors. Moreover, FYN is involved in the regulation of multiple cancer-related signaling pathways, including interactions with ERK, COX-2, STAT5, MET and AKT. FYN is therefore an attractive therapeutic target for various tumor types, and suppressing FYN can improve the prognosis and prolong the life of patients. The purpose of this review is to provide an overview of FYN's structure, expression, upstream regulators, downstream substrate molecules, and biological functions in tumors.

Preliminary exploration of the co-regulation of Alzheimer's disease pathogenic genes by microRNAs and transcription factors

Background

Alzheimer's disease (AD) is the most common form of age-related neurodegenerative disease. Unfortunately, due to the complexity of pathological types and clinical heterogeneity of AD, there is a lack of satisfactory treatment for AD. Previous studies have shown that microRNAs and transcription factors can modulate genes associated with AD, but the underlying pathophysiology remains unclear.

Methods

The datasets GSE1297 and GSE5281 were downloaded from the gene expression omnibus (GEO) database and analyzed to obtain the differentially expressed genes (DEGs) through the "R" language "limma" package. The GSE1297 dataset was analyzed by weighted correlation network analysis (WGCNA), and the key gene modules were selected. Next, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis for the key gene modules were performed. Then, the protein-protein interaction (PPI) network was constructed and the hub genes were identified using the STRING database and Cytoscape software. Finally, for the GSE150693 dataset, the "R" package "survivation" was used to integrate the data of survival time, AD transformation status and 35 characteristics, and the key microRNAs (miRNAs) were selected by Cox method. We also performed regression analysis using least absolute shrinkage and selection operator (Lasso)-Cox to construct and validate prognostic features associated with the four key genes using different databases. We also tried to find drugs targeting key genes through DrugBank database.

Results

GO and KEGG enrichment analysis showed that DEGs were mainly enriched in pathways regulating chemical synaptic transmission, glutamatergic synapses and Huntington's disease. In addition, 10 hub genes were selected from the PPI network by using the algorithm Between Centrality. Then, four core genes (TBP, CDK7, GRM5, and GRIA1) were selected by correlation with clinical information, and the established model had very good prognosis in different databases. Finally, hsa-miR-425-5p and hsa-miR-186-5p were determined by COX regression, AD transformation status and aberrant miRNAs.

Conclusion

In conclusion, we tried to construct a network in which miRNAs and transcription factors jointly regulate pathogenic genes, and described the process that abnormal miRNAs and abnormal transcription factors TBP and CDK7 jointly regulate the transcription of AD central genes GRM5 and GRIA1. The insights gained from this study offer the potential AD biomarkers, which may be of assistance to the diagnose and therapy of AD.

Identification of hub biomarkers of myocardial infarction by single-cell sequencing, bioinformatics, and machine learning

Background

Myocardial infarction (MI) is one of the first cardiovascular diseases endangering human health. Inflammatory response plays a significant role in the pathophysiological process of MI. Messenger RNA (mRNA) has been proven to play a key role in cardiovascular diseases. Single-cell sequencing (SCS) technology is a new technology for high-throughput sequencing analysis of genome, transcriptome, and epigenome at the single-cell level, and it also plays an important role in the diagnosis and treatment of cardiovascular diseases. Machine learning algorithms have a wide scope of utilization in biomedicine and have demonstrated superior efficiency in clinical trials. However, few studies integrate these three methods to investigate the role of mRNA in MI. The aim of this study was to screen the expression of mRNA, investigate the function of mRNA, and provide an underlying scientific basis for the diagnosis of MI.

Methods

In total, four RNA microarray datasets of MI, namely, GSE66360, GSE97320, GSE60993, and GSE48060, were downloaded from the Gene Expression Omnibus database. The function analysis was carried out by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Disease Ontology (DO) enrichment analysis. At the same time, inflammation-related genes (IRGs) were acquired from the GeneCards database. Then, 52 co-DEGs were acquired from differentially expressed genes (DEGs) in differential analysis, IRGs, and genes from SCS, and they were used to construct a protein-protein interaction (PPI) network. Two machine learning algorithms, namely, (1) least absolute shrinkage and selection operator and (2) support vector machine recursive feature elimination, were used to filter the co-DEGs. Gene set enrichment analysis (GSEA) was performed to screen the hub-modulating signaling pathways associated with the hub genes. The results were validated in GSE97320, GSE60993, and GSE48060 datasets. The CIBERSORT algorithm was used to analyze 22 infiltrating immune cells in the MI and healthy control (CON) groups and to analyze the correlation between these immune cells. The Pymol software was used for molecular docking of hub DEGs and for potential treatment of MI drugs acquired from the COREMINE.

Results

A total of 126 DEGs were in the MI and CON groups. After screening two machine learning algorithms and key co-DEGs from a PPI network, two hub DEGs (i.e., IL1B and TLR2) were obtained. The diagnostic efficiency of IL1B, TLR2, and IL1B + TLR2 showed good discrimination in the four cohorts. GSEA showed that KEGG enriched by DEGs were mainly related to inflammation-mediated signaling pathways, and GO biological processes enriched by DEGs were linked to biological effects of various inflammatory cells. Immune analysis indicated that IL1B and TLR2 were correlated with various immune cells. Dan shen, san qi, feng mi, yuan can e, can sha, san qi ye, san qi hua, and cha shu gen were identified as the potential traditional Chinese medicine (TCM) for the treatment of MI. 7-hydroxyflavone (HF) had stable combinations with IL1B and TLR2, respectively.

Conclusion

This study identified two hub DEGs (IL1B and TLR2) and illustrated its potential role in the diagnosis of MI to enhance our knowledge of the underlying molecular mechanism. Infiltrating immune cells played an important role in MI. TCM, especially HF, was a potential drug for the treatment of MI.

Neurons: The Interplay between Cytoskeleton, Ion Channels/Transporters and Mitochondria

Neurons are permanent cells whose key feature is information transmission via chemical and electrical signals. Therefore, a finely tuned homeostasis is necessary to maintain function and preserve neuronal lifelong survival. The cytoskeleton, and in particular microtubules, are far from being inert actors in the maintenance of this complex cellular equilibrium, and they participate in the mobilization of molecular cargos and organelles, thus influencing neuronal migration, neuritis growth and synaptic transmission. Notably, alterations of cytoskeletal dynamics have been linked to alterations of neuronal excitability. In this review, we discuss the characteristics of the neuronal cytoskeleton and provide insights into alterations of this component leading to human diseases, addressing how these might affect excitability/synaptic activity, as well as neuronal functioning. We also provide an overview of the microscopic approaches to visualize and assess the cytoskeleton, with a specific focus on mitochondrial trafficking.

Down-Regulation of Insulin Like Growth Factor 1 Involved in Alzheimer’s Disease via MAPK, Ras, and FoxO Signaling Pathways

The inability to halt or even delay the course of Alzheimer's disease (AD) forces the development of new molecular signatures and therapeutic strategies. Insulin like growth factor 1 (IGF1) is a promising target for AD treatment, yet exact mechanisms of AD ascribed to IGF1 remain elusive. Herein, gene expression profiles of 195 samples were analyzed and 19,245 background genes were generated, among which 4,424 differentially expressed genes (DEGs) were overlapped between AD/control and IGF1-low/high groups. Based on such DEGs, seven co-expression modules were established by weight gene correlation network analysis (WGCNA). The turquoise module had the strongest correlation with AD and IGF1-low expression, the DEGs of which were enriched in GABAergic synapse, long-term potentiation, mitogen-activated protein kinase (MAPK), Ras, and forkhead box O (FoxO) signaling pathways. Furthermore, cross-talking pathways of IGF1, including MAPK, Ras, and FoxO signaling pathways were identified in the protein-protein interaction network. According to the area under the curve (AUC) analysis, down-regulation of IGF1 exhibited good diagnostic performance in AD prediction. Collectively, our findings highlight the involvement of low IGF1 in AD pathogenesis via MAPK, Ras, and FoxO signaling pathways, which might advance strategies for the prevention and therapy of AD based on IGF1 target.

GNG2 acts as a tumor suppressor in breast cancer through stimulating MRAS signaling

G-protein gamma subunit 2 (GNG2) is involved in several cell signaling pathways, and is essential for cell proliferation and angiogenesis. However, the role of GNG2 in tumorigenesis and development remains unclear. In this study, 1321 differentially expressed genes (DEGs) in breast cancer (BC) tissues were screened using the GEO and TCGA databases. KEGG enrichment analysis showed that most of the enriched genes were part of the PI3K-Akt signaling pathway. We identified GNG2 from the first five DEGs, its expression was markedly reduced in all BC subtype tissues. Cox regression analysis showed that GNG2 was independently associated with overall survival in patients with luminal A and triple-negative breast cancers (TNBC). GNG2 over-expression could significantly block the cell cycle, inhibit proliferation, and promote apoptosis in BC cells in vitro. In animal studies, GNG2 over-expression inhibited the growth of BC cells. Further, we found that GNG2 significantly inhibited the activity of ERK and Akt in an MRAS-dependent manner. Importantly, GNG2 and muscle RAS oncogene homolog (MRAS) were co-localized in the cell membrane, and the fluorescence resonance energy transfer (FRET) experiment revealed that they had direct interaction. In conclusion, the interaction between GNG2 and MRAS likely inhibits Akt and ERK activity, promoting apoptosis and suppressing proliferation in BC cells. Increasing GNG2 expression or disrupting the GNG2-MRAS interaction in vivo could therefore be a potential therapeutic strategy to treat BC.

Proteomic mechanisms for the regulation of growth, photosynthetic activity and nitrogen fixation in Nostoc sp. PCC 7120 exposed to three antibiotic contaminants

This study investigated the influences of three frequently detected antibiotics in surface waters, ciprofloxacin, tetracycline and sulfamethoxazole, on the growth, photosynthetic activity, nitrogen-fixing capacity and proteomic expression profiles of Nostoc sp. PCC 7120, through a 15-day exposure test at environmentally relevant exposure doses of 50-200 ng/L. Cyanobacterial growth was stimulated by 100 ng/L and 200 ng/L of ciprofloxacin and sulfamethoxazole as well as 50-200 ng/L of tetracycline. The nitrogenase synthesis ability in each cyanobacterial cell was stimulated by 50-200 ng/L of ciprofloxacin while inhibited by 100 ng/L and 200 ng/L of tetracycline and sulfamethoxazole. At the exposure dose of 100 ng/L for each antibiotic, the variation of total nitrogen in the culture medium indicated that the nitrogen-fixing capacity of Nostoc sp. was determined by total nitrogenase concentration calculated by cell density × nitrogenase synthesis ability. Therefore, ciprofloxacin enhanced nitrogen fixation through the stimulation of both cyanobacterial growth and nitrogenase synthesis, while tetracycline and sulfamethoxazole enhanced nitrogen fixation merely through growth stimulation. At the exposure dose of 100 ng/L, only two downregulated proteins, a phosphonate ABC transporter and a methionine aminopeptidase, as well as one upregulated protein, the phenylalanine-tRNA ligase alpha subunit, were commonly shared by three antibiotic-treated groups. Ciprofloxacin upregulated proteins related to nitrogen fixation, carbon catabolism and biosynthesis, but downregulated photosynthesis-related proteins. In contrast, tetracycline and sulfamethoxazole increased the photosynthetic activity of Nostoc sp. through upregulating photosynthesis-related proteins, but downregulated proteins related to nitrogen fixation, carbon catabolism and biosynthesis. The resistance of Nostoc sp. PCC 7120 to three target antibiotics were related with the responses of RNA synthesis regulatory proteins. Stimulation of cyanobacterial nitrogen fixation by antibiotic contaminants could aggravate eutrophication in aquatic environments.

Transcriptome profiling of five brain regions in a 6-hydroxydopamine rat model of Parkinson's disease

BACKGROUND

Parkinson's disease (PD) is a neurodegenerative disease, and its pathogenesis is unclear. Previous studies mainly focus on the lesions of substantia nigra (SN) and striatum (Str) in PD. However, lesions are not limited. The olfactory bulb (OB), subventricular zone (SVZ), and hippocampus (Hippo) are also affected in PD.

AIM

To reveal gene expression changes in the five brain regions (OB, SVZ, Str, SN, and Hippo), and to look for potential candidate genes and pathways that may be correlated with the pathogenesis of PD.

MATERIALS AND METHODS

We established control group and 6-hydroxydopamine (6-OHDA) PD model group, and detected gene expressions in the five brain regions using RNA-seq and real-time quantitative polymerase chain reaction (RT-qPCR). We further analyzed the RNA-seq data by bioinformatics.

RESULTS

We identified differentially expressed genes (DEGs) in all five brain regions. The DEGs were significantly enriched in the "dopaminergic synapse" and "retrograde endocannabinoid signaling," and Gi/o-GIRK is the shared cascade in the two pathways. We further identified Ephx2, Fam111a, and Gng2 as the potential candidate genes in the pathogenesis of PD for further studies.

CONCLUSION

Our study suggested that gene expressions change in the five brain regions following exposure to 6-OHDA. The "dopaminergic synapse," "retrograde endocannabinoid signaling," and Gi/o-GIRK may be the key pathways and cascade of the synaptic damage in 6-OHDA PD rats. Ephx2, Fam111a, and Gng2 may play critical roles in the pathogenesis of PD.

Profiling Basal Forebrain Cholinergic Neurons Reveals a Molecular Basis for Vulnerability Within the Ts65Dn Model of Down Syndrome and Alzheimer's Disease

Basal forebrain cholinergic neuron (BFCN) degeneration is a hallmark of Down syndrome (DS) and Alzheimer's disease (AD). Current therapeutics have been unsuccessful in slowing disease progression, likely due to complex pathological interactions and dysregulated pathways that are poorly understood. The Ts65Dn trisomic mouse model recapitulates both cognitive and morphological deficits of DS and AD, including BFCN degeneration. We utilized Ts65Dn mice to understand mechanisms underlying BFCN degeneration to identify novel targets for therapeutic intervention. We performed high-throughput, single population RNA sequencing (RNA-seq) to interrogate transcriptomic changes within medial septal nucleus (MSN) BFCNs, using laser capture microdissection to individually isolate ~500 choline acetyltransferase-immunopositive neurons in Ts65Dn and normal disomic (2N) mice at 6 months of age (MO). Ts65Dn mice had unique MSN BFCN transcriptomic profiles at ~6 MO clearly differentiating them from 2N mice. Leveraging Ingenuity Pathway Analysis and KEGG analysis, we linked differentially expressed gene (DEG) changes within MSN BFCNs to several canonical pathways and aberrant physiological functions. The dysregulated transcriptomic profile of trisomic BFCNs provides key information underscoring selective vulnerability within the septohippocampal circuit. We propose both expected and novel therapeutic targets for DS and AD, including specific DEGs within cholinergic, glutamatergic, GABAergic, and neurotrophin pathways, as well as select targets for repairing oxidative phosphorylation status in neurons. We demonstrate and validate this interrogative quantitative bioinformatic analysis of a key dysregulated neuronal population linking single population transcript changes to an established pathological hallmark associated with cognitive decline for therapeutic development in human DS and AD.

Neuronal Cytoskeleton in Intellectual Disability: From Systems Biology and Modeling to Therapeutic Opportunities

Intellectual disability (ID) is a pathological condition characterized by limited intellectual functioning and adaptive behaviors. It affects 1–3% of the worldwide population, and no pharmacological therapies are currently available. More than 1000 genes have been found mutated in ID patients pointing out that, despite the common phenotype, the genetic bases are highly heterogeneous and apparently unrelated. Bibliomic analysis reveals that ID genes converge onto a few biological modules, including cytoskeleton dynamics, whose regulation depends on Rho GTPases transduction. Genetic variants exert their effects at different levels in a hierarchical arrangement, starting from the molecular level and moving toward higher levels of organization, i.e., cell compartment and functions, circuits, cognition, and behavior. Thus, cytoskeleton alterations that have an impact on cell processes such as neuronal migration, neuritogenesis, and synaptic plasticity rebound on the overall establishment of an effective network and consequently on the cognitive phenotype. Systems biology (SB) approaches are more focused on the overall interconnected network rather than on individual genes, thus encouraging the design of therapies that aim to correct common dysregulated biological processes. This review summarizes current knowledge about cytoskeleton control in neurons and its relevance for the ID pathogenesis, exploiting in silico modeling and translating the implications of those findings into biomedical research.

Investigating the effects of chronic perinatal alcohol and combined nicotine and alcohol exposure on dopaminergic and non-dopaminergic neurons in the VTA

The ventral tegmental area (VTA) is the origin of dopaminergic neurons and the dopamine (DA) reward pathway. This pathway has been widely studied in addiction and drug reinforcement studies and is believed to be the central processing component of the reward circuit. In this study, we used a well-established rat model to expose mother dams to alcohol, nicotine-alcohol, and saline perinatally. DA and non-DA neurons collected from the VTA of the rat pups were used to study expression profiles of miRNAs and mRNAs. miRNA pathway interactions, putative miRNA-mRNA target pairs, and downstream modulated biological pathways were analyzed. In the DA neurons, 4607 genes were differentially upregulated and 4682 were differentially downregulated following nicotine-alcohol exposure. However, in the non-DA neurons, only 543 genes were differentially upregulated and 506 were differentially downregulated. Cell proliferation, differentiation, and survival pathways were enriched after the treatments. Specifically, in the PI3K/AKT signaling pathway, there were 41 miRNAs and 136 mRNAs differentially expressed in the DA neurons while only 16 miRNAs and 20 mRNAs were differentially expressed in the non-DA neurons after the nicotine-alcohol exposure. These results depicted that chronic nicotine and alcohol exposures during pregnancy differentially affect both miRNA and gene expression profiles more in DA than the non-DA neurons in the VTA. Understanding how the expression signatures representing specific neuronal subpopulations become enriched in the VTA after addictive substance administration helps us to identify how neuronal functions may be altered in the brain.

Transcriptomic analysis delineates potential signature genes and miRNAs associated with the pathogenesis of asthma

Asthma is a multifarious disease affecting several million people around the world. It has a heterogeneous risk architecture inclusive of both genetic and environmental factors. This heterogeneity can be utilised to identify differentially expressed biomarkers of the disease, which may ultimately aid in the development of more localized and molecularly targeted therapies. In this respect, our study complies with meta-analysis of microarray datasets containing mRNA expression profiles of both asthmatic and control patients, to identify the critical Differentially Expressed Genes (DEGs) involved in the pathogenesis of asthma. We found a total of 30 DEGs out of which 13 were involved in the pathway and functional enrichment analysis. Moreover, 5 DEGs were identified as the hub genes by network centrality-based analysis. Most hub genes were involved in protease/antiprotease pathways. Also, 26 miRNAs and 20 TFs having an association with these hub genes were found to be intricated in a 3-node miRNA Feed-Forward Loop. Out of these, miR-34b and miR-449c were identified as the key miRNAs regulating the expression of SERPINB2 gene and SMAD4 transcription factor. Thus, our study is suggestive of certain miRNAs and unexplored pathways which may pave a way to unravel critical therapeutic targets in asthma.

Cystathionine-β-Synthase: Molecular Regulation and Pharmacological Inhibition

Cystathionine-β-synthase (CBS), the first (and rate-limiting) enzyme in the transsulfuration pathway, is an important mammalian enzyme in health and disease. Its biochemical functions under physiological conditions include the metabolism of homocysteine (a cytotoxic molecule and cardiovascular risk factor) and the generation of hydrogen sulfide (H2S), a gaseous biological mediator with multiple regulatory roles in the vascular, nervous, and immune system. CBS is up-regulated in several diseases, including Down syndrome and many forms of cancer; in these conditions, the preclinical data indicate that inhibition or inactivation of CBS exerts beneficial effects. This article overviews the current information on the expression, tissue distribution, physiological roles, and biochemistry of CBS, followed by a comprehensive overview of direct and indirect approaches to inhibit the enzyme. Among the small-molecule CBS inhibitors, the review highlights the specificity and selectivity problems related to many of the commonly used "CBS inhibitors" (e.g., aminooxyacetic acid) and provides a comprehensive review of their pharmacological actions under physiological conditions and in various disease models.

Role of 3-Mercaptopyruvate Sulfurtransferase in the Regulation of Proliferation and Cellular Bioenergetics in Human Down Syndrome Fibroblasts

Down syndrome (trisomy of human chromosome 21) is a common genetic disorder. Overproduction of the gaseous mediator hydrogen sulfide (H₂S) has been implicated in the pathogenesis of neurological and metabolic deficits associated with Down syndrome. Several lines of data indicate that an important enzyme responsible for H₂S overproduction in Down syndrome is cystathionine-β-synthase (CBS), an enzyme localized on chromosome 21. The current study explored the possibility that a second H₂S-producing enzyme, 3-mercaptopyruvate sulfurtransferase (3-MST), may also contribute to the development of functional deficits of Down syndrome cells. Western blotting analysis demonstrated a significantly higher level of 3-MST protein expression in human Down syndrome fibroblasts compared to cells from healthy control individuals; the excess 3-MST was mainly localized to the mitochondrial compartment. Pharmacological inhibition of 3-MST activity improved mitochondrial electron transport and oxidative phosphorylation parameters (but did not affect the suppressed glycolytic parameters) and enhanced cell proliferation in Down syndrome cells (but not in healthy control cells). The findings presented in the current report suggest that in addition to the indisputable role of CBS, H₂S produced from 3-MST may also contribute to the development of mitochondrial metabolic and functional impairments in Down syndrome cells.

Cognitive impairment and transcriptomic profile in hippocampus of young mice after multiple neonatal exposures to sevoflurane

Children with repeated inhalational anesthesia may develop cognitive disorders. This study aimed to investigate the transcriptome-wide response of hippocampus in young mice that had been exposed to multiple sevoflurane in the neonatal period. Mice received 3% sevoflurane for 2 h on postnatal day (PND) 6, 8, and 10, followed by arterial blood gas test on PND 10, behavioral experiments on PND 31–36, and RNA sequencing (RNA-seq) of hippocampus on PND 37. Functional annotation and protein-protein interaction analyses of differentially expressed genes (DEGs) and quantitative reverse transcription polymerase chain reaction (qPCR) were performed. Neonatal sevoflurane exposures induced cognitive and social behavior disorders in young mice. RNA-seq identified a total of 314 DEGs. Several enriched biological processes (ion channels, brain development, learning, and memory) and signaling pathways (oxytocin signaling pathway and glutamatergic, cholinergic, and GABAergic synapses) were highlighted. As hub-proteins, Pten was involved in nervous system development, synapse assembly, learning, memory, and behaviors, Nos3 and Pik3cd in oxytocin signaling pathway, and Cdk16 in exocytosis and phosphorylation. Some top DEGs were validated by qPCR. This study revealed a transcriptome-wide profile in mice hippocampus after multiple neonatal exposures to sevoflurane, promoting better understanding of underlying mechanisms and investigation of preventive strategies.