An integrative Bayesian approach to matrix-based analysis in neuroimaging

Structural connectivity of thalamic subnuclei in major depressive disorder: An ultra-high resolution diffusion MRI study at 7-Tesla

BACKGROUND

The thalamus serves as a central relay station within the brain, and thalamic connectional anomalies are increasingly thought to be present in major depressive disorder (MDD). However, the use of conventional MRI scanners and acquisition techniques has prevented a thorough examination of the thalamus and its subnuclear connectional profile. We combined ultra-high field diffusion MRI acquired at 7.0 Tesla to map the white matter connectivity of thalamic subnuclei.

METHODS

Fifty-three MDD patients and 12 healthy controls (HCs) were involved in the final analysis. FreeSurfer was used to segment the thalamic subnuclei, and MRtrix was used to perform the preprocessing and tractography. Fractional anisotropy, axial diffusivity, mean diffusivity, radial diffusivity, and streamline count of thalamic subnuclear tracts were measured as proxies of white matter microstructure. Bayesian multilevel model was used to assess group differences in white matter metrics for each thalamic subnuclear tract and the association between these white matter metrics and clinical features in MDD.

RESULTS

Evidence was found for reduced whiter matter metrics of the tracts spanning from all thalamic subnuclei among MDD versus HC participants. Moreover, evidence was found that white matter in various thalamic subnuclear tracts is related to medication status, age of onset and recurrence in MDD.

CONCLUSIONS

Structural connectivity was generally reduced in thalamic subnuclei in MDD participants. Several clinical characteristics are related to perturbed subnuclear thalamic connectivity with cortical and subcortical circuits that govern sensory processing, emotional function, and goal-directed behavior.

Sensitivity of Functional Arterial Spin Labelling in Detecting Cerebral Blood Flow Changes

Aims/Background Arterial spin labelling (ASL) is a non-invasive magnetic resonance imaging (MRI) method. ASL techniques can quantitatively measure cerebral perfusion by fitting a kinetic model to the difference between labelled images (tag images) and ones which are acquired without labelling (control images). ASL functional MRI (fMRI) provides quantitative perfusion maps by using arterial water as an endogenous tracer instead of depending on vascular blood oxygenation level.This study aimed to assess the number of pulsed ASL blocks that were needed to provide accurate and reliable regional estimates of cerebral blood flow (CBF) changes when participants engaged in visually guided saccade and fixation task; evaluate the localization to cortical control saccade versus fixation; investigate the relationship between the sensitivity of ASL fMRI and the number of blocks; and compare the sensitivity of blood oxygen level-dependent (BOLD) fMRI and ASL fMRI. Methods The experiment was a block-design paradigm consisting of two conditions: fixation and saccade. No response other than the eye movements of the participants was recorded during the scans. ASL and BOLD fMRI scans were conducted on all participants during the same session. The fMRI study consisted of two functional experiments: a CBF contrast was provided using the ASL sequence, and an optimized BOLD contrast was provided using the BOLD sequence. Results From group analysis in all divided blocks of ASL sessions (4, 6, 8...... 14, 16, 18......26, 28, 30), ASL yielded significant activation clusters in the visual cortex of the bilateral hemisphere from block 4. There was no false activation from block 4. No activation cluster was found by reversing analysis of block 2. Robust and consistent activation in the visual cortex was observed in each of the 14 divided blocks group analysis, and no activation was found in the eye field of the brain. The sensitivity of 4-block was found to be better than that of 8-block. More significant activation clusters of the visual cortex were found in BOLD than in ASL. No activation cluster of parietal eye field (PEF), frontal eye field (FEF) and supplementary eye field (SEF) was detected in ASL. The voxel size of the activation cluster increased with the increasing number of blocks, and the percent signal change in the activation cluster decreased with the escalating block number. The voxel size was positively correlated with the number of blocks (correlation coefficient = 0.98, p < 0.0001), and the percent signal change negatively correlated with the number of blocks (correlation coefficient = -0.90, p < 0.0001). Conclusion The 4-block pulsed functional ASL (fASL) presents accurate and reliable activation, with minimal time-on-task effect and little adverse impact of time, in participants engaging in visually guided saccade and fixation tasks. Despite having lower sensitivity than BOLD fMRI, ASL can determine accurate activation location. Although the time-on-task effects affect the observation for the sensitivity of ASL over task time, it is suggested that ASL fMRI may provide a powerful method for pinpointing the time-on-task effect over a long period of time.

