Patients
Twenty-nine patients with MB (21 male and 8 female) were recruited at the Ophthalmology Department of the First Affiliated Hospital of Nanchang University. These subjects satisfied the following criteria: 1) blind in one eye; 2) contralateral eye is normal without cataract, optic neuritis, or other eye diseases.
In addition, 29 healthy controls (21 male and 8 female) were recruited and the two groups were similar in gender balance (P > 0.99), age (P = 0.792), and weight (P = 0.881). Control subjects were included if they satisfied the following criteria: 1) normal naked eye or normal corrected vision ; 2) no neurological diseases; 3) no mental disorder; 4) able to have an MRI scan (for example, they did not have pacemaker or implanted metal device).
The research was authorized by the Human Research Ethics Committee of the First Affiliated Hospital of Nanchang University. Each participant understood the aim, methods and possible risks of the research, and signed a declaration of informed consent, and all the experiments were performed in accordance with the Declaration of Helsinki.
MRI data collection
The Trio 3-Tesla MR scanner (Siemens, Munich, Germany) was used. Before scanning, each participant was asked to relax, close their eyes, and minimize movement[29]. To obtain functional data, a 3D metamorphic gradient echo pulse sequence was used. The following parameters were used for a 176-image scan: acquisition matrix 256×256; field of view 250×250 mm; echo time 2.26 ms; repetition time 1,900 ms; thickness 1.0 mm; gap 0.5 mm; flip angle 9°. For a 240-image scan, parameters were as follows: acquisition matrix 64×64; field of view 220×220 mm; thickness 4.0 mm; gap 1.2 mm; repetition time 2,000 ms; echo time 30 ms; flip angle. 90°, 29 axial.
fMRI processing
MRIcro software (Nottingham University, Nottingham, UK) was used to sort the data, and to identify and exclude incomplete or flawed data. Remaining data were processed, including space standardization, head movement correction, slice time, and digital image format conversion using DPARSFA (http://rfmri.org/DPARSF). Linear regression was used to eliminate the influence of factors such as signals originating from white matter.
Because excessive head movement may have a significant impact on the fMRI sequence, participants with head movements > 3mm and the data were excluded. Due to inter-individual variations in brain size and structure, each brain image was standardized[30]. We used regions of interests (ROI) of the central white matter region to deal with irrelevant variables[31].
fMRI data were processed using the PerAF method, a relatively reliable and direct measurement of brain activity. First, the average BOLD signal value was calculated, then the signal strength at a range of time points was normalized to this value. This process resulted in an amplitude at each time point as a percentage of the average across the time series, and a signal change percentage similarity index, referred to as percentage amplitude fluctuation (PerAF). The formula used to calculate the PerAF value of a single voxel is as follows:
PerAF = 1n∑i = 1n|Xi − µµ|×100% (1)
µ = 1n∑i = 1nXi (2)
Where Xi represents the signal strength, n is the total number of time points, and µ is the mean value of the time series[21].
Correlation analysis
We obtained the anxiety scores (AS) and depression scores (DS) of MB patients by doing the Hospital Anxiety and Depression Score (HADS). We looked for correlations between each score and the PerAF values of the following brain regions: Frontal_Sup_Orb_L/ Frontal_Inf_Orb_L, and Frontal_Inf_Oper_L using Pearson’s correlation analysis (P<0.05 was considered significant). GraphPad Prism 8.0 software was used to plot linear correlations.
Statistical analysis
For between-group comparisons, SPSS software, version 20.0 (IBM Corp., Armonk, NY, USA) was used to conduct independent sample t tests, and P<0.05 was considered significant. The REST software was used to conduct independent sample t tests comparing PerAF values between the two groups. Gaussian random field theory was used for multiple comparison correction, and the voxel level threshold was p<0.001. AlphaSim, part of the REST toolbox, was used for correction, the cluster size was set at >49 voxels, and the level was p<0.05. Receiver operating characteristic (ROC) curves were used to compare the average PerAF values of the relevant brain areas between MB and HC groups and to obtain estimates of diagnostic accuracy based on the area under the curve (AUC). As explained above, Pearson’s correlation was used to evaluate the relationship between PerAF and anxiety/depression scores. All averaged data are presented in the form of mean ± standard deviation. The regions were defined using automatic anatomic labeling based on the Montreal Neurological Institute data set.