Age-related changes in neural volume and microstructure associated with interleukin-6 are ameliorated by a calorie-restricted diet in old rhesus monkeys

Peripheral and central production of proinflammatory cytokines (PIC), including interleukin-6 (IL-6), upregulates oxidative stress in brain and increases age-related neural atrophy and risk for neurodegenerative pathologies ( Ershler, 1993 ; Kagiwada et al., 2004 ; Baranowska-Bik et al., 2008 ; Qin et al., 2008 ; Chen et al., 2008 ; Baune et al., 2009 ). IL-6 concentrations in sera and spleen typically increase with age ( Ferrucci et al., 1999 ; Ershler & Keller, 2000 ; Bruunsgaard et al., 2001 ) and age-related chronic inflammatory states ( Brouqui et al. 1994 ; Devaux et al. 1997 ; De Keyser et al. 1998 ), which may prime microglia to overproduce PIC and contribute to neuronal degradation ( Godbout et al., 2005 ; Deng et al., 2006; Cunningham et al., 2007 ; Kumagai et al., 2007 ; Perry et al., 2007 ; Teeling and Perry, 2008). In humans, peripheral IL-6 is related to lower global and regional gray and white matter volume (GM, WM; Hauss-Wegrzyniak et al., 2000 ; Jefferson et al., 2007 ; Marsland et al., 2008 ). Aging is similarly related to diffuse regional atrophy ( Sowell et al., 2003 ), and therefore it was of interest to further investigate the relationship between neural indices and systemic IL-6 using a whole brain voxel-based approach in an aged animal model, the rhesus monkey ( Macaca mulatta ), that undergoes some of the age-related brain atrophy seen in humans ( Roth et al., 2004 ). IL-6 is uniquely suited for conducting these analyses and translating results to humans. The protein structure of IL-6 shares greater than 98% homology between monkeys and humans, and there is a strong relationship between peripheral and central IL-6 levels in rhesus monkeys ( Reyes and Coe, 1996 ). In addition, because IL-6 is produced by many types of cells in the body, it is a stable biomarker of many physiological systems ( Rao et al., 1994 ; Willette et al., 2007 ). By contrast, other PIC such as interleukin-1 and tumor necrosis factor-alpha are quite low in the blood stream in the absence of infection. IL-6 is also more readily measured in the healthy individual, and is widely used as a biomarker for chronic, low-grade systemic inflammation typically associated with aging (Bruunsgaard, 2002) and many disease states (Fonseca et al., 2009). Neural atrophy related to IL-6 levels may be reduced by a calorically restricted diet (CR), which downregulates mRNA and protein expression of PIC ( Spaulding et al., 1997 ; Arvidsson et al., 2004 ; You et al., 2007 ). CR also decreases Alzheimer’s disease-like neural sequelae in animal models ( Qin et al., 2006 ; Halagappa et al., 2007 ), and can lower proinflammatory or oxidative stress factors in aged murine neocortex ( Lee et al., 2000 ) and hippocampus ( Wu et al., 2008 ). CR-induced reductions in visceral adiposity, body weight, and inflammatory pathologies may mediate decreases in peripheral IL-6 and were determined as possible mediating health factors ((Bruunsgaard, 2002; Bulcão et al., 2006 ; Colman et al., 2009 ). Finally, to confirm that IL-6 is a stable biomarker over time ( Rao et al., 1994 ), this protein was assessed twice: when the monkey reached 20 years of age and later at the time of the scan. This cross-sectional study used animals from a longitudinal CR project initiated in 1989 ( Ramsey et al., 2000 ). Voxel-wise analyses of segmented T1-weighted and diffusion tensor imaging (DTI) data were performed. We hypothesized that higher IL-6 levels acquired proximal to the time of the MRI scan would be associated with: 1) decreased global GM and WM ( Pang et al., 2003 ; Jefferson et al., 2007 ); and 2) decreased regional volumes, as well as lower fractional anisotropy (FA) and increased mean diffusivity (MD), in diffuse areas and tracts sensitive to age-related atrophy including: cortical and intrahemispheric parietal areas, superior temporal sulcus, corpus callosum, thalamus, the dorsal convexity post central sulcus, cingulate cortex, insula, and dorsal or medial portions of prefrontal cortex (PFC; Hauss-Wegrzyniak et al., 2000 ; Good et al., 2001 ; Sowell et al., 2003 ; Pfefferbaum et al., 2005 ; Sullivan et al., 2006 ; Kochunov et al., 2007 ; Pagani et al., 2008 ; Alexander et al., 2008 ; Marsland et al., 2008 ). FA is typically used to examine intravoxel tract coherence sensitive to changes in myelin and axons, while MD can broadly reflect changes in WM and GM tissue density. CR was expected to lower IL-6 and reduce associated neural atrophy. Further, Colman and colleagues (2009) recently showed that CR interacts with age to protect against atrophy. A complementary set of interaction analyses was thus conducted to assess the extent to which CR monkeys showed less atrophy per picogram/milliliter (pg/ml) increase in IL-6.

