Alterations of the default mode network and cognitive impairments in patients with chronic obstructive pulmonary disease.

Background and objectives: Cognitive impairment is a common extrapulmonary comorbidity in COPD patients. The default mode network (DMN) plays a critical role in maintaining the normal activities of humans, and its function can be evaluated by resting state functional magnetic resonance imaging. The aim of this study was to investigate the correlations between cognition and function changes of the DMN in COPD patients.
Methods: One hundred and thirteen eligible participants including 30 control subjects and 83 COPD patients matched for demographic characteristics were recruited. All participants performed cognitive function tests and underwent resting state functional magnetic resonance imaging.
Results: The total cognitive function scores of COPD patients were significantly different from those of control subjects (P<0.05) and worsened with the degree of airflow obstruction. The activated brain regions in the DMN of COPD patients were less than those of normal controls. Six activated brain regions in the DMN were found to develop significantly different functional connectivity (FC) values among the subjects. Meanwhile, the FC values of the left posterior cingulate cortex and left hippocampus correlated well with cognitive functions and pulmonary function.
Conclusion: COPD patients have cognitive impairments that correlate well with disease severity. FC changes in activated brain regions in the DMN may predict cognitive impairment, and the left posterior cingulate cortex and left hippocampus may be important brain regions related to cognitive impairment in COPD patients.

Introduction

Chronic obstructive pulmonary disease is characterized by irreversible airflow limitation and is often accompanied by many extrapulmonary comorbidities. Dysfunction of the central nervous system has been verified as a distinct extrapulmonary comorbidity of COPD patients in numerous published studies.1,2 Meanwhile, studies have demonstrated that COPD patients with cognitive impairment have lower levels of daily activity and an increased risk of acute exacerbation,3–5 which may contribute to the overall disease severity in individuals.

Neuroimaging studies can display changes in brain metabolism, structure, and functions in COPD patients. Using magnetic resonance spectroscopy, Shim et al found a significant change in brain metabolism in severe COPD patients.6 With single photon emission computed tomography, Ortapamuk and Naldoken showed significantly lower brain perfusion in frontal and parietal lobes in COPD patients.7 Dodd et al displayed decreased white matter integrity and disturbed functional activation of gray matter in COPD patients.8

Resting state functional magnetic resonance imaging (rs-fMRI) is an effective way to investigate brain functions characterized by spatially coherent spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, and it is performed under resting and waking states without any thinking activities or cognitive tasks. The brain is quite active under this relatively resting state.9 Functional connectivity (FC) analyses of rs-fMRI can demonstrate the dynamic interactions between different brain regions; it can be defined as the temporal correlation of neurophysiologic index measured in different brain regions.10 Currently, the rs-fMRI is widely used in the areas of neurological and psychiatric diseases.11–13 Most of the brain’s energy budget is consumed by the spontaneous on-going activity when the brain is at rest, and this spontaneous activity is most prominent in a distributed cortical network commonly called the default mode network (DMN), which plays a critical role in maintaining the metabolic activity of brain tissue, as well as in the process of spatial orientation, internal monitoring, and memory in the resting state.14

Previous studies indicated that cognitive function of COPD patients is closely related to pulmonary function and varied in different COPD grades.15 Our previous study demonstrated that a decreased volume of bilateral hippocampal tissue was related to cognitive impairment in COPD patients.10,16 However, little is known about the correlations between brain function changes (especially the DMN) and cognitive impairment in COPD patients. So, in this study, we evaluated cognition and function changes in the DMN in both control subjects and COPD patients. We hypothesized that COPD patients of different disease grades may experience varied cognitive impairment, as well as brain function changes. We aimed to investigate the correlation of brain function with cognitive impairment in COPD patients.

Methods

Study subjects

