Early onset MJ use was associated with different patterns of cortical architecture.
•
Early vs. late onset divergence was in brain regions underlying higher-order cognition.
•
Findings were above and beyond effects of alcohol and current age.
Abstract
Background
As the most commonly used illicit substance during early adolescence, long-term or latent effects of early adolescent marijuana use across adolescent developmental processes remain to be determined.
Methods
We examined cortical thickness, gray/white matter border contrast (GWR) and local gyrification index (LGI) in 42 marijuana (MJ) users. Voxelwise regressions assessed early-onset (age <16) vs. late-onset (≥16 years-old) differences and relationships to continued use while controlling for current age and alcohol use.
Results
Although groups did not differ by onset status, groups diverged in their correlations between cannabis use and cortical architecture. Among early-onset users, continued years of MJ use and current MJ consumption were associated with thicker cortex, increased GWR and decreased LGI. Late-onset users exhibited the opposite pattern. This divergence was observed in all three morphological measures in the anterior dorsolateral frontal cortex (p < .05, FWE-corrected).
Conclusions
Divergent patterns between current MJ use and elements of cortical architecture were associated with early MJ use onset. Considering brain development in early adolescence, findings are consistent with disruptions in pruning. However, divergence with continued use for many years thereafter suggests altered trajectories of brain maturation during late adolescence and beyond.
With more than 25% of high school seniors reporting recent use and 6.5% of 12th graders being daily users (Johnston et al., 2014), marijuana (MJ) is the most frequently used illicit substance among adolescents. Across all age groups over 70% of new drug initiates start with using MJ at an average age of 18 years (SAMHSA, 2014). Indeed, the scope of MJ use prevalence is of great public interest, as MJ use in early adolescence is associated with increased risk of greater substance use, legal problems, disrupting education, injuries/medical problems, developing psychopathology, cognitive changes and chronic psychosocial struggles (CASA, 2011; Fergusson and Horwood, 1997; Fergusson et al., 1996; Patton et al., 2002). Taken together, rates of MJ use are suggestive of an epidemic based in adolescence, which is concerning not just due to societal cost, but also due to the potential to offset sensitive brain development during this period.
Adolescent MJ use has the potential to cause structural and functional changes in the brain by altering cannabinoid signaling. One possible mechanism would be blunt neurotoxic influence. For example, delta9-tetrahydrocannabinol (THC), the primary psychoactive component in MJ that binds CB1 receptors, is reported to cause cell shrinkage and damage DNA strands in THC-treated neuron cultures (Chan et al., 1998). This may be the mechanism by which smaller volumes have been observed in individuals exposed to cannabis during adolescence (Battistella et al., 2014). However, it is more likely that MJ exerts its influence on brain development indirectly. The cannabinoid system plays a role in modulating other neurotransmitters, including gamma-aminobutyric acid (GABA), glutamate and monoamines (Lopez-Moreno et al., 2008). Specifically, activation of CB1 receptors is associated with down-regulating inhibitory GABAergic transmission in cortical interneurons during adolescence (Caballero and Tseng, 2012; Cass et al., 2014). In addition, CB signaling inhibits microglia function (Walter et al., 2003). These two points are important because cortical pruning processes involve glial-mediated synaptic elimination and altering the excitatory/inhibitory balance is liable to disrupt the selective tagging and preserving synapses (Selemon, 2013). The impact of this indirect influence on the developing brain may be in the observations of abnormal connectivity in those who began MJ use in adolescence (Jacobus et al., 2009). Evidence from human neuroimaging studies lends greater support to MJ-related disruptions to brain development.
Structural neuroimaging studies have indicated that volumes of several brain areas are smaller in heavy adult MJ users especially in areas enriched with cannabinoid 1 (CB1) receptors, such as medial temporal lobe, and prefrontal cortex (Lorenzetti et al., 2010). Studies of adult chronic MJ users note brain volume reductions in temporal lobe, insula, and prefrontal cortex, amygdala and hippocampus (Battistella et al., 2014; Cousijn et al., 2012; Filbey et al., 2014; Matochik et al., 2005; Yucel et al., 2008). Among different characteristics of MJ involvement (e.g., dependence symptoms, use frequency, consumption), the age of initial MJ use is a robust factor that has been associated with smaller brain volumes in users. For example, Battistella et al. (2014) observed left parahippocampal gyrus and right temporal pole structural differences in 25 regular MJ users compared to 22 occasional users, however, even the occasional users who began smoking MJ during adolescence (before age 18) demonstrated similar brain changes as the regular users. Our group has also found links with early MJ use onset (Bava et al., 2009) and structural connectivity with orbitofrontal cortex in a cohort of daily MJ users, suggesting complex neuroadaptive processes related to MJ use in the context of adolescent brain development (Filbey et al., 2014). These findings underscore the potential for significant heterogeneity in brain changes among adult MJ users, especially those who began using MJ during neurodevelopment.
Studies comparing early adolescent MJ use to users initiating MJ use in later adolescence provide further evidence for the potential of MJ to cause enduring change. The few studies that have directly investigated the timing of the effects of MJ during adolescence have noted divergent neurodevelopment effects. For example, in an fMRI study by Gruber and colleagues, functional and behavioral differences during an interference task were reported between early (before age 16) and late (after age 16) MJ users (Gruber et al., 2012) (Sagar et al., 2015). The same group also reported decreased white matter integrity in early onset vs. late onset MJ users (mean age 14.46 vs. 17.93) (Gruber et al., 2014). Similar differential effects have also been noted in parietal lobe activation between early and late adolescent binge drinkers during a spatial working memory task (Tapert et al., 2004). These studies highlight the importance of clarifying the differential neural effects of early- and late-adolescent onset use.
To that end, in the current study, we compared daily MJ users who were early onset users (<16 years old) versus late onset users (≥16 years old) on measures of cortical morphology that are sensitive to developmental changes. We aimed to characterize both the effect of early onset status on cortical morphology as well as assess for morphological patterns linked to the continued use of MJ after early and late adolescent MJ initiation. We expected early onset users to show a morphological pattern consistent with disruption of early adolescent brain development (e.g., increased cortical thickness, greater gray/white definition of the cortical ribbon via disruptions to adolescent pruning processes) that may be more consistent with indirect impact of MJ of brain development. While gray matter decline has been shown to be associated with marijuana use, particularly in areas rich in CB1 receptors, increased cortical thickness and greater gray/white definition in the cortical ribbon point to potential disruption in neurodevelopment (i.e. synaptic pruning) that may result from MJ use at key developmental stages (i.e. earlier as opposed to later in adolescent neuronal development). Such disruptions may extend to gyrification as well. While this process begins in utero, there is evidence that gyrification is ongoing into adolescence (Armstrong et al., 1995; Alemán-Gómez et al., 2013; Klein et al., 2014) and may also display aberrant developmental patterns in the presence of MJ use.
2. Methods
This study was approved by the University of Texas at Dallas (UTD) and University of Texas Southwestern Medical Center (UTSW) Institutional Review Boards. All participants were recruited from the Dallas-Ft.Worth metro area via flyers and advertisements. Following informed consent, MJ users completed two sessions – a baseline appointment for collecting demographic, psychosocial and behavioral measures and a thorough substance use history. Three days later the participants returned for a neuroimaging appointment. Prior to their scanning session, participants were asked to be abstinent from MJ use for 72 h, from alcohol for 24 h, and from caffeine and cigarettes for the preceding 2 h. These were confirmed by self-report (MJ, alcohol, caffeine and cigarettes), quantitative THC urinalysis (MJ), and by breath alcohol level of .000 (alcohol) at the start of their session.
2.1. Participants
We scanned 45 regular heavy MJ users as part of the parent project. Inclusion criteria were: right-handedness, English as the primary language and no histories of psychosis, traumatic brain injury, and MRI contraindications (e.g., pregnancy, non-removal metallic implants, claustrophobia). One subject reported a history of anxiety and depression and one other reported a history of ADHD as a child. Additional exclusions for the current study included: Axis I diagnosis (via SCID) other than cannabis use disorder, unusable sMRI due to motion artifact or poor signal-to-noise ratio that precluded accurate tissue segmentation (n = 1) and incomplete drug use histories (n = 2). Of the 42 remaining cases, 22 were early onset users (onset of first use before age 16). Group categorization using onset of regular use as opposed to onset of first use maintained the same grouping (mean early onset of regular use = 16.5, mean late onset of regular use = 19.0). Regular use was defined as at least one time per week. To determine how age of onset of regular MJ use influenced our reported effects, we performed these analyses while covarying for age of onset of regular use (see Supplement). Table 1 summarizes demographic and substance use information according to onset status. Table 2 summarizes the correlation between age and identified marijuana use variables. Only MJ years of use and current age showed a statistically significant correlation. Participants were recruited based on self-reported daily MJ use and a positive urinalysis for THC metabolites at their baseline visit. All of the participants were screened via urinalysis for other drugs of abuse and were excluded if drugs (other than MJ) were detected. Participants were required to have used MJ for a minimum of 5000 lifetime occasions and self-report daily use (without >24 h abstinence) for the last 60 days.
2.2. MRI acquisition and analysis
2.2.1. Image acquisition
Scanning sessions took place at the Advanced Imaging Research Center at the University of Texas, Southwestern Medical Center three days following their initial visit. Another verification of THC metabolites via urinalysis was also performed before the scan. MRI images were collected using a 3T Philips whole-body scanner equipped with Quasar gradient subsystem (40 mT/m amplitude, a slew rate of 220 mT/m/ms). High-resolution T1-weighted anatomical scans were collected using a MPRAGE sequence: TR/TE/TI = 2100/3.70/1100 ms; flip angle = 12°; field of view = 256 mm × 256 mm; slab thickness = 160 mm (along left-right direction); voxel size = 1 mm × 1 mm × 1 mm, Total scan time = 3 m 57 s.
2.2.2. Image processing
MPRAGE anatomical scans were pre-processed for surface-based analyses using FreeSurfer v5.3 semi-automated pipeline (http://surfer.nmr.mgh.harvard.edu). This semi-automated pipeline included spatial (Talairach) and signal intensity normalization of images, volumetric segmentation and subcortical labeling (Dale et al., 1999; Fischl et al., 2002). Outer gray matter and white matter boundaries were then identified and reconstructed into a mesh of over 150,000 tessellated vertices to allow point-to-point surface measures (Fischl et al., 1999). Next, gyral anatomy is aligned to a standard spherical template using surface convexity and curvature measures. Resulting surfaces were inspected, blind to MJ onset status, to identify and correct any errors made during cortical reconstruction. Modifications to the volumes were made as necessary to correct for tissue misclassifications according to FreeSurfer's wiki manual (Schmansky et al., 2010). In preparation for analysis, each morphological measure for each case was co-registered to a standard template (fsaverage). Anatomical labels in FreeSurfer (Desikan et al., 2006) were used for interpretation of results.
2.3. Morphological measures
2.3.1. Cortical thickness
The width of the cortical ribbon was measured as the distance between corresponding vertices of the white matter and gray matter surfaces at each vertex in the cortical mantel (Fischl and Dale, 2000).
2.3.2. Gray–white matter ratio (GWR)
To assess the quality of cortical ribbon definition, a tissue contrast between gray and white matter signal intensities was computed as a percent ratio (W − G)/(.5*(W + G)) (from pctsurfcon v1.11.2.1, inbuilt component of FreeSurfer pipeline v5.3, 2011). White matter signal intensities were measured at an absolute length of 1 mm below the gray–white border surface and gray matter signal was measured 30% into the cortical ribbon (Salat et al., 2009).
2.3.3. Local gyrification index
The cortical surface from FreeSurfer's main pipeline is further processed to create an outer surface that encapsulates the gyral and sulcal curvature for each hemisphere, which serves as a basis for calculating a local gyrification index (Schaer et al., 2012). LGI is measured as the amount of cortex within the sulcal folds beneath the outer surface compared to the amount of visible cortex that touches the outer surface. Cortical maps are generated from repeated iterations of delineating a 25 mm radius sphere on the outer surface and its corresponding point on the cortical surface using a matching algorithm.
2.4. Background and premorbid characteristics
2.4.1. Sample characteristics
Age, gender, education level, ethnicity, along with other background information, was obtained using a standard demographics questionnaire. The two-subtest administration of the Wechsler Abbreviated Scale of Intelligence (Vocabulary and Matrix Reasoning) provided estimates of intellect (Wechsler, 1999).
2.4.2. Substance use
The Substance Use Disorder modules of the Structured Clinical Interview for DSM-IV (SCID) (First et al., 2002) were administered by a trained research assistant to assess for lifetime and current symptoms of abuse and dependence for alcohol, nicotine, MJ and other substances. The SCID interview also provided the onset of use information. A Time Line Follow-Back (TLFB) approach was used to quantify alcohol, nicotine, and MJ use patterns for 90 days prior to study participation (Sobell and Sobell, 1992). Marijuana use in grams was obtained via self-report in response to probes aimed at quantifying their regular use.
2.5. Statistical analyses
Statistical analyses were conducted in SPSS 18.0 for behavioral and psychosocial measures whereas general linear model group comparisons on surfaced-based morphology measures were carried out FreeSurfer's built-in application QDEC (v1.5). Independent samples t-tests, Mann–Whitney U-tests or chi-square tests, compared groups on background and demographic variables (see Table 1). Before statistical analysis was conducted, the dependent measures of cortical thickness, GWR and LGI were smoothed using a FWHM Gaussian filter with a width of either 10 or 15 mm. Separate univariate general linear model (GLM) was then used to model cortical thickness, GWR and LGI with onset status of MJ use as a between groups factor. The dependent variables were thickness, gray–white ratio or local gyrification index and the independent variables were either recent monthly MJ use in grams (MJ grams) or duration of MJ use (MJ years). Age and total drinks in the past 2 months were treated as nuisance covariates in the model. Using MJ years of use and MJ grams as independent predictors of interest allowed us to characterize and differentiate the latent developmental effects from cumulative and current effects of MJ use. The variable “marijuana years of use” was based on the participants’ response to the question “For how many years have you been using marijuana regularly?” Of note, an outlier in the early onset group was removed before the statistical comparisons were performed.
3. Results
3.1. Cortical thickness
There were no regions of group differences in cortical thickness by early onset status alone, controlled for age and alcohol use. However, MJ use characteristics were correlated with anterior dorsolateral prefrontal cortex thickness based on onset status. Early onset users showed increased thickness with increased MJ grams while late onset users showed thinner cortex with increased MJ grams (p < 0.05 uncorrected) (Table 3). The same pattern emerged with more years of MJ use being associated with thicker region of the right medial temporal lobe in the early onset users and the reverse for the late onset users (p < 0.05 uncorrected) (Fig. 1).
Table 3. Clusters of significant age of onset × marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.
The effect sizes were derived from Freesurfer's tool in the significant region of interest using mri_segstats. This was also confirmed manually by using the F-value reported by Freesurfer.
3.2. Gray–white matter contrast
There were no regions of group differences in gray–white matter contrast by early onset status alone, controlled for age and alcohol use. However, current MJ consumption (grams) and onset status were differentially correlated with gray–white matter contrast in a left anterior dorsal frontal region (p < 0.05, FWE corrected). Increased gray–white contrast with heavier MJ use was seen in the early onset users and the opposite was seen in later onset users (heavier current use linked to decreasing GWR). The same pattern was seen between duration of MJ use in two prefrontal cortex clusters of the right dorsal frontal and medial orbitofrontal area p < 0.05, FWE corrected – more years of MJ use were linked to greater GWR among early users (Fig. 1).
3.3. Gyrification
MJ use onset status alone showed no significant main effects above age and alcohol covariates. However, onset status was correlated with divergent patterns between local gyrification and MJ use, whereby early onset users showed decreasing LGI with increasing MJ consumption and longer duration of use in prefrontal cortex regions p < 0.05, FWE corrected. The left hemisphere clusters encompassed the majority of the length of the middle lateral surface of the left cortex, including motor cortices, parietal lobe and multimodal integration areas (Fig. 1).
4. Discussion
The present study was designed to characterize the cortical architecture in adolescent onset MJ users by comparing early adolescent onset users to late adolescent onset in MJ use on measures of cortical thickness, gray/white matter contrast and gyrification. The primary finding was that early versus late onset MJ users showed a divergent pattern in cortical thickness, definition of the cortical ribbon and local gyrification with continued use through and beyond adolescent years. Specifically, early onset users showed cortical thickening, enhanced gray/white matter contrast, and decreased gyrification in association with more years of MJ use and current consumption of MJ in grams in frontal and temporal regions – areas that underlie higher order cognition including executive functioning, learning and memory. Findings were above and beyond effects of alcohol and current age, therefore, results are less likely to reflect morphological trends due to aging.
Our findings did not find the expected age of onset differences previously reported in marijuana users (Gruber et al., 2012, 2014). This inconsistency suggests that the age of onset effects may be more robust in brain white matter connectivity (Gruber et al., 2014) and function (Gruber et al., 2012) than brain surface morphometry. To date, the few studies that have described altered cortical morphology in MJ users have led to mixed findings. Mata et al. (2010) identified brain regions with decreased sulcal depth suggestive of lower gyrification in a study of adult MJ users. Jacobus and Tapert (2014) recently reported increased cortical thickness in the entorhinal cortex among 24 adolescent MJ users (mean age = 17.7, mean MJ onset age = 15.4) relative to peer controls. However, the authors also reported a negative relationship between cortical thickness and total MJ use in the right paracentral gyrus, and they observed consistent positive relationships in various brain regions between age of MJ onset and thickness. In the only other known adolescent study of cortical thickness and MJ, Lopez-Larson and colleagues studied 18 adolescent heavy MJ users (similar in age and MJ onset as Jacobus and Tapert, 2014) and reported mixed findings of increased thickness in prefrontal/insula regions and decreased thickness in posterior/temporal lobe areas in the MJ users compared to controls. In contrast to Jacobus and Tapert (2014), Lopez-Larson et al. (2011) found areas of the frontal lobe and insula that were thinner with increased urine THC metabolites and thicker with earlier age of onset. Select findings from the current study align with aspects of both of these studies, with a consensus supporting findings of a negative dose-dependent relationship between MJ use and cortical thickness. Given the low availability of studies to compare, this consensus is very limited. Although Jacobus et al. and Lopez-Larson et al. found the opposite effect of age of onset on thickness, the pattern of divergence among early vs. late onset users in the current study is more consistent with the latter study, whereby we saw early onset users exhibit thicker cortex with continued MJ use. Taken together, findings of increased thickness related to early MJ onset accompanied by negative dose-dependent relationships with MJ exposure may reflect two distinct processes. One process may be specific to the interactions with cortical development during early adolescence, likely leading to a disruption in pruning, and, the other, specific to the pharmacological effect with heavy chronic MJ use.
In the only known study to examine the curvature-morphology of the cortex in adult MJ users, Mata et al. (2010) identified decreased sulcal concavity and thinner sulci in 23 MJ users compared to controls (n = 44), also in prefrontal areas. However, they did not observe significant relationships with age, MJ onset age, or cumulative MJ use. It is interesting that the authors detected group level differences (MJ vs. controls) but no correlations with MJ use characteristics such as dose or age of onset, whereas our primary findings are the consistent effects of continued MJ use differing after early or late adolescent onset. There are substantial methodological explanations for this disparity. For example, the current study did not compare morphology in MJ users to a normative control sample, therefore, it is feasible that group-level differences may emerge with such a comparison. Likewise, we deliberately covaried for current age in order to control for brain changes with aging and thus optimize our interrogation of developmental effects of early onset age and of aspects of continued use.
The heterogeneity of MJ effects clearly suggests a multifactorial system of neurobiological processes involved. The primary results uphold that age of onset is a robust variable that differentiates heavy MJ users based on early versus late MJ onset. However, this group distinction relied on current use characteristics. Therefore, in the absence of group-level differences, the interactions between onset age and current use indicates that continued cannabis exposure and early adolescent developmental factors both contribute to a dynamic and sustained departure from what is expected based on developmental studies.
Typical synaptic refinement processes during early adolescence are in the context of long-term depression and potentiation of cortical neurons in order to facilitate neuronal remodeling. Thus, the normal course of early adolescent development is uniquely vulnerable to disruption by MJ due to the electrochemical conditions and maturity of brain processes that would not present together again. Cass and colleagues tested the sensitivity of early adolescence cannabinoid exposure in an animal model (Cass et al., 2014). They found that acute administration of cannabinoid agonists in early, middle and late adolescent rats led to a state of frequency-dependent disinhibition of neurons in the frontal cortex in the early-to-middle adolescent rats, but not in the late adolescent rats. Moreover, the authors also noted that adult rats previously exposed to cannabinoid agonists in adolescence displayed comparable neuronal disinhibition. Thus, by changing the inhibitory/excitatory landscape during adolescence, MJ can influence lasting changes to typical cortical remodeling during sensitive early adolescent years.
The sequence of pruning and myelination likely plays a formative role in lasting changes from early adolescent onset MJ use. With decreased synaptic elimination, our findings of greater GW border contrast may reflect greater proliferation of myelin at the boundary of the cortical ribbon where non-pruned synapses remained with linked axons. Findings of altered white matter tissue qualities are mixed in adolescent and adult MJ user samples. Some report both increases and decreases in fractional anisotropy (FA) and average water diffusion (Bava et al., 2009) whereas others report consistent decreases in FA among adolescent MJ users (Ashtari et al., 2009; Jacobus et al., 2009) or null findings (Delisi et al., 2006). Two studies of diffusion tensor imaging in adult MJ users reported reduced FA in users compared to controls (Gruber et al., 2011, 2014). In addition to equivocal findings, research is needed to address the microstructural changes that could result in altered definition of the cortical ribbon. For example, rather than whole brain techniques that assess diffusion measures along major white matter tracts, indices assessing axonal organization along radial and interneuronassociation fibers along the cortical ribbon are needed. This scenario played out could result in increased gray matter (thicker cortex from disrupted pruning) and the myelination of connections to these spared terminals would result in increased density of white matter at the cortical boundary. Without any known studies of adolescent development of the gray/white tissue contrast at the cortical border to serve as a point of comparison, we speculate that early adolescent disruption of pruning and subsequent myelination of connections at the cortical boundary would be reflected by increased GWR as we saw in the current study.
5. Limitations and conclusions
The cross-sectional nature of this study limits causal attributions in terms of what we can infer to be directly related to the effects of MJ. Although a longitudinal design is optimal for addressing brain changes directly due to MJ, cross-sectional studies facilitate data-driven hypotheses that can be assessed directly in prospective studies.
It is important to keep in mind that the participants were not explicitly asked for possible years of abstinence during their period of regular use, which may have created possible inflation in reported duration of regular use. However, because the participants provided number of years of “regular” marijuana use, this inherently suggests continued, uninterrupted years of use. Concurrent nicotine use could have also influenced our reported results. But in the absence of a larger sample size and the presence of huge variance in nicotine use in the current sample, we were unable to verify the effect of nicotine use in the reported results.
Interpretation of these findings is also limited by the lack of behavioral anchors for the observed morphological effects and lack of information on other aspects of developmental history that could further characterize the effects of marijuana during neurodevelopment. This is further limited by the absence of “expected” patterns based on normative data. Given the varied directions of effects and the small sample size, these findings should be replicated and be viewed as preliminary.
To conclude, early MJ use was linked to altered neurodevelopmental patterns in brain regions sub-serving higher-order cognitive process. Clinical implications include need for early, targeted intervention. Given that the most robust results were related to interactions between onset age and continued use through emerging adulthood, harm reduction approaches may be effective in moderating adolescent MJ use to levels that are less likely to cause long-term developmental changes.
Conflict of interest
The authors report no conflicts of interest.
Acknowledgements
This research was funded by the National Institute on Drug Abuse (R01 DA030344, Filbey). We would like to thank all the participants who volunteered for this study. We are also very grateful to Talha Alvi, Sina Aslan, Jessica Baine, Collette Bice, Vicki Germer, Ariel Ketcherside, Alison King, Brittany Kuhn, Tyler Rhinehardt, Wing Ting To and the team of lab interns for their assistance with recruitment, running participants and data management.
Appendix A. Supplementary data
The following are Supplementary data to this article:
M. Ashtari, K. Cervellione, J. Cottone, B.A. Ardekani, S. Sevy, S. KumraDiffusion abnormalities in adolescents and young adults with a history of heavy cannabis use
D.K. Cass, E. Flores-Barrera, D.R. Thomases, W.F. Vital, A. Caballero, K.Y. TsengCB1 cannabinoid receptor stimulation during adolescence impairs the maturation of GABA function in the adult rat prefrontal cortex
J. Cousijn, R.W. Wiers, K.R. Ridderinkhof, W. van den Brink, D.J. Veltman, A.E. GoudriaanGrey matter alterations associated with cannabis use: results of a VBM study in heavy cannabis users and healthy controls
L.E. Delisi, H.C. Bertisch, K.U. Szulc, M. Majcher, K. Brown, A. Bappal, B.A. ArdekaniA preliminary DTI study showing no brain structural change associated with adolescent cannabis use
R.S. Desikan, F. Segonne, B. Fischl, B.T. Quinn, B.C. Dickerson, D. Blacker, ..., R.J. KillianyAn automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest
M.B. First, R.L. Spitzer, G. Miriam, J.B.W. WilliamsStructured Clinical Interview for DSM-IV-TR Axis I Disorders, Research Version, Non-patient Edition (SCID-I/NP)
Biometrics Research, New York State Psychiatric Institute, New York (2002)
B. Fischl, D.H. Salat, E. Busa, M. Albert, M. Dieterich, C. Haselgrove, ..., A.M. DaleWhole brain segmentation: automated labeling of neuroanatomical structures in the human brain
A. Giorgio, K.E. Watkins, M. Chadwick, S. James, L. Winmill, G. Douaud, ..., A.C. JamesLongitudinal changes in grey and white matter during adolescence
N. Gogtay, J.N. Giedd, L. Lusk, K.M. Hayashi, D. Greenstein, ..., P.M. ThompsonDynamic mapping of human cortical development during childhood through early adulthood
Proc. Natl. Acad. Sci. U.S.A., 101 (21) (2004), pp. 8174-8179
S.A. Gruber, M.M. Silveri, M.K. Dahlgren, D. Yurgelun-ToddWhy so impulsive? White matter alterations are associated with impulsivity in chronic marijuana smokers
Exp. Clin. Psychopharmacol., 19 (3) (2011), pp. 231-242, 10.1037/a0023034
K.M. Hasan, A. Sankar, C. Halphen, L.A. Kramer, M.E. Brandt, J. Juranek, ..., L. Ewing-CobbsDevelopment and organization of the human brain tissue compartments across the lifespan using diffusion tensor imaging
J. Jacobus, T. McQueeny, S. Bava, B.C. Schweinsburg, L.R. Frank, T.T. Yang, S.F. TapertWhite matter integrity in adolescents with histories of marijuana use and binge drinking
J.A. Lopez-Moreno, G. Gonzalez-Cuevas, G. Moreno, M. NavarroThe pharmacology of the endocannabinoid system: functional and structural interactions with other neurotransmitter systems and their repercussions in behavioral addiction
D.H. Salat, S.Y. Lee, A.J. van der Kouwe, D.N. Greve, B. Fischl, H.D. RosasAge-associated alterations in cortical gray and white matter signal intensity and gray to white matter contrast
M. Schaer, M.B. Cuadra, N. Schmansky, B. Fischl, J.P. Thiran, S. EliezHow to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index
S.F. Tapert, A.D. Schweinsburg, V.C. Barlett, S.A. Brown, L.R. Frank, G.G. Brown, M.J. MeloyBlood oxygen level dependent response and spatial working memory in adolescents with alcohol use disorders
Alcohol Clin. Exp. Res., 28 (10) (2004), pp. 1577-1586
M. Verdurand, V. Nguyen, D. Stark, D. Zahra, M.C. Gregoire, I. Greguric, K. ZavitsanouComparison of cannabinoid CB(1) receptor binding in adolescent and adult rats: a positron emission tomography study using [F]MK-9470
M. Yucel, N. Solowij, C. Respondek, S. Whittle, A. Fornito, C. Pantelis, D.I. LubmanRegional brain abnormalities associated with long-term heavy cannabis use
Early onset MJ use was associated with different patterns of cortical architecture.
+
+
+
+
Early vs. late onset divergence was in brain regions underlying higher-order cognition.
+
+
+
+
Findings were above and beyond effects of alcohol and current age.
+
+
+
+
+
+
+ Background
+
As the most commonly used illicit substance during early adolescence, long-term or latent effects of early adolescent marijuana use across adolescent developmental processes remain to be determined.
+
+
+ Methods
+
We examined cortical thickness, gray/white matter border contrast (GWR) and local gyrification index (LGI) in 42 marijuana (MJ) users. Voxelwise regressions assessed early-onset (age <16) vs. late-onset (≥16 years-old) differences and relationships to continued use while controlling for current age and alcohol use.
+
+
+ Results
+
Although groups did not differ by onset status, groups diverged in their correlations between cannabis use and cortical architecture. Among early-onset users, continued years of MJ use and current MJ consumption were associated with thicker cortex, increased GWR and decreased LGI. Late-onset users exhibited the opposite pattern. This divergence was observed in all three morphological measures in the anterior dorsolateral frontal cortex (p < .05, FWE-corrected).
+
+
+ Conclusions
+
Divergent patterns between current MJ use and elements of cortical architecture were associated with early MJ use onset. Considering brain development in early adolescence, findings are consistent with disruptions in pruning. However, divergence with continued use for many years thereafter suggests altered trajectories of brain maturation during late adolescence and beyond.
With more than 25% of high school seniors reporting recent use and 6.5% of 12th graders being daily users (Johnston et al., 2014), marijuana (MJ) is the most frequently used illicit substance among adolescents. Across all age groups over 70% of new drug initiates start with using MJ at an average age of 18 years (SAMHSA, 2014). Indeed, the scope of MJ use prevalence is of great public interest, as MJ use in early adolescence is associated with increased risk of greater substance use, legal problems, disrupting education, injuries/medical problems, developing psychopathology, cognitive changes and chronic psychosocial struggles (CASA, 2011, Fergusson and Horwood, 1997, Fergusson et al., 1996, Patton et al., 2002). Taken together, rates of MJ use are suggestive of an epidemic based in adolescence, which is concerning not just due to societal cost, but also due to the potential to offset sensitive brain development during this period.
+
Despite its prevalence, the impact of MJ use on adolescent brain development is not fully known. Important neuromaturational processes during adolescence through young adulthood are believed to bring about improved higher-order cognition by refining neural systems locally and globally through white and gray matter development (Casey et al., 2005, Giedd, 2008, Paus, 2005). In general, gray matter reductions and cortical thinning coincide with increased white matter volume and organization through adolescence and young adulthood, suggestive of synaptic pruning and axonal myelination (Giorgio et al., 2010, Gogtay et al., 2004, Hasan et al., 2007, Lebel et al., 2010, Shaw et al., 2008). The endogenous cannabinoid (CB) system is also immature during adolescence (Anavi-Goffer and Mulder, 2009, Verdurand et al., 2011). In an animal model (Verdurand et al., 2011) imaged CB1 receptor binding using PET and found relatively lower activation of CB1 receptors in adolescent rats compared to adult rats in brain areas including those in the frontal cortex, temporal lobe (hippocampus and amygdala) and sub-cortical regions including striatal regions, thalamus, hypothalamus, superior colliculus. Thus, adolescence represents a developmental period with vulnerability to structural and functional changes due to exogenous MJ exposure.
+
Adolescent MJ use has the potential to cause structural and functional changes in the brain by altering cannabinoid signaling. One possible mechanism would be blunt neurotoxic influence. For example, delta9-tetrahydrocannabinol (THC), the primary psychoactive component in MJ that binds CB1 receptors, is reported to cause cell shrinkage and damage DNA strands in THC-treated neuron cultures (Chan et al., 1998). This may be the mechanism by which smaller volumes have been observed in individuals exposed to cannabis during adolescence (Battistella et al., 2014). However, it is more likely that MJ exerts its influence on brain development indirectly. The cannabinoid system plays a role in modulating other neurotransmitters, including gamma-aminobutyric acid (GABA), glutamate and monoamines (Lopez-Moreno et al., 2008). Specifically, activation of CB1 receptors is associated with down-regulating inhibitory GABAergic transmission in cortical interneurons during adolescence (Caballero and Tseng, 2012, Cass et al., 2014). In addition, CB signaling inhibits microglia function (Walter et al., 2003). These two points are important because cortical pruning processes involve glial-mediated synaptic elimination and altering the excitatory/inhibitory balance is liable to disrupt the selective tagging and preserving synapses (Selemon, 2013). The impact of this indirect influence on the developing brain may be in the observations of abnormal connectivity in those who began MJ use in adolescence (Jacobus et al., 2009). Evidence from human neuroimaging studies lends greater support to MJ-related disruptions to brain development.
+
Structural neuroimaging studies have indicated that volumes of several brain areas are smaller in heavy adult MJ users especially in areas enriched with cannabinoid 1 (CB1) receptors, such as medial temporal lobe, and prefrontal cortex (Lorenzetti et al., 2010). Studies of adult chronic MJ users note brain volume reductions in temporal lobe, insula, and prefrontal cortex, amygdala and hippocampus (Battistella et al., 2014, Cousijn et al., 2012, Filbey et al., 2014, Matochik et al., 2005, Yucel et al., 2008). Among different characteristics of MJ involvement (e.g., dependence symptoms, use frequency, consumption), the age of initial MJ use is a robust factor that has been associated with smaller brain volumes in users. For example, Battistella et al. (2014) observed left parahippocampal gyrus and right temporal pole structural differences in 25 regular MJ users compared to 22 occasional users, however, even the occasional users who began smoking MJ during adolescence (before age 18) demonstrated similar brain changes as the regular users. Our group has also found links with early MJ use onset (Bava et al., 2009) and structural connectivity with orbitofrontal cortex in a cohort of daily MJ users, suggesting complex neuroadaptive processes related to MJ use in the context of adolescent brain development (Filbey et al., 2014). These findings underscore the potential for significant heterogeneity in brain changes among adult MJ users, especially those who began using MJ during neurodevelopment.
+
Studies comparing early adolescent MJ use to users initiating MJ use in later adolescence provide further evidence for the potential of MJ to cause enduring change. The few studies that have directly investigated the timing of the effects of MJ during adolescence have noted divergent neurodevelopment effects. For example, in an fMRI study by Gruber and colleagues, functional and behavioral differences during an interference task were reported between early (before age 16) and late (after age 16) MJ users (Gruber et al., 2012) (Sagar et al., 2015). The same group also reported decreased white matter integrity in early onset vs. late onset MJ users (mean age 14.46 vs. 17.93) (Gruber et al., 2014). Similar differential effects have also been noted in parietal lobe activation between early and late adolescent binge drinkers during a spatial working memory task (Tapert et al., 2004). These studies highlight the importance of clarifying the differential neural effects of early- and late-adolescent onset use.
+
To that end, in the current study, we compared daily MJ users who were early onset users (<16 years old) versus late onset users (≥16 years old) on measures of cortical morphology that are sensitive to developmental changes. We aimed to characterize both the effect of early onset status on cortical morphology as well as assess for morphological patterns linked to the continued use of MJ after early and late adolescent MJ initiation. We expected early onset users to show a morphological pattern consistent with disruption of early adolescent brain development (e.g., increased cortical thickness, greater gray/white definition of the cortical ribbon via disruptions to adolescent pruning processes) that may be more consistent with indirect impact of MJ of brain development. While gray matter decline has been shown to be associated with marijuana use, particularly in areas rich in CB1 receptors, increased cortical thickness and greater gray/white definition in the cortical ribbon point to potential disruption in neurodevelopment (i.e. synaptic pruning) that may result from MJ use at key developmental stages (i.e. earlier as opposed to later in adolescent neuronal development). Such disruptions may extend to gyrification as well. While this process begins in utero, there is evidence that gyrification is ongoing into adolescence (Armstrong et al., 1995, Alemán-Gómez et al., 2013, Klein et al., 2014) and may also display aberrant developmental patterns in the presence of MJ use.
+
+
+
+ Methods
+
This study was approved by the University of Texas at Dallas (UTD) and University of Texas Southwestern Medical Center (UTSW) Institutional Review Boards. All participants were recruited from the Dallas-Ft.Worth metro area via flyers and advertisements. Following informed consent, MJ users completed two sessions – a baseline appointment for collecting demographic, psychosocial and behavioral measures and a thorough substance use history. Three days later the participants returned for a neuroimaging appointment. Prior to their scanning session, participants were asked to be abstinent from MJ use for 72 h, from alcohol for 24 h, and from caffeine and cigarettes for the preceding 2 h. These were confirmed by self-report (MJ, alcohol, caffeine and cigarettes), quantitative THC urinalysis (MJ), and by breath alcohol level of .000 (alcohol) at the start of their session.
+
+
+ Participants
+
We scanned 45 regular heavy MJ users as part of the parent project. Inclusion criteria were: right-handedness, English as the primary language and no histories of psychosis, traumatic brain injury, and MRI contraindications (e.g., pregnancy, non-removal metallic implants, claustrophobia). One subject reported a history of anxiety and depression and one other reported a history of ADHD as a child. Additional exclusions for the current study included: Axis I diagnosis (via SCID) other than cannabis use disorder, unusable sMRI due to motion artifact or poor signal-to-noise ratio that precluded accurate tissue segmentation (n = 1) and incomplete drug use histories (n = 2). Of the 42 remaining cases, 22 were early onset users (onset of first use before age 16). Group categorization using onset of regular use as opposed to onset of first use maintained the same grouping (mean early onset of regular use = 16.5, mean late onset of regular use = 19.0). Regular use was defined as at least one time per week. To determine how age of onset of regular MJ use influenced our reported effects, we performed these analyses while covarying for age of onset of regular use (see Supplement). Table 1 summarizes demographic and substance use information according to onset status. Table 2 summarizes the correlation between age and identified marijuana use variables. Only MJ years of use and current age showed a statistically significant correlation. Participants were recruited based on self-reported daily MJ use and a positive urinalysis for THC metabolites at their baseline visit. All of the participants were screened via urinalysis for other drugs of abuse and were excluded if drugs (other than MJ) were detected. Participants were required to have used MJ for a minimum of 5000 lifetime occasions and self-report daily use (without >24 h abstinence) for the last 60 days.
Sample characteristics. MJ, marijuana.
Measure
Early onset (n = 20)
Late onset (n = 22)
p-Value
Effect size***
Statistic
Mean
(SD)
Min–Max
Mean
(SD)
Min–Max
|t/U
Age
32.50
(8.01)
21–50
30.25
(7.19)
21–47
0.316
0.302
t = 1.01
Education (years)
12.91
(2.54)
8–18
13.26
(2.40)
10–19
0.651
0.144
t = 0.456
Gender (male)
55%
73%
0.241
0.034
χ2 = 1.41
Ethnicity (% Caucasian)
50%
50%
0.566
0.008
χ2 = 0.336
IQ*
108
(9.99)
88–124
105
(13.54)
83–129
0.351
0.298
t = 0.94
Age of first MJ use**
13.18
(1.89)
9–15
16.90
(1.48)
16–21
<0.001
0.866
U = 0
Age of regular MJ use**
16.50
(3.57)
9–25
19.00
(4.29)
16–36
0.004
0.439
U = 108
Substance use in the last 60 days
MJ grams (daily)
2.14
(1.79)
0.50–7.50
1.65
(1.21)
0.46–4.23
0.338
0.083
U = 182
# EtOH drinks
44.09
(76.30)
0–310
38.85
(58.61)
0–183
0.588
0.039
U = 198.05
Max # EtOH drinks
6.62
(7.21)
0–31
6.25
(5.80)
0–21
0.8
0.095
U = 210
# EtOH drinking days
11.09
(16.69)
0–59
9.84
(13.51)
0–60
0.537
0.006
U = 185.5
# Binge EtOH drinking days
4.36
(12.50)
0–59
2.90
(5.53)
0–19
0.968
0.035
U = 218.5
# EtOH drinks per day
2.99
(2.27)
0–7.40
3.56
(3.29)
0–14.00
0.82
0.289
U = 211
# Cigarette days
1.18
(3.72)
0–17
2.95
(1.17)
0–21
0.06
0.294
U = 159
# Cigarettes per day
0.22
(0.55)
0–2.00
0.78
(1.23)
0–4.50
0.057
0.296
U = 158
Max # cigarettes
0.25
(0.61)
0–2
0.96
(1.60)
0–6
0.054
0.290
U = 157.5
Illicit drug use/past 90 days
14%
5%
Lifetime illicit drug use
73%
75%
IQ scores derived from Wechsler Abbreviated Scale of Intelligence Vocabulary and Matrix Reasoning subtests.
p < .05; SS, standard score; |t|, absolute value of student's t, U is the Mann–Whitney U's score.
The effect sizes of the above table were calculated either based on mean differences if normally distributed, correlation coefficient or F-value score using the default Cohen's effect size formula for respective metrics.
The correlations between current age and all MJ use variables.
Scanning sessions took place at the Advanced Imaging Research Center at the University of Texas, Southwestern Medical Center three days following their initial visit. Another verification of THC metabolites via urinalysis was also performed before the scan. MRI images were collected using a 3T Philips whole-body scanner equipped with Quasar gradient subsystem (40 mT/m amplitude, a slew rate of 220 mT/m/ms). High-resolution T1-weighted anatomical scans were collected using a MPRAGE sequence: TR/TE/TI = 2100/3.70/1100 ms; flip angle = 12°; field of view = 256 mm × 256 mm; slab thickness = 160 mm (along left-right direction); voxel size = 1 mm × 1 mm × 1 mm, Total scan time = 3 m 57 s.
+
+
+
+ Image processing
+
MPRAGE anatomical scans were pre-processed for surface-based analyses using FreeSurfer v5.3 semi-automated pipeline (http://surfer.nmr.mgh.harvard.edu). This semi-automated pipeline included spatial (Talairach) and signal intensity normalization of images, volumetric segmentation and subcortical labeling (Dale et al., 1999, Fischl et al., 2002). Outer gray matter and white matter boundaries were then identified and reconstructed into a mesh of over 150,000 tessellated vertices to allow point-to-point surface measures (Fischl et al., 1999). Next, gyral anatomy is aligned to a standard spherical template using surface convexity and curvature measures. Resulting surfaces were inspected, blind to MJ onset status, to identify and correct any errors made during cortical reconstruction. Modifications to the volumes were made as necessary to correct for tissue misclassifications according to FreeSurfer's wiki manual (Schmansky et al., 2010). In preparation for analysis, each morphological measure for each case was co-registered to a standard template (fsaverage). Anatomical labels in FreeSurfer (Desikan et al., 2006) were used for interpretation of results.
The width of the cortical ribbon was measured as the distance between corresponding vertices of the white matter and gray matter surfaces at each vertex in the cortical mantel (Fischl and Dale, 2000).
+
+
+
+ Gray–white matter ratio (GWR)
+
To assess the quality of cortical ribbon definition, a tissue contrast between gray and white matter signal intensities was computed as a percent ratio (W − G)/(.5*(W + G)) (from pctsurfcon v1.11.2.1, inbuilt component of FreeSurfer pipeline v5.3, 2011). White matter signal intensities were measured at an absolute length of 1 mm below the gray–white border surface and gray matter signal was measured 30% into the cortical ribbon (Salat et al., 2009).
+
+
+
+ Local gyrification index
+
The cortical surface from FreeSurfer's main pipeline is further processed to create an outer surface that encapsulates the gyral and sulcal curvature for each hemisphere, which serves as a basis for calculating a local gyrification index (Schaer et al., 2012). LGI is measured as the amount of cortex within the sulcal folds beneath the outer surface compared to the amount of visible cortex that touches the outer surface. Cortical maps are generated from repeated iterations of delineating a 25 mm radius sphere on the outer surface and its corresponding point on the cortical surface using a matching algorithm.
Age, gender, education level, ethnicity, along with other background information, was obtained using a standard demographics questionnaire. The two-subtest administration of the Wechsler Abbreviated Scale of Intelligence (Vocabulary and Matrix Reasoning) provided estimates of intellect (Wechsler, 1999).
+
+
+
+ Substance use
+
The Substance Use Disorder modules of the Structured Clinical Interview for DSM-IV (SCID) (First et al., 2002) were administered by a trained research assistant to assess for lifetime and current symptoms of abuse and dependence for alcohol, nicotine, MJ and other substances. The SCID interview also provided the onset of use information. A Time Line Follow-Back (TLFB) approach was used to quantify alcohol, nicotine, and MJ use patterns for 90 days prior to study participation (Sobell and Sobell, 1992). Marijuana use in grams was obtained via self-report in response to probes aimed at quantifying their regular use.
+
+
+
+
+ Statistical analyses
+
Statistical analyses were conducted in SPSS 18.0 for behavioral and psychosocial measures whereas general linear model group comparisons on surfaced-based morphology measures were carried out FreeSurfer's built-in application QDEC (v1.5). Independent samples t-tests, Mann–Whitney U-tests or chi-square tests, compared groups on background and demographic variables (see Table 1). Before statistical analysis was conducted, the dependent measures of cortical thickness, GWR and LGI were smoothed using a FWHM Gaussian filter with a width of either 10 or 15 mm. Separate univariate general linear model (GLM) was then used to model cortical thickness, GWR and LGI with onset status of MJ use as a between groups factor. The dependent variables were thickness, gray–white ratio or local gyrification index and the independent variables were either recent monthly MJ use in grams (MJ grams) or duration of MJ use (MJ years). Age and total drinks in the past 2 months were treated as nuisance covariates in the model. Using MJ years of use and MJ grams as independent predictors of interest allowed us to characterize and differentiate the latent developmental effects from cumulative and current effects of MJ use. The variable “marijuana years of use” was based on the participants’ response to the question “For how many years have you been using marijuana regularly?” Of note, an outlier in the early onset group was removed before the statistical comparisons were performed.
+
+
+
+
+ Results
+
+
+ Cortical thickness
+
There were no regions of group differences in cortical thickness by early onset status alone, controlled for age and alcohol use. However, MJ use characteristics were correlated with anterior dorsolateral prefrontal cortex thickness based on onset status. Early onset users showed increased thickness with increased MJ grams while late onset users showed thinner cortex with increased MJ grams (p < 0.05 uncorrected) (Table 3). The same pattern emerged with more years of MJ use being associated with thicker region of the right medial temporal lobe in the early onset users and the reverse for the late onset users (p < 0.05 uncorrected) (Fig. 1).
Clusters of significant age of onset × marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.
Measure
Label@Max, Extended coverage
Side
Max-log(p)
VtxMax
Size (mm2)
x
y
z
Correlate
P (corr)
F-value
Effect Size**
Thickness
Lingual
R
−2.488
127,927
1110
19
−69
−7
MJ Years
0.016
10.07
0.063
GWR
Rostral middle frontal,
L
−2.668
42,505
1730
−23
53
10
MJ Grams
0.001
11.09
1.969
Rostral middle frontal
R
−3.565
94,896
2661
39
30
18
MJ Years
0.0002
16.6
0.744
Medial orbitofrontal
R
−3.304
84,773
1368
6
47
−20
MJ Years
0.013
14.92
0.796
LGI
Inferior parietal
L
3.456
122,169
2565
−44
−62
44
MJ Grams
0.015
15.89
0.131
p(corr), family-wise error fully corrected.
The effect sizes were derived from Freesurfer's tool in the significant region of interest using mri_segstats. This was also confirmed manually by using the F-value reported by Freesurfer.
Early vs. late onset marijuana users show divergent morphological patterns based on current marijuana use (measured in grams; MJ grams) in overlapping areas of anterior prefrontal cortex. GWR, gray/white matter border ratio; LGI, local gyrification index.
+
+
+
+ Gray–white matter contrast
+
There were no regions of group differences in gray–white matter contrast by early onset status alone, controlled for age and alcohol use. However, current MJ consumption (grams) and onset status were differentially correlated with gray–white matter contrast in a left anterior dorsal frontal region (p < 0.05, FWE corrected). Increased gray–white contrast with heavier MJ use was seen in the early onset users and the opposite was seen in later onset users (heavier current use linked to decreasing GWR). The same pattern was seen between duration of MJ use in two prefrontal cortex clusters of the right dorsal frontal and medial orbitofrontal area p < 0.05, FWE corrected – more years of MJ use were linked to greater GWR among early users (Fig. 1).
+
+
+
+ Gyrification
+
MJ use onset status alone showed no significant main effects above age and alcohol covariates. However, onset status was correlated with divergent patterns between local gyrification and MJ use, whereby early onset users showed decreasing LGI with increasing MJ consumption and longer duration of use in prefrontal cortex regions p < 0.05, FWE corrected. The left hemisphere clusters encompassed the majority of the length of the middle lateral surface of the left cortex, including motor cortices, parietal lobe and multimodal integration areas (Fig. 1).
+
+
+
+
+ Discussion
+
The present study was designed to characterize the cortical architecture in adolescent onset MJ users by comparing early adolescent onset users to late adolescent onset in MJ use on measures of cortical thickness, gray/white matter contrast and gyrification. The primary finding was that early versus late onset MJ users showed a divergent pattern in cortical thickness, definition of the cortical ribbon and local gyrification with continued use through and beyond adolescent years. Specifically, early onset users showed cortical thickening, enhanced gray/white matter contrast, and decreased gyrification in association with more years of MJ use and current consumption of MJ in grams in frontal and temporal regions – areas that underlie higher order cognition including executive functioning, learning and memory. Findings were above and beyond effects of alcohol and current age, therefore, results are less likely to reflect morphological trends due to aging.
+
Our findings did not find the expected age of onset differences previously reported in marijuana users (Gruber et al., 2012, Gruber et al., 2014). This inconsistency suggests that the age of onset effects may be more robust in brain white matter connectivity (Gruber et al., 2014) and function (Gruber et al., 2012) than brain surface morphometry. To date, the few studies that have described altered cortical morphology in MJ users have led to mixed findings. Mata et al. (2010) identified brain regions with decreased sulcal depth suggestive of lower gyrification in a study of adult MJ users. Jacobus and Tapert (2014) recently reported increased cortical thickness in the entorhinal cortex among 24 adolescent MJ users (mean age = 17.7, mean MJ onset age = 15.4) relative to peer controls. However, the authors also reported a negative relationship between cortical thickness and total MJ use in the right paracentral gyrus, and they observed consistent positive relationships in various brain regions between age of MJ onset and thickness. In the only other known adolescent study of cortical thickness and MJ, Lopez-Larson and colleagues studied 18 adolescent heavy MJ users (similar in age and MJ onset as Jacobus and Tapert, 2014) and reported mixed findings of increased thickness in prefrontal/insula regions and decreased thickness in posterior/temporal lobe areas in the MJ users compared to controls. In contrast to Jacobus and Tapert, 2014, Lopez-Larson et al., 2011 found areas of the frontal lobe and insula that were thinner with increased urine THC metabolites and thicker with earlier age of onset. Select findings from the current study align with aspects of both of these studies, with a consensus supporting findings of a negative dose-dependent relationship between MJ use and cortical thickness. Given the low availability of studies to compare, this consensus is very limited. Although Jacobus et al. and Lopez-Larson et al. found the opposite effect of age of onset on thickness, the pattern of divergence among early vs. late onset users in the current study is more consistent with the latter study, whereby we saw early onset users exhibit thicker cortex with continued MJ use. Taken together, findings of increased thickness related to early MJ onset accompanied by negative dose-dependent relationships with MJ exposure may reflect two distinct processes. One process may be specific to the interactions with cortical development during early adolescence, likely leading to a disruption in pruning, and, the other, specific to the pharmacological effect with heavy chronic MJ use.
+
In the only known study to examine the curvature-morphology of the cortex in adult MJ users, Mata et al. (2010) identified decreased sulcal concavity and thinner sulci in 23 MJ users compared to controls (n = 44), also in prefrontal areas. However, they did not observe significant relationships with age, MJ onset age, or cumulative MJ use. It is interesting that the authors detected group level differences (MJ vs. controls) but no correlations with MJ use characteristics such as dose or age of onset, whereas our primary findings are the consistent effects of continued MJ use differing after early or late adolescent onset. There are substantial methodological explanations for this disparity. For example, the current study did not compare morphology in MJ users to a normative control sample, therefore, it is feasible that group-level differences may emerge with such a comparison. Likewise, we deliberately covaried for current age in order to control for brain changes with aging and thus optimize our interrogation of developmental effects of early onset age and of aspects of continued use.
+
The heterogeneity of MJ effects clearly suggests a multifactorial system of neurobiological processes involved. The primary results uphold that age of onset is a robust variable that differentiates heavy MJ users based on early versus late MJ onset. However, this group distinction relied on current use characteristics. Therefore, in the absence of group-level differences, the interactions between onset age and current use indicates that continued cannabis exposure and early adolescent developmental factors both contribute to a dynamic and sustained departure from what is expected based on developmental studies.
+
Typical synaptic refinement processes during early adolescence are in the context of long-term depression and potentiation of cortical neurons in order to facilitate neuronal remodeling. Thus, the normal course of early adolescent development is uniquely vulnerable to disruption by MJ due to the electrochemical conditions and maturity of brain processes that would not present together again. Cass and colleagues tested the sensitivity of early adolescence cannabinoid exposure in an animal model (Cass et al., 2014). They found that acute administration of cannabinoid agonists in early, middle and late adolescent rats led to a state of frequency-dependent disinhibition of neurons in the frontal cortex in the early-to-middle adolescent rats, but not in the late adolescent rats. Moreover, the authors also noted that adult rats previously exposed to cannabinoid agonists in adolescence displayed comparable neuronal disinhibition. Thus, by changing the inhibitory/excitatory landscape during adolescence, MJ can influence lasting changes to typical cortical remodeling during sensitive early adolescent years.
+
The sequence of pruning and myelination likely plays a formative role in lasting changes from early adolescent onset MJ use. With decreased synaptic elimination, our findings of greater GW border contrast may reflect greater proliferation of myelin at the boundary of the cortical ribbon where non-pruned synapses remained with linked axons. Findings of altered white matter tissue qualities are mixed in adolescent and adult MJ user samples. Some report both increases and decreases in fractional anisotropy (FA) and average water diffusion (Bava et al., 2009) whereas others report consistent decreases in FA among adolescent MJ users (Ashtari et al., 2009, Jacobus et al., 2009) or null findings (Delisi et al., 2006). Two studies of diffusion tensor imaging in adult MJ users reported reduced FA in users compared to controls (Gruber et al., 2011, Gruber et al., 2014). In addition to equivocal findings, research is needed to address the microstructural changes that could result in altered definition of the cortical ribbon. For example, rather than whole brain techniques that assess diffusion measures along major white matter tracts, indices assessing axonal organization along radial and interneuron association fibers along the cortical ribbon are needed. This scenario played out could result in increased gray matter (thicker cortex from disrupted pruning) and the myelination of connections to these spared terminals would result in increased density of white matter at the cortical boundary. Without any known studies of adolescent development of the gray/white tissue contrast at the cortical border to serve as a point of comparison, we speculate that early adolescent disruption of pruning and subsequent myelination of connections at the cortical boundary would be reflected by increased GWR as we saw in the current study.
+
+
+
+ Limitations and conclusions
+
The cross-sectional nature of this study limits causal attributions in terms of what we can infer to be directly related to the effects of MJ. Although a longitudinal design is optimal for addressing brain changes directly due to MJ, cross-sectional studies facilitate data-driven hypotheses that can be assessed directly in prospective studies.
+
It is important to keep in mind that the participants were not explicitly asked for possible years of abstinence during their period of regular use, which may have created possible inflation in reported duration of regular use. However, because the participants provided number of years of “regular” marijuana use, this inherently suggests continued, uninterrupted years of use. Concurrent nicotine use could have also influenced our reported results. But in the absence of a larger sample size and the presence of huge variance in nicotine use in the current sample, we were unable to verify the effect of nicotine use in the reported results.
+
Interpretation of these findings is also limited by the lack of behavioral anchors for the observed morphological effects and lack of information on other aspects of developmental history that could further characterize the effects of marijuana during neurodevelopment. This is further limited by the absence of “expected” patterns based on normative data. Given the varied directions of effects and the small sample size, these findings should be replicated and be viewed as preliminary.
+
To conclude, early MJ use was linked to altered neurodevelopmental patterns in brain regions sub-serving higher-order cognitive process. Clinical implications include need for early, targeted intervention. Given that the most robust results were related to interactions between onset age and continued use through emerging adulthood, harm reduction approaches may be effective in moderating adolescent MJ use to levels that are less likely to cause long-term developmental changes.
+
+
+ Conflict of interest
+
The authors report no conflicts of interest.
+
+
+
+
+ References
+
+
+
+
+ Alemán-Gómez
+ Y.
+
+
+ The human cerebral cortex flattens during adolescence
+
+ 33
+ 38
+ 2013
+ 15004
+ 15010
+ 24048830
+
+
+
+
+
+
+ Anavi-Goffer
+ S.
+
+
+ Mulder
+ J.
+
+
+ The polarised life of the endocannabinoid system in CNS development
+
+ 10
+ 10
+ 2009
+ 1591
+ 1598
+ 19533710
+
+
+
+
+
+
+ Armstrong
+ E.
+
+
+ The ontogeny of human gyrification
+
+ 5
+ 1
+ 1995
+ 56
+ 63
+ 7719130
+
+
+
+
+
+
+ Ashtari
+ M.
+
+
+ Cervellione
+ K.
+
+
+ Cottone
+ J.
+
+
+ Ardekani
+ B.A.
+
+
+ Sevy
+ S.
+
+
+ Kumra
+ S.
+
+
+ Diffusion abnormalities in adolescents and young adults with a history of heavy cannabis use
+
+ 43
+ 3
+ 2009
+ 189
+ 204
+ 19111160
+
+
+
+
+
+
+ Battistella
+ G.
+
+
+ Fornari
+ E.
+
+
+ Annoni
+ J.M.
+
+
+ Chtioui
+ H.
+
+
+ Dao
+ K.
+
+
+ Fabritius
+ M.
+
+ …
+
+ Giroud
+ C.
+
+
+ Long-term effects of cannabis on brain structure
+
+ 39
+ 9
+ 2014
+ 2041
+ 2048
+ 24633558
+
+
+
+
+
+
+ Bava
+ S.
+
+
+ Frank
+ L.R.
+
+
+ McQueeny
+ T.
+
+
+ Schweinsburg
+ B.C.
+
+
+ Schweinsburg
+ A.D.
+
+
+ Tapert
+ S.F.
+
+
+ Altered white matter microstructure in adolescent substance users
+
+ 173
+ 3
+ 2009
+ 228
+ 237
+ pii:S0925-4927(09)00089-4
+ 19699064
+
+
+
+
+
+
+ Caballero
+ A.
+
+
+ Tseng
+ K.Y.
+
+
+ Association of cannabis use during adolescence, prefrontal cb1 receptor signaling, and schizophrenia
+
+ 3
+ 2012
+ 101
+ 22654759
+
+
+
+
+
+
+ CASA
+
+
+ Adolescent Substance Use: America's #1 Public Health Problem
+ 2011
+ The National Center on Addiction and Substance Abuse (CASA) at Columbia University
+ New York
+
+
+
+
+
+
+ Casey
+ B.J.
+
+
+ Tottenham
+ N.
+
+
+ Liston
+ C.
+
+
+ Durston
+ S.
+
+
+ Imaging the developing brain: what have we learned about cognitive development?
+
+ 9
+ 3
+ 2005
+ 104
+ 110
+ pii:S1364-6613(05)00030-6
+ 15737818
+
+
+
+
+
+
+ Cass
+ D.K.
+
+
+ Flores-Barrera
+ E.
+
+
+ Thomases
+ D.R.
+
+
+ Vital
+ W.F.
+
+
+ Caballero
+ A.
+
+
+ Tseng
+ K.Y.
+
+
+ CB1 cannabinoid receptor stimulation during adolescence impairs the maturation of GABA function in the adult rat prefrontal cortex
+
+ 19
+ 5
+ 2014
+ 536
+ 543
+ 24589887
+
+
+
+
+
+
+ Chan
+ G.C.
+
+
+ Hinds
+ T.R.
+
+
+ Impey
+ S.
+
+
+ Storm
+ D.R.
+
+
+ Hippocampal neurotoxicity of Delta9-tetrahydrocannabinol
+
+ 18
+ 14
+ 1998
+ 5322
+ 5332
+ Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/9651215, http://www.jneurosci.org/content/18/14/5322.full.pdf
+ 9651215
+
+
+
+
+
+
+ Cousijn
+ J.
+
+
+ Wiers
+ R.W.
+
+
+ Ridderinkhof
+ K.R.
+
+
+ van den Brink
+ W.
+
+
+ Veltman
+ D.J.
+
+
+ Goudriaan
+ A.E.
+
+
+ Grey matter alterations associated with cannabis use: results of a VBM study in heavy cannabis users and healthy controls
+
+ 59
+ 4
+ 2012
+ 3845
+ 3851
+ 21982932
+
+
+
+
+
+
+ Dale
+ A.M.
+
+
+ Fischl
+ B.
+
+
+ Sereno
+ M.I.
+
+
+ Cortical surface-based analysis I. Segmentation and surface reconstruction
+
+ 9
+ 2
+ 1999
+ 179
+ 194
+ 9931268
+
+
+
+
+
+
+ Delisi
+ L.E.
+
+
+ Bertisch
+ H.C.
+
+
+ Szulc
+ K.U.
+
+
+ Majcher
+ M.
+
+
+ Brown
+ K.
+
+
+ Bappal
+ A.
+
+
+ Ardekani
+ B.A.
+
+
+ A preliminary DTI study showing no brain structural change associated with adolescent cannabis use
+
+ 3
+ 2006
+ 17
+ 16684342
+
+
+
+
+
+
+ Desikan
+ R.S.
+
+
+ Segonne
+ F.
+
+
+ Fischl
+ B.
+
+
+ Quinn
+ B.T.
+
+
+ Dickerson
+ B.C.
+
+
+ Blacker
+ D.
+
+ …
+
+ Killiany
+ R.J.
+
+
+ An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest
+
+ 31
+ 3
+ 2006
+ 968
+ 980
+ 16530430
+
+
+
+
+
+
+ Fergusson
+ D.M.
+
+
+ Horwood
+ L.J.
+
+
+ Early onset cannabis use and psychosocial adjustment in young adults
+
+ 92
+ 3
+ 1997
+ 279
+ 296
+ Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/9219390
+ 9219390
+
+
+
+
+
+
+ Fergusson
+ D.M.
+
+
+ Lynskey
+ M.T.
+
+
+ Horwood
+ L.J.
+
+
+ The short-term consequences of early onset cannabis use
+
+ 24
+ 4
+ 1996
+ 499
+ 512
+ Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/8886945
+ 8886945
+
+
+
+
+
+
+ Filbey
+ F.M.
+
+
+ Aslan
+ S.
+
+
+ Calhoun
+ V.D.
+
+
+ Spence
+ J.S.
+
+
+ Damaraju
+ E.
+
+
+ Caprihan
+ A.
+
+
+ Segall
+ J.
+
+
+ Long-term effects of marijuana use on the brain
+
+ 2014
+
+
+
+
+
+
+ First
+ M.B.
+
+
+ Spitzer
+ R.L.
+
+
+ Miriam
+ G.
+
+
+ Williams
+ J.B.W.
+
+
+ Structured Clinical Interview for DSM-IV-TR Axis I Disorders, Research Version, Non-patient Edition (SCID-I/NP)
+ 2002
+ Biometrics Research, New York State Psychiatric Institute
+ New York
+
+
+
+
+
+
+ Fischl
+ B.
+
+
+ Dale
+ A.M.
+
+
+ Measuring the thickness of the human cerebral cortex from magnetic resonance images
+
+ 97
+ 20
+ 2000
+ 11050
+ 11055
+ 10984517
+
+
+
+
+
+
+ Fischl
+ B.
+
+
+ Salat
+ D.H.
+
+
+ Busa
+ E.
+
+
+ Albert
+ M.
+
+
+ Dieterich
+ M.
+
+
+ Haselgrove
+ C.
+
+ …
+
+ Dale
+ A.M.
+
+
+ Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain
+
+ 33
+ 3
+ 2002
+ 341
+ 355
+ Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/11832223
+ 11832223
+
+
+
+
+
+
+ Fischl
+ B.
+
+
+ Sereno
+ M.I.
+
+
+ Dale
+ A.M.
+
+
+ Cortical surface-based analysis II: Inflation, flattening, and a surface-based coordinate system
+
+ 9
+ 2
+ 1999
+ 195
+ 207
+ 9931269
+
+
+
+
+
+
+ Giedd
+ J.N.
+
+
+ The teen brain: insights from neuroimaging
+
+ 42
+ 4
+ 2008
+ 335
+ 343
+ 18346658
+
+
+
+
+
+
+ Giorgio
+ A.
+
+
+ Watkins
+ K.E.
+
+
+ Chadwick
+ M.
+
+
+ James
+ S.
+
+
+ Winmill
+ L.
+
+
+ Douaud
+ G.
+
+ …
+
+ James
+ A.C.
+
+
+ Longitudinal changes in grey and white matter during adolescence
+
+ 49
+ 1
+ 2010
+ 94
+ 103
+ pii:S1053-8119(09)00864-7
+ 19679191
+
+
+
+
+
+
+ Gogtay
+ N.
+
+
+ Giedd
+ J.N.
+
+
+ Lusk
+ L.
+
+
+ Hayashi
+ K.M.
+
+
+ Greenstein
+ D.
+
+ …
+
+ Thompson
+ P.M.
+
+
+ Dynamic mapping of human cortical development during childhood through early adulthood
+
+ 101
+ 21
+ 2004
+ 8174
+ 8179
+ 15148381
+
+
+
+
+
+
+ Gruber
+ S.A.
+
+
+ Dahlgren
+ M.K.
+
+
+ Sagar
+ K.A.
+
+
+ Gonenc
+ A.
+
+
+ Killgore
+ W.D.
+
+
+ Age of onset of marijuana use impacts inhibitory processing
+
+ 511
+ 2
+ 2012
+ 89
+ 94
+ 22306089
+
+
+
+
+
+
+ Gruber
+ S.A.
+
+
+ Dahlgren
+ M.K.
+
+
+ Sagar
+ K.A.
+
+
+ Gonenc
+ A.
+
+
+ Lukas
+ S.E.
+
+
+ Worth the wait: effects of age of onset of marijuana use on white matter and impulsivity
+
+ 231
+ 8
+ 2014
+ 1455
+ 1465
+ 24190588
+
+
+
+
+
+
+ Gruber
+ S.A.
+
+
+ Silveri
+ M.M.
+
+
+ Dahlgren
+ M.K.
+
+
+ Yurgelun-Todd
+ D.
+
+
+ Why so impulsive? White matter alterations are associated with impulsivity in chronic marijuana smokers
+
+ 19
+ 3
+ 2011
+ 231
+ 242
+ 21480730
+
+
+
+
+
+
+ Hasan
+ K.M.
+
+
+ Sankar
+ A.
+
+
+ Halphen
+ C.
+
+
+ Kramer
+ L.A.
+
+
+ Brandt
+ M.E.
+
+
+ Juranek
+ J.
+
+ …
+
+ Ewing-Cobbs
+ L.
+
+
+ Development and organization of the human brain tissue compartments across the lifespan using diffusion tensor imaging
+
+ 18
+ 16
+ 2007
+ 1735
+ 1739
+ 17921878
+
+
+
+
+
+
+ Jacobus
+ J.
+
+
+ McQueeny
+ T.
+
+
+ Bava
+ S.
+
+
+ Schweinsburg
+ B.C.
+
+
+ Frank
+ L.R.
+
+
+ Yang
+ T.T.
+
+
+ Tapert
+ S.F.
+
+
+ White matter integrity in adolescents with histories of marijuana use and binge drinking
+
+ 31
+ 6
+ 2009
+ 349
+ 355
+ pii:S0892-0362(09)00145-7
+ 19631736
+
+
+
+
+
+
+ Jacobus
+ J.
+
+
+ Tapert
+ S.F.
+
+
+ Effects of cannabis on the adolescent brain
+
+ 20
+ 13
+ 2014
+ 2186
+ 2193
+ 23829363
+
+
+
+
+
+
+ Johnston
+ L.D.
+
+
+ O’Malley
+ P.M.
+
+
+ Miech
+ R.A.
+
+
+ Bachman
+ J.G.
+
+
+ Schulenberg
+ J.E.
+
+
+ Monitoring the Future National Survey Results on Drug Use, 1975–2013
+ 2014
+ Institute for Social Research, The University of Michigan
+ Ann Arbor, MI
+
+
+
+
+
+
+ Klein
+ D.
+
+
+ Adolescent brain maturation and cortical folding: evidence for reductions in gyrification
+
+ 9
+ 1
+ 2014
+ e84914
+ 24454765
+
+
+
+
+
+
+ Lebel
+ C.
+
+
+ Caverhill-Godkewitsch
+ S.
+
+
+ Beaulieu
+ C.
+
+
+ Age-related variations of white matter tracts
+
+ 52
+ 1
+ 2010
+ 20
+ 31
+ pii:S1053-8119(10)00359-9
+ 20362683
+
+
+
+
+
+
+ Lopez-Larson
+ M.P.
+
+
+ Altered prefrontal and insular cortical thickness in adolescent marijuana users
+
+ 220
+ 1
+ 2011
+ 164
+ 172
+ 21310189
+
+
+
+
+
+
+ Lopez-Moreno
+ J.A.
+
+
+ Gonzalez-Cuevas
+ G.
+
+
+ Moreno
+ G.
+
+
+ Navarro
+ M.
+
+
+ The pharmacology of the endocannabinoid system: functional and structural interactions with other neurotransmitter systems and their repercussions in behavioral addiction
+
+ 13
+ 2
+ 2008
+ 160
+ 187
+ 18422831
+
+
+
+
+
+
+ Lorenzetti
+ V.
+
+
+ Lubman
+ D.I.
+
+
+ Whittle
+ S.
+
+
+ Solowij
+ N.
+
+
+ Yucel
+ M.
+
+
+ Structural MRI findings in long-term cannabis users: what do we know?
+
+ 45
+ 11
+ 2010
+ 1787
+ 1808
+ 20590400
+
+
+
+
+
+
+ Mata
+ I.
+
+
+ Gyrification brain abnormalities associated with adolescence and early-adulthood cannabis use
+
+ 1317
+ 2010
+ 297
+ 304
+ 20045399
+
+
+
+
+
+
+ Matochik
+ J.A.
+
+
+ Eldreth
+ D.A.
+
+
+ Cadet
+ J.L.
+
+
+ Bolla
+ K.I.
+
+
+ Altered brain tissue composition in heavy marijuana users
+
+ 77
+ 1
+ 2005
+ 23
+ 30
+ Retrieved from http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T63-4DCN07C-1-1&_cdi=5019&_user=2629161&_pii=S0376871604002066&_orig=search&_coverDate=01%2F07%2F2005&_sk=999229998&view=c&wchp=dGLbVzb-zSkWA&md5=f38ff3ce0e5bc6e93f01659d96437f10&ie=/sdarticle.pdf
+ 15607838
+
+
+
+
+
+
+ Patton
+ G.C.
+
+
+ Coffey
+ C.
+
+
+ Carlin
+ J.B.
+
+
+ Degenhardt
+ L.
+
+
+ Lynskey
+ M.
+
+
+ Hall
+ W.
+
+
+ Cannabis use and mental health in young people: cohort study
+
+ 325
+ 7374
+ 2002
+ 1195
+ 1198
+ 12446533
+
+
+
+
+
+
+ Paus
+ T.Å.
+
+
+ Mapping brain maturation and cognitive development during adolescence
+
+ 9
+ 2
+ 2005
+ 60
+ 68
+ 15668098
+
+
+
+
+
+
+ Sagar
+ K.A.
+
+
+ Dahlgren
+ M.K.
+
+
+ Gonenc
+ A.
+
+
+ Racine
+ M.T.
+
+
+ Dreman
+ M.W.
+
+
+ Gruber
+ S.A.
+
+
+ The impact of initiation: early onset marijuana smokers demonstrate altered Stroop performance and brain activation
+
+ 2015
+
+
+
+
+
+
+ Salat
+ D.H.
+
+
+ Lee
+ S.Y.
+
+
+ van der Kouwe
+ A.J.
+
+
+ Greve
+ D.N.
+
+
+ Fischl
+ B.
+
+
+ Rosas
+ H.D.
+
+
+ Age-associated alterations in cortical gray and white matter signal intensity and gray to white matter contrast
+
+ 48
+ 1
+ 2009
+ 21
+ 28
+ 19580876
+
+
+
+
+
+
+ SAMHSA
+
+
+ Results from the 2013 National Survey on Drug Use and Health: Summary of National Findings NSDUH Series H-48, HHS Publication No. (SMA) 14-4863
+ 2014
+ Substance Abuse and Mental Health Services Administration
+ Rockville, MD
+
+
+
+
+
+
+ Schaer
+ M.
+
+
+ Cuadra
+ M.B.
+
+
+ Schmansky
+ N.
+
+
+ Fischl
+ B.
+
+
+ Thiran
+ J.P.
+
+
+ Eliez
+ S.
+
+
+ How to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index
+
+ 59
+ 2012
+ e3417
+ 22230945
+
+
+
+
+
+
+ Schmansky
+ N.
+
+
+ Stevens
+ A.
+
+
+ Subramaniam
+ K.
+
+
+ Greve
+ D.N.
+
+
+ Kakunoori
+ S.
+
+
+ Pacheco
+ J.
+
+
+ Troubleshooting Your Output
+ 2010
+ FsTutorial
+ Retrieved from http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/TroubleshootingData
+
+
+
+
+
+
+ Selemon
+ L.D.
+
+
+ A role for synaptic plasticity in the adolescent development of executive function
+
+ 3
+ 2013
+ e238
+ 23462989
+
+
+
+
+
+
+ Shaw
+ P.
+
+
+ Kabani
+ N.J.
+
+
+ Lerch
+ J.P.
+
+
+ Eckstrand
+ K.
+
+
+ Lenroot
+ R.
+
+
+ Gogtay
+ N.
+
+ …
+
+ Wise
+ S.P.
+
+
+ Neurodevelopmental trajectories of the human cerebral cortex
+
+ 28
+ 14
+ 2008
+ 3586
+ 3594
+ pii:28/14/3586
+ 18385317
+
+
+
+
+
+
+ Sobell
+ L.C.
+
+
+ Sobell
+ M.B.
+
+
+ Timeline follow-back: a technique for assessing self-reported alcohol consumption
+
+
+ Litten
+ R.Z.
+
+
+ Allen
+ J.P.
+
+
+
+ 1992
+ Humana Press
+ Totowa, NJ, USA
+ 41
+ 72
+
+
+
+
+
+
+ Tapert
+ S.F.
+
+
+ Schweinsburg
+ A.D.
+
+
+ Barlett
+ V.C.
+
+
+ Brown
+ S.A.
+
+
+ Frank
+ L.R.
+
+
+ Brown
+ G.G.
+
+
+ Meloy
+ M.J.
+
+
+ Blood oxygen level dependent response and spatial working memory in adolescents with alcohol use disorders
+
+ 28
+ 10
+ 2004
+ 1577
+ 1586
+ Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/15597092
+ 15597092
+
+
+
+
+
+
+ Verdurand
+ M.
+
+
+ Nguyen
+ V.
+
+
+ Stark
+ D.
+
+
+ Zahra
+ D.
+
+
+ Gregoire
+ M.C.
+
+
+ Greguric
+ I.
+
+
+ Zavitsanou
+ K.
+
+
+ Comparison of cannabinoid CB(1) receptor binding in adolescent and adult rats: a positron emission tomography study using [F]MK-9470
+
+ 2011
+ 2011
+ 548123
+ 22187642
+
+
+
+
+
+
+ Walter
+ L.
+
+
+ Franklin
+ A.
+
+
+ Witting
+ A.
+
+
+ Wade
+ C.
+
+
+ Xie
+ Y.
+
+
+ Kunos
+ G.
+
+ …
+
+ Stella
+ N.
+
+
+ Nonpsychotropic cannabinoid receptors regulate microglial cell migration
+
+ 23
+ 4
+ 2003
+ 1398
+ 1405
+ pii:23/4/1398
+ 12598628
+
+
+
+
+
+
+ Wechsler
+ D.
+
+
+ Manual for the Wechsler Abbreviated Scale of Intelligence
+ 1999
+ The Psychological Corporation
+ San Antonio, TX
+
+
+
+
+
+
+ Yucel
+ M.
+
+
+ Solowij
+ N.
+
+
+ Respondek
+ C.
+
+
+ Whittle
+ S.
+
+
+ Fornito
+ A.
+
+
+ Pantelis
+ C.
+
+
+ Lubman
+ D.I.
+
+
+ Regional brain abnormalities associated with long-term heavy cannabis use
+
+ 65
+ 6
+ 2008
+ 694
+ 701
+ 18519827
+
+
+
+
+
+ Supplementary data
+
The following are Supplementary data to this article:
+
+
+ Acknowledgements
+
This research was funded by the National Institute on Drug Abuse (R01 DA030344, Filbey). We would like to thank all the participants who volunteered for this study. We are also very grateful to Talha Alvi, Sina Aslan, Jessica Baine, Collette Bice, Vicki Germer, Ariel Ketcherside, Alison King, Brittany Kuhn, Tyler Rhinehardt, Wing Ting To and the team of lab interns for their assistance with recruitment, running participants and data management.
+
+
+
+
+
Supplementary material related to this article can be found, in the online version, at http://dx.doi.org/10.1016/j.dcn.2015.10.001.
+
+
+
+
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000.csv b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000.csv
new file mode 100644
index 0000000..f023fbb
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000.csv
@@ -0,0 +1,21 @@
+Measure,Early onset (n = 20),Early onset (n = 20),Early onset (n = 20),Late onset (n = 22),Late onset (n = 22),Late onset (n = 22),p-Value,Effect size???,Statistic
+Unnamed: 0_level_1,Mean,(SD),Min–Max,Mean,(SD),Min–Max,Unnamed: 7_level_1,Unnamed: 8_level_1,|t/U
+Age,32.50,(8.01),21–50,30.25,(7.19),21–47,0.316,0.302,t = 1.01
+Education (years),12.91,(2.54),8–18,13.26,(2.40),10–19,0.651,0.144,t = 0.456
+Gender (male),55%,,,73%,,,0.241,0.034,χ2 = 1.41
+Ethnicity (% Caucasian),50%,,,50%,,,0.566,0.008,χ2 = 0.336
+IQ???,108,(9.99),88–124,105,(13.54),83–129,0.351,0.298,t = 0.94
+Age of first MJ use???,13.18,(1.89),9–15,16.90,(1.48),16–21,<0.001,0.866,U = 0
+Age of regular MJ use???,16.50,(3.57),9–25,19.00,(4.29),16–36,0.004,0.439,U = 108
+Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days,Substance use in the last 60 days
+MJ grams (daily),2.14,(1.79),0.50–7.50,1.65,(1.21),0.46–4.23,0.338,0.083,U = 182
+# EtOH drinks,44.09,(76.30),0–310,38.85,(58.61),0–183,0.588,0.039,U = 198.05
+Max # EtOH drinks,6.62,(7.21),0–31,6.25,(5.80),0–21,0.8,0.095,U = 210
+# EtOH drinking days,11.09,(16.69),0–59,9.84,(13.51),0–60,0.537,0.006,U = 185.5
+# Binge EtOH drinking days,4.36,(12.50),0–59,2.90,(5.53),0–19,0.968,0.035,U = 218.5
+# EtOH drinks per day,2.99,(2.27),0–7.40,3.56,(3.29),0–14.00,0.82,0.289,U = 211
+# Cigarette days,1.18,(3.72),0–17,2.95,(1.17),0–21,0.06,0.294,U = 159
+# Cigarettes per day,0.22,(0.55),0–2.00,0.78,(1.23),0–4.50,0.057,0.296,U = 158
+Max # cigarettes,0.25,(0.61),0–2,0.96,(1.60),0–6,0.054,0.290,U = 157.5
+Illicit drug use/past 90 days,14%,,,5%,,,,,
+Lifetime illicit drug use,73%,,,75%,,,,,
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000_info.json b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000_info.json
new file mode 100644
index 0000000..950f279
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_000_info.json
@@ -0,0 +1 @@
+{"table_id": "tbl0005", "table_label": "Table 1", "table_caption": "Sample characteristics. MJ, marijuana.", "table_foot": "*IQ scores derived from Wechsler Abbreviated Scale of Intelligence Vocabulary and Matrix Reasoning subtests.", "n_header_rows": 2, "table_data_file": "table_000.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001.csv b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001.csv
new file mode 100644
index 0000000..abbbe2d
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001.csv
@@ -0,0 +1,5 @@
+Measure,Early onset,Late onset
+First MJ use,r = 0.038,r = 0.189
+Regular MJ use,r = 0.289,r = 0.203
+MJ years of use,r = 0.898???,r = 0.623???
+MJ grams,r = 0.123,r = 0.206
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001_info.json b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001_info.json
new file mode 100644
index 0000000..c6bb046
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_001_info.json
@@ -0,0 +1 @@
+{"table_id": "tbl0010", "table_label": "Table 2", "table_caption": "The correlations between current age and all MJ use variables.", "table_foot": "*p\u00a0<\u00a00.001.", "n_header_rows": 1, "table_data_file": "table_001.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002.csv b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002.csv
new file mode 100644
index 0000000..7445b91
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002.csv
@@ -0,0 +1,6 @@
+Measure,"Label@Max, Extended coverage",Side,Max-log(p),VtxMax,Size (mm2),x,y,z,Correlate,P (corr),F-value,Effect Size???
+Thickness,Lingual,R,−2.488,"127,927",1110,19,−69,−7,MJ Years,0.016,10.07,0.063
+GWR,"Rostral middle frontal,",L,−2.668,"42,505",1730,−23,53,10,MJ Grams,0.001,11.09,1.969
+,Rostral middle frontal,R,−3.565,"94,896",2661,39,30,18,MJ Years,0.0002,16.6,0.744
+,Medial orbitofrontal,R,−3.304,"84,773",1368,6,47,−20,MJ Years,0.013,14.92,0.796
+LGI,Inferior parietal,L,3.456,"122,169",2565,−44,−62,44,MJ Grams,0.015,15.89,0.131
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002_info.json b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002_info.json
new file mode 100644
index 0000000..bfd3e11
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/table_002_info.json
@@ -0,0 +1 @@
+{"table_id": "tbl0015", "table_label": "Table 3", "table_caption": "Clusters of significant age of onset\u00a0\u00d7\u00a0marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.", "table_foot": "p(corr), family-wise error fully corrected.", "n_header_rows": 1, "table_data_file": "table_002.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/tables.xml b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/tables.xml
new file mode 100644
index 0000000..35835ea
--- /dev/null
+++ b/tests/data/sample_inputs/3qT3nzK9bLZ7/source/pubget/tables/tables.xml
@@ -0,0 +1,2 @@
+
+4691364265074334691364S1878-9293(15)00091-210.1016/j.dcn.2015.10.001tbl0005Table 1Sample characteristics. MJ, marijuana.*IQ scores derived from Wechsler Abbreviated Scale of Intelligence Vocabulary and Matrix Reasoning subtests.
Sample characteristics. MJ, marijuana.
Measure
Early onset (n = 20)
Late onset (n = 22)
p-Value
Effect size***
Statistic
Mean
(SD)
Min–Max
Mean
(SD)
Min–Max
|t/U
Age
32.50
(8.01)
21–50
30.25
(7.19)
21–47
0.316
0.302
t = 1.01
Education (years)
12.91
(2.54)
8–18
13.26
(2.40)
10–19
0.651
0.144
t = 0.456
Gender (male)
55%
73%
0.241
0.034
χ2 = 1.41
Ethnicity (% Caucasian)
50%
50%
0.566
0.008
χ2 = 0.336
IQ*
108
(9.99)
88–124
105
(13.54)
83–129
0.351
0.298
t = 0.94
Age of first MJ use**
13.18
(1.89)
9–15
16.90
(1.48)
16–21
<0.001
0.866
U = 0
Age of regular MJ use**
16.50
(3.57)
9–25
19.00
(4.29)
16–36
0.004
0.439
U = 108
Substance use in the last 60 days
MJ grams (daily)
2.14
(1.79)
0.50–7.50
1.65
(1.21)
0.46–4.23
0.338
0.083
U = 182
# EtOH drinks
44.09
(76.30)
0–310
38.85
(58.61)
0–183
0.588
0.039
U = 198.05
Max # EtOH drinks
6.62
(7.21)
0–31
6.25
(5.80)
0–21
0.8
0.095
U = 210
# EtOH drinking days
11.09
(16.69)
0–59
9.84
(13.51)
0–60
0.537
0.006
U = 185.5
# Binge EtOH drinking days
4.36
(12.50)
0–59
2.90
(5.53)
0–19
0.968
0.035
U = 218.5
# EtOH drinks per day
2.99
(2.27)
0–7.40
3.56
(3.29)
0–14.00
0.82
0.289
U = 211
# Cigarette days
1.18
(3.72)
0–17
2.95
(1.17)
0–21
0.06
0.294
U = 159
# Cigarettes per day
0.22
(0.55)
0–2.00
0.78
(1.23)
0–4.50
0.057
0.296
U = 158
Max # cigarettes
0.25
(0.61)
0–2
0.96
(1.60)
0–6
0.054
0.290
U = 157.5
Illicit drug use/past 90 days
14%
5%
Lifetime illicit drug use
73%
75%
IQ scores derived from Wechsler Abbreviated Scale of Intelligence Vocabulary and Matrix Reasoning subtests.
p < .05; SS, standard score; |t|, absolute value of student's t, U is the Mann–Whitney U's score.
The effect sizes of the above table were calculated either based on mean differences if normally distributed, correlation coefficient or F-value score using the default Cohen's effect size formula for respective metrics.
Table 1
Sample characteristics. MJ, marijuana.
Measure
Early onset (n = 20)
Late onset (n = 22)
p-Value
Effect size???
Statistic
Mean
(SD)
Min–Max
Mean
(SD)
Min–Max
|t/U
Age
32.50
(8.01)
21–50
30.25
(7.19)
21–47
0.316
0.302
t = 1.01
Education (years)
12.91
(2.54)
8–18
13.26
(2.40)
10–19
0.651
0.144
t = 0.456
Gender (male)
55%
73%
0.241
0.034
χ2 = 1.41
Ethnicity (% Caucasian)
50%
50%
0.566
0.008
χ2 = 0.336
IQ???
108
(9.99)
88–124
105
(13.54)
83–129
0.351
0.298
t = 0.94
Age of first MJ use???
13.18
(1.89)
9–15
16.90
(1.48)
16–21
<0.001
0.866
U = 0
Age of regular MJ use???
16.50
(3.57)
9–25
19.00
(4.29)
16–36
0.004
0.439
U = 108
Substance use in the last 60 days
MJ grams (daily)
2.14
(1.79)
0.50–7.50
1.65
(1.21)
0.46–4.23
0.338
0.083
U = 182
# EtOH drinks
44.09
(76.30)
0–310
38.85
(58.61)
0–183
0.588
0.039
U = 198.05
Max # EtOH drinks
6.62
(7.21)
0–31
6.25
(5.80)
0–21
0.8
0.095
U = 210
# EtOH drinking days
11.09
(16.69)
0–59
9.84
(13.51)
0–60
0.537
0.006
U = 185.5
# Binge EtOH drinking days
4.36
(12.50)
0–59
2.90
(5.53)
0–19
0.968
0.035
U = 218.5
# EtOH drinks per day
2.99
(2.27)
0–7.40
3.56
(3.29)
0–14.00
0.82
0.289
U = 211
# Cigarette days
1.18
(3.72)
0–17
2.95
(1.17)
0–21
0.06
0.294
U = 159
# Cigarettes per day
0.22
(0.55)
0–2.00
0.78
(1.23)
0–4.50
0.057
0.296
U = 158
Max # cigarettes
0.25
(0.61)
0–2
0.96
(1.60)
0–6
0.054
0.290
U = 157.5
Illicit drug use/past 90 days
14%
5%
Lifetime illicit drug use
73%
75%
*IQ scores derived from Wechsler Abbreviated Scale of Intelligence Vocabulary and Matrix Reasoning subtests.**p < .05; SS, standard score; |t|, absolute value of student's t, U is the Mann–Whitney U's score.***The effect sizes of the above table were calculated either based on mean differences if normally distributed, correlation coefficient or F-value score using the default Cohen's effect size formula for respective metrics.tbl0010Table 2The correlations between current age and all MJ use variables.*p < 0.001.
The correlations between current age and all MJ use variables.
Measure
Early onset
Late onset
First MJ use
r = 0.038
r = 0.189
Regular MJ use
r = 0.289
r = 0.203
MJ years of use
r = 0.898*
r = 0.623**
MJ grams
r = 0.123
r = 0.206
p < 0.001.
p < 0.005.
Table 2
The correlations between current age and all MJ use variables.
Measure
Early onset
Late onset
First MJ use
r = 0.038
r = 0.189
Regular MJ use
r = 0.289
r = 0.203
MJ years of use
r = 0.898???
r = 0.623???
MJ grams
r = 0.123
r = 0.206
*p < 0.001.**p < 0.005.tbl0015Table 3Clusters of significant age of onset × marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.p(corr), family-wise error fully corrected.
Clusters of significant age of onset × marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.
Measure
Label@Max, Extended coverage
Side
Max-log(p)
VtxMax
Size (mm2)
x
y
z
Correlate
P (corr)
F-value
Effect Size**
Thickness
Lingual
R
−2.488
127,927
1110
19
−69
−7
MJ Years
0.016
10.07
0.063
GWR
Rostral middle frontal,
L
−2.668
42,505
1730
−23
53
10
MJ Grams
0.001
11.09
1.969
Rostral middle frontal
R
−3.565
94,896
2661
39
30
18
MJ Years
0.0002
16.6
0.744
Medial orbitofrontal
R
−3.304
84,773
1368
6
47
−20
MJ Years
0.013
14.92
0.796
LGI
Inferior parietal
L
3.456
122,169
2565
−44
−62
44
MJ Grams
0.015
15.89
0.131
p(corr), family-wise error fully corrected.
The effect sizes were derived from Freesurfer's tool in the significant region of interest using mri_segstats. This was also confirmed manually by using the F-value reported by Freesurfer.
Table 3
Clusters of significant age of onset × marijuana use interactions. GWR, gray/white matter border ratio; LGI, local gyrification index.
Understanding the meaning of abstract mathematical symbols is a cornerstone of arithmetic learning in children. Studies have long focused on the role of spatial intuitions in the processing of numerals. However, it has been argued that such intuitions may also underlie symbols that convey fundamental arithmetic concepts, such as arithmetic operators. In the present cross-sectional study, we used fMRI to investigate how and when associations between arithmetic operators and brain regions processing spatial information emerge in children from 3rd to 10th grade. We found that the mere perception of a ‘+’ sign elicited grade-related increases of spatial activity in the right hippocampus. That is, merely perceiving ‘+’ signs – without any operands – elicited enhanced hippocampal activity after around 7th grade (12–13 years old). In these children, hippocampal activity in response to a ‘+’ sign was further correlated with the degree to which calculation performance was facilitated by the preview of that sign before an addition problem, an effect termed operator-priming. Grade-related increases of hippocampal spatial activity were operation-specific because they were not observed with ‘×’ signs, which might evoke rote retrieval rather than numerical manipulation. Our study raises the possibility that hippocampal spatial mechanisms help build associations between some arithmetic operators and space throughout age and/or education.
Humans are unique in their ability to represent abstract mathematical concepts by culturally invented symbols, such as Arabic numerals and arithmetic signs. Because these symbols are arbitrary, learning the relationship between their identity and the concept they represent is a challenge during early math education in children. Most prior studies have focused on the mechanisms supporting the acquisition of symbols representing numerical quantities (Piazza et al., 2007, Ansari, 2008, Holloway and Ansari, 2009, Lyons and Ansari, 2009, Mundy and Gilmore, 2009). However, efficient processing of symbols that convey fundamental arithmetic concepts (i.e., operators) may be an important and largely neglected aspect of arithmetic skills. This is suggested by the operator-priming effect (Roussel et al., 2002, Fayol and Thevenot, 2012, Mathieu et al., 2017), whereby the anticipated presentation of a ‘+’ or ‘-’ sign 150 ms before a single-digit addition or subtraction problem facilitates problem-solving in adults.
What aspect of the processing of an operator may cause the operator-priming effect in adults? A first possibility is that an arithmetic sign may automatically evoke a network of facts. For example, the perception of a ‘+’ or ‘-’ sign might pre-activate a network of additive or subtractive facts that would have been built in declarative memory after years of practice (Campbell and Xue, 2001, Ashcraft, 1992). Pre-activating such a network would facilitate the retrieval of the answer from memory when operands are presented. A second possibility is that an arithmetic sign may prime a specific procedure that would have been “automatized” after its repeated practice during arithmetic learning. For instance, Fayol and Thevenot argued that perceiving a ‘+’ or ‘-’ sign might trigger an automatized procedure that could be “linked to the convocation of the mental number line and could correspond to a preparation for a quick left-to-right or right-to-left browsing of this mental line” (Fayol and Thevenot, 2012). This proposal echoes the idea that adding or subtracting numbers involves rightward and leftward shifts of attention from a source to a target number along a mental map of numbers oriented from left to right, i.e., the mental number line (MNL) (Hubbard et al., 2005, Masson and Pesenti, 2014, Mathieu et al., 2016, Pinheiro-Chagas et al., 2017). Pre-activating such a procedure would result in a facilitation of subsequent calculation when operands are presented, thereby explaining the operator-priming effect.
Interestingly, two lines of evidence favor the procedural over the declarative interpretation of the operator-priming effect. First, the effect is not observed with the ‘×’ sign and multiplication problems (Roussel et al., 2002, Mathieu et al., 2017). Multiplication problems, however, are explicitly learned by rote in school and multiplication is unanimously viewed as the operation having the strongest association with a network of facts in memory (Campbell and Xue, 2001, Galfano et al., 2003, Thibodeau et al., 1996). Therefore, the lack of operator-priming effect for multiplication problems is difficult to reconcile with the idea that the effect is due to associations between operators and networks of stored facts. Second, in line with Fayol and Thevenot’s proposal that ‘+’ and ‘−’ signs may prime a spatial scanning of the MNL, a recent study suggests that ‘+’ and ‘−’ signs do evoke spatial intuitions. Specifically, Pinhas et al. (2014) found that, when instructed to categorize ‘+’ and ‘−’ signs with left-hand or right-hand responses, adults tend to respond faster to ‘+’ signs with the right hand than with the left hand, whereas they tend to respond faster to ‘-’ signs with the left hand than with the right hand (Pinhas et al., 2014). Thus, ‘+’ and ‘−’ signs appear to have some automatic associations with the right and left sides of space, respectively.
Using fMRI, we recently found that such spatial associations may stem from the fact that some arithmetic operators are automatically processed in brain regions involved in spatial attention in adults. We showed that the mere perception of a ‘+’ sign elicits greater activity than the mere perception of a ‘×’ sign in brain regions underlying overt spatial attention. These included the frontal eye fields (FEF) and the posterior superior parietal lobule (PSPL) (Mathieu et al., 2017). Thus, perceiving a ‘+’ sign (but not a ‘×’ sign) may be associated with a deployment of spatial attention in educated adults. Therefore, the rightward shifts of attention that have been posited to underlie addition problem-solving (Hubbard et al., 2005, Masson and Pesenti, 2014, Mathieu et al., 2016) might be primed by the mere preview of the addition sign (but not by the preview of a multiplication sign because multiplication is typically learned by rote and unlikely to be associated with movements along the MNL). Overall, there is mounting evidence that at least some arithmetic operators (e.g., ‘+’ but not ‘×’ signs) evoke spatial intuitions in adults, and that these intuitions may relate to the operator-priming effect.
However, associations between operators and space are arguably not innate. Therefore, a fundamental outstanding question is how and when such associations emerge in the developing brain. To answer that question, we studied 34 children from 3rd to 10th grade while they performed 3 tasks. First, fMRI activity was measured while children were instructed to make eye saccades towards visually presented targets. This allowed us to precisely localize several regions of interest (ROIs) involved in spatial attention across children. Second, fMRI activity was measured in these spatial attention ROIs while children were presented with trials in which a ‘+’ sign was displayed without any operands (hereafter addition sign-only trials). As in our previous study in adults (Mathieu et al., accepted), activity during the perception of addition sign-only trials was compared to activity associated with trials in which a ‘×’ sign was displayed without any operands (hereafter multiplication sign-only trials) because these do not appear to evoke any specific intuitions in adults (Fayol and Thevenot 2012). This allowed us to identify the spatial attention ROIs in which activity in response to a ‘+’ sign (as compared to a ‘×’ sign) increases with age and/or education, as well as the developmental time course of these effects.1 Third, outside of the scanner, we asked subjects to perform an operator-priming task and measured the correlation between inter-individual differences in the size of the operator-priming effect and inter-individual differences in sign-related activity in spatial attention ROIs as a function of grade. This allowed us to evaluate when sign-related activity in spatial attention ROIs leads to an operator-priming effect in children.
2. Material and methods
2.1. Participants
Forty-two right-handed children from 3rd to 10th grade participated in the study. All were native French speakers. Participants did not have prior history of neurological disease, psychiatric disorders, learning disabilities or attention deficits. All children and parents provided written informed consent to participate in the study, which was approved by the local ethics committee (CPP Sud-Est-II). Families received 80€ for their participation. Data from 8 subjects were excluded because of excessive head-movement in the scanner (see criteria in the Section 2.7., n = 3), poor whole-brain coverage (i.e. susceptibility artefacts from dental braces, n = 3) and unacceptably low performance during the task (i.e., lower than 50% accuracy on the sign-plus-operand trials, n = 2). Therefore, the final sample consisted of 34 children (20 males) from 3rd to 10th grade (age range: 8–15, mean age = 11.37, SD = 1.84). For each child, a continuous measure of grade was calculated by taking into account the specific date within the grade year when that child was scanned. The whole sample (n = 34) was evenly split into three groups as a function of grade: 11 children were from the ‘lower grades’ group (grade 3.2–5.4; mean = 4.4), 11 children were from the ‘intermediate grades’ group (grade 5.6–6.9; mean = 6.2), and 12 children were from the ‘higher grades’ group (grade 7.6–10.2; mean = 8.5).
2.2. Standardized measures
Children were administered standardized tests of intellectual and arithmetic abilities to ensure that there were no age differences with respect to those measures. Full-scale IQ was measured using the NEMI-2 (Cognet, 2006). Basic arithmetic knowledge was evaluated with the Math-Fluency subtest of the Woodcock-Johnson-III Tests of Achievement (WJ-III) (Woodcock et al., 2001). Across all participants, standardized (i.e., age-normalized) scores on IQ (mean = 112; SD = 10) and Math Fluency (mean = 106; SD = 16) tests were within the normal range. One-way ANOVAs with the between-subject factor group (lower, intermediate, higher grades) revealed no main effect of group on IQ (F(2,31) = 0.591, p= 0.560, BF10 = 0.29), indicating that age-normalized intellectual abilities were similar across groups. However, there was a main effect of group on Math Fluency (F(2,31) = 5.867, p= 0.007, BF10 = 7.24): Children from intermediate grades had a higher age-normalized score (mean = 118; SD = 18) than children from lower (mean = 100; SD = 11) and higher grades (mean = 100; SD = 13). Therefore, we included standardized Math-Fluency scores as nuisance covariate in all of our analyses.
2.3. Behavioral session
After standardized testing, children participated in a behavioral session during which they performed an operator-priming task adapted from Fayol and Thevenot (2012) and Roussel et al. (2002). Children were asked to evaluate 56 single-digit addition and 56 multiplication problems composed of operands between 2 and 9. Problems were presented in both commutative orders. Tie problems were excluded. Problems with a sum smaller than or equal to 11 and a product smaller or equal to 24 were considered small. Other problems were considered large.
In each trial, a problem was presented with an answer (Fig. 1a). The arithmetic sign was presented either 150 ms before (Negative SOA condition) or at the same time (Null SOA condition) as the operands (Fig. 1a). All problems were presented once in both SOA condition with a valid answer. Twenty-eight addition and 28 multiplication problems were also presented in both SOA condition with an invalid answer (obtained by adding or subtracting 1 to or from the valid answer). Trials were pseudorandomly ordered so that no more than three problems of the same type appeared consecutively. Problems with an invalid answer were randomly chosen across subjects and the order of blocks was counter-balanced between subjects. The experiment started with 8 practice trials.
The experiment was controlled by Presentation software (Neurobehavioral Systems, Albany, CA). Problems were displayed in white Arial 60-point font on a black background. All trials started with the presentation of a white central fixation dot for 1500 ms, immediately followed by a red central fixation dot for 1000 ms signaling that the problem was about to be presented, either in the negative SOA condition or in the null SOA condition (Fig. 1a). Subjects had a maximum of 5000 ms to evaluate whether the response was valid or invalid as quickly as possible by pressing one of two keys on the computer keyboard.
2.4. fMRI session
During fMRI scanning, children performed a spatial attention localizer task and an arithmetic task. The spatial attention localizer task consisted in alternating blocks of fixation and saccades. During saccade blocks (n = 9), participants were asked to make saccades towards several successive target dots. Each saccade block contained 16 target dots (0.2° visual angle) that appeared at random positions with an eccentricity of 3°, 3.5°, 4°, 4.5°, 5° or 5.5° in the left or right visual field for an average of 800 ms (with a jitter of ±200 ms). During fixation blocks (n = 9), participants were asked to maintain fixation on a central dot for 12,800 ms. Block order was counterbalanced across children.
During the arithmetic task, children were presented with sign-only and sign-plus-operands versions of addition and multiplication trials (Fig. 1b). Each trial started with the presentation of either a ‘+’ or a ‘×’ sign at the center of the screen for 150 ms. In sign-only trials (n = 30), the trial ended with the presentation of the sign and was simply followed by the inter-trial period of fixation (see below). These sign-only trials were our trials of interest and allowed us to isolate neural activity due to the presentation of a sign alone. We also included in the experiment sign-plus-operands trials (n = 50). In those filler trials, the ‘+’ or ‘×’ sign was immediately followed by a single-digit addition or multiplication problem (respectively) presented with an answer. Participants were asked to indicate whether the answer was true or false. The goal of these filler trials (for which associated activity would be difficult to interpret because any effects could be attributable to the anticipatory presentation of the operator, the appearance of the operands, or a combination of both of these factors) was only to keep children engaged and attentive in the scanner. They also induced an arithmetic context, thereby ensuring that the ‘+’ and ‘×’ signs presented in sign-only trials were perceived as arithmetic signs. Problems in sign-plus-operand trials were constructed following the same criteria as in the behavioral session. Finally, the baseline consisted in trials in which the arithmetic sign was replaced by an abstract non-arithmetic sign (i.e., ‘◊’). We included 30 baseline sign-only trials (in which the ‘◊’ sign was presented in isolation) and 50 baseline sign-plus-operand trials (in which the ‘◊’ sign was followed by 3 letters and participants had to indicate whether one of these letters was a B). All trials were followed by a variable period of visual fixation ranging from 3000 ms to 3800 ms. That period consisted in a central white fixation dot that turned red 1000 ms before the onset of the next trial. The arithmetic task was decomposed in 4 functional runs. All trials were intermixed and the timing and order of trial presentation within each run was optimized for estimation efficiency using optseq2 (http://surfer.nmr.mgh.harvard.edu/optseq/). Behavioral responses were recorded using an MR-compatible response device.
Stimuli were generated using Presentation software (Neurobehavioral Systems, Albany, CA). Prior scanning, children were familiarized with the fMRI environment during a practice session that took place after the standardized testing and the behavioral session. During this practice session, children learned to minimize head movement in a mock fMRI scanner. The actual scanning session took place no more than 3 weeks after the practice session.
2.5. Behavioral analyses
RT data associated with the operator-priming task were normalized using a logarithmic transformation prior all analyses to improve the conformity of the data to the standard assumptions of parametric testing. Following Fayol and Thevenot (2012), mean RT was analyzed using planned comparisons that followed from a within-subject ANOVA with the factors Operation (Addition/Multiplication) and SOA (Negative/Null), conducted separately for each group. We report for all effects the corresponding Bayes factors (BF10), indicating the strength of evidence for the alternative hypothesis (H1) relative to the null hypothesis (H0). Substantial evidence in favor of the alternative hypothesis is typically suggested by a BF10 greater than 3 (Jeffreys, 1961, Dienes, 2011).
2.6. fMRI data acquisition
Images were collected with a Siemens Prisma 3T MRI scanner (Siemens Healthcare, Erlangen, Germany) at the CERMEP Imagerie du vivant in Lyon, France. The BOLD signal was measured with a susceptibility weighted single-shot EPI sequence. Imaging parameters were as follows: TR = 2000 ms, TE = 24 ms, flip angle = 80°, matrix size = 128 × 120, field of view = 220 × 206 mm, slice thickness = 3 mm (0.48 mm gap), number of slices = 32. A high-resolution T1-weighted whole-brain anatomical volume was also collected for each participant. Parameters were as follows: TR = 3500 ms, TE = 2.24 ms, flip angle = 8°, matrix size = 256 × 256, field of view = 224 × 224 mm, slice thickness = 0.9 mm, number of slices = 192.
2.7. fMRI preprocessing
Data analysis was performed using SPM12 (http://www.fil.ion.ucl.ac.uk/spm). Functional images were corrected for slice acquisition delays and spatially realigned to the first image of the first run. Images were then spatially smoothed with a Gaussian filter equal to twice the voxel size. ArtRepair was used to help remove motion from the functional images prior to normalization (Mazaika et al., 2009). Volumes with rapid scan-to-scan movements of greater than 1.5 mm were repaired by interpolation of the two nearest non-repaired scans. Each run with more than 5% of the total number of volumes replaced was removed from the analyses. A subject was excluded from further analysis if more than one run was removed. The number of volumes replaced did not differ between grade groups (F(2,31) = 2.20; p = 0.13). Finally, functional images were normalized into the standard MNI space (normalized voxel size, 2 × 2 × 3.5 mm3).
2.8. fMRI processing
Event-related statistical analysis was performed according to the general linear model (GLM). For the localizer task, brain activity associated with saccades and fixation blocks was modeled as epochs with onsets and offsets time-locked to the beginning and the end of each block. Each epoch was convolved with a canonical hemodynamic response function (HRF) and the time series data from each run were high-pass filtered (1/128 Hz). Finally, serial correlations were corrected using an autoregressive AR(1) model. Following our previous study using the same task in adults (Mathieu et al., 2017), brain activity associated with sign-only trials during the arithmetic task was estimated using a finite impulse response (FIR) model. We modeled 8 time points with an interval of 2 s (corresponding to one TR) ranging from the onset of the sign to 16 s after the sign. The magnitude of the fMRI response for each type of sign-only trial was calculated by subtracting activity at the onset of the sign (i.e., 1st bin, or 0 s after the onset) from the peak activity (i.e., 4th bin, or ∼8 s after the onset). The time series data from each run were high-pass filtered (1/128 Hz), and serial correlations were corrected using an autoregressive AR(1) model.
2.9. Region of interest (ROI) definition and analyses
The present study used a Region-of-Interest (ROI) approach to analyze brain activity associated with sign-only trials in brain regions involved in the orienting of spatial attention in children. All ROIs were independently defined using the contrast of saccades versus fixation blocks in the spatial attention localizer task. All subject-specific contrasts were entered into a random effect (RFX) one-sample t-test across subjects. The RFX contrast map was then thresholded across the whole-brain using an uncorrected voxel-level threshold of p < 0.001 and a false-discovery-rate (FDR) corrected cluster-level threshold of p < 0.05 (Chumbley and Friston 2009). Using the SPM toolbox Marsbar (http://marsbar.sourceforge.net/), ROIs were defined as 6-mm radius spheres around the peak coordinate of each region.
Within each ROI and for each participant, we calculated the average response (parameter estimates) for ‘+’ signs using the contrast of addition sign-only trials versus baseline sign-only trials. Similarly, we calculated the average response for ‘×’ signs using the contrast of multiplication sign-only trials versus baseline sign-only trials. Two analyses were performed in each ROI. First, we identified the ROI(s) in which a difference in fMRI activity between ‘+’ and ‘×’ signs emerged throughout age and/or education, using a 3 × 2 ANOVA with the between-subject factor group (lower/intermediate/higher grades) and the within-subject factor sign (‘+’/‘×’). Second, we tested in these ROIs whether inter-individual differences in the fMRI response to an arithmetic sign were correlated with inter-individual differences in the size of the operator-priming effect (measured in the behavioral session) for each group. In all analyses, we report uncorrected P values as well as P values corrected for multiple comparisons across all identified ROIs using the Bonferroni procedure. Bayes factor are also reported.
3. Results
3.1. The spatial localizer task activates a brain network encompassing frontal, parietal, occipital and hippocampal regions
Contrasting saccades to fixation blocks in the spatial attention localizer task, we first identified 10 clusters supporting the orienting of spatial attention across subjects: the bilateral Frontal Eye Field (FEF), bilateral Posterior Superior Parietal Lobule (PSPL), bilateral Middle Temporal Gyri (MTG), bilateral Middle Occipital Gyri (MOG), right dorsal Inferior Frontal Gyrus (dIFG), and right Hippocampus (see Table 1 and Fig. 2). Therefore, a large brain network was involved in the orienting of spatial attention across subjects. Each of these regions served as an ROI in subsequent analyses.
Table 1. Brain regions that were activated during the spatial attention localizer task. Each of these regions constituted an ROI.
Anatomical location
∼BA
Cluster size (mm3)
MNI coordinates
Z-score
X
Y
Z
L. Middle Occipital Gyrus
17
24640
−24
−96
10
5.37
R. Middle Occipital Gyrus/Calcarine
18
–
14
−92
2
4.77
L. Frontal Eye Field
6
5824
−34
−8
55
5.34
R. Middle Temporal Gyrus
21
9534
54
−52
13
4.68
R. Frontal Eye Field
6
5166
28
−4
48
4.65
R. Posterior Superior Parietal Lobule
7
5250
18
−66
66
4.61
L. Posterior Superior Parietal Lobule
7
3178
−22
−64
66
4.58
L. dorsal Inferior Frontal Gyrus
6/44
1316
56
10
20
4.52
L. Middle Temporal Gyrus
21
3626
−50
−54
13
4.33
R. Hippocampus
–
630
24
−16
−18
3.48
L. = left; R. = right; ∼BA = approximate Brodmannʼs area; MNI = Montreal Neurological Institute.
3.2. The right hippocampus region identified by the spatial localizer task is increasingly activated in response to a ‘+’ sign but not to a ‘×’ sign throughout age and/or education
A 3 × 2 ANOVA with the between-subject factor group (lower/intermediate/higher grades) and the within-subject factor sign (‘+’/‘×’) was then conducted in each of the 10 ROIs identified by the spatial attention localizer. We found an interaction between group and sign in the right hippocampus (F(2.30) = 6.75, p= 0.0038, pcorr = 0.038; Fig. 3a and b), but not in any other ROIs (all Fs < 3.44, all ps > 0.046, all pscorr > 0.46). Bayes Factor analysis indicated substantial evidence for this interaction in the right hippocampus (BF10 = 11.53), while no or anecdotal evidence for such an interaction was found in the other ROIs (BF10 < 1.63). Follow-up t-tests in the hippocampus ROI revealed that children from higher grades (average grade = 8.5) exhibited greater activity for addition than multiplication sign-only trials (t11 = 3.02, p = 0.012), whereas there was no difference between signs in children from intermediate grades (average grade = 6.2) (t10 = 0.87, p = 0.41) and even a trend for less activity for addition than multiplication sign-only trials in children from lower grades (average grade = 4.4) (t10 = 2.22, p = 0.051). Bayes Factor analysis indicated substantial evidence for a difference of activity between addition and multiplication sign-only trials in children from higher grades (BF10 = 5.21), but evidence for this difference was absent in intermediate grades (BF10 = 0.41) and anecdotal in lower grades (BF10 = 1.69). Finally, across all groups, addition sign-only trials were associated with greater activity than baseline sign-only trials in children from higher grades (t11 = 2.63, p = 0.023), but not in any other groups (all ts < 1.32, all ps > 0.21). Multiplication sign-only trials were associated with greater activity than baseline sign-only trials in none of the groups (all ts < 1.38, all ps > 0.20). Bayes Factor analysis indicated substantial evidence for a difference between addition sign-only trials and baseline sign-only trials in children from higher grades (BF10 = 3.00), but no or anecdotal evidence in the other groups (all BF10s < 0.60) and for multiplication sign-only trials (all BF10s < 0.63). Overall, then, a difference in fMRI response to a ‘+’ and a ‘×’ sign emerged throughout age and/or education in the right hippocampus.
The findings above were confirmed by an additional correlation analysis in which grade was treated as a continuous predictor across all subjects. The difference in activity between addition and multiplication sign-only trials was positively correlated with grade in the right hippocampus (r= 0.53, p= 0.001, pcorr = 0.01; Fig. 3c). No other significant grade-related changes were found in any other regions (all rs < 0.38, all ps > 0.03, all pscorr > 0.30). Bayes Factor analysis indicated substantial evidence for this correlation in the right hippocampus (BF10 = 20.97), but no or anecdotal evidence in any other ROIs (all BF10s < 2.03). In the right hippocampus, the correlation between grade and the contrast of addition sign-only trials versus baseline sign-only trials, was also near significance (r= 0.33, p= 0.056; BF10 = 1.22).2 The correlation between grade and the contrast of multiplication sign-only trials versus baseline sign-only trials, however, was not significant (r= −0.18, p= 0.32; BF10 = 0.35).
3.3. Spatial hippocampal activity in response to a ‘+’ sign relates to an addition-priming effect in children from higher grades
We then tested whether the hippocampal response to a ‘+’ sign observed in children from higher grades was related to the operator-priming effect. To this aim, each child performed a version of the operator-priming task outside of the scanner (see Fig. 1a).
First, we tested whether the results obtained by Fayol and Thevenot (2012) in adults (i.e., an operator-priming effect for addition but not for multiplication across subjects) could be extended to our children participants. Because children from lower grades had a performance close to chance on large problems (58%), we exclusively focused our analyses on small problems for which accuracy was significantly above chance in all groups (lower grades: 80%, intermediate grades: 92%, higher grades: 96%). Planned comparisons revealed an operator-priming effect for addition in children from higher grades (1491 ms versus 1577 ms; F(1,11) = 8.11, p= 0.016), but not in children from lower grades (2289 ms versus 2417 ms; F(1,10) = 2.66, p= 0.134) and intermediate grades (1530 versus 1509 ms; F(1,10) = 0.01, p= 0.941). No operator-priming effect for multiplication was observed in any groups (lower grades: F(1,10) = 3.50, p= 0.091; intermediate grades: F(1,10) = 1.52, p= 0.246; higher grades: F(1,11) = 0.14, p= 0.715). Bayes Factor analysis indicated substantial evidence for an operator-priming effect with addition problems in children from higher grades (BF10 = 4.08), but no evidence in children from intermediate (BF10 = 0.30) and lower (BF10 = 0.83) grades. There was also no or anecdotal evidence for an operator-priming effect with multiplication problems in any group (higher grades: BF10 = 0.31; intermediate grades: BF10 = 0.55; lower grades: BF10 = 1.09).
Second, we tested whether the size of the operator-priming effect in children (measured on small problems) was correlated to the magnitude of the response for addition sign-only trials versus baseline sign-only trials in the right hippocampus. Such a correlation was found to be highly significant in children from higher grades (r= 0.82, p= 0.0012, Fig. 4), surviving Bonferroni correction for multiple comparisons between the two conditions and across the three groups (pcorr = 0.007). That is, children from higher grades who show greater responses to ‘+’ signs in the right hippocampus are those who show larger operator-priming effect with addition problems. No significant correlation was found in children from lower (r= 0.15, p= 0.66 Fig. 4) and intermediate (r= 0.24, p= 0.48, Fig. 4) grades. Bayes Factor analysis indicated substantial evidence for the correlation in children from higher grades (BF10 = 38.15), but no evidence in children from lower (BF10 = 0.40) and intermediate (BF10 = 0.46) grades. There was also no significant (and anecdotal evidence for a) correlation between the operator-priming effect for addition problems and the fMRI response to multiplication sign-only trials (compared to baseline sign-only trials) in the right hippocampus, in any of the groups (lower grades: r= 0.06, p= 0.87, BF10 = 0.37; intermediate grades: r= 0.32, p= 0.34, BF10 = 0.56; higher grades: r= 0.51, p= 0.09, BF10 = 1.28). Therefore, not only did we observe an operator-priming effect for addition in the only group in which we also observed a greater hippocampal response to ‘+’ than ‘×’ signs (i.e., children from higher grades), but inter-individual differences in the size of the operator-priming effect in that group was also related to hippocampal activity.
Third, we tested whether the correlation between the operator-priming effect and the contrast of addition sign-only trials versus baseline sign-only trials increased over grade. This was done by transforming the correlation coefficient in each group to a Fischer’s z score before comparing the groups using the cocor package (Diedenhofen and Musch, 2015). Although the correlation was not greater in children from intermediate than lower grades (z = 0.19, p = 0.43, one-tailed), it was significantly greater in children from higher than lower grades (z = 2.07, p = 0.019, one-tailed) and in children from higher than intermediate grades (z = 1.88, p = 0.030, one-tailed). Therefore, this brain-behavior correlation increased over grade.
4. Discussion
In the present study, we used fMRI and a cross-sectional design to investigate (i) how and when spatial processing related to the perception of an addition sign emerges in the developing brain, and (ii) to what extent it contributes to the emergence of an operator-priming effect.
4.1. The mere perception of a ‘+’ sign is associated with increased hippocampal spatial activity throughout age and/or education
It has been shown that the processing of a ‘+’ sign is associated with the right side of space (Pinhas et al., 2014) and activates brain regions involved in overt spatial attention in adults (Mathieu et al., 2017). Therefore, we expected arithmetic learning to be associated with increased recruitment of brain regions involved in spatial attention in response to the perception of a ‘+’ sign throughout age and/or education in children. This was the case in a region of the right hippocampus that we identified in our spatial attention localizer task. Therefore, it is possible that hippocampal spatial mechanisms may scaffold the progressive association between an arithmetic operator (i.e., a ‘+’ sign) and spatial intuitions throughout age and/or education. There is increasing evidence that the hippocampal formation, and particularly the right hippocampus, may house a ‘sense of space’ (Buffalo, 2015). Specifically, the right hippocampus has been extensively reported to support spatial representation and navigation in humans (Maguire et al., 1998, Burgess et al., 2002) as well as in non-human primates and rodents (O'keefe and Nadel, 1978, Bird and Burgess, 2008). For example, the hippocampus is typically activated when human participants learn to navigate through a mental representation of space (i.e., mental scanning) (Mellet et al., 2002, Spiers and Maguire, 2006). Interestingly, a recent study in monkeys demonstrated that neurons in the hippocampal formation may encode the direction of overt (Killian et al., 2015) as well as covert (Wilming et al., 2015) shifts of attention. Therefore, the hippocampal formation is likely a critical region for both representing a mental map of space and navigating along that map (Killian et al., 2012, Meister and Buffalo, 2016).
Why would such a hippocampal spatial navigation mechanism be increasingly recruited by the mere perception of a’+’ sign throughout age and/or education? One possibility is that this mechanism might enable children to construct a detailed representation of numbers in mental space, as well as to navigate along that mental representation. Indeed, there is overwhelming evidence that numbers of increasing size are organized along a left-to-right mental map (i.e., the MNL) in adults (Fischer and Shaki 2014). This spatial representation may enable individuals to add or subtract numbers by navigating from a source to a target number to the left or right of that MNL. This is supported by behavioral studies showing that addition and subtraction problem-solving is associated with rightward and leftward shifts of attention (Masson and Pesenti, 2014, Mathieu et al., 2016), as well as by a neuroimaging study indicating an overlap between the brain regions involved in overt shifts of attention and those involved in arithmetic calculation in adults (Knops et al., 2009). Such strategies may be acquired early by children, sometimes even explicitly in the classroom where addition and subtraction is often demonstrated on visual number lines. Yet, it is only with practice that they might become progressively attached to and evoked by an arithmetic operator such as a ‘+’, which might explain the grade-related increases of activity in this region in response to the ‘+’ sign (and the fact that it is only by 7th grade that children exhibit significant activity in response to that sign).
4.2. Hippocampal spatial activity in response to a ‘+’ sign relates to the operator-priming effect in children from higher grades
A critical question is to what extent this automatic processing of a ‘+’ sign in hippocampal spatial mechanisms is associated with children’s behavior. To answer this question, we asked all children to perform a version of the operator-priming task developed by Fayol and Thevenot (2012) and Roussel et al. (2002). First, we replicated the operator-priming effect observed in adults with addition problems (i.e., a facilitation of problem-solving when the operator is presented 150 ms before the operands), but only in children from higher grades (after around 7th grade). Like in adults, this effect was specific to addition problems and not observed with multiplication problems. Thus, the perception of a ‘+’ sign (but not that of a ‘×’ sign) appears to pre-activate a process that is likely used to solve the subsequent problem in children from higher grades. More central to our current interest, we found that the size of the operator-priming effect in these children was highly correlated with the degree of activation of hippocampal spatial mechanisms in response to a ‘+’ sign. This indicates that hippocampal spatial activity may be at the source of the operator priming-effect in older children, perhaps because these children might prepare for an attentional movement along the MNL as soon as a ‘+’ sign is presented. Because no brain-behavior correlation was observed in younger children, extensive practice might be needed before such mechanisms are triggered by the mere perception of the sign.
4.3. Hippocampal spatial activity in response to a ‘+’ sign is transient in development
Strikingly, the spatial brain mechanisms that respond to the mere perception of a ‘+’ sign appear to be different in children and adults. That is, albeit we found increased hippocampal spatial activity throughout age and/or education in the present study, we did not identify these mechanisms in our previous study in adult participants using the exact same task (Mathieu et al., 2017). Rather, we found increased activity in response to a ‘+’ sign in neocortical regions of the FEF and PSPL in adults. Therefore, the contribution of the hippocampus to the automatic processing of a ‘+’ sign is likely transient. Such a transient involvement of the hippocampus is consistent with a wealth of studies that have demonstrated that the spatial representations initially supported by the hippocampus during learning become independent from this brain structure over experience and transferred to neocortical regions (Rosenbaum et al., 2004, Hirshhorn et al., 2012b). For example, longitudinal studies demonstrate that right hippocampal activity associated with learning to mentally navigate through a new environment disappears and is replaced by neocortical activity when individuals become familiar with that environment (Spiers and Maguire, 2007, Hirshhorn et al., 2012a). It is possible that the same phenomenon is at play here: The hippocampus may be involved in the early representation of (and navigation along) the MNL before that representation is transferred to neocortical regions of the fronto-parietal cortex. Future investigations with a wider age sample than in the present study are needed to test this hypothesis.
4.4. Can right hippocampal involvement in the present study reflect mnemonic operations involved in learning arithmetic?
Although there is no doubt that the hippocampus supports spatial processing (Burgess et al., 2002, Spiers and Maguire, 2007), this brain structure is also well known to support the encoding and consolidation of verbal declarative knowledge into long-term memory (Eichenbaum, 2004). In fact, previous developmental studies have largely explained the involvement of the hippocampus during arithmetic learning by referring to its role in declarative memory rather than spatial processing (Rivera et al., 2005, De Smedt et al., 2011, Cho et al., 2011, Cho et al., 2012, Qin et al., 2014). This interpretation relies on the claim that results of well-practiced arithmetic facts (e.g., 2 + 3 or 4 × 2) might become progressively retrieved from memory (rather than calculated) over the course of learning and development (Campbell and Xue, 2001). The hippocampus might thus support the encoding and consolidation of networks of arithmetic facts in children.
Can the role of the hippocampus in declarative memory explain the operator-specific activity over grade (and correlation with the operator-priming effect) observed in the region of the right hippocampus identified by our spatial localizer task? We acknowledge that we did not have a task identifying processes involved in declarative memory. Thus, even if the right hippocampus is usually more associated with spatial than mnemonic processes (Burgess et al., 2002), it is possible that the hippocampal cluster that we identified as being involved in spatial processing may also be involved in some aspects of declarative memory. One might thus argue that grade-related increases of activity in relation to ‘+’ signs reflect the progressive association between a ‘+’ and a network of additive facts. This explanation, however, can be ruled out by an examination of activity related to ‘×’ signs. Because single-digit multiplication problems are almost exclusively learned by rote in school, multiplication is the operation that is perhaps the most associated with a network of stored facts in the literature (Campbell and Xue, 2001). Thus, if increased hippocampal activity in relation to ‘+’ signs were due to the progressive building of a network of additive facts, increased activity in that same region should have been observed during the perception of ‘×’ signs (perhaps even more so for the perception of ‘+’ signs). Yet, this is not the case. Not only did we not find any grade-related increase of activity for ‘×’ signs in the hippocampal cluster identified by our spatial localizer task, but activity was significantly greater for ‘+’ than ‘×’ signs in higher graders (who are as proficient in single-digit multiplication as addition). Similarly, no operator-priming effect was observed for multiplication problems in higher graders, indicating that the operator-priming effect observed for addition is likely to have little to do with the pre-activation of a network of stored facts (because this should be also observed for multiplication). Therefore, the specificity of our results to ‘+’ signs (as compared to ‘×’ signs) in the right hippocampus ROI makes it very unlikely that our results are related to mnemonic operations. In our view, emerging associations between ‘+’ signs and spatial intuitions related to the MNL are the best explanation of the effects reported here.
Of course, the fact that the role of the hippocampus in declarative memory is unlikely to explain our operator-specific findings in the right hippocampus ROI does not mean that hippocampal mechanisms supporting mnemonic operations do not contribute to arithmetic learning. Instead, they indicate that the hippocampus might contribute to arithmetic learning through its role in both declarative memory and spatial processing. Interestingly, the operator-specific activity observed in our (spatially localized) right hippocampal cluster is not observed in a mirror (left lateralized) cluster that is not activated in the localizer contrast (see Supplementary information). In that mirror region, no difference was observed between activity related to ‘+’ and ‘×’ signs in any group of children (and left hippocampal activity was not related to the operator-priming effect). Thus, the developmental effect reported here appears to be restricted to the right hippocampus. This specificity suggests that the observed developmental changes in the right hippocampus may not simply reflect general brain maturation but rather mechanisms that are specific to arithmetic learning.
4.5. Limitations
It is worth acknowledging here 2 potential limitations of the present work. First, as is the case for any cross-sectional fMRI studies, our study is correlational in nature. Thus, although our findings are consistent with the idea that the right hippocampus might scaffold the progressive association between (at least some) arithmetic operators and space throughout age and/or education, future studies might specifically investigate the causal role of these hippocampal mechanisms. Second, our finding of a correlation between grade and the processing of an addition sign in the right hippocampus (see Fig. 3C) relies on a relatively large sample size of 34 children. However, other findings involve subgroups of participants and therefore rely on smaller sample sizes. In particular, null findings in relation to these subgroups might be difficult to interpret because of potential lack of power. For example, whereas we found an operator-priming effect in children from higher grades and no effect in children from intermediate grades, there was no significant difference between these groups in terms of response times in negative SOA trials (1491 ms versus 1530 ms; t21 = 0.21; p = 0.84; BF10 = 0.39). Behavioral studies focusing on the operator-priming effect in children might test whether this difference emerges with larger sample sizes. More generally, future studies are needed to improve our understanding of the present results.
5. Conclusion
In sum, our findings suggest that the right hippocampus might contribute to the progressive association between (at least some) arithmetic operators and space throughout age and/or education. Therefore, our study raises the possibility that increased hippocampal activity during arithmetic learning in children may be explained by the role of this structure in spatial representations as well as in declarative memory.
Conflict of interest
The authors declare no competing financial interests.
Acknowledgments
This research was supported by a grant from the European Union (Marie Curie Career Integration Grant n° PCIG12-GA-2012-333602) to J.P. and a grant from the French Ministry of Higher Education and Research to R.M. We thank the Hospices Civils de Lyon for sponsoring the research, as well as Flora Schwartz and the MRI engineers (Franck Lamberton and Danielle Ibarrola) at the CERMEP-Lyon platform for their assistance in collecting the fMRI data. Finally, we are grateful to Pr. Christian Scheiber for his help with the pre-MRI medical exams.
Appendix A. Supplementary data
The following is Supplementary data to this article:
S. Cho, A.W. Metcalfe, C.B. Young, S. Ryali, D.C. Geary, V. MenonHippocampal–prefrontal engagement and dynamic causal interactions in the maturation of children's fact retrieval
B. De Smedt, I.D. Holloway, D. AnsariEffects of problem size and arithmetic operation on brain activation during calculation in children with varying levels of arithmetical fluency
M. Hirshhorn, C. Grady, R.S. Rosenbaum, G. Winocur, M. MoscovitchThe hippocampus is involved in mental navigation for a recently learned, but not a highly familiar environment: a longitudinal fMRI study
M. Hirshhorn, C. Grady, R.S. Rosenbaum, G. Winocur, M. MoscovitchBrain regions involved in the retrieval of spatial and episodic details associated with a familiar environment: an fMRI study
I.D. Holloway, D. AnsariMapping numerical magnitudes onto symbols: the numerical distance effect and individual differences in children’s mathematics achievement
R. Mathieu, J. Epinat-Duclos, M. Sigovan, A. Breton, A. Cheylus, M. Fayol, C. Thevenot, J. PradoWhat’s behind a ‘ + ’ sign? Perceiving an arithmetic operator recruits brain circuits for spatial orienting
E. Mellet, S. Bricogne, F. Crivello, B. Mazoyer, M. Denis, N. Tzourio-MazoyerNeural basis of mental scanning of a topographic representation built from a text
S. Qin, S. Cho, T. Chen, M. Rosenberg-Lee, D.C. Geary, V. MenonHippocampal-neocortical functional reorganization underlies children's cognitive development
S.M. Rivera, A. Reiss, M.A. Eckert, V. MenonDevelopmental changes in mental arithmetic: evidence for increased functional specialization in the left inferior parietal cortex
R.S. Rosenbaum, M. Ziegler, G. Winocur, C.L. Grady, M. MoscovitchI have often walked down this street before: fMRI studies on the hippocampus and other structures during mental navigation of an old environment
J.-L. Roussel, M. Fayol, P. BarrouilletProcedural vs. direct retrieval strategies in arithmetic: a comparison between additive and multiplicative problem solving
Note that to induce an arithmetic context and disguise the goal of the experiment, we also included trials in which a ‘ + ’ or a ‘ × ’ sign was followed 150 ms later by operands and participants were asked to solve the problem. The low temporal resolution of fMRI, however, makes it impossible to dissociate activity associated with the sign from activity associated with operands in these problems. Therefore, they were simply designed to be filler trials.
There was a tendency for a correlation between grade and activity associated with addition sign-only trials (versus fixation) in the right hippocampus (r = 0.29, p = 0.09; BF10 = 0.82), but no correlation for baseline sign-only trials (versus fixation) (r = −0.10, p = 0.58; BF10 = 0.25). Thus, the correlation between grade and the contrast of addition sign-only trials versus baseline sign-only trials was more likely driven by changes of activity in addition sign-only trials than in baseline sign-only trials.
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/identifiers.json b/tests/data/sample_inputs/8EVW7TUtC9cx/identifiers.json
new file mode 100644
index 0000000..4b13178
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/identifiers.json
@@ -0,0 +1 @@
+{"pmid": "26878057", "doi": "10.1523/ENEURO.0107-15.2016", "pmcid": "PMC4745181"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/coordinates.csv b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/coordinates.csv
new file mode 100644
index 0000000..ba82d07
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/coordinates.csv
@@ -0,0 +1,51 @@
+table_id,table_label,table_caption,table_number,x,y,z,p_value,region,size,statistic,groups
+T3,Table 3:,,,-16.0,-96.0,24.0,,,,,
+T3,Table 3:,,,-30.0,-76.0,34.0,,,,,
+T3,Table 3:,,,-52.0,22.0,-10.0,,,,,
+T3,Table 3:,,,18.0,-98.0,16.0,,,,,
+T3,Table 3:,,,8.0,-82.0,38.0,,,,,
+T3,Table 3:,,,-30.0,-72.0,34.0,,,,,
+T3,Table 3:,,,-32.0,-32.0,6.0,,,,,
+T3,Table 3:,,,28.0,-74.0,46.0,,,,,
+T3,Table 3:,,,-24.0,-66.0,-38.0,,,,,
+T3,Table 3:,,,30.0,58.0,12.0,,,,,
+T3,Table 3:,,,54.0,-42.0,54.0,,,,,
+T3,Table 3:,,,10.0,10.0,2.0,,,,,
+T3,Table 3:,,,-34.0,6.0,-10.0,,,,,
+T3,Table 3:,,,18.0,-76.0,-38.0,,,,,
+T3,Table 3:,,,52.0,-40.0,54.0,,,,,
+T3,Table 3:,,,6.0,22.0,66.0,,,,,
+T3,Table 3:,,,6.0,6.0,4.0,,,,,
+T3,Table 3:,,,30.0,58.0,14.0,,,,,
+T3,Table 3:,,,-30.0,58.0,10.0,,,,,
+T3,Table 3:,,,32.0,58.0,14.0,,,,,
+T3,Table 3:,,,-34.0,52.0,4.0,,,,,
+T3,Table 3:,,,30.0,58.0,12.0,,,,,
+T3,Table 3:,,,-24.0,56.0,6.0,,,,,
+T3,Table 3:,,,32.0,56.0,10.0,,,,,
+T3,Table 3:,,,-30.0,58.0,8.0,,,,,
+T3,Table 3:,,,-16.0,-96.0,24.0,,,,,
+T3,Table 3:,,,-30.0,-76.0,34.0,,,,,
+T3,Table 3:,,,-52.0,22.0,-10.0,,,,,
+T3,Table 3:,,,18.0,-98.0,16.0,,,,,
+T3,Table 3:,,,8.0,-82.0,38.0,,,,,
+T3,Table 3:,,,-30.0,-72.0,34.0,,,,,
+T3,Table 3:,,,-32.0,-32.0,6.0,,,,,
+T3,Table 3:,,,28.0,-74.0,46.0,,,,,
+T3,Table 3:,,,-24.0,-66.0,-38.0,,,,,
+T3,Table 3:,,,30.0,58.0,12.0,,,,,
+T3,Table 3:,,,54.0,-42.0,54.0,,,,,
+T3,Table 3:,,,10.0,10.0,2.0,,,,,
+T3,Table 3:,,,-34.0,6.0,-10.0,,,,,
+T3,Table 3:,,,18.0,-76.0,-38.0,,,,,
+T3,Table 3:,,,52.0,-40.0,54.0,,,,,
+T3,Table 3:,,,6.0,22.0,66.0,,,,,
+T3,Table 3:,,,6.0,6.0,4.0,,,,,
+T3,Table 3:,,,30.0,58.0,14.0,,,,,
+T3,Table 3:,,,-30.0,58.0,10.0,,,,,
+T3,Table 3:,,,32.0,58.0,14.0,,,,,
+T3,Table 3:,,,-34.0,52.0,4.0,,,,,
+T3,Table 3:,,,30.0,58.0,12.0,,,,,
+T3,Table 3:,,,-24.0,56.0,6.0,,,,,
+T3,Table 3:,,,32.0,56.0,10.0,,,,,
+T3,Table 3:,,,-30.0,58.0,8.0,,,,,
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/metadata.json b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/metadata.json
new file mode 100644
index 0000000..cdf6c67
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/metadata.json
@@ -0,0 +1,11 @@
+{
+ "title": "Testosterone Modulates Altered Prefrontal Control of Emotional Actions in Psychopathic Offenders123",
+ "authors": "Volman, Inge; von Borries, Anna Katinka Louise; Bulten, Berend Hendrik; Verkes, Robbert Jan; Toni, Ivan; Roelofs, Karin; Volman, Inge; von Borries, Anna Katinka Louise; Bulten, Berend Hendrik; Verkes, Robbert Jan; Toni, Ivan; Roelofs, Karin",
+ "journal": "eNeuro",
+ "keywords": "amygdala\nconnectivity\nemotion\nfMRI\nprefrontal\npsychopathy\n",
+ "abstract": " \nPsychopathic individuals are notorious for their controlled goal-directed aggressive behavior. Yet, during social challenges, they often show uncontrolled emotional behavior. Healthy individuals can control their social emotional behavior through anterior prefrontal cortex (aPFC) downregulation of neural activity in the amygdala, with testosterone modulating aPFC\u2013amygdala coupling. This study tests whether individual differences in this neuroendocrine system relate to the paradoxical lack of emotional control observed in human psychopathic offenders. Emotional control was operationalized with an fMRI-adapted approach\u2013avoidance task requiring rule-driven control over rapid emotional responses. Fifteen psychopathic offenders and 19 matched healthy control subjects made approaching and avoiding movements in response to emotional faces. Control of social emotional behavior was required during affect-incongruent trials, when participants had to override affect-congruent, automatic action tendencies and select the opposite response. Psychopathic offenders showed less control-related aPFC activity and aPFC\u2013amygdala coupling during trials requiring control of emotional actions, when compared with healthy control subjects. This pattern was particularly pronounced in psychopathic individuals with high endogenous testosterone levels. These findings suggest that reduced prefrontal coordination underlies reduced behavioral control in psychopathic offenders during emotionally provoking situations. Even though the modest sample size warrants replication, the modulatory role of endogenous testosterone on the aPFC\u2013amygdala circuit suggests a neurobiological substrate of individual differences that is relevant for the advancement of treatment and the reduction of recidivism.\n \n ",
+ "publication_year": 2016,
+ "coordinate_space": "MNI",
+ "license": "http://creativecommons.org/licenses/by/4.0/",
+ "text": true
+}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/text.txt b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/text.txt
new file mode 100644
index 0000000..b96f0ed
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/processed/pubget/text.txt
@@ -0,0 +1,154 @@
+
+## Significance Statement
+
+Psychopathic criminals are commonly seen as instrumentally abusive and emotionally callous, yet social challenges often trigger uncontrolled emotional behavior in those individuals. This study shows how this paradoxical aspect of psychopathy relates to altered neuroendocrine interactions between testosterone and the cerebral circuit coordinating emotional action tendencies. The anterior prefrontal cortex, a region necessary for controlling emotional behavior, showed blunted responses and reduced connectivity with the amygdala in psychopathic criminals engaged in controlling their emotional action tendencies. This cerebral pattern was strongest in psychopathic individuals with high endogenous testosterone levels. This neuroendocrine signature of altered emotional control highlights the relevance of considering the testosterone level of individual psychopathic patients during treatment of their impulsive behavior.
+
+
+## Introduction
+
+Psychopathy is a disorder often associated with blunted emotional responding and increased goal-directed behavior ( ; ). On the other hand, offenders with psychopathy also show a paradoxical increase in impulsive behavior and uncontrolled aggression after emotional provocations ( ; ; ; ; ; ), which may be related to heightened testosterone levels ( ; ). These two aspects of psychopathy are also distinguished within the most commonly used psychopathy checklist, the Psychopathy Check List-Revised (PCL-R), potentially reflecting differing traits among psychopathic individuals ( ; ). Importantly, enhanced difficulty in controlling emotional impulses, a crucial component of criminal psychopathy associated with PCL-R factor 2, has been largely neglected by cognitive neuroscience. Yet, the clinical relevance of this cognitive trait is large: reduced behavioral control and increased impulsivity predict recidivism in psychopathic offenders ( ), and behavioral control in psychopathic offenders appears particularly fragile when dealing with emotionally relevant behavior ( ; , chapter 7; ). Accordingly, understanding the neurobiological systems underlying the altered control of social emotional behavior in psychopathic individuals is relevant for improving currently available interventions, which are plagued by low treatment response and high recidivism ( ). Here we study those neuroendocrine systems in a group of psychopathic offenders engaged in an experimental paradigm that requires rule-driven control of emotional behavior.
+
+Previous investigations of psychopathy showed altered reactivity to emotional material in several brain regions that include the anterior part of the PFC (aPFC) and the amygdala ( ; ; ). Furthermore, individuals with psychopathy showed decreased functional and anatomical connectivity between the PFC and amygdala at rest ( ; ), an indication that these brain regions might have a reduced ability to interact effectively. Studies in healthy participants have shown that this cerebral circuit is necessary for implementing the control of emotionally relevant actions ( ). Namely, aPFC downregulates neural processing in the amygdala during emotional control ( ), while high levels of endogenous testosterone reduce such control-related connectivity between aPFC and amygdala ( ). Those findings raise the possibility that aPFC–amygdala connectivity is altered when psychopathic offenders need to control emotionally relevant actions, with high levels of endogenous testosterone exacerbating that altered connectivity.
+
+This study tests these hypotheses by measuring brain activity with functional magnetic resonance imaging (fMRI) in 15 psychopathic criminals and 19 matched healthy control subjects dealing with a challenge to control their emotional behavior. The psychopathy sample was obtained by focused and comprehensive screening excluding confounds that are frequently associated with random criminal sampling (e.g., medication use, comorbidity). The social approach–avoidance (AA) task was used to provide reliable indexes of control over social emotional behavior ( ; ; , ). Behaviorally, psychopathic participants previously showed altered AA behavior to explicitly approaching and avoiding emotional faces ( ). Similar findings occurred after testosterone administration in healthy participants ( ). Interestingly, a more subtle version of the AA task has been shown to be sensitive to testosterone-related alterations and genetic variations in the aPFC–amygdala pathway, while keeping behavior constant across experimental groups ( ), opening the way for isolating neural vulnerability factors ( ) in psychopathy. During this task, participants respond to affective faces (happy, angry) presented for a short time with approach and avoidance movements. Automatic emotional tendencies (approach–happy and avoid–angry faces; affect-congruent response conditions) need to be controlled during affect-incongruent response conditions in order to apply the counterintuitive action of approaching angry and avoiding happy faces ( ; ). Healthy participants respond more slowly and rely more strongly on the aPFC when emotional control is required, operationalized by the differences evoked between affect-incongruent and affect-congruent trials ( ; ). Accordingly, this study tests whether exerting control over emotionally relevant actions is reflected by reduced functionality of the aPFC–amygdala circuit in psychopathic individuals, suggesting less prefrontal regulation of emotional actions. In addition, it sets out to test whether this alteration is intensified by high levels of endogenous testosterone.
+
+The emotional control AA task. The AA task involved the presentation of happy and angry faces, and the performance of approach and avoidance responses. During the AA task, the participants had to select their response according to the perceived emotion of the face. At the beginning of each block of 12 trials, the participants received instructions on whether to pull the joystick toward themselves (approach) or push it away (avoid) when seeing a face with a particular emotion. When viewing happy or angry faces, automatic stimulus–response tendencies trigger corresponding approach or avoidance actions. These tendencies could be followed during the affect-congruent condition (approach–happy, avoid–angry). In contrast, when task instructions required participants to avoid happy faces or to approach angry faces, automatic tendencies needed to be controlled and overridden with the instructed response (affect-incongruent condition). Participants saw the faces and moved the joystick while lying in a MR scanner (top left corner of the table). Figure adapted from ).
+
+
+## Materials and Methods
+
+### Participants
+
+The psychopathic group was recruited from in-patient populations of the Pompestichting and Oldenkotte, forensic psychiatric institutes (TBS-clinics) in the Netherlands. TBS-clinics are facilities for criminal offenders with a mental disorder treated on behalf of the state.
+
+Seventeen male psychopathic violent offenders (age range, 23-56 years) participated; all had received a diagnosis with a PCL-R score of ≥26, according to European standards ( ; ; ). PCL-R consensus scores were obtained by trained clinicians based on a structured PCL-R interview, clinical status, and history. After the independent scoring, the two raters compared their scores and came to the consensus score. When no consensus could be found, a third independent rater was included in the process. Dutch versions of the National Adult Reading Test and Edinburgh Handedness Inventory were used to assess IQ levels and right-handedness ( ; ). Twenty-one healthy male control subjects (HCs) matched for age, right-handedness, and IQ, without criminal records or history of psychiatric disorders, were recruited from staff of the clinics. All participants received oral and written information about the experiment and gave written informed consent according to guidelines of the local ethics committee (Commissie Mensengebonden Onderzoek region Arnhem-Nijmegen). Psychiatric exclusion criteria consisted of neurological, axis-I, and axis-II disorders, besides antisocial personality disorder for the psychopathic group. They were screened for these exclusion criteria by trained psychologists using Dutch versions of the Structured Clinical Interview (SCID; ) and Mini-International Neuropsychiatric Interview (MINI; ) for Diagnostic and Statistical Manual of Mental Disorders , 4th edition, disorders. All participants were asked about drug use and medical/neurological history to exclude the following: alcohol use of >3 units/day, cannabis, or other illicit drug use 1 week before, psychotropic medication other than oxazepam 5 d before, 1 unit of alcohol or oxazepam use within 24 h before the experiment; history of trauma capitis; visual and auditive disorder; and neurological disorder. Furthermore, general exclusion criteria for MRI experiments were applied. Two psychopathic patients (PPs) and two HCs were excluded from the analyses, due to incomplete scanning procedures (1 PP, 1 HC) or too many errors on the task (>16%, representing the outlier with a z -score >3). The final groups did not differ in age, IQ, and handedness (see ).
+
+Demographical data
+
+
+### Procedure
+
+Two test sessions took place. During the first session, right-handedness, IQ, MINI, and SCID were assessed. During the second session, participants completed several questionnaires upon arrival in the laboratory, including the State-Trait Anxiety Inventory (STAI) to measure anxiety levels ( ). Next, they provided saliva for the testosterone measurement. Afterward, participants were positioned in the 1.5 T MR scanner and familiarized with the task setup. Immediately after this, the fMRI session started with the AA task (duration, 30 min) followed by another task (not included in this report). After a short break outside the scanner, the anatomical scan (duration, 5 min) and an unrelated task were acquired in the side-by-side 3 T MR scanner.
+
+
+### Experimental task
+
+The AA task consisted of 24 blocks (with 12 trials per block and a baseline period of 21-24 s) during which participants had to respond to visually presented faces either by pulling a joystick toward themselves (approach) or by pushing it away from themselves (avoid; ). The participants had to categorize faces as happy, angry, and neutral (filler items), based on their affective expressions. During each block, two of the three affective expressions were presented as stimuli, because only two responses could be given to categorize the stimulus. This resulted in six different block types each used four times, representing the affect (happy–angry, happy–neutral, angry–neutral) × movement (approach–avoid) combinations. At the start of each block, participants received written instructions regarding the required response mapping. The affect × movement combinations were pseudorandomly and evenly distributed (with no affect combination repetition), and the combination of the first block was counterbalanced across participants. Within each block, affective expressions and gender types were pseudorandomly presented, avoiding three or more sequential presentations of the same expression/gender, and two presentations of the same facial model. Each face was presented for 100 ms, preceded by a 300 ms blank screen, and followed by the participant’s response, a blank screen, and by a pseudorandom intertrial interval (ITI; 1-3 s). A baseline period of 21-24 s preceded each block. The faces were from 36 models (18 male) obtained from several databases ( ; ; ; ), each showing all expressions. The pictures were in grayscale, matched for brightness and contrast values, and displayed against a black background. To exclude influence from hair and nonfacial contours, the faces were trimmed. Joystick displacements of >80% along the sagittal plane within 2 s from stimulus presentation were marked as valid responses. Invalid responses were signaled for 1 s with written feedback stating “you did not move your joystick far enough.” After moving the joystick, participants had to return to the starting position (defined as the central area extending 20% along the sagittal plane) before the end of the ITI. Otherwise, visual feedback indicated “return the joystick to the starting position,” and the ITI was repeated after participants returned the joystick. The training at the beginning consisted of six blocks; one block of eight trials for each of the six affect × movement combinations. Different visual stimuli were used during the training and scanning blocks.
+
+
+### Materials and apparatus
+
+The fMR images were acquired on a 1.5 T MRI scanner (Avanto, Siemens Medical Systems) with an eight-channel head coil using a multiecho generalized autocalibrating partially parallel acquisitions (GRAPPA) sequence [ ; repetition time (TR), 2.14 ms; five echo times (TEs), 9.4/21/33/44/56 ms; 34 transversal slices; ascending acquisition; distance factor, 17%; effective voxel size, 3.3 × 3.3 × 3.5 mm; field of view (FOV), 212 mm]. High-resolution anatomical images were acquired on a 3 T MRI scanner with a 32-channel head coil using a magnetization prepared rapid gradient echo sequence (TR, 2300 ms; TE, 3.03 ms; 192 sagittal slices; voxel size, 1.0 × 1.0 × 1.0 mm; FOV, 256 mm).
+
+An MR-compatible joystick (Fiber Optic Joystick, Current Designs; sampling rate, 550 Hz) was placed on participants’ abdomens to ensure comfortable push-and-pull movements ( ). Participants wore MR-compatible headphones to reduce scanner noise (Commander XG MRI Audio System, Resonance Technologies). Stimuli were projected at the center of a screen, viewed via a mirror above the participant’s head, with a visual angle of 4° × 6° (width × height). Stimuli presentation and acquisition of joystick positions were controlled by a PC running Presentation version 13 ( ).
+
+
+### Salivary measurements
+
+Participants filled two Salicaps (IBL) with saliva for testosterone measurement, which were stored at −25°C. Testosterone concentration was measured using competitive chemiluminescence immunoassay with a sensitivity of 0.0025 ng/ml (IBL International, Tecan). Intra-assay and interassay coefficients are between 10% and 12%. To control variables influencing testosterone levels, participants were instructed to refrain from any food, cigarettes, and drinks (except water) for 1 h before the experiment.
+
+
+### Behavioral analysis
+
+Behavioral data was analyzed using MATLAB version 7.9 (MathWorks) and PASW Statistics 18 (SPSS Inc.). First, to obtain a precise measure of movement onset [reaction time (RT)], the joystick movement for each trial was reconstructed using the joystick displacement measurements. Excluded trials showed a joystick movement in the wrong direction, an extreme RT (<150 or >1500 ms), peak velocity (<0.1 cm/s), or movement time (>400 ms); or an error rate of above chance level in a block (in that case, the whole block was excluded). RTs and testosterone levels were log transformed to obtain a normal distribution. Second, following previous studies ( ; ), we conducted three-way repeated-measures ANOVA (ANCOVArm) on the mean RT and error rates, with factors group (PP, HC), movement (approach, avoid), and valence (happy, angry), including standardized testosterone and STAI state as covariate. A measure of anxiety (STAI) was included to account for the effects of psychopathy type (e.g., primary vs secondary); and the possible effects on emotional behavior, hormonal levels, amygdala, and prefrontal cortex functioning ( ; ; ; ). The α-level was set at p < 0.05.
+
+
+### Functional MRI data
+
+#### Single-subject analyses
+
+Imaging data were preprocessed and analyzed using SPM8 (Statistical Parametric Mapping; ). The first four volumes of each participant’s dataset were discarded to allow for T equilibration. Given the multiecho GRAPPA MR sequence (Poser et al., 2006), head motion parameters were estimated on MR images with the shortest TE (9.4 ms), since these are least affected by possible artifacts. These motion correction parameters, estimated using a least-squares approach with six rigid body transformation parameters (translations, rotations), were applied to the five echo images collected for each excitation. After spatial realignment, the five echo images were combined into a single MR volume using an optimized echo weighting method (Poser et al., 2006). The time series for each voxel was temporally realigned to the first slice in time. The T -weighted image was spatially coregistered to the mean of the functional images. The fMRI time series were transformed and resampled at an isotropic voxel size of 2 mm into standard Montreal Neurological Institute (MNI) space by unified segmentation and normalization using the coregistered T -weighted image ( ). The normalized functional images were spatially smoothed using an isotropic 8 mm full-width at half-maximum Gaussian kernel.
+
+The fMRI time series of each subject were further analyzed using an event-related approach in the context of general linear model, including the following effects: approach–happy, approach–neutral, approach–angry, avoid–happy, avoid–neutral, and avoid–angry. Trials excluded from behavioral analyses and periods of instructions or feedback were modeled as regressors. Vectors describing the time of picture presentation (onset) and RT of each event (duration) were convolved with the canonical hemodynamic response function. Potential confounding effects of residual head movement were modeled using original, squared, cubic, first-order, and second-order derivatives of the movement correction parameters ( ). Three further regressors, describing the time course of signal intensities of white matter, CSF, and the portion of the MR image outside the skull were also added. This procedure accounts for image intensity shifts due to hand movements within or near the magnetic field of the scanner ( ). Finally, fMRI time series were high-pass filtered (cutoff 120 s). Temporal autocorrelation was modeled as a first-order autoregressive process.
+
+
+#### Group analyses
+
+Consistent effects across participants and between groups were tested using a random-effects multiple regression analysis that included six contrast images (approach–happy, approach–neutral, approach–angry, avoid–happy, avoid–neutral, avoid–angry) per participant. Together, these images represented the estimated cerebral effects from 12 conditions of the experimental design [group (PP, HC) × valence (happy, neutral, angry) × response (approach, avoid)]. Standardized log-transformed testosterone and standardized STAI state levels were included in the multiple regression analysis as condition-specific [group (PP, HC) × valence (happy, neutral, angry) × response (approach, avoid)] regressors, generating another 12 regressors per variable.
+
+All analyses assessed the congruency effect, reflecting task-related differences of affect-incongruent (approach–angry, avoid–happy) versus affect-congruent trials (approach–happy, avoid–angry; ; ). We considered two effects. First, to test for general effects of congruency, we performed an analysis on the congruency effect over both groups and for each group separately. When assessing the effects of one group explicitly, we also tested whether those effects were specific to that group and were significantly weaker in the other group (at p < 0.05 uncorrected) by masking the statistical map describing the congruency effect in the first group (using multiple comparisons correction, see below) with the statistical map describing the group × congruency contrast. Second, to test whether testosterone differentially modulated the control of emotionally relevant actions in the groups, we performed a group × congruency contrast on the regressor parametrizing interindividual differences in testosterone on task-related conditions. If such an interaction is present, the testosterone modulation on the congruency effect of each group separately is considered. In addition to whole-brain analyses, we used a volume of interest (VOI) on coordinates previously found to be modulated by testosterone during the congruency effect in healthy students (two 8-mm-radius spheres centered on the following MNI coordinates: x , −30; y , 58; and z, 2; and x , 32; y , 54; and z , 8; ).
+
+The reported activations are corrected for multiple comparisons using familywise error (FWE) correction. For whole-brain analyses, we made inferences at cluster level (FWE: p < 0.05, corresponding to a cluster size of >140 on the basis of intensity threshold, p < 0.001). For VOI analyses, we made inferences at voxel-level (FWE corrected, p < 0.05; ; ). Anatomical inference is drawn by superimposing SPM showing significant signal changes on structural images of participants. For anatomical accuracy, we report only activation peaks in gray matter. Anatomical landmarks were identified using the atlas of . Brodmann areas (BAs) were assigned by superimposing significant SPM on the SPM anatomy toolbox ( ) and MRIcron template ( /).
+
+
+
+### Connectivity analyses
+
+The aim of the following analysis was to test whether inter-regional coupling of the aPFC (see Results) with the amygdala and other brain regions during the congruency effect was different between the groups and modulated by testosterone. To test for these effects, we used the psychophysiological interactions (PPIs) method ( ). More specifically, we tested for significant differences between the regression coefficients of each voxel over the right aPFC during the affect-incongruent versus the affect-congruent conditions. To select voxels to be included in the VOI, we used the following anatomical constraints ( ): for each participant, selected voxels fell within a sphere with a radius of 4 mm around the peak voxel corresponding to the activated cluster of the congruency effect over both groups (coordinates: x , 30; y , 58; z , 14; see Results). Participant specific contrast images were generated describing the PPI between the time courses of the right aPFC VOI and affect-incongruent versus affect-congruent conditions. Group differences and testosterone modulations on task-related coupling between the aPFC and other regions were then assessed using a multiple regression design on participant-specific contrast images with their corresponding testosterone (log-transformed, standardized) and STAI state (standardized) levels as subject- and group-specific regressors. In addition to whole-brain analyses, we assessed significant voxel-level effects (FWE corrected for multiple comparisons, p < 0.05) within the amygdala, defined on the Automated Anatomical Labeling atlas ( ) using the WFU PickAtlas tool ( ).
+
+
+
+## Results
+
+### Behavioral results
+
+Fifteen psychopathic criminals (PPs; PCL-R score of ≥26, according to European standards ( ; ; ) and 19 HCs (for demographics, see ) were included in the analyses. Participants performed the task accurately and consistently (error rates: PPs, 7.9%; HCs, 7.3%; omissions: PPs, 1.6%; HCs, 1.5%; undefined responses: PPs, 0.9%; HCs, 0.3%; ).
+
+RTs and error rates for each group and factor of the AA task
+
+A significant movement × valence interaction for the RTs indicated that, over groups, participants responded more slowly during affect-incongruent (approach–angry, avoid–happy) than during affect-congruent trials (approach–happy, avoid–angry; F = 10.4, p = 0.003; ). This congruency effect replicates the behavioral results from previous fMRI studies ( ; ). Furthermore, there were main effects of movement ( F = 26.3, p < 0.001) and valence ( F = 28.7, p < 0.001), reflecting the slowing of avoidance movements and responses to angry faces in general ( ). There were no significant effects involving group, including no main effect ( p > 0.3). The congruency effect correlated positively (without corrections for multiple comparisons) with the PCL-R total score ( p = 0.048, R = 0.517, respectively). Excluding anxiety from the analyses did not affect the outcomes. Moreover, when including the neutral conditions in the analyses, the movement × valence (happy, neutral, angry) interaction for RTs remained significant ( F = 5.5, p = 0.010), showing that neutral approach–avoidance effects are intermediary compared with happy and angry ( ).
+
+Behavioral results. Mean RTs (±SEM) for the affect-congruent and affect-incongruent conditions of the AA task for the healthy control subjects and psychopathic offenders. The groups were significantly slower to provide affect-incongruent responses (approach–angry; avoid–happy) than affect-congruent responses (approach–happy; avoid–angry), with no significant group differences.
+
+For the error rates, the three-way ANCOVArm showed main effects of movement ( F = 27.5, p < 0.001), valence ( F = 25.9, p < 0.001), and testosterone ( F = 4.6, p = 0.040), and a valence × testosterone interaction ( F = 4.3, p = 0.047). There were no other significant effects for the error rates ( p > 0.15).
+
+Endogenous testosterone levels [median (SD): PPs, 101 pg/ml (70 pg/ml); HCs, 90 pg/ml (46 pg/ml)] and state anxiety levels [STAI mean (SD): PPs, 32 (8); HCs, 32 (5)] did not differ between groups ( p > 0.4), and showed no correlations with psychopathy (PCL-R) scores or with each other ( p > 0.1).
+
+
+### fMRI results
+
+#### Multiple regression analyses
+
+To assess the two main questions of this study, we isolated cerebral structures showing stronger responses during affect-incongruent than affect-congruent trials (congruency effect), and cerebral structures in which the congruency effect was modulated by testosterone levels.
+
+The results showed a significant congruency effect across groups in the aPFC [ROI analysis: MNI coordinates ( x , y , z ): (30, 58, 14) and (−30 58 10); p = 0.001 and 0.036; t = 4.46 and 3.43; for further details, see ]. As expected, this effect was driven by the healthy control group, and it was significantly weaker in the psychopathic offenders [ p = 0.001 and 0.040; t = 4.58 and 3.40, on the congruency effect in healthy control subjects masked implicitly by group (HC > PP) × congruency interaction]. The implicit masking demonstrates that the group × congruency interaction is also significant at p < 0.05 within the significant voxels corrected for multiple comparisons on the HC congruency effect. The psychopathy group showed no significant congruency effect in this region ( p > 0.3). There was also a significant congruency effect across groups in the right superior parietal lobule (whole-brain analysis); this effect was driven mainly by the psychopathy group ( ).
+
+Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)
+
+Critically, testosterone modulated the congruency effect in the aPFC differently in psychopathic offenders and healthy control subjects (whole-brain analysis on testosterone × group × congruency: MINI coordinates ( x , y , z ): (30, 58, 12); p < 0.001; t = 5.10; for all details, see ). Post hoc analyses revealed that, in the psychopathy group, congruency effects decreased as testosterone levels increased [MNI coordinates ( x , y , z ): (32, 56, 10) and (−30, 58, 8); p = 0.002 and 0.015; t = 4.34 and 3.74]. The modulatory effect of testosterone on congruency was absent in the healthy control subjects ( p ≥ 0.05; ). The whole-brain analysis also showed an effect in the right caudate nucleus and right inferior supramarginal gyrus, driven by reduced congruency effects as a function of testosterone in the psychopathy group ( ; ).
+
+Testosterone modulations of the cerebral congruency effect in psychopathic offenders and healthy control subjects. A , D , Brain image showing testosterone-modulated congruency effects (affect-incongruent−affect-congruent) in the psychopathic offenders in the bilateral aPFC ( A ) and right supramarginal gyrus ( D ). B , E , Bar graphs showing the mean activation (±SEM) of the active voxels within the yellow circles per group. * p < 0.05. ns, Not significant. C , F , Scatterplots showing the correlation of the mean activation of active voxels within the yellow circles with testosterone (log-transformed and standardized) for the healthy control group and the psychopathy group. The ROI activations are presented at p < 0.05, uncorrected for visualization purposes. There are no outliers [Mahalanobis distances D < 4.2 (cutoff at p < 0.05; D = 7.74); ; ]. Healthy control subjects show an increased aPFC activity for the congruency effect and no modulation by testosterone, while in psychopathic offenders endogenous testosterone levels modulate the activity of the aPFC and right supramarginal gyrus.
+
+
+#### Effective connectivity analyses
+
+Given the relevance of aPFC–amygdala connectivity for implementing emotional control as evoked by the AA task ( ), we assessed whether psychopathy also resulted in altered connectivity along that neural pathway. Connectivity analyses using the right aPFC [4-mm-radius sphere; central voxel from main analysis (MNI coordinates: x , 30; y , 58; z , 14)] as the seed region on the congruency effect indicated a significant group difference (PP > HC) with the right amygdala ( ; ROI analysis; extent, 3 voxels; t = 3.82; p = 0.027; MNI coordinates of local maxima: x , 32; y , 0; z , −16). When testing effects for both groups separately, healthy control subjects showed a significant negative coupling between the right aPFC and amygdala (ROI analysis; extent: 3 voxels, t = 3.70; p = 0.036; MNI coordinates of local maxima: x , 32; y , 0; z , −16), while psychopathic offenders showed no differential connectivity effect. Post hoc testing on right amygdala voxels showing the group interaction (threshold, p < 0.05 FWE) indicated a significant positive correlation with testosterone over both groups (ROI analysis; extent, 1 voxel; t = 2.29; p = 0.029; MNI coordinates of local maxima: x , 32; y , 2; z , −16). There was no correlation between aPFC–amygdala connectivity and the PCL-R scores ( p > 0.2).
+
+Group difference on congruency-related aPFC–amygdala connectivity. A , Brain images illustrating the congruency-related modulation of connectivity between the right aPFC (yellow circle, axial slice) and the right amygdala (coronal slice) for the congruency contrast. The activations are presented at p < 0.05, uncorrected for visualization purposes. B , Bar graph visualizing the strength of the congruency-specific change (±SEM) in aPFC–amygdala connectivity for the healthy control subjects and psychopathic offenders. There is a significant negative aPFC–amygdala coupling in the healthy control subjects, which is not present in the psychopathic offenders.
+
+
+
+
+## Discussion
+
+This study indicates that psychopathic offenders show reduced aPFC activity as well as less aPFC–amygdala connectivity during the control of emotional behavior. Emotional control was measured by comparing affect-incongruent and affect-congruent approach–avoidance responses to emotional faces (congruency effect on the AA task; ). When healthy control subjects exerted emotional control, reaction times, aPFC activity, and aPFC–amygdala anticorrelations increased, confirming previous observations ( ). In contrast, psychopathic offenders did not show this typical control-related pattern of aPFC activity and connectivity. In addition, these effects were significantly modulated by endogenous testosterone. Namely, psychopathic individuals with relatively lower testosterone levels showed a neural activity and connectivity pattern that resembled the findings in healthy control subjects, while this pattern was absent in those with higher testosterone levels. This indicates that especially psychopathic individuals with high testosterone levels have less prefrontal regulation of amygdala-driven emotional actions when the control of emotional behavior is required.
+
+### Emotional control in psychopathy
+
+Imaging studies have illustrated an association between psychopathy and altered processing of fear, including altered amygdala responses ( ; ; ), attentional deficits for peripheral stimuli ( ), and moral/empathic insensitivity ( ; ). However, psychopathic offenders also show clear impulsivity problems ( ), for example, when control is required during emotionally provoking situations. To address this relatively unexplored but crucial component of criminal psychopathy, we used a paradigm requiring rule-driven control of emotional actions. With this paradigm, it was possible to move beyond simple motor inhibition and to target the flexible control of emotionally driven action tendencies.
+
+First, the aPFC (also called BA 10) was less active in psychopathic offenders as a function of testosterone. The aPFC is a region crucial for the control of social emotional behavior. When aPFC functioning is temporarily disrupted, participants have increased difficulty in overriding emotional tendencies with rule-driven behavior ( ). Moreover, the aPFC seems especially important for integrating and coordinating multiple cognitive processes to facilitate response selection ( ; ). For example, transcranial magnetic stimulation-induced reduction of aPFC functioning during the control of emotional behavior decreased activity in brain areas associated with rule selection (posterior parietal cortex), while both amygdala activity and automatic action tendencies increased ( ). The current study indicates that psychopathic individuals with especially high testosterone levels recruited the aPFC less when the control of emotional responses was needed. This finding suggests that they have reduced coordination of rule-based behavior with emotional information.
+
+Second, connectivity between the aPFC and amygdala also differed significantly between groups. Healthy control subjects showed a negative aPFC–amygdala coupling during the control of social emotional behavior, whereas psychopathic individuals showed no significant coupling between these regions. Evidence of anatomical connectivity alterations between these regions in psychopathic individuals and the relation of that tract to social emotional behavior modifications support these findings ( ). Although these results cannot resolve the direction of these connectivity effects, a previous study ( ) using this paradigm showed an effective connectivity modulation of emotional control on the connection from aPFC to amygdala. Also, animal studies ( ) suggest strong prefrontal inhibitory connections that control automatic amygdala responses. The absence of this aPFC–amygdala coupling in psychopathic offenders suggests that in this group the aPFC has a reduced ability to inhibit amygdala-driven responses. This study used subtle emotional provocations, but stronger emotional events result in stronger amygdala responses, increasing the bias for automatic emotional behavior ( ). A lack of prefrontal control likely reduces the ability to inhibit these biases and lead to an increased expression of automatic emotional actions even when they are not beneficial ( ; ).
+
+Testosterone administration studies also illustrated a decoupling between the prefrontal cortex and the amygdala, suggesting that testosterone reduces the communication between the PFC and amygdala ( ; ; ) and, within the AA task, reduces top-down control. The association between testosterone levels and enhanced social aggression and dominance seeking, and reduced impulse control in the general population ( ; ; ) supports the relevance of testosterone in this process. Even amygdala responses to angry faces have recently been found to be enhanced after testosterone administration and in psychopathic individuals ( ; ; ). There is a clear association between testosterone and aggression after provocation, which has been related to reduced activity in the orbital frontal cortex, a region just ventral of the aPFC ( ). Interestingly, psychopathic offenders with lower testosterone levels displayed a pattern similar to that in healthy control subjects, while the psychopathic individuals with high testosterone levels showed less aPFC activity and aPFC–amygdala coupling. This could provide a potential vulnerability factor explaining the difference between the goal-directed “successful” psychopath and the “unsuccessful” psychopath with reduced impulse control ( ; ). We hypothesize that especially psychopathic individuals with high testosterone levels fail to inhibit amygdala-driven action tendencies using the aPFC during the control of emotional behavior.
+
+
+Endogenous testosterone levels also modulated control-related activity in the supramarginal gyrus and caudate nucleus of the psychopathy group. The supramarginal gyrus was previously found to be involved during emotional control on the AA task in a healthy student sample ( ). Previous work indicated that it plays an important role in action organization ( ), and that psychopathic individuals show reduced supramarginal gyrus activity compared with control subjects when reasoning about other people’s emotional state ( ). The current findings, emphasizing the role of supramarginal gyrus during emotional control in psychopathic offenders with low testosterone levels, could indicate the facilitation of action preparation in trials with affect-incongruent stimulus–response mapping. The caudate nucleus is important for incorporating predicted action outcome, when selecting the most beneficial behavioral goal ( ), and has previously found to be larger in psychopathy ( ). In light of these findings, our results suggest that psychopathic offenders with low endogenous testosterone levels, as opposed to those with high testosterone levels, have more interference of automatic action tendencies and outcomes associated with the facial emotions (e.g., approach–happy) that are opposite to the required actions during affect-incongruent trials ( ).
+
+
+### Interpretational issues
+
+Individuals with psychopathy have been suggested to have difficulty recognizing emotional expressions. However, this impairment seems quite specific to fear, rather than the emotional expressions used here (anger and happiness; ; ). Furthermore, the groups assessed in this study made comparable numbers of errors, suggesting that psychopathic offenders had no special difficulty in recognizing the clear emotional expressions used in this study.
+
+This study used a relatively subtle manipulation to target the emotional control system. The rationale of this choice was to detect neural vulnerability markers without affecting behavioral performance. Psychopathic offenders performing a more salient behavioral version of the AA task showed reduced avoidance of angry faces ( ). In this study, angry faces evoked numerically similar behavioral effects ( ) and, additionally, aPFC effects ( post hoc inspection of extracted parameters). Although these observations could be interpreted as a sign that psychopathic offenders have a tendency to approach angry faces, those observations were not statistically significant between groups [behavioral and aPFC group effects on angry faces: p > 0.2; p = 0.271; z = 2.54, on the angry–congruency effect in healthy control subjects masked implicitly by group (HC > PP) × angry–congruency interaction]. Future investigation is needed to directly test whether more provocative paradigms induce specific effects for angry faces. A previous study ( ) using this fMRI task in participants with genetic susceptibility for developing aggressive disorders, also found no group-specific behavioral effects. That study suggested that alterations of the aPFC–amygdala pathway might reflect a vulnerability factor for psychopathologies.
+
+Previously, endogenous testosterone modulated the aPFC and aPFC–amygdala coupling in a sample of healthy students ( ). In that study, a different demographic group of healthy control subjects similarly showed a testosterone modulation of aPFC–amygdala coupling, but no testosterone modulation of aPFC activity. This difference in the strength of testosterone-modulatory effects might be related to between-group differences in age (mean healthy control subjects, 41; mean students, 22; ), educational level (staff of forensic psychiatric institute vs university students), or general anxiety [STAI trait, lower in healthy control subjects of the current study; mean (SD): 29 (4.4) and 34 (6.9), respectively; t = −2.605; p = 0.014]. A limitation of this study is the modest sample size. Our focus to exclude moderating factors of comorbid disorders (except antisocial personality disorder) and recent drug use has the advantage that the sample is relatively homogeneous, but future studies using larger samples are needed for replication and to define subsamples.
+
+
+### Conclusion
+
+Psychopathic offenders showed reduced aPFC activity and aPFC–amygdala connectivity during control of emotional actions, suggesting a decreased coordination of emotional information during rule-driven behavior. Moreover, endogenous testosterone modulated the involvement of these neural mechanisms. Psychopathic offenders with high testosterone levels showed less involvement of the aPFC, aPFC–amygdala connectivity, supramarginal gyrus, and caudate nucleus, whereas psychopathic individuals with low testosterone levels recruited the aPFC in a fashion similar to that of healthy control subjects. These findings suggest that a lack of prefrontal control during emotional actions may explain enhanced impulsivity in psychopathic offenders during emotionally provoking situations. They outline a neuroendocrine model underlying impulsive emotional behavior in psychopathy and support the relevance of assessing a potential imbalance in testosterone function to guide treatment. It remains to be seen whether these neuroendocrine alterations of emotional control are also present in highly impulsive or antisocial individuals.
+
+
+
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/26878057.xml b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/26878057.xml
new file mode 100644
index 0000000..10d0863
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/26878057.xml
@@ -0,0 +1,1650 @@
+
+
+
+
+
+ eNeuro
+ eNeuro
+ eneuro
+ eneuro
+ eNeuro
+
+ eNeuro
+
+ 2373-2822
+
+ Society for Neuroscience
+
+
+
+ 26878057
+ 4745181
+ eN-NWR-0107-15
+ 10.1523/ENEURO.0107-15.2016
+
+
+ 1
+
+
+ New Research
+
+ Cognition and Behavior
+
+
+
+
+ Testosterone Modulates Altered Prefrontal Control of Emotional Actions in Psychopathic Offenders123
+ Prefrontal–Amygdala Connectivity in Psychopathy
+
+
+
+ http://orcid.org/0000-0003-3467-1841
+
+ Volman
+ Inge
+
+
+ 1
+
+
+ 2
+
+
+ 3
+
+
+
+
+ von Borries
+ Anna Katinka Louise
+
+
+ 2
+
+
+ 3
+
+
+ 4
+
+
+ 5
+
+
+
+
+ Bulten
+ Berend Hendrik
+
+
+ 5
+
+
+
+
+ Verkes
+ Robbert Jan
+
+
+ 3
+
+
+ 4
+
+
+ 5
+
+
+
+ http://orcid.org/0000-0003-0936-3601
+
+ Toni
+ Ivan
+
+
+ 3
+
+
+
+ http://orcid.org/0000-0002-8863-8978
+
+ Roelofs
+ Karin
+
+
+ 2
+
+
+ 3
+
+
+ Sobell Department of Motor Neuroscience and Movement Disorders, UCL Institute of Neurology, University College London, London WC1N 3BG, United Kingdom
+ Behavioural Science Institute, Radboud University Nijmegen, 6525 HR, Nijmegen, The Netherlands
+ Donders Institute for Brain, Cognition and Behavior, Radboud University Nijmegen, 6525 EN, Nijmegen, The Netherlands
+ Department of Psychiatry, UMC Sint Radboud, 6525 GA, Nijmegen, The Netherlands
+ Pompestichting, 6532 CN, Nijmegen, The Netherlands
+
+
+
+
+
The authors declare no competing financial interests.
+
+
+
+
Author contributions: I.V., A.K.L.v.B., B.H.B., R.J.V., I.T., and K.R. designed research; I.V. and A.K.L.v.B. performed research; I.V., A.K.L.v.B., I.T., and K.R. analyzed data; and I.V., A.K.L.v.B., B.H.B., R.J.V., I.T., and K.R. wrote the paper.
+
+
+
+
This work was supported by VIDI Grant 452-07-008 from the Netherlands Organization for Scientific Research (NWO) awarded to K.R. supporting I.V., Marie Curie Individual Fellowship MSCA-IF-2014-EF_660397 within the European Union's Horizon 2020 Framework Programme awarded to I.V., VICI Grant 453-08-002 from the NWO awarded to I.T., and Starting Grant ERC_StG2012_313749 from the European Research Council and VICI Grant 453-12-001 from the NWO awarded to K.R.
Psychopathic individuals are notorious for their controlled goal-directed aggressive behavior. Yet, during social challenges, they often show uncontrolled emotional behavior. Healthy individuals can control their social emotional behavior through anterior prefrontal cortex (aPFC) downregulation of neural activity in the amygdala, with testosterone modulating aPFC–amygdala coupling. This study tests whether individual differences in this neuroendocrine system relate to the paradoxical lack of emotional control observed in human psychopathic offenders. Emotional control was operationalized with an fMRI-adapted approach–avoidance task requiring rule-driven control over rapid emotional responses. Fifteen psychopathic offenders and 19 matched healthy control subjects made approaching and avoiding movements in response to emotional faces. Control of social emotional behavior was required during affect-incongruent trials, when participants had to override affect-congruent, automatic action tendencies and select the opposite response. Psychopathic offenders showed less control-related aPFC activity and aPFC–amygdala coupling during trials requiring control of emotional actions, when compared with healthy control subjects. This pattern was particularly pronounced in psychopathic individuals with high endogenous testosterone levels. These findings suggest that reduced prefrontal coordination underlies reduced behavioral control in psychopathic offenders during emotionally provoking situations. Even though the modest sample size warrants replication, the modulatory role of endogenous testosterone on the aPFC–amygdala circuit suggests a neurobiological substrate of individual differences that is relevant for the advancement of treatment and the reduction of recidivism.
+
+
+
+ amygdala
+ connectivity
+ emotion
+ fMRI
+ prefrontal
+ psychopathy
+
+
+
+ Netherlands organisation for scientific research
+ 452-07-008
+
+
+ Netherlands organisation for scientific research
+ 453-08-002
+
+
+ European research council
+ ERC_StG2012_313749
+
+
+
+
+
+
+
+
+
+
+
+
+ cover-date
+ January/February 2016
+
+
+
+
+
+
+ Significance Statement
+
Psychopathic criminals are commonly seen as instrumentally abusive and emotionally callous, yet social challenges often trigger uncontrolled emotional behavior in those individuals. This study shows how this paradoxical aspect of psychopathy relates to altered neuroendocrine interactions between testosterone and the cerebral circuit coordinating emotional action tendencies. The anterior prefrontal cortex, a region necessary for controlling emotional behavior, showed blunted responses and reduced connectivity with the amygdala in psychopathic criminals engaged in controlling their emotional action tendencies. This cerebral pattern was strongest in psychopathic individuals with high endogenous testosterone levels. This neuroendocrine signature of altered emotional control highlights the relevance of considering the testosterone level of individual psychopathic patients during treatment of their impulsive behavior.
+
+
+ Introduction
+
Psychopathy is a disorder often associated with blunted emotional responding and increased goal-directed behavior (Blair, 2010; Anderson and Kiehl, 2012). On the other hand, offenders with psychopathy also show a paradoxical increase in impulsive behavior and uncontrolled aggression after emotional provocations (Cornell et al., 1996; Hare, 2003; Patrick et al., 2005; Malterer et al., 2008; Blair, 2010; Anderson and Kiehl, 2012), which may be related to heightened testosterone levels (Stålenheim et al., 1998; Dolan et al., 2001). These two aspects of psychopathy are also distinguished within the most commonly used psychopathy checklist, the Psychopathy Check List-Revised (PCL-R), potentially reflecting differing traits among psychopathic individuals (Hare, 2003; Anderson and Kiehl, 2012). Importantly, enhanced difficulty in controlling emotional impulses, a crucial component of criminal psychopathy associated with PCL-R factor 2, has been largely neglected by cognitive neuroscience. Yet, the clinical relevance of this cognitive trait is large: reduced behavioral control and increased impulsivity predict recidivism in psychopathic offenders (Walters, 2003), and behavioral control in psychopathic offenders appears particularly fragile when dealing with emotionally relevant behavior (Hare, 2003; Blair et al., 2005, chapter 7; Malterer et al. 2008). Accordingly, understanding the neurobiological systems underlying the altered control of social emotional behavior in psychopathic individuals is relevant for improving currently available interventions, which are plagued by low treatment response and high recidivism (Hare, 2003). Here we study those neuroendocrine systems in a group of psychopathic offenders engaged in an experimental paradigm that requires rule-driven control of emotional behavior.
+
Previous investigations of psychopathy showed altered reactivity to emotional material in several brain regions that include the anterior part of the PFC (aPFC) and the amygdala (Anderson and Kiehl, 2012; Blair, 2013; Decety et al., 2015). Furthermore, individuals with psychopathy showed decreased functional and anatomical connectivity between the PFC and amygdala at rest (Craig et al., 2009; Motzkin et al., 2011), an indication that these brain regions might have a reduced ability to interact effectively. Studies in healthy participants have shown that this cerebral circuit is necessary for implementing the control of emotionally relevant actions (Volman et al., 2011a). Namely, aPFC downregulates neural processing in the amygdala during emotional control (Volman et al., 2011a, 2013), while high levels of endogenous testosterone reduce such control-related connectivity between aPFC and amygdala (Volman et al., 2011b). Those findings raise the possibility that aPFC–amygdala connectivity is altered when psychopathic offenders need to control emotionally relevant actions, with high levels of endogenous testosterone exacerbating that altered connectivity.
+
This study tests these hypotheses by measuring brain activity with functional magnetic resonance imaging (fMRI) in 15 psychopathic criminals and 19 matched healthy control subjects dealing with a challenge to control their emotional behavior. The psychopathy sample was obtained by focused and comprehensive screening excluding confounds that are frequently associated with random criminal sampling (e.g., medication use, comorbidity). The social approach–avoidance (AA) task was used to provide reliable indexes of control over social emotional behavior (Fig. 1; Roelofs et al., 2009; Volman et al., 2011a,b). Behaviorally, psychopathic participants previously showed altered AA behavior to explicitly approaching and avoiding emotional faces (Von Borries et al., 2012). Similar findings occurred after testosterone administration in healthy participants (Enter et al., 2014). Interestingly, a more subtle version of the AA task has been shown to be sensitive to testosterone-related alterations and genetic variations in the aPFC–amygdala pathway, while keeping behavior constant across experimental groups (Volman et al., 2011b, 2013), opening the way for isolating neural vulnerability factors (Price and Friston, 1999) in psychopathy. During this task, participants respond to affective faces (happy, angry) presented for a short time with approach and avoidance movements. Automatic emotional tendencies (approach–happy and avoid–angry faces; affect-congruent response conditions) need to be controlled during affect-incongruent response conditions in order to apply the counterintuitive action of approaching angry and avoiding happy faces (Chen and Bargh, 1999; Roelofs et al., 2009). Healthy participants respond more slowly and rely more strongly on the aPFC when emotional control is required, operationalized by the differences evoked between affect-incongruent and affect-congruent trials (Roelofs et al., 2009; Volman et al., 2011b). Accordingly, this study tests whether exerting control over emotionally relevant actions is reflected by reduced functionality of the aPFC–amygdala circuit in psychopathic individuals, suggesting less prefrontal regulation of emotional actions. In addition, it sets out to test whether this alteration is intensified by high levels of endogenous testosterone.
+
+
+
+
The emotional control AA task. The AA task involved the presentation of happy and angry faces, and the performance of approach and avoidance responses. During the AA task, the participants had to select their response according to the perceived emotion of the face. At the beginning of each block of 12 trials, the participants received instructions on whether to pull the joystick toward themselves (approach) or push it away (avoid) when seeing a face with a particular emotion. When viewing happy or angry faces, automatic stimulus–response tendencies trigger corresponding approach or avoidance actions. These tendencies could be followed during the affect-congruent condition (approach–happy, avoid–angry). In contrast, when task instructions required participants to avoid happy faces or to approach angry faces, automatic tendencies needed to be controlled and overridden with the instructed response (affect-incongruent condition). Participants saw the faces and moved the joystick while lying in a MR scanner (top left corner of the table). Figure adapted from Volman et al. (2011a, 2013).
The psychopathic group was recruited from in-patient populations of the Pompestichting and Oldenkotte, forensic psychiatric institutes (TBS-clinics) in the Netherlands. TBS-clinics are facilities for criminal offenders with a mental disorder treated on behalf of the state.
+
Seventeen male psychopathic violent offenders (age range, 23-56 years) participated; all had received a diagnosis with a PCL-R score of ≥26, according to European standards (Hare et al., 1991; Rasmussen et al., 1999; Hildebrand et al., 2004). PCL-R consensus scores were obtained by trained clinicians based on a structured PCL-R interview, clinical status, and history. After the independent scoring, the two raters compared their scores and came to the consensus score. When no consensus could be found, a third independent rater was included in the process. Dutch versions of the National Adult Reading Test and Edinburgh Handedness Inventory were used to assess IQ levels and right-handedness (Oldfield, 1971; Schmand et al., 1991). Twenty-one healthy male control subjects (HCs) matched for age, right-handedness, and IQ, without criminal records or history of psychiatric disorders, were recruited from staff of the clinics. All participants received oral and written information about the experiment and gave written informed consent according to guidelines of the local ethics committee (Commissie Mensengebonden Onderzoek region Arnhem-Nijmegen). Psychiatric exclusion criteria consisted of neurological, axis-I, and axis-II disorders, besides antisocial personality disorder for the psychopathic group. They were screened for these exclusion criteria by trained psychologists using Dutch versions of the Structured Clinical Interview (SCID; Groenestijn et al., 1999) and Mini-International Neuropsychiatric Interview (MINI; Van Vliet et al., 2000) for Diagnostic and Statistical Manual of Mental Disorders, 4th edition, disorders. All participants were asked about drug use and medical/neurological history to exclude the following: alcohol use of >3 units/day, cannabis, or other illicit drug use 1 week before, psychotropic medication other than oxazepam 5 d before, 1 unit of alcohol or oxazepam use within 24 h before the experiment; history of trauma capitis; visual and auditive disorder; and neurological disorder. Furthermore, general exclusion criteria for MRI experiments were applied. Two psychopathic patients (PPs) and two HCs were excluded from the analyses, due to incomplete scanning procedures (1 PP, 1 HC) or too many errors on the task (>16%, representing the outlier with a z-score >3). The final groups did not differ in age, IQ, and handedness (see Table 1).
+
+
+
+
Demographical data
+
+
+
+
+
+
+
+
+
+
+
+
Psychopathic offenders(n = 15)
+
HCs(n = 19)
+
p value
+
+
+
+
+
Age
+
37.8 (7.9)
+
40.7 (10.3)
+
0.368
+
+
+
IQ
+
101 (10)
+
102 (9)
+
0.761
+
+
+
Handedness
+
50.7 (81)
+
59.2 (62)
+
0.729
+
+
+
PCL-R total
+
30.4 (3.5)
+
+
+
+
+
PCL-R F1
+
12.1 (2.6)
+
+
+
+
+
PCL-R F2
+
14.1 (2.3)
+
+
+
+
+
+
+
+
Values are presented as the mean (SD), unless otherwise indicated. F1, factor 1; F2, factor 2.
+
+
+
+
+
+ Procedure
+
Two test sessions took place. During the first session, right-handedness, IQ, MINI, and SCID were assessed. During the second session, participants completed several questionnaires upon arrival in the laboratory, including the State-Trait Anxiety Inventory (STAI) to measure anxiety levels (Spielberger, 1983). Next, they provided saliva for the testosterone measurement. Afterward, participants were positioned in the 1.5 T MR scanner and familiarized with the task setup. Immediately after this, the fMRI session started with the AA task (duration, 30 min) followed by another task (not included in this report). After a short break outside the scanner, the anatomical scan (duration, 5 min) and an unrelated task were acquired in the side-by-side 3 T MR scanner.
+
+
+ Experimental task
+
The AA task consisted of 24 blocks (with 12 trials per block and a baseline period of 21-24 s) during which participants had to respond to visually presented faces either by pulling a joystick toward themselves (approach) or by pushing it away from themselves (avoid; Fig. 1). The participants had to categorize faces as happy, angry, and neutral (filler items), based on their affective expressions. During each block, two of the three affective expressions were presented as stimuli, because only two responses could be given to categorize the stimulus. This resulted in six different block types each used four times, representing the affect (happy–angry, happy–neutral, angry–neutral) × movement (approach–avoid) combinations. At the start of each block, participants received written instructions regarding the required response mapping. The affect × movement combinations were pseudorandomly and evenly distributed (with no affect combination repetition), and the combination of the first block was counterbalanced across participants. Within each block, affective expressions and gender types were pseudorandomly presented, avoiding three or more sequential presentations of the same expression/gender, and two presentations of the same facial model. Each face was presented for 100 ms, preceded by a 300 ms blank screen, and followed by the participant’s response, a blank screen, and by a pseudorandom intertrial interval (ITI; 1-3 s). A baseline period of 21-24 s preceded each block. The faces were from 36 models (18 male) obtained from several databases (Ekman and Friesen, 1976; Matsumoto and Ekman, 1988; Lundqvist et al., 1998; Martinez and Benavente, 1998), each showing all expressions. The pictures were in grayscale, matched for brightness and contrast values, and displayed against a black background. To exclude influence from hair and nonfacial contours, the faces were trimmed. Joystick displacements of >80% along the sagittal plane within 2 s from stimulus presentation were marked as valid responses. Invalid responses were signaled for 1 s with written feedback stating “you did not move your joystick far enough.” After moving the joystick, participants had to return to the starting position (defined as the central area extending 20% along the sagittal plane) before the end of the ITI. Otherwise, visual feedback indicated “return the joystick to the starting position,” and the ITI was repeated after participants returned the joystick. The training at the beginning consisted of six blocks; one block of eight trials for each of the six affect × movement combinations. Different visual stimuli were used during the training and scanning blocks.
+
+
+ Materials and apparatus
+
The fMR images were acquired on a 1.5 T MRI scanner (Avanto, Siemens Medical Systems) with an eight-channel head coil using a multiecho generalized autocalibrating partially parallel acquisitions (GRAPPA) sequence [Poser et al., 2006; repetition time (TR), 2.14 ms; five echo times (TEs), 9.4/21/33/44/56 ms; 34 transversal slices; ascending acquisition; distance factor, 17%; effective voxel size, 3.3 × 3.3 × 3.5 mm; field of view (FOV), 212 mm]. High-resolution anatomical images were acquired on a 3 T MRI scanner with a 32-channel head coil using a magnetization prepared rapid gradient echo sequence (TR, 2300 ms; TE, 3.03 ms; 192 sagittal slices; voxel size, 1.0 × 1.0 × 1.0 mm; FOV, 256 mm).
+
An MR-compatible joystick (Fiber Optic Joystick, Current Designs; sampling rate, 550 Hz) was placed on participants’ abdomens to ensure comfortable push-and-pull movements (Fig. 1). Participants wore MR-compatible headphones to reduce scanner noise (Commander XG MRI Audio System, Resonance Technologies). Stimuli were projected at the center of a screen, viewed via a mirror above the participant’s head, with a visual angle of 4° × 6° (width × height). Stimuli presentation and acquisition of joystick positions were controlled by a PC running Presentation version 13 (http://www.neurobs.com).
+
+
+ Salivary measurements
+
Participants filled two Salicaps (IBL) with saliva for testosterone measurement, which were stored at −25°C. Testosterone concentration was measured using competitive chemiluminescence immunoassay with a sensitivity of 0.0025 ng/ml (IBL International, Tecan). Intra-assay and interassay coefficients are between 10% and 12%. To control variables influencing testosterone levels, participants were instructed to refrain from any food, cigarettes, and drinks (except water) for 1 h before the experiment.
+
+
+ Behavioral analysis
+
Behavioral data was analyzed using MATLAB version 7.9 (MathWorks) and PASW Statistics 18 (SPSS Inc.). First, to obtain a precise measure of movement onset [reaction time (RT)], the joystick movement for each trial was reconstructed using the joystick displacement measurements. Excluded trials showed a joystick movement in the wrong direction, an extreme RT (<150 or >1500 ms), peak velocity (<0.1 cm/s), or movement time (>400 ms); or an error rate of above chance level in a block (in that case, the whole block was excluded). RTs and testosterone levels were log transformed to obtain a normal distribution. Second, following previous studies (Roelofs et al., 2009; Volman et al., 2011b), we conducted three-way repeated-measures ANOVA (ANCOVArm) on the mean RT and error rates, with factors group (PP, HC), movement (approach, avoid), and valence (happy, angry), including standardized testosterone and STAI state as covariate. A measure of anxiety (STAI) was included to account for the effects of psychopathy type (e.g., primary vs secondary); and the possible effects on emotional behavior, hormonal levels, amygdala, and prefrontal cortex functioning (Freitas-Ferrari et al., 2010; Koenigs et al., 2011; Giltay et al., 2012; Fouche et al., 2013). The α-level was set at p < 0.05.
Imaging data were preprocessed and analyzed using SPM8 (Statistical Parametric Mapping; http://www.fil.ion.ucl.ac.uk/spm). The first four volumes of each participant’s dataset were discarded to allow for T1 equilibration. Given the multiecho GRAPPA MR sequence (Poser et al., 2006), head motion parameters were estimated on MR images with the shortest TE (9.4 ms), since these are least affected by possible artifacts. These motion correction parameters, estimated using a least-squares approach with six rigid body transformation parameters (translations, rotations), were applied to the five echo images collected for each excitation. After spatial realignment, the five echo images were combined into a single MR volume using an optimized echo weighting method (Poser et al., 2006). The time series for each voxel was temporally realigned to the first slice in time. The T1-weighted image was spatially coregistered to the mean of the functional images. The fMRI time series were transformed and resampled at an isotropic voxel size of 2 mm into standard Montreal Neurological Institute (MNI) space by unified segmentation and normalization using the coregistered T1-weighted image (Ashburner and Friston, 2005). The normalized functional images were spatially smoothed using an isotropic 8 mm full-width at half-maximum Gaussian kernel.
+
The fMRI time series of each subject were further analyzed using an event-related approach in the context of general linear model, including the following effects: approach–happy, approach–neutral, approach–angry, avoid–happy, avoid–neutral, and avoid–angry. Trials excluded from behavioral analyses and periods of instructions or feedback were modeled as regressors. Vectors describing the time of picture presentation (onset) and RT of each event (duration) were convolved with the canonical hemodynamic response function. Potential confounding effects of residual head movement were modeled using original, squared, cubic, first-order, and second-order derivatives of the movement correction parameters (Lund et al., 2005). Three further regressors, describing the time course of signal intensities of white matter, CSF, and the portion of the MR image outside the skull were also added. This procedure accounts for image intensity shifts due to hand movements within or near the magnetic field of the scanner (Verhagen et al., 2006). Finally, fMRI time series were high-pass filtered (cutoff 120 s). Temporal autocorrelation was modeled as a first-order autoregressive process.
+
+
+ Group analyses
+
Consistent effects across participants and between groups were tested using a random-effects multiple regression analysis that included six contrast images (approach–happy, approach–neutral, approach–angry, avoid–happy, avoid–neutral, avoid–angry) per participant. Together, these images represented the estimated cerebral effects from 12 conditions of the experimental design [group (PP, HC) × valence (happy, neutral, angry) × response (approach, avoid)]. Standardized log-transformed testosterone and standardized STAI state levels were included in the multiple regression analysis as condition-specific [group (PP, HC) × valence (happy, neutral, angry) × response (approach, avoid)] regressors, generating another 12 regressors per variable.
+
All analyses assessed the congruency effect, reflecting task-related differences of affect-incongruent (approach–angry, avoid–happy) versus affect-congruent trials (approach–happy, avoid–angry; Roelofs et al., 2009; Volman et al., 2011b). We considered two effects. First, to test for general effects of congruency, we performed an analysis on the congruency effect over both groups and for each group separately. When assessing the effects of one group explicitly, we also tested whether those effects were specific to that group and were significantly weaker in the other group (at p < 0.05 uncorrected) by masking the statistical map describing the congruency effect in the first group (using multiple comparisons correction, see below) with the statistical map describing the group × congruency contrast. Second, to test whether testosterone differentially modulated the control of emotionally relevant actions in the groups, we performed a group × congruency contrast on the regressor parametrizing interindividual differences in testosterone on task-related conditions. If such an interaction is present, the testosterone modulation on the congruency effect of each group separately is considered. In addition to whole-brain analyses, we used a volume of interest (VOI) on coordinates previously found to be modulated by testosterone during the congruency effect in healthy students (two 8-mm-radius spheres centered on the following MNI coordinates: x, −30; y, 58; and z, 2; and x, 32; y, 54; and z, 8; Volman et al., 2011b).
+
The reported activations are corrected for multiple comparisons using familywise error (FWE) correction. For whole-brain analyses, we made inferences at cluster level (FWE: p < 0.05, corresponding to a cluster size of >140 on the basis of intensity threshold, p < 0.001). For VOI analyses, we made inferences at voxel-level (FWE corrected, p < 0.05; Worsley et al., 1996; Friston, 1997). Anatomical inference is drawn by superimposing SPM showing significant signal changes on structural images of participants. For anatomical accuracy, we report only activation peaks in gray matter. Anatomical landmarks were identified using the atlas of Duvernoy et al. (1991). Brodmann areas (BAs) were assigned by superimposing significant SPM on the SPM anatomy toolbox (Eickhoff et al., 2005) and MRIcron template (http://www.mccauslandcenter.sc.edu/mricro/mricron/).
+
+
+
+ Connectivity analyses
+
The aim of the following analysis was to test whether inter-regional coupling of the aPFC (see Results) with the amygdala and other brain regions during the congruency effect was different between the groups and modulated by testosterone. To test for these effects, we used the psychophysiological interactions (PPIs) method (Friston et al., 1997). More specifically, we tested for significant differences between the regression coefficients of each voxel over the right aPFC during the affect-incongruent versus the affect-congruent conditions. To select voxels to be included in the VOI, we used the following anatomical constraints (Stephan et al., 2010): for each participant, selected voxels fell within a sphere with a radius of 4 mm around the peak voxel corresponding to the activated cluster of the congruency effect over both groups (coordinates: x, 30; y, 58; z, 14; see Results). Participant specific contrast images were generated describing the PPI between the time courses of the right aPFC VOI and affect-incongruent versus affect-congruent conditions. Group differences and testosterone modulations on task-related coupling between the aPFC and other regions were then assessed using a multiple regression design on participant-specific contrast images with their corresponding testosterone (log-transformed, standardized) and STAI state (standardized) levels as subject- and group-specific regressors. In addition to whole-brain analyses, we assessed significant voxel-level effects (FWE corrected for multiple comparisons, p < 0.05) within the amygdala, defined on the Automated Anatomical Labeling atlas (Tzourio-Mazoyer et al., 2002) using the WFU PickAtlas tool (Maldjian et al., 2003).
+
+
+
+ Results
+
+ Behavioral results
+
Fifteen psychopathic criminals (PPs; PCL-R score of ≥26, according to European standards (Rasmussen et al., 1999; Hare, 2003; Hildebrand et al., 2004) and 19 HCs (for demographics, see Table 1) were included in the analyses. Participants performed the task accurately and consistently (error rates: PPs, 7.9%; HCs, 7.3%; omissions: PPs, 1.6%; HCs, 1.5%; undefined responses: PPs, 0.9%; HCs, 0.3%; Table 2).
+
+
+
+
RTs and error rates for each group and factor of the AA task
+
+
+
+
+
+
+
+
+
+
+
+
+
Psychopathic offenders
+
HCs
+
+
+
+
Approach
+
Avoid
+
Approach
+
Avoid
+
+
+
+
+
Errors (%)
+
+
+
+
+
+
+
Happy
+
3.2 (0.9)
+
8.9 (1.8)
+
2.4 (0.8)
+
7.7 (1.1)
+
+
+
Neutral
+
6.1 (1.3)
+
5.8 (1.1)
+
7.1 (1.4)
+
5.2 (1.0)
+
+
+
Angry
+
10.1 (2.2)
+
13.1 (2.1)
+
9.6 (1.8)
+
11.6 (1.8)
+
+
+
RT (ms)
+
+
+
+
+
+
+
Happy
+
554 (25)
+
625 (35)
+
553 (23)
+
603 (25)
+
+
+
Neutral
+
666 (28)
+
687 (31)
+
639 (21)
+
668 (24)
+
+
+
Angry
+
630 (25)
+
665 (33)
+
620 (24)
+
630 (23)
+
+
+
+
+
+
Values are presented as the mean (SE).
+
+
+
+
A significant movement × valence interaction for the RTs indicated that, over groups, participants responded more slowly during affect-incongruent (approach–angry, avoid–happy) than during affect-congruent trials (approach–happy, avoid–angry; F(1,29) = 10.4, p = 0.003; Fig. 2). This congruency effect replicates the behavioral results from previous fMRI studies (Roelofs et al., 2009; Volman et al., 2011b, 2013). Furthermore, there were main effects of movement (F(1,29) = 26.3, p < 0.001) and valence (F(1,29) = 28.7, p < 0.001), reflecting the slowing of avoidance movements and responses to angry faces in general (Table 2). There were no significant effects involving group, including no main effect (p > 0.3). The congruency effect correlated positively (without corrections for multiple comparisons) with the PCL-R total score (p = 0.048, R = 0.517, respectively). Excluding anxiety from the analyses did not affect the outcomes. Moreover, when including the neutral conditions in the analyses, the movement × valence (happy, neutral, angry) interaction for RTs remained significant (F(1,28) = 5.5, p = 0.010), showing that neutral approach–avoidance effects are intermediary compared with happy and angry (Table 2).
+
+
+
+
Behavioral results. Mean RTs (±SEM) for the affect-congruent and affect-incongruent conditions of the AA task for the healthy control subjects and psychopathic offenders. The groups were significantly slower to provide affect-incongruent responses (approach–angry; avoid–happy) than affect-congruent responses (approach–happy; avoid–angry), with no significant group differences.
+
+
+
+
For the error rates, the three-way ANCOVArm showed main effects of movement (F(1,29) = 27.5, p < 0.001), valence (F(1,29) = 25.9, p < 0.001), and testosterone (F(1,29) = 4.6, p = 0.040), and a valence × testosterone interaction (F(1,29) = 4.3, p = 0.047). There were no other significant effects for the error rates (p > 0.15).
+
Endogenous testosterone levels [median (SD): PPs, 101 pg/ml (70 pg/ml); HCs, 90 pg/ml (46 pg/ml)] and state anxiety levels [STAI mean (SD): PPs, 32 (8); HCs, 32 (5)] did not differ between groups (p > 0.4), and showed no correlations with psychopathy (PCL-R) scores or with each other (p > 0.1).
To assess the two main questions of this study, we isolated cerebral structures showing stronger responses during affect-incongruent than affect-congruent trials (congruency effect), and cerebral structures in which the congruency effect was modulated by testosterone levels.
+
The results showed a significant congruency effect across groups in the aPFC [ROI analysis: MNI coordinates (x, y, z): (30, 58, 14) and (−30 58 10); pFWE = 0.001 and 0.036; t = 4.46 and 3.43; for further details, see Table 3]. As expected, this effect was driven by the healthy control group, and it was significantly weaker in the psychopathic offenders [pFWE = 0.001 and 0.040; t = 4.58 and 3.40, on the congruency effect in healthy control subjects masked implicitly by group (HC > PP) × congruency interaction]. The implicit masking demonstrates that the group × congruency interaction is also significant at puncorrected < 0.05 within the significant voxels corrected for multiple comparisons on the HC congruency effect. The psychopathy group showed no significant congruency effect in this region (pFWE > 0.3). There was also a significant congruency effect across groups in the right superior parietal lobule (whole-brain analysis); this effect was driven mainly by the psychopathy group (Table 3).
+
+
+
+
Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Anatomical region
+
Putative BA
+
Side
+
x
+
+
y
+
+
z
+
+
Voxels (n)
+
p value
+
t value
+
+
+
+
+
+ Whole-brain effects
+
+
+
+
+
+
+
+
+
+
+
+
+ Congruency effect over groups
+
+
+
+
+
+
+
+
+
+
+
+
Cuneus
+
18
+
R/L
+
−16
+
−96
+
24
+
634
+
<0.001
+
4.85
+
+
+
SPL/Superior occipital gyrus
+
7/19
+
L
+
−30
+
−76
+
34
+
254
+
0.004
+
4.37
+
+
+
IFG
+
45/47
+
L
+
−52
+
22
+
−10
+
223
+
0.007
+
4.37
+
+
+
Cuneus
+
18
+
R
+
18
+
−98
+
16
+
190
+
0.016
+
4.52
+
+
+
+ Congruency effect for psychopathy group
+
+
+
+
+
+
+
+
+
+
+
+
SPL/superior occipital gyrus
+
7/19
+
R/L
+
8
+
−82
+
38
+
925
+
<0.001
+
4.98
+
+
+
Angular gyrus
+
39/19
+
L
+
−30
+
−72
+
34
+
337
+
0.001
+
4.35
+
+
+
Superior temporal gyrus
+
42
+
L
+
−32
+
−32
+
6
+
214
+
0.009
+
4.51
+
+
+
SPL
+
7
+
R
+
28
+
−74
+
46
+
158
+
0.034
+
4.58
+
+
+
Cerebellum
+
+
L
+
−24
+
−66
+
−38
+
146
+
0.046
+
4.63
+
+
+
+ Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
+
+
+
+
+
+
+
+
+
+
+
+
aPFC
+
10
+
R
+
30
+
58
+
12
+
391
+
<0.001
+
5.10
+
+
+
Supramarginal gyrus
+
40
+
R
+
54
+
−42
+
54
+
325
+
0.001
+
4.66
+
+
+
Caudate nucleus
+
+
R
+
10
+
10
+
2
+
273
+
0.002
+
4.69
+
+
+
Putamen/Insula
+
+
L
+
−34
+
6
+
−10
+
176
+
0.022
+
5.22
+
+
+
Cerebellum
+
+
R
+
18
+
−76
+
−38
+
188
+
0.016
+
5.09
+
+
+
+ Negative testosterone modulation of Congruency effect in psychopathy
+
+
+
+
+
+
+
+
+
+
+
+
Supramarginal gyrus
+
40
+
R
+
52
+
−40
+
54
+
657
+
<0.001
+
5.32
+
+
+
Precentral/superior frontal gyrus
+
6
+
R/L
+
6
+
22
+
66
+
471
+
<0.001
+
5.65
+
+
+
Caudate nucleus
+
+
R
+
6
+
6
+
4
+
228
+
0.006
+
4.28
+
+
+
+ VOI on bilateral aPFC
+
+
+
+
+
+
+
+
+
+
+
+
+ Congruency effect over groups
+
+
10
+
R
+
30
+
58
+
14
+
31
+
0.001
+
4.46
+
+
+
+
10
+
L
+
−30
+
58
+
10
+
5
+
0.036
+
3.43
+
+
+
+ Congruency effect in healthy control subjects
+
+
10
+
R
+
32
+
58
+
14
+
12
+
0.001
+
4.58
+
+
+
+
10
+
L
+
−34
+
52
+
4
+
2
+
0.040
+
3.40
+
+
+
+ Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
+
+
10
+
R
+
30
+
58
+
12
+
145
+
<0.001
+
5.10
+
+
+
10
+
L
+
−24
+
56
+
6
+
15
+
0.010
+
3.87
+
+
+
+ Negative testosterone modulation of congruency effect in psychopathy
+
+
10
+
R
+
32
+
56
+
10
+
77
+
0.002
+
4.34
+
+
+
10
+
L
+
−30
+
58
+
8
+
17
+
0.015
+
3.74
+
+
+
+
+
+
Coordinates are defined in MNI (x, y, z) space. The p values represent the FWE cluster-level corrected values for the whole-brain analyses and FWE voxel-level corrected values for the VOI analyses. IFG, Inferior frontal gyrus; L, left; R, right; SPL, superior parietal lobule.
+
+
+
+
Critically, testosterone modulated the congruency effect in the aPFC differently in psychopathic offenders and healthy control subjects (whole-brain analysis on testosterone × group × congruency: MINI coordinates (x, y, z): (30, 58, 12); pFWE < 0.001; t = 5.10; for all details, see Table 3). Post hoc analyses revealed that, in the psychopathy group, congruency effects decreased as testosterone levels increased [MNI coordinates (x, y, z): (32, 56, 10) and (−30, 58, 8); pFWE = 0.002 and 0.015; t = 4.34 and 3.74]. The modulatory effect of testosterone on congruency was absent in the healthy control subjects (pFWE ≥ 0.05; Fig. 3A–C). The whole-brain analysis also showed an effect in the right caudate nucleus and right inferior supramarginal gyrus, driven by reduced congruency effects as a function of testosterone in the psychopathy group (Fig. 3D–F; Table 3).
+
+
+
+
Testosterone modulations of the cerebral congruency effect in psychopathic offenders and healthy control subjects. A, D, Brain image showing testosterone-modulated congruency effects (affect-incongruent−affect-congruent) in the psychopathic offenders in the bilateral aPFC (A) and right supramarginal gyrus (D). B, E, Bar graphs showing the mean activation (±SEM) of the active voxels within the yellow circles per group. *pFWE < 0.05. ns, Not significant. C, F, Scatterplots showing the correlation of the mean activation of active voxels within the yellow circles with testosterone (log-transformed and standardized) for the healthy control group and the psychopathy group. The ROI activations are presented at p < 0.05, uncorrected for visualization purposes. There are no outliers [Mahalanobis distances D
+2i < 4.2 (cutoff at p < 0.05; D = 7.74); Barnett and Lewis, 1978; Stevens, 1996]. Healthy control subjects show an increased aPFC activity for the congruency effect and no modulation by testosterone, while in psychopathic offenders endogenous testosterone levels modulate the activity of the aPFC and right supramarginal gyrus.
+
+
+
+
+
+ Effective connectivity analyses
+
Given the relevance of aPFC–amygdala connectivity for implementing emotional control as evoked by the AA task (Volman et al., 2011a, 2013), we assessed whether psychopathy also resulted in altered connectivity along that neural pathway. Connectivity analyses using the right aPFC [4-mm-radius sphere; central voxel from main analysis (MNI coordinates: x, 30; y, 58; z, 14)] as the seed region on the congruency effect indicated a significant group difference (PP > HC) with the right amygdala (Fig. 4A,B; ROI analysis; extent, 3 voxels; t = 3.82; pFWE = 0.027; MNI coordinates of local maxima: x, 32; y, 0; z, −16). When testing effects for both groups separately, healthy control subjects showed a significant negative coupling between the right aPFC and amygdala (ROI analysis; extent: 3 voxels, t = 3.70; pFWE = 0.036; MNI coordinates of local maxima: x, 32; y, 0; z, −16), while psychopathic offenders showed no differential connectivity effect. Post hoc testing on right amygdala voxels showing the group interaction (threshold, p < 0.05 FWE) indicated a significant positive correlation with testosterone over both groups (ROI analysis; extent, 1 voxel; t = 2.29; pFWE = 0.029; MNI coordinates of local maxima: x, 32; y, 2; z, −16). There was no correlation between aPFC–amygdala connectivity and the PCL-R scores (p > 0.2).
+
+
+
+
Group difference on congruency-related aPFC–amygdala connectivity. A, Brain images illustrating the congruency-related modulation of connectivity between the right aPFC (yellow circle, axial slice) and the right amygdala (coronal slice) for the congruency contrast. The activations are presented at p < 0.05, uncorrected for visualization purposes. B, Bar graph visualizing the strength of the congruency-specific change (±SEM) in aPFC–amygdala connectivity for the healthy control subjects and psychopathic offenders. There is a significant negative aPFC–amygdala coupling in the healthy control subjects, which is not present in the psychopathic offenders.
+
+
+
+
+
+
+
+ Discussion
+
This study indicates that psychopathic offenders show reduced aPFC activity as well as less aPFC–amygdala connectivity during the control of emotional behavior. Emotional control was measured by comparing affect-incongruent and affect-congruent approach–avoidance responses to emotional faces (congruency effect on the AA task; Roelofs et al. 2009). When healthy control subjects exerted emotional control, reaction times, aPFC activity, and aPFC–amygdala anticorrelations increased, confirming previous observations (Volman et al., 2011b, 2013). In contrast, psychopathic offenders did not show this typical control-related pattern of aPFC activity and connectivity. In addition, these effects were significantly modulated by endogenous testosterone. Namely, psychopathic individuals with relatively lower testosterone levels showed a neural activity and connectivity pattern that resembled the findings in healthy control subjects, while this pattern was absent in those with higher testosterone levels. This indicates that especially psychopathic individuals with high testosterone levels have less prefrontal regulation of amygdala-driven emotional actions when the control of emotional behavior is required.
+
+ Emotional control in psychopathy
+
Imaging studies have illustrated an association between psychopathy and altered processing of fear, including altered amygdala responses (Blair, 2010; Moul et al., 2012; Decety et al., 2015), attentional deficits for peripheral stimuli (Baskin-Sommers et al., 2011), and moral/empathic insensitivity (Decety et al., 2013; Marsh and Cardinale, 2012). However, psychopathic offenders also show clear impulsivity problems (Hare, 2003), for example, when control is required during emotionally provoking situations. To address this relatively unexplored but crucial component of criminal psychopathy, we used a paradigm requiring rule-driven control of emotional actions. With this paradigm, it was possible to move beyond simple motor inhibition and to target the flexible control of emotionally driven action tendencies.
+
First, the aPFC (also called BA 10) was less active in psychopathic offenders as a function of testosterone. The aPFC is a region crucial for the control of social emotional behavior. When aPFC functioning is temporarily disrupted, participants have increased difficulty in overriding emotional tendencies with rule-driven behavior (Volman et al., 2011a). Moreover, the aPFC seems especially important for integrating and coordinating multiple cognitive processes to facilitate response selection (Ramnani and Owen, 2004; Haggard, 2008). For example, transcranial magnetic stimulation-induced reduction of aPFC functioning during the control of emotional behavior decreased activity in brain areas associated with rule selection (posterior parietal cortex), while both amygdala activity and automatic action tendencies increased (Volman et al., 2011a). The current study indicates that psychopathic individuals with especially high testosterone levels recruited the aPFC less when the control of emotional responses was needed. This finding suggests that they have reduced coordination of rule-based behavior with emotional information.
+
Second, connectivity between the aPFC and amygdala also differed significantly between groups. Healthy control subjects showed a negative aPFC–amygdala coupling during the control of social emotional behavior, whereas psychopathic individuals showed no significant coupling between these regions. Evidence of anatomical connectivity alterations between these regions in psychopathic individuals and the relation of that tract to social emotional behavior modifications support these findings (Von Der Heide et al., 2013). Although these results cannot resolve the direction of these connectivity effects, a previous study (Volman et al., 2013) using this paradigm showed an effective connectivity modulation of emotional control on the connection from aPFC to amygdala. Also, animal studies (Quirk and Gehlert, 2003) suggest strong prefrontal inhibitory connections that control automatic amygdala responses. The absence of this aPFC–amygdala coupling in psychopathic offenders suggests that in this group the aPFC has a reduced ability to inhibit amygdala-driven responses. This study used subtle emotional provocations, but stronger emotional events result in stronger amygdala responses, increasing the bias for automatic emotional behavior (Quirk and Gehlert, 2003). A lack of prefrontal control likely reduces the ability to inhibit these biases and lead to an increased expression of automatic emotional actions even when they are not beneficial (Quirk and Gehlert, 2003; Volman et al., 2011a).
+
Testosterone administration studies also illustrated a decoupling between the prefrontal cortex and the amygdala, suggesting that testosterone reduces the communication between the PFC and amygdala (Eisenegger et al., 2011; Van Wingen et al., 2011; Bos et al., 2012) and, within the AA task, reduces top-down control. The association between testosterone levels and enhanced social aggression and dominance seeking, and reduced impulse control in the general population (Van Wingen et al., 2011; Montoya et al., 2012; Carré et al., 2013) supports the relevance of testosterone in this process. Even amygdala responses to angry faces have recently been found to be enhanced after testosterone administration and in psychopathic individuals (Van Wingen et al., 2011; Carré et al., 2013; Radke et al., 2015). There is a clear association between testosterone and aggression after provocation, which has been related to reduced activity in the orbital frontal cortex, a region just ventral of the aPFC (Mehta and Beer, 2010). Interestingly, psychopathic offenders with lower testosterone levels displayed a pattern similar to that in healthy control subjects, while the psychopathic individuals with high testosterone levels showed less aPFC activity and aPFC–amygdala coupling. This could provide a potential vulnerability factor explaining the difference between the goal-directed “successful” psychopath and the “unsuccessful” psychopath with reduced impulse control (Gao and Raine, 2010; Anderson and Kiehl, 2012). We hypothesize that especially psychopathic individuals with high testosterone levels fail to inhibit amygdala-driven action tendencies using the aPFC during the control of emotional behavior.
+
+
Endogenous testosterone levels also modulated control-related activity in the supramarginal gyrus and caudate nucleus of the psychopathy group. The supramarginal gyrus was previously found to be involved during emotional control on the AA task in a healthy student sample (Volman et al., 2011b). Previous work indicated that it plays an important role in action organization (Jubault et al., 2007), and that psychopathic individuals show reduced supramarginal gyrus activity compared with control subjects when reasoning about other people’s emotional state (Sommer et al., 2010). The current findings, emphasizing the role of supramarginal gyrus during emotional control in psychopathic offenders with low testosterone levels, could indicate the facilitation of action preparation in trials with affect-incongruent stimulus–response mapping. The caudate nucleus is important for incorporating predicted action outcome, when selecting the most beneficial behavioral goal (Grahn et al., 2008), and has previously found to be larger in psychopathy (Glenn et al., 2010). In light of these findings, our results suggest that psychopathic offenders with low endogenous testosterone levels, as opposed to those with high testosterone levels, have more interference of automatic action tendencies and outcomes associated with the facial emotions (e.g., approach–happy) that are opposite to the required actions during affect-incongruent trials (Grahn et al., 2008).
+
+
+ Interpretational issues
+
Individuals with psychopathy have been suggested to have difficulty recognizing emotional expressions. However, this impairment seems quite specific to fear, rather than the emotional expressions used here (anger and happiness; Marsh and Blair, 2008; Von Borries et al., 2012). Furthermore, the groups assessed in this study made comparable numbers of errors, suggesting that psychopathic offenders had no special difficulty in recognizing the clear emotional expressions used in this study.
+
This study used a relatively subtle manipulation to target the emotional control system. The rationale of this choice was to detect neural vulnerability markers without affecting behavioral performance. Psychopathic offenders performing a more salient behavioral version of the AA task showed reduced avoidance of angry faces (Von Borries et al., 2012). In this study, angry faces evoked numerically similar behavioral effects (Table 2) and, additionally, aPFC effects (post hoc inspection of extracted parameters). Although these observations could be interpreted as a sign that psychopathic offenders have a tendency to approach angry faces, those observations were not statistically significant between groups [behavioral and aPFC group effects on angry faces: p > 0.2; pFWE = 0.271; z = 2.54, on the angry–congruency effect in healthy control subjects masked implicitly by group (HC > PP) × angry–congruency interaction]. Future investigation is needed to directly test whether more provocative paradigms induce specific effects for angry faces. A previous study (Volman et al., 2013) using this fMRI task in participants with genetic susceptibility for developing aggressive disorders, also found no group-specific behavioral effects. That study suggested that alterations of the aPFC–amygdala pathway might reflect a vulnerability factor for psychopathologies.
+
Previously, endogenous testosterone modulated the aPFC and aPFC–amygdala coupling in a sample of healthy students (Volman et al., 2011b). In that study, a different demographic group of healthy control subjects similarly showed a testosterone modulation of aPFC–amygdala coupling, but no testosterone modulation of aPFC activity. This difference in the strength of testosterone-modulatory effects might be related to between-group differences in age (mean healthy control subjects, 41; mean students, 22; Peper et al., 2011), educational level (staff of forensic psychiatric institute vs university students), or general anxiety [STAI trait, lower in healthy control subjects of the current study; mean (SD): 29 (4.4) and 34 (6.9), respectively; t(37) = −2.605; p = 0.014]. A limitation of this study is the modest sample size. Our focus to exclude moderating factors of comorbid disorders (except antisocial personality disorder) and recent drug use has the advantage that the sample is relatively homogeneous, but future studies using larger samples are needed for replication and to define subsamples.
+
+
+ Conclusion
+
Psychopathic offenders showed reduced aPFC activity and aPFC–amygdala connectivity during control of emotional actions, suggesting a decreased coordination of emotional information during rule-driven behavior. Moreover, endogenous testosterone modulated the involvement of these neural mechanisms. Psychopathic offenders with high testosterone levels showed less involvement of the aPFC, aPFC–amygdala connectivity, supramarginal gyrus, and caudate nucleus, whereas psychopathic individuals with low testosterone levels recruited the aPFC in a fashion similar to that of healthy control subjects. These findings suggest that a lack of prefrontal control during emotional actions may explain enhanced impulsivity in psychopathic offenders during emotionally provoking situations. They outline a neuroendocrine model underlying impulsive emotional behavior in psychopathy and support the relevance of assessing a potential imbalance in testosterone function to guide treatment. It remains to be seen whether these neuroendocrine alterations of emotional control are also present in highly impulsive or antisocial individuals.
+
+
+
+
+
+ References
+
+
+Anderson
+NE, Kiehl
+KA (2012) The psychopath magnetized: insights from brain imaging.
+
+16:52-60. 10.1016/j.tics.2011.11.008
+22177031
+
+
+
+Ashburner
+J, Friston
+KJ (2005) Unified segmentation.
+
+26:839-851. 10.1016/j.neuroimage.2005.02.018
+15955494
+
+
+
+Barnett
+V, Lewis
+T (1978) . New York: Wiley.
+
+
+
+Baskin-Sommers
+AR, Curtin
+JJ, Newman
+JP (2011) Specifying the attentional selection that moderates the fearlessness of psychopathic offenders.
+
+22:226-234. 10.1177/0956797610396227
+21245494
+
+
+
+Blair
+RJR (2010) Neuroimaging of psychopathy and antisocial behavior: a targeted review.
+
+12:76-82. 10.1007/s11920-009-0086-x
+20425314
+
+
+
+Blair
+RJR (2013) Psychopathy: cognitive and neural dysfunction.
+15:181–190.24174892
+
+
+
+Blair
+RJR, Mitchell
+D, Blair
+K (2005) . Oxford, UK: Blackwell.
+
+
+
+Bos
+PA, Panksepp
+J, Bluthé
+RM, van Honk
+J (2012) Acute effects of steroid hormones and neuropeptides on human social-emotional behavior: a review of single administration studies.
+
+33:17-35. 10.1016/j.yfrne.2011.01.002
+21256859
+
+
+
+Carré
+JM, Hyde
+LW, Neumann
+CS, Viding
+E, Hariri
+AR (2013) The neural signatures of distinct psychopathic traits.
+
+8:122-135. 10.1080/17470919.2012.703623
+22775289
+
+
+
+Chen
+M, Bargh
+JA (1999) Consequences of automatic evaluation: immediate behavioral predispositions to approach or avoid the stimulus.
+
+25:215-224. 10.1177/0146167299025002007
+
+
+
+Cornell
+DG, Warren
+J, Hawk
+G, Stafford
+E, Oram
+G, Pine
+D (1996) Psychopathy in instrumental and reactive violent offenders.
+
+64:783-790. 8803369
+
+
+
+Craig
+MC, Catani
+M, Deeley
+Q, Latham
+R, Daly
+E, Kanaan
+R, Picchioni
+M, McGuire
+PK, Fahy
+T, Murphy
+DGM (2009) Altered connections on the road to psychopathy.
+
+14:946-953. 10.1038/mp.2009.40
+19506560
+
+
+
+Decety
+J, Skelly
+LR, Kiehl
+KA (2013) Brain response to empathy-eliciting scenarios involving pain in incarcerated individuals with psychopathy.
+
+70:638-645. 10.1001/jamapsychiatry.2013.27
+23615636
+
+
+
+Decety
+J, Chen
+C, Harenski
+CL, Kiehl
+KA (2015) Socioemotional processing of morally-laden behavior and their consequences on others in forensic psychopaths.
+
+36:2015-2026. 10.1002/hbm.22752
+25641358
+
+
+
+Dolan
+M, Anderson
+IM, Deakin
+JF (2001) Relationship between 5-HT function and impulsivity and aggression in male offenders with personality disorders.
+178:352-359.11282815
+
+
+
+Duvernoy
+HM, Cabanis
+EA, Vannson
+JL (1991) . Vienna: Springer.
+
+
+
+Eickhoff
+SB, Stephan
+KE, Mohlberg
+H, Grefkes
+C, Fink
+GR, Amunts
+K, Zilles
+K (2005) A new SPM toolbox for combining probabilistic cytoarchitectonic maps and functional imaging data.
+
+25:1325-1335. 10.1016/j.neuroimage.2004.12.034
+15850749
+
+
+
+Eisenegger
+C, Haushofer
+J, Fehr
+E (2011) The role of testosterone in social interaction.
+
+15:263-271. 10.1016/j.tics.2011.04.008
+21616702
+
+
+
+Ekman
+P, Friesen
+WV (1976) . Palo Alto, CA: Consulting Psychologist.
+
+
+
+Enter
+D, Spinhoven
+P, Roelofs
+K (2014) Alleviating social avoidance: effects of single dose testosterone administration on approach-avoidance action.
+
+65:351-354. 10.1016/j.yhbeh.2014.02.001
+24530652
+
+
+
+Fouche
+JP, van Der Wee
+NJ, Roelofs
+K, Stein
+DJ (2013) Recent advances in the brain imaging of social anxiety disorder.
+
+28:102-105. 10.1002/hup.2281
+23239106
+
+
+
+Freitas-Ferrari
+MC. (2010) Neuroimaging in social anxiety disorder: a systematic review of the literature.
+
+34:565-580. 10.1016/j.pnpbp.2010.02.028
+20206659
+
+
+
+Friston
+KJ (1997) Testing for anatomically specified regional effects.
+
+5:133-136. 10096418
+
+
+
+Friston
+KJ, Buechel
+C, Fink
+GR, Morris
+J, Rolls
+E, Dolan
+RJ (1997) Psychophysiological and modulatory interactions in neuroimaging.
+
+6:218-229. 10.1006/nimg.1997.0291
+9344826
+
+
+
+Gao
+Y, Raine
+A (2010) Successful and unsuccessful psychopaths: a neurobiological model.
+
+28:194-210. 10.1002/bsl.924
+20422645
+
+
+
+Giltay
+EJ, Enter
+D, Zitman
+FG, Penninx
+BW, van Pelt
+J, Spinhoven
+P, Roelofs
+K (2012) Salivary testosterone: associations with depression, anxiety disorders, and antidepressant use in a large cohort study.
+
+72:205-213. 10.1016/j.jpsychores.2011.11.01422325700
+
+
+
+Glenn
+AL, Raine
+A, Yaralian
+PS, Yang
+Y (2010) Increased volume of the striatum in psychopathic individuals.
+
+67:52-58. 10.1016/j.biopsych.2009.06.018
+19683706
+
+
+
+Grahn
+JA, Parkinson
+JA, Owen
+AM (2008) The cognitive functions of the caudate nucleus.
+
+86:141-155. 10.1016/j.pneurobio.2008.09.004
+18824075
+
+
+
+Groenestijn
+MAC, Akkerhuis
+GW, Kupka
+RW, Schneider
+N, Nolen
+WA (1999) . Lisse, The Netherlands: Swets Test.
+
+
+
+Haggard
+P (2008) Human volition: towards a neuroscience of will.
+
+9:934-946. 10.1038/nrn2497
+19020512
+
+
+
+Hare
+RD (2003) , Ed 2.
+Toronto, Canada: Multi-Health Systems.
+
+
+
+Hare
+RD, Hart
+SD, Harpur
+TJ (1991) Psychopathy and the DSM-IV criteria for antisocial personality disorder.
+
+100:391-398. 1918618
+
+
+
+Hildebrand
+M, De Ruiter
+C, Nijman
+H (2004) PCL-R psychopathy predicts disruptive behavior among male offenders in a Dutch forensic psychiatric hospital.
+
+19:13-29. 10.1177/088626050325904714680527
+
+
+
+Jubault
+T, Ody
+C, Koechlin
+E (2007) Serial organization of human behavior in the inferior parietal cortex.
+
+27:11028-11036. 10.1523/JNEUROSCI.1986-07.2007
+17928444
+
+
+
+Koenigs
+M, Baskin-Sommers
+A, Zeier
+J, Newman
+JP (2011) Investigating the neural correlates of psychopathy: a critical review.
+
+16:792-799. 10.1038/mp.2010.124
+21135855
+
+
+
+Lund
+TE, Nørgaard
+MD, Rostrup
+E, Rowe
+JB, Paulson
+OB (2005) Motion or activity: their role in intra- and inter-subject variation in fMRI.
+
+26:960-964. 10.1016/j.neuroimage.2005.02.021
+15955506
+
+
+
+Lundqvist
+D, Flykt
+A, Öhman
+A (1998) . Solna, Sweden: Department of Clinical Neuroscience, Psychology Section, Karolinska Institute.
+
+
+
+Maldjian
+JA, Laurienti
+PJ, Kraft
+RA, Burdette
+JH (2003) An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets.
+
+19:1233-1239. 12880848
+
+
+
+Malterer
+MB, Glass
+SJ, Newman
+JP (2008) Psychopathy and trait emotional intelligence.
+
+44:735-745. 10.1016/j.paid.2007.10.007
+18438451
+
+
+
+Marsh
+AA, Blair
+RJ (2008) Deficits in facial affect recognition among antisocial populations: a meta-analysis.
+
+32:454-465. 10.1016/j.neubiorev.2007.08.003
+17915324
+
+
+
+Marsh
+AA, Cardinale
+EM (2012) When psychopathy impairs moral judgments: neural responses during judgments about causing fear.
+
+9:3-11
+10.1093/scan/nss097
+22956667
+
+
+
+Martinez
+AM, Benavente
+R (1998) The AR face database. Barcelona, Spain: Computer Vision Center, Universidad Autònoma de Barcelona 24.
+
+
+
+Matsumoto
+D, Ekman
+P (1988) . San Francisco: University of California, Human Interaction Laboratory.
+
+
+
+Mehta
+PH, Beer
+J (2010) Neural mechanisms of the testosterone-aggression relation: the role of orbitofrontal cortex.
+
+22:2357-2368. 10.1162/jocn.2009.21389
+19925198
+
+
+
+Montoya
+ER, Terburg
+D, Bos
+PA, van Honk
+J (2012) Testosterone, cortisol, and serotonin as key regulators of social aggression: a review and theoretical perspective.
+
+36:65-73. 10.1007/s11031-011-9264-3
+22448079
+
+
+
+Motzkin
+JC, Newman
+JP, Kiehl
+KA, Koenigs
+M (2011) Reduced prefrontal connectivity in psychopathy.
+
+31:17348-17357. 10.1523/JNEUROSCI.4215-11.2011
+22131397
+
+
+
+Moul
+C, Killcross
+S, Dadds
+MR (2012) A model of differential amygdala activation in psychopathy.
+
+119:789-806. 10.1037/a0029342
+22800411
+
+
+
+Oldfield
+RC (1971) The assessment and analysis of handedness: the Edinburgh inventory.
+
+9:97-113. 5146491
+
+
+
+Patrick
+CJ, Hicks
+BM, Krueger
+RF, Lang
+AR (2005) Relations between psychopathy facets and externalizing in a criminal offender sample.
+19:339-356.16178678
+
+
+
+Poser
+BA, Versluis
+MJ, Hoogduin
+JM, Norris
+DG (2006) BOLD contrast sensitivity enhancement and artifact reduction with multiecho EPI: parallel-acquired inhomogeneity-desensitized fMRI.
+55:1227-1235.
+16680688
+
+
+
+Peper
+JS, van den Heuvel
+MP, Mandl
+RC, Pol
+HE, van Honk
+J (2011) Sex steroids and connectivity in the human brain: a review of neuroimaging studies.
+
+36:1101-1113. 10.1016/j.psyneuen.2011.05.004
+21641727
+
+
+
+Price
+CJ, Friston
+KJ (1999) Scanning patients with tasks they can perform.
+
+8:102-108. 10524600
+
+
+
+Quirk
+GJ, Gehlert
+DR (2003) Inhibition of the amygdala: key to pathological states?
+
+985:263-272. 12724164
+
+
+
+Radke
+S, Volman
+I, Mehta
+P, van Son
+V, Enter
+D, Sanfey
+A, Toni
+I, de Bruijn
+ERA, Roelofs
+K (2015) Testosterone biases the amygdala toward social threat approach.
+
+1:e1400074. 10.1126/sciadv.1400074
+26601187
+
+
+
+Ramnani
+N, Owen
+AM (2004) Anterior prefrontal cortex: insights into function from anatomy and neuroimaging.
+
+5(3):184-194. 10.1038/nrn1343
+14976518
+
+
+
+Rasmussen
+K, Storsaeter
+O, Levander
+S (1999) Personality disorders, psychopathy, and crime in a Norwegian prison population.
+
+22:91-97. 10086294
+
+
+
+Roelofs
+K, Minelli
+A, Mars
+RB, van Peer
+J, Toni
+I (2009) On the neural control of social emotional behavior.
+
+4:50-58. 10.1093/scan/nsn036
+19047074
+
+
+
+Schmand
+B, Bakker
+D, Saan
+R, Louman
+J (1991) The Dutch Reading Test for Adults: a measure of premorbid intelligence level.
+22:15-19.1877068
+
+
+
+Sommer
+M, Sodian
+B, Döhnel
+K, Schwerdtner
+J, Meinhardt
+J, Hajak
+G (2010) In psychopathic patients emotion attribution modulates activity in outcome-related brain areas.
+
+182:88-95. 10.1016/j.pscychresns.2010.01.007
+20417065
+
+
+
+Spielberger
+CD (1983) . Palo Alto, CA: Consulting Psychologists.
+
+
+
+Stålenheim
+EG, Eriksson
+E, von Knorring
+L, Wide
+L (1998) Testosterone as a biological marker in psychopathy and alcoholism.
+
+77:79-88. 9541143
+
+
+
+Stephan
+KE, Penny
+WD, Moran
+RJ, den Ouden
+HEM, Daunizeau
+J, Friston
+KJ (2010) Ten simple rules for dynamic causal modeling.
+
+49:3099-3109. 10.1016/j.neuroimage.2009.11.015
+19914382
+
+
+
+Stevens
+J (1996) , Ed 3
+Mahwah, NJ: Lawrence Erlbaum Associates.
+
+
+
+Tzourio-Mazoyer
+N, Landeau
+B, Papathanassiou
+D, Crivello
+F, Etard
+O, Delcroix
+N, Mazoyer
+B, Joliot
+M (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain.
+
+15:273-289. 10.1006/nimg.2001.0978
+11771995
+
+
+
+Van Vliet
+IM, Leroy
+H, Van Megen
+HJM (2000) . Leiden, The Netherlands: LUMC.
+
+
+
+Van Wingen
+G, Ossewaarde
+L, Bäckström
+T, Hermans
+EJ, Fernández
+G (2011) Gonadal hormone regulation of the emotion circuitry in humans.
+
+191:38-45. 10.1016/j.neuroscience.2011.04.042
+21540080
+
+
+
+Verhagen
+L, Grol
+MJ, Dijkerman
+HC, Toni
+I (2006) Studying visually-guided reach-to-grasp movements in an MR-environment.
+31:S45.
+
+
+
+Volman
+I, Roelofs
+K, Koch
+S, Verhagen
+L, Toni
+I (2011a) Anterior prefrontal cortex inhibition impairs control over social emotional actions.
+
+21:1766-1770. 10.1016/j.cub.2011.08.050
+22000109
+
+
+
+Volman
+I, Toni
+I, Verhagen
+L, Roelofs
+K (2011b) Endogenous testosterone modulates prefrontal-amygdala connectivity during social emotional behavior.
+21:2282-2290.21339377
+
+
+
+Volman
+I, Verhagen
+L, Den Ouden
+HEM, Fernandez
+G, Rijpkema
+M, Franke
+B, Toni
+I, Roelofs
+R (2013) Reduced serotonin transporter availability decreases prefrontal control of the amygdala.
+
+33:8974-8979. 10.1523/JNEUROSCI.5518-12.201323699508
+
+
+
+Von Borries
+AKL, Volman
+I, De Bruijn
+ERA, Bulten
+BH, Verkes
+RJ, Roelofs
+K (2012) Psychopaths lack the automatic avoidance of social threat: relation to instrumental aggression.
+
+200:761-766. 10.1016/j.psychres.2012.06.026
+22819277
+
+
+
+Von Der Heide
+RJ, Skipper
+LM, Klobusicky
+E, Olson
+IR (2013) Dissecting the uncinate fasciculus: disorders, controversies and a hypothesis.
+
+136:1692-1707. 10.1093/brain/awt094
+23649697
+
+
+
+Walters
+GD (2003) Predicting institutional adjustment and recidivism with the psychopathy checklist factor scores: a meta-analysis.
+
+27:541-558. 14593797
+
+
+
+Worsley
+KJ, Marrett
+S, Neelin
+P, Vandal
+AC, Friston
+KJ, Evans
+AC (1996) A unified statistical approach for determining significant signals in images of cerebral activation.
+
+4:58-73. 10.1002/(SICI)1097-0193(1996)4:1<58::AID-HBM4>3.0.CO;2-O
+20408186
+
+
+
+
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000.csv b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000.csv
new file mode 100644
index 0000000..b4a7369
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000.csv
@@ -0,0 +1,7 @@
+Unnamed: 0,Psychopathic offenders(n = 15),HCs(n = 19),p value
+Age,37.8 (7.9),40.7 (10.3),0.368
+IQ,101 (10),102 (9),0.761
+Handedness,50.7 (81),59.2 (62),0.729
+PCL-R total,30.4 (3.5),,
+PCL-R F1,12.1 (2.6),,
+PCL-R F2,14.1 (2.3),,
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000_info.json b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000_info.json
new file mode 100644
index 0000000..4d4ae78
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_000_info.json
@@ -0,0 +1 @@
+{"table_id": "T1", "table_label": "Table 1:", "table_caption": "Demographical data", "table_foot": "Values are presented as the mean (SD), unless otherwise indicated. F1, factor 1; F2, factor 2.", "n_header_rows": 1, "table_data_file": "table_000.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001.csv b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001.csv
new file mode 100644
index 0000000..e75ba72
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001.csv
@@ -0,0 +1,10 @@
+Unnamed: 0_level_0,Psychopathic offenders,Psychopathic offenders,HCs,HCs
+Unnamed: 0_level_1,Approach,Avoid,Approach,Avoid
+Errors (%),,,,
+Happy,3.2 (0.9),8.9 (1.8),2.4 (0.8),7.7 (1.1)
+Neutral,6.1 (1.3),5.8 (1.1),7.1 (1.4),5.2 (1.0)
+Angry,10.1 (2.2),13.1 (2.1),9.6 (1.8),11.6 (1.8)
+RT (ms),,,,
+Happy,554 (25),625 (35),553 (23),603 (25)
+Neutral,666 (28),687 (31),639 (21),668 (24)
+Angry,630 (25),665 (33),620 (24),630 (23)
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001_info.json b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001_info.json
new file mode 100644
index 0000000..cb4eedf
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_001_info.json
@@ -0,0 +1 @@
+{"table_id": "T2", "table_label": "Table 2:", "table_caption": "RTs and error rates for each group and factor of the AA task", "table_foot": "Values are presented as the mean (SE).", "n_header_rows": 2, "table_data_file": "table_001.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002.csv b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002.csv
new file mode 100644
index 0000000..20751bb
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002.csv
@@ -0,0 +1,32 @@
+Anatomical region,Putative BA,Side,x,y,z,Voxels (n),p value,t value
+Whole-brain effects,,,,,,,,
+Congruency effect over groups,,,,,,,,
+Cuneus,18,R/L,−16,−96,24,634.0,<0.001,4.85
+SPL/Superior occipital gyrus,7/19,L,−30,−76,34,254.0,0.004,4.37
+IFG,45/47,L,−52,22,−10,223.0,0.007,4.37
+Cuneus,18,R,18,−98,16,190.0,0.016,4.52
+Congruency effect for psychopathy group,,,,,,,,
+SPL/superior occipital gyrus,7/19,R/L,8,−82,38,925.0,<0.001,4.98
+Angular gyrus,39/19,L,−30,−72,34,337.0,0.001,4.35
+Superior temporal gyrus,42,L,−32,−32,6,214.0,0.009,4.51
+SPL,7,R,28,−74,46,158.0,0.034,4.58
+Cerebellum,,L,−24,−66,−38,146.0,0.046,4.63
+Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction,,,,,,,,
+aPFC,10,R,30,58,12,391.0,<0.001,5.1
+Supramarginal gyrus,40,R,54,−42,54,325.0,0.001,4.66
+Caudate nucleus,,R,10,10,2,273.0,0.002,4.69
+Putamen/Insula,,L,−34,6,−10,176.0,0.022,5.22
+Cerebellum,,R,18,−76,−38,188.0,0.016,5.09
+Negative testosterone modulation of Congruency effect in psychopathy,,,,,,,,
+Supramarginal gyrus,40,R,52,−40,54,657.0,<0.001,5.32
+Precentral/superior frontal gyrus,6,R/L,6,22,66,471.0,<0.001,5.65
+Caudate nucleus,,R,6,6,4,228.0,0.006,4.28
+VOI on bilateral aPFC,,,,,,,,
+Congruency effect over groups,10,R,30,58,14,31.0,0.001,4.46
+,10,L,−30,58,10,5.0,0.036,3.43
+Congruency effect in healthy control subjects,10,R,32,58,14,12.0,0.001,4.58
+,10,L,−34,52,4,2.0,0.040,3.4
+Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction,10,R,30,58,12,145.0,<0.001,5.1
+Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction,10,L,−24,56,6,15.0,0.010,3.87
+Negative testosterone modulation of congruency effect in psychopathy,10,R,32,56,10,77.0,0.002,4.34
+Negative testosterone modulation of congruency effect in psychopathy,10,L,−30,58,8,17.0,0.015,3.74
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002_info.json b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002_info.json
new file mode 100644
index 0000000..70d5cb8
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/table_002_info.json
@@ -0,0 +1 @@
+{"table_id": "T3", "table_label": "Table 3:", "table_caption": "Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)", "table_foot": "Coordinates are defined in MNI (x, y, z) space. The p values represent the FWE cluster-level corrected values for the whole-brain analyses and FWE voxel-level corrected values for the VOI analyses. IFG, Inferior frontal gyrus; L, left; R, right; SPL, superior parietal lobule.", "n_header_rows": 1, "table_data_file": "table_002.csv"}
\ No newline at end of file
diff --git a/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/tables.xml b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/tables.xml
new file mode 100644
index 0000000..7904b9d
--- /dev/null
+++ b/tests/data/sample_inputs/8EVW7TUtC9cx/source/pubget/tables/tables.xml
@@ -0,0 +1,2 @@
+
+4745181268780574745181eN-NWR-0107-1510.1523/ENEURO.0107-15.2016T1Table 1:Demographical dataValues are presented as the mean (SD), unless otherwise indicated. F1, factor 1; F2, factor 2.
Demographical data
Psychopathic offenders(n = 15)
HCs(n = 19)
p value
Age
37.8 (7.9)
40.7 (10.3)
0.368
IQ
101 (10)
102 (9)
0.761
Handedness
50.7 (81)
59.2 (62)
0.729
PCL-R total
30.4 (3.5)
PCL-R F1
12.1 (2.6)
PCL-R F2
14.1 (2.3)
Values are presented as the mean (SD), unless otherwise indicated. F1, factor 1; F2, factor 2.
Table 1:
Demographical data
Psychopathic offenders(n = 15)
HCs(n = 19)
p value
Age
37.8 (7.9)
40.7 (10.3)
0.368
IQ
101 (10)
102 (9)
0.761
Handedness
50.7 (81)
59.2 (62)
0.729
PCL-R total
30.4 (3.5)
PCL-R F1
12.1 (2.6)
PCL-R F2
14.1 (2.3)
Values are presented as the mean (SD), unless otherwise indicated. F1, factor 1; F2, factor 2.T2Table 2:RTs and error rates for each group and factor of the AA taskValues are presented as the mean (SE).
RTs and error rates for each group and factor of the AA task
Psychopathic offenders
HCs
Approach
Avoid
Approach
Avoid
Errors (%)
Happy
3.2 (0.9)
8.9 (1.8)
2.4 (0.8)
7.7 (1.1)
Neutral
6.1 (1.3)
5.8 (1.1)
7.1 (1.4)
5.2 (1.0)
Angry
10.1 (2.2)
13.1 (2.1)
9.6 (1.8)
11.6 (1.8)
RT (ms)
Happy
554 (25)
625 (35)
553 (23)
603 (25)
Neutral
666 (28)
687 (31)
639 (21)
668 (24)
Angry
630 (25)
665 (33)
620 (24)
630 (23)
Values are presented as the mean (SE).
Table 2:
RTs and error rates for each group and factor of the AA task
Psychopathic offenders
HCs
Approach
Avoid
Approach
Avoid
Errors (%)
Happy
3.2 (0.9)
8.9 (1.8)
2.4 (0.8)
7.7 (1.1)
Neutral
6.1 (1.3)
5.8 (1.1)
7.1 (1.4)
5.2 (1.0)
Angry
10.1 (2.2)
13.1 (2.1)
9.6 (1.8)
11.6 (1.8)
RT (ms)
Happy
554 (25)
625 (35)
553 (23)
603 (25)
Neutral
666 (28)
687 (31)
639 (21)
668 (24)
Angry
630 (25)
665 (33)
620 (24)
630 (23)
Values are presented as the mean (SE).T3Table 3:Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)Coordinates are defined in MNI (x, y, z) space. The p values represent the FWE cluster-level corrected values for the whole-brain analyses and FWE voxel-level corrected values for the VOI analyses. IFG, Inferior frontal gyrus; L, left; R, right; SPL, superior parietal lobule.
Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)
Anatomical region
Putative BA
Side
x
y
z
Voxels (n)
p value
t value
Whole-brain effects
Congruency effect over groups
Cuneus
18
R/L
−16
−96
24
634
<0.001
4.85
SPL/Superior occipital gyrus
7/19
L
−30
−76
34
254
0.004
4.37
IFG
45/47
L
−52
22
−10
223
0.007
4.37
Cuneus
18
R
18
−98
16
190
0.016
4.52
Congruency effect for psychopathy group
SPL/superior occipital gyrus
7/19
R/L
8
−82
38
925
<0.001
4.98
Angular gyrus
39/19
L
−30
−72
34
337
0.001
4.35
Superior temporal gyrus
42
L
−32
−32
6
214
0.009
4.51
SPL
7
R
28
−74
46
158
0.034
4.58
Cerebellum
L
−24
−66
−38
146
0.046
4.63
Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
aPFC
10
R
30
58
12
391
<0.001
5.10
Supramarginal gyrus
40
R
54
−42
54
325
0.001
4.66
Caudate nucleus
R
10
10
2
273
0.002
4.69
Putamen/Insula
L
−34
6
−10
176
0.022
5.22
Cerebellum
R
18
−76
−38
188
0.016
5.09
Negative testosterone modulation of Congruency effect in psychopathy
Supramarginal gyrus
40
R
52
−40
54
657
<0.001
5.32
Precentral/superior frontal gyrus
6
R/L
6
22
66
471
<0.001
5.65
Caudate nucleus
R
6
6
4
228
0.006
4.28
VOI on bilateral aPFC
Congruency effect over groups
10
R
30
58
14
31
0.001
4.46
10
L
−30
58
10
5
0.036
3.43
Congruency effect in healthy control subjects
10
R
32
58
14
12
0.001
4.58
10
L
−34
52
4
2
0.040
3.40
Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
10
R
30
58
12
145
<0.001
5.10
10
L
−24
56
6
15
0.010
3.87
Negative testosterone modulation of congruency effect in psychopathy
10
R
32
56
10
77
0.002
4.34
10
L
−30
58
8
17
0.015
3.74
Coordinates are defined in MNI (x, y, z) space. The p values represent the FWE cluster-level corrected values for the whole-brain analyses and FWE voxel-level corrected values for the VOI analyses. IFG, Inferior frontal gyrus; L, left; R, right; SPL, superior parietal lobule.
Table 3:
Clusters showing significantly larger activity for the affect-incongruent vs the affect-congruent conditions (emotion-control effect)
Anatomical region
Putative BA
Side
x
y
z
Voxels (n)
p value
t value
Whole-brain effects
Congruency effect over groups
Cuneus
18
R/L
−16
−96
24
634
<0.001
4.85
SPL/Superior occipital gyrus
7/19
L
−30
−76
34
254
0.004
4.37
IFG
45/47
L
−52
22
−10
223
0.007
4.37
Cuneus
18
R
18
−98
16
190
0.016
4.52
Congruency effect for psychopathy group
SPL/superior occipital gyrus
7/19
R/L
8
−82
38
925
<0.001
4.98
Angular gyrus
39/19
L
−30
−72
34
337
0.001
4.35
Superior temporal gyrus
42
L
−32
−32
6
214
0.009
4.51
SPL
7
R
28
−74
46
158
0.034
4.58
Cerebellum
L
−24
−66
−38
146
0.046
4.63
Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
aPFC
10
R
30
58
12
391
<0.001
5.10
Supramarginal gyrus
40
R
54
−42
54
325
0.001
4.66
Caudate nucleus
R
10
10
2
273
0.002
4.69
Putamen/Insula
L
−34
6
−10
176
0.022
5.22
Cerebellum
R
18
−76
−38
188
0.016
5.09
Negative testosterone modulation of Congruency effect in psychopathy
Supramarginal gyrus
40
R
52
−40
54
657
<0.001
5.32
Precentral/superior frontal gyrus
6
R/L
6
22
66
471
<0.001
5.65
Caudate nucleus
R
6
6
4
228
0.006
4.28
VOI on bilateral aPFC
Congruency effect over groups
10
R
30
58
14
31
0.001
4.46
10
L
−30
58
10
5
0.036
3.43
Congruency effect in healthy control subjects
10
R
32
58
14
12
0.001
4.58
10
L
−34
52
4
2
0.040
3.40
Negative testosterone modulation of group (psychopathic offenders > healthy control subjects) × congruency interaction
10
R
30
58
12
145
<0.001
5.10
10
L
−24
56
6
15
0.010
3.87
Negative testosterone modulation of congruency effect in psychopathy
10
R
32
56
10
77
0.002
4.34
10
L
−30
58
8
17
0.015
3.74
Coordinates are defined in MNI (x, y, z) space. The p values represent the FWE cluster-level corrected values for the whole-brain analyses and FWE voxel-level corrected values for the VOI analyses. IFG, Inferior frontal gyrus; L, left; R, right; SPL, superior parietal lobule.
\ No newline at end of file
diff --git a/tests/test_dataset.py b/tests/test_dataset.py
new file mode 100644
index 0000000..783e894
--- /dev/null
+++ b/tests/test_dataset.py
@@ -0,0 +1,7 @@
+from ns_pipelines.dataset import Dataset
+
+
+def test_dataset(sample_data):
+ dataset = Dataset(sample_data)
+
+ assert len(dataset) == 3
diff --git a/tests/test_participant_demographics.py b/tests/test_participant_demographics.py
new file mode 100644
index 0000000..8678bdc
--- /dev/null
+++ b/tests/test_participant_demographics.py
@@ -0,0 +1,20 @@
+import pytest
+from pathlib import Path
+
+from ns_pipelines import ParticipantDemographicsExtractor
+from ns_pipelines.dataset import Dataset
+
+
+@pytest.mark.vcr(record_mode="once", filter_headers=["authorization"])
+def test_ParticipantDemographicsExtractor(sample_data, tmp_path):
+ """Test the word count extraction pipeline."""
+ pde = ParticipantDemographicsExtractor(
+ extraction_model="gpt-4o-mini-2024-07-18",
+ prompt_set="ZERO_SHOT_MULTI_GROUP_FTSTRICT_FC",
+ env_variable="API_CLIENT_OPENAI_KEY",
+ env_file=str(Path(__file__).parents[1] / ".keys"),
+ )
+ dataset = Dataset(sample_data)
+ output_dir = tmp_path / "participant_demographics"
+ pde.run(dataset, output_dir)
+ assert True
diff --git a/tests/test_word_count.py b/tests/test_word_count.py
new file mode 100644
index 0000000..438dfb6
--- /dev/null
+++ b/tests/test_word_count.py
@@ -0,0 +1,33 @@
+from pathlib import Path
+from ns_pipelines import WordCountExtractor, WordDevianceExtractor
+from ns_pipelines.dataset import Dataset
+
+
+def test_WordCountExtractor(sample_data, tmp_path):
+ """Test the word count extraction pipeline."""
+ wce = WordCountExtractor()
+ dataset = Dataset(sample_data)
+ output_dir = tmp_path / "word_count"
+ wce.run(dataset, output_dir)
+ # rerun where the output directory already exists
+ # no ouputs generated
+ wce.run(dataset, output_dir)
+ # rerun with preference of ace
+ wce_ace = WordCountExtractor(input_sources=("ace", "pubget"))
+ wce_ace.run(dataset, output_dir)
+ assert True
+
+
+def test_WordDevianceExtractor(sample_data, tmp_path):
+ """Test the word deviance extraction pipeline."""
+ wde = WordDevianceExtractor()
+ dataset = Dataset(sample_data)
+ output_dir = tmp_path / "word_deviance"
+ wde.run(dataset, output_dir)
+ # rerun where the output directory already exists
+ # no ouputs generated
+ wde.run(dataset, output_dir)
+ # rerun with preference of ace
+ wde_ace = WordDevianceExtractor(input_sources=("ace", "pubget"))
+ wde_ace.run(dataset, output_dir)
+ assert True