Animal Model
Five-days-old male Sprague Dawley (SD) rats (Sibeifu, Beijing, China, license NO. SCXK-Beijing-2016-0002) were purchased and housed in animal cages (35 cm x 25 cm x 20 cm; before ablactation: a litter/cage; from ablactation at third week to adult: 3 rats/cages) under the following conditions: room temperature of 23 ± 2℃, relative humidity of 60% ± 5% artificially lighted from 7:00 a.m. to 7:00 p.m. with light intensity similar to daylight. The rats had free access to conventional rodent lab chow and water. The cages and bedding were changed every day. The person in charge of the animal’s welfare carried out the experiments. We used only male rats for the study to avoid possible confounding effects of gender due to differences in hormones. This study was conducted in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The experimental protocol was approved by the Animal Care and Use Committee in the Institute of Acupuncture and Moxibustion, China Academy of Chinese Medical Sciences (D2019-11-14-1). All efforts were made to minimize suffering.
The modeling, intervention and experimental protocol are shown in Figure 1. Except for the control group, all animals were gavaged with iodoacetamide (IA, 0.2 mL 0.1% IA in 2% sucrose solution) at the age of 10 days daily for 6 days21. This process has been utilized in previous studies to bulit the model of FD rats that also induced depression-like behaviors in a rodent model of FD.22, 23 Rats in the control group did not received any intervention.
Intervention
Transcutaneous auricular vagus nerve stimulation (taVNS)
The rats in the taVNS group were placed in a fixator, permitting relatively free movement. After 30 min of adaptive fixation, two opposite magnetic electrodes (+/-) were placed on both sides (inside and outside) of bilateral auricula concha so that electrical current could conduct through the tissue. Saline was applied between the electrode and the skin to improve electrical current conduction. At age of 9 weeks, taVNS was performed 30 min once daily for 14 days using the following parameters: 100 Hz, 0.5ms, 0.1s-on, 0.4s-off and an intensity of 0.5 mA via a watch-size electrical stimulator (Transcutaneous Electrical Acustimulation, SNM-FDCM01, Medkinetic, Ningbo, China).
Behavioral Tests
Open field test (OFT)
Rats were acclimatized to the laboratory environment one hour before the test. A black square opaque box (100 cm x 100 cm x 40 cm) was utilized for the OFT test. The box was placed in the center of the room with a camera fixed above the box on the ceiling. Four light sources were used to eliminate shadows. Each rat was placed in the center of the box, and the camera recorded animal movement. The horizontal and vertical scores were recorded for 5 minutes. The horizontal score was the number of horizontal crossing grids of the rat (all four paws enter the square to be counted). The vertical score was the number of vertical standing (the two front paws of the rat were off the cage bottom or against the side wall of the open box). The room was kept quiet during the recording, and the animals had no visual contact with the investigator. The chamber was thoroughly cleaned with 75% ethanol and water to eradicate any remnant animal scent.24
Forced swimming test (FST)
Rats were placed individually in a transparent plexiglass cylinder (60 cm in height, 30 cm in diameter) filled with fresh water (24 -26℃) up to a height of 45 cm from the bottom. The immobility time of each rat, a total of 5 minutes, was individually recorded.25
MRI Experiments
All MRI experiments were conducted using a 7.0 T animal MRI scanner (Bruker ClinScan Germany). During the MRI scans, the animals were anesthetized with isoflurane (1.0–1.5%, 1.0 L/ min). The rodent’s body temperature was maintained using a warm water circulation system. BOLD rs-fMRI images were acquired with the following parameters: field of view (FOV) = 25×20 mm2, thickness = 1.0 mm, slice gap=0.2 mm, repetition time (TR)/echo time (TE) = 2691.4 ms/27.6 ms, flip angle = 90°. 26
The fMRI image data pre-processing was performed using Statistical parametric mapping 12 (SPM12) software on the Matlab R2022b. The data were converted to Neuroimaging informatics technology initiative (NIFIT) format for recognition and processing within the software. The image proctor data from the first 10 time points of each subject were removed to exclude the effect of the signal onset instability of the magnetic field and to achieve the adaptation of the experimental animals to the magnetic field. Slice timing: interpolation was used to make the results of each scan converge to the actual results at a certain time. Realign: removing data where the head moved more than 3 mm in each direction (X, Y, Z) or rotated more than 3° in any angle to reduce the effect of head movement on the image. Spatial normalization: the brain images of experimental rats were aligned to the standard brain to exclude the differences in anatomical structures of individual brain tissues and overcome the influence of different brain structures on the data. Smooth: spatial smoothing was performed by using a 6mm half-orbit full-width (FWHM) Gaussian smoothing kernel to blur the differences between the brain structures of different subjects, improve the signal-to-noise ratio of the images, and reduce the remaining inter-individual differences after normalization, to make the images more in line with the assumption of Gaussian distribution. Detrend: used to reduce the effects on the data due to the increased temperature of the MRI machine and the subjects' adaptation to the environment. Low-frequency filtering: extracting their low-frequency signals in the frequency range of 0.01-0.08 Hz to remove the effects of high-frequency respiratory heartbeats and low-frequency noise. The ALFF and ReHo were calculated.27
Statistical Analysis
SPSS 26.0 software was used for statistical analysis. A Shapiro-Wilk test was used to examine the normality of the data. If the data fit a normal distribution, one-way ANOVA was used for comparison among multiple groups. The Fisher's least significant difference (LSD) post hoc test was used for data with a homogeneity of variance. Tamhane's T2 (M) test was used after ANOVA if the data were not by homogeneity of variance. If the data did not conform to normality, the Kruskal-Wallis H test was applied. Analyses and graphing were conducted by GraphPad Prism 8.0. Differences were considered statistically significant at P<0.05.
Differences of ALFF and ReHo among the three groups were used One-way ANOVA with F > 0 representing a significant brain region of difference among the three groups combined with significance criteria of voxel P<0.005 and cluster P<0.05.
To examine the associations between the evaluated behaviors and the cerebral alterations of taVNS for FD, we performed Pearson (data fit a normal distribution) or Spearman (data did not have a normal distribution) correlation analyses between the behavioral test indices and ALLF/ReHo values in the brain regions that exhibited significant differences between the control and IA groups, as well as between the IA and taVNS groups. Differences were considered statistically significant at P<0.05.