Altered microstructure of brain white matter in females with anorexia nervosa: a diffusion tensor imaging study.

OBJECTIVE: Structural studies have reported anorexia nervosa (AN) patients with abnormal gray matter in several brain regions and dysfunction in some connected neural circuits. However, the role of white matter (WM) in AN patients has rarely been investigated. The present study aimed to assess alterations in WM microstructure of the entire brain in females with AN using a voxel-based method on diffusion tensor imaging (DTI) data.
MATERIALS AND METHODS: The study enrolled 8 female patients with AN and 14 age-matched females as controls (CW). The DTI data was collected from each subject to calculate the fractional anisotropy (FA) maps of the whole brain by the DTI-Studio software. Subsequently, a 2-sample t-test (P<0.05, corrected) was performed to detect the difference in FA maps of AN and CW group, and a Pearson's correlation analyzed the relationship between mean FA value of brain regions and body mass index (BMI).
RESULTS: Compared with CW, AN patients revealed a significant decrease in FA maps in the left superior frontal gyrus, medial frontal gyrus, anterior cingulate cortex, middle frontal gyrus, inferior frontal gyrus, thalamus, and bilateral insula. Moreover, significantly positive correlations were established between the mean FA value of the left inferior frontal gyrus, insula as well as thalamus and BMI in AN patients.
CONCLUSIONS: Our findings supported the presence of WM abnormality in patients with AN. The significant differences of FA maps, in patients with AN, were associated with their aberrant BMI. The results further improved our understanding of the pathophysiological mechanisms underlying AN.

Introduction

Anorexia nervosa (AN) is a severe psychiatric eating disorder of unknown etiology related to self-driven food refusal and emaciation, altered body perception, and preoccupations with weight and shape.1 It is primarily characterized by disrupted eating patterns or weight-control behavior resulting in heightened anxiety and significantly increased damage to health, causing impaired physical and psychosocial functions.2–4 Because of its low therapeutic compliance and high expulsion rate, it has been a crucial factor of adolescent deaths in recent years.5

Recent studies demonstrated that the structural abnormality of the brain was a prerequisite for the development of AN.6 Several studies have shown alterations of gray matter volume (GMV), using voxel-wise analysis, in AN patients compared with the control group. For instance, Mühlau et al found a decreased global GMV and regional GM in the bilateral anterior cingulate cortex (ACC) in AN patients by ~1% and 5%, respectively, which remarkably correlated with the lowest body mass index (BMI).2 Boghi et al found a significant reduction in total white matter (WM) volume and local GM atrophy in cerebellum, hypothalamus, caudate nucleus, and frontal, parietal and temporal areas. In addition, an association between BMI and GM was also found in the hypothalamus.7 Fonville et al detected smaller global GMV and regional decrease in the cerebellum, temporal, frontal, and occipital lobes in adult participants with AN, and a relation between GMV loss and duration of illness in the cerebellum and mesencephalon.8 Moreover, some of these volume depletions may return to normal in recovered patients,9 while other studies have shown their persistence.2,10,11

Structural neuroimaging studies of AN primarily concentrate on alterations of gray matter. To date, study focused on the abnormality of WM is rare. Brain WM axons connect between cortical and subcortical brain structures and thus have a critical impact on cognitive and emotional processing.12 Diffusion tensor imaging (DTI), an advanced magnetic resonance (MR) technique, can reveal the microstructural integrity of WM pathways.13 One of the most important DTI metrics, fractional anisotropy (FA) is capable of measuring water diffusion along the WM tracts. It has been previously applied in other psychiatric disorders to investigate the change of WM microstructure.14 Lower FA inversely affects the WM integrity because of either fragile axonal coherence, density or myelination.15,16 A recent DTI study suggested that WM integrity of AN patients was destroyed in limbic and association pathways, including fimbria-fornix, fronto-occipital, and posterior WM cingulum. The fimbria-fornix WM integrity can also predict the harm avoidance.12 Another tract-based study found that AN patients showed significant FA decrease in the parietal part of the left superior longitudinal fasciculus, with increased mean diffusion (MD) and radial diffusion (RD) but unaltered axial diffusion (AD). Patients with AN also revealed remarkably increased MD in the fornix, accompanied by low FA, and high RD and AD.17 However, these studies mainly focused on the altered WM integrity in several fiber tracts related to AN, but studies investigating WM microstructure of the whole brain using voxel-based method are scarce.