How Distributed Subcortical Integration of Reward and Threat May Inform Subsequent Approach-Avoidance Decisions

Healthy and successful living involves carefully navigating rewarding and threatening situations by balancing approach and avoidance behaviors. Excessive avoidance to evade potential threats often leads to forfeiting potential rewards. However, little is known about how reward and threat information is integrated neurally to inform approach or avoidance. In this preregistered study, participants (N behavior = 31, 17F; N MRI = 28, 15F) made approach-avoidance decisions under varying reward (monetary gains) and threat (electrical stimulations) prospects during functional MRI scanning. In contrast to theorized parallel subcortical processing of reward and threat information before cortical integration, Bayesian multivariate multilevel analyses revealed subcortical reward and threat integration prior to indicating approach-avoidance decisions. This integration occurred in the ventral striatum, thalamus, and bed nucleus of the stria terminalis (BNST). When reward was low, risk-diminishing avoidance decisions dominated, which was linked to more positive tracking of threat magnitude prior to indicating avoidance than approach decisions. In contrast, the amygdala exhibited dual sensitivity to reward and threat. While anticipating outcomes of risky approach decisions, we observed positive tracking of threat magnitude within the salience network (dorsal anterior cingulate cortex, thalamus, periaqueductal gray, BNST). Conversely, after risk-diminishing avoidance, characterized by reduced reward prospects, we observed more negative tracking of reward magnitude in the ventromedial prefrontal cortex and ventral striatum. These findings shed light on the temporal dynamics of approach-avoidance decision-making. Importantly, they demonstrate the role of subcortical integration of reward and threat information in balancing approach and avoidance, challenging theories positing predominantly separate subcortical processing prior to cortical integration.

Estimation of fMRI responses related to epileptic discharges using Bayesian hierarchical modeling

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a unique and noninvasive method for epilepsy presurgical evaluation. When selecting voxels by null-hypothesis tests, the conventional analysis may overestimate fMRI response amplitudes related to interictal epileptic discharges (IEDs), especially when IEDs are rare. We aimed to estimate fMRI response amplitudes represented by blood oxygen level dependent (BOLD) percentage changes related to IEDs using a hierarchical model. It involves the local and distributed hemodynamic response homogeneity to regularize estimations. Bayesian inference was applied to fit the model. Eighty-two epilepsy patients who underwent EEG-fMRI and subsequent surgery were included in this study. A conventional voxel-wise general linear model was compared to the hierarchical model on estimated fMRI response amplitudes and on the concordance between the highest response cluster and the surgical cavity. The voxel-wise model overestimated fMRI responses compared to the hierarchical model, evidenced by a practically and statistically significant difference between the estimated BOLD percentage changes. Only the hierarchical model differentiated brief and long-lasting IEDs with significantly different BOLD percentage changes. Overall, the hierarchical model outperformed the voxel-wise model on presurgical evaluation, measured by higher prediction performance. When compared with a previous study, the hierarchical model showed higher performance metric values, but the same or lower sensitivity. Our results demonstrated the capability of the hierarchical model of providing more physiologically reasonable and more accurate estimations of fMRI response amplitudes induced by IEDs. To enhance the sensitivity of EEG-fMRI for presurgical evaluation, it may be necessary to incorporate more appropriate spatial priors and bespoke decision strategies.

Reduced vmPFC-insula functional connectivity in generalized anxiety disorder: a Bayesian confirmation study

Differences in the correlated activity of networked brain regions have been reported in individuals with generalized anxiety disorder (GAD) but an overreliance on null-hypothesis significance testing (NHST) limits the identification of disorder-relevant relationships. In this preregistered study, we applied both a Bayesian statistical framework and NHST to the analysis of resting-state fMRI scans from females with GAD and matched healthy comparison females. Eleven a-priori hypotheses about functional connectivity (FC) were evaluated using Bayesian (multilevel model) and frequentist (t-test) inference. Reduced FC between the ventromedial prefrontal cortex (vmPFC) and the posterior-mid insula (PMI) was confirmed by both statistical approaches and was associated with anxiety sensitivity. FC between the vmPFC-anterior insula, the amygdala-PMI, and the amygdala-dorsolateral prefrontal cortex (dlPFC) region pairs did not survive multiple comparison correction using the frequentist approach. However, the Bayesian model provided evidence for these region pairs having decreased FC in the GAD group. Leveraging Bayesian modeling, we demonstrate decreased FC of the vmPFC, insula, amygdala, and dlPFC in females with GAD. Exploiting the Bayesian framework revealed FC abnormalities between region pairs excluded by the frequentist analysis and other previously undescribed regions in GAD, demonstrating the value of applying this approach to resting-state FC data in clinical investigations.

