PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(7): e11449.
Published online 2010 July 6. doi:  10.1371/journal.pone.0011449
PMCID: PMC2897842

Structural Modifications of the Brain in Acclimatization to High-Altitude

Pedro Antonio Valdes-Sosa, Editor

Abstract

Adaptive changes in respiratory and cardiovascular responses at high altitude (HA) have been well clarified. However, the central mechanisms underlying HA acclimatization remain unclear. Using voxel-based morphometry (VBM) and diffusion tensor imaging (DTI) with fractional anisotropy (FA) calculation, we investigated 28 Han immigrant residents (17–22 yr) born and raised at HA of 2616–4200 m in Qinghai-Tibetan Plateau for at least 17 years and who currently attended college at sea-level (SL). Their family migrated from SL to HA 2–3 generations ago and has resided at HA ever since. Control subjects were matched SL residents. HA residents (vs. SL) showed decreased grey matter volume in the bilateral anterior insula, right anterior cingulate cortex, bilateral prefrontal cortex, left precentral cortex, and right lingual cortex. HA residents (vs. SL) had significantly higher FA mainly in the bilateral anterior limb of internal capsule, bilateral superior and inferior longitudinal fasciculus, corpus callosum, bilateral superior corona radiata, bilateral anterior external capsule, right posterior cingulum, and right corticospinal tract. Higher FA values in those regions were associated with decreased or unchanged radial diffusivity coinciding with no change of longitudinal diffusivity in HA vs. SL group. Conversely, HA residents had lower FA in the left optic radiation and left superior longitudinal fasciculus. Our data demonstrates that HA acclimatization is associated with brain structural modifications, including the loss of regional cortical grey matter accompanied by changes in the white matter, which may underlie the physiological adaptation of residents at HA.

Introduction

According to WHO (1996), there were approximately 140 million people living at high altitude (HA) over 2500 m. As of 2006, approximately 12 million Tibetan natives and Han lowland immigrants permanently reside between 2200 to 5200 m on the Qinghai–Tibetan Plateau, and every year hundreds of thousands of lowlanders traveled up to the Tibetan plateau [1]. A large amount of evidences have shown these natives and immigrants in HA environment (hypoxia, hypobaric press, UV rays from the sun, cold, and dehydration) have developed adaptive changes in respiratory and cardiovascular regulations, which directly related to oxygen transport [2][5].

The brain is the control centre of the body. At HA, through afferent feedback, the adaptation in the cardiovascular and respiratory systems may act on the control centers in the brain. On the other hand, as the central nervous system is highly oxidative, it inevitably suffers from hypoxic stress. The mental disturbances of chronic mountain sickness may be the strongest indicator of central nervous system failing in acclimatization to HA hypoxia [6]. Several researches on HA residents mainly focused on cerebral glucose metabolic rates [7] and cerebral autoregulation [8], [9]. Up to now, how the brain of residents in structure acclimatized to HA remains unclear.

Hypoxia is one of the most important environmental factors at HA. Laboratory observations on hypoxic animals and patients suffered from chronic obstructive pulmonary diseases (COPD) or obstructive sleep apnea (OSA) revealed increases of cerebral microvessel density in the hippocampus, the cerebellum, and the motor and the somatosensory cortices [10][12], the impairments of grey matter (GM) in the cerebellum, the parietal cortex, the anterior cingulate, the caudate, the putamen, the thalamus, the hippocampus, and the parahippocampus [13][16] and the lesions of white matter (WM) in the pons, the frontal, the temporal and the parietal cortices, the corpus callosum, projections to and from the cerebellum, and within the limbic system [17], [18]. All those studies suggested multiple regional brain changes could have been occurred in chronic hypoxia-exposed HA residents.

The prefrontal cortex, the insular cortex, and the cingulate cortex have been proved involved in cardiovascular control [19][21], and among those brain areas, anterior insular and anterior cingulate cortices play an important role in the unpleasantness of dyspnea [22], which often occurred during the adaptation to HA hypoxia. Recently, Paulus et al. [23] proposed the hypotheses that the anterior insular and cingulate cortices should be needed to process perturbation of the homeostatic balance in extreme environments. Aerobic capacity was strongly correlated with the right anterior insular GM density [24] and decrease in aerobic capacity has been well known in sea-level (SL) residents who were acclimatized to chronic HA hypoxia during the developmental period [25], [26]. Thus, a decreased GM in anterior insula at residents acclimatized to HA was suggested. Based on the above data, we expected that, in order to adapt changed peripheral physiology at HA environment, the modifications in brain may include these areas.