Although a recent voxel-based DTI study on AN has been reported, the outcomes presented only reduced FA value in the bilateral posterior thalamic radiation, while the correlation analysis between decreased FA value and BMI of AN patients was not explored.6 Previous studies have demonstrated that the common cognitive dysfunctions, such as executive functions, visuospatial processes, self-image perception, emotional regulation, and reward processing, were found in AN patients.18 Therefore, to explore the mechanism underlying the etiology of AN, the voxel-based analysis of the entire brain on DTI data was performed between AN and controls in the present study. We hypothesized that the altered FA would be found in widespread regions, such as prefrontal cortex (PFC), parietal lobe, and subcortical regions, that are closely associated with the impaired functions in AN patients. Concurrently, significant correlations would be exhibited between altered FA and BMI in AN patients.

Materials and methods

Participants

Eight right-handed female patients with AN (average age 17.6±2.2 years, average height 162±3.0 cm), fulfilling Diagnostic and Statistical Manual of Mental Disorders fourth edition text revision criteria for the restricted subtype of AN,19 from the Department of Psychiatry, First Affiliated Hospital of Zhejiang University School of Medicine, China, were consecutively recruited. In addition, 14 controls (females, average age 19.1±3.1 years, average height 161±3.9 cm) matched in age and height (P=0.239, P=0.600, respectively, 2-sample 2-tailed Student’s t-test), gender, and dexterity, were recruited via proximal advertisements. Comorbid psychiatric disorders, neurologic conditions, and abuse of any substance except nicotine comprised the exclusion criteria. Before scanning, all patients were required to have had at least 1 week of supervised meals and hydration during their regular hospital admission. This minimized the alleged biases caused by dehydration in brain measurements.20 All controls were screened using the General Health Questionnaire and a clinical semi-structured interview to exclude any psychiatric or other medical condition.21 We also ensured that the controls had a BMI within the healthy range (average BMI, 20.1±1.7 kg*m−2) and that they did not present unhealthy eating behaviors (eg, constant dieting) or sub-threshold symptoms of any eating disorder. The demographic characteristics of each patient are displayed in Table 1. This study was approved by Institutional Ethics Committee of the First Affiliated Hospital of Zhejiang University School of Medicine. Informed consent was obtained from all participants or legal guardians for those under the legal age of consent.

Table 1

Demographic characteristics and clinical assessment of AN patients

ID	Age (years)	Gender	Handedness	Duration (month)	Height (cm)	Weight (kg)	BMI (kg/m2)
01	16	F	R	6	162	32.5/50a	12.7
02	18	F	R	12	164	40/54	14.9
03	17	F	R	12	160	42.5/48	16.6
04	22	F	R	6	168	37/51	13.1
05	15	F	R	6	158.5	33.5/42	13.3
06	16	F	R	6	163	38/44	14.4
07	19	F	R	24	160	40/44	15.6
08	18	F	R	12	161	36/51	13.9

Note:

Represents the standard weight for each patient.

Abbreviations: AN, anorexia nervosa; BMI, body mass index; F, female; R, right-handed.

Image acquisition and processing

Participants were scanned using a 3T MR imaging system (GE sign HDXM) with an 8-channel phased array head coil. Cushions restricted head movements, and the participants wore earplugs. DTI data was obtained using a single-shot echo-planar imaging sequence in 44 axial planes with 25 nonlinear diffusion sensitization gradients (b=1,000 s/mm2), and a reference image without diffusion weighting (b0-image). Imaging parameters were: repetition time/echo time (TR/TE), 12,000/89 ms; slice thickness, 3 mm (no slice gap); number of excitations, 2; matrix, 128×128; field of view (FOV), 240×240 mm2; voxel size, 1×1×3 mm3. DTI was performed using axial sections parallel to the anterior–posterior commissural line to cover the entire brain. As an anatomical reference, high-resolution structural images were acquired using a magnetized rapid gradient-echo (MPRAGE) 3-dimensional T1-weighted sequence: TR/TE, 7.1/2.9 ms; 146 slices with thickness 1 mm; flip angle, 8°; matrix, 256×256; FOV, 240×240 mm2; voxel size, 1×1×1 mm3. DTI data processing and analysis were performed using DTI-Studio program SPM8 (Johns Hopkins Medical Institute, Baltimore, ML, USA), and software tools developed in-house. Subsequently, the echo planar distortions induced by eddy-current were corrected using an affine transformation algorithm, following which the parametric maps of FA were created to perform the voxel-wise analysis and region of interest (ROI) analysis.

Voxel-wise analysis of FAmaps between AN and CW

The voxel-wise of FA maps was conducted by using statistical parametric mapping (SPM8). First, b0 images of each CW were normalized to Montreal Neurological Institute space by an echo-planar imaging template, and then averaged. Next, the averaged images were spatially leveled using an anisotropic Gaussian filter with a full width at half maximum of 6 mm. Thus, we obtained a new template of b0 image in the Montreal Neurological Institute space. Subsequently, individual b0 images of all patients and CW were registered to the new template (the normalization used the nonlinear registration method in SPM, including affine transformation and nonlinear transformation) and corresponding transformation matrices were also obtained. FA maps of each participant were then normalized to Montreal Neurological Institute space by using their corresponding transformation matrix and leveled by an anisotropic Gaussian filter with a full width at half maximum of 6 mm. Finally, a two-sample Student’s t-test was performed on the leveled FA of patients and CW. The threshold was set at a combined cut-off value of P<0.05 (false discovery rate corrected) and a minimal cluster size of 50 voxels.

