Anatomical and diffusion MRI of deep gray matter in pediatric spina bifida.

Individuals with spina bifida myelomeningocele (SBM) exhibit brain abnormalities in cortical thickness, white matter integrity, and cerebellar structure. Little is known about deep gray matter macro- and microstructure in this population. The current study utilized volumetric and diffusion-weighted MRI techniques to examine gray matter volume and microstructure in several subcortical structures: basal ganglia nuclei, thalamus, hippocampus, and amygdala. Sixty-six children and adolescents (ages 8-18; M = 12.0, SD = 2.73) with SBM and typically developing (TD) controls underwent T1- and diffusion-weighted neuroimaging. Microstructural results indicated that hippocampal volume was disproportionately reduced, whereas the putamen volume was enlarged in the group with SBM. Microstructural analyses indicated increased mean diffusivity (MD) and fractional anisotropy (FA) in the gray matter of most examined structures (i.e., thalamus, caudate, hippocampus), with the putamen exhibiting a unique pattern of decreased MD and increased FA. These results provide further support that SBM differentially disrupts brain regions whereby some structures are volumetrically normal whereas others are reduced or enlarged. In the hippocampus, volumetric reduction coupled with increased MD may imply reduced cellular density and aberrant organization. Alternatively, the enlarged volume and significantly reduced MD in the putamen suggest increased density.

Introduction

Spina bifida myelomeningocele (SBM), identified by the herniated protrusion of the spinal cord and meninges from the vertebrae at birth (Detrait et al., 2005), is the most common congenital birth defect affecting the central nervous system compatible with survival. It occurs in approximately .3–.5 per 1000 live births (Au et al., 2010). Almost all individuals with SBM incur structural brain abnormalities, including variable anomalies of the midbrain and corpus callosum. Most consistently observed is the Chiari II malformation of the hindbrain, which occurs in over 90% of cases with SBM (Juranek and Salman, 2010). Resultant obstruction at the level of the fourth ventricle leads to hydrocephalus in the majority (roughly 80%) of people with SBM, with many receiving shunt diversion as treatment.

The high comorbidity of hydrocephalus and shunt treatment in individuals with SBM disrupts brain macro- and microstructure (Del Bigio, 2010), leading to variable neural and cognitive outcomes (Hampton et al., 2011). Although intellectual disability is infrequent in SBM, the neurobehavioral profile of SBM consists of strengths in lower-order cognitive and motor domains, such as association and motor learning. Weaknesses are in higher-order, contextually based cognitive processes, including problem solving and on-line information integration (Dennis and Barnes, 2010). There are often relative strengths in motor learning and adaptation to environmental stimuli, but weaknesses in dynamic motor regulation such as impaired movement coordination and fine motor control (reviewed in: Dennis and Barnes, 2010).

Given the characteristic features of the Chiari II malformation and motor impairments, initial neuroimaging studies of SBM targeted cerebellar abnormalities, revealing generally reduced cerebellar volumes in these individuals relative to typically developing (TD) peers (Dennis et al., 2004; Juranek and Salman, 2010). However, despite overall volume reductions, individuals with SBM show significant volumetric enlargement in anterior cerebellar portions and reductions in posterior–inferior regions (Juranek et al., 2008).

At the level of the cerebrum, individuals with SBM also exhibit anterior to posterior patterning in the cortex as well as in the corpus callosum (Hannay et al., 2009; Hannay et al., 2008; Juranek et al., 2008; Treble et al., 2013). Cortical thickness is significantly increased in frontal regions, but decreased in posterior regions (Juranek et al., 2008). When examining cortical complexity (i.e., gyrification) a similar pattern emerges whereby frontal regions generally have significantly less complexity (except the middle frontal region) despite significantly greater complexity in parietal and temporal cortices relative to TD controls (Treble et al., 2013). Posterior regions of the corpus callosum may also be more affected in SBM; the genu is generally well preserved whereas the splenium shows the greatest dysgenesis (Hannay et al., 2009).

White matter (WM) and gray matter (GM) volumetric and microstructural alterations also exist in individuals with SBM. WM integrity is particularly reduced in long association fiber tracts connecting posterior and anterior brain regions (Hasan et al., 2008a; Ou et al., 2011). Most studies have reported decreased fractional anisotropy (FA) and increased mean diffusivity (MD) in SBM relative to controls (Ou et al., 2011; Hasan et al., 2008a; Williams et al., 2013). Lower FA has also been observed in the genu of the corpus callosum, but not in the anterior commissure (Herweh et al., 2009). In a recent probabilistic tractography study of attention pathways in SBM, posterior pathways showed greater reduction in WM integrity (decreased FA) than frontal pathways (Williams et al., 2013). Tectal beaking occurs with high frequency in people with SBM and was further associated with WM alterations in both anterior and posterior attention pathways. Limbic tract integrity, as previously examined through fiber tracking and region of interest (ROI) techniques, is also reduced in SBM, particularly in the fornix and cingulum (Kumar et al., 2010; Ou et al., 2011; Vachha et al., 2006).