Effects of schema on the relationship between post-encoding brain connectivity and subsequent durable memory

Schemas can facilitate memory consolidation. Studies have suggested that interactions between the hippocampus and the ventromedial prefrontal cortex (vmPFC) are important for schema-related memory consolidation. However, in humans, how schema accelerates the consolidation of new information and relates to durable memory remains unclear. To address these knowledge gaps, we used a human analogue of the rodent spatial schema task and resting-state fMRI to investigate how post-encoding brain networks can predict long-term memory performance in different schema conditions. After participants were trained to obtain schema-consistent or schema-inconsistent object-location associations, they learned new object-location associations. The new associations were tested after the post-encoding rest in the scanner and 24 h later outside the scanner. The Bayesian multilevel modelling was applied to analyse the post-encoding brain networks. The results showed that during the post-encoding, stronger vmPFC- anterior hippocampal connectivity was associated with durable memory in the schema-consistent condition, whereas stronger object-selective lateral occipital cortex (LOC)-ventromedial prefrontal connectivity and weaker connectivity inside the default mode network were associated with durable memory in the schema inconsistent condition. In addition, stronger LOC-anterior hippocampal connectivity was associated with memory in both schema conditions. These results shed light on how schemas reconfigure early brain networks, especially the prefrontal-hippocampal and stimuli-relevant cortical networks and influence long-term memory performance.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

Emotion dynamics as hierarchical Bayesian inference in time

What fundamental property of our environment would be most valuable and optimal in characterizing the emotional dynamics we experience in daily life? Empirical work has shown that an accurate estimation of uncertainty is necessary for our optimal perception, learning, and decision-making. However, the role of this uncertainty in governing our affective dynamics remains unexplored. Using Bayesian encoding, decoding and computational modeling, on a large-scale neuroimaging and behavioral data on a passive movie-watching task, we showed that emotions naturally arise due to ongoing uncertainty estimations about future outcomes in a hierarchical neural architecture. Several prefrontal subregions hierarchically encoded a lower-dimensional signal that highly correlated with the evolving uncertainty. Crucially, the lateral orbitofrontal cortex (lOFC) tracked the temporal fluctuations of this uncertainty and was predictive of the participants' predisposition to anxiety. Furthermore, we observed a distinct functional double-dissociation within OFC with increased connectivity between medial OFC and DMN, while with that of lOFC and FPN in response to the evolving affect. Finally, we uncovered a temporally predictive code updating an individual's beliefs spontaneously with fluctuating outcome uncertainty in the lOFC. A biologically relevant and computationally crucial parameter in the theories of brain function, we propose uncertainty to be central to the definition of complex emotions.

Neural similarity between mentalizing and live social interaction during the transition to adolescence

Social interactions are essential for human development, yet little neuroimaging research has examined their underlying neurocognitive mechanisms using socially interactive paradigms during childhood and adolescence. Recent neuroimaging research has revealed activity in the mentalizing network when children engage with a live social partner, even when mentalizing is not required. While this finding suggests that social‐interactive contexts may spontaneously engage mentalizing, it is not a direct test of how similarly the brain responds to these two contexts. The current study used representational similarity analysis on data from 8‐ to 14‐year‐olds who made mental and nonmental judgments about an abstract character and a live interaction partner during fMRI. A within‐subject, 2 (Mental/Nonmental) × 2 (Peer/Character) design enabled us to examine response pattern similarity between conditions, and estimate fit to three conceptual models of how the two contexts relate: (1) social interaction and mentalizing about an abstract character are represented similarly; (2) interactive peers and abstract characters are represented differently regardless of the evaluation type; and (3) mental and nonmental states are represented dissimilarly regardless of target. We found that the temporal poles represent mentalizing and peer interactions similarly (Model 1), suggesting a neurocognitive link between the two in these regions. Much of the rest of the social brain exhibits different representations of interactive peers and abstract characters (Model 2). Our findings highlight the importance of studying social‐cognitive processes using interactive approaches, and the utility of pattern‐based analyses for understanding how social‐cognitive processes relate to each other.

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.

A comprehensive macaque fMRI pipeline and hierarchical atlas

