fMRI Signal Research Articles

Learning Robust and Sparse Principal Components With the α-Divergence.

In this paper, novel robust principal component analysis (RPCA) methods are proposed to exploit the local structure of datasets. The proposed methods are derived by minimizing the α -divergence between the sample distribution and the Gaussian density model. The α- divergence is used in different frameworks to represent variants of RPCA approaches including orthogonal, non-orthogonal, and sparse methods. We show that the classical PCA is a special case of our proposed methods where the α- divergence is reduced to the Kullback-Leibler (KL) divergence. It is shown in simulations that the proposed approaches recover the underlying principal components (PCs) by down-weighting the importance of structured and unstructured outliers. Furthermore, using simulated data, it is shown that the proposed methods can be applied to fMRI signal recovery and Foreground-Background (FB) separation in video analysis. Results on real world problems of FB separation as well as image reconstruction are also provided.

Uncovering Hemispheric Asymmetry and Directed Oscillatory Brain-Heart Interplay in Anxiety Processing: An fMRI Study.

Brain-heart interactions (BHI) are critical for generating and processing emotions, including anxiety. Understanding specific neural correlates would be instrumental for greater comprehension and potential therapeutic interventions of anxiety disorders. While prior work has implicated the pontine structure as a central processor in cardiac regulation in anxiety, the distributed nature of anxiety processing across the cortex remains elusive. To address this, we performed a whole-brain-heart analysis using the full frequency directed transfer function to study resting-state spectral differences in BHI between high and low anxiety groups undergoing fMRI scans. Our findings revealed a hemispheric asymmetry in low-frequency interplay (0.05 Hz - 0.15 Hz) characterized by ascending BHI to the left insula and descending BHI from the right insula. Furthermore, we provide evidence supporting the "pacemaker hypothesis", highlighting the pons' function in regulating cardiac activity. Higher frequency interplay (0.2 Hz - 0.4Hz) demonstrate a preference for ascending interactions, particularly towards ventral prefrontal cortical activity in high anxiety groups, suggesting the heart's role in triggering a cognitive response to regulate anxiety. These findings highlight the impact of anxiety on BHI, contributing to a better understanding of its effect on the resting-state fMRI signal, with further implications for potential therapeutic interventions in treating anxiety disorders.

A hybrid spatiotemporal deep belief network and sparse representation-based framework reveals multilevel core functional components in decoding multitask fMRI signals.

Decoding human brain activity on various task-based functional brain imaging data is of great significance for uncovering the functioning mechanism of the human mind. Currently, most feature extraction model-based methods for brain state decoding are shallow machine learning models, which may struggle to capture complex and precise spatiotemporal patterns of brain activity from the highly noisy fMRI raw data. Moreover, although decoding models based on deep learning methods benefit from their multilayer structure that could extract spatiotemporal features at multiscale, the relatively large populations of fMRI datasets are indispensable, and the explainability of their results is elusive. To address the above problems, we proposed a computational framework based on hybrid spatiotemporal deep belief network and sparse representations to differentiate multitask fMRI (tfMRI) signals. Using a relatively small cohort of tfMRI data as a test bed, our framework can achieve an average classification accuracy of 97.86% and define the multilevel temporal and spatial patterns of multiple cognitive tasks. Intriguingly, our model can characterize the key components for differentiating the multitask fMRI signals. Overall, the proposed framework can identify the interpretable and discriminative fMRI composition patterns at multiple scales, offering an effective methodology for basic neuroscience and clinical research with relatively small cohorts.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training.

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

Novel inductively coupled ear-bars (ICEs) to enhance restored fMRI signal from susceptibility compensation in rats.