Little is known about deep GM structure and integrity in SBM. This knowledge is relevant given disruptions in periventricular white matter, the fornix, and long association fibers, which have connections with subcortical GM structures such as basal ganglia nuclei and the hippocampus. However, no studies to our knowledge have examined specific limbic system structures and few have investigated the thalamus and basal ganglia structures. Of the three studies that have investigated specific basal ganglia nuclei and periventricular GM regions using quantitative methods, the most consistent finding is reduced GM integrity in the caudate nucleus of the basal ganglia (Hasan et al., 2008b; Kumar et al., 2010; Kumar et al., 2011). However, discrepancies exist between these three studies of GM integrity. While decreased FA and increased MD values were observed in the caudate nucleus of children with SBM relative to controls in one study (Kumar et al., 2010), two other studies reported increased FA in this region relative controls (Hasan et al., 2008b; Kumar et al., 2011).

Similarly, some discrepancies exist with regard to published reports of putamen integrity. One study reported no significant differences in SBM compared to controls in DTI metrics of the putamen (Hasan et al., 2008b) and another study reported increased FA in this structure (Kumar et al., 2011). Kumar et al. (2011) also investigated GM integrity of the thalamus, but reported no significant differences between patients and controls (Kumar et al., 2011), despite its proximal location to periventricular WM, which has abnormal features in SBM (Hasan et al., 2008a). To date, no studies have examined whether volumetric differences in these regions (i.e., limbic, basal ganglia, thalamus) are evident in SBM. The clinical implications of elucidating subcortical GM volume and integrity differences are relevant since most GM structures are clinically read as “normal” during qualitative radiological review (Fletcher et al., 2005).

To further understand subcortical GM in SBM and clarify current discrepancies in existing literature, we examined the volumetric and microstructural properties of deep GM regions in SBM and a typically developing comparison group.

Hypotheses

Volumetric measures (corrected for total brain volume) of deep GM structures will be systematically reduced in participants with SBM relative to TD controls.

Measures of microstructure in deep GM will reflect reduced (indicated by increased mean diffusivity and reduced fractional anisotropy) integrity of these regions in participants with SBM relative to typically developing controls.

Materials and methods

Participants

Participants included children and adolescents (N = 66) between 8 and 18 years (M = 12.0, SD = 2.73) of age with SBM (n = 48) and typically developing controls (TD; n = 18). Participants were drawn from the Houston cohort of an international project examining neurobehavioral outcomes of spina bifida (Fletcher et al., 2005). While participants with SBM were recruited from clinics at Texas Children's Hospital and Shriners Hospital for Children, TD participants were recruited from advertisements and community contacts and had no previous or current neurological difficulties. Participants with SBM had hydrocephalus (20.8% mild at time of testing) and all but one had been treated with a shunt near the time of birth. Most participants with SBM had lumbar and sacral spinal lesions (83.3%), Chiari II malformation (87.5%), callosal hypogenesis or hypoplasia (thinning) (95.8%), less than 5 shunt revisions (81.2%), and infrequent seizure history (none: 66.7%; past: 8.3%; current: 4.2%).

Participants in both groups had verbal and/or nonverbal intelligence quotient (IQ) scores between 70 and 120 per the Stanford–Binet Intelligence Scale: Fourth Edition (Thorndike et al., 1986). Full-scale IQ estimates ranged between 63 and 112 for the group with SBM and between 84 and 120 for the TD group. Some children with lower IQ scores were eliminated from the parent study because of concerns about their capacity for following directions on neuropsychological assessment procedures and are not discussed below. All study procedures, including MRI acquisitions, were conducted between 2005 and 2010 and were in compliance with IRB approval at The University of Houston and The University of Texas Health Science Center at Houston.

Demographic characteristics of the sample are displayed for each group in Table 1. As expected, participants with SBM had significantly lower verbal and nonverbal IQ scores (ps ≤ .001) than controls and significantly lower socioeconomic status (SES) per the Four Factor Index of Social Status (Hollingshead, 1975). The group with SBM also had a greater number of non-right handed individuals, which is common in samples of SBM (Fletcher et al., 2005) and other samples involving congenital brain injury (e.g., Ross et al., 1987).