In this study, to investigate structural modifications in brain, residents in the Qinghai-Tibetan Plateau were examined using magnetic resonance imaging (MRI). Voxel-based morphometry (VBM) provides a quantitative and comprehensive assessment of anatomical differences throughout the brain [27]. It had been widely applied in numerous clinical and neuroscience researches. In the present study, we used VBM to reveal possible changes in GM volume. Diffusion tensor imaging (DTI) produces in vivo images of biological tissue. The anisotropy of the fibrous microstructure reflected from restricted water diffusion is measured with fractional anisotropy (FA) value. FA can reflect neural axons of the WM in the brain and used as a measure for WM integrity [28]. FA values have successful clinical application [29][31]. Current DTI data analyses of FA values include region-of-interest (ROI) analysis and voxel-wise analysis. We first used conventional regional analyses to characterize the pattern of anatomic alterations of HA effects by examining multiple, separate ROIs that spanned all brain structures. This anatomic structure based ROI analysis of FA values has been successful used to study clinical diseases [32][34]. Tract-Based Spatial Statistics (TBSS) is a new voxel-wise method proposed recently. It alleviates the alignment-related problems by applying both linear and nonlinear alignment to data and additionally projecting the FA values of individual subjects in given spatial locations to the thinned common “FA-skeleton” of major WM structures, thus improving sensitivity, objectivity, and interpretability of analysis of multi-subject DTI data [35], [36]. Moreover, diffusion tensor eigenvalues (longitudinal and radial diffusivities) 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 vs. obliquely aligned myelinated fibers [37]. TBSS has also been successful used to study clinical diseases [37]. In the present study, using TBSS, we measured FA and longitudinal and radial diffusivities to examine alterations in alignment of myelinated fibers in WM tracts in the whole brain.

The influence of hypoxic stress, and the organism's response to it, are greater during growth than during adulthood. The differences between the highland and lowland natives in their physiological performance and morphology are mostly due to adaptations acquired during the developmental period [38]. The study on the Bolivians of foreign ancestry acclimatized to high altitude since birth or during growth attained suggested that developmental acclimatization is important in the attainment of normal physiological functions at HA [39]. Therefore, HA young adult residents (range 17–22 yr) born and raised at HA were studied in the present study. All of their ancestors migrated from SL.

Results

Physiological characteristics

There were no significant differences in hemoglobin levels, circulating red blood cell count, blood pressure, and pulse rate between HA residents and SL controls. Body height of both males and females in HA residents were significantly larger than that of their SL controls (p = 0.006; p = 0.015, respectively). HA females had a higher diastolic pressure (p = 0.017) compared with SL females. No significant differences in pulmonary function and hematological measurements were found between HA residents and SL subjects (Table 1).

Table 1
Physiological Characteristics.

Average volume of GM, WM and cerebrospinal fluid (CSF)

No subject from either group showed visible abnormalities on T1-weighted structural images. There were also no significant differences in average volumes of the whole brain, the GM or the WM between the two groups; but HA residents showed a significant increase in CSF total volume (t = 2.302, p = 0.025) (Table 2).

Table 2
Average volume (mm3) of GM, WM and CSF in HA and SL groups (mean ± SD).

GM

VBM analysis showed that HA residents had decreased GM volume compared with SL controls in the bilateral anterior insula, the right anterior cingulate cortex, the bilateral prefrontal cortex, the left precentral cortex, and the right lingual cortex. (two-sample t-test, |t|>2.70, p<0.01, FEW corrected) (Figure 1). Coordinate information is shown in Table 3.

Figure 1
Gray matter volume decrease in HA residents vs. SL controls.
Table 3
Regional information of decreased gray matter volume in HA subjects compared with SL controls.

FA

ROI analysis

ROI analysis revealed significant increase of FA values in both the right and left anterior limb of internal capsule (ALIC) and significant decrease at the right posterior cingulum in HA residents compared with SL controls. No significant differences were detected in other areas. Mean FA values of each ROI for HA residents and SL controls are shown in Table 4 (two-sample t-test, p<0.05).

Table 4
Detailed regional FA value (mean ± SD) changes in HA subjects compared with SL controls.

Whole brain voxel-wise statistic analysis

Whole brain voxel-wise statistic analysis showed HA residents had significantly higher FA in a broad range of brain areas compared with SL controls (p<0.05, uncorrected) (Figure 2, Table 5). The significant regions (clusters size>100 voxels) included the right and left superior longitudinal fasciculus (SLF), the right and the left inferior longitudinal fasciculus (ILF), the corpus callosum (genu and body), the right and left superior corona radiata (SCR), the right and left ALIC, the right and left anterior external capsule, the right posterior cingulum and the right corticospinal tract. Conversely, compared with SL subjects, there were some regions with lower FA values in HA residents at the left optic radiation and left SLF (the inferior frontal lobe, the pars opercularis, the BA 44) (clusters size >8 voxels, Table 5).

Figure 2
Statistical maps of group comparison of FA values on a voxel-wise basis (results of TBSS).
Table 5
Main regions showing greater (clusters >100 voxels) and reduced (clusters >8 voxels) FA in HA subjects and SL controls.

Radial diffusivity

Higher FA (clusters >100 voxels) in the regions of the right and left SLF (frontal precentral cortex, frontal opercular cortex, parietal postcentral gyrus, inferior frontal pars triangularis), the right ILF (temporal supramarginal gyrus), the right and left SCR (frontal pole, BA 8, frontal precentral gyrus, parietal precuncus lobe, paracingulate gyrus, and postcentral gyrus), the right ALIC, the right and left external capsule (the anterior insula), and the corpus callosum (genu and body) were associated with decreased radial diffusivity in HA residents compared with SL controls. No significant differences were detected in the higher FA regions of the left SLF (the frontal lobe, BA 11 and the frontal Sub-Gyral lobe), the right and left ILF (the posterior thalamus, the temporal lobe and the temporal pianum temporle lobe), the right SCR (the frontal pole, BA 9), the left ALIC, the right posterior cingulatum, and the right corticospinal tract between the HA and SL groups. Lower FA values (clusters size >8 voxels) in the both left optic radiation (V5) and left SLF (inferior frontal lobe, BA 44) were associated with no changed radial diffusivity in HA subjects compared with SL controls. (Table 5)