Functional magnetic resonance imaging faces inherent challenges when applied to deep-brain areas in rodents, e.g. entorhinal cortex, due to the signal loss near the ear cavities induced by susceptibility artifacts and reduced sensitivity induced by the long distance from the surface array coil. Given the pivotal roles of deep brain regions in various diseases, optimized imaging techniques are needed. To mitigate susceptibility-induced signal losses, we introduced baby cream into the middle ear. To enhance the detection sensitivity of deep brain regions, we implemented inductively coupled ear-bars, resulting in approximately a 2-fold increase in sensitivity in entorhinal cortex. Notably, the inductively coupled ear-bar can be seamlessly integrated as an add-on device, without necessitating modifications to the scanner interface. To underscore the versatility of inductively coupled ear-bars, we conducted echo-planner imaging-based task functional magnetic resonance imaging in rats modeling Alzheimer's disease. As a proof of concept, we also demonstrated resting-state-functional magnetic resonance imaging connectivity maps originating from the left entorhinal cortex-a central hub for memory and navigation networks-to amygdala hippocampal area, Insular Cortex, Prelimbic Systems, Cingulate Cortex, Secondary Visual Cortex, and Motor Cortex. This work demonstrates an optimized procedure for acquiring large-scale networks emanating from a previously challenging seed region by conventional magnetic resonance imaging detectors, thereby facilitating improved observation of functional magnetic resonance imaging outcomes.

Breath-hold BOLD fMRI without CO2 sampling enables estimation of venous cerebral blood volume: potential use in normalization of stimulus-evoked BOLD fMRI data

BOLD fMRI signal has been used in conjunction with vasodilatory stimulation as a marker of cerebrovascular reactivity (CVR): the relative change in cerebral blood flow (CBF) arising from a unit change in the vasodilatory stimulus. Using numerical simulations, we demonstrate that the variability in the relative BOLD signal change induced by vasodilation is strongly influenced by the variability in deoxyhemoglobin-containing cerebral blood volume (CBV), as this source of variability is likely to be more prominent than that of CVR. It may, therefore, be more appropriate to describe the relative BOLD signal change induced by an isometabolic vasodilation as a proxy of deoxygenated CBV (CBVdHb) rather than CVR. With this in mind, a new method was implemented to map a marker of CBVdHb, termed BOLD-CBV, based on the normalization of voxel-wise BOLD signal variation by an estimate of the intravascular venous BOLD signal from voxels filled with venous blood. The intravascular venous BOLD signal variation, recorded during repeated breath-holding, was extracted from the superior sagittal sinus in a cohort of 27 healthy volunteers and used as a regressor across the whole brain, yielding maps of BOLD-CBV. In the same cohort, we demonstrated the potential use of BOLD-CBV for the normalization of stimulus-evoked BOLD fMRI by comparing group-level BOLD fMRI responses to a visuomotor learning task with and without the inclusion of voxel-wise vascular covariates of BOLD-CBV and the BOLD signal change per mmHg variation in end-tidal carbon dioxide (BOLD-CVR). The empirical measure of BOLD-CBV accounted for more between-subject variability in the motor task-induced BOLD responses than BOLD-CVR estimated from end-tidal carbon dioxide recordings. The new method can potentially increase the power of group fMRI studies by including a measure of vascular characteristics and has the strong practical advantage of not requiring experimental measurement of end-tidal carbon dioxide, unlike traditional methods to estimate BOLD-CVR. It also more closely represents a specific physiological characteristic of brain vasculature than BOLD-CVR, namely blood volume.

FMRI-based neurofeedback strategies and the way forward to treating phasic psychiatric symptoms.

Auditory verbal hallucinations (AVH) are the perfect illustration of phasic symptoms in psychiatric disorders. For some patients and in some situations, AVH cannot be relieved by standard therapeutic approaches. More advanced treatments are needed, among which neurofeedback, and more specifically fMRI-based neurofeedback, has been considered. This paper discusses the different possibilities to approach neurofeedback in the specific context of phasic symptoms, by highlighting the strengths and weaknesses of the available neurofeedback options. It concludes with the added value of the recently introduced information-based neurofeedback. Although requiring an online fMRI signal classifier, which can be quite complex to implement, this neurofeedback strategy opens a door toward an alternative treatment option for complex phasic symptomatology.

High-temporal resolution functional PET/MRI reveals coupling between human metabolic and hemodynamic brain response