Functional neuroimaging research in the non-human primate (NHP) has been advancing at a remarkable rate. The increase in available data establishes a need for robust analysis pipelines designed for NHP neuroimaging and accompanying template spaces to standardize the localization of neuroimaging results. Our group recently developed the NIMH Macaque Template (NMT), a high-resolution population average anatomical template and associated neuroimaging resources, providing researchers with a standard space for macaque neuroimaging . Here, we release NMT v2, which includes both symmetric and asymmetric templates in stereotaxic orientation, with improvements in spatial contrast, processing efficiency, and segmentation. We also introduce the Cortical Hierarchy Atlas of the Rhesus Macaque (CHARM), a hierarchical parcellation of the macaque cerebral cortex with varying degrees of detail. These tools have been integrated into the neuroimaging analysis software AFNI to provide a comprehensive and robust pipeline for fMRI processing, visualization and analysis of NHP data. AFNI's new @animal_warper program can be used to efficiently align anatomical scans to the NMT v2 space, and afni_proc.py integrates these results with full fMRI processing using macaque-specific parameters: from motion correction through regression modeling. Taken together, the NMT v2 and AFNI represent an all-in-one package for macaque functional neuroimaging analysis, as demonstrated with available demos for both task and resting state fMRI.

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.

Cross-sectional exploration of brain functional connectivity in the triadic development model of adolescents

Adolescence represents a transitional stage with increased risk taking and mood dysregulation. These vulnerabilities are accountable by developmental dynamics in the triadic functional brain networks underlying reward seeking (REW), emotional avoidance (EMO), and cognitive regulation (COG). However, these triadic dynamics, though conceptually established, have yet been investigated directly. Capitalizing on public database of resting-state fMRI from 222 adolescents (8-18 years old, 89F133M), this study examined cross-sectional development profiles of functional connectivity (FC) by jointly considering bilateral seeds of the ventral striatum, amygdala, and dorsal lateral prefrontal cortex in probing the networks of REW, EMO, and COG, respectively. Positive and negative FCs were considered separately for clarification of synergetic and suppressive interactions. While the REW and EMO mostly exhibited quadratic FC changes across age, suggesting reduced reward sensitivity and risk avoidance, the COG exhibited both linear and quadratic FC changes, suggesting both protracted maturation of cognitive ability and lowered top-down regulation. Additional age × gender effects were identified in the precentral gyrus and superior medial prefrontal cortex, which may associate risky action and emotion dysregulation to boys and girls, respectively. These results provide network evidence in substantiating the "triadic model" and deepening existing insights into neurodevelopmental mechanisms associated with adolescent behavior.

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.

Theta-burst TMS to the posterior superior temporal sulcus decreases resting-state fMRI connectivity across the face processing network

Humans process faces by using a network of face-selective regions distributed across the brain. Neuropsychological patient studies demonstrate that focal damage to nodes in this network can impair face recognition, but such patients are rare. We approximated the effects of damage to the face network in neurologically normal human participants by using theta burst transcranial magnetic stimulation (TBS). Multi-echo functional magnetic resonance imaging (fMRI) resting-state data were collected pre- and post-TBS delivery over the face-selective right superior temporal sulcus (rpSTS), or a control site in the right motor cortex. Results showed that TBS delivered over the rpSTS reduced resting-state connectivity across the extended face processing network. This connectivity reduction was observed not only between the rpSTS and other face-selective areas, but also between nonstimulated face-selective areas across the ventral, medial, and lateral brain surfaces (e.g., between the right amygdala and bilateral fusiform face areas and occipital face areas). TBS delivered over the motor cortex did not produce significant changes in resting-state connectivity across the face processing network. These results demonstrate that, even without task-induced fMRI signal changes, disrupting a single node in a brain network can decrease the functional connectivity between nodes in that network that have not been directly stimulated.

Human behavior is dependent on brain networks that perform different cognitive functions. We combined theta burst transcranial magnetic stimulation (TBS) with resting-state fMRI to study the face processing network. Disruption of the face-selective right posterior superior temporal sulcus (rpSTS) reduced fMRI connectivity across the face network. This impairment in connectivity was observed not only between the rpSTS and other face-selective areas, but also between nonstimulated face-selective areas on the ventral and medial brain surfaces (e.g., between the right amygdala and bilateral fusiform face areas and occipital face areas). Thus, combined TBS/fMRI can be used to approximate and measure the effects of focal brain damage on brain networks, and suggests such an approach may be useful for mapping intrinsic network organization.

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