The subject composition of the two dietary conditions did not differ by age or sex. CR monkeys had lower basal IL-6 at the scan time point than did the age- and sex-matched controls [Z(45) = -2.201, p < .05]. There was a median difference of 1.12 pg/ml between the CR and control groups (Median =1.43 pg/ml, IQR=2.05, versus 2.55 pg/ml, IQR=2.43, respectively) ( Figure 1 ). For the healthy monkeys without pathologies at the time of scan, IL-6 was correlated with early values obtained at 20 years of age [r(29) = 0.47, p < .05]. By contrast, monkeys with some type of inflammatory pathology near time of scan did not show this association, which also negated the stability of the IL-6 values over time for all subjects [r(45) = -0.01, N.S.]. The two diet conditions had statistically similar IL-6 values at 20 years of age (Z(45) = -.754, N.S.). Volumetric and microstructural result maps related to IL-6 were comparable whether cytokine data were used from 20 years of age and at time of scan independently or averaged together (data not shown). The presence of inflammation-related conditions was similar in both CR (n=8) and controls (n = 7) in this sub-sample from the larger project. A binary variable representing this marker of health was correlated with higher IL-6 concentrations near the time of scan [r(45) = 0.35, p < .05] but not at 20 years of age [r(45) = 0.13, N.S]. The two monkeys with the highest IL-6 levels at the time of the scan had been diagnosed with chronic endometriosis. This gynecological condition is common in older female monkeys; therefore, they were included in all neural analyses as representative of normative aging in this species. Central adiposity and total body mass were not correlated with IL-6 sampled near the time of scan. The TBV of the control and CR monkeys did not differ significantly. Higher peripheral IL-6 near scanning was moderately correlated with lower TBV [r = -.41, p < .01]. This relationship appeared stronger for global GM [r = -.52, p <.001] as compared to WM [r = -.34, p < .05]. IL-6 was associated with reduced GM volume in parietal and temporal regions sensitive to age-related atrophy. Several clusters were found in association cortices along the intraparietal sulcus; medial parietal cortex; and superior temporal sulcus in lateral, medial, and superior temporal areas ( Figure 2a-c ; Table 1 ). GM regions typically not associated with age-related atrophy included primary visual cortex and extrastriate areas. There was no significant relationship between IL-6 levels and regional WM volume after correction for multiple comparisons. An IL-6*CR Status model examined if CR monkeys had more volume than controls across subjects per pg/ml increase in IL-6, suggesting a protective effect against IL-6 related atrophy. For GM volume, no clusters were significant after Type 1 error correction. For WM volumes ( Figure 2d-g ), CR monkeys with higher IL-6 showed mitigated atrophy relative to controls within intraparietal areas along the inferior longitudinal fasciculus [-22, 2, 14; t(38) = 3.41; voxel cluster size = 1,636] and a local maxima in medial longitudinal fasciculus [-17, 4, 22; t(38)= 3.07]; cingulum, uncinate fasciculus, and other sub-cortical bundles along the frontal forceps [2, 38, 16; peak t-value = 3.36; voxel cluster size = 2,525], as well as a local maxima in uncinate fasciculus in deep WM [-8, 37, 14; peak t-value = 3.04]; and extreme capsule and striatal bundle proximal to orbital PFC [14, 32, 14; peak t-value = 3.34; voxel cluster size = 1,108]. Progressively higher IL-6 across all animals was related to lower FA values along the genu, body, and splenium of the corpus callosum, as well as posterior cingulum ( Figure 3a ; Table 2 ). Clusters were also found within thalamic bundles at the levels of the ventral-posterior and pulvinar nuclei, as well as WM proximal to the somatosensory areas consistent with the superior longitudinal fasciculus ( Figure 3b-c ). Other associations were noted in arcuate fasciculus along an intraparietal extent of lateral sulcus, as well as anterior commissure and uncinate fasciculus. A small anterior portion of this latter tract showed overlap with the WM volume result map generated by the IL-6*CR contrast. To determine whether WM volume influenced the FA result, a peak voxel (2, 44, 16) in this area of overlap was taken. A full correlation indicated a mild-moderate negative relationship between FA and IL-6 [r(32) = -.392, p = .013], which was not influenced when partially out WM voxel values [r(29) = -.388, p = .016]. An IL-6*CR Status model was used to test where CR monkeys had higher FA relative to controls across subjects per pg/ml increase in IL-6, suggesting a protective effect. After correction, the CR monkeys showed a mitigated relationship with IL-6 associated atrophy in superior longitudinal fasciculus along the dorsal convexity anterior to the central sulcus ( Figure 3d-e ; Table 2 ). Removal of individual IL-6 values at the high or low end of the range did not undermine the associations or group differences beyond an expected loss of statistical power. Split level partial correlations indicated that the positive association of IL-6 and FA for CR monkeys was non-significant (data not shown). To determine whether FA results were driven primarily by myelin or axonal changes, post hoc analyses of radial and axial diffusivity maps were conducted within the space of the FA results map. Increased peripheral IL-6 was related to greater radial diffusion along the corpus callosum, thalamic bundles, and other areas highly consistent in cluster size and significance with the FA map ( Supplementary Figure 1 ). By contrast, higher axial diffusivity was restricted to part of the splenium and corticocortical fibers adjacent to the somatosensory areas. Comparable to FA, higher levels of IL-6 at scan were associated with increased MD of WM throughout an extensive cluster encompassing the corpus callosum and cingulum ( Figure 4 ; Table 3 ). MD in GM tissue was also higher in several brain areas. The main cluster ranged from the global maxima adjoining right hemisphere insula and somatosensory areas out to local maxima caudally in the splenium and rostrally to ventral PFC. The caudal extent in part included occipito-splenial fibers; intraparietal and parietal cortex areas; temporal regions generally along the superior temporal sulcus; posterior cingulate cortex; and longitudinal fasciculi in parietal lobe. Association regions extended anteriorly to dorsal premotor cortex, supplementary motor area (SMA), pre-SMA, frontal eyefield, and then dorsal PFC. This rostral portion of the cluster consisted of anterior cingulate cortex; dorsal, infralimbic, orbitofrontal and ventral PFC, as well as gyrus rectus; precentral operculum; optic tract; agranular insula; and the basal forebrain. Medial and lateral extents incorporated striatum; dysgranular anterior insula and claustrum; uncinate, arcuate, and longitudinal fasciculi; superior temporal sulcus in medial temporal lobe; and medial temporal lobe structures including parahippocampus, hippocampus, and amygdala. Another separate cluster included occipito-splenial fibers and extrastriate cortex. An IL-6*CR Status model tested where CR monkeys had lower MD relative to controls for every pg/ml increase in IL-6 across all subjects. A striking cluster extended approximately 25 mm along the frontal convexity in WM and GM, spanning from dorsal premotor cortex to dorsal PFC and additionally incorporating superior longitudinal fasciculus, SMA, pre-SMA, frontal eye field, and anterior cingulate cortex ( Figure 4 ; Table 3 ). This region bilaterally overlapped with part of the IL-6 main effect. Clusters were also found in primary and extrastriate visual cortices, as well as hippocampus and uncinate fasciculus near the superior temporal sulcus. Removal of individual IL-6 values at the high or low end of the range did not undermine the associations or group differences beyond an expected loss of statistical power. Further, a similar result map with comparable maxima and sub-maxima was produced when the average of IL-6 taken at 20 years of age and the time of scan was used as the independent variable [-6, 33, 34; t(29) = 4.52; voxel cluster size = 5,369]. This averaged IL-6 was normally distributed and had no especially high or low values.