PurposePositron emission tomography (PET) provides precise molecular information on physiological processes, but its low temporal resolution is a major obstacle. Consequently, we characterized the metabolic response of the human brain to working memory performance using an optimized functional PET (fPET) framework at a temporal resolution of 3 s.MethodsThirty-five healthy volunteers underwent fPET with [18F]FDG bolus plus constant infusion, 19 of those at a hybrid PET/MRI scanner. During the scan, an n-back working memory paradigm was completed. fPET data were reconstructed to 3 s temporal resolution and processed with a novel sliding window filter to increase signal to noise ratio. BOLD fMRI signals were acquired at 2 s.ResultsConsistent with simulated kinetic modeling, we observed a constant increase in the [18F]FDG signal during task execution, followed by a rapid return to baseline after stimulation ceased. These task-specific changes were robustly observed in brain regions involved in working memory processing. The simultaneous acquisition of BOLD fMRI revealed that the temporal coupling between hemodynamic and metabolic signals in the primary motor cortex was related to individual behavioral performance during working memory. Furthermore, task-induced BOLD deactivations in the posteromedial default mode network were accompanied by distinct temporal patterns in glucose metabolism, which were dependent on the metabolic demands of the corresponding task-positive networks.ConclusionsIn sum, the proposed approach enables the advancement from parallel to truly synchronized investigation of metabolic and hemodynamic responses during cognitive processing. This allows to capture unique information in the temporal domain, which is not accessible to conventional PET imaging.

ASD-GResTM: Deep Learning Framework for ASD classification using Gramian Angular Field.

Autism Spectrum Disorder (ASD) is a heterogeneous disorder in children, and the current clinical diagnosis is accomplished using behavioral, cognitive, developmental, and language metrics. These clinical metrics can be imperfect measures as they are subject to high test-retest variability, and are influenced by assessment factors such as environment, social structure, or comorbid disorders. Advances in neuroimaging coupled with machine-learning provides an opportunity to develop methods that are more quantifiable, and reliable than existing clinical techniques. In this paper, we design and develop a deep-learning model that operates on functional magnetic resonance imaging (fMRI) data, and can classify between ASD and neurotypical brains. We introduce a novel strategy to transform time-series data extracted from fMRI signals into Gramian Angular Field (GAF) while locking in the temporal and spatial patterns in the data. Our motivation is to design and develop a novel framework that could encode the time-series, acquired from fMRI data, into images that can be used by deep-learning architectures that have been successful in computer vision. In our proposed framework called ASD-GResTM, we used a Convolutional Neural Network (CNN) to extract useful features from GAF images. We then used a Long Short-Term Memory (LSTM) layer to learn the activities between the regions. Finally, the output representations of the last LSTM layer are applied to a single-layer perceptron (SPL) to get the final classification. Our extensive experimentation demonstrates high accuracy across 4 centers, and outperforms state-of-the-art models on two centers with an increase in the accuracy of 17.58% and 6.7%, respectively as compared to the state of the art. Our model achieved the maximum accuracy of 81.78% with high degree of sensitivity and specificity. All training, validation, and testing was accomplished using openly available ABIDE-I benchmarking dataset.

Iron content affects age group differences in associative learning-related fMRI activity

Brain regions accumulate different amounts of iron with age, with older adults having higher iron in the basal ganglia (globus pallidus, putamen, caudate) relative to the hippocampus. This has important implications for functional magnetic resonance imaging (fMRI) studies in aging as the presence of iron may influence both neuronal functioning as well as the measured fMRI (BOLD) signal, and these effects will vary across age groups and brain regions. To test this hypothesis, the current study examined the effect of iron on age group differences in task-related activity within each basal nuclei and the hippocampus. Twenty-eight younger and 22 older adults completed an associative learning task during fMRI acquisition. Iron content (QSM, R2*) was estimated from a multi-echo gradient echo sequence. As previously reported, older adults learned significantly less than younger adults and age group differences in iron content were largest in the basal ganglia (putamen, caudate). In the hippocampus (early task stage) and globus pallidus (late task stage), older adults had significantly higher learning-related activity than younger adults both before and after controlling for iron. In the putamen (late task stage), however, younger adults had significantly higher learning-related activity than older adults that was only seen after controlling for iron. These findings support the notion that age-related differences in iron influence both neuronal functioning and the measured fMRI signal in select basal nuclei. Moreover, previous fMRI studies in aging populations may have under-reported age group differences in task-related activity by not accounting for iron within these regions.