This study involved 117 participants, including 30 control subjects and 87 COPD patients. The control subjects were all volunteers, and the COPD patients were recruited from the pulmonary department of the First Affiliated Hospital of Anhui Medical University. The diagnoses and classification of COPD followed the Global Initiative for Chronic Obstructive Lung Disease (GOLD) guidelines from 2013.17 COPD patients were all in a stable state (with no exacerbations during the past 8 weeks). Exclusion criteria were as follows: 1) comorbidities such as diabetes, liver failure, cardiovascular disease, neurological disorders, and malignant tumor; 2) sleep apnea syndrome; 3) mental disorders; 4) history of alcoholism and/or drugs; 5) dementia; 6) visual and/or auditory impairment; 7) education levels below 5 years; 8) receiving medicines that may cause cognitive impairment. All participants were required to be right-handed.18 All of them received a complete physical examination by a respiratory physician, and their cognitive status was evaluated by a neuropsychologist. All subjects provided written consent to participate and were informed of the possible risks of the study. The study was performed according to the Declaration of Helsinki and approved by the ethics committee of the First Affiliated Hospital of Anhui Medical University (NO Quick-PJ-2016-12-11).

Arterial blood gas analysis

Arterial blood gas analysis was performed while patients breathed room air in the supine position. Power of hydrogen (PH), arterial partial pressure of oxygen (PaO2), arterial partial pressure of carbon dioxide (PaCO2), and blood oxygen saturation (SaO2) were evaluated with a Stat Profile Critical Care Xpress (Nova Biomedical, Waltham, MA, USA).

Pulmonary function tests

A standardized pulmonary function test was performed with a dry spirometer device (Erich Jaeger GmbH, Hoechberg, Germany) at 15 minutes after inhaling salbutamol 400 µg (Ventolin; GlaxoSmithKline, London, UK); forced vital capacity (FVC), forced expiratory volume in 1 second (FEV1), and the FEV1/FVC ratio were recorded. The subjects with FEV1/FVC<0.7 and FEV1≥80% predicted were classified as mild COPD, 50%≤FEV1≤80% predicted were classified as moderate COPD, 30%≤FEV1≤50% predicted were classified as severe COPD.

Cognitive function test

Cognitive function was assessed using the Montreal Cognitive Assessment (MoCA).19 The MoCA is a rapid screening instrument for mild cognitive impairment. The test assesses multiple cognitive functions including attention, concentration, executive functions, memory, language visuoconstructional skills, conceptual thinking, orientation, and calculation. The highest score is 30, and a score >26 means normal. All subjects completed the cognitive function tests according to Liesker et al.20

Image acquisition

Imaging data were acquired using a 3.0-T MRI scanner (Signa HD 3.0T; General Electric, Boston, MA, USA) with a 16-channel head coil in the Department of Radiology, First Affiliated Hospital of Anhui Medical University. Each participant was instructed to remain stationary, close his or her eyes and not think of anything while scanning. Subjects lay supine in the scanner and wore earplugs to minimize machine noise. Head movement was minimized by inserting soft foam padding between the subject’s head and the head coil. To capture the whole brain structure and eliminate organic brain disease, we acquired axis T2-weighted imaging paralleled with the anterior–posterior commissure. The parameters were as follows: slice thickness =5 mm, gap =1.5 mm, 19 slices, repetition time (TR) =5,087 ms, echo time (TE) =120 ms, field of view (FOV) =24×19 mm, matrix =384×256. At the same position, the resting state fMRI was performed with axis spin echo-echo planar imaging scanning with the following parameters: slice thickness =3.6 mm, gap =0.4 mm, 38 slices, TR =2,000 ms, flip angle =90°, TE =35 ms, FOV =240×240 mm, matrix =64×64, and 240 time points in total. A three-dimensional anatomic structure image using speed phase gradient reverse recovery sequence sagittal scanning was performed with the following parameters: slice thickness =1.2 mm, no gap, 166 slices, TR =8.5 ms, flip angle =12°, TE =3.3 ms, FOV =240×240 mm, matrix =256×256, number of excitations =2. All images were transferred to the ADM4.4 imaging workstation for further analysis.

Image analysis

Data pre-processing

The image format was converted from DICOM to NIFTI in advance; the rs-fMRI data preprocessing comprised slice-timing correction, realignment, spatial normalization, and spatial smoothing as implemented in MATLAB 7.12.0 (R2011a; The MathWorks, Natick, MA, USA), Statistical Parametric Mapping (SPM12; Wellcome Department of Cognitive Neurology, London, UK), and Data Processing Assistant for Resting State fMRI (DPARSF3.2). Firstly, all subjects required time for aligning and motion correction; we computed frame-wise head motion for each subject, and one-way analysis of variance (ANOVA) indicated no significant difference among the four groups; head motion translation values >2 mm in any direction of the x, y, and z parameters or rotation values more than 2° of any angle were excluded. Secondly, the first 10 seconds of data were discarded for scanner calibration and subject adaptation. Then we performed the coregistration between functional and anatomical images and normalization to the standard Montreal Neurological Institute (MNI) 152 brain atlas; figure data were resampled based on a volume of 3×3×3 mm3 as a unit. Spatial smoothing was used with a 4-mm full-width at half maximum Gaussian smoothing kernel.

