Fighting or embracing multiplicity in neuroimaging? neighborhood leverage versus global calibration

Age of onset of obsessive-compulsive disorder differentially affects white matter microstructure

Previous diffusion MRI studies have reported mixed findings on white matter microstructure alterations in obsessive-compulsive disorder (OCD), likely due to variation in demographic and clinical characteristics, scanning methods, and underpowered samples. The OCD global study was created across five international sites to overcome these challenges by harmonizing data collection to identify consistent brain signatures of OCD that are reproducible and generalizable. Single-shell diffusion measures (e.g., fractional anisotropy), multi-shell Neurite Orientation Dispersion and Density Imaging (NODDI) and fixel-based measures, were extracted from skeletonized white matter tracts in 260 medication-free adults with OCD and 252 healthy controls. We additionally performed structural connectome analysis. We compared cases with controls and cases with early (<18) versus late (18+) OCD onset using mixed-model and Bayesian multilevel analysis. Compared with healthy controls, adult OCD individuals showed higher fiber density in the sagittal stratum (B[SE] = 0.10[0.05], P = 0.04) and credible evidence for higher fiber density in several other tracts. When comparing early (n = 145) and late-onset (n = 114) cases, converging evidence showed lower integrity of the posterior thalamic radiation -particularly radial diffusivity (B[SE] = 0.28[0.12], P = 0.03)-and lower global efficiency of the structural connectome (B[SE] = 15.3[6.6], P = 0.03) in late-onset cases. Post-hoc analyses indicated divergent direction of effects of the two OCD groups compared to healthy controls. Age of OCD onset differentially affects the integrity of thalamo-parietal/occipital tracts and the efficiency of the structural brain network. These results lend further support for the role of the thalamus and its afferent fibers and visual attentional processes in the pathophysiology of OCD.

Language access differentially alters functional connectivity during emotion perception across cultures

Introduction

It is often assumed that the ability to recognize the emotions of others is reflexive and automatic, driven only by observable facial muscle configurations. However, research suggests that accumulated emotion concept knowledge shapes the way people perceive the emotional meaning of others' facial muscle movements. Cultural upbringing can shape an individual's concept knowledge, such as expectations about which facial muscle configurations convey anger, disgust, or sadness. Additionally, growing evidence suggests that access to emotion category words, such as "anger," facilitates access to such emotion concept knowledge and in turn facilitates emotion perception.

Methods

To investigate the impact of cultural influence and emotion concept accessibility on emotion perception, participants from two cultural groups (Chinese and White Americans) completed a functional magnetic resonance imaging scanning session to assess functional connectivity between brain regions during emotion perception. Across four blocks, participants were primed with either English emotion category words ("anger," "disgust") or control text (XXXXXX) before viewing images of White American actors posing facial muscle configurations that are stereotypical of anger and disgust in the United States.

Results

We found that when primed with "disgust" versus control text prior to seeing disgusted facial expressions, Chinese participants showed a significant decrease in functional connectivity between a region associated with semantic retrieval (the inferior frontal gyrus) and regions associated with semantic processing, visual perception, and social cognition. Priming the word "anger" did not impact functional connectivity for Chinese participants relative to control text, and priming neither "disgust" nor "anger" impacted functional connectivity for White American participants.

Discussion

These findings provide preliminary evidence that emotion concept accessibility differentially impacts perception based on participants' cultural background.

Associations of prior concussion severity with brain microstructure using mean apparent propagator magnetic resonance imaging

