Data Collection and potential drugs screening
The gene expression profile (GSE2871) was obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/), which was sequenced using the GPL85 platform. At early post-injury timepoint, animals will be sacrificed, brain regions (parietal cortex and hippocampus, ipsilateral and contralateral to injury) will be dissected and RNA isolated. RNA will be used to synthesize cRNA probes for microarray hybridization. We identified a number of potential drugs that may act on some of the TBI-associated genes identified in this study in HERB database (http://herb.ac.cn/) and the other public database (https://drugcentral.org/).
Identifification of DEGs and GSEA
The “Limma” R package was used to screen differentially expressed genes (DEGs) between TBI and normal samples, and genes with P < 0.05 and |log2FC| >1 were regarded as DEGs. GSEA-4.1.0 was used to input the expression data and phenotypic data, and the five most significantly up-regulated pathways and the most significantly down-regulated pathways were plotted, respectively.
Screening of the Critical Genes
To find out the core genes that were altered after TBI, the downloaded dataset was used to construct a weighted gene co-expression network using the “WGCNA” R package. Then, we performed a cluster analysis of the samples and calculated the pearson correlation coefficient between each pair of genes to evaluate the expression similarity of genes and acquire a correlation matrix. We further used a soft threshold function to transform the correlation matrix into a weighted neighborhood matrix, and a soft join algorithm was used to select the optimal soft threshold to ensure that gene correlations fit the scale-free distribution to the greatest extent possible. Subsequently, the neighborhood matrix was transformed into a topological overlap matrix. After obtaining the co-expression modules, the key modules were screened out by correlation analysis, and the genes of the key modules were regarded as the key genes. Based on WGCNA screening, DEGs and key genes were intersected to obtain the target genes.
Identification of TBI hub genes based on machine-learning algorithms
Least absolute shrinkage and selection operator (LASSO) is a regression analysis method that performs both gene selection and classification[13]. First, the R package glmnet (Version4.1.2) was used to fit the logistic LASSO regression model. Next, the SVM-RFE algorithm was used to screen potential genes using the “e1071” R package. In addition, the random forest (RF) algorithm was also conducted to screen potential genes using the “randomForest” R package. Finally, the intersection of the genes obtained by LASSO, SVM-RFE and RF machine-learning algorithms was taken by Veen graph as the hub genes of TBI.
Experimental animals and grouping
The male ICR mice (6–8 weeks, 20-30g) used in this experiment were purchased from Hangzhou Ziyuan Experimental Animal Technology Co., LTD., and fed in the mouse house for 7 days for the experiment. The temperature of the mouse house is controlled at 22 ~ 25◦C, the humidity is 70% ~ 75%, and the water and food are adequately supplied. All animal experiments were conducted in accordance with the animal welfare policy of The Affiliated Huai'an Hospital of Xuzhou Medical University and were approved by The Affiliated Huai'an Hospital of Xuzhou Medical University. The experimental animals were randomly divided into Sham group, TBI group and TBI + Cal group (Gallic acid, 40mg/kg, ip). In Sham group, the bone window was opened but not traumatized. The TBI group was traumatized but not treated. Gallic acid (MedChemExpress, USA) is administered once 30 minutes after TBI and then once daily. The experimental section of this article fully complies with the ARRIVE guidelines.
The construction of TBI model
The Feeney free-fall impact method was used to construct the TBI model[29]. First, the mice were anesthetized with isoflurane by an inhalation anesthesia machine (M5209, Changsha Maiyue Biotechnology Co., Ltd.) and fixed in a stereoscope, then the local hair was removed and disinfected with iodophor. Then the scalp was cut about 1.5cm along the median line of the mouse skull, exposing the skull and stripping the periosteum. Drilling was then done with a 2mm diameter cranial rotatory to create a bone window of about 5mm diameter, exposing the parietal lobe with as little damage to the dura as possible. Parameters for setting up the craniocerebral injury percussion apparatus are as follows: weight mass 20g, drop height 15cm, impact depth 1.5mm. After the attack, the mice were hemostatic and the scalp was tightly sutured, and then returned to the original cage for feeding. All operations should pay attention to the principle of asepsis, and the instruments should be autoclaved and disinfected in advance during the operation.
Quantitative Real-Time PCR
Total RNA extraction and cDNA reverse transcription of mouse brain tissue were performed by the Vazyme (Vazyme Biotech, Nanjing, China) kit. Using GAPDH as internal parameter, 2−ΔΔCt method was used to calculate the expression levels. Primer sequences for target genes:
STK39 (Forward)5’ - CAAACCCAGGCAAGAACGC − 3’ ;
STK39 (Reverse)5’ - GCCACAGCTCATCTTTGACCAC − 3’ ;
Kcnd3 (Forward)5’ - CACCAGTCGCTCCAGCCTTAA − 3’ ;
Kcnd3 (Reverse)5’ - GACGACATTGCTGGTTATGGAAG − 3’ ;
Apoc3 (Forward)5’ - GAGTCCGATATAGCTGTGGTGG − 3’ ;
Apoc3 (Reverse)5’ - GTTGGTTGGTCCTCAGGGTTAG − 3’ ;
FOXE3 (Forward)5’ - CGACTGTTTCGTCAAGGTGC − 3’ ;
FOXE3 (Reverse)5’ - CGTTGTCGAACATGTCAGCG − 3’ ;
CHRNB1 (Forward)5’ - CCGTTATCCTTAGTGTTGTGGTC − 3’ ;
CHRNB1 (Reverse)5’ - AGTGATGTGGTTCAGGGAGTTG − 3’ ;
NPW (Forward)5’ - CTGCTAGAGCCTTCGGAGAGAC − 3’ ;
NPW (Reverse)5’ - ATCGGTTCTTGGGCCTGACA − 3’ ;
GAPDH (Forward)5’ - CCTCGTCCCGTAGACAAAATG − 3’ ;
GAPDH (Reverse)5’ - TGAGGTCAATGAAGGGGTCGT − 3’ .
Western blotting
The tissue blocks were washed with pre-cooled PBS for 2–3 times to remove the blood stain, cut into small pieces and placed in a homogenizing tube. 2 homogenizing beads of 4mm were added, and the lysate of 10 times the tissue volume was added (protease inhibitor was added before use), and homogenizing procedure was set for homogenizing. Take out the homogenated tube, place the ice lysate for 30min, and shake every 5min to ensure complete tissue cracking; The supernatant was collected by centrifugation at 12000rpm at 4℃ for 10min. Protein concentration was quantified using the protein quantification kit (Servicebio, China). Add a certain amount of loading buffer and boil at 95℃ for 10min to complete sample preparation. Then the sample is applied and electrophoresis is performed successively until the sample reaches the lower edge of the gel. Then the gel was transferred and closed for 1 ~ 2h. After closure, the primary antibody was incubated overnight (4℃). The next day, the secondary antibody was washed 3 times with PBS buffer for 5min each time, and then added to the incubator and incubated at room temperature for 30min; And then clean it again with PBS buffer for 3 times, 10min each time. Luminescent solution (Servicebio, China) was configured and protein band detection and gray scale analysis were performed by imaging system (CLINX, China). The main reagents include primary antibody TfR1 (1:1000, Proteintech, USA), NOX2 (1:1000, Proteintech, USA) and GPX4 (1:1000, Proteintech, USA), secondary antibody Polyclonal Goat Anti-Mouse IgG labeled by HRP (1:2000, Servicebio, China), HRP labeled goat anti-Rabbit secondary antibody (1:2000, Servicebio, China).
Transmission electron microscopy
2mm brain tissue around the site of brain injury and normal brain tissue at the same site were cut into 1mm3 and immersed in electron microscope fixative (G1102-100ML, Servicebio, China). The tissue blocks were fixed at room temperature for 2h away from light, and then stored in a refrigerator at 4℃. After the tissue blocks were removed, the surface fixing solution was washed with PBS (PH = 7.4) for 3 times, and the washing time was 15min each time. The rinsed tissue blocks were fixed with 0.1mol/L PBS (PH = 7.4) prepared with 1% osmic acid and then placed in a dark room at room temperature for 2h. After that, the tissue blocks were rinsed with PBS (PH = 7.4) for 3 times, and each rinsing time was 15min. Finally, different concentrations of ethanol and 100% acetone were dehydrated. After completing the above operations, the tissue was sliced with a thickness of 60-80nm. The tissue sections were double-stained with uranium lead and then left to dry overnight at room temperature. Finally, the sections were analyzed under transmission electron microscopy.
Immunohistochemistry
We first dewaxed the tissue sections according to the previous experimental procedures[30]. After rinsing with PBS for three times, it was incubated in hydrogen peroxide solution in the dark, and then rinsed with PBS solution for three times. The tissue was uniformly covered with 3% BSA in the tissue chemical circle and closed at room temperature for 30 min. The main reagents include primary antibody TfR1 (1:1000, Proteintech, USA), NOX2 (1:1000, Proteintech, USA) and GPX4 (1:1000, Proteintech, USA). After adding the primary antibody, place it in a refrigerator at 4°C overnight and rinse it the next day. After rinsing the primary antibody and adding the second antibody, incubate in room temperature for 50 minutes. Secondary antibody Polyclonal Goat Anti-Mouse IgG labeled by HRP (1:2000, Servicebio, China), HRP labeled goat anti-Rabbit secondary antibody (1:2000, Servicebio, China). After completing the above steps, retain the core by DAB color development method for 3 minutes and then rinse. Then wash with hematoxylin differentiation solution and then rinse with water, and finally, hematoxylin is rinsed after returning to blue. In the last step, the sections were dehydrated in anhydrous ethanol and closed under a microscope.
Perl’s staining
First, the paraffin sections of the brain tissue were dewaxed to water. Prepare the Prussian blue dye, add it to the section and dye for 1 hour. Then use distilled water to wash the excess dye on the surface of the section. The sections were stained with nuclear fast red staining solution for 3 minutes and then the excess dye solution was rinsed. Finally, the slices are dehydrated and sealed with neutral gum. Image J 1.53 software was used for image processing.
Nissl staining
The paraffin sections were dewaxed to water and then washed 3 times with PBS. The slices were then incubated by Nissl staining solution (Beyotime Biotechnology, China) for 10 min. The analysis results were photographed with an optical microscope (Nikon 80 i, Japan).
Fluoro-Jade B Staining
First, the paraffin sections are dewaxed to water, and then the FJB working liquid is added. The specific process was as follows: the pen circle was organized, 50% ice acetic acid was used as solvent, FJB working liquid (Merck, Germany) was configured at 1: 400, the diluted FJB green fluorescent probe was added, and the nucleus was re-dyed with DAPI after 4◦C overnight. Images were collected after sealing.
Lesion degree assessment
We first stained the tissue with HE. The simple process is as follows: 1. Dewaxing paraffin sections to water; 2. Hematoxylin staining; 3. eosin staining; 4. dehydration seal. We started the section from the defect edge to the normal tissue edge with a thickness of 30 µm. Then the brain tissue defect volume was calculated using NIH Image J software (Bethesda, MD, USA).
Evans blue extravasation assay
We measured the amount of EB dye in mouse brain tissue 3 days after TBI to assess the extent of BBB destruction. In simple terms, after injecting EB dye (2%, 2 µL/g) through the tail vein, the mice were anesthetized and then injected with PBS through the left ventricle of the heart to eliminate the localized dye from the sinus bleeding. The brain tissue was then weighed, then the sample was soaked in formamide solution and homogenized at a concentration of 200 mg of tissue per milliliter. 37◦C warm bath for 48h, centrifuge at 6000 rpm for 20 mins and remove the supernatant. The absorbance of the mixture at 632 nm was measured using a spectrophotometer (BioTek, Winooski, VT, USA).
Behavioral experiments
The perception and memory ability of mice were investigated by Novel object recognition (NOR) experiment[31]. Within 3 days of the training stage, an identical object (yellow cube) was placed in a box (100 cm in diameter and 50 cm in height) at a position 20 cm away from the wall, and then the mouse was placed in the middle of the two objects for it to explore by itself. In the experiment phase, a familiar object (yellow cube) was replaced with a new object (white cylinder), and then the mouse was placed in the middle of the two objects and allowed to explore freely for 5 min. New object recognition rate (NORI): the proportion of old object recognition time to all object recognition time.
The Morris water maze experiment is divided into training and experiment stages[30]. The whole process was recorded and evaluated using a video tracking system (Anhui Zhenghua Biological Instrument Equipment Co., LTD. : Huaibei, China). Specific experiments are as follows:During the training phase (the first 7 days), 4 training sessions are performed daily. The mice were randomly placed into the water from a quadrant facing the wall of the pool. If the mouse successfully found the platform within 90 seconds, it was allowed to stay on the platform for 10 seconds, and if the mouse failed to find the platform within the specified time, it was guided to the platform for 10 seconds with a guide stick. The water temperature was kept at 24◦C during the experiment. In the experimental stage (day 8 of the water maze test), the platform was removed from the water, and the mice were put into the pool from the opposite quadrant of the platform for a single test. Each mouse swam for 90 seconds. The movement track of the mice, the number of times they crossed the platform and the time they stayed in the quadrant where the platform was located were recorded.
Statistic analysis
Statistical analyses were performed and plotted using GraphPad Prism 8.4.2 software (San Diego, CA, USA) and SPSS 20.0. Data were presented as means ± SD. (*P < 0.05, **P < 0.01, ***P < 0.001).