To ensure that each rs-fMRI data point was the best representative of spontaneous neural activities, the effects of physiological sources were minimized by regressing out estimated predictors, which included white matter, cerebral spinal fluid, the global signal, and six head motion parameters. The waveform of each voxel was finally used to remove the linear trends of time courses and for temporal band-pass filtering (0.01–0.08 Hz) to reduce low-frequency drift and high-frequency physiological noise. The individual data were used for further FC analyses.

Functional connectivity analyses

In the FC analyses, the DMN is anchored in the PCC, and left posterior cingulate cortex (L-PCC) is an important functional region in the DMN.21,22 We chose L-PCC as the seed in the Anatomical Automatic Labeling (AAL) model; other brain regions of the DMN were set as regions of interest (ROIs), and the parameters of MNI coordinates were calculated and expressed as (x, y, z). The average values of the BOLD time series of L-PCC in the AAL model were then acquired, and we computed the Pearson’s correlation coefficients between the average BOLD time series of L-PCC and those from the whole brain. The correlation coefficients were transformed to better fit a normal distribution using Fisher’s Z transformation, and thus, the FC values of the different brain regions in the DMN were obtained for each participant.

Statistical analysis

Measurement data were expressed as the means ± standard deviation, and enumeration data were expressed as a percentage. SPSS Statistics version 16.0 (SPSS, Inc., Chicago, IL, USA) was used for statistical analysis. Group comparisons were tested using one-way ANOVA with a Student–Newman–Keuls post hoc correction for multiple comparisons or the chi-square test.

Voxel-wise k-sample t-tests with one covariate (age) were performed to assess the differences in the FC of different brain regions in the DMN; the Monte Carlo method correction was performed using a whole brain mask of MNI template. A cluster >10 voxels with P-value <0.001 was considered as significant. Pearson’s correlation analysis was used to evaluate the relationships between the FC values of different brain regions in the DMN and PaO2, PaCO2, MoCA, and FEV1. P<0.05 was considered as significant.

Results

Demographic characteristics of the subjects

Four COPD patients were excluded for head movements, and 113 eligible participants including 30 control subjects, 29 mild COPD patients, 30 moderate COPD patients, and 24 severe COPD patients were finally recruited. Among them, 55 COPD patients received blood gas analysis. The 113 subjects were classified into four groups as shown in Table 1, where the four groups were statistically similar with respect to age, sex, BMI, smoking status, and education level (P>0.05). However, FEV1, PaO2, SaO2 and PaCO2 were significantly different among the four groups (P<0.05).

Table 1

Demographic characteristics of participants

Variables	Control group (n=30)	Mild COPD (n=29)	Moderate COPD (n=30)	Severe COPD (n=24)	F/χ2-value	P-value
Age (years)	57.43±8.66	59.38±8.27	60.90±10.63	63.42±9.53	1.96	0.12a
Sex (M/F)	(21/9)	(13/16)	(16/14)	(17/7)	7.01	0.07b
Education (years)	9.20±2.29	9.86±2.66	9.60±2.06	9.25±2.80	0.46	0.71a
BMI (kg/m2)	23.9±3.02	24.20±4.40	23.25±4.20	22.16±3.41	1.48	0.22a
Smoking history (Y/N)	(19/11)	(12/17)	(15/15)	(16/8)	4.63	0.20b
Smoking index (packs/year)	22.79±6.91	21.08±9.06	24.27±7.04	25.00±5.05	0.84	0.47a
FEV1 (%)	106.62±22.86	85.47±4.17*	66.23±8.21*,#	38.18±7.39*,#,+	132.87	0.001a
Blood gas analysisc
PH	7.40±0.01	7.40±0.02	7.40±0.01	7.40±0.02	0.104	0.958a
PaO2 (mmHg)	115.65±8.92	101.27±8.64*	76.90±5.79*,#	66.97±3.75*,#,+	163.19	0.001a
SaO2 (%)	38.88±2.73	38.81±3.99	39.19±2.54	40.05±4.21+	5.15	0.003a
PaCO2 (mmHg)	98.22±0.51	97.09±0.79*	94.59±0.97*,#	92.04±1.38*,#,+	120.29	0.001a