The response time paradox in functional magnetic resonance imaging analyses.

Response times (RTs) are often the main signal of interest in cognitive psychology but are often ignored in functional MRI (fMRI) analyses. In fMRI analysis the intensity of the signal serves as a proxy for the intensity of local neuronal activity, but changes in either the intensity or the duration of neuronal activity can yield identical fMRI signals. Therefore, if RTs are ignored and pair with neuronal durations, fMRI results claiming intensity differences may be confounded by RTs. We show how ignoring RTs goes beyond this confound, where longer RTs are paired with larger activation estimates, to lesser-known issues where RTs become confounds in group-level analyses and, surprisingly, how the RT confound can induce other artificial group-level associations with variables that are not related to the condition contrast or RTs. We propose a new time-series model to address these issues and encourage increasing focus on what the widespread RT-based signal represents.

The impact of emotional valence and stimulus habituation on fMRI signal reliability during emotion generation

BackgroundThe emotional domain is often impaired across many neurological diseases, for this reason it represents a relevant target of rehabilitation interventions. Functional changes in neural activity related to treatment can be assessed with functional MRI (fMRI) using emotion-generation tasks in longitudinal settings. Previous studies demonstrated that within-subject fMRI signal reliability can be affected by several factors such as repetition suppression, type of task and brain anatomy. However, the differential role of repetition suppression and emotional valence of the stimuli on the fMRI signal reliability and reproducibility during an emotion-generation task involving the vision of emotional pictures is yet to be determined. MethodsSixty-two healthy subjects were enrolled and split into two groups: group A (21 subjects, test-retest reliability on same-day and with same-task-form), group B (30 subjects, test-retest reproducibility with 4-month-interval using two equivalent-parallel forms of the task). Test-retest reliability and reproducibility of fMRI responses and patterns were evaluated separately for positive and negative emotional valence conditions in both groups. The analyses were performed voxel-wise, using the general linear model (GLM), and via a region-of-interest (ROI)-based approach, by computing the intra-class correlation coefficient (ICC) on the obtained contrasts. ResultsThe voxel-wise GLM test yielded no significant differences for both conditions in reliability and reproducibility analyses. As to the ROI-based approach, across all areas with significant main effects of the stimuli, the reliability, as measured with ICC, was poor (<0.4) for the positive condition and ranged from poor to excellent (0.4–0.75) for the negative condition. The ICC-based reproducibility analysis, related to the comparison of two different parallel forms, yielded similar results. DiscussionThe voxel-wise GLM analysis failed to capture the poor reliability of fMRI signal which was instead highlighted using the ROI-based ICC analysis. The latter showed higher signal reliability for negative valence stimuli with respect to positive ones. The implementation of two parallel forms allowed to exclude neural suppression as the predominant effect causing low signal reliability, which could be instead ascribed to the employment of different neural strategies to cope with emotional stimuli over time. This is an invaluable information for a better assessment of treatment and rehabilitation effects in longitudinal studies of emotional neural processing.

Modulatory effects of fMRI acquisition time of day, week and year on adolescent functional connectomes across spatial scales: Implications for inference

Metabolic, hormonal, autonomic and physiological rhythms may have a significant impact on cerebral hemodynamics and intrinsic brain synchronization measured with fMRI (the resting-state connectome). The impact of their characteristic time scales (hourly, circadian, seasonal), and consequently scan timing effects, on brain topology in inherently heterogeneous developing connectomes remains elusive. In a cohort of 4102 early adolescents with resting-state fMRI (median age = 120.0 months; 53.1 % females) from the Adolescent Brain Cognitive Development Study, this study investigated associations between scan time-of-day, time-of-week (school day vs weekend) and time-of-year (school year vs summer vacation) and topological properties of resting-state connectomes at multiple spatial scales. On average, participants were scanned around 2 pm, primarily during school days (60.9 %), and during the school year (74.6 %). Scan time-of-day was negatively correlated with multiple whole-brain, network-specific and regional topological properties (with the exception of a positive correlation with modularity), primarily of visual, dorsal attention, salience, frontoparietal control networks, and the basal ganglia. Being scanned during the weekend (vs a school day) was correlated with topological differences in the hippocampus and temporoparietal networks. Being scanned during the summer vacation (vs the school year) was consistently positively associated with multiple topological properties of bilateral visual, and to a lesser extent somatomotor, dorsal attention and temporoparietal networks. Time parameter interactions suggested that being scanned during the weekend and summer vacation enhanced the positive effects of being scanned in the morning. Time-of-day effects were overall small but spatially extensive, and time-of-week and time-of-year effects varied from small to large (Cohen's f ≤ 0.1, Cohen's d<0.82, p < 0.05). Together, these parameters were also positively correlated with temporal fMRI signal variability but only in the left hemisphere. Finally, confounding effects of scan time parameters on relationships between connectome properties and cognitive task performance were assessed using the ABCD neurocognitive battery. Although most relationships were unaffected by scan time parameters, their combined inclusion eliminated associations between properties of visual and somatomotor networks and performance in the Matrix Reasoning and Pattern Comparison Processing Speed tasks. Thus, scan time of day, week and year may impact measurements of adolescent brain's functional circuits, and should be accounted for in studies on their associations with cognitive performance, in order to reduce the probability of incorrect inference.