Table 1

Demographic data for the spina bifida myelomeningocele (SBM) and typically developing control (TD) groups.

Variable	SBM(n = 42)	TD(n = 18)	p-Value	Pairwise
Sex (n [% male])	24 (57.1)	11 (61.1)	.421
Handedness (n [% right])	35 (72.9)	16 (88.9)	.051	TD > SBM
Age at MRI (M [SD])	12.34 (2.9)	11.24 (2.2)	.148
Ethnicity (n [% Hispanic])	28 (58.3)	14 (77.8)	.282
Socioeconomic status (M [SD])*	31.36 (13.2)	39.83 (11.3)	.018	TD > SBM
Full-scale IQ (M [SD])	85.23 (12.4)	104.18 (10.3)	<.001	TD > SBM
Verbal IQ (M [SD])	85.54 (16.0)	99.56 (13.0)	.001	TD > SBM
Nonverbal IQ (M [SD])	92.00 (14.1)	107.44 (14.1)	<.001	TD > SBM

Note: IQ and socioeconomic status estimates obtained from the Stanford–Binet Intelligence Scales and Four Factor Index of Social Status, respectively (Thorndike et al., 1986; Hollingshead, 1975).

Magnetic resonance image acquisition

All MRI acquisitions were acquired at the University of Texas Medical School in Houston using a research-dedicated Philips 3 T scanner with SENSE (Sensitivity Encoding) technology and an 8-channel phased array head coil.

T1-weighted imaging

After a conventional scout sequence, high-resolution T1-weighted anatomical images were acquired in the coronal plane using a 3D turbo fast echo sequence with the following parameters: voxel dimensions = .94 × .94, slice thickness = 1.5 mm, TR = 6.50–6.70 ms, TE = 3.04–3.14 ms, flip angle = 8°, DFOV = 240 mm, matrix = 256 × 256.

Diffusion tensor imaging (DTI)

DTI images were acquired in the axial plane using a single-shot spin-echo diffusion sensitized echo-planar imaging sequence. Diffusion sensitizing gradients were applied in 21 directions (weighting: b = 1000 s/mm2) with one reference image (b = 0 s/mm2) and the following parameters: voxel dimensions = .94 × .94 mm, slice thickness = 3 mm, TR = 6500 ms, TE = 65 ms, flip angle = 90°, DFOV = 240 mm, matrix = 256 × 256.

Magnetic resonance data analysis

All scans underwent initial radiological evaluation. The basal ganglia structures (i.e., caudate, putamen, globus pallidus) were all read as “normal” for both groups and the thalamus was read as normal for all participants in the TD group and for 95.8% of participants in the group with SBM (2 showed unilateral abnormality and 1 had massa intermedia < 1 cm).

Prior to data preprocessing, anatomical MRI and DTI data were reviewed for image quality by an expert user with extensive neuroimaging knowledge and experience (JJ) who was blind to age, gender, and neuropsychological performance outcomes. Morphometric analyses were conducted in FreeSurfer v4.0.5 () (Fischl, 2012) using a 64-bit Linux computer. Skull-stripping and regional segmentation were executed using a fully automated command (e.g., recon-all –all) for each T1-weighted image, yielding 3 classes of voxels: gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) (Fischl et al., 2002; Fischl et al., 2004a; Fischl et al., 2004b). Subsequently, within FreeSurfer, a fully automated routine determined delimiting boundaries of deep GM structures (i.e., caudate, putamen, pallidum, thalamus, hippocampus, amygdala). The brains of people with SBM are dysmorphic and boundaries are not always reliably indicated by the FreeSurfer output. The output images from FreeSurfer's segmentation and classification algorithms were visually inspected and extensively manually edited by a single rater (JJ) with substantial expertise using FreeSurfer's software as well its predecessor, Cardviews (Filipek et al., 1989). An example of finalized segmentation after manual editing is shown in Fig. 1. At the time of MRI acquisitions (collected between 2005 and 2010), FreeSurfer v4.0.5 (released June 2008) was the most recent version available for analyses. Although subsequent revisions of FreeSurfer were released, FreeSurfer v5.0.0 (released in 2010) was the first version that retained manual edits (i.e., persistent edits) and the older version of FreeSurfer was maintained to ensure consistency across study subjects.

Fig. 1

Examples of final segmentations following extensive manual edits. A) SBM participant's T1-weighted MRI (left panel) with color-coded segmentation boundaries of deep gray matter structures (right panel). B) TD participant's T1-weighted MRI (left panel) with color-coded segmentation boundaries of deep gray matter structures (right panel).