ROI analysis of FA maps

According to the significant differences in FA maps detected by voxel-wise analysis, we identified and analyzed the ROIs to quantify further the degree of secondary degeneration in patients compared with CW. First, we made binary masks by using the clusters with significant group differences of FA maps between groups. Then, the respective binary masks were applied to the normalized FA maps of all the participants. Finally, a value for each ROI was obtained by averaging all voxel values within each mask. Furthermore, correlation analysis was carried out between BMI and each ROI value across all patients and a P-value <0.05 was considered to be statistically significant.

Clinical assessment

The BMI is one of the most important markers of AN, and its relevance has also been recognized from a general health perspective. A low BMI at the beginning of treatment has been identified as pivotal risk factors for poor prognosis in AN patients. BMI in the present study is used to evaluate the illness severity of AN, which is calculated as the body mass divided by the square of the height in each subject.

Results

Voxel-wise analysis of FA maps

Compared with CW, the significantly decreased FA was found in the left superior frontal gyrus (SFG), medial frontal gyrus (medFG), ACC, middle frontal gyrus (MFG), inferior frontal gyrus (IFG), thalamus, and bilateral insula in AN groups (Figure 1; Table 2).

Figure 1

Brain regions showing decreased FA in the AN group relative to healthy controls.

Note: The color bars indicate the Student’s t-test value and the statistical at P<0.05 (false discovery rate-corrected).

Abbreviations: AN, anorexia nervosa; FA, fractional anisotropy; L, left; R, right.

Table 2

Regions showing decreased FA in AN group

Regions	BA	MNI coordinates			Cluster	T
		X	Y	Z
Left superior frontal gyrus	9	0	3	51	193	−7.50
Left medial frontal gyrus	8	0	42	24	184	−7.50
Left anterior cingulate cortex	32				118	−4.83
Left middle frontal gyrus	9/46	−36	24	27	57	−6.67
Left inferior frontal gyrus	44/45	−51	12	27	92	−5.66
Left thalamus	–	−3	−3	6	112	−5.20
Left insula	13	−36	12	−18	105	−6.02
Right insula	13	42	9	−9	52	−6.09

Abbreviations: AN, anorexia nervosa; BA, Broca’s area; FA, fractional anisotropy; MNI, Montreal Neurological Institute.

AN group showed reduced FA in the left SFG (P<0.01), medFG (P<0.01), ACC (P<0.01), MFG (P<0.01), IFG (P<0.01), thalamus (P<0.01), and bilateral insula (P<0.01) (Figure 2). More importantly, we found significantly positive correlations between the mean FA value of the left IFG (r=0.916; P=0.001), insula (r=0.898; P=0.002) as well as thalamus (r=0.815; P=0.014) and BMI of AN patients (Figure 3).

Figure 2

The mean FA of region of interests with differences between groups were displayed.

Note: **P<0.01.

Abbreviations: ACC, anterior cingulate cortex; AN, anorexia nervosa; FA, fractional anisotropy; HC, healthy controls; IFG, inferior frontal gyrus; L, left; MFG, middle frontal gyrus; R, right.

Figure 3

The clusters (yellow circles) were associated with significant differences between the AN and CW groups, and differences within which the FA index significantly correlated with BMI.

Abbreviations: AN, anorexia nervosa; BMI, body mass index; CW, control women; FA, fractional anisotropy; IFG, inferior frontal gyrus; L, left; R, right.

Reproducibility

The sample size was relatively small in the present study. The reproducibility and robustness of these findings on FA were validated by a leave-one-out validation. Specifically, 1 AN patient was left out of the sample, and the group comparisons were performed based on the permutated samples (ie, 8 AN vs 14 CW). This led to a total of 8 2-sample t-test images; for each voxel, the number of tests where this voxel exhibited significant differences between groups was calculated as the reproducibility of FA differences between AN patients and CW. Figure 4 indicates the highly reproducible altered FA patterns across these tests as well as the original test reported.

Figure 4

The reproducible patterns of FA variations in AN group.

Abbreviations: AN, anorexia nervosa; FA, fractional anisotropy; L, left; R, right.

Discussion

The findings of the present study supported our hypothesis that AN patients exhibit widespread changes of FA values in the regions related to AN compared with CW. Specifically, we found reduced FA values in the left SFG, medFG, ACC, MFG, IFG, thalamus, and bilateral insula in the patients with AN. Furthermore, significantly positive correlations were found between the mean FA value of the left IFG, insula as well as thalamus and BMI in AN patients, demonstrating that the altered FA value may reflect the degeneration of WM integrity in the particular regions of AN patients, which can be associated with the underlying AN mechanism.