Magnetic resonance imaging (MRI) diffusion studies have shown chronic microstructural tissue abnormalities in athletes with history of concussion, but with inconsistent findings. Concussions with post-traumatic amnesia (PTA) and/or loss of consciousness (LOC) have been connected to greater physiological injury. The novel mean apparent propagator (MAP) MRI is expected to be more sensitive to such tissue injury than the conventional diffusion tensor imaging. This study examined effects of prior concussion severity on microstructure with MAP-MRI. Collegiate-aged athletes (N = 111, 38 females; ≥6 months since most recent concussion, if present) completed semistructured interviews to determine the presence of prior concussion and associated injury characteristics, including PTA and LOC. MAP-MRI metrics (mean non-Gaussian diffusion [NG Mean], return-to-origin probability [RTOP], and mean square displacement [MSD]) were calculated from multi-shell diffusion data, then evaluated for associations with concussion severity through group comparisons in a primary model (athletes with/without prior concussion) and two secondary models (athletes with/without prior concussion with PTA and/or LOC, and athletes with/without prior concussion with LOC only). Bayesian multilevel modeling estimated models in regions of interest (ROI) in white matter and subcortical gray matter, separately. In gray matter, the primary model showed decreased NG Mean and RTOP in the bilateral pallidum and decreased NG Mean in the left putamen with prior concussion. In white matter, lower NG Mean with prior concussion was present in all ROI across all models and was further decreased with LOC. However, only prior concussion with LOC was associated with decreased RTOP and increased MSD across ROI. Exploratory analyses conducted separately in male and female athletes indicate associations in the primary model may differ by sex. Results suggest microstructural measures in gray matter are associated with a general history of concussion, while a severity-dependent association of prior concussion may exist in white matter.

Increased neural differentiation after a single session of aerobic exercise in older adults

Aging is associated with decreased cognitive function. One theory posits that this decline is in part due to multiple neural systems becoming dedifferentiated in older adults. Exercise is known to improve cognition in older adults, even after only a single session. We hypothesized that one mechanism of improvement is a redifferentiation of neural systems. We used a within-participant, cross-over design involving 2 sessions: either 30 minutes of aerobic exercise or 30 minutes of seated rest (n = 32; ages 55-81 years). Both functional Magnetic Resonance Imaging (fMRI) and Stroop performance were acquired soon after exercise and rest. We quantified neural differentiation via general heterogeneity regression. There were 3 prominent findings following the exercise. First, participants were better at reducing Stroop interference. Second, while there was greater neural differentiation within the hippocampal formation and cerebellum, there was lower neural differentiation within frontal cortices. Third, this greater neural differentiation in the cerebellum and temporal lobe was more pronounced in the older ages. These data suggest that exercise can induce greater neural differentiation in healthy aging.

Layer-specific, retinotopically-diffuse modulation in human visual cortex in response to viewing emotionally expressive faces

Viewing faces that are perceived as emotionally expressive evokes enhanced neural responses in multiple brain regions, a phenomenon thought to depend critically on the amygdala. This emotion-related modulation is evident even in primary visual cortex (V1), providing a potential neural substrate by which emotionally salient stimuli can affect perception. How does emotional valence information, computed in the amygdala, reach V1? Here we use high-resolution functional MRI to investigate the layer profile and retinotopic distribution of neural activity specific to emotional facial expressions. Across three experiments, human participants viewed centrally presented face stimuli varying in emotional expression and performed a gender judgment task. We found that facial valence sensitivity was evident only in superficial cortical layers and was not restricted to the retinotopic location of the stimuli, consistent with diffuse feedback-like projections from the amygdala. Together, our results provide a feedback mechanism by which the amygdala directly modulates activity at the earliest stage of visual processing.

The impact of neighborhood disadvantage on amygdala reactivity: Pathways through neighborhood social processes

Youth growing up in disadvantaged neighborhoods are more likely than their advantaged peers to face negative behavioral and mental health outcomes. Although studies have shown that adversity can undermine positive development via its impact on the developing brain, few studies have examined the association between neighborhood disadvantage and neural function, and no study has investigated potential social mechanisms within the neighborhood that might link neighborhood disadvantage to altered neural function. The current study evaluated the association between neighborhood disadvantage and amygdala reactivity during socioemotional face processing. We also assessed whether and which neighborhood-level social processes were related to amygdala reactivity, and whether these social processes mediated or moderated the association between neighborhood disadvantage and altered amygdala reactivity. We examined these aims in a registered report, using a sample of twins aged 7-19 years (N = 354 families, 708 twins) recruited from birth records with enrichment for neighborhood disadvantage. Twins completed a socioemotional face processing fMRI task and a sample of unrelated participants from the twins' neighborhoods were also recruited to serve as informants on neighborhood social processes. We found that neighborhood disadvantage was associated with greater right amygdala reactivity to threat, but only when neighborhood informants perceived norms in the neighborhood to be more permissive regarding general safety and management. The findings from this research add to the growing literature highlighting the influence of neighborhood disadvantage on amygdala function and the ways that supportive social processes may buffer the impact of adversity on brain function.

