Grey and white matter abnormalities in chronic obstructive pulmonary disease: a case–control study
- Haiyan Zhang1,
- Xiaochuan Wang2,
- Jianzhong Lin3,
- Yinchuan Sun1,
- Yongxia Huang1,
- Tianhe Yang3,
- Shili Zheng4,
- Ming Fan5,
- Jiaxing Zhang1
- 1Department of Physiology, Medical College of Xiamen University, Xiamen, China
- 2Department of Neurology, the Affiliated Hospital of Hangzhou Normal University, Hangzhou, China
- 3Magnetic Resonance Center, Zhongshan Hospital, Medical College of Xiamen University, Xiamen, China
- 4Department of Respiratory, Zhongshan Hospital, Medical College of Xiamen University, Xiamen, China
- 5Department of Brain Protection and Plasticity, Institute of Basic Medical Sciences, Beijing, China
- Correspondence to Dr Jiaxing Zhang;
- Received 12 January 2012
- Accepted 6 March 2012
- Published 24 April 2012
Objectives The irreversible airflow limitation characterised by chronic obstructive pulmonary disease (COPD) causes a decrease in the oxygen supply to the brain. The aim of the present study was to investigate brain structural damage in COPD.
Design Retrospective case–control study. Patients with COPD and healthy volunteers were recruited. The two groups were matched in age, gender and educational background.
Setting A hospital and a number of communities: they are all located in southern Fujian province, China.
Participants 25 stable patients and 25 controls were enrolled from December 2009 to May 2011.
Methods Using voxel-based morphometry and tract-based spatial statistics based on MRI to analyse grey matter (GM) density and white matter fractional anisotropy (FA), respectively, and a battery of neuropsychological tests were performed.
Results Patients with COPD (vs controls) showed decreased GM density in the limbic and paralimbic structures, including right gyrus rectus, left precentral gyrus, bilateral anterior and middle cingulate gyri, bilateral superior temporal gyri, bilateral anterior insula extending to Rolandic operculum, bilateral thalamus/pulvinars and left caudate nucleus. Patients with COPD (vs controls) had decreased FA values in the bilateral superior corona radiata, bilateral superior and inferior longitudinal fasciculus, bilateral optic radiation, bilateral lingual gyri, left parahippocampal gyrus and fornix. Lower FA values in these regions were associated with increased radial diffusivity and no changes of longitudinal diffusivity. Patients with COPD had poor performances in the Mini-Mental State Examination, figure memory and visual reproduction. GM density in some decreased regions in COPD had positive correlations with arterial blood Po2, negative correlations with disease duration and also positive correlations with visual tasks.
Conclusion The authors demonstrated that COPD exhibited loss of regional GM accompanied by impairment of white matter microstructural integrity, which was associated with disease severity and may underlie the pathophysiological and psychological changes of COPD.
Decreased oxygen supply to brain may cause neuronal damage in COPD. However, the damage remains largely uninvestigated.
We found that COPD extends to the brain, with the loss of regional cortical grey matter accompanied by impairment in the white matter microstructural integrity.
Our findings would be help for clinical therapy of COPD.
Strengths and limitations of this study
Multiple analyses were used based on MR images. The statistic power for FA analysis was weak.
Chronic obstructive pulmonary disease (COPD) continues to be a major cause of morbidity and mortality. It is increasingly recognised that COPD extends beyond the lung.1 The irreversible airflow limitation characterised by COPD usually develops arterial oxygen desaturation, which could subsequently result in a decrease in oxygen transport to the brain. Hypoxia during COPD has been previously proven to induce cerebral perfusion decline2 and metabolic changes.3–6 Moreover, systematic inflammation1 may also cause neuronal damages in the brain of patients with COPD. In patients with COPD, clinical symptoms such as neuropsychological deficits,7 depression and anxiety,1 and physical disability1 have been well documented. Taken together, these data suggest the presence of brain structural alteration. However, until now, it remains largely uninvestigated.
Voxel-based morphometry (VBM) and Tract-Based Spatial Statistics (TBSS)8 based on MRI were adopted to measure grey matter (GM) density and white matter (WM) fibrous microstructure properties in tracts, respectively. VBM is an automatic quantitative volumetric technique over the whole brain using voxel by voxel analysis without prior speciﬁcation of regions-of-interest for analysis, and it does not rely on arbitrarily predeﬁned structures. Recently, the preprocessing steps of VBM have been improved with the Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra (DARTEL) registration method,9 which can achieve more accurate intersubject registration of brain images. TBSS is a recently introduced method, which uses diffusion tensor MR imaging (DTI) to measure differences in fractional anisotropy (FA) between groups. TBSS increases the sensitivity and the interpretability of the results compared with voxel-based approaches based purely on non-linear registration.8 Moreover, diffusion tensor eigenvalues were also included in the analysis since they can help interpret FA changes in WM tracts by providing information regarding likely alterations in the proportion of longitudinally versus obliquely aligned myelinated fibres. The VBM and TBSS methods have been extensively applied in clinical researches, including the evaluation of morphological characteristics of high-altitude residents in our previous study.10
Dyspnoea is the most common complaint and most disabling symptom in patients with COPD. Functional MRI studies on breathlessness, air hunger and inspiratory loaded breathing have revealed that a large number of brain regions, including the frontal cortex, parietal cortex, temporal cortex, limbic cortex, cerebellar cortex and brainstem were activated by dyspnoea.11 These dyspnoea-activated brain regions have been shown to be impaired in patients with congenital central hypoventilation syndrome,12 in patients with obstructive sleep apnoea13 and in high-altitude residents.10 We therefore hypothesised that patients with COPD would have similar cerebral impairment.
Twenty-five patients were enrolled from December 2009 to May 2011. All patients had undergone a period of 30–45 days of in-hospital rehabilitation following an acute exacerbation of COPD. At the time of data collection, patients were in stable condition. Among these patients, 12 discharged patients were recruited during their rest at home and 13 patients were recruited when they were awaiting discharge from hospital. Patients were diagnosed in Zhongshan Hospital (Xiamen, China) according to the diagnostic criteria of Global Initiative for Chronic Obstructive Lung Disease.14 Twenty-five healthy volunteers, with comparable age, gender and educational background, comprised the control group. All the subjects were free from a known history of cerebrovascular accident, heart failure, neurological disorders, obstructive sleep apnoea, coronary artery disease, diabetes or other diseases known to affect cognition. Patients were provided with therapy including inhaled ipratropium bromide, bricanyl, ventoline and budesonide. Demographic characteristics of the patients and healthy volunteers were listed in table 1. Procedures were fully explained, and all subjects were provided written informed consent before participating in the study. The experimental protocol was approved by the Research Ethics Review Board of Xiamen University.
Physiological and neuropsychological tests
Physiological and neuropsychological tests and activities of daily living (ADL) (score range 14–56)15 were conducted 1 day before the MRI scan. Physiological tests include pulse rate and arterial blood pressure measures, arterial blood gas analysis and pulmonary function measure. Blood samples were taken in the morning between 07:00 and 07:30 h. The neuropsychological tests include (1) the Chinese version of the Mini-Mental State Examination (MMSE) that measured the general cognitive function and (2) the visual reproduction test, figure memory test and digital span forward and backward tasks, which, taken from the Chinese revised version of the Wechsler Memory Scale, were used to measure visual construction ability, visuospatial memory and short-term working memory, respectively.16 All data were analysed using SPSS V.19.0. Independent t test measured between-group differences. Statistical significance was set at p<0.05.
MRI data acquisition
Images were acquired on a Siemens Trio Tim 3.0T (Erlangen, Germany) at MRI Research Center (Zhongshan Hospital). A three-dimensional (3D) structural MRI was acquired from each subject using a T1-weighted MPRAGE sequence (TR/TE=1900 ms/2.48 ms, FOV=25×25 cm2, NEX=1, matrix=512×256, slice thickness=1.0 mm). Conventional two-dimensional T1 and T2 images were also acquired and examined for any incidental findings. A DTI pulse sequence with single shot diffusion-weighted echo planar imaging (TR/TE=3600/95 ms, FOV=24×24 cm2, NEX=2, matrix=128×128, slice thickness=4.5 mm) was applied sequentially in 30 non-collinear directions (b-value=1000 s/mm2) with one scan without diffusion weighting (b=0s/mm2). The following data analyses were conducted by two researchers who were blinded to the status of subjects.
VBM analysis of 3D T1 images
The 3D T1 images were used for GM analysis using VBM8 toolbox implemented in SPM8 (Wellcome Department of Imaging Neuroscience, University College London, London, UK). The following processing steps were carried out: (1) the images were inspected and set at the anterior commissure. Each reorientated image was segmented into GM, WM and cerebrospinal fluid in native space, and procrustes aligned GM images were generated by a rigid transformation. (2) The DARTEL registration method was used to create a study-specific template by using the aligned images from all the patients and controls to improve intersubject registration of structural images.9 The procedure implicated in six iterations, which began with the averaging of aligned data to generate an original template. Then, the first iteration of the registration was completed on each subject and a new template was created. After this, the second iteration began. When six iterations were finished, the template was generated, which was the average of the DARTEL registered data. During iterations, all images were warped to the template yielding a series of flow fields that parameterised deformations. (3) The normalised images were transformed into Montreal Neurological Institute space. These GM images were then smoothed using a Gaussian kernel of 8 mm full-width at half-maximum. Independent t tests were performed to examine between-group differences. The statistical parametric map was generated with the voxel level threshold at t>3.7734 p<0.01 (false discovery rate FDR correction with gender, age, education and total intracranial volume as covariates).
To analysed the correlation of GM image values with cognitive or physiological measurement, the following steps were first taken: (1) regions-of-interests were created for clusters showing differences between groups and (2) using these regions-of-interests masks, the GM values were extracted from each individual's normalised and smoothed GM maps. Then, the correlations were analysed using SPSS. Statistical significance was set at p<0.05, with gender, age and education as covariates.
TBSS analysis of DTI
DCM2MII was used to convert diffusion tensor images from the proprietary scanner format to the NIFTI format. Then, the images were processed using the FSL 4.1.5 software package (http://www.fmrib.ox.ac.uk/fsl/). The images were realigned to the b-value (b0) image by affine transformations using FMRIB's diffusion toolbox17 to minimise distortions and reduce head motion artefacts. In order to remove non-brain tissue components and background noise, a brain mask was created from the first b0 image and then applied in the DTI to extract brain voxels using Brain Extraction Tool. After these processes, using DTIFit within the FMRIB's diffusion toolbox, the images were calculated to get the FA, λ1 (longitudinal diffusivity) (the magnitude of diffusion along the principal diffusion direction) and λ23 (radial diffusivity) (the magnitude of diffusion in the two orthogonal directions perpendicular to the principal diffusion direction) maps. The whole-brain voxel-wise statistic analysis of the FA images was performed using TBSS in FSL.8 TBSS processing includes the following steps: (1) align the FA images of all subjects to a template that was arbitrarily selected from those FA images by non-linear registrations; (2) transform all the aligned FA images into 1×1×1 mm3 MNI152 space by affine registrations to remove the effect of cross-subject spatial variability that remains after the non-linear registration; (3) create the mean FA image and filter to retain only the centre of the WM tracts, with the threshold FA≥0.20, and successfully exclude voxels, which consisted of GM or cerebrospinal fluid in the majority of subjects, so as to create the mean FA skeleton8; (4) project individual subjects' FAs onto mean FA skeleton; (5) following these steps, data were fed into voxel-wise cross-subject statistical analyses. In all cases, the null distribution was built up over 5000 permutations, with significance analysed using independent t tests at p<0.05 levels, uncorrected for multiple comparisons. We determined the anatomic localisation of each cluster by means of the FSL atlas tool, which incorporates several anatomic templates, including the Harvard-Oxford Cortical Structural Atlas, Harvard-Oxford Subcortical Structural Atlas, Talairach Daemon Labels and MNI Structural Atlas.
Within the cluster of changed FA, mean λ1 and λ23 values were extracted from each individual's λ1 and λ23 maps. Values were analysed using SPSS. Analysis of variance statistic was used to identify the group differences for these distinct brain locations. Statistical significance was set at p<0.05.
Physiological and behavioural findings
Compared with the controls, independent t test showed that patients with COPD had significant decreases in arterial blood Sao2 and Po2, and increases in arterial blood Pco2 and heart rate. Patients with COPD had significantly lower values in one second over forced vital capacity (FVC), forced expiratory volume (FEV) and FEV1/FVC values and higher respiratory rate. The disease staging categories of patients with COPD based on FEV1% predicted were as follows: FEV1 =82% predicted, n=1; 54.9%≤FEV1<78% (64.9±7.6) predicted, n=8; 31.9%≤FEV1<48.4% (42.0±5.4) predicted, n=7; FEV1<29.9% (26.2±2.9) predicted, n=9. Patients with COPD had significantly lower scores in ADL, MMSE test, visual reproduction and figure memory (table 2).
No subject from either group showed visible abnormalities on T1-weighted structural images. VBM analysis showed that patients with COPD had decreased GM densities compared with healthy controls in the right gyrus rectus, left precentral gyrus, bilateral anterior and middle cingulate gyri, bilateral superior temporal gyri, bilateral anterior insula extending to Rolandic operculum (base of the pre- and post-central gyri), bilateral thalamus/pulvinars and left caudate nucleus (cluster size >100 voxels) (figures 1–4, table 3).
FA, longitudinal diffusivity and radial diffusivity in relation to COPD
Whole-brain voxel-wise statistic analysis showed that patients with COPD had significantly lower FA in a broad range of brain regions compared with controls (figure 5, table 4). The significantly affected regions (clusters size >40 voxels) included the superior corona radiata (corresponding to bilateral precuneus and bilateral superior parietal lobules), superior longitudinal fasciculus (bilateral supramarginal gyri), inferior longitudinal fasciculus (left superior temporal gyrus, right middle temporal gyrus and fusiform gyrus), bilateral optic radiation, bilateral lingual gyri, left parahippocampal gyrus and fornix.
Lower FA values in these regions were associated with increased radial diffusivity and no changes of longitudinal diffusivity in patients with COPD versus controls (table 4).
Correlations between MRI measurement and disease severity
The correlations were listed in table 5. In patients with COPD, partial correlation (controlling for disease duration, FEV1/FVC, age, education and gender) revealed that the GM density in the bilateral anterior cingulate cortex, left superior temporal cortex, bilateral insula/superior temporal/Rolandic operculum, bilateral thalamus/pulvinar and left caudate nucleus had positive correlations with arterial blood Po2. Partial correlation (controlling for Po2, FEV1/FVC, age, education and gender) revealed that the GM density in the bilateral anterior cingulate cortex, right insula/superior temporal/Rolandic operculum and right thalamus/pulvinar had negative correlations with disease duration. Partial correlation (controlling for age, education and gender) analysis showed that the GM density in the left superior temporal lobes and left insula/superior temporal/Rolandic operculum in patients with COPD was significantly correlated with figure memory score and the GM density in left precentral gyrus and left thalamus/pulvinar in patients with COPD correlated significantly with visual reproduction.
Our present study revealed that patients with COPD had decreased regional GM density confined to the limbic and paralimbic structures. GM density in impaired regions in patients with COPD had significant positive correlation with arterial blood Po2 and negative correlation with disease duration. The decreased WM FA value with increased radial diffusivity value was detected mainly in the visual cortex of the occipital lobe, the posterior parietal lobe as well as the temporal lobe. Decreased FA was associated with compromised myelin structure, changes in axonal morphologic structure and altered interaxonal spacing of fibre bundles.18 Radial diffusivity measures motion of water molecules perpendicular to fibres, and an increase of radial diffusivity is interpreted as abnormalities in myelinated membranes.19 Consequently, decreased FA and increased radial diffusivity in COPD indicated the impairment of WM microstructural integrity.
Previously, Borson et al6 only measured the volume of hippocampus in patients with COPD using region-of-interest analysis and did not find any significant change. The impaired brain regions in COPD have also been found in other chronic hypoxic diseases. For example, decrease in GM volume/concentration in the gyrus rectus, precentral gyrus, anterior cingulate cortex, multiple sites within the temporal lobes, insular cortex, thalamus and caudate nucleus were detected in patients with obstructive sleep apnoea.13 Impairments of WM microstructure in the temporal lobe, parietal lobe, fornix and corona radiata were found in patients with congenital central hypoventilation syndrome.12 In our previous study, the decrease in GM volume in the anterior insula, anterior cingulate cortex and precentral cortex were found in high-altitude residents.10
GM density in impaired regions in patients with COPD had a strongly positive correlation with the arterial blood Po2, which suggested that the impairment in GM may result from low blood oxygen. Moreover, the GM density in some impaired regions showed negative correlations with disease duration. It is already known that hypoxia can induce metabolic decreases3–5 and cerebral perfusion decline2 in COPD. In addition, patients with COPD often suffer from systemic inflammation, which can exacerbate neuronal injury.1 A greater proportion of regions showing GM loss located in limbic/paralimbic cortex in patients with COPD may be due to the fact that phylogenetically older regions of the brain showed sharper vascular responses to hypoxia than evolutionary younger regions.20
A larger cortical network including the anterior insula and anterior cingulate cortex underlie the perception of dyspnoea,11 and these regions play an important role in regulating the cardiovascular system.21 Posterior thalamus was implicated in suppressing the ventilatory response to hypoxia.22 Hippocampus has been proved to control arterial pressure and heart rate.23 Thus, the morphological impairments in these regions may play a role in respiratory and cardiovascular disturbances, such as higher heart rate and higher respiratory rate, in patients with COPD tested in our study.
In the present study, patients with COPD had poorer performance in MMSE, visuospatial memory and visual construction task. These results were consistent with those found in previous studies in patients with COPD.7 In line with the present findings of COPD, we previously found that long-time living in mild high-altitude hypoxic environment only impaired cognitive performances confined to visual reproduction and short-time complex figure memory.24 In our present study, the decreases in GM density in the frontal precentral cortex, insula/superior temporal cortex/Rolandic operculum and thalamus/pulvinar may be responsible for the performance deficit in visual-related tasks since the GM density in these areas showed a significant positive correlation with figure memory or visual reproduction score. The following data support our findings: (1) recent research has identified the inferior frontal cortex served as a source of top–down modulation underlying attention to visual features25; (2) studies on patients using fMRI and positron emission tomography demonstrated Rolandic operculum as one of the visual structures26; (3) the pulvinar region of the thalamus is known to project to posterior parietal lobe and inferior temporal lobe. The pulvinar has been implicated in various visual functions in lesion studies.27 Declines in memory and executive function make contributions to declines in ADL.28 Visual construction tasks reflect executive function. Therefore, the decreases in GM density in the above regions that relate to visual construction may also be responsible for ADL deficits.
Our present study found the impairments of WM limited to the pathways of visual processing, including optic radiation, posterior parietal lobe (superior parietal lobule, supramarginal gyrus and precuneus) and the inferior temporal fusiform and lingual gyri. Visual information enters the primary visual cortex via optic radiation to the visual cortex. Cortical areas along the posterior parietal ‘dorsal stream’ are primarily concerned with spatial localisation and directing attention, while cortical areas along the inferior temporal ‘ventral stream’ are mainly concerned with the recognition and identification of visual stimuli.29 COPD also showed impaired WM in the middle temporal gyrus. Middle temporal cortex is important for the long-term build-up of perceptual memory for ambiguous motion stimuli.30 Based on the above data, our findings in WM may also clarify the mechanisms underlying the deficit in visual-related tasks. In addition, impaired WM in input and output fibres of hippocampus (fornix) may be related to the deficit in MMSE. Previous study on patients with Alzheimer's disease found that the volumes of hippocampus were significantly reduced and the volumes of the left hippocampus correlated significantly with the MMSE score.31 The limitation of our study is the weak statistical power of FA value analysis because the results obtained in the TBSS analysis could not survive multiple comparison correction.
In summary, we first demonstrated that COPD extended beyond the lung to the brain, with the decrease of regional GM density accompanied by impairment in the WM microstructural integrity. Our ﬁndings suggest signiﬁcant participation of these structures in responding to hypoxic challenges, which include cardiovascular and air-hunger components. The brain structural changes may also underlie the psychological and mood changes of COPD.
HZ and XW contributed equally to this work.
To cite: Zhang H, Wang X, Lin J, et al. Grey and white matter abnormalities in chronic obstructive pulmonary disease: a case–control study. BMJ Open 2012;2:e000844. doi:10.1136/bmjopen-2012-000844
Contributors HZ contributed to conception and design, the accuracy of the data analysis, drafting the article and final approval of the version to be published. XW contributed to the study design, revising the article critically for important intellectual content and final approval of the version to be published. JL, YS, YH, TY, SZ and MF contributed to acquisition of data, revising the article critically for important intellectual content and final approval of the version to be published. JZ contributions to conception and design, acquisition of data and interpretation of data; drafting the article and revising it critically for important intellectual content; and final approval of the version to be published.
Funding This work was supported by National Science Foundation of China (Project No. 31071041; 81171324), Bureau of Science and Technology of Hangzhou (20090833B10) and Health Bureau of Zhejiang Province (2009A168).
Competing interests None.
Patient consent Obtained.
Ethics approval The experimental protocol was approved by the Research Ethics Review Board of Xiamen University.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement We are pleased to share our data with other researchers. Jiaxing Zhang et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.