Characterizing systemic physiological effects on the blood oxygen level dependent signal of resting-state fMRI in time-frequency space using wavelets.

Systemic physiological dynamics, such as heart rate variability (HRV) and respiration volume per time (RVT), are known to account for significant variance in the blood oxygen level dependent (BOLD) signal of resting-state functional magnetic resonance imaging (rsfMRI). However, synchrony between these cardiorespiratory changes and the BOLD signal could be due to neuronal (i.e., autonomic activity inducing changes in heart rate and respiration) or vascular (i.e., cardiorespiratory activity facilitating hemodynamic changes and thus the BOLD signal) effects and the contributions of these effects may differ spatially, temporally, and spectrally. In this study, we characterize these brain-body dynamics using a wavelet analysis in rapidly sampled rsfMRI data with simultaneous pulse oximetry and respiratory monitoring of the Human Connectome Project. Our time-frequency analysis across resting-state networks (RSNs) revealed differences in the coherence of the BOLD signal and heartbeat interval (HBI)/RVT dynamics across frequencies, with unique profiles per network. Somatomotor (SMN), visual (VN), and salience (VAN) networks demonstrated the greatest synchrony with both systemic physiological signals when compared to other networks; however, significant coherence was observed in all RSNs regardless of direct autonomic involvement. Our phase analysis revealed distinct frequency profiles of percentage of time with significant coherence between BOLD and systemic physiological signals for different phase offsets across RSNs, suggesting that the phase offset and temporal order of signals varies by frequency. Lastly, our analysis of temporal variability of coherence provides insight on potential influence of autonomic state on brain-body communication. Overall, the novel wavelet analysis enables an efficient characterization of the dynamic relationship between cardiorespiratory activity and the BOLD signal in spatial, temporal, and spectral dimensions to inform our understanding of autonomic states and improve our interpretation of the BOLD signal.

Functional connectivity of sensorimotor network is enhanced in spastic diplegic cerebral palsy: A multimodal study using fMRI and MEG

ObjectiveTo assess the effects to functional connectivity (FC) caused by lesions related to spastic diplegic cerebral palsy (CP) in children and adolescents using multiple imaging modalities. MethodsWe used resting state magnetoencephalography (MEG) envelope signals in alpha, beta and gamma ranges and resting state functional magnetic resonance imaging (fMRI) signals to quantify FC between selected sensorimotor regions of interest (ROIs) in 11 adolescents with spastic diplegic cerebral palsy and 24 typically developing controls. Motor performance of the hands was quantified with gross motor, fine motor and kinesthesia tests. ResultsIn fMRI, participants with CP showed enhanced FC within posterior parietal regions; in MEG, they showed enhanced interhemispheric FC between sensorimotor regions and posterior parietal regions both in alpha and lower beta bands. There was a correlation between the kinesthesia score and fronto-parietal connectivity in the control population. ConclusionsCP is associated with enhanced FC in sensorimotor network. This difference is not correlated with hand coordination performance. The effect of the lesion is likely not fully captured by temporal correlation of ROI signals. SignificanceBrain lesions can show as increased temporal correlation of activity between remote brain areas. We suggest this effect is likely separate from typical physiological correlates of functional connectivity.