Reliability of task-evoked neural activation during face-emotion paradigms: Effects of scanner and psychological processes

Assessing and improving test–retest reliability is critical to efforts to address concerns about replicability of task‐based functional magnetic resonance imaging. The current study uses two statistical approaches to examine how scanner and task‐related factors influence reliability of neural response to face‐emotion viewing. Forty healthy adult participants completed two face‐emotion paradigms at up to three scanning sessions across two scanners of the same build over approximately 2 months. We examined reliability across the main task contrasts using Bayesian linear mixed‐effects models performed voxel‐wise across the brain. We also used a novel Bayesian hierarchical model across a predefined whole‐brain parcellation scheme and subcortical anatomical regions. Scanner differences accounted for minimal variance in temporal signal‐to‐noise ratio and task contrast maps. Regions activated during task at the group level showed higher reliability relative to regions not activated significantly at the group level. Greater reliability was found for contrasts involving conditions with clearly distinct visual stimuli and associated cognitive demands (e.g., face vs. nonface discrimination) compared to conditions with more similar demands (e.g., angry vs. happy face discrimination). Voxel‐wise reliability estimates tended to be higher than those based on predefined anatomical regions. This work informs attempts to improve reliability in the context of task activation patterns and specific task contrasts. Our study provides a new method to estimate reliability across a large number of regions of interest and can inform researchers' selection of task conditions and analytic contrasts.

Neural correlates of working memory and compensation at different stages of cognitive impairment in Parkinson’s disease

•

PD patients have two types of compensatory mechanisms for working memory (WM).

•

First, hyperactivation for different WM load tasks depending on cognitive status.

•

PD-CogNL has hyperactivation for moderate and heavy working memory load tasks.

•

PD-MCI and PDD have hyperactivation for control and light working memory load tasks.

•

Second, bilateral recruitment of WM-related areas improves WM performance.

Working memory (WM) impairment is one of the most frequent cognitive deficits in Parkinson's disease (PD). However, it is not known how neural activity is altered and compensatory responses eventually fail during progression. We aimed to elucidate neural correlates of WM and compensatory mechanisms in PD. Eighteen cognitively normal PD patients (PD-CogNL), 16 with PD with mild cognitive impairment (PD-MCI), 11 with PD with dementia (PDD), and 17 healthy controls (HCs) were evaluated. Subjects performed an n-back task. Functional MRI data were analyzed by event-related analysis for correct responses. Brain activations were evaluated by comparing them to fixation cross or 0-back task, and correlated with n-back task performance. When compared to fixation cross, PD-CogNL patients had more activation in WM areas than HCs for both the 2- and 3-back tasks. PD-MCI and PDD patients had more activation in WM areas than HCs for the 0- and 1-back task. 2-back task performance was correlated with brain activations (vs. 0-back task) in the bilateral dorsolateral prefrontal cortex and frontal eye field (FEF) and left rostral prefrontal cortex, caudate nucleus, inferior/superior parietal lobule (IPL/SPL), and anterior insular cortex as well as anterior cingulate cortex. 3-back task performance was correlated with brain activations (vs. 0-back task) in the left FEF, right caudate nucleus, and bilateral IPL/SPL. Additional activations on top of the 0-back task, rather than fixation cross, are the neural correlates of WM. Our results suggest PD patients have two types of compensatory mechanisms: (1) Hyperactivation for different WM load tasks depending on their cognitive status. PD-CogNL have hyperactivation for moderate and heavy working memory load tasks while maintaining normal working memory performance. In contrast, PD-MCI and PDD have hyperactivation for control task and light working memory load task, leaving less neural resources to further activate for more demanding tasks and resulting in impaired working memory performance. (2) Bilateral recruitment of WM-related areas, in particular the DLPFC, FEF, IPL/SPL and caudate nucleus, to improve WM performance.

Providing Evidence for the Null Hypothesis in Functional Magnetic Resonance Imaging Using Group-Level Bayesian Inference