PFC, involved in regulating emotion and motivation, is a hub of the cognitive control network.22 On the other hand, the left anterior and posterior insula activation is associated with gastric distension, and thus, is mediated by BMI;23 right anterior insula is connected with self-recognition, the “abstract representation of oneself,”24 and interoceptive awareness,25 and hence, a fixed perception of being fat while severely underweight presenting in AN patients.26 The altered signal was discovered in several structural and functional studies in patients with AN. In GM studies, reduced GMV was consistently observed in the prefrontal lobe and insula of AN patients compared with the control group.27–29 The task-based functional magnetic resonance imaging study showed increased activations in the PFC, insula, and amygdala of AN patients, using a distorted fat body image of the patients as a visual task. These regions were composed of a network related to the affective component.30–32 In a study with AN and bulimia nervosa, the patterns of connectivity between the insula and frontal brain regions distinguished the groups during a visual food-cue task.33 A similar study of AN indicated reduced resting-state connectivity in cognitive control contributing area in the frontal cortex.34 On the study of WM, FA values in adolescents with AN were lower in the fornix, posterior frontal, and parietal areas, but higher in anterior frontal, orbitofrontal, and temporal lobes.35 Another study manifested positive correlations between BMI and FA, suggesting that weight loss could be related to the changes in FA.36 In agreement with previous studies, decreased FA values in AN patients were revealed in the PFC and insula in the present study, which positively correlated with BMI. This may indicate that the altered WM integrity in PFC and insula is associated with dysfunctional emotions and cognition control in AN patients.

The primary function of the thalamus is to filter and relay sensory information from the spinothalamic tract to other brain regions.37 Using graph theory method, a recent network analysis in patients with acute AN, found locally decreased strength and degree, while increased path length in the middle/posterior insula and the thalamus. Reduced network efficiency in the thalamus and the middle/posterior insula may reflect the altered filtering and calibration of nociceptive, thermal, and other interoceptive signals that may contribute to abnormal representations of body size as well as hunger.37 The posterior thalamic radiation, involved in the processing of the body image revealed decreased FA values in AN patients.6 Voxel-wise analysis also found reduced GMV in the extrastriate body area (EBA) in females suffering from acute AN.38 The parietal cortex region containing the EBA showed alteration of the diffusivity via the apparent diffusion coefficient-analysis.6 The left mediodorsal thalamic nucleus revealed reduced FA value, which is connected to the ACC, the PFC, and the ventral striatum.6 All of these regions have been implicated in the pathology of AN, contributing to the impairments in cognitive domains, especially set-shifting ability, executive control and habit learning, and in reward processing. Decreased FA of the left thalamus in AN patients in the present study may suggest the degeneration of WM integrity in this region. Moreover, the significant positive correlation between mean FA in this region and BMI may indicate the structural basis of functional impairments of emotion and cognition induced by AN in the patients.

ACC involvement in the manifestation of psychiatric symptoms is well established,39,40 which appears to be a recent brain adaptation that is linked to cognitive and emotional processes unique to humans.41 Findings from functional neuroimaging studies as well as positron emission tomography highlight the role of ACC in AN patients.42–43 Recently, these outcomes were supported by several structural studies using voxel-wise analysis. For example, patients recovered from AN showed a reduced GMV in the anterior cingulum.2 A recent study using voxel-based morphometry in anorectic patients illustrated volume depletion in right anterior insular cortex, bilateral parahippocampal gyrus, left fusiform gyrus, left cerebellum, and right posterior cingulate while an increased volume of the right dorsolateral PFC was observed.28 Additionally, ACC plays a major role in the reward mechanism,44,45 which could be associated with set-shifting deficit among AN patients.46 The reduced FA value of ACC found in AN patients in the present study indicated that the damage of WM in the ACC may provide the structural mechanism underlying the dysfunction of emotions and reward processes.

Limitations and future directions

The present study uncovered abnormal FA in specific AN-related regions in female patients with AN using the voxel-wise analysis method on DTI data. Nevertheless, some limitations remain. First, the sample size was relatively small, and thus, further studies with larger sample size are indispensable. Second, body image distortion and other disorders with specific behavioral or cognitive phenomena and attitudes were not assessed; thus, we cannot evaluate whether there is a direct association between the observed microstructural abnormalities related to AN and the perceived body image distortion in the patients.

Conclusion

In this study, we have demonstrated that WM microstructural alterations are present in the patients with AN. By a voxel-based approach on DTI data, we identified regions with prominent FA changes and significant correlations with BMI. This is closely associated with current hypotheses about the pathology of AN.