Notes: Data are presented as the mean ± standard deviation, except where indicated otherwise.

ANOVA;

Chi-square test;

Control group (n=15), mild COPD (n=16), moderate COPD (n=22), severe COPD (n=17).

P<0.05 versus control group;

P<0.05 versus mild COPD;

P<0.05 versus moderate COPD.

Abbreviations: ANOVA, analysis of variance; BMI, body mass index; F, female; FEV1, forced expiratory volume in 1 second; M, male; N, no; PaCO2, arterial partial pressure of carbon dioxide; PaO2, arterial partial pressure of oxygen; PH, power of hydrogen; SaO2, oxygen saturation; Y, yes.

Cognitive function

The total MoCA scores were significantly different among the four groups (P<0.05, shown in Table 2). Within single score comparisons, the scores of executive function, attention, delayed recall, abstraction, and orientation were significantly different among the four groups (P<0.05).

Table 2

Cognitive function test scores of participants

Variables (score)	Control group	Mild COPD	Moderate COPD	Severe COPD	F-value	P-value
MoCA	26.17±1.29	24.34±2.39*	22.20±3.44*,#	21.75±2.83*,#	17.44	0.001
Visuospatial executive	4.53±0.57	4.03±0.86	3.40±1.04*	3.12±0.99*,#	14.43	0.001
Naming function	3.00±0.00	2.90±0.31	2.77±0.50	2.83±0.38	2.38	0.07
Attention	5.40±0.56	5.24±0.69	4.73±0.91*	4.88±0.74*	5.19	0.002
Verbal	1.90±0.48	1.79±0.82	1.60±0.62	1.46±0.78	2.26	0.08
Abstraction	1.80±0.41	1.45±0.78	1.33±0.76*	1.46±0.59	2.83	0.04
Delay recall	3.33±0.96	2.90±1.24	2.40±1.25*	2.58±0.83*	4.10	0.008
Orientation	5.97±0.18	5.90±0.41	5.77±0.43	5.46±0.59*,#	7.51	0.001

Notes: Data are presented as the mean ± standard deviation, except where indicated otherwise.

P<0.05, versus control group;

P<0.05, versus mild COPD.

Abbreviation: MoCA, Montreal Cognitive Assessment.

Functional connectivity

The activated brain regions in the DMN of the four groups are displayed in Figure 1. We found that the activated brain regions in the DMN of the severe COPD patients were obviously fewer than in the other three groups. The post hoc multiple comparisons are shown in Table 3.

Figure 1

The cross-sectional MRI images of brain activation regions (yellow areas) related to DMN (the left posterior cingulate cortex as the seed) among four groups.

Notes: (A) Control subjects; (B) mild COPD group; (C) moderate COPD group; (D) severe COPD group.

Abbreviations: DMN, default mode network; MRI, magnetic resonance imaging.

Table 3

The post hoc multiple comparisons of FC values of DMN among four groups

Brain region	Voxel	Z-value	P-value	MNI coordinates			Brodmann area
				X	Y	Z
