the nodes can be brain regions or voxels, or even
neurons, while edges can be identified by anatomical
connectivity, functional connectivity depending on
the characteristics of the data set. Functional
connectivity indicates a correlation in time between
the neurophysiological activity of two points. Based
on graph theory, Chen found that the whole brain
network of OSA patients exhibited decreased
smallworld topological property (Chen, 2018).
However, Huang found that the smallworld property
of the OSA brain network was not significantly
different to HC’s, but OSA performed significantly
lower in Cp and higher in Lp (Huang, 2019). Hence,
it follows that the functional network organization of
OSA patients are not clear yet.
In this paper, we intend to explore the impact of
OSA disease on patients' functional brain networks
from a complex network perspective.
2 MATERIAL METHOD
2.1 Data Acquisition
The data of this experiment were divided into two
groups, the OSA patient group and the healthy
controls group (HC), which included 24 OSA patients
with typical indications of OSA disease and met the
diagnostic criteria for the disease in the relevant draft;
21 healthy volunteers were also recruited as the
control group.
In this study, MRI data were acquired using a GE
Signa HDx magnetic resonance scanner with a field
strength of 3.0 T. Scans and resting-state fMRI data
were acquired using Gradient Recalled Echo (GRE)
single excitation planar echo imaging (EPI) sequence
with the following parameters: TR=2000ms,
TE=30ms, FA=90°, FOV was 240*240mm², a 64*64
matrix was used, thickness =3mm, the layer spacing
(gap) was 1mm, a total of 38 layers were divided, and
180 time points were acquired in each volume.
2.2 Data Pre-Processing
The DPARSF (Data Processing Assistant for Resting-
State fMRI) software based on MATLAB was used
to preprocess the functional MRI data of both OSA
and HC related to the following steps.
The images of the first 10 time points were
excluded to avoid the potential noise and instability.
Slice timing was applied to correct this time
difference. Head movement correction was used to
avoid a slight head movement. The unified standard
spatial EPI template with a voxel size of 3×3×3mm3
was used for transformation in order to facilitate the
later study. The whole brain average signal,
cerebrospinal fluid and white matter signal and
motion signal were regressed out. The linear drift was
removed. Then a band-pass filtering (0.01-0.08Hz)
was used to remove the interference of low frequency
and high frequency signals.
2.3 Construction of Brain Networks
and Calculation of Topological
Properties
Brain networks were constructed and topological
properties were calculated. Based on the human
Brainnetome Atlas template (Fan, 2016), the
cerebrum of each subject was divided into 246 brain
regions (nodes). The Pearson correlation coefficient
of the average time series of each two brain regions
was calculated as the functional connectivity (edges),
then followed by the Fisher r-to-z transformation to
normalize it. The z-scored 246×246 correlation
coefficient matrix was obtained.
The correlation matrices were binarized by a pre-
selected value of sparsity K (0.05 ≤ K ≤ 0.5). For a
specific K, we got an undirected binarized network
for each subject, then we applied graph theory to
calculate the topological properties of each brain
network. The global network properties contain
clustering coefficient (C
p
), characteristic path length
(L
p
), normalized clustering coefficient (γ),
normalized characteristic path length (λ), and small-
world property (σ). The nodal properties contain the
betweenness centrality (BC) and degree centrality
(DC) of each brain regions.
2.4 Statistical Analysis
Two-samples t-test was performed for each parameter
corresponding to the two groups of subjects. p<0.05
is considered to be statistically different.
3 RESULTS AND DISCUSSIONS
3.1 Global Properties
The area under curve (AUC) of each global parameter
(C
p
, L
p
, γ,λ, and σ) in OSA group and HC group is
shown in Figure 1(a). The Cp, Lp,λ of OSA patients
were significantly lower than healthy controls (p <
0.05), while theγand σ of OSA patients were
significantly higher than healthy control (p < 0.05).