Abstract
Psychopathy is a personality disorder characterized by callous antisocial behavior and criminal recidivism. Here we examine whether psychopathy is associated with alterations in functional connectivity in three large-scale cortical networks. Using fMRI in 142 adult male prison inmates, we computed resting-state functional connectivity using seeds from the default mode network, frontoparietal network, and cingulo-opercular network. To determine the specificity of our findings to these cortical networks, we also calculated functional connectivity using seeds from two comparison primary sensory networks: visual and auditory networks. Regression analyses related network connectivity to overall psychopathy scores and to subscores for the “factors” and “facets” of psychopathy: Factor 1, interpersonal/affective traits; Factor 2, lifestyle/antisocial traits; Facet 1, interpersonal; Facet 2, affective; Facet 3, lifestyle; Facet 4, antisocial. Overall psychopathy severity was associated with reduced functional connectivity between lateral parietal cortex and dorsal anterior cingulate cortex. The two factor scores exhibited contrasting relationships with functional connectivity: Factor 1 scores were associated with reduced functional connectivity in the three cortical networks, whereas Factor 2 scores were associated with heightened connectivity in the same networks. This dissociation was evident particularly in the functional connectivity between anterior insula and dorsal anterior cingulate cortex. The facet scores also demonstrated distinct patterns of connectivity. We found no associations between psychopathy scores and functional connectivity within visual or auditory networks. These findings provide novel evidence on the neural correlates of psychopathy and suggest that connectivity between cortical association hubs, such as the dorsal anterior cingulate cortex, may be a neurobiological marker of the disorder.
- antisocial personality disorder
- fMRI
- fronto-parietal network
- functional connectivity
- psychopathy
- resting-state
Introduction
Psychopathy is a personality disorder characterized by a callous lack of empathy and impulsive antisocial behavior. Present in approximately one-quarter of adult prison inmates, psychopathy is a significant predictor of violent crime and recidivism (Harris et al., 1991; Hare, 2003). Thus, identifying the psychological and neurobiological mechanisms underlying this disorder could have profound implications for the clinical and legal management of psychopathic criminals, as well as for the basic understanding of human social behavior. Initial research on the neural correlates of psychopathy has demonstrated abnormalities in the structure and function of numerous cortical and subcortical brain regions, with the predominance of findings highlighting frontal, temporal, and limbic areas involved in social and affective processing (Yang and Raine, 2009; Blair, 2010; Koenigs et al., 2011; Anderson and Kiehl, 2012). However, the neural correlates of the disorder have yet to be identified comprehensively.
One major knowledge gap in this area of research concerns the role that alterations in cortical network connectivity may play in psychopathy. In recent years, psychiatric illnesses have come to be understood as disorders of neural circuits (Morris and Cuthbert, 2012). This perspective, which underscores large-scale network connectivity as a critical biological parameter related to mental health, is supported by a rapidly expanding body of neuroimaging research that associates a variety of psychiatric disorders with altered resting-state connectivity in distinct cortical networks (Menon, 2011; Sylvester et al., 2012; Whitfield-Gabrieli and Ford, 2012). In particular, Menon (2011) suggests that dysfunction in connectivity within and between three major cortical networks—the default mode network, central executive/frontoparietal network, and salience/cingulo-opercular network—may underlie various psychiatric and neurological disorders. However, to date, no study has examined systematically whether psychopathy is associated with alterations of resting-state functional connectivity within and between these three major cortical association networks. Given that psychopathy involves severe deficits in the integration of social and affective information into decision-making (Cleckley, 1976; Newman and Lorenz, 2003)—a function that would presumably depend critically on patent connectivity among these same cortical networks—there is strong reason to hypothesize alterations in the resting-state connectivity in these networks. This study is the first to test this hypothesis in a large sample of criminal offenders with a wide range of psychopathy severity.
Furthermore, the large and diverse inmate sample affords us a unique opportunity to examine the neural correlates of the two primary “factors” and the four “facets,” or dimensions, of psychopathy. Factor 1 corresponds to the interpersonal/affective traits of psychopathy (e.g., callousness, egocentrism), whereas Factor 2 corresponds to the lifestyle/antisocial features (e.g., impulsivity, irresponsibility; Hare et al., 1990). Factor 1 can be further subdivided into Facet 1 (interpersonal traits) and Facet 2 (affective traits), whereas Factor 2 can be further subdivided into Facet 3 (lifestyle traits) and Facet 4 (antisocial traits; Hare and Neumann, 2005). Although several previous neuroimaging studies have identified neural correlates of individual psychopathy factors and/or facets (Yang et al., 2005, 2009; de Oliveira-Souza et al., 2008; Ermer et al., 2012; Contreras-Rodríguez et al., 2014), none of these studies have examined resting-state functional connectivity across such a broad range of psychopathy severity. In summary, this study uses resting-state fMRI in a large sample of adult prison inmates to identify novel neural correlates of psychopathy.
Materials and Methods
Participants.
One hundred fifty-five adult male inmates, recruited from a medium-security Wisconsin correctional facility, participated in the present study. Thirteen subjects were excluded because of excessive head motion during the scan (see below, Preprocessing of functional data). Participants were selected based on the following inclusion criteria: (1) age <45 years; (2) intelligence quotient (IQ) >70; (3) no history of psychosis or bipolar disorder; (4) no history of significant head injury or postconcussion symptoms; and (5) no current use of psychotropic medications. Informed consent was obtained both orally and in writing.
Psychopathy was assessed for all participants using the Psychopathy Checklist Revised (PCL-R; Hare, 2003). The PCL-R is a 20-item scale. Scores for each item were rated as 0, 1, or 2 based on the severity of each trait determined from a semi-structured interview and file review administered by trained research assistants. Inter-rater reliability (intraclass correlation) for total PCL-R score was 0.98 based on 10 dual ratings. Total PCL-R, Factor 1 (interpersonal/affective traits), Factor 2 (lifestyle/antisocial traits), Facet 1 (interpersonal traits), Facet 2 (affective traits), Facet 3 (lifestyle traits), and Facet 4 (antisocial traits) scores were used for separate regression analyses. For the group analyses, participants were classified as psychopathic if their PCL-R scores were ≥30 (n = 46) and nonpsychopathic if their PCL-R scores were ≤20 (n = 49; Hare, 2003). For demographic and neuropsychological information, see Table 1.
Participant characteristics
MRI data acquisition.
All MRI data were acquired on correctional facility grounds with the Siemens 1.5T Avanto Mobile MRI System of the Mind Research Network with advanced SQ gradients (maximum slew rate, 200 T · m−1 · s−1; 346 T · m−1 · s−1 vector summation; rise time, 200 μs) equipped with a 12-element head coil. Head motion was limited using padding and restraint.
Resting-state functional images were collected while subjects lay still and awake, passively viewing a fixation cross for 5.5 min (158 volumes). T2*-weighted gradient-echo functional echo planar images (EPIs) were acquired with the following parameters: repetition time (TR), 2000 ms; echo time (TE), 39 ms; flip angle, 75°; field of view (FOV), 24 × 24 cm; matrix, 64 × 64; slice thickness, 4 mm; gap, 1 mm; voxel size, 3.75 × 3.75 × 5 mm; 27 sequential axial oblique slices. A T1-weighted structural image was acquired using a four-echo MPRAGE sequence (TR, 2530; TE, 1.64, 3.5, 5.36, 7.22 ms; flip angle, 7°; FOV, 256 × 256 mm; matrix, 128 × 128; slice thickness, 1.33 mm; voxel size, 1 × 1 × 1.33 mm; 128 interleaved sagittal slices). All four echoes were averaged into a single high-resolution image that was used to aid in the spatial normalization of EPI volumes and visualization of group statistics.
Preprocessing of functional data.
All fMRI data analysis was performed using AFNI (Automated Functional Neuro-Imaging; Cox, 1996) and FSL (FMRIB Software Library; http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/). The following preprocessing steps were performed: (1) EPI volumes were slice time corrected using the first slice as a reference and (2) motion corrected by rigid body alignment to the first EPI acquisition, (3) images were deobliqued, (4) the first three volumes were omitted, (5) data were then motion corrected (3dvolreg function in AFNI) and (6) despiked to remove extreme time series outliers and then (7) bandpass filtered (0.009 < f < 0.08) and spatially smoothed with a 6 mm full-width at half-maximum Gaussian kernel (Fox et al., 2005). The skull-stripped anatomical scan for each participant was rigidly coregistered with the EPI and diffeomorphically aligned to Montreal Neurological Institute (MNI)-152 template space using a symmetric normalization algorithm (Avants and Gee, 2004). The transformation matrix from this registration procedure was then used to align the EPI scans for each participant to MNI-152 space. Finally, the EPI scans were resampled to 3 mm cubic voxels for subsequent functional connectivity analyses.
Because individual differences in subject motion can contribute to resting-state correlations (Power et al., 2012; Satterthwaite et al., 2013; Yan et al., 2013), we examined motion for each subject. We excluded subjects with mean framewise motion displacement (i.e., volume to volume movement across the time series) >2 mm and/or total scan time <4 min after censoring all time points with framewise motion displacement >0.2 mm and extreme time series displacement (i.e., time points in which >10% of voxels were outliers; Power et al., 2012; Satterthwaite et al., 2013; Yan et al., 2013). These specific thresholds for motion (mean framewise motion displacement <2 mm and >4 min remaining after censoring time points for both framewise motion >0.2 mm and extreme time series displacement >10% of voxels) were selected based on previous research to provide the most contemporary and conservative criteria for motion correction. Thus, 13 subjects were excluded from additional analysis (n = 3 with mean framewise motion >2 mm, n = 13 with total scan time <4 min after motion censoring). In the remaining sample (after censoring), there were no significant differences between groups in root mean squared head-position change [psychopathic, 0.138 ± 0.08 (mean ± SD); nonpsychopathic, 0.122 ± 0.06; p = 0.25; Power et al., 2012] or in the number of frames remaining (psychopathic, 147.83 ± 8.97; nonpsychopathic, 148.43 ± 8.28, p = 0.73).
Normalized T1 anatomical images were segmented into gray matter, white matter, and CSF segments using FAST (FMRIB Automated Segmentation Tool) in FSL (Zhang et al., 2001). White matter and CSF segments were used as masks to extract a representative time series from each tissue type.
Region of interest selection for cortical networks.
To assess functional connectivity for three main cortical networks of interest [default mode (DMN), frontoparietal (FPN), and cingulo-opercular (CON)], we selected a priori seed regions of interest (ROIs) that have been identified from previous studies (Andrews-Hanna et al., 2007; Dosenbach et al., 2007; Table 2). To determine the specificity of our findings to these cortical networks, we also selected two frequently used primary sensory comparison networks [visual (VN) and auditory (AN); Andrews-Hanna et al., 2007; Boveroux et al., 2010]. Given the aim of the study on connectivity within and between DMN, FPN, and CON, we focus our results on clusters that fall within these networks. We defined the regions within these networks based on a recent large-scale resting-state study (n = 1000 subjects) that demonstrated the reliability and reproducibility of these same cortical networks (Yeo et al., 2011). Specifically, we used masks for the DMN (network 7), FPN (network 6), dorsal attention network (network 3), and CON (network 4) (http://www.freesurfer.net/fswiki/CorticalParcellation; Yeo et al., 2011). Each network mask was aligned to MNI-152 space and resampled to 3 mm3 resolution. Based on these network parcellations, the following regions were included in each of our cortical networks of interest if they overlapped with the network masks mentioned above: (1) DMN: the medial prefrontal cortex (mPFC), the superior frontal gyrus extending into superior portions of the middle frontal gyrus, the posterior cingulate cortex (PCC)/retrosplenial cortex, the inferior and medial lateral temporal gyri, and the inferior parietal lobule (angular gyrus); (2) FPN: the superior parietal lobule [including the intraparietal sulcus (IPS)], the superior, middle, and inferior frontal gyri [including the dorsolateral PFC (dlPFC)], the precuneus, and the inferior temporal cortex; (3) dorsal attention network: the precuneus (the superior portion bilaterally); and (4) CON: the dorsal anterior cingulate cortex (dACC)/mesial frontal cortex (msFC), the anterior insula/frontal operculum (aIfO), the mid-insula, and the supramarginal gyrus. The regions included in the comparison networks were the following: (1) VN: the lingual gyrus and the cuneus; and (2) AN: the primary auditory cortex/superior temporal gyrus.
Cortical networks evaluated for resting-state functional connectivity
Functional connectivity analysis.
To evaluate functional connectivity for each of the seed ROIs, we generated a spherical 6 mm radius seed centered on the coordinates for each ROI using the Wake Forest University PickAtlas (Maldjian et al., 2003). All seeds were aligned to MNI-152 space.
Functional connectivity was assessed for each network seed ROI using the mean resting-state BOLD time series, extracted for each participant. The mean time series from each ROI was included in a GLM with 15 regressors of no interest: 1–12, six motion parameters (three translations and three rotations) obtained from the rigid-body alignment of EPI volumes and their six derivatives; 13, the white matter time series; 14, the ventricular (CSF) time series; and 15, a second-order polynomial to model baseline signal and slow drift. To further control for subject motion during the GLM, volumes were censored for extreme time series displacement (i.e., time points in which >10% of voxels were outliers) and framewise motion displacement (i.e., volume to volume movement) >0.2 mm (Power et al., 2012; Yan et al., 2013).
To convert the output of R2 values from the GLM to correlation coefficients (r), the following equation implemented in AFNI was used: 3dcalc -a Corr_subject1.nii ‘[sub-brick for R2]’ -b Corr_subject1.nii “[sub-brick for Beta]” -expr ‘ispositive(b) × sqrt(a) − isnegative(b) × sqrt(a)’. Then, correlation coefficients for each network seed ROI were converted to z scores via Fisher's r to z transform and corrected for degrees of freedom using the following equation: z = log ((1 + r)/(1 − r)) × √((n − 3))/2, where n is degrees of freedom (i.e., time points in the fMRI resting-state scan after censoring). The resulting z-score maps were entered into second-level statistical analyses.
Regression analyses for seed-based functional connectivity.
We performed linear regression analyses (3dRegAna function in AFNI) to examine the relationship between PCL-R scores and functional connectivity for all cortical network seed ROIs. We performed separate regressions for each of the following PCL-R scores: total PCL-R, Factor 1 (covarying for Factor 2), Factor 2 (covarying for Factor 1), Facet 1, Facet 2, Facet 3, and Facet 4. To correct for multiple comparisons, we used familywise error (FWE) correction at the cluster level using a whole-brain mask (3dClustSim function in AFNI; Forman et al., 1995; Carp, 2012) and applied cluster extent thresholding to our regression results. The cluster extent threshold corresponded to the statistical probability (α = 0.05, or 5% chance) of identifying a random noise cluster at a predefined voxelwise (i.e., whole-brain) threshold of p < 0.005 (uncorrected). Using this whole-brain FWE cluster correction, a cluster-corrected size of ≥72 voxels was significant at pFWE < 0.05 in the regression analyses reported below for PCL-R scores. Regression results were overlaid on the normalized mean anatomical image.
Follow-up regression analyses were performed for all significant results to account statistically for potential differences in age, substance use disorder (SUD), race, and IQ. We first extracted the average z score from all significant clusters identified for each subject, and then we ran separate regressions (in R) to examine the association between the average z score for each significant cluster and the PCL-R score after controlling for age, SUD, race, and IQ. Regressions with PCL-R factor scores also included partial correlations (e.g., models with Factor 1 also contained Factor 2, along with age, SUD, race, and IQ). Regressions with PCL-R facet scores included the facet score after controlling for age, SUD, race, and IQ. Importantly, all results were corrected for multiple tests using the Holm–Bonferroni correction method (Holm, 1979). Using the Holm–Bonferroni method, we performed the following steps: (1) ordered the p values for all regression results from smallest to largest, (2) multiplied the smallest p value by [k − (n − 1)], where k is the number of null hypothesis tests or contrasts and n is the test number (i.e., n = 1 for the first p value in the ordered list), (3) multiplied the next smallest p value by [k − (n − 1)] (i.e., n = 2 for the second p value), and (4) continued this iterative process for all remaining p values until the Holm–Bonferroni corrected p value exceeded 0.05. Below we only report results that remained significant after the inclusion of relevant covariates (age, SUD, race, IQ, and Factor 1/Factor 2) and Holm–Bonferroni correction, at p < 0.05.
Group comparisons for seed-based functional connectivity.
As a follow-up to the seed-based regression analyses, we investigated differences in functional connectivity between psychopathic and nonpsychopathic inmates. We performed unpaired two-sample t tests on the z-score maps derived from each seed ROI from the seed-based functional connectivity analysis. As described above for the regressions, we used an FWE correction approach at the cluster level using a whole-brain mask to correct for multiple comparisons at the voxelwise level (Forman et al., 1995; Carp, 2012). In this group comparison, a cluster size of ≥89 voxels was significant at pFWE < 0.05 in the two-sample t tests reported. Group correlation maps were overlaid on the normalized mean anatomical image.
Results
Regression analysis: total PCL-R
For the FPN (Fig. 1, Table 3), higher total PCL-R scores were associated with significantly diminished connectivity between the right (R) IPS and the bilateral dACC of the CON and right precuneus of the FPN. There were no significant relationships between total PCL-R scores and functional connectivity for any seeds of the DMN, CON, or either of the comparison networks.
Total PCL-R scores were related inversely to functional connectivity between the R IPS and the bilateral dACC and the right precuneus. Top row, Seed ROIs (red–orange) for the R IPS displayed significant regression results (pFWE < 0.05) for total PCL-R scores and R IPS connectivity. All results and seed ROIs are displayed on the group average structural MRI in MNI space. Bottom row, Scatter plots show the relationship between total PCL-R scores and the connectivity values for the R IPS and the bilateral dACC (A) and the R IPS and the right precuneus (pC; B). Each regression line (plotted in black) depicts the predicted model for total PCL-R and R IPS connectivity including covariates (age, SUD, race, and IQ); confidence intervals are plotted in blue lines.
Resting-state functional connectivity results: total PCL-R and factor scores
Between-group analysis: psychopathic versus nonpsychopathic inmates
Similar to the regression analyses with total PCL-R, in the between-group analysis, we found significant differences for only FPN seeds. Psychopathic inmates exhibited significantly reduced functional connectivity between the R IPS and right precuneus (within the FPN) and between the R IPS and left (L) dACC (extending bilaterally) and the right insula of the CON (Fig. 2, Table 3). There were no significant group differences for either of the two comparison networks.
Group differences in functional connectivity between the R IPS, the L dACC, the right precuneus, and the right insula. Top row, R IPS connectivity for the psychopathic group (Psych(P)). Middle row, R IPS connectivity for the nonpsychopathic group (Non-psych(NP)). Bottom row, Group difference maps (Diff(P-NP)) show significantly diminished connectivity between the R IPS and the left dACC extending bilaterally (x = −1; i), the right precuneus (x = 8; ii), and the right insula (x = 38; iii) in psychopathic offenders compared with nonpsychopathic offenders (pFWE < 0.05). Results are displayed on the group average structural MRI in MNI space. Color bars indicate uncorrected t values.
Regression analysis: PCL-R factor scores
For the DMN (Fig. 3, Table 3), the individual factor scores were associated with distinct patterns of functional connectivity. Higher Factor 1 scores were associated with significantly decreased connectivity between the medial and lateral regions of the DMN [PCC seed: right superior frontal gyrus; L inferior parietal lobule (IPL) seed: left superior frontal gyrus; R IPL seed: left posterior cingulate]. In contrast, higher Factor 2 scores were associated with significantly increased functional connectivity between the medial prefrontal and parietal regions of the DMN (mPFC seed: left dorsomedial PFC, left superior frontal gyrus, right posterior cingulate; PCC seed: left dorsomedial PFC, left superior frontal gyrus extending to middle frontal gyrus). In addition, higher Factor 2 scores were associated with significant functional connectivity between the mPFC seed and left precuneus (within the FPN).
Distinct patterns of functional connectivity for PCL-R factor scores. Left column, Seed ROIs (red–orange) with significant regression results are displayed for each network. All results and seed ROIs are displayed on the group average structural MRI in MNI space. Right column, Conjunction maps are presented to display the overlap of significant regression results (pFWE < 0.05) for both PCL-R factor scores for the DMN (rows 1–4), FPN (rows 5–7), and CON (rows 8–10) seeds. Color codes for PCL-R factor scores: light blue, Factor 1 (F1, negative relationship); yellow, Factor 2 (F2, positive relationship); green, overlap between F1 and F2 results.
For the FPN (Fig. 3, Table 3), higher Factor 1 scores, similar to total PCL-R, were associated with significantly decreased connectivity between the lateral parietal, medial parietal, and lateral prefrontal regions of the FPN (L IPS seed: right superior frontal gyrus; R IPS seed: left precuneus; L dlPFC seed: left superior parietal lobule). In contrast, higher Factor 2 scores were associated with significantly increased functional connectivity between lateral prefrontal regions of the FPN (L dlPFC seed: right middle frontal gyrus) and right supramarginal gyrus of the CON (L dlPFC seed).
For the CON (Fig. 3, Table 3), as in the other two networks, the individual factor scores were again associated with distinct patterns of functional connectivity. Higher Factor 1 scores were associated with decreased connectivity between the insula, anterior cingulate, and lateral parietal regions of the CON (L aIfO seed: right dorsal anterior cingulate extending bilaterally, right anterior insula, left supramarginal gyrus; R aIfO seed: right dorsal anterior cingulate, right inferior frontal gyrus/insula, left inferior frontal gyrus/insula, right postcentral gyrus extending to supramarginal gyrus), whereas higher Factor 2 scores were associated with increased connectivity between the insula, anterior cingulate, and lateral parietal regions of the CON (L aIfO seed: right mid-insula; R aIfO seed: left dorsal anterior cingulate extending bilaterally, right supramarginal gyrus, right inferior frontal gyrus/insula, right medial frontal gyrus, left inferior frontal gyrus). Higher Factor 2 scores were also associated with increased connectivity between the dACC/msFC seed and right precuneus of the FPN.
Similar to total PCL-R scores, we found no significant associations between factor scores and functional connectivity for the comparison primary sensory networks.
Regression analysis: PCL-R facet scores
For the DMN (Fig. 4, Table 4), higher Facet 1 (interpersonal) scores were associated with significantly diminished functional connectivity between the R IPL seed and medial prefrontal and medial parietal regions of the DMN (right ventromedial PFC, right posterior cingulate). There were no significant associations between the other three facet scores and functional connectivity for the DMN.
Distinct patterns of functional connectivity for PCL-R facet scores. A, Left, R IPL seed. Right, Higher Facet 1 (interpersonal) scores were associated with reduced connectivity between the R IPL and regions of the DMN. B, Left, R IPS seed. Right, Higher Facet 1 (interpersonal) scores were associated with reduced connectivity between the R IPS, pC, and dACC. C, Left, R IPS seed. Right, Higher Facet 3 (lifestyle) scores were associated with reduced connectivity between the R IPS and the superior frontal gyrus. D, Left, R IPS seed. Right, Higher Facet 4 (antisocial) scores were associated with greater connectivity between R IPS and regions of the FPN but reduced connectivity with insula. E, Left, dACC/msFC seed. Right, Higher Facet 1 (interpersonal) scores were associated with reduced functional connectivity between the dACC/msFC and IPL. F, Left, R aIfO seed. Right, Higher Facet 2 (affective) scores were associated with diminished connectivity between the R aIfO and supramarginal gyrus. Significant regression results (pFWE < 0.05) are displayed on the group average structural MRI in MNI space. Color bar indicates uncorrected t values for A–F.
Resting-state functional connectivity results: facet scores
For the FPN (Fig. 4, Table 4), individual facet scores were associated with distinct patterns of connectivity. Higher Facet 1 (interpersonal) scores were associated with significantly decreased connectivity between the R IPS seed and bilateral precuneus of the FPN, and the right medial frontal cortex/dorsal anterior cingulate of the CON. Similarly, higher Facet 3 (lifestyle) scores were associated with significantly decreased functional connectivity between the R IPS seed and right superior frontal gyrus within the CON. In contrast, higher Facet 4 (antisocial) scores were associated with significantly increased functional connectivity between the R IPS seed and prefrontal regions of the FPN (R IPS seed: left middle frontal gyrus, right anterior medial frontal gyrus) but significantly decreased functional connectivity with right mid-insula of the CON.
For the CON (Fig. 4, Table 4), higher Facet 1 (interpersonal) scores were associated with significantly decreased functional connectivity between the dACC/msFC seed and R IPL extending into the superior parietal lobule of the FPN. Similarly, higher Facet 2 (affective) scores predicted significantly decreased functional connectivity between the R aIfO seed and right postcentral gyrus/supramarginal gyrus within the CON. There were no significant associations between the other two Facet scores and functional connectivity for the CON.
Consistent with the other psychopathy scores, we found no significant associations between facet scores and comparison primary sensory network connectivity. Whole-brain results for clusters outside the networks of interest are reported in Table 5.
Whole-brain results for rest-state functional connectivity outside networks of interest
Discussion
In this study, we used resting-state fMRI in a large sample of prison inmates to examine whether psychopathy is linked to alterations in functional connectivity within and between the major cortical networks of the brain. Our results demonstrate multiple significant associations between psychopathic traits and cortical connectivity. Overall psychopathy severity (total PCL-R score) was associated principally with diminished connectivity between frontal and parietal regions, whereas the major subsets of psychopathic traits, including the two factors and four facets, exhibited distinct patterns of correlations. Factor 1 was linked exclusively to decreased cortical connectivity, whereas Factor 2 was linked exclusively to increased cortical connectivity. Similarly, Facets 1 and 2 were linked exclusively to decreased cortical connectivity, whereas Facet 4 was associated primarily with increased cortical connectivity. We did not find evidence for associations between connectivity and psychopathy scores for the comparison sensory networks. Here we consider each of the major findings in turn.
With respect to overall psychopathy severity, both the regression analysis and the between-group analysis revealed reduced functional connectivity between the right superior parietal cortex and the dACC, which are regions of the FPN and CON, respectively. The regression analysis included data from 142 inmates with a wide range of psychopathy severity (total PCL-R scores range, 5–36), whereas the between-groups analysis included 46 psychopathic and 49 nonpsychopathic offenders, respectively. This is the largest sample size ever reported in a study of resting-state functional connectivity of psychopathy. The convergent results between analyses indicate that diminished parietal–dACC connectivity is a robust neurobiological correlate of psychopathy. Although this functional connectivity result is novel, a number of previous studies have linked psychopathy with dACC dysfunction (Kiehl et al., 2001; Müller et al., 2003; Birbaumer et al., 2005; Koenigs, 2012; Ly et al., 2012; Meffert et al., 2013; Contreras-Rodríguez et al., 2014; Cope et al., 2014). For example, task-based fMRI studies have associated psychopathy with decreased dACC activity across a variety of affective paradigms (Kiehl et al., 2001; Müller et al., 2003; Birbaumer et al., 2005; Meffert et al., 2013; Cope et al., 2014). Furthermore, two resting-state fMRI studies have reported diminished dACC–insula connectivity in psychopathy (Ly et al., 2012; Contreras-Rodríguez et al., 2014). Although the functional significance of decreased parietal–dACC connectivity in psychopathy is unknown, some theories suggest the interaction between these regions, as portions of the FPN and CON or as part of the unitary frontoparietal control system, could be important for integrating salient emotional and cognitive information to facilitate goal-directed action (Vincent et al., 2008; Menon, 2011). Thus, in psychopathy, disrupted communication between regions of the FPN and CON (e.g., parietal–dACC connectivity) could lead to deficits in directing attention to salient stimuli. Behaviorally, this could account for impairments in psychopathy in the allocation of attention to salient external stimuli (e.g., threat-related cues; Newman et al., 1997, 2010; Hiatt et al., 2004; Glass and Newman, 2009). However, future neuroimaging and behavioral work will be crucial for understanding the relationship between parietal–dACC connectivity and psychopathic traits.
Although previous research has identified diminished DMN connectivity in psychopathy (Glenn et al., 2009; Motzkin et al., 2011), we did not observe this relationship for overall psychopathy severity in the present study. However, we did find distinct DMN connectivity associations with individual factor and facet scores—with diminished DMN connectivity related to Factor 1 and Facet 1 but enhanced connectivity related to Factor 2. Thus, one potential explanation of our null findings for overall psychopathy is that the distribution of factor/facet scores was different from the previous studies (e.g., greater number of participants with high Factor 1/Facet 1 scores; Glenn et al., 2009; Motzkin et al., 2011). Additional work in an independent sample is necessary to determine the precise relationships among overall psychopathy, psychopathy factors, and DMN connectivity.
With respect to individual factor and facet scores, our study demonstrates a striking neural dissociation for the major components of psychopathic personality. Factor 1, which includes the interpersonal and affective features of psychopathy (e.g., callous, egocentric, manipulative, and emotionally shallow) was linked to diminished connectivity in each of the three main cortical networks of interest (DMN, FPN, and CON). Similarly, Facet 1, the interpersonal traits also included in Factor 1, was associated with decreased connectivity in the same main networks. In contrast, Factor 2, which includes the lifestyle and antisocial features of psychopathy (e.g., impulsive, irresponsible, parasitic, and criminally versatile) was linked to heightened connectivity in the same networks. Facet 4, the antisocial traits also included in Factor 2, was associated with primarily increased connectivity in the FPN. Despite the fact that the two factor scores and four facet scores are generally positively correlated (factor scores, r = 0.68 in this inmate sample), these brain imaging findings suggest that the neurobiological correlates of each factor and facet are clearly separable. These results are also consistent with neuroimaging studies of psychopathy reporting distinct neural correlates for psychopathy factor and facet scores (Glenn et al., 2009; Contreras-Rodríguez et al., 2014). Specifically, Contreras-Rodríguez et al. (2014) found negative correlations between Factor 1 scores and gray matter volume but positive correlations between Factor 2 scores and gray matter volume for several brain regions, including the medial dorsal PFC, frontal operculum, and lateral PFC. During a moral decision-making fMRI task, Glenn et al. (2009) found an association between Facet 1 and decreased activity in the mPFC, PCC, and IPL of the DMN but no DMN associations with other facets. It is important to note that the dissociations observed here were not an artifact attributable to the partial correlations of factor scores (i.e., including Factor 2 score as a covariate for the Factor 1 analysis and vice versa). Even the zero-order correlations between functional connectivity values and each factor score (without covarying for the other) supports the same pattern of results obtained from the factor score regression analyses (Table 6). This partial correlation method of parsing the unique relationship between each factor score and a dependent measure of interest is standard practice in psychopathy research (Patrick et al., 1993; Patrick, 1994; Frick et al., 1999; Verona et al., 2001).
Zero-order correlations between factor scores and rest-state connectivity
Perhaps the most remarkable dissociation for the two factors was observed in the connectivity between the anterior insula and dACC regions of the CON. The CON, also referred to as the “salience network” (Seeley et al., 2007), may play an important role in several interrelated functions germane to psychopathy, including empathy (Engen and Singer, 2013), salience detection (Seeley et al., 2007), and cognitive control (Dosenbach et al., 2007; Menon and Uddin, 2010). Recent work has also suggested that a dACC–anterior insula network may function to coordinate relative levels of activity in the DMN (which has been linked to internally directed attention, such as mind wandering and self-reflection) and FPN (which has been linked to externally directed attention, such as cognitively demanding tasks and external stimulus processing; Sridharan et al., 2008; Goulden et al., 2014). Future work could examine whether Factor 1 and Factor 2 traits exhibit similar dissociations with regard to the flexible reallocation of attention between internal thoughts (Newman and Lorenz, 2003; Koenigs et al., 2010) and salient external stimuli (e.g., threat/emotion-related cues; Newman et al., 1997, 2010; Hiatt et al., 2004; Baskin-Sommers et al., 2009; Glass and Newman, 2009). In addition, task-based fMRI studies using emotion-related, cognitive, or social cognition tasks could be used to examine whether insula–dACC dysfunction is linked to altered DMN–FPN dynamics in psychopathy. More broadly, future studies could explore the relationship between cortical network connectivity (i.e., resting state) and relevant cognitive and emotional measures (e.g., executive function, self-related, and empathy) in relation to factors and facets of psychopathy.
We note that the observed relationships between factors scores and insula–dACC connectivity bear striking resemblance to the relationships between factor scores and anxiety. Like insula–dACC connectivity, anxiety is related negatively to Factor 1 scores but related positively to Factor 2 scores (Derefinko, 2014). Interestingly, CON connectivity/activity (in particular the anterior insula) is correlated positively with anxiety in healthy subjects (Seeley et al., 2007) and has been implicated in anxiety disorders (Menon, 2011). Although debate remains regarding the precise relationship between anxiety and psychopathy (in terms of both factor scores and psychopathy subtypes; Skeem et al., 2011), this combination of findings suggests that the delineation of the psychobiological mechanisms of psychopathic traits could benefit from future research on the link between insula–dACC connectivity, anxiety, and factor scores, as well as psychopathy subtypes. If these dissociable functional connectivity findings reflect distinct etiopathological mechanisms for Factor 1 and Factor 2 traits, then the development of more effective treatment and prevention strategies for psychopathy may benefit from targeting these traits separately. Furthermore, it is possible that neuroimaging measures could offer more accurate prognostic information than the standard interview assessment of psychopathic traits and behaviors (Aharoni et al., 2013).
In summary, here we have reported novel findings of functional brain abnormalities in psychopathic criminals (using both continuous and between-group comparison approaches). More specifically, we have demonstrated evidence of altered interactions within multiple cortical networks, including distinct network connectivity characteristics for individual psychopathy factors and facets. These results bolster an emerging literature on the neuropsychological basis of psychopathy.
Footnotes
This work was supported by National Institutes of Health Grants MH070539 (principle investigator, KAK), DA026505 (principle investigator, K.A.K.), MH078980 (principle investigator, J.N.), MH087525 (principle investigator, J. Decety), MH090169 (principle investigator, D.S. Kosson). We thank the many individuals at the Wisconsin Department of Corrections for making this research possible and are especially indebted to Deputy Warden Tom Nickel and Dr. Kevin Kallas. We also thank the Mind Research Network staff Keith Harenski, Dr. Nate Anderson, and Mike Maurer for their assistance with MRI data collection.
The authors declare no competing financial interests.
- Correspondence should be addressed to Carissa L. Philippi and Michael Koenigs, Department of Psychiatry, University of Wisconsin–Madison, 6001 Research Park Boulevard, Madison, WI 53719. cphilippi{at}wisc.edu; mrkoenigs{at}wisc.edu