Topological Data Analysis Captures Task-Driven fMRI Profiles in Individual Participants: A Classification Pipeline Based on Persistence.

BOLD-based fMRI is the most widely used method for studying brain function. The BOLD signal while valuable, is beset with unique vulnerabilities. The most notable of these is the modest signal to noise ratio, and the relatively low temporal and spatial resolution. However, the high dimensional complexity of the BOLD signal also presents unique opportunities for functional discovery. Topological Data Analyses (TDA), a branch of mathematics optimized to search for specific classes of structure within high dimensional data may provide particularly valuable applications. In this investigation, we acquired fMRI data in the anterior cingulate cortex (ACC) using a basic motor control paradigm. Then, for each participant and each of three task conditions, fMRI signals in the ACC were summarized using two methods: a) TDA based methods of persistent homology and persistence landscapes and b) non-TDA based methods using a standard vectorization scheme. Finally, using machine learning (with support vector classifiers), classification accuracy of TDA and non-TDA vectorized data was tested across participants. In each participant, TDA-based classification out-performed the non-TDA based counterpart, suggesting that our TDA analytic pipeline better characterized task- and condition-induced structure in fMRI data in the ACC. Our results emphasize the value of TDA in characterizing task- and condition-induced structure in regional fMRI signals. In addition to providing our analytical tools for other users to emulate, we also discuss the unique role that TDA-based methods can play in the study of individual differences in the structure of functional brain signals in the healthy and the clinical brain.

Mindfulness-oriented recovery enhancement in opioid use disorder: Extended emotional regulation and neural effects and immediate effects of guided meditation in a pilot sample

ObjectiveMindfulness-Oriented Recovery Enhancement (MORE) is an efficacious intervention to aid recovery from substance use disorder. This study in a pilot sample of individuals in treatment for opioid use disorder (OUD) characterizes longer-term changes after the MORE intervention and immediate effects of a brief MORE guided meditation session. DesignTwelve female participants in residential treatment for OUD completed an 8-week MORE intervention. Participants completed two sessions: one before and one after the 8-week MORE intervention. Each session included an emotional regulation questionnaire outside an MRI scanner first and then a 10-minute guided MORE meditation inside the scanner during which functional magnetic resonance imaging (fMRI) data were collected. Emotional regulation was measured after 8-weeks of MORE intervention. In addition, functional connectivity (i.e. correlated fMRI signal) between regions in a hypothesized affect regulation network was measured during the meditation state to assess change in brain network function due to 8-weeks of MORE. For each 10-min guided meditation, we also assessed their mood and opioid craving. ResultsNine participants completed all measurements. Participants’ emotional regulation difficulty significantly decreased after 8-weeks of MORE intervention. Furthermore, after 8-weeks of MORE, there was significantly increased connectivity between left ventromedial prefrontal cortex and left amygdala and between left ventrolateral prefrontal cortex and left nucleus accumbens captured during a meditation state. In both sessions, positive mood significantly increased after 10-min of guided mediation, however opioid craving was not significantly influenced. ConclusionsThis pilot study characterizes potential benefits of 8-week MORE intervention in improving emotional regulation difficulty and brain function. A 10-min guided MORE meditation may immediately improve mood, with potential to reduce acute stress- or cue-provoked craving. These results warrant future studies with larger sample size.

Auditory and reward structures reflect the pleasure of musical expectancies during naturalistic listening

