Abstract
Neurofibromatosis type 1 (NF1) has regularly been associated with cognitive, social, and behavioral problems. The fact that many different cognitive and behavioral impairments have been observed in NF1 suggests that networks of brain regions are involved rather than specific brain regions. Here, we examined whether functional connectivity was different in NF1 and, if so, whether associations were present with cognitive, social, and behavioral outcomes. Fourteen NF1 patients (8 male, age: M=12.49, SD=2.65) and 30 healthy controls (HC; 23 male, age: M=12.30, SD=2.94; p=0.835) were included. Functional connectivity was assessed using functional resting-state scanning. We analyzed brain regions that have been associated with cognitive and social functions: the bilateral ventral anterior cingulate cortex (vACC), the bilateral amygdala, the bilateral orbitofrontal cortex (OFC), and the posterior cingulate cortex (PCC). For NF1 patients, connection strengths between brain regions showing HC-NF1 differences were correlated with parent reports of cognitive, social, and behavioral functioning. Compared to HC, patients showed differences in functional connectivity between the left vACC and the frontal cortex, insula, and subcortical areas (caudate, putamen), between the left amygdala and the frontal cortex, insula, supramarginal gyrus, and PCC/precuneus, and between the left OFC and frontal and subcortical areas (caudate, pallidum). In patients, indications were found for associations between increased frontofrontal and temporofrontal functional connectivity with cognitive, social, and behavioral deficits (r-range=0.536–0.851). NF1 patients showed differences in functional connectivity between areas associated with cognitive and social functioning when compared to controls. This, plus the fact that connectivity strengths in these networks were associated with worse cognitive, social, and behavioral outcomes, suggests a neuropathological basis for the widespread deficits observed in NF1.
Introduction
Neurofibromatosis Type 1 (NF1) is an autosomal genetic disorder, affecting one in 3500 newborns (National Institutes of Health, 1988), characterized by many different symptoms such as glial neoplasms, café-au-lait spots, skin fold freckling, optic nerve glioma, osseous (bone) abnormalities, iris Lisch nodules, and cutaneous and plexiform neurofibromas (Payne et al., 2010). Brain manifestations include tumors, megalencephaly, and anomalies of the corpus callosum (Payne et al., 2010).
Core behavioral consequences include cognitive impairment, social malfunction, and behavioral abnormalities. Cognitive abnormalities have been expressed in lower intelligence quotient (IQ), learning disabilities, visuospatial processing impairments, language problems, and executive deficits (Billingsley et al., 2004; Clements-Stephens et al., 2008; Coudé et al., 2006; Cutting et al., 2002; Hyman et al., 2005, 2006; North et al., 1994; Payne et al., 2011; Rowbotham et al., 2009; Roy et al., 2010). Social and behavioral problems include difficulties with social skills and communication, problems with forming friendships, rejection by peers, and social isolation. Although social problems have been associated with cognitive difficulties, low academic achievement, and attention deficit hyperactivity disorder (Barton and North, 2004; Gresham et al., 1996; Huijbregts et al., 2010), it was shown that social malfunction remained, although attenuated, after controlling for cognitive status (Huijbregts et al., 2010).
Previous functional magnetic resonance imaging (fMRI) studies in other diseases, examining alterations in networks of brain regions associated with cognitive and social functioning (Salmi et al., 2013; Wheaton et al., 2014), highlighted several pivotal areas being associated with measured malfunction, including the amygdala, orbitofrontal cortex (OFC), and anterior (ACC) and posterior cingulate cortex (PCC) (Gasquoine, 2013; Jackowski et al., 2012; Leech and Sharp, 2014; Murray et al., 2014; Salmi et al., 2013). The amygdala has been associated with determining the emotional significance of visual, auditory, and olfactory signals, and is involved in coping with the social environment and making social judgments (Adolphs et al., 1998), especially when evaluating facial expressions. Lesion studies revealed that amygdala damage results in deficits in retrieving socially relevant knowledge on the basis of facial appearance or an impairment in the ability to recognize fear from facial expressions (Adolphs, 1994, 1998). Deficits in face recognition, identification of facial emotions, especially fearful expressions, and in matching facial emotions have already been reported in NF1 (Huijbregts et al., 2010). The amygdala shares extensive reciprocal connections with the OFC, damage to which has also been associated with social deficits (Beer et al., 2006). The OFC is involved in social and antisocial behavior (Passamonti et al., 2010), self-monitoring (Beer et al., 2006), and social adaptation (Beyer et al., 2014). However, imaging studies of the OFC in NF1 are sparse. Violante et al. (2013) reported on local structural orbitofrontal differences between NF1 patients and control subjects and suggested these might underlie the executive impairments in NF1. The ACC and the PCC, both part of the default mode network (DMN, an assembly of synchronized spatially remote areas shown to be deactivated during goal-directed behavior), have been reported to be fundamentally involved in both social and cognitive processes (Greicius et al., 2003; Mars et al., 2012). The ACC comprises a ventral section (vACC), associated specifically with emotional aspects, especially the assessment of salience of emotional and motivational information and the regulation of emotional responses, and a dorsal section (dACC), which has been associated with modulation of attention and executive function (Bush et al., 2000).
In NF1, despite the disease-related impact on the brain [e.g., unidentified bright objects (Ferraz-Filho et al., 2012), macrocephaly (Steen et al., 2001), microstructural alterations (Ferraz-Filho et al., 2011)] and the established behavioral, cognitive, and social deficits, fMRI studies are scarce. A number of task-based functional imaging studies showed abnormal activation of brain regions associated with phonological and visual processing, as well as the effect of medication on the DMN (Billingsley et al., 2004; Chabernaud et al., 2012; Clements-Stephens et al., 2008; Violante et al., 2012), all reporting on task-related functional network changes.
Resting-state functional imaging has emerged as a valuable method of assessing integrity of brain functional connectivity networks in patient cohorts in the past decade, as data acquisition is less dependent on patient compliance, which appears to be particularly important in patient cohorts suffering from attention deficits, problems with task comprehension, or lying still during task execution. During resting-state fMRI, spontaneous fluctuations of the BOLD signal within remote brain regions are used to infer temporal coherencies (Biswal et al., 1995; Greicius et al., 2003). The temporal synchronization between remote brain regions is commonly interpreted as functional connectivity, which in turn is suggested to mirror brain functions (Smith et al., 2009).
Given the executive, social, and behavioral impairments in NF1, we hypothesized that these deficits would be associated with abnormal functional brain connectivity within brain networks related to these specific domains. Furthermore, we investigated whether any IQ differences were also associated with the functional connectivity differences of interest (and, if so, intended to correct for them). A seed-based connectivity approach, using seven seeds related to NF1 pathophysiology and sociocognitive problems (bilateral amygdala, bilateral OFC, bilateral vACC, and PCC), was used to provide a comprehensive overview of functional connectivity abnormalities in NF1.
Materials and Methods
Participants
In this cross-sectional study, 14 NF1 patients (8 male, age: M=12.49, SD=2.65) and 30 healthy controls (HCs, 23 male, age: M=12.30, SD=2.94; p=0.835, sex: χ2=1.748, df=1, p=0.186) underwent structural and functional MRI. All NF1 subjects fulfilled the diagnostic criteria specified by the National Institutes of Health (1988). Written informed consent was obtained.
The experiment was approved by the Ethics Committee of the Leiden University Medical Center (Leiden, The Netherlands) and the Department of Education and Child Studies, Leiden University (Leiden, The Netherlands).
Imaging protocol
All subjects underwent scanning at the Leiden University Medical Center. Imaging was performed on a Philips 3 Tesla Achieva MRI scanner using an 8-channel SENSE receiver head coil (Philips Healthcare, Best, The Netherlands). In each subject, a three-dimensional (3D) T1-weighted anatomical scan was acquired with the following scan parameters for registration purposes: 3D T1 TFE sequence, 140 axial slices, TR 9.8 msec, TE 4.6 msec, flip angle 8°, 1.16×0.92 mm, 1.2 mm slice thickness, and no slice gap. The scanning parameters of the resting-state fMRI scan were as follows: GE-EPI sequence, 160 volumes, 38 axial slices, TR 2.2 sec, TE 30 msec, flip angle 80°, 2.75 mm isotropic voxels, and 0.25 mm slice gap (duration=5.8 min). During the resting-state scan, participants were instructed to lie still with their eyes closed and not to fall asleep. All anatomical scans were reviewed by a neuroradiologist.
Preprocessing resting-state data
FMRI data preprocessing was performed using FMRIBs Software Library (FSL 5.0.4, Oxford, United Kingdom;
Resting-state data analyses
The functional connectivity analysis was performed using a dual regression method (part of FSL) (Filippini et al., 2009). We chose 4 mm spherical regions of interest in the vACC left (Montreal Neurologic Institute [MNI] coordinates: −2, 38, 10) and right (4, 36, 12), the amygdala left (−22, −6, −16) and right (28, −2, −20), the OFC left (−16, 12, −20) and right (20, 24, −20), and the PCC (0, −40, 30). The center coordinate of each seed was extracted from the Harvard-Oxford Cortical or Subcortical Atlases. Each seed was paired with two confound variables (white matter and cerebrospinal fluid) and used in the spatial and temporal dual regression technique [as described in Filippini et al. (2009)] to perform voxelwise between group comparisons. In the dual regression analysis, for each patient, each seed (and white matter and cerebrospinal fluid) was used as a spatial regressor in a general linear model to identify individual temporal dynamics (regression 1). The time courses were then used as a set of temporal regressors in a second general linear model (regression 2), resulting in beta-weights for each voxel. These were used in the subsequent higher level analysis. Group-level inferences were carried out using voxelwise nonparametric permutation testing (5000 permutations, Randomise v2.9, Nichols and Holmes, 2002). Results were corrected for multiple comparisons across space using cluster-based correction (threshold z=2.3, FDR-corrected p=0.05). Group comparisons were performed using two-tailed independent samples t-tests. To identify the direction of functional connectivity changes, we extracted the mean beta-values for each individual within clusters of significant difference of group comparisons.
Questionnaires
The neuropsychological assessment involved estimation of IQ based on the Wechsler Intelligence Scale for Children (WISC, Wechsler et al., 2011). The following subtests were obtained: verbal comprehension (vocabulary, comprehension), visual-spatial abilities (block design), fluid reasoning (picture completion), and processing speed (symbol search, coding). Verbal IQ was estimated based on the subscales vocabulary and comprehension and performance IQ was estimated on the basis of all other subtests. For other aspects of cognitive and social functioning, several questionnaires were administered. For social skills, the parent version of the Social Skills Rating System (SSRS) (Gresham and Elliott, 1990) was used. The four SSRS-dimensions (cooperation, assertion, responsibility, self-control) as well as the total SSRS score were used as indicators of social skills. The Social Responsiveness Scale (SRS) (Constantino, 2005) is a standard measure of autism spectrum. It includes five interrelated dimensions: social awareness, social cognition, reciprocal social communication, social motivation, and autistic mannerisms. The Child Behavior Checklist (CBCL) (Achenbach, 1992) was used to determine behavioral problems. Executive functioning was assessed using the parent-rated Behavior Rating Inventory of Executive Function (BRIEF) (Gioia et al., 2000) and the Dysexecutive Questionnaire (DEX, Wilson et al., 1996). The DEX is part of the Behavioral Assessment of the Dysexecutive Syndrome (Wilson et al., 1996) and is designed to assess a range of cognitive, behavioral, and emotional problems associated with executive dysfunction.
The BRIEF total score and the DEX cognition were used as indicators of executive function. Due to a high correlation between scores on the two questionnaires (r=0.836, p<0.001), principal component analyses for dimensionality reduction were used to extract one fixed component. To account for the reported influence of cognition on social functioning, this component was used as a covariate in the group comparison of social psychometrics to control for the effect. Test scores were available for all patients and nine control children.
Statistical analyses
SPSS 21 for Windows (SPSS, Chicago, IL) was used for statistical analyses. Whenever appropriate, tested with the Kolmogorov–Smirnov test for normal distribution of data, parametric tests were used. If data were not normally distributed, nonparametric tests were applied. For the NF1 patients, Pearson correlations were performed (uncorrected for multiple comparisons) between the beta-scores of areas of functional connectivity differences and IQ, executive, social, and behavioral outcome measures.
Results
IQ, cognitive and social profile
When compared to controls, NF1 patients showed a significantly lower performance IQ (MHC=11.36, SDHC=2.05; MNF1=8.91, SDNF1=2.29; p=0.016) and total IQ scores (MHC=11.17, SDHC=1.88; MNF1=8.92, SDNF1=1.77; p=0.009). Only a trend was found for verbal IQ (MHC=10.78, SDHC=2.43; MNF1=8.93, SDNF1=2.23; p=0.074). Also, NF1 patients showed poorer cognitive and social scores in all but one examined domains (Table 1). As social functioning was reported to be influenced by cognitive control deficits (Huijbregts et al., 2010), subscales of social questionnaires were corrected for executive function. The SRS total score, including the subscale social motivation and autistic mannerisms, as well as SSRS assertion and CBCL social problems, remained significantly different after correction (Table 1).
Group Differences Between NF1 Patients and Healthy Controls on Cognitive/Executive Function and Social Questionnaires
As executive functions can impact on social skills, p-corr shows group differences, in social scales only, after correction for cognition/executive function (uncorrected for multiple comparisons).
Higher scores indicate better performance.
Mann–Whitney U-Test, otherwise t-test.
DEX, Dysexecutive Questionnaire; BRIEF, Behavior Rating Inventory of Executive Functioning; SRS, Social Responsiveness Scale; SSRS, Social Skills Rating System; CBCL, Child Behavior Checklist; uncorrected for multiple comparisons; NF1, neurofibromatosis type 1.
Seed-based functional connectivity
Mean functional connectivity maps for the seed regions and their associated areas can be seen in Figure 1 (1a, 2a, 3a, 4–7). For a precise group-specific network description, see Supplementary Data; Supplementary Data are available online at

Mean functional connectivity maps of the seed regions in the vACC (
To determine whether functional connectivity is higher or lower in NF1 patients when compared to controls between the respective areas (e.g., between the vACC left and the frontal pole right), beta-scores served to indicate on directionality (Figs. 1 and 2). Functional connectivity in the abovementioned areas of group differences was not significantly different from zero in controls, highlighting that within the patient group, positive or negative connectivity exists between areas where no significant connectivity on average in controls is observed. Therefore, all regions with a difference show an increased connection (either positive or negative interaction) in NF1.

Mean beta-values of NF1 children and controls in areas of connectivity difference (peak voxel) seen in Figure 1. NF1 patients show positive coupling between the vACC left and the frontal pole
To examine the association between differences in functional connectivity and executive, social, and behavioral outcomes, as well as IQ, correlations between betas of areas of significant differences (i.e., frontal pole right, OFC right, PCC/precuneus, and the pallidum left, Fig. 2) and questionnaire scores were conducted (Table 2). These revealed that higher functional connectivity between (1) the left vACC and the right frontal pole correlated positively with the DEX, BRIEF, and SRS scores, and subscores of the CBCL and (2) a higher functional connectivity between the left amygdala and the right OFC correlated with subscores of the CBCL. These results indicate that the higher the functional connectivity between the stated areas, the worse the executive, social, and behavioral performance of NF1 patients. No correlation was found between areas of lower functional connectivity (PCC/precunus, pallidum left) and questionnaire outcome measures. Also, no correlation was found with IQ.
Significant Correlations Between Areas of Changed Functional Connectivity and Social/Executive Measures in NF1 Patients Only (Uncorrected for Multiple Comparisons)
Significant after correction for multiple comparisons (FDR, q=0.05).
L, left; r, right; vACC, ventral anterior cingulate cortex; FP frontal pole; AMG, amygdala; OFC, orbitofrontal cortex.
Discussion
Using resting-state fMRI to investigate functional connectivity alterations within networks associated with cognitive, social, and behavioral performance, we observed several significant functional alterations in the connectivity patterns of NF1 children when compared to controls. These differences revealed increases of interactions between certain brain regions in NF1 patients.
Positive coupling was found between left vACC and areas of the frontal cortex, insula cortex, and subcortical areas, and between the left amygdala and the frontal cortex, parietal cortex, and insula cortex. Negative coupling was found between the left amygdala and the PCC/precuneous, and between left OFC and frontal and subcortical areas. These changes were partly associated with behavioral outcome measures in NF1 patients. The positive coupling between the left vACC and the frontal pole and between the left amygdala and the right OFC was correlated with worse executive, social, and behavioral performances in NF1 patients. Importantly, equivalent correlations between functional connectivity and IQ have not been found, corroborating that deficits in executive functions, social skills, and behavior cannot be ascribed to global deficits only. We conclude that the increased functional connectivity between frontofrontal areas and temporofrontal areas is maladaptive rather than a failed attempt of the brain to compensate by functional reorganization.
In addition to the functional connectivity changes, we found in line with previous literature (Huijbregts and de Sonneville, 2011; Rowbotham et al., 2009), widespread impairments in executive, social, and behavioral function in the patients. The social deficits partly remained even after controlling for executive dysfunction (Huijbregts and de Sonneville, 2011; Huijbregts et al., 2010). Together, these findings further inform us on the neuropathological basis for the executive, social, and behavioral phenotype that has consistently been associated with NF1.
The seed-based approach revealed NF1-related connectivity deviations of the left vACC, left amygdala, and left OFC. The most prominent deviations were found in both the left vACC and left amygdala showing functional increases in connectivity with the insular cortex and frontoparietal operculum. The insular cortex has been associated with a variety of different functions, such as cognitive control (Cole and Schneider, 2007), individual emotions (Nieuwenhuys, 2012; Phan et al., 2002), and social emotion, including empathy and compassion, as well as interpersonal phenomena such as fairness and cooperation (Lamm and Singer, 2010). Also, the importance of the insular cortex has already been shown in several diseases associated with social-cognitive dysfunction such as mood disorders, post-stroke depression, panic disorders, post-traumatic stress disorders, obsessive compulsive disorder, and schizophrenia (Nagai et al., 2007). Hence, a shift in functional connectivity of the vACC and the amygdala with the insular cortex might be a strong indicator of underlying pathophysiological change in the context of social and cognitive status in NF1.
In addition, differences in functional connectivity were identified between the vACC and the dACC and the frontal pole. Previous literature highlighted the segregated functional role of the ventral and dorsal cingulate cortex, ascribing them an emotional and cognitive function, respectively (Bush et al., 2000). A reciprocal suppression between these two cingulate areas was reported, that is, deactivation of the vACC during cognitive tasks and deactivation of the dACC during emotional tasks (Bush et al., 2000). Furthermore, the dACC has strong connections to the frontal pole, which has been associated with social cognition, attention, episodic memory, and multitasking (Burgess et al., 2007; Gilbert et al., 2006, 2010). The imbalance between dorsal and ventral suppression within the anterior cingulate in adolescents with NF1, seen in their increased functional connectivity when compared to controls, together with the increased connectivity to the frontal pole, might explain the cognitive influence on social functioning in NF1 children, and hence, the attenuation of social deficits when controlling for cognitive dysfunction. Additionally, a divergence in the functional connectivity of the ACC, together with the medial prefrontal cortex and the retrosplenial cortex, all part of the DMN, has already been reported in a prior study examining visual processing of NF1 patients (Violante et al., 2012). The authors reported that patients failed to deactivate these areas when compared to controls. As the DMN has been associated with daydreaming and self-referential thought, they concluded that patients had a higher tendency toward task-irrelevant thoughts. Although we cannot corroborate this conclusion with our results, the repeated finding of functional differences between patients and controls associated with the ACC indicates its significant role in NF1-related brain pathology.
Finally, a negative coupling in functional connectivity between the amygdala and the PCC/precuneus was found. Both the PCC and precuneus are main hubs of the DMN (Damoiseaux et al., 2006). Previously, altered connectivity between regions of the DMN has been associated with a variety of mental and developmental disorders, including schizophrenia, autism, attention deficit, and anxiety disorders (Leech and Sharp, 2014). The fact that patients show more negative coupling between the amygdala and the PCC/precuneus could be indicative of a frontoparietal pathophysiological process. The network associated with the PCC seed did not reveal any group differences. Although this seems counterintuitive, considering the wide range of abnormalities in NF1 and the disease-related sensibility of the DMN, seen in other diseases (Leech and Sharp, 2014), outcomes of networks depend on the precise location of the seed. We therefore cannot exclude that various other seed coordinates within the PCC might reveal connectivity differences between groups. Further studies are needed to examine its potential role in NF1.
A limitation of the current study is the modest size of the patient group, although similar group sizes have been used in other fMRI studies in NF1 (Chabernaud et al., 2012; Violante et al., 2012). Behavioral data were available for only a subset of control subjects, although group differences were very much in line with those observed in earlier studies with larger control samples. Another limitation concerns the fact that results are uncorrected for multiple comparisons. Differences in functional connectivity strongly indicate on deviant brain organization in NF1 patients, but due to limited power, results become nonsignificant after correction. Also, the correlation results remained uncorrected. It therefore needs to be emphasized that the results have to be interpreted with care and need further verification in bigger samples.
Despite these limitations, this study provides comprehensive insights into the role of abnormal functional connectivity and its association to executive, social, and behavioral deficits observed in NF1.
Footnotes
Acknowledgments
This study was funded by the Austrian Science Fund (FWF: J3500 B23) and the Dutch Neurofibromatosis Association. S.A.R. is supported by a grant of the Netherlands Organization for Scientific Research (NWO).
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