Higher peripheral levels of IL-6 in these aged monkeys were associated with smaller regional brain volume and less microstructural density primarily in parietal, frontal, callosal, and prefrontal areas. These age-sensitive regions have been implicated in voxel- or ROI based analyses in aged humans, rats, and monkeys ( Hauss-Wegrzyniak et al., 2000 ; Baron et al., 2001 ; Rabinovici et al., 2007; Alexander et al., 2008 ; Falangola et al., 2008 ). Higher IL-6 was also correlated with moderately less global GM and WM. Similar relationships have been found when correlating age or IL-6 to TBV and global GM or WM among normally aging humans ( Jernigan et al., 2001 ; Liu et al., 2003 ; Jefferson et al., 2007 ). Importantly, the monkeys on the CR diet had significantly lower IL-6 levels at the time of the scan and would manifest less neural atrophy specifically related to IL-6 versus controls. Due to the main effect of CR on GM that has been found previously ( Colman et al., 2009 ), dietary condition was included as a covariate to isolate the association between IL-6 and the current neural indices. By examining CR*IL-6 interactions, it was further possible to show that the relationship between IL-6 and atrophy in several neural regions was attenuated among CR monkeys. There are several mechanisms through which CR may confer protection against IL-6 associated atrophy beyond the main effect of lowering IL-6 itself. In mouse epididymal adipose tissue, CR downregulates expression of 55 proinflammatory genes related to the complement system, acute phase response, interferons, and lysosomal proteases ( Higami et al., 2006 ). Attenuated synthesis of any of these proteins may inhibit the ability of IL-6 to adversely affect cellular function and subsequently moderate neural atrophy. CR also raises anti-inflammatory factors in circulation such as Interleukin-10, which is related to decreased IL-6 synthesis ( Jung et al., 2008 ) and may centrally protect parenchymal tissue sensitive to the influence of PIC. The lack of IL-6 association with lower regional versus global WM volume was unexpected, but may be explained by DTI. Aging increases redundant myelin lamellae and oligodendrocytes in the brain ( Peters et al., 1991 ; Wozniak & Lim, 2006 ), which could falsely indicate intact regional WM even as axonal or myelin loss occurs. Indeed, voxel-wise regression analyses of systemic IL-6 with reduced FA and increased MD were highly significant in regions and tracts that are the first to show myelin changes with age, such as the genu of the corpus callosum and PFC ( Hayakawa et al., 1991 ; Pfefferbaum et al., 2005 ; Pagani et al., 2008 ). Indeed, increased radial diffusivity of the corpus callosum and other regions accurately corresponded with decreased FA, suggesting that increased systemic inflammation represented by peripheral IL-6 was primarily related to demyelination ( Song et al., 2002 ). WM volume did not appear to influence these results when included as a covariate. Thus, IL-6 showed cross-modal correspondence of volumetric and microstructural association clusters in the parietal and temporal lobes, whereas prefrontal and callosal correlations were significant only for DTI indices. Peripheral IL-6 concentrations among clinically healthy old monkeys were significantly related over time, which is comparable to results in human elders ( Rao et al., 1994 ). The IL-6 levels were not related to central adiposity or body mass, but did appear to be associated with chronic inflammatory pathologies in a subset of the animals in both dietary conditions. Although it is now well accepted that adipose tissue produces IL-6 ( Bulcão et al., 2006 ), and that IL-6 covaries with BMI in humans ( Marsland et al., 2008 ), the current sample of monkeys did not include morbidly obese animals in either group. It has been established that the CR monkeys weigh much less than do controls and exhibit many metabolic benefits, including better glucoregulation ( Colman et al., 2009 ). Yet, in concert with Marsland and coworkers (2008) , it is thus unlikely in this primate model that central adiposity is the primary mediator of the relationship between IL-6 and neural atrophy. By examining CR*IL-6 interactions, monkeys on CR diets had reduced IL-6 related volumetric atrophy and lower microstructural density in frontal and prefrontal areas. These same areas manifest a decreased multi-chemical MRI spectroscopy profile with human aging, including sensorimotor cortex and dorsal PFC ( Grachev and Apkarian, 2001 ). The most striking CR*IL-6 moderation of reduced FA and increased MD incorporated GM and WM tissue along the dorsal convexity, including superior longitudinal fasciculus, dorsal premotor cortex, SMA, pre-SMA, frontal eyefield, and dorsal PFC. This cluster shared substantial overlap with a comparable region found in the IL-6 association analysis for MD. This correspondence suggests that individuals on a CR diet, with different physiological profiles, may show either moderated IL-6 related atrophy or stronger abatement of its correlational relationship with neural regions that are critical for daily function. Integrity of pre-SMA tissue, for example, is important for performing voluntary saccades ( Boxer et al., 2006 ). The protective effect of CR, beyond the main effect reduction of IL-6 near scan, was also related to preserved WM and GM tissue density in hippocampus and uncinate fasciculus, which are critical structures that are relevant for normal and pathological aging ( Uylings et al., 2002 ). There are several limitations to the current cross-sectional study. It was not feasible to collect CSF or histological tissue from this longitudinally studied monkey cohort to establish causality between cytokine activity and induction of neural atrophy. Such samples will eventually be collected at time of death to permit a direct examination of the degree of neurodegeneration at mortality. Due to the correlational nature of the analyses, it was also not possible to determine whether systemic inflammation impacted the brain or, alternatively, the neural changes were occurring in parallel or perhaps even antecedent to the shift in inflammatory physiology. However, most studies have been based on the conceptual premise that peripheral PIC play a causative role in inducing subsequent upregulation in neural parenchyma ( Dantzer, 2004 ). To be sure, Offner and colleagues (2009) found that ischemic lesions lead to marked immunosuppression of circulating PIC like interferon gamma and tumor necrosis factor-alpha ( Offner et al., 2009 ). Finally, it must be emphasized that IL-6 is likely not the sole biomarker of inflammation related to atrophy in age-sensitive regions. The PIC tend to act as a unit by co-activating one another in the periphery and brain, consequently mediating neural and behavioral changes ( Dantzer, 2004 ). However, only IL-6 is reliably detected in healthy individuals, and it is a stable biomarker of low-grade inflammation in individuals across time ( Rao et al., 1994 ; Reyes and Coe, 1996 ; Willette et al., 2007 ). In summary, systemic IL-6 appears to provide a valuable and sensitive index of health and aging in the nonhuman primate, which is correlated with the magnitude of brain aging globally and in a number of sensitive regions.