The FC values decreased gradually
Control group > mild COPD
L-inferior temporal gyrus	45	3.62	0.001	−42	−9	−42	20
R-hippocampus	16	4.07	0.001	24	−6	−21	–
R-middle occipital gyrus	23	3.35	0.001	30	−90	15	18
R-fusiform gyrus	16	2.91	0.001	33	−36	−24	37
Control group > moderate COPD
L-middle temporal gyrus	20	3.41	0.001	−54	0	−30	21
Control group > severe COPD
R-medial orbitofrontal cortex	293	4.78	0.001	3	51	−12	11
R-middle orbitofrontal cortex	68	3.82	0.001	24	66	−9	11
R-supper medial frontal gyrus	16	3.32	0.001	6	51	36	9
R-opercular part of inferior frontal gyrus	16	3.14	0.001	51	24	36	44
L-anterior cingulate cortex	293	4.53	0.001	−3	45	3	10
L-posterior cingulate cortex	111	4.54	0.001	−3	−48	27	30
R-middle cingulum cortex	37	3.46	0.001	3	−15	39	23
L-middle temporal gyrus	35	3.35	0.001	−66	−15	−9	21
L-inferior temporal gyrus	240	4.20	0.001	−48	−21	−30	20
L-hippocampus	240	4.72	0.001	−15	−9	−21	–
L-precuneus gyrus	50	3.68	0.001	−3	−45	57	5
R-rectus gyrus	293	4.29	0.001	6	39	−18	11
R-posterior central gyrus	87	3.85	0.001	39	−24	54	4
Mild COPD > moderate COPD
None
Mild COPD > severe COPD
L-medial orbitofrontal cortex	95	3.60	0.001	6	42	−12	11
R-middle frontal gyrus	16	3.23	0.001	33	24	39	46
L-superior frontal gyrus	32	3.19	0.001	−21	24	42	0
L-opercularis frontalis inferioris gyrus	33	2.91	0.001	−42	15	33	44
L-anterior cingulate cortex	95	3.71	0.001	−3	45	3	10
L-posterior cingulate cortex	120	4.32	0.001	−3	−45	27	26
L-precuneus gyrus	120	3.68	0.001	−6	−54	21	23
L-inferior parietal gyrus	41	3.54	0.001	−51	−27	39	2
L-posterior central gyrus	41	3.56	0.001	−54	−18	42	3
Moderate COPD > severe COPD
R-medial orbitofrontal cortex	51	4.40	0.001	6	45	−12	11
L-superior frontal gyrus	19	3.41	0.001	−24	15	48	8
L-anterior cingulate cortex	26	4.02	0.001	−6	30	24	32
R-inferior temporal gyrus	12	3.97	0.001	48	−33	−24	20
R-superior parietal gyrus	24	3.73	0.001	15	−54	63	5
R-rectus gyrus	51	4.00	0.001	3	36	−15	11
L-precuneus gyrus	12	3.95	0.001	−9	−48	63	5
L-fusiform gyrus	17	4.37	0.001	−24	−39	−15	30
R-fusiform gyrus	34	4.25	0.001	30	−9	−33	36
L-cuneus gyrus	11	3.77	0.001	0	−81	30	18
The FC values increased gradually
Severe COPD > moderate COPD
R-rolandic operculum	10	4.34	0.001	42	0	15	48
Severe COPD > mild COPD
None
Severe COPD > control group
R-rolandic operculum	15	4.45	0.001	42	−3	15	48
Moderate COPD > mild COPD
R-orbit inferior frontal gyrus	16	3.91	0.001	30	33	−15	47
R-supramarginal gyrus	17	3.60	0.001	63	−33	42	40
R-parahippocampus	10	4.15	0.001	27	−12	−36	–
Moderate COPD > control group
L-triangle of inferior frontal gyrus	18	3.96	0.001	−42	39	9	45
R-middle frontal gyrus	13	3.81	0.001	27	9	48	6
L-supramarginal gyrus	39	4.08	0.001	−57	−39	30	48
L-inferior parietal gyrus	11	3.94	0.001	−48	−39	42	40
L-lingual gyrus	23	4.17	0.001	−15	−78	−51	18
Mild COPD > control group
None

Abbreviations: DMN, default mode network; FC, functional connectivity; L, left; MNI, Montreal Neurological Institute; R, right.

Figure 2 shows the activated brain regions in the DMN that have significantly different FC values among the four groups from three different planes, as well as the three-dimensional images. The results are also presented by parameters of the MNI coordinates and are shown in Table 4. The meaningfully activated brain regions in the DMN were L-PCC, left anterior cingulated cortex (L-ACC), left supra-marginal gyrus (L-SMG), left hippocampus (L-HP), right parahippocampus (R-PHP), and right medial orbitofrontal cortex (R-MOFC).

Figure 2

Brain regions of DMN that had significantly different FC values among the four groups.

Notes: (A) Horizontal axial; (B) coronal; (C) sagittal; (D) the three-dimensional brain map of significant group effect (red areas). The statistic t-score bar is shown in the maps.

Abbreviations: DMN, default mode network; FC, functional connectivity; L, left; R, right.

Table 4

Brain regions in DMN that have significantly different FC

Brain region	Voxel	Z-value	P-value	MNI coordinates			Brodmann area
				X	Y	Z