Classical null hypothesis significance testing is limited to the rejection of the point-null hypothesis; it does not allow the interpretation of non-significant results. This leads to a bias against the null hypothesis. Herein, we discuss statistical approaches to ‘null effect’ assessment focusing on the Bayesian parameter inference (BPI). Although Bayesian methods have been theoretically elaborated and implemented in common neuroimaging software packages, they are not widely used for ‘null effect’ assessment. BPI considers the posterior probability of finding the effect within or outside the region of practical equivalence to the null value. It can be used to find both ‘activated/deactivated’ and ‘not activated’ voxels or to indicate that the obtained data are not sufficient using a single decision rule. It also allows to evaluate the data as the sample size increases and decide to stop the experiment if the obtained data are sufficient to make a confident inference. To demonstrate the advantages of using BPI for fMRI data group analysis, we compare it with classical null hypothesis significance testing on empirical data. We also use simulated data to show how BPI performs under different effect sizes, noise levels, noise distributions and sample sizes. Finally, we consider the problem of defining the region of practical equivalence for BPI and discuss possible applications of BPI in fMRI studies. To facilitate ‘null effect’ assessment for fMRI practitioners, we provide Statistical Parametric Mapping 12 based toolbox for Bayesian inference.

Beyond linearity in neuroimaging: Capturing nonlinear relationships with application to longitudinal studies

The ubiquitous adoption of linearity for quantitative predictors in statistical modeling is likely attributable to its advantages of straightforward interpretation and computational feasibility. The linearity assumption may be a reasonable approximation especially when the variable is confined within a narrow range, but it can be problematic when the variable's effect is non-monotonic or complex. Furthermore, visualization and model assessment of a linear fit are usually omitted because of challenges at the whole brain level in neuroimaging. By adopting a principle of learning from the data in the presence of uncertainty to resolve the problematic aspects of conventional polynomial fitting, we introduce a flexible and adaptive approach of multilevel smoothing splines (MSS) to capture any nonlinearity of a quantitative predictor for population-level neuroimaging data analysis. With no prior knowledge regarding the underlying relationship other than a parsimonious assumption about the extent of smoothness (e.g., no sharp corners), we express the unknown relationship with a sufficient number of smoothing splines and use the data to adaptively determine the specifics of the nonlinearity. In addition to introducing the theoretical framework of MSS as an efficient approach with a counterbalance between flexibility and stability, we strive to (a) lay out the specific schemes for population-level nonlinear analyses that may involve task (e.g., contrasting conditions) and subject-grouping (e.g., patients vs controls) factors; (b) provide modeling accommodations to adaptively reveal, estimate and compare any nonlinear effects of a predictor across the brain, or to more accurately account for the effects (including nonlinear effects) of a quantitative confound; (c) offer the associated program 3dMSS to the neuroimaging community for whole-brain voxel-wise analysis as part of the AFNI suite; and (d) demonstrate the modeling approach and visualization processes with a longitudinal dataset of structural MRI scans.

To pool or not to pool: Can we ignore cross-trial variability in FMRI?

In this work, we investigate the importance of explicitly accounting for cross-trial variability in neuroimaging data analysis. To attempt to obtain reliable estimates in a task-based experiment, each condition is usually repeated across many trials. The investigator may be interested in (a) condition-level effects, (b) trial-level effects, or (c) the association of trial-level effects with the corresponding behavior data. The typical strategy for condition-level modeling is to create one regressor per condition at the subject level with the underlying assumption that responses do not change across trials. In this methodology of complete pooling, all cross-trial variability is ignored and dismissed as random noise that is swept under the rug of model residuals. Unfortunately, this framework invalidates the generalizability from the confine of specific trials (e.g., particular faces) to the associated stimulus category ("face"), and may inflate the statistical evidence when the trial sample size is not large enough. Here we propose an adaptive and computationally tractable framework that meshes well with the current two-level pipeline and explicitly accounts for trial-by-trial variability. The trial-level effects are first estimated per subject through no pooling. To allow generalizing beyond the particular stimulus set employed, the cross-trial variability is modeled at the population level through partial pooling in a multilevel model, which permits accurate effect estimation and characterization. Alternatively, trial-level estimates can be used to investigate, for example, brain-behavior associations or correlations between brain regions. Furthermore, our approach allows appropriate accounting for serial correlation, handling outliers, adapting to data skew, and capturing nonlinear brain-behavior relationships. By applying a Bayesian multilevel model framework at the level of regions of interest to an experimental dataset, we show how multiple testing can be addressed and full results reported without arbitrary dichotomization. Our approach revealed important differences compared to the conventional method at the condition level, including how the latter can distort effect magnitude and precision. Notably, in some cases our approach led to increased statistical sensitivity. In summary, our proposed framework provides an effective strategy to capture trial-by-trial responses that should be of interest to a wide community of experimentalists.