Forty-five rhesus monkeys between 19 and 31 years of age were utilized from a longitudinal CR project at the Wisconsin National Primate Research Center (WNPRC). 18 animals were fed a normal diet (mean age = 23.84 ± 2.79 (SD); 12 females, 6 males), while 27 CR subjects had been on a moderately restricted diet (30% reduction of intake) for approximately 12 to 17 years (mean age = 24.32 ± 2.77 (SD); 15 females, 12 males). Details of the CR manipulation have been described previously ( Kemnitz et al., 1993 ; Ramsey et al., 2000 ). Animals were maintained at 21 °C; humidity 50-65%; 12h/12h light-dark cycle. Food was present for 6-8 hours per day; water was available ad libitum . Animals were individually housed and had continuous visual and auditory contact with other monkeys in the room. Research and care staff regularly inspected and interacted with the cohort approximately three times per day. Manipulanda and other environmental enrichment were provided to alleviate boredom. Scan acquisition, image processing, and IL-6 assays were conducted by personnel blind to dietary condition. The study was approved by the Institutional Animal Care and Use Committee. Blood was collected at two time points: 1) within 6 months of the neuroimaging scan date to determine basal levels of circulating IL-6, and 2) when each animal reached 20 years of age. IL-6 was assayed using a commercial Enyzme-Linked Immunosorbent Assay (ELISA; R&D Systems, Minneapolis, MN) using methods similar to those described previously ( Willette et al., 2007 ). Body weight (kg) was recorded within 48 hours prior to the MRI scan. Central adipose tissue in the peritoneum was assessed by MRI using the following procedure. First, a region of interest (ROI) was traced around the viscera of the thoracic cavity subjacent to subcutaneous adipose tissue. Subject-specific intensity thresholding was then used to isolate adipose tissue, followed by summing all voxels to obtain total central adiposity. It has been reported that Body Mass Index (BMI) does not appear to influence the relationship of IL-6 and neural atrophy ( Marsland et al., 2008 ). Thus, for the purposes of this report, only visceral adiposity was evaluated as a correlate of IL-6 values. Finally, the presence of veterinarian-diagnosed inflammatory or illness conditions near the time of scan was considered, because such pathologies might markedly raise PIC ( Brouqui et al. 1994 ; Devaux et al. 1997 ; De Keyser et al. 1998 ). The small sample size for glucoregulatory pathologies such as diabetes precluded statistical analysis (controls, n = 4; CR, n = 0). Inflammatory ailments included: arthritis (3), peritonitis (1), endometriosis (3), adenocarcinoma (4), gastric bloat (1), and bloat co-morbid with endometriosis (1). No acute or severe clinical illnesses due to infectious pathogens were noted near the time of sampling or scan. Monkeys were thus grouped into a binary variable as clinically healthy or having an inflammatory condition that could potentially affect IL-6 levels. Animals were sedated and images acquired as previously described on a General Electric 3.0 T Signa MR unit (GE Medical Systems, Milwaukee, WI, USA) using a quadrature Tx/Rx volume coil with an 18cm diameter ( McLaren et al., 2009 ). Briefly, monkeys were maintained immobile in an optimal plane of anesthesia. A three-dimensional head-only coronal T1-weighted inversion recovery-prepped spoiled gradient echo (IR-prepped SPGR) sequence was collected with the following parameters: repetition time (TR) = 8.772 ms; echo time (TE) = 1.876 ms; inversion time (TI) = 600 ms; flip angle = 10°; number of excitations (NEX) = 2; matrix = 256×224; field of view = 160mm. 124 coronal slices were acquired with a thickness of 0.7 mm and no gap resulting in near isotropic 0.625 × 0.625 × 0.7 mm voxels. Head-only DTI data were acquired using a single-shot, spin-echo, diffusion-weighted echo-planar imaging (DW-EPI) sequence with diffusion gradients in twelve optimal directions and was employed using minimum energy criterion ( Hasan et al., 2001 ). A diffusion-weighting of b = 816 s/mm 2 was used. Other imaging parameters for this sequence were: axial acquisition; TR/TE = 10000/77.2ms; flip angle = 90°; NEX = 6, field of view = 160 mm; in-plane matrix = 120×120; slice thickness = 2.5mm; slice gap = 2.5mm; 124 slices. Image distortion in the DW-EPI data was minimized by a higher-order shimming protocol that was run before the DTI scan. Voxel size was resampled to 0.5 × 0.5 × 0.5 mm for both scan sequences. Finally, for the torso, a three-dimensional coronal T1-weighted fast gradient echo (FGRE) sequence was acquired during the same scan session. Parameters were: TR/TE = 15.7/1.404 ms; 27 coronal sections; flip angle = 20°; slice thickness = 4 mm; slice gap = 6 mm; NEX = 2; matrix =256 × 128; field of view = 200 mm; 0.781 × 0.781 × 6 mm voxels. T1-weighted images were manually deskulled, rotated, normalized to the 112RM-SL monkey atlas, and segmented into GM and WM probability maps in SPM5 as previously described ( McLaren et al., 2009 , 2010 ). Subject data were normalized to a common space for VBM analyses through Diffeomorphic Anatomical Registration using Exponentiated Lie Algebra (DARTEL; Ashburner, 2007 ; McLaren et al., 2010 ). These warped and modulated GM and WM images were thus moved into the same space, coordinates, and orientation as the D99-SL brain ( Saleem et al., 2002 ). A 4mm FWHM Gaussian kernel was applied to smooth the images for analysis. An absolute threshold of 0.1 was applied for GM or WM volume to restrict analyses to the tissue class of interest. The DTI data were processed as follows: image distortions caused by eddy currents were corrected using FSL’s FDT Diffusion Toolbox ( http://www.fmrib.ox.ac.uk/fsl/fdt/index.htm ). Three-dimensional maps of the diffusion tensor and derived measures were calculated using DTIFIT in FSL, including: FA, MD, axial diffusivity (λ 1 ), and radial diffusivity ([λ 2 + λ 3 ] / 2]). In order to bring DTI maps into the 112RM-SL atlas for voxel-wise analysis, each subject’s non-diffusion weighted image (B0 map) was coregistered and normalized to a T2-weighted template through the combination of two transforms, where T2 acquisition has been described elsewhere ( McLaren et al., 2009 ). The normalized images were smoothed using a 4mm FWHM isotropic Gaussian kernel. FA analyses were restricted to WM by using a binary prior probability map, in order discern changes in major fiber tracts; an absolute threshold of 0.1 was applied for each monkey to reduce the influence of morphology on results. MD maps were masked using a binary version of the 112RM-SL rhesus atlas, because MD is not overtly biased toward GM or WM tissue classes relative to FA. MD results did not change when binary GM and WM prior probability maps were separately used as masks and the resulting MD result maps combined. Images were inspected by experts (ASF, GX) in order to locate scanner artifacts that could influence volumetric or microstructural results. One T1-weighted scan and nine DTI scans were excluded from analyses due to motion, acquisition, respiration, and other distortion artifacts. Due to wide variation in IL-6 levels for each group, a two-tailed Mann-Whitney ranks test using the Z-statistic was used to assess group differences in peripheral IL-6 at time of scan, which is reported as the median and interquartile range (IQR). This nonparametric test is less sensitive to the influence of very high or low values than a comparison of means. Two-tailed Pearson’s r correlations were used to evaluate if IL-6 levels were related to body weight, visceral adiposity, and health status, as well as to TBV, global GM and WM. IL-6 levels at 20 years of age and at time of scan were also compared. Alpha was set at .05 for these analyses. For neural association analyses, covariates included age, gender, and dietary condition; TBV was additionally incorporated when testing for regional volume associations. GM or WM probability maps were entered into regional analyses of volume as the dependent variable ( Ashburner and Friston, 2000 ). The independent variable for association analyses was IL-6 collected near the time of scan--although the average of IL-6 at both time points was used to confirm the stability of result maps. Of particular interest was if CR monkeys showed incrementally less atrophy per pg/ml unit increase in IL-6 relative to controls, suggesting a protective effect. Thus, an IL-6*CR interaction term was entered into separate general linear models to test for this relationship; IL-6 and dietary condition vectors, in addition to standard covariates, were included to remove their main effects. The t-statistic threshold at the voxel level was set at p < 0.005 (uncorrected). Analyses were corrected for multiple comparisons at the cluster level using Monte Carlo simulations, in the Analysis of Functional NeuroImages (AFNI) package AlphaSim, to derive a cluster corrected p-value less than .05 (Forman et al., 1995). This correction required GM volume clusters to contain at least 875 edge-connected voxels in order to be significant. For WM volume areas, as well as FA maps, the necessary cluster size after correction was 854 voxels. MD analyses examining both GM and WM required a cluster size of at least 1,106 voxels. Although violation of non-stationarity is a concern for cluster correction of structural imaging data, the use of a 4mm FWHM kernel does not appreciably bias the true alpha when using this permutational estimate of cluster size ( Hayasaka et al., 2004 ). To determine if alterations in myelin or axons predominantly drove FA results ( Song et al., 2002 ), the IL-6 association FA map was used as an explicit mask for post hoc voxel-wise analyses of radial and axial diffusivity thresholded at a voxel alpha of .05 (uncorrected) over 50 contiguous voxels. Coordinates of significant clusters correspond to the space of the Saleem-Logothetis atlas and are displayed on the 112RM-SL underlays ( Saleem et al., 2002 ). Major fiber tracts were designated based on autoradiography conducted in rhesus monkeys ( Schmahmann and Pandya, 2006 ).

Post hoc association of IL-6 and radial and axial diffusivity in areas of significant FA associations. The FA map was used as an explicit mask for these voxel-wise analyses. The radial diffusion and FA maps closely corresponded in the location and size of clusters, whereas axial diffusivity is related to FA in part of the splenium and somatosensory corticocortical fibers. The ‘warm’ and green color schemes for FA and post hoc analyses represent t-values. Brains are oriented in neurological space.