L-hippocampus	10	4.18	0.001	−15	−9	−21	–
L-anterior cingulate cortex	10	4.11	0.001	−3	−48	27	30
R-medial orbitofrontal cortex	54	4.08	0.001	3	51	−12	11
L-posterior cingulate cortex	54	3.90	0.001	−3	45	3	10
R-parahippocampus	12	4.01	0.001	27	−12	−36	36
L-supramarginal cortex	15	3.45	0.001	−60	−39	30	48

Abbreviations: DMN, default mode network; FC, functional connectivity; L, left; MNI, Montreal Neurological Institute; R, right.

Group comparisons of FC analyses (shown in Figure 3) revealed that compared with the other three groups, the severe COPD patients had significantly lower FC values in L-ACC, R-MOFC, and L-HP (P<0.01), while for L-SMG, both moderate and severe COPD patients had significantly different FC values when compared with the mild COPD patients and the control group (all P<0.05). For L-PCC, the severe COPD patients had significantly lower FC values than mild COPD patients and the control group (P<0.01). For R-PHP, the severe COPD patients had significantly lower FC values than moderate COPD patients and the control group (P<0.01), and the mild COPD patients had significantly lower FC values compared with moderate COPD patients and the control group (P<0.01).

Figure 3

Comparisons of FC values of brain regions in DMN among the four groups.

Notes: For L-ACC, R-MOFC, and L-HP, severe COPD patients had significantly lower FC than the other three groups (P<0.01), and mild COPD patients had a lower FC of L-HP compared with the control group (P<0.05). Whereas, for L-SMG, both moderate and severe COPD patients had significantly different FC values when compared with the mild COPD patients and the control group (all P<0.05). For L-PCC, severe COPD patients had significantly lower FC values than mild COPD patients and the control group (P<0.01). For R-PHP, severe COPD patients had significantly lower FC values than moderate COPD patients and the control group (P<0.01), and mild COPD patients had significantly different FC values compared with moderate COPD patients and the control group (P<0.01). **P<0.01; *P<0.05.

Abbreviations: DMN, default mode network; FC, functional connectivity; L-ACC, left anterior cingulated cortex; L-HP, left hippocampus; L-SMG, left supra marginal gyrus; L-PCC, left posterior cingulate cortex; R-MOFC, right medial orbitofrontal cortex; R-PHP, right parahippocampus.

Correlation analysis

Fifty-five participants including 16 mild, 22 moderate, and 17 severe COPD patients underwent blood gas analysis. Pearson’s correlation analysis showed that the MoCA score was significantly correlated with PaO2 (r=0.368, P=0.006) and FEV1 (r=0.295, P=0.029), but not with PaCO2 (r=−0.056, P=0.685). The FC values of L-HP and L-PCC were positively correlated with PaO2 (r=0.276 and 0.433, P=0.042 and 0.001, respectively) but not with PaCO2 (all P>0.05).

The FEV1 of all COPD patients was positively correlated with the FC values of L-ACC (r=0.410, P=0.001), L-PCC (r=0.457, P=0.001), L-HP (r=0.248, P=0.024), and R-MOFC (r=0.394, P=0.001). The MoCA score was positively correlated with the FC values of L-HP (r=0.233, P=0.034) and L-PCC (r=0.215, P=0.050). The visuospatial executive function score was positively correlated with the FC value of L-PCC (r=0.387, P=0.001), while it was negatively related to L-SMG (r=−0.323, P=0.003). Meanwhile, the naming function score was negatively correlated with the FC value of the R-PHP (r=−0.225, P=0.041), and the orientation score was significantly correlated with the FC value of the R-MOFC (r=0.263, P=0.018).

Discussion

In the present study, we demonstrated that COPD patients developed cognitive impairment, which was consistent with numerous published studies.15,23 After classifying the COPD patients into three groups (mild, moderate, and severe) by FEV1, cognitive function was significantly different among the three groups and declined with worsening of the FEV1. Cognitive impairment mainly manifested in visuospatial executive function, attention, abstraction, delayed recall, and orientation. Using rs-fMRI, we found that activated brain regions in the DMN of COPD patients were obviously less than that of normal control subjects, especially the severe COPD group. Meanwhile, six activated brain regions in the DMN were found to attain significantly different FC values among the four groups. The MoCA scores correlated well with the FC values of L-HP and L-PCC.