Longitudinal diffusivity

All the regions have no significant changes in the longitudinal diffusivity between HA and SL groups (Table 5).

Discussion

Our present study first revealed that HA acclimatization was associated with brain structural modifications, which included the loss of regional cortical GM volume and WM structures. As we anticipated, the changed GMs were mainly confined in the prefrontal cortex, the anterior insular cortex, the anterior cingulate cortex, and the lingual cortex. Importantly, we found significant changes in anisotropy and diffusivity that reflect widespread alterations in fiber pathways. No significant differences were found in the average volume of GM, WM, and the whole brain of HA residents compared to SL controls, indicating that no global brain changes occurred during the acclimatization to HA. Higher CSF in HA residents may be a sign of a larger ventricle space in HA than in SL brain. Since our subjects had lived at SL for over one year, the present study suggested that the structural changes persisted even after HA residents had relocated at normoxia environment for a long time.

The implications for the loss of GM in the anterior insula, the anterior cingulate cortex, and the prefrontal cortex

Our present findings confirmed the hypotheses proposed by Paulus et al. [23] that the anterior insular and cingulate cortices should be needed to process perturbation of the homeostatic balance in extreme environments and agreed with the results of MRI studies on hypoxic OSA patients who showed a loss of anterior cingulate GM [13], [18].

The insular cortex is connected with the hypothalamic, the midbrain, the pontine, and the medullary brain regions that are involved in cardiovascular control and receives visceral sensory information arising from baroreceptors and chemoreceptors within the cardiovascular system [19]. The information is then relayed to the anterior cingulate cortex [23]. Electrical cortical stimulation and human neuroimaging studies have demonstrated that the activation of the anterior insula, the dorsal anterior cingulate cortex and the prefrontal cortex, modulated sympathetic nerve activity, heart rate, and blood pressure [21], [24], [40][42]. Dyspnea, often occurred during the adaptation to HA, activated the insular cortex, which had been shown in the inspiratory loading and volitional breathing tasks [2], [22], [43], [44]. Therefore, elevations in the sympathoadrenal system in the SL natives temporarily exposed to HA may be due to the activity of the anterior insula and its connected brain areas. HA hypoxia was also reported to have induced activation of the carotid baroreceptors, but it was not sufficient to completely counteract the catecholamine-induced increase in blood pressure and peripheral vasoconstriction [45]. However, after some years of residence the blood pressure tends to gradually decline, even falling below that observed at SL [46]. The decreased GM volume in the anterior cingulate and the insular cortices that found in the present study might be responsible for the decreased blood pressure. As the lesion to the insula has been proved to disrupt the representation of internal states that underpinning motivation [47], reduce in the anterior insular GM may be a reason why HA residents had a blunted hypoxic ventilatory response [48].

The HA residents have not shown significant differences with their SL peers in the peripheral physiology after they had relocated SL for over one year, which suggested pulmonary and cardiovascular functions had readapted to SL after acclimatization to HA [2], [3]. However, possibly due to failure in neuronal regenerate in SL normoxia, the loss of brain GM persisted. This may explain why HA residents persisted blunting of the ventilatory response to hypoxia [49] and had the lesser degree of sympathetic activation [45] despite long-term acclimatization at SL. One shortcoming of the present study was that we did not test hypoxic ventilatory response.

The loss of GM in the anterior insula, the prefrontal cortex, and the anterior cingulate cortex found in the present study may clarify the mechanisms underlining the decreased aerobic capacity [24][26] and decreased appetite [50] and may suggest changed cognitive functions [51][53] in HA acclimatized immigrants.

The implications for the changes in WM

In the present study, although anatomic structure based ROI analysis only detected changed FA confined to a few areas, changed FA values were found in a broad areas of WM using TBSS analysis. One of the most interesting findings of our study is the changes in FA values were symmetric in most regions between the left and right hemispheres in HA residents except the superior longitudinal fasciculus. The symmetric FA increase in both hemispheres implies the equal impacts of HA on the left and right sides of the brains. The discrepant findings between ROI and TBSS analyses may be due to that ROI analysis is difficult to objectively and reproducibly place ROIs on small or thin tracts on the images of individual subjects, when the slice orientation and anatomical details may show variation between individuals and the boundaries of the WM tracts are not easily identified. Moreover, the FA value obtained from ROI analysis is largely affected by location and size of the ROIs; both being limiting factors in accuracy of ROI method since the ROI selection is often conducted without a prior knowledge about the exact location [54]. However, both two analysis methods consistently showed significantly higher FA in the bilateral ALIC.

Although the voxel-based analysis provides the ability to map the FA changes within the whole brain at once without a priori define ROIs it has been shown that the interpretability of the results from voxel-based DTI may sometimes be questionable [55], [56]. The TBSS method used here overcomes the problems in the voxel-based DTI analysis in regards of alignment and spatial smoothing, and thus seems to be a far more accurate voxelwise analysis method when comparing multi-subject diffusion data [35]. Greater FA may reflect greater myelination of WM fibers, increased number of myelinated fibers, smaller axonal diameter, or reduced neural branches within MRI voxel [37], [57], [58]. Whereas reduced FA was associated with local cerebral edema, cerebrospinal fluid, compromised myelin structure, changes in axonal morphologic structure, and altered interaxonal spacing of fiber bundles [57], . In our study, further analysis found that most of those regions were associated with reduced radial diffusivity coinciding with no change of longitudinal diffusivity. Radial diffusivity was positively correlated with the mean axon diameter while the longitudinal diffusivity was negatively correlated with it [63]. Greater FA values associated with reduced radial diffusivity coinciding with no change of longitudinal diffusivity can result from increased myelin thickness, smaller axonal diameter or extracellular space [64].