Neuroimaging faces the daunting challenge of multiple testing - an instance of multiplicity - that is associated with two other issues to some extent: low inference efficiency and poor reproducibility. Typically, the same statistical model is applied to each spatial unit independently in the approach of massively univariate modeling. In dealing with multiplicity, the general strategy employed in the field is the same regardless of the specifics: trust the local "unbiased" effect estimates while adjusting the extent of statistical evidence at the global level. However, in this approach, modeling efficiency is compromised because each spatial unit (e.g., voxel, region, matrix element) is treated as an isolated and independent entity during massively univariate modeling. In addition, the required step of multiple testing "correction" by taking into consideration spatial relatedness, or neighborhood leverage, can only partly recoup statistical efficiency, resulting in potentially excessive penalization as well as arbitrariness due to thresholding procedures. Moreover, the assigned statistical evidence at the global level heavily relies on the data space (whole brain or a small volume). The present paper reviews how Stein's paradox (1956) motivates a Bayesian multilevel (BML) approach that, rather than fighting multiplicity, embraces it to our advantage through a global calibration process among spatial units. Global calibration is accomplished via a Gaussian distribution for the cross-region effects whose properties are not a priori specified, but a posteriori determined by the data at hand through the BML model. Our framework therefore incorporates multiplicity as integral to the modeling structure, not a separate correction step. By turning multiplicity into a strength, we aim to achieve five goals: 1) improve the model efficiency with a higher predictive accuracy, 2) control the errors of incorrect magnitude and incorrect sign, 3) validate each model relative to competing candidates, 4) reduce the reliance and sensitivity on the choice of data space, and 5) encourage full results reporting. Our modeling proposal reverberates with recent proposals to eliminate the dichotomization of statistical evidence ("significant" vs. "non-significant"), to improve the interpretability of study findings, as well as to promote reporting the full gamut of results (not only "significant" ones), thereby enhancing research transparency and reproducibility.

Functional and resting-state characterizations of a periventricular heterotopic nodule associated with epileptogenic activity

The object of this study was to extensively characterize a region of periventricular nodular heterotopia (PVNH) in an epilepsy patient to reveal its possible neurocognitive functional role(s). The authors used 3-T MRI approaches to exhaustively characterize a single, right hemisphere heterotopion in a high-functioning adult male with medically responsive epilepsy, which had manifested during late adolescence. The heterotopion proved to be spectroscopically consistent with a cortical-like composition and was interconnected with nearby ipsilateral cortical fundi, as revealed by fiber tractography (diffusion-weighted imaging) and resting-state functional connectivity MRI (rsfMRI). Moreover, the region of PVNH demonstrated two novel characterizations for a heterotopion. First, functional MRI (fMRI), as distinct from rsfMRI, showed that the heterotopion was significantly modulated while the patient watched animated video scenes of biological motion (i.e., cartoons). Second, rsfMRI, which demonstrated correlated brain activity during a task-negative state, uniquely showed directionality within an interconnected network, receiving positive path effects from patent cortical and cerebellar foci while outputting only negative path effects to specific brain foci.These findings are addressed in the context of the impact on noninvasive presurgical brain mapping strategies for adult and pediatric patient workups, as well as the impact of this study on an understanding of the functional cortical architecture underlying cognition from a neurodiversity and evolutionary perspective.

An integrative Bayesian approach to matrix-based analysis in neuroimaging

Understanding the correlation structure associated with brain regions is a central goal in neuroscience, as it informs about interregional relationships and network organization. Correlation structure can be conveniently captured in a matrix that indicates the relationships among brain regions, which could involve electroencephalogram sensors, electrophysiology recordings, calcium imaging data, or functional magnetic resonance imaging (FMRI) data-We call this type of analysis matrix-based analysis, or MBA. Although different methods have been developed to summarize such matrices across subjects, including univariate general linear models (GLMs), the available modeling strategies tend to disregard the interrelationships among the regions, leading to "inefficient" statistical inference. Here, we develop a Bayesian multilevel (BML) modeling framework that simultaneously integrates the analyses of all regions, region pairs (RPs), and subjects. In this approach, the intricate relationships across regions as well as across RPs are quantitatively characterized. The adoption of the Bayesian framework allows us to achieve three goals: (a) dissolve the multiple testing issue typically associated with seeking evidence for the effect of each RP under the conventional univariate GLM; (b) make inferences on effects that would be treated as "random" under the conventional linear mixed-effects framework; and (c) estimate the effect of each brain region in a manner that indexes their relative "importance". We demonstrate the BML methodology with an FMRI dataset involving a cognitive-emotional task and compare it to the conventional GLM approach in terms of model efficiency, performance, and inferences. The associated program MBA is available as part of the AFNI suite for general use.