Enjoying music consistently engages key structures of the neural auditory and reward systems such as the right superior temporal gyrus (R STG) and ventral striatum (VS). Expectations seem to play a central role in this effect, as preferences reliably vary according to listeners’ uncertainty about the musical future and surprise about the musical past. Accordingly, VS activity reflects the pleasure of musical surprise, and exhibits stronger correlations with R STG activity as pleasure grows. Yet the reward value of musical surprise – and thus the reason for these surprises engaging the reward system – remains an open question. Recent models of predictive neural processing and learning suggest that forming, testing, and updating hypotheses about one’s environment may be intrinsically rewarding, and that the constantly evolving structure of musical patterns could provide ample opportunity for this procedure. Consistent with these accounts, our group previously found that listeners tend to prefer melodic excerpts taken from real music when it either validates their uncertain melodic predictions (i.e., is high in uncertainty and low in surprise) or when it challenges their highly confident ones (i.e., is low in uncertainty and high in surprise). An independent research group (Cheung et al., 2019) replicated these results with musical chord sequences, and identified their fMRI correlates in the STG, amygdala, and hippocampus but not the VS, raising new questions about the neural mechanisms of musical pleasure that the present study seeks to address. Here, we assessed concurrent liking ratings and hemodynamic fMRI signals as 24 participants listened to 50 naturalistic, real-world musical excerpts that varied across wide spectra of computationally modeled uncertainty and surprise. As in previous studies, liking ratings exhibited an interaction between uncertainty and surprise, with the strongest preferences for high uncertainty/low surprise and low uncertainty/high surprise. FMRI results also replicated previous findings, with music liking effects in the R STG and VS. Furthermore, we identify interactions between uncertainty and surprise on the one hand, and liking and surprise on the other, in VS activity. Altogether, these results provide important support for the hypothesized role of the VS in deriving pleasure from learning about musical structure.

Statistical Approaches to Identify Pairwise and High-Order Brain Functional Connectivity Signatures on a Single-Subject Basis.

Keeping up with the shift towards personalized neuroscience essentially requires the derivation of meaningful insights from individual brain signal recordings by analyzing the descriptive indexes of physio-pathological states through statistical methods that prioritize subject-specific differences under varying experimental conditions. Within this framework, the current study presents a methodology for assessing the value of the single-subject fingerprints of brain functional connectivity, assessed both by standard pairwise and novel high-order measures. Functional connectivity networks, which investigate the inter-relationships between pairs of brain regions, have long been a valuable tool for modeling the brain as a complex system. However, their usefulness is limited by their inability to detect high-order dependencies beyond pairwise correlations. In this study, by leveraging multivariate information theory, we confirm recent evidence suggesting that the brain contains a plethora of high-order, synergistic subsystems that would go unnoticed using a pairwise graph structure. The significance and variations across different conditions of functional pairwise and high-order interactions (HOIs) between groups of brain signals are statistically verified on an individual level through the utilization of surrogate and bootstrap data analyses. The approach is illustrated on the single-subject recordings of resting-state functional magnetic resonance imaging (rest-fMRI) signals acquired using a pediatric patient with hepatic encephalopathy associated with a portosystemic shunt and undergoing liver vascular shunt correction. Our results show that (i) the proposed single-subject analysis may have remarkable clinical relevance for subject-specific investigations and treatment planning, and (ii) the possibility of investigating brain connectivity and its post-treatment functional developments at a high-order level may be essential to fully capture the complexity and modalities of the recovery.

Multifractal long-range dependence pattern of functional magnetic resonance imaging in the human brain at rest.

Long-range dependence is a prevalent phenomenon in various biological systems that characterizes the long-memory effect of temporal fluctuations. While recent research suggests that functional magnetic resonance imaging signal has fractal property, it remains unknown about the multifractal long-range dependence pattern of resting-state functional magnetic resonance imaging signals. The current study adopted the multifractal detrended fluctuation analysis on highly sampled resting-state functional magnetic resonance imaging scans to investigate long-range dependence profile associated with the whole-brain voxels as specific functional networks. Our findings revealed the long-range dependence's multifractal properties. Moreover, long-term persistent fluctuations are found for all stations with stronger persistency in whole-brain regions. Subsets with large fluctuations contribute more to the multifractal spectrum in the whole brain. Additionally, we found that the preprocessing with band-pass filtering provided significantly higher reliability for estimating long-range dependence. Our validation analysis confirmed that the optimal pipeline of long-range dependence analysis should include band-pass filtering and removal of daily temporal dependence. Furthermore, multifractal long-range dependence characteristics in healthy control and schizophrenia are different significantly. This work has provided an analytical pipeline for the multifractal long-range dependence in the resting-state functional magnetic resonance imaging signal. The findings suggest differential long-memory effects in the intrinsic functional networks, which may offer a neural marker finding for understanding brain function and pathology.