The ALIC contains the anterior thalamic peduncle, which connects the dorsomedial and anterior thalamic nuclei with the prefrontal cortex and the cingulate cortex. From a functional neuroanatomic perspective, ALIC is involved in the medial limbic circuit (composed of the hippocampal formation, the mammillary bodies, the anterior thalamic nuclei, and the cingulate gyrus) and the basolateral limbic circuit (interconnecting the orbitofrontal cortex, the dorsomedial thalamic nucleus, the amygdala, and the anterior temporal cortex) [65]. The increased FA values with reduced radial diffusivity and no change of longitudinal diffusivity in ALIC described in our study indicated that HA acclimatization was involved with improvement of the fronto-thalamic or cingulate-thalamic structural connectivity, which may result in higher functional connectivity between the cortical cortex and subcortical regions partially linked by these circuits.

The changed FA in the posterior cingulate fiber tract may alter communication between components of the Papez circuit in HA residents [66]. The body of the corpus callosum contains fibres important for connecting motor and sensory cortices. The radiating bundle of corona radiata fibres generates ascending paths to the motor cortices and descending paths such as the corticospinal tract to spinal motor neurons. Therefore, increases of FA at those locations may be linked to improvement of motor skills in acclimatization to HA.

The mechanisms that involved in brain structural modulations

The structural changes in the brain of HA residents may be related to the decrease in glucose metabolism [7] and/or impaired cerebral autoregulation [8], [9]. Moreover, GM loss mostly occurred in new cortex in HA residents found in the present study may be due to the unevenly distributed cerebral blood flow and cerebral blood volume during baseline conditions. Previous study revealed that phylogenetically older regions of the brain, which receives greater than average increases in cerebral blood flow, showed larger vascular responses to hypoxia than evolutionary younger regions that generally received below average increases [67].

What determined the brain structural modifications: gene or developmental environment?