For DTI data, within-participant co-registration of the GM masks (obtained from the FreeSurfer analyses) in T1-space with DTI space (B0 or nodif image) were performed using FSL v4.1.0 (FMRIB's Software Library, ) (Smith et al., 2004). Minor head motion and eddy currents were corrected in the DTI series with the Eddy Current Correction tool included in FSL's Diffusion Toolbox v2.0. Skull-stripping and removal of non-brain tissue were performed using BET v2.1. FSL's Linear Image Registration Tool v5.5 (FLIRT) (Jenkinson et al., 2002) was used to perform a within-participant inter-modal linear registration between each participant's T1-weighted series (reference image) and corresponding non-diffusion-weighted images (b = 0) from the DTI dataset. The resultant transformation matrix (and calculated inverse) was used to transform the binary segmentation masks of each GM structure (obtained from high-resolution T1-weighted images using FreeSurfer) to corresponding diffusion-weighted space (detailed in Juranek et al., 2012). The diffusion tensors were reconstructed using FSL's DTIFIT tool within the Diffusion Toolbox. This transformation matrix was then applied to each participant to mask individual maps of MD and FA for each GM structure. To be as conservative as possible, each participant's binary segmentation mask was eroded using a 2 × 2 × 2 mm3 kernel to reduce contamination from neighboring voxels of CSF and/or WM. Final quantitative data were obtained using command line utilities to compute average MD and FA circumscribed by each eroded deep GM mask in each participant.

Data analyses

Preliminary analyses included examination of the distributions of the primary volumetric and DTI variables within and across groups and evaluation of hemispheric differences by group in addition to potential covariates and moderators (e.g., age, sex). IQ was not considered as an appropriate covariate. A covariate that is an attribute or an intrinsic characteristic of the condition cannot be meaningfully adjusted to “control” for the attribute, instead creating groups with artificial mean differences reflecting the association of the covariate with the disorder (Dennis et al., 2009). In addition, handedness was not included as a covariate given the lack of variance in the TD group (Kutner et al., 2005).

Across groups, distributions of volumetric data were generally normal. Within groups, large deviations from group means (˜2 SD) were rare, but where they occurred, they were specifically evaluated for impact. The group with SBM had significantly smaller total brain volume (TBV) relative to the TD group, F(1,64) = 24.40, p < .001. Therefore volumetric measures were evaluated using corrected ratios (GM structure volume to TBV) for all analyses.

To examine potential hemispheric differences between groups, separate 2 (hemisphere: left, right) × 2 (group: SBM, TD) repeated measures ANOVAs were examined for each corrected GM volume and DTI value (MD, FA). With the exception of thalamus volumes (p = .031), none of the group × hemisphere interactions were significant for GM volumes or DTI metrics (p > .05). Across groups, there were significant hemisphere effects (p < .015) found for volumes of pallidum and putamen, FA in all basal ganglia structures (caudate, pallidum, putamen), and MD in the amygdala, pallidum, and putamen. However, since none of the group × hemisphere interaction effects were significant, indicating similar lateralization of volumetric and DTI metric values for participants in the SBM and TD groups, hemisphere was not included in the main analyses below. There were also no significant group × hemisphere × sex interactions, indicating hemispheric similarity for males and females. Given these findings, whole corrected GM structure volumes (average of left and right hemispheres) and DTI values were computed for each GM structure and examined in subsequent analyses, with the exception of thalamus volumes.

Across groups, age was not significantly correlated (r < .190) with GM volumes, although when groups were examined separately age was associated with putamen volumes in the TD, but not in the group with SBM. However, age was significantly associated with microstructure values both across and within groups separately. Since the literature has shown age effects on GM volume and microstructure values across development (e.g., Lebel et al., 2008), it was included as a covariate in all analyses.

To address the study hypotheses, two sets of results were examined, first comparing deep GM volumes and then DTI values between the SBM and TD groups on each structure. Each deep GM structure was examined using separate ANCOVAs with group (SBM, TD) entered as the independent variable. For thalamus volume, for which there was a significant group × hemisphere interaction, a 2-way repeated measure ANCOVA with one repeated factor (hemisphere) and one between-subjects factor (group) was utilized. For all models, age was included.

Results

Macrostructure

Table 2 presents means and standard deviations for each deep GM structure volume (with and without correction for TBV, which differed between groups). ANCOVA results (with covariate age) indicated significant effects (shown in Fig. 2) of group for volumes of the putamen, F(1,63) = 8.55, p = .005, η2partial = .119, and hippocampus, F(1,63) = 39.73, p < .001, η2partial = .387, but not for the pallidum, amygdala, or caudate, F(1,63) < 1, p ≥ .335. Participants with SBM had smaller hippocampi and larger putamen volumes. Age was only significantly positively associated with putamen volume, F(1,63) = 5.05, p = .028, η2partial = .074.

Table 2

Mean (standard deviations) for raw and corrected (raw to total brain volume) deep gray matter structures for participants with spina bifida myelomeningocele (SBM; n = 48) and typically developing controls (TD; n = 18). Values reflect averages between left and right hemispheres.

	Raw volume (mm3)		Corrected volume (mm3)
Brain region	SBM	TD	SBM	TD	Pairwise(corrected)
Total brain volume	1,183,298.86 (97,723.81)	1,291,996.14 (83,270.50)			SBM < TD ***
Hippocampus	4845.14 (1060.48)	7034.00 (812.60)	.004 (.0009)	.005 (.0007)	SBM < TD ***
Amygdala	3699.86 (684.60)	3961.27 (419.55)	.003 (.0006)	.003 (.0003)
Caudate	7434.11 (1147.82)	7912.05 (812.64)	.006 (.0008)	.006 (.0006)
Putamen	12,302.23 (1442.31)	12,731.45 (1007.44)	.0104 (.0010)	.010 (.0006)	SBM > TD **
Pallidum	3682.62 (594.91)	3944.68 (586.11)	.003 (.0005)	.003 (.0004)	SBM > TB *+
Thalamus++	14,910.35 (2430.92)	15,743.18 (1856.24)	.013 (.0021)	.012 (.0012)

Note: reported values reflect the average volume for left and right hemispheres.

Group differences for pallidum volume were no longer significant after demographic variables (sex, socioeconomic status) were included as covariates.

Between-group differences on thalamus volumes were examined using a 2 (group) × 2 (hemisphere) within-participants ANCOVA.

p ≤ .05.

p ≤ .01.

p ≤ .001.

Fig. 2

Group differences on corrected putamen (p = .005) and hippocampus (p < .001) volumes between participants with spina bifida myelomeningocele (SBM; n = 48) and typically developing controls (TD; n = 18). Values reflect averages between left and right hemispheres. Note: bars represent standard error.

For the thalamus, the effects of hemisphere and hemisphere × age were not significant (ps ≥ .593), although the hemisphere × group effect was significant, F(1,63) = 4.25, p = .043, η2partial = .063. For the between-group effects, neither age nor group was statistically significant (ps ≥ .274). Follow-up examination of the hemisphere × group interaction revealed a greater, though nonsignificant (ps ≥ .139), discrepancy between groups (SBM > TD) in the left (M[SD]: SBM = .0063 [.0011]; TD = .0059 [.0005]) relative to the right (M [SD]; SBM = .0061 [.0009]; TD = .0060 [.0005]) thalamus.

Microstructure

Results for microstructure (MD and FA) are presented in Table 3. Only overall models and effects that met critical level of alpha (p < .05) are discussed below.

Table 3

Means (standard deviations) of DTI values for the spina bifida myelomeningocele (SBM, n = 49) and typically developing control (TD, n = 18) groups.

	Mean diffusivity			Fractional anisotropy
Brain region	SBM	TD	Pairwise	SBM	TD	Pairwise
Hippocampus	1.089 (.127)	0.954 (.056)	SBM > TD ***	0.178 (.044)	0.168 (.048)
Amygdala	0.894 (.014)	0.855 (.139)		0.209 (.030)	0.186 (.016)	SBM > TD **
Caudate	0.789 (.049)	0.751 (.015)	SBM > TD ***	0.163 (.030)	0.125 (.016)	SBM > TD ***
Putamen	0.745 (.015)	0.756 (.011)	SBM < TD *	0.158 (.017)	0.143 (.012)	SBM > TD ***
Pallidum	0.765 (.029)	0.774 (.019)		0.301 (.033)	0.244 (.027)	SBM > TD ***
Thalamus	0.799 (.035)	0.775 (.022)	SBM > TD **	0.302 (.020)	0.275 (.014)	SBM > TD ***

Note: MD M (SD) reported with values ×10−3; reported values reflect the average MD/FA for left and right hemispheres.

p < .05.

p ≤ .01.

p ≤ .001.

Mean diffusivity

There was a significant effect of group on MD values in the thalamus, F(1,63) = 9.26, p = .003, η2partial = .128, caudate, F(1,63) = 12.57, p = .001, η2partial = .166, putamen, F(1,63) = 5.12, p = .027, η2partial = .075, and hippocampus, F(1,63) = 19.67, p < .001, η2partial = .238. MD was significantly higher for the group with SBM in the thalamus, caudate and hippocampus and significantly lower in the putamen relative to the TD group. Age was significantly negatively associated with MD in the putamen, F(1,63) = 6.08, p = .016, η2partial = .088, and pallidum, F(1,63) = 21.19, p < .001, η2partial = .252.

Fractional anisotropy

There was a significant effect of group on FA values in the thalamus, F(1,63) = 24.34, p < .001, η2partial = .279, putamen, F(1,63) = 11.22, p = .001, η2partial = .151, caudate, F(1,63) = 25.34, p < .001, η2partial = .287, pallidum, F(1,63) = 39.46, p < .001, η2partial = .385, and amygdala, F(1,63) = 10.92, p = .002, η2partial = .148. FA was significantly higher for the group with SBM in all structures relative to the TD group. Age was significantly positively associated with FA in the pallidum, F(1,63) = 4.57, p = .036, η2partial = .068, and in the thalamus, F(1,63) = 15.25, p < .001, η2partial = .195.

Sex and socioeconomic status as covariates

Sex

Given the potential sex differences in GM volumes and microstructural metric values, sex was added to the original model as a covariate along with age. With covariates sex and age, volumetric results remained similar except for reduced group differences in pallidum volumes (p = .551); groups remained significantly different on hippocampus (p < .001) and putamen (p = .004) volumes. Group differences on MD values did not meet the critical value of alpha for the pallidum and amygdala, but remained significant in the putamen (p = .047), caudate (p = .001), thalamus (p = .004), and hippocampus (p < .001). The addition of sex as a model covariate also did not change results for FA values in any region; all deep GM structures, except the thalamus, had increased FA in the group with SBM relative to the TD group (p ≤ .001).

Socioeconomic status

Results were also not substantially altered by the addition of SES as a model covariate. When both age and SES were controlled, group differences in GM region volume remained similar, with the exception of pallidum volume (p = .717). The hippocampus remained significantly reduced and the putamen remained significantly enlarged in the group with SBM compared to controls. The additional covariate also did not substantially alter any of the findings for MD or FA; significant group differences in MD persisted in all regions except the amygdala and pallidum and groups remained significantly different on FA values in all structures except the hippocampus.

Examination of shunt status in the group with SBM.

Given that all but one participant in the group with SBM had received shunt intervention as treatment for hydrocephalus, the number of shunt revisions was examined as a predictor for GM volumes and MD and FA. Linear regression was used to predict whether more revisions resulted in smaller corrected volumes and worse integrity compared to fewer shunt revisions for all examined GM regions. For participants with SBM, regression results indicated that number of shunt revisions was not strongly associated with GM volumes (p ≥ .264), FA values (p ≥ .148), or MD values (p ≥ .080).

Discussion

Failure of neural tube closure in utero results in an altered pattern of programmed fetal development leading to gross anatomical brain dysmorphologies. In SBM, the considerable inter-individual variability in measures of both structure and function is documented (Fletcher et al., 2005). The present study probed distinct structural properties (both macro- and microstructural indices) of deep GM in children and adolescents with SBM relative to a TD comparison group.

Contrary to initial expectations of ubiquitous disproportionate reduction in volume and microstructural integrity of subcortical GM structures following SBM, findings (summarized in Tables 2 and 3) supported a model of regionally specific GM areas being differentially impacted by this early gestational insult. Whereas the hippocampus showed volumetric reduction, the putamen was significantly enlarged. The majority of other GM structures (i.e., thalamus, amygdala, caudate, globus pallidus) were volumetrically scaled to TBV in both groups. Similarly, results from the DTI metrics (i.e., MD, FA) suggested that deep GM microstructural integrity was not uniformly disrupted in SBM. Significant group differences across structures typically followed a pattern where the group with SBM showed increased MD values. Three exceptions were in the putamen, where MD was decreased, and in the amygdala and pallidum, which had similar values between groups. With respect to FA values, all evaluated structures (except the hippocampus) demonstrated significant group differences with increased FA values in the group with SBM relative to controls.

These findings are consistent with the overall concept that SBM results in a complex pattern of disproportionate enlargement in some regions (e.g., volume of the putamen and anterior cerebellum and prefrontal cortical thickness and complexity), disproportionate reduction of volume in other regions (e.g., the hippocampus and inferior–posterior cerebellum and parietal cortical thickness), in conjunction with typical structure of some brain regions (Juranek and Salman, 2010). The anomalous structural patterning in SBM may reflect complex compensatory mechanisms of neural plasticity in this population. Additionally, the mechanical effects of hydrocephalus likely complicate the protracted brain development in SBM.

The hippocampus may be particularly affected in SBM given the disproportionately reduced volume and increased MD observed relative to controls. Volumetric reduction coupled with increased MD could imply lower cellular density and/or aberrant cellular structure and organization in SBM. The hippocampus develops early in gestation (with individual subfields distinguishable around 15–19 weeks) and is particularly sensitive to early gestational factors associated with the occurrence of SBM, including malnutrition (Arnold and Trojanowski, 1996; Morgane et al., 2002) and mechanical injury due to hydrocephalus (Del Bigio, 2010). Previous studies in clinical and preclinical populations indicate thinning and stretching of the fimbria and fornix (Del Bigio et al., 2003; Naidich et al., 1982), which likely disrupts subcortical and collosal connectivity of the hippocampus (Del Bigio and Zhang, 1998; Mataro et al., 2006). Microstructural abnormalities of cellular deterioration in axons and synapses in hydrocephalic animal models further support the current interpretation of diminished cell density and structural organization of the hippocampus in the group with SBM (Kriebel and McAllister, 2000). Thus, the macro- and microstructural integrity of the hippocampus in SBM may reflect initial damage resulting from failure of neural tube closure coupled with prolonged damage resulting from hydrocephalus.

The putamen reflected a unique profile of increased volume, increased FA, and decreased MD in the group with SBM. The disproportionately increased volume and significantly reduced MD values observed in the putamen are not indicative of decreased density, but rather increased density. The literature of typical development from childhood to late adolescence and adulthood supports negative relations between age and MD values in the putamen (e.g., Abe et al., 2008; Lebel et al., 2008; Mukherjee et al., 2002). However, typical development of subcortical gray matter volumes generally peaks in early adolescence and decreases into young adulthood (Giedd et al., 1999a; Giedd et al., 1999b; Good et al., 2001; Ostby et al., 2009). When taken together with the larger volumes observed in the putamen, potential explanations for group differences in MD include atypical structural organization and aberrant apoptosis, or programmed cell death that typically occurs across adolescence. Mukerjee et al. (2002) suggested that deviations from isotropy could reflect deviations in structural patterning during neuronal and glial arborization, which supports our interpretation of increased FA and decreased MD in the current sample with SBM (Mukherjee et al., 2002). However, lack of research in early brain structure of infants and young children with SBM complicates elucidation of group differences. Knowledge of early brain macro- and microstructure in SBM would increase current understanding regarding structural and maturational anomalies in SBM.

Surprisingly, the current analyses also indicated distinct volumetric and microstructural patterning between the putamen and caudate in SBM. Structural properties of the caudate and the putamen (or dorsal striatum) were expected to be similar given the common embryonic derivation of these structures from the transient ventricular eminence (e.g., telencephalic origin) early in gestation (Marín and Rubenstein, 2003). In contrast, the globus pallidum emerges from diencephalic origins (Müller and O'Rahilly, 2006) and migrates to its final telencephalic position in the median eminence which together with the lateral eminence forms the floor of the third ventricle (Müller and O'Rahilly, 2006). Thus, structural properties of the caudate and the putamen are expected to be very similar since neurons within their boundaries are derived from the same origin. In contrast, structural properties of the pallidum are not expected to share structural similarities with the striatal nuclei since those neurons are derived from a different source of progenitor cells. Dissimilarity of macro- and microstructure in the caudate and putamen may reflect differences in their subcortical location, particularly to the ventricles. Decreased WM integrity within the posterior limb of the internal capsule in SBM may result in increased interstitial space in the extracellular neuropil (Kriebel et al., 1993; Ou et al., 2011). The putamen is surrounded by the internal and external capsules, which may anchor the putamen and reduce displacement following hydrocephalus. However, compromised WM integrity of the internal capsule in individuals with SBM may allow for the putamen to expand.

While no group differences were observed for the thalamus or caudate volumes, these periventricular structures exhibited a similar pattern of increased MD and FA in the group with SBM. Given the proximity of these structures to the ventricles, increased MD and FA may suggest lower density and cellular degeneration following the mechanical effects of hydrocephalus. Periventricular WM disruption (stretching) is the most frequent consequence of hydrocephalus and results in atrophy and gliosis of periventricular WM, increased extracellular space, and neuronal degeneration (Del Bigio, 2010; Del Bigio et al., 2003; Kriebel et al., 1993). Neuropathological studies in hydrocephalic rats indicate that hydrocephalus severity is positively related to thalamic apoptosis following disrupted periventricular WM (Del Bigio, 2010; Mori et al., 2002). Potentially, our observations of increased MD and FA in the thalamus and caudate of children and adolescents with SBM may be explained by less densely packed cells resulting from induced cell death from hydrocephalus. Our finding of increased FA in the caudate in SBM corroborates preliminary findings from a subset of this cohort (Hasan et al., 2008b), which was interpreted as reflecting reduced dendritic branching.

Limitations

Functional outcomes across individuals with SBM vary greatly and may be complicated by hydrocephalus and associated shunting. As previously mentioned, hydrocephalus disrupts neural structure, regardless of etiology (Del Bigio, 2010; Scott et al., 1998; Yeates and Enrile, 2005). The majority of children and adolescents with SBM in the current sample had hydrocephalus and had histories of shunt surgery. The current findings of nonsignificant relations between macro- and microstructural values of GM and number of shunt revisions should be interpreted cautiously. Shunting procedures were not uniform across individuals with SBM in the current sample and can disrupt brain structure during insertion and, when applicable, revision. Both the location and number of revisions could have impacted brain structure. Future research could benefit from an examination of how, specifically, shunt insertion (and revision) alters neural structure in this population and whether shunt location influences these outcomes. Therefore, current results represent a complex, and potentially synergistic, interplay between neural tube defect (i.e., SBM), hydrocephalus (resulting from the characteristic Chiari II malformation), and shunting history.

The current study examined a fairly broad age range, although statistical attempts to account for age were made. Although age was not strongly associated with regional subcortical GM volumes, it was significantly associated with the DTI metric measurements in the majority of analyses. However, for most models, group accounted for a greater amount of variance than age. Additionally, the relations observed between age and GM microstructure (negative for MD and positive for FA) support previous literature (Lebel et al., 2008; Lenroot et al., 2007; Ostby et al., 2009). Methodological weaknesses include the sample size of the comparison group and the use of a less advanced neuroimaging protocol relative to contemporary standards. A larger comparison group of TD individuals would be desirable especially given the broad age-range assessed. Future research may also benefit from a more advanced imaging protocol, including the utilization of greater angular resolution by increasing the number of diffusion directions for the DTI acquisition. However, we feel that the current data acquisition and processing were sufficient to address the current analytic questions and were state of the art at the time of data collection. Lastly, while we did not assess for puberty, subsequent developmental changes could have influenced findings, as subcortical GM development peaks after puberty (Lenroot et al., 2007). Future research would benefit from an investigation into the timing and hormonal effects of puberty on brain maturation in patients with SBM, particularly given reports of precocious puberty in this population (Meyer and Landau, 1984; Oakeshott and Hunt, 2003).

Implications and future directions

The current study examined the structural properties of deep GM in order to elucidate potential neural correlates underlying disrupted neurobehavioral outcomes. These findings are especially interesting given that these areas were largely read as normal in radiological review. Disruptions in higher-order cognition and motor functioning following SBM may be accounted for by aberrant development of subcortical GM. For instance, frontostriatal circuitry, which comprises the prefrontal, basal ganglia, and thalamic regions, is involved in multiple aspects of motor learning and control, including motor initiation, inhibition, and planning (Alexander et al., 1986). Additionally, the thalamus and certain basal ganglia nuclei, particularly the caudate and putamen, are shown to facilitate working memory and attentional control, which are impaired in individuals with SBM (Baier et al., 2010; Brown et al., 2008; Burmeister et al., 2005; Cohen and Frank, 2009; Rose and Holmbeck, 2007; Vachha and Adams, 2005). Furthermore, limbic system dysfunction, especially hippocampal, in SBM supports diminished memory abilities in this population (Dennis et al., 2007; Dennis et al., 2010; Hampton et al., 2011; Vachha and Adams, 2005; Vachha et al., 2006).

Overall, the present results support our concept of regionally specific brain areas being differentially impacted by this early gestational insult. However, the specific functional impairments associated with such structural abnormalities have yet to be identified. Examining neurobehavioral correlates of specific subcortical structures and the underlying behavioral differences observed in SBM would promote our understanding of these relations. Furthermore, while informative, macrostructural anomalies may only represent CNS damage associated with SBM, and do not explain complex interactions between regional microstructure in relation to global brain circuitry. Future research would significantly benefit from longitudinal investigations of subcortical GM development in children with SBM. Such investigations may illuminate mechanisms of plasticity and provide a foundation for effective intervention in this population.

Conflicts of interest

None to report.