Among the various brain mode networks in the resting state, the DMN is the most active24 and is mainly composed of the ACC, PCC, medial prefrontal cortex, dorsal thalamus, cuneus, precuneus, HP, and temporal lobe. The seed-based method to analyze FC values in the DMN is now widely used in the study of neuropsychiatric diseases.22,25 In this study, we chose L-PCC as the seed to analyze the FC values because of its critical role in the DMN connection.21 A recent study showed that changes in FC values in brain regions occurred earlier than the structural changes of brain gray matter.26 In this study, we found that there were six activated brain regions in the DMN that had significantly different FC values among the four groups. Meanwhile, compared with normal control subjects, the severe COPD patients had significant lower FC values in L-ACC, R-MOFC, and L-HP, while mild COPD patients had significantly lower FC values in L-HP and R-PHP; however, in moderate COPD patients, the FC values of some brain regions (R-PHP and L-HP) were higher than those in the mild and severe groups (shown in Figure 3). We speculate that the brain function of COPD patients possibly has the ability to compensate, and this ability may be the strongest in moderate COPD patients.

Both the HP and PHP belong to the medial temporal lobe system, and they receive cortex information from the frontal lobe, parietal lobe, and occipital lobe. Previous studies have shown that the HP was obviously related to recent memory and this study mainly refers to spatial learning and situational memory.27,28 The cingulate gyrus is an important component of the limbic system, which is closely related to the higher-level activities of humans, such as emotion, study, and memory. Nestor et al found that patients with mild cognitive impairment may have a severe reduction of metabolism of the HP and PCC.29 Zhang et al reported significant decreases in the amplitude of low-frequency fluctuation in the bilateral PCC of COPD patients.10 In our study, we found that the FC values of L-HP and L-PCC were positively correlated with both FEV1 and cognitive function scores, which may indicate that brain function of COPD patients varied with the disease severity and that L-HP and L-PCC may play an important role in the process of COPD patients’ cognitive impairment. A single score analysis revealed that the visuospatial executive function score was related to the FC values of L-PCC and L-SMG, as well as the relations between naming function and R-PHP. We speculate that the FC values of activated brain regions in the DMN may predict cognitive impairment of COPD patients; however, further studies are needed to gain more insight into the underlying mechanisms.

COPD is a chronic pulmonary disease, and most patients exhibit hypoxia and/or hypercapnia with the development of the disease. Oxygen deficiency can cause increased blood viscosity, which may directly lead to a decrease in the perfusion of the brain’s prefrontal lobe, reducing the blood supply to some brain regions. Moreover, oxygen deficiency can increase the production of oxygen free radicals and damage the neurons.30,31 Some researchers have considered that decreased cellular metabolism and brain perfusion caused by oxygen deficiency are the main causes of damage to neurons.32,33 In our study, we demonstrated that the cognitive function of COPD patients correlated well with PaO2, but not PaCO2. However, cognitive impairment can also be found in COPD patients without hypoxemia,20 and cognitive functions can partially recover after receiving long-term oxygen therapy.1,34 We speculate that the possible mechanism may lie in that the cognitive impairment of COPD patients could be caused by various independent factors, including abnormalities of brain structure, decreased lung function, hypoxemia, and hypercapnia.35–38 More research is needed to gain more insight into this phenomenon. Meanwhile, we found that the FC values of L-PCC and L-HP were positively correlated with PaO2; this is reasonable as oxygen deficiency can increase the amount of deoxygenated hemoglobin, thus affecting the BOLD signal.

There are several limitations in this study. Firstly, this study was a cross-sectional study to better interpret the role of hypoxia in the change of brain function and neural mechanisms of DMN; further longitudinal studies about cognitive function and disease progression of COPD are necessary. Secondly, the very severe COPD patients were not included in our study because they could not undergo rs-fMRI. Other limitations included the relatively small sample size and the fact that we did not analyze the effects of disease duration. All of these factors may affect the representativeness of our results. We will expand the sample size in further studies to improve the reliability of our conclusion.

Conclusions

COPD patients have cognitive impairment, which correlates well with disease severity. FC changes in activated brain regions in the DMN may predict cognitive impairment, and L-PCC and L-HP may be the important brain regions that are related to cognitive impairment of COPD patients.