Animals
Adult male Sprague‒Dawley rats (4 months old, weighing 350-490 × g) were purchased from Jiangsu Aniphe Biolaboratory (Jiangsu, China) and were maintained under specific-pathogen-free (SPF) conditions with free access to food and water and a 12-h light/12-h dark cycle in the Laboratory Animal Center (Guangzhou Medical University, Guangzhou, China). All procedures used in the present study were performed in accordance with the guidelines on animal care. All experiments were approved by the Animal Experiments and Ethics Committee of Guangzhou Medical University (No. 2020143, 21 August 2020). The rats were assigned randomly to five groups: the control group without injection, the sham group with injection of vehicle only, and the 3 6-OHDA groups injected with 6-OHDA for observation times of 2, 3 and 5 weeks, separately.
6-OHDA lesions
The neurotoxic rat model of PD was constructed by a unilateral injection of 6-OHDA into the medial forebrain bundle (MFB) as reported previously16. Using the rat brain atlas of Paxinos and Watson (Paxinos and Watson, 2007) as a guide, 6-OHDA (Sigma‒Aldrich, Shanghai, China) was injected into the MFB based on the following coordinates: anterior-posterior, -4.4 mm; medial-lateral, -1.2 mm; dorsal-ventral, -7.8 mm relative to bregma. The rats were placed in a stereotaxic device under 1% pentobarbital sodium anesthesia. Here, 6-OHDA was dissolved at 2 μg per 1 μL in vehicle solution that contained 0.2% ascorbic acid in an equal volume of 0.9% NaCl. In the PD rats, 4 μL of the 6-OHDA solution was administered at a rate of 400 nl/min using a Hamilton syringe equipped with a 30-gauge needle. The needle was left after injection for an additional 10 min and withdrawn at a speed of 1 mm/min. The same coordinates were used in rats undergoing sham operation, but only vehicle was delivered.
Apomorphine-induced rotation test
To evaluate the presence of 6-OHDA lesions, the apomorphine-induced turning rate was assessed according to a previously described protocol (0.5 mg/kg, i.p. apomorphine) 16. At 2, 3, and 5 weeks after 6-OHDA injection, the rats were injected with 0.5 mg/kg apomorphine (Sigma‒Aldrich, Shanghai, China) at 0.5 mg/ml dissolved in NaCl by intraperitoneal injection. Then, the rats were transferred individually to a transparent plastic cylinder, and the rotations were recorded for 30 min in a self-constructed rotometer according to a previously reported method16, 17. The total number of complete 360° anti-clockwise rotations was counted, and the mean number per minute was calculated and expressed in the Results section.
Immunohistochemical staining (ICH)
Rats were anesthetized with chloral hydrate and then perfusion-fixed with precooled 4% paraformaldehyde (Sigma‒Aldrich, Shanghai, China) followed by 0.9% NaCl. The brains were separated, postfixed in 4% PFA for 24 h and cryoprotected in 30% w/v sucrose in PBS for dehydration. After embedding in OCT medium (Tissue-Tek, Sakura Finetek, USA), the brains were sectioned at 15 or 20 μm using a freezing microtome. The brain sections were incubated with primary antibodies against tyrosine hydroxylase (TH, Millipore, MAB318, 1:400) and Iba1 (Proteintech, 10904, 1:200). For immunofluorescence staining, the sections were counterstained with DAPI. The morphometric analyses were performed using the ImageJ program (NIH, Bethesda, MA).
The number of TH+ neurons in the substantia nigra pars compacta (SNpc) was counted using an optical fractionator. Every three sections from a total of 9-12 sections per animal throughout the entire SNpc were counted in the lesion side under a ×20 objective view. The SNpc was identified by defined anatomic landmarks. TH+ cells with optimally visualized nuclei in one view were considered a valid count. To count the number of Iba1+ cells based on immunofluorescence staining, three different sections containing 1,200 µm × 1,400 µm standardized areas in the SNpc were selected randomly as described previously18. The total area of Iba1+ cells was calculated as the whole image fluorescence. After the mean background fluorescence was subtracted, the program ImageJ was used to identify and count the number of Iba1+ cells. The density of Iba1+ cells was determined based on the number of cells and area. Three sections were quantified per rat.
Quantification of TH expression in the striatum was performed using the Fiji image analysis system in ImageJ as previously described16, 18. The whole striatum region in both hemispheres was identified as a region of interest (ROI). TH immunoreactivity was measured in both hemispheres. After the sections were scanned, the images were converted to binary mode. All the rats shared the same threshold. The percentage was considered as the TH intensity of the lesion side compared with the contralateral side in the same section. A total of 4-5 striatal sections were assessed per rat.
Western Blot (WB)
The midbrain was separated, and its proteins were extracted as previously reported16, 18. Protein concentrations were estimated using a Nanodrop™ 8,000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The protein extracts were run on a 4-12% SDS‒PAGE gel and then transferred to PVDF membranes. The membranes were then incubated with primary antibodies against TH (Millipore, MAB318, 1:2000) or β-actin (Bioworld, AP0060, 1:5000). After washing with Tris-buffered saline containing 0.1% Tween-20, the membranes were incubated with HRP-conjugated secondary antibodies. An enhanced chemiluminescence system (ECL, Millipore) was used for visualization. The blots were assessed using the Image Lab program, and the results were normalized to the intensity of β-actin.
UPLC‒MS/MS analysis of CSF
After anesthesia with 1% pentobarbital sodium (40 mg/kg, i.p.), the rats were placed in a stereotaxic device. The head was position at a 135 degree angle to the body. After shaving and disinfection, a 2-cm longitudinal incision along the back midline in the occipital was generated to expose the foramen magnum. The insulin needle was used directly to puncture the dura mater to the cerebellomedullary cistern. Then, the CSF was drawn for UPLC‒MS/MS analysis.
Lipid extraction from CSF was performed as previously reported19. Fifty microliters of CSF sample was mixed with methyl tert-butyl ether and lipid standards to 1 ml. After incubation in an ultrasonic water bath with a frequency of 40 kHz and a power of 100 W for 5 min, 500 μl Milli-Q water was added, and the sample was centrifuged at 12,000 r/min at 10 °C. Then, 500 μl supernatant was collected and mixed with 100 μl mobile phase B. Then, 20 µl of each sample was added to a pool as a quality control. The detection system included a Shim-pack UFLC SHIMADZU CBM30A ultra-performance liquid chromatograph (UPLC) (SHIMADZU, JAPAN) coupled with a QTRAP5500 tandem mass spectrometer (MS/MS) (SCIEX, USA). Chromatographic separation was performed on a Thermo Accucore™ C30 column (2.1 mm × 100 mm, 2.6 μm). The column temperature was maintained at 45 °C. Mobile Phases A (60/40 acetonitrile/water) and B (10/90 acetonitrile/water) both contained 0.1% acetic acid and 10 mmol/L ammonium formate. The injection volume was 2.0 μL, and the flow rate was 0.35 mL/min. The gradient elution program was set as follows: 20% B at the beginning; 30-60% B between 2 to 9 min, 85-90% B between 9 to 15.5 min, 95% B between 15.5 to 17.5 min; 20% B for 17.5 to 20 min. The desolvation temperature was set as 500 °C. The MS voltage in positive and negative modes was set to 5500 V and -4500 V, respectively. The pressures of the curtain gas, gas II and ion source gas were 35 psi, 55 psi, and 45 psi, respectively. The ion pairs were detected by optimized collision energy and declustering potential in the triple quadrupole LC‒MS/MS.
Lipidomic analysis in CSF
The UPLC‒MS/MS analysis used a widely targeted metabolome method based on the Metware database (Metware Biotechnology Co., Ltd, Wuhan, Hubei, China)20. The retention time, ion pair information, and secondary spectrum data were used for qualitative analysis. The quantitative analysis was performed by the multiple reaction monitoring modes (MRM) of triple quadrupole mass spectrometry, and the data were analyzed using Analyst 1.6.3 software (AB SCIEX, Framingham, MA, USA) and MultiQuant (version 3.0, AB SCIEX). For peak area determinations, an individual area of the same metabolite was normalized to the integrated area of all peaks. All sample extracts were mixed as the quality control (QC), and one QC sample was inserted into every ten samples. Total ion flow diagrams of various QC samples were used to evaluate the repeatability of metabolite extraction and detection. The lipid species were classified as eicosanoids, triglycerides (TGs), phosphatidylcholine (PC), sphingomyelin (SM), free fatty acids (FFAs), diglyceride (DG), lysophosphatidylcholine (LPC), phosphatidylethanolamine (PE), ceramide (Cer), carnitine (CAR), monoglyceride (MG), cholesterol ester (CE), phosphatidylserine (PS), hexosylceramide (HexCer), sphingomyelin (SPH), butyric acid (BA), and phosphatidylglycerol (PG).
Unsupervised and supervised methods were used to evaluate whether the lipidomic signature in CSF could be determined after the induction of lesions using 6-OHDA. In the comparison between the early stage after 6-OHDA treatment and the control, the PCA results showed that PC-1 and PC-2 accounted for 26.2% and 16.5%, respectively (Supplement data 1 A). The PLS-DA results indicate that 22% of the total variance can be explained by component-1, and 12% of the total variance can be explained by component-2 (Supplement data 1 B). Using the OPLS-DA method, the orthogonal T score and T score were 22.4% and 9.3%, respectively (Supplement data 1 C). OPLS-DA data were subject to S-plot analysis (Supplement data 1 D). The Q2 intercept value of - 0.237 in the permutation test (n = 200) suggested that the OPLS-DA model had statistically effective quality and robustness (Supplement data 1 E). The lipids with the top 15 VIP values were identified and are shown in the VIP score plot (Supplement data 1 F).
The same methods were selected to compare the control with early- and late-stage 6-OHDA lesions. The PCA results showed that PC-1 and PC-2 accounted for 23.8% and 14.9%, respectively, in the comparison of the late stage with the control (Supplement data 2 A, Supplement data 3 A) and 30.7% and 18.5%, respectively, in the comparison between the early and late stages. The PLS-DA results showed that 15.1% of the total variance can be explained component-1 and 14.2% of the total variance can be explained by component-2 in the comparison of the late stage with the control. These values were 20.1% for component-1 and 19.3% for component-2 in the comparison of the two stages (Supplement data 2 B, 3 B). Using the OPLS-DA method, the orthogonal T score and T score were 22.1% and 6% as well as 25.7% and 13.7% in the two comparisons, respectively (Supplement data 2 C, 3 C). The S-plot and permutation test were performed, and the Q2 intercept values were - 0.301 and -0.275 in two comparisons (Supplement data 2 D, E; 3 D, E). The lipids with the top 15 VIP values are shown in the VIP score plot (Supplement data 2 F; 3 F).
Statistical analysis
The data are presented as the means ± standard errors of the means. Comparisons between two groups accepted the unpaired t test or Mann‒Whitney test after the normality test. Comparisons among multiple groups were analyzed by one-way ANOVA followed by the LSD post hoc test. Data analyses were performed using SPSS 16.0 (IBM, USA), and p < 0.05 indicated a statistically significant difference.
In lipidomic analysis, SIMCA version 15 software (Umetrics, Umeå, Sweden) and MetaboAnalyst software (Version 5.0, https://www.metaboanalyst.ca/) were used for the principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), and orthogonal partial least-squares discrimination analysis (OPLS-DA). The volcano plot, box plot and heatmap were drawn using GraphPad Prism (GraphPad Prism® Software version 8.0.2 for Windows; La Jolla, CA) and Origin software (Version 2022, OriginLab Inc., USA). The metabolites with variable importance in the projection (VIP) ≥ 1 and absolute log2Fold change ≤ 1 were used for cluster analysis with Origin software (Version 2022, OriginLab Inc., USA). Violin plots were generated using Hiplot (https://hiplot.com.cn). Pearson correlations were analyzed by Origin software (Version 2022, OriginLab Inc., USA). In the Pearson correlation and linear mixed-effects model, the significance of each variable was assessed at a level of 0.05.