Controllability over stressor decreases responses in key threat-related brain areas

Controllability over stressors has major impacts on brain and behavior. In humans, however, the effect of controllability on responses to stressors is poorly understood. Using functional magnetic resonance imaging (fMRI), we investigated how controllability altered responses to a shock-plus-sound stressor with a between-group yoked design, where participants in controllable and uncontrollable groups experienced matched stressor exposure. Employing Bayesian multilevel analysis at the level of regions of interest and voxels in the insula, and standard voxelwise analysis, we found that controllability decreased stressor-related responses across threat-related regions, notably in the bed nucleus of the stria terminalis and anterior insula. Posterior cingulate cortex, posterior insula, and possibly medial frontal gyrus showed increased responses during control over stressor. Our findings support the idea that the aversiveness of stressors is reduced when controllable, leading to decreased responses across key regions involved in anxiety-related processing, even at the level of the extended amygdala.

Untangling the relatedness among correlations, part III: Inter-subject correlation analysis through Bayesian multilevel modeling for naturalistic scanning

While inter-subject correlation (ISC) analysis is a powerful tool for naturalistic scanning data, drawing appropriate statistical inferences is difficult due to the daunting task of accounting for the intricate relatedness in data structure as well as handling the multiple testing issue. Although the linear mixed-effects (LME) modeling approach (Chen et al., 2017a) is capable of capturing the relatedness in the data and incorporating explanatory variables, there are a few challenging issues: 1) it is difficult to assign accurate degrees of freedom for each testing statistic, 2) multiple testing correction is potentially over-penalizing due to model inefficiency, and 3) thresholding necessitates arbitrary dichotomous decisions. Here we propose a Bayesian multilevel (BML) framework for ISC data analysis that integrates all regions of interest into one model. By loosely constraining the regions through a weakly informative prior, BML dissolves multiplicity through conservatively pooling the effect of each region toward the center and improves collective fitting and overall model performance. In addition to potentially achieving a higher inference efficiency, BML improves spatial specificity and easily allows the investigator to adopt a philosophy of full results reporting. A dataset of naturalistic scanning is utilized to illustrate the modeling approach with 268 parcels and to showcase the modeling capability, flexibility and advantages in results reporting. The associated program will be available as part of the AFNI suite for general use.

Brain function and clinical characterization in the Boston adolescent neuroimaging of depression and anxiety study

We present a Human Connectome Project study tailored toward adolescent anxiety and depression. This study is one of the first studies of the Connectomes Related to Human Diseases initiative and is collecting structural, functional, and diffusion-weighted brain imaging data from up to 225 adolescents (ages 14-17 years), 150 of whom are expected to have a current diagnosis of an anxiety and/or depressive disorder. Comprehensive clinical and neuropsychological evaluations and longitudinal clinical data are also being collected. This article provides an overview of task functional magnetic resonance imaging (fMRI) protocols and preliminary findings (N = 140), as well as clinical and neuropsychological characterization of adolescents. Data collection is ongoing for an additional 85 adolescents, most of whom are expected to have a diagnosis of an anxiety and/or depressive disorder. Data from the first 140 adolescents are projected for public release through the National Institutes of Health Data Archive (NDA) with the timing of this manuscript. All other data will be made publicly-available through the NDA at regularly scheduled intervals. This article is intended to serve as an introduction to this project as well as a reference for those seeking to clinical, neurocognitive, and task fMRI data from this public resource.