The brain structural modifications found in our subjects who ancestors migrated from SL and have resident at HA at least two generations may occur in genetically adapted HA natives, such as Tibetan, and in adolescents who immigrated to HA with parent after their birth. A recently study of Hogan et al. [68] may support our hypothesis. Hogan et al.([68] tested a reduced psychomotor speed in children (6–10 years) and adolescents (13–16 years) of mixed-ethnic background (including Native American, European or African ancestry whose ancestor immigrate to HA during different periods) who were born and raised at HA, they found the reduced psychomotor speed was well correlated with the reduced cerebral metabolism and blood flow. For the proportion of European, Native American and African genetic admixture was comparable across altitude groups, therefore, chronic hypoxic exposure rather than genetic inheritance appears to affect the neurocognitive development at HA. However, Brutsaert et al. [69] found that pulmonary function (forced vital capacity and forced expiratory volume) measures were larger in HA natives compared with low altitude natives born and raised at HA, suggesting a genetic effect. In contrast, forced vital capacity and forced expiratory volume were similar in HA natives and SL natives at low altitude, suggesting that the genetic potential for larger lung volumes at HA depended upon developmental exposure to HA. In summary, results from those studies emphasized the importance of developmental adaptation to HA.

As our present findings in white matter change regions are in agreement with most of regions showing FA increases during adolescence [70], we speculate that higher FA with decreased radial diffusivity coinciding with no change of longitudinal diffusivity is a developmental modification in fiber microstructure in acclimatization to HA. Recently, through longitudinal studies, Bava et al. [70] have documented linear increases in FA with decreases in radial diffusivity across typical adolescent development continuing through the second decade of life, showing FA increases in the bilateral SLF, SCR, thalamic radiations, posterior internal capsule, corticospinal tract, arcuate fasciculus, superior and mid-temporal white matter, inferior parietal white matter, and the corpus callosum. Future study should be done to explore how the brain acclimatized to HA in residents who immigrated from SL to HA at adult and have been resided at HA for many years.

In summary, we demonstrated that HA acclimatization was associated with regional brain structural modifications, which mainly related to cardiovascular or respiratory regulations. These regional brain changes may underlie the physiological functions of HA residents at high altitude and indicate some cognition and motor skill deficits or enhancement. These changes persisted even after HA inhabitants relocated to a SL residence. Future study is needed to explore the functional modifications of the brain in acclimatization to HA and structural brain modulations in relation to the behaviors and physiological functions of HA residents.

Materials and Methods

Subjects

Subjects in the current study consisted of 28 HA immigrant residents (mean age 20.4 yr; range 17–22 yr) born and raised at the altitude of 2616–4200 m in Qinghai-Tibetan Plateau for at least 17 years and currently attended college at Chendu (<400 m), and their ancestors migrated from SL and have resident at HA at least two generations (Table 6). The 28 matched control subjects (mean age 20.9 yr; range 17–23 yr) were their schoolmates, all of whom were lowlanders, born and living at sea level below 400 m and without any prior exposure to high altitude. HA and SL students did not differ in IQ or college enrollment scores and all of them were selected from Han populations to avoid possible racial differences. Physiological and MRI studies were examined at Chendu. Subjects were excluded if they had: (1) chronic mountain sickness, (2) a documented neurological disorder, or (3) a past history of head injury with loss of consciousness. Procedures were fully explained, and all subjects provided written informed consent before participating in the study. The experimental protocol was approved by the Research Ethics Review Board of the Institute of Psychology, Chinese Academy of Sciences.

Table 6
Demographic information of HA subjects and SL controls (based on self-reported).

Physical and physiological assessments

Before MRI scanning, subjects underwent physical and physiological examinations, including body weight, hematological measurements, blood pressure, pulse rate, and pulmonary function. Hematological measurements were tested using a hematology analyzer (Sysmex XE-2100, TOA Medical Electronics, Kobe, Japan). Pulmonary function was tested using a pulmonary function testing device (Master Screen Body, JAEGER, German).

MRI data acquisition

Structural images were acquired on a GE 3.0 T Signa Excite Gemse MRI system (GE Medical, Milwaukee, WI, USA) at Huaxi Magnetic Resonance Research Center (West China Hospital, Chengdu, China).

A 3D structural MRI was acquired from each subject using a T1-weighted MPRAGE sequence (TR/TE = 8.5 ms/3.4 ms, TI = 400 ms, FOV = 280×280 mm2, in-plane resolution = 0.547×1.094 mm2, flip angle = 12°), yielding 156 contiguous axial slices (1 mm thick) covering the whole brain. A DTI pulse sequence with single shot diffusion-weighted echo planar imaging (TR/TE = 10000/70.8 ms, FOV = 240×240 mm2, in-plane resolution = 1.875×1.875 mm2) was applied sequentially in 16 different directions. We acquired 42 contiguous 3-mm thick slices covering the whole brain.

VBM analysis of 3D T1 images

The 3D T1 images were used for GM volume analysis using VBM implemented in SPM2 toolbox (Wellcome Department of Imaging Neuroscience, London). For each section of GM, WM and cerebral spinal fluid compartments, we constructed an anatomical template map as well as a probability map specific to the current study based on the 3D images from all the subjects. Using these study-specific templates, the 3D images for each individual were spatially normalized to Talairach space [71] and segmented, and then were smoothed using a Gaussian kernel of 10 mm full-width at half-maximum (FWHM). Random-effect two-sample t-tests were performed to examine between-group differences. The statistical parametric map was generated with threshold at |t|>2.70, p<0.01 (FWE corrected).

Whole brain voxel-wise statistic analysis and ROI analyses of DTI

ROI analysis

All analysis were conducted with the DTI Studio software version 2.40 provided by Johns Hopkins University [72] (https://www.dtistudio.org/) on the FA image of each participant. The appropriate slice on which the ROI was identified was chosen by using the anatomical landmarks from the FA maps. The ROIs were then drawn manually by an experienced operator blind to subjects' status, using standardized guidelines based on location and size. ROIs were placed in a total of 10 areas. The ROIs were depicted in Figure 3 and defined as Table 7. Data were further analyzed using SPSS. ANOVA statistic identified the differences between the HA and SL groups. Data were presented as mean±SD. Statistical significance was set at p<0.05.

Figure 3
Examples of ROI demarcation on FA images.
Table 7
Regional definitions for FA values analyses (Depicted in Figure 3).

Whole brain voxel-wise statistic analysis

We used DCM2MII to convert diffusion tensor images from the proprietary scanner format to the NIFTI format. Then images were processed using FSL 4.1.5 software package (http://www.fmrib.ox.ac.uk/fsl/). Images were realigned to the b-value (b0) image to remove eddy current distortions and motion artifacts using FDT (FMRIB's diffusion toolbox) [73]. Brain mask was created from the first b0 image using BET (Brain Extraction Tool). After those processes images were calculated with the FDT for FA and longitudinal diffusivity and radial diffusivity maps. The analysis of FA images was performed using the TBSS package in FSL [35], [36]. TBSS processing includes the following steps: (1) Align the FA images of all subjects to a template which was arbitrarily selected from those FA images by nonlinear registrations; (2) Transform all the aligned FA images into 1×1×1 mm3 MNI152 space by affine registrations; (3) Create the mean FA image and filter to retain only the center of the WM tracts so as to create the mean FA skeleton; (4) Project individual subjects' FA was put onto the skeleton. (5) Following these steps, data were fed into voxel-wise cross-subject statistical analyses with the following group comparisons: HA vs. SL and SL vs. HA. In all cases, the null distribution was built up over 5000 permutations, and significance was tested at p<0.05 levels, uncorrected for multiple comparisons. We determined the anatomic localization of each cluster by means of the FSL atlas tool, which incorporates several anatomic templates, including the Talairach atlas, MNI structural atlas, Julich histological atlas, Harvard-Oxford cortical and subcortical structural atlases, and the Johns Hopkins University DTI-based WM atlases.

Longitudinal (principal diffusion direction, λ1) and radial (transverse diffusion component, [(λ2+λ3)/2] diffusivity values were computed for clusters showing a significant FA change (increase or decrease) between the HA residents and SL controls. Data were analyzed using SPSS. ANOVA statistic identified the differences between the HA and SL groups. Data were presented as mean ± SD. Statistical significance was set at p<0.05.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was funded by National Science Foundation of China (Project No. 30425008, No. 60628101), Ministry of Science and Technology (Project No. 2007CB512300) and China Postdoctoral Science Foundation (Project No. 20060390129). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Wu T, Kayser B. High altitude adaptation in Tibetans. High Alt Med Biol. 2006;7:193–208. [PubMed]
2. Zhuang J, Droma T, Sun S, Janes C, McCullough RE, et al. Hypoxic ventilatory responsiveness in Tibetan compared with Han residents of 3,658 m. J Appl Physiol. 1993;74:303–311. [PubMed]
3. Curran LS, Zhuang J, Sun SF, Moore LG. Ventilation and hypoxic ventilatory responsiveness in Chinese-Tibetan residents at 3,658 m. J Appl Physiol. 1997;83:2098–2104. [PubMed]
4. Beall CM. Two routes to functional adaptation: Tibetan and Andean high-altitude natives. Proc Natl Acad Sci USA. 2007;104 S1:8655–8660. [PubMed]
5. Penaloza D, Arias-Stella J. The heart and pulmonary circulation at high altitudes: healthy highlanders and chronic mountain sickness. Circulation. 2007;115:1132–1146. [PubMed]
6. Moore LG, Niermeyer S, Vargas E. Does chronic mountain sickness (CMS) have perinatal origins? Respir Physiol Neurobiol. 2007;158:180–189. [PubMed]
7. Hochachka PW, Clark CM, Brown WD, Stanley C, Stone CK, et al. The brain at high altitude: hypometabolism as a defense against chronic hypoxia? J Cereb Blood Flow Metab. 1994;14:671. [PubMed]
8. Jansen GF, Krins A, Basnyat B, Odoom JA, Ince C. The role of the altitude level on cerebral autoregulation in man resident at high altitude. J Appl Physiol. 2007;103:518–523. [PubMed]
9. Claydon VE, Gulli G, Slessarev M, Appenzeller O, Zenebe G, et al. Cerebrovascular responses to hypoxia and hypocapnia in Ethiopian high altitude dwellers. Stroke. 2008;39:336–342. [PubMed]
10. LaManna JC, Vendel LM, Farrell RM. Brain adaptation to chronic hypobaric hypoxia in rats. J Appl Physiol. 1992;72:2238–2243. [PubMed]
11. Boero JA, Ascher J, Arregui A, Rovainen C, Woolsey TA. Increased brain capillaries in chronic hypoxia. J Appl Physiol. 1999;86:1211–1219. [PubMed]
12. Kanaan A, Farahani R, Douglas RM, LaManna JC, Haddad GG. Effect of chronic continuous or intermittent hypoxia and reoxygenation on cerebral capillary density and myelination. Am J Physiol Regul Integr Comp Physiol. 2006;290:R1105–1114. [PubMed]
13. Macey PM, Henderson LA, Macey KE, Alger JR, Frysinger RC, et al. Brain morphology associated with obstructive sleep apnea. Am J Respir Crit Care Med. 2002;166:1382–1387. [PubMed]
14. Morrell MJ, McRobbie DW, Quest RA, Cummin AR, Ghiassi R, et al. Changes in brain morphology associated with obstructive sleep apnea. Sleep Med. 2003;4:451–454. [PubMed]
15. Morrell MJ, Twigg G. Neural consequences of sleep disordered breathing: the role of intermittent hypoxia. Adv Exp Med Biol. 2006;588:75–88. [PubMed]
16. Huang BY, Castillo M. Hypoxic-ischemic brain injury: imaging findings from birth to adulthood. Radiographics. 2008;28:417–439. [PubMed]
17. Wideroe M, Olsen O, Pedersen TB, Goa PE, Kavelaars A, et al. Manganese-enhanced magnetic resonance imaging of hypoxic-ischemic brain injury in the neonatal rat. Neuroimage. 2009;45:880–890. [PubMed]
18. Macey PM, Kumar R, Woo MA, Valladares EM, Yan-Go FL, et al. Brain Structural Changes in Obstructive Sleep Apnea. Sleep. 2008;31:967–977. [PubMed]
19. Verberne AJ, Owens NC. Cortical modulation of the cardiovascular system. Prog Neurobiol. 1998;54:149–168. [PubMed]
20. Green AL, Paterson DJ. Identification of neurocircuitry controlling cardiovascular function in humans using functional neurosurgery: implications for exercise control. Exp Physiol. 2008;93:1022–1028. [PubMed]
21. Wager TD, Waugh CE, Lindquist M, Noll DC, Fredrickson BL, et al. Brain mediators of cardiovascular responses to social threat: part I: Reciprocal dorsal and ventral sub-regions of the medial prefrontal cortex and heart-rate reactivity. Neuroimage. 2009;47:821–835. [PMC free article] [PubMed]
22. Davenport PW, Vovk A. Cortical and subcortical central neural pathways in respiratory sensations. Respir Physiol Neurobiol. 2009;167:72–86. [PubMed]
23. Paulus MP, Potterat EG, Taylor MK, Van Orden KF, Bauman J, et al. A neuroscience approach to optimizing brain resources for human performance in extreme environments. Neurosci Biobehav Rev. 2009;33:1080–1088. [PMC free article] [PubMed]
24. Peters J, Dauvermann M, Mette C, Platen P, Franke J, et al. Voxel-based morphometry reveals an association between aerobic capacity and grey matter density in the right anterior insula. 2009. Neuroscience In Press. [PubMed]
25. Frisancho AR, Martinez C, Velasquez T, Sanchez J, Montoye H. Influence of developmental adaptation on aerobic capacity at high altitude. J Appl Physiol. 1973;34:176–180. [PubMed]
26. Marconi C, Marzorati M, Grassi B, Basnyat B, Colombini A, et al. Second generation Tibetan lowlanders acclimatize to high altitude more quickly than Caucasians. . J Physiol. 2004;556(Pt 2):661–671. [PubMed]
27. Ashburner J, Friston KJ. Voxel-based morphometry-the methods. Neuroimage. 2000;11:805–821. [PubMed]
28. Gulani V, Sundgren PC. Diffusion tensor magnetic resonance imaging. J Neuroophthalmol. 2006;26:51–60. [PubMed]
29. Beppu T, Inoue T, Shibata Y, Yamada N, Kurose A, et al. Fractional anisotropy value by diffusion tensor magnetic resonance imaging as a predictor of cell density and proliferation activity of glioblastomas. Surg Neurol. 2005;63:56–61. [PubMed]
30 . Sidaros A, Engberg AW, Sidaros K, Liptrot MG, Herning M, et al. Diffusion tensor imaging during recovery from severe traumatic brain injury and relation to clinical outcome: a longitudinal study. Brain. 2008;131:559. [PubMed]
31. Chua TC, Wen W, Slavin MJ, Sachdev PS. Diffusion tensor imaging in mild cognitive impairment and Alzheimer's disease: a review. Curr Opin Neurol. 2008;21:83–92. [PubMed]
32. Head D, Buckner RL, Shimony JS, Williams LE, Akbudak E, et al. Differential vulnerability of anterior white matter in nondemented aging with minimal acceleration in dementia of the Alzheimer type: evidence from diffusion tensor imaging. Cereb Cortex. 2004;14:410–423. [PubMed]
33. Zhang Y, Schuff N, Jahng GH, Bayne W, Mori S, et al. Diffusion tensor imaging of cingulum fibers in mild cognitive impairment and Alzheimer disease. Neurology. 2007;68:13–19. [PMC free article] [PubMed]
34. Chang L, Wong V, Nakama H, Watters M, Ramones D, et al. Greater than age-related changes in brain diffusion of HIV patients after 1 year. J Neuroimmune Pharmacol. 2008;3:265–274. [PMC free article] [PubMed]
35. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage. 2006;31:1487–1505. [PubMed]
36. Smith SM, Johansen-Berg H, Jenkinson M, Rueckert D, Nichols TE, et al. Acquisition and voxelwise analysis of multi-subject diffusion data with tract-based spatial statistics. Nat Protocols. 2007;2:499–503. [PubMed]
37. Versace A, Almeida JR, Hassel S, Walsh ND, Novelli M, et al. Elevated left and reduced right orbitomedial prefrontal fractional anisotropy in adults with bipolar disorder revealed by tract-based spatial statistics. Arch Gen Psychiatry. 2008;65:1041–1052. [PMC free article] [PubMed]
38. Frisancho AR. Developmental adaptation to high altitude hypoxia. Int J Biometeorol. 1977;21:135–146. [PubMed]
39. Frisancho AR, Frisancho HG, Milotich M, Brutsaert T, Albalak R, et al. Developmental, genetic, and environmental components of aerobic capacity at high altitude. Am J Phys Anthropol. 1995;96:431–442. [PubMed]
40. Burns SM, Wyss JM. The involvement of the anterior cingulate cortex in blood pressure control. Brain Res. 1985;340:71–77. [PubMed]
41. Critchley HD, Mathias CJ, Josephs O, O'Doherty J, Zanini S, et al. Human cingulate cortex and autonomic control: converging neuroimaging and clinical evidence. Brain. 2003;126:2139–2152. [PubMed]
42. Kimmerly DS, Wong SW, Salzer D, Menon R, Shoemaker JK. Forebrain regions associated with postexercise differences in autonomic and cardiovascular function during baroreceptor unloading. Am J Physiol Heart Circ Physiol. 2007;293:H299–306. [PubMed]
43. Frisancho AR, Juliao PC, Barcelona V, Kudyba CE, Amayo G, et al. Developmental components of resting ventilation among high- and low-altitude Andean children and adults. Am J Phys Anthropol. 1999;109:295–301. [PubMed]
44. von Leupoldt A, Dahme B. Cortical Substrates for the Perception of Dyspnea. Chest. 2005;128:345–354. [PubMed]
45. Bernardi L, Passino C, Spadacini G, Calciati A, Robergs R, et al. Cardiovascular autonomic modulation and activity of carotid baroreceptors at altitude. Clin Sci (Lond) 1998;95:565–573. [PubMed]
46. Tripathy V, Satapathy KC, Gupta R. Nutritional Status and Hypertension among Tibetan Adults in India. Human Ecology. 2006;14:77–82.
47. Dani JA, Montague PR. Disrupting addiction through the loss of drug-associated internal states. Nat Neurosci. 2007;10:403–404. [PubMed]
48. Frisancho AR. Functional adaptation to high altitude hypoxia. Science. 1975;187(4174):313–319. [PubMed]
49. Rivera-Ch M, Gamboa A, Leon-Velarde F, Palacios JA, O'Connor DF, et al. Selected contribution: High-altitude natives living at sea level acclimatize to high altitude like sea-level natives. J Appl Physiol. 2003;94:1263–1268. [PubMed]
50. Westerterp-Plantenga MS, Westerterp KR, Rubbens M, Verwegen CR, Richelet JP, et al. Appetite at “high altitude” [Operation Everest III (Comex-'97)]: a simulated ascent of Mount Everest. J Appl Physiol. 1999;87:391–399. [PubMed]
51. Constantinidis C, Procyk E. The primate working memory networks. Cogn Affect Behav Neurosci. 2004;4:444–465. [PubMed]
52. Silton RL, Heller W, Towers DN, Engels AS, Spielberg JM, et al. The time course of activity in dorsolateral prefrontal cortex and anterior cingulate cortex during top-down attentional control. Neuroimage. 2010;50:1292–1302. [PubMed]
53. Hyafil A, Summerfield C, Koechlin E. Two mechanisms for task switching in the prefrontal cortex. J Neurosci. 2009;29:5135–5142. [PubMed]
54. Giuliani NR, Calhoun VD, Pearlson GD, Francis A, Buchanan RW. Voxel-based morphometry versus region of interest: a comparison of two methods for analyzing gray matter differences in schizophrenia. Schizophr Res. 2005;74:135–147. [PubMed]
55. Bookstein FL. “Voxel-based morphometry” should not be used with imperfectly registered images. Neuroimage. 2001;14:1454–1462. [PubMed]
56. Davatzikos C. Why voxel-based morphometric analysis should be used with great caution when characterizing group differences. Neuroimage. 2004;23:17–20. [PubMed]
57. Beaulieu C. The basis of anisotropic water diffusion in the nervous system- a technical review. NMR Biomed. 2002;15:435–455. [PubMed]
58. Hoeft F, Barnea-Goraly N, Haas BW, Golarai G, Ng D, et al. More is not always better: increased fractional anisotropy of superior longitudinal fasciculus associated with poor visuospatial abilities in Williams syndrome. J Neurosci. 2007;27:11960–11965. [PubMed]
59. Mukherjee P, Miller JH, Shimony JS, Philip JV, Nehra D, et al. Diffusion-tensor MR imaging of gray and white matter development during normal human brain maturation. AJNR Am J Neuroradiol. 2002;23:1445–1456. [PubMed]
60. Schwartz ED, Hackney DB. Diffusion-weighted MRI and the evaluation of spinal cord axonal integrity following injury and treatment. Exp Neurol. 2003;184:570–589. [PubMed]
61. Concha L, Beaulieu C, Gross DW. Bilateral limbic diffusion abnormalities in unilateral temporal lobe epilepsy. Ann Neurol. 2005;57:188–196. [PubMed]
62. Thomalla G, Glauche V, Koch MA, Beaulieu C, Weiller C, et al. Diffusion tensor imaging detects early Wallerian degeneration of the pyramidal tract after ischemic stroke. Neuroimage. 2004;22:1767–1774. [PubMed]
63. Barazany D, Basser PJ, Assaf Y. In vivo measurement of axon diameter distribution in the corpus callosum of rat brain. Brain. 2009;132:1210–1220. [PMC free article] [PubMed]
64. Gao W, Lin W, Chen Y, Gerig G, Smith JK, et al. Temporal and spatial development of axonal maturation and myelination of white matter in the developing brain. AJNR. Am J Neuroradiol. 2009;30:290–296. [PMC free article] [PubMed]
65. Wobrock T, Kamer T, Roy A, Vogeley K, Schneider-Axmann T, et al. Reduction of the internal capsule in families affected with schizophrenia. Biol Psychiatry. 2008;63:65–71. [PubMed]
66. Papez JW. A proposed mechanism of emotion. 1937 J Neuropsychiatry Clin. Neurosci. 1995;7:103–112. [PubMed]
67. Binks AP, Cunningham VJ, Adams L, Banzett RB. Gray matter blood flow change is unevenly distributed during moderate isocapnic hypoxia in humans. J Appl Physiol. 2008;104:212–217. [PubMed]
68. Hogan AM, Virues-Ortega J, Botti AB, Bucks R, Holloway JW, et al. Development of aptitude at altitude. Dev Sci. 2010;13:533–544. [PubMed]
69. Brutsaert TD, Soria R, Caceres E, Spielvogel H, Haas JD. Effect of developmental and ancestral high altitude exposure on chest morphology and pulmonary function in Andean and European/North American natives. Am J Hum Biol. 1999;1:383–395. [PubMed]
70. Bava S, Thayer R, Jacobus J, Ward M, Jernigan TL, et al. Longitudinal characterization of white matter maturation during adolescence. Brain Res. 2010;1327:38–46. [PMC free article] [PubMed]
71. Talairach J, Tournoux PA. New York: Thieme Press; 1988. A co-planar stereotactic atlas of the human brain.
72. Jiang H, van Zijl PC, Kim J, Pearlson GD, Mori S. DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking. Comput Methods Prog Biomed. 2006;81:106–116. [PubMed]
73. Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal. 2001;5:143–156. [PubMed]
74. Witelson SF. Hand and sex differences in the isthmus and genu of the human corpus callosum. A postmortem morphological study. Brain. 1989;112:799–835. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science