Effects of spatial smoothing on inter-subject correlation based analysis of FMRI

Recommended Resting-State fMRI Acquisition and Preprocessing Steps for Preoperative Mapping of Language and Motor and Visual Areas in Adult and Pediatric Patients with Brain Tumors and Epilepsy

Resting-state (rs) fMRI has been shown to be useful for preoperative mapping of functional areas in patients with brain tumors and epilepsy. However, its lack of standardization limits its widespread use and hinders multicenter collaboration. The American Society of Functional Neuroradiology, American Society of Pediatric Neuroradiology, and the American Society of Neuroradiology Functional and Diffusion MR Imaging Study Group recommend specific rs-fMRI acquisition approaches and preprocessing steps that will further support rs-fMRI for future clinical use. A task force with expertise in fMRI from multiple institutions provided recommendations on the rs-fMRI steps needed for mapping of language, motor, and visual areas in adult and pediatric patients with brain tumor and epilepsy. These were based on an extensive literature review and expert consensus.Following rs-fMRI acquisition parameters are recommended: minimum 6-minute acquisition time; scan with eyes open with fixation; obtain rs-fMRI before both task-based fMRI and contrast administration; temporal resolution of ≤2 seconds; scanner field strength of 3T or higher. The following rs-fMRI preprocessing steps and parameters are recommended: motion correction (seed-based correlation analysis [SBC], independent component analysis [ICA]); despiking (SBC); volume censoring (SBC, ICA); nuisance regression of CSF and white matter signals (SBC); head motion regression (SBC, ICA); bandpass filtering (SBC, ICA); and spatial smoothing with a kernel size that is twice the effective voxel size (SBC, ICA).The consensus recommendations put forth for rs-fMRI acquisition and preprocessing steps will aid in standardization of practice and guide rs-fMRI program development across institutions. Standardized rs-fMRI protocols and processing pipelines are essential for multicenter trials and to implement rs-fMRI as part of standard clinical practice.

Spatial Smoothing Effect on Group-Level Functional Connectivity during Resting and Task-Based fMRI

Spatial smoothing is a preprocessing step applied to neuroimaging data to enhance data quality by reducing noise and artifacts. However, selecting an appropriate smoothing kernel size can be challenging as it can lead to undesired alterations in final images and functional connectivity networks. However, there is no sufficient information about the effects of the Gaussian kernel size on group-level results for different cases yet. This study investigates the influence of kernel size on functional connectivity networks and network parameters in whole-brain rs-fMRI and tb-fMRI analyses of healthy adults. The analysis includes {0, 2, 4, 6, 8, 10} mm kernels, commonly used in practical analyses, covering all major brain networks. Graph theoretical measures such as betweenness centrality, global/local efficiency, clustering coefficient, and average path length are examined for each kernel. Additionally, principal component analysis (PCA) and independent component analysis (ICA) parameters, namely kurtosis and skewness, are evaluated for the functional images. The findings demonstrate that kernel size directly affects node connections, resulting in modifications to functional network structures and PCA/ICA parameters. However, network metrics exhibit greater resilience to these changes.

Functional MRI responses to naturalistic stimuli are increasingly typical across early childhood

While findings show that throughout development, there are child- and age-specific patterns of brain functioning, there is also evidence for significantly greater inter-individual response variability in young children relative to adults. It is currently unclear whether this increase in functional "typicality" (i.e., inter-individual similarity) is a developmental process that occurs across early childhood, and what changes in BOLD response may be driving changes in typicality. We collected fMRI data from 81 typically developing 4-8-year-old children during passive viewing of age-appropriate television clips and asked whether there is increasing typicality of brain response across this age range. We found that the "increasing typicality" hypothesis was supported across many regions engaged by passive viewing. Post hoc analyses showed that in a priori ROIs related to language and face processing, the strength of the group-average shared component of activity increased with age, with no concomitant decline in residual signal or change in spatial extent or variability. Together, this suggests that increasing inter-individual similarity of functional responses to audiovisual stimuli is an important feature of early childhood functional brain development.

Thalamo-cortical inter-subject functional correlation during movie watching across the adult lifespan

An increasing number of studies have shown that the functional interactions between the thalamus and cerebral cortices play an important role in cognitive function and are influenced by age. Previous studies have revealed age-related changes in the thalamo-cortical system within individuals, while neglecting differences between individuals. Here, we characterized inter-subject functional correlation (ISFC) between the thalamus and several cortical brain networks in 500 healthy participants aged 18–87 years old from the Cambridge Centre for Aging and Neuroscience (Cam-CAN) cohort using movie-watching state fMRI data. General linear models (GLM) were performed to assess age-related changes in ISFC of thalamo-cortical networks and the relationship between ISFC and fluid intelligence. We found significant age-related decreases in ISFC between the posterior thalamus (e.g., ventral posterior nucleus and pulvinar) and the attentional network, sensorimotor network, and visual network (FDR correction with p < 0.05). Meanwhile, the ISFC between the thalamus (mainly the mediodorsal nucleus and ventral thalamic nuclei) and higher-order cortical networks, including the default mode network, salience network and control network, showed complex changes with age. Furthermore, the altered ISFC of thalamo-cortical networks was positively correlated with decreased fluid intelligence (FDR correction with p < 0.05). Overall, our results provide further evidence that alterations in the functional integrity of the thalamo-cortical system might play an important role in cognitive decline during aging.

Evaluating denoising strategies in resting-state functional magnetic resonance in traumatic brain injury (EpiBioS4Rx)

Resting-state functional MRI is increasingly used in the clinical setting and is now included in some diagnostic guidelines for severe brain injury patients. However, to ensure high-quality data, one should mitigate fMRI-related noise typical of this population. Therefore, we aimed to evaluate the ability of different preprocessing strategies to mitigate noise-related signal (i.e., in-scanner movement and physiological noise) in functional connectivity (FC) of traumatic brain injury (TBI) patients. We applied nine commonly used denoising strategies, combined into 17 pipelines, to 88 TBI patients from the Epilepsy Bioinformatics Study for Anti-epileptogenic Therapy clinical trial. Pipelines were evaluated by three quality control (QC) metrics across three exclusion regimes based on the participant's head movement profile. While no pipeline eliminated noise effects on FC, some pipelines exhibited relatively high effectiveness depending on the exclusion regime. Once high-motion participants were excluded, the choice of denoising pipeline becomes secondary - although this strategy leads to substantial data loss. Pipelines combining spike regression with physiological regressors were the best performers, whereas pipelines that used automated data-driven methods performed comparatively worse. In this study, we report the first large-scale evaluation of denoising pipelines aimed at reducing noise-related FC in a clinical population known to be highly susceptible to in-scanner motion and significant anatomical abnormalities. If resting-state functional magnetic resonance is to be a successful clinical technique, it is crucial that procedures mitigating the effect of noise be systematically evaluated in the most challenging populations, such as TBI datasets.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

Inattentive and hyperactive traits differentially associate with inter-individual functional synchrony during video viewing in young children without ADHD

Inattention and hyperactivity present on a spectrum and may influence the way children perceive and interact with the world. We investigated whether normative variation in inattentive and hyperactive traits was associated with differences in brain function, while children watched clips from an age-appropriate television program. Functional magnetic resonance imaging (fMRI) data and parent reports of inattention and hyperactivity traits were collected from 81 children 4–7 years of age with no parent-reported diagnoses. Data were analyzed using intersubject correlations (ISCs) in mixed effects models to determine if inattentive and hyperactive traits were associated with idiosyncrasy of fMRI response to the video. We hypothesized that pairs of children with higher average inattention and hyperactivity scores would show less interindividual brain synchrony to one another than pairs with lower average scores on these traits. Video watching engaged widespread visual, auditory, default mode and dorsal prefrontal regions. Inattention and hyperactivity were separably associated with ISC in many of these regions. Our findings suggest that the spectrum of inattention and hyperactivity traits in children without ADHD are differentially associated with neural processing of naturalistic video stimuli, which may have implications for understanding how children with different levels of these traits process audiovisual information in unconstrained conditions.

Effects of Parent Emotion Socialization on the Neurobiology Underlying Adolescent Emotion Processing: A Multimethod fMRI Study

Parents' emotion socialization (ES) practices impact socioemotional development throughout adolescence. Little is known, however, regarding the neurobiology underlying these effects. This study used functional magnetic resonance imaging (fMRI) to examine how parent ES practices relate to adolescent brain function during emotion processing. Thirty-three adolescents (ages 14-16) reported on ES practices of a focal parent (primarily mothers) using the Emotions as a Child (EAC) Scale. Adolescents also completed a conflict discussion task with this parent, and parents' statements were coded for emotional valence. Adolescents performed two fMRI tasks: a standard emotion processing (EP) task (n = 32) and the Testing Emotional Attunement and Mutuality (TEAM) task (n = 27). The EP task consisted of viewing emotional pictures and either reacting naturally or using cognitive reappraisal to regulate emotional responses. The TEAM task was performed with the parent and included trials during which adolescents were shown that their parent made an error, costing the dyad $5. Parent negative verbalizations during the conflict discussion were associated with greater activity in the thalamus during the emotion reactivity condition of the EP task and in the thalamus, superior medial and superior frontal gyri, anterior insula, and dorsolateral prefrontal cortex during the costly error condition of the TEAM task. Unsupportive ES was associated with greater activity in the supplementary motor area and less activity in the paracentral gyrus and amygdala during the costly error condition of the TEAM task. This study supports the premise that ES influences adolescents' emotion-related neural processing, particularly when using ecologically valid tasks in social contexts.

Connectome spatial smoothing (CSS): Concepts, methods, and evaluation

Structural connectomes are increasingly mapped at high spatial resolutions comprising many hundreds-if not thousands-of network nodes. However, high-resolution connectomes are particularly susceptible to image registration misalignment, tractography artifacts, and noise, all of which can lead to reductions in connectome accuracy and test-retest reliability. We investigate a network analogue of image smoothing to address these key challenges. Connectome Spatial Smoothing (CSS) involves jointly applying a carefully chosen smoothing kernel to the two endpoints of each tractography streamline, yielding a spatially smoothed connectivity matrix. We develop computationally efficient methods to perform CSS using a matrix congruence transformation and evaluate a range of different smoothing kernel choices on CSS performance. We find that smoothing substantially improves the identifiability, sensitivity, and test-retest reliability of high-resolution connectivity maps, though at a cost of increasing storage burden. For atlas-based connectomes (i.e. low-resolution connectivity maps), we show that CSS marginally improves the statistical power to detect associations between connectivity and cognitive performance, particularly for connectomes mapped using probabilistic tractography. CSS was also found to enable more reliable statistical inference compared to connectomes without any smoothing. We provide recommendations for optimal smoothing kernel parameters for connectomes mapped using both deterministic and probabilistic tractography. We conclude that spatial smoothing is particularly important for the reliability of high-resolution connectomes, but can also provide benefits at lower parcellation resolutions. We hope that our work enables computationally efficient integration of spatial smoothing into established structural connectome mapping pipelines.

Topography of hippocampal connectivity with sensorimotor cortex revealed by optimizing smoothing kernel and voxel size

Studies of the hippocampus use smaller voxel sizes and smoothing kernels than cortical activation studies, typically using a multivoxel seed with specified radius for connectivity analysis. This study identified optimal processing parameters for evaluating hippocampal connectivity with sensorimotor cortex (SMC), comparing effectiveness by varying parameters during both activation and connectivity analysis. Using both 3mm and 4mm isovoxels, smoothing kernels of 0-10mm were evaluated on the amplitude and extent of motor activation and hippocampal connectivity with SMC. Psychophysiological interactions (PPI) identified hippocampal connectivity with SMC during volitional movements, and connectivity effects from multivoxel seeds were compared with alternate methods; a structural seed represented the mean connectivity map from all voxels within a region, whereas a functional seed represented the regional voxel with maximal SMC connectivity. With few exceptions, the same parameters were optimal for activation and connectivity. Larger isovoxels showed larger activation volumes in both SMC and the hippocampus; connectivity volumes from structural seeds were also larger, except from the posterior hippocampus. Regardless of voxel size, the 10mm smoothing kernel generated larger activation and connectivity volumes from structural seeds, as well as larger beta estimates at connectivity maxima; structural seeds also produced larger connectivity volumes than multivoxel seeds. Functional seeds showed lesser effects from voxel size and smoothing kernels. Optimal parameters revealed topography in structural seed connectivity along both the longitudinal axis and mediolateral axis of the hippocampus. These results indicate larger voxels and smoothing kernels can improve sensitivity for detecting both cortical activation and hippocampal connectivity.

The "Narratives" fMRI dataset for evaluating models of naturalistic language comprehension

The "Narratives" collection aggregates a variety of functional MRI datasets collected while human subjects listened to naturalistic spoken stories. The current release includes 345 subjects, 891 functional scans, and 27 diverse stories of varying duration totaling ~4.6 hours of unique stimuli (~43,000 words). This data collection is well-suited for naturalistic neuroimaging analysis, and is intended to serve as a benchmark for models of language and narrative comprehension. We provide standardized MRI data accompanied by rich metadata, preprocessed versions of the data ready for immediate use, and the spoken story stimuli with time-stamped phoneme- and word-level transcripts. All code and data are publicly available with full provenance in keeping with current best practices in transparent and reproducible neuroimaging.

rsfMRI based evidence for functional connectivity alterations in adults with developmental stuttering

Persistent developmental stuttering (PDS) is defined as a speech disorder mainly characterized by intermittent involuntary disruption in normal fluency, time patterning, and rhythm of speech. Although extensive functional neuroimaging studies have explored brain activation alterations in stuttering, the main affected brain regions/networks in PDS still remain unclear. Here, using functional magnetic resonance imaging (fMRI), we investigated resting-state whole-brain functional connectivity of 15 adults who stutter (PDS group) and 15 age-matched control individuals to reveal the connectivity abnormalities associated with stuttering. We were also interested in exploring how the severity of stuttering varies across individuals to understand the compensatory mechanism of connectivity pattern in patients showing less symptoms. Our results revealed decreased connectivity of left frontal pole and left middle frontal gyrus (MidFG) with right precentral/postcentral gyrus in stuttering individuals compared with control participants, while less symptomatic PDS individuals showed greater functional connectivity between left MidFG and left caudate. Additionally, our finding indicated reduced connectivity in the PDS group between the left superior temporal gyrus (STG) and several brain regions including the right limbic lobe, right fusiform, and right cerebellum, as well as the left middle temporal gyrus (MTG). We also observed that PDS individuals with less severe symptoms had stronger connectivity between right MTG and several left hemispheric regions including inferior frontal gyrus (IFG) and STG. The connectivity between right fronto-orbital and right MTG was also negatively correlated with stuttering severity. These findings may suggest the involvement of right MTG and left MidFG in successful compensatory mechanisms in more fluent stutterers.

Cortical and subcortical contributions to interference resolution and inhibition - An fMRI ALE meta-analysis

Interacting with our environment requires the selection of appropriate responses and the inhibition of others. Such effortful inhibition is achieved by a number of interference resolution and global inhibition processes. This meta-analysis including 57 studies and 73 contrasts revisits the overlap and differences in brain areas supporting interference resolution and global inhibition in cortical and subcortical brain areas. Activation likelihood estimation was used to discern the brain regions subserving each type of cognitive control. Individual contrast analysis revealed a common activation of the bilateral insula and supplementary motor areas. Subtraction analyses demonstrated the voxel-wise differences in recruitment in a number of areas including the precuneus in the interference tasks and the frontal pole and dorsal striatum in the inhibition tasks. Our results display a surprising lack of subcortical involvement within these types of cognitive control, a finding that is likely to reflect a systematic gap in the field of functional neuroimaging.

Functional Connectivity of the Human Paraventricular Thalamic Nucleus: Insights From High Field Functional MRI

The paraventricular thalamic nucleus (PVT) is a small but highly connected nucleus of the dorsal midline thalamus. The PVT has garnered recent attention as a context-sensitive node within the thalamocortical arousal system that modulates state-dependent motivated behaviors. Once considered related to generalized arousal responses with non-specific impacts on behavior, accumulating evidence bolsters the contemporary view that discrete midline thalamic subnuclei belong to specialized corticolimbic and corticostriatal circuits related to attention, emotions, and cognition. However, the functional connectivity patterns of the human PVT have yet to be mapped. Here, we combined high-quality, high-resolution 7T and 3T resting state MRI data from 121 young adult participants from the Human Connectome Project (HCP) and thalamic subnuclei atlas masks to investigate resting state functional connectivity of the human PVT. The 7T results demonstrated extensive positive functional connectivity with the brainstem, midbrain, ventral and dorsal medial prefrontal cortex (mPFC), anterior and posterior cingulate, ventral striatum, hippocampus, and amygdala. These connections persist upon controlling for functional connectivity of the rest of the thalamus. Whole-brain contrasts provided further evidence that, compared to three nearby midline thalamic subnuclei, functional connectivity of the PVT is strong with the hippocampus, amygdala, ventral and dorsal mPFC, and middle temporal gyrus. These findings suggest that, even during rest, the human PVT is functionally coupled with many regions known to be structurally connected to rodent and non-human primate PVT. Further, cosine similarity analysis results suggested the PVT is integrated into the default mode network (DMN), an intrinsic connectivity network associated with episodic memory and self-referential thought. The current work provides a much-needed foundation for ongoing and future work examining the functional roles of the PVT in humans.

Effects of different smoothing on global and regional resting functional connectivity

PURPOSE

Spatial smoothing is an essential pre-processing step in the process of analysing functional magnetic resonance imaging (fMRI) data, both during an experimental task or during resting-state fMRI (rsfMRI). The main benefit of this spatial smoothing step is to artificially increase the signal-to-noise ratio of the fMRI signal. Previous fMRI studies have investigated the impact of spatial smoothing on task fMRI data, while rsfMRI studies usually apply the same analytical process used for the task data. However, this study investigates changes in different rsfMRI analyses, such as ROI-to-ROI, seed-to-voxels and ICA analyses.

METHODS

Nineteen healthy volunteers were scanned using rsfMRI with three applied smoothing kernels: 0 mm, 4 mm and 8 mm. Appropriate statistical comparisons were made.

RESULTS

The findings showed that spatial smoothing has a greater effect on rsfMRI data when analysed using seed-to-voxel-based analysis. The effect was less pronounced when analysing data using ROI-ROI or ICA analyses. The results demonstrated that even when analysing the data without the application of spatial smoothing, the results were significant compared with data analysed using a typical smoothing kernel. However, data analysed with lower-smoothing kernels produced greater negative correlations, particularly with the ICA analysis.

CONCLUSION

The results suggest that a medium smoothing kernel (around 4 mm) may be preferable, as it is comparable with the 8 mm kernel in all of the analyses performed. It is also recommended that the researchers consider analysing the data using two different smoothing kernels, as this will help to confirm the significance of the results and avoid overestimating the findings.

Effects of spatial smoothing on group-level differences in functional brain networks

Brain connectivity with functional magnetic resonance imaging (fMRI) is a popular approach for detecting differences between healthy and clinical populations. Before creating a functional brain network, the fMRI time series must undergo several preprocessing steps to control for artifacts and to improve data quality. However, preprocessing may affect the results in an undesirable way. Spatial smoothing, for example, is known to alter functional network structure. Yet, its effects on group-level network differences remain unknown. Here, we investigate the effects of spatial smoothing on the difference between patients and controls for two clinical conditions: autism spectrum disorder and bipolar disorder, considering fMRI data smoothed with Gaussian kernels (0–32 mm). We find that smoothing affects network differences between groups. For weighted networks, incrementing the smoothing kernel makes networks more different. For thresholded networks, larger smoothing kernels lead to more similar networks, although this depends on the network density. Smoothing also alters the effect sizes of the individual link differences. This is independent of the region of interest (ROI) size, but varies with link length. The effects of spatial smoothing are diverse, nontrivial, and difficult to predict. This has important consequences: The choice of smoothing kernel affects the observed network differences.

Spatial smoothing is a preprocessing tool commonly applied to reduce the amount of noise in functional magnetic resonance imaging (fMRI) data. However, smoothing is known to affect the outcomes of functional brain network analysis at the level of individual subjects in undesired ways. Here, we investigate how spatial smoothing affects the observed differences in brain network structure between subject groups. Using fMRI data from two clinical populations and healthy controls, we show that the between-group differences in network structure depend on the amount of spatial smoothing applied during preprocessing in a nontrivial way. The optimal level of spatial smoothing is difficult to define and probably depends on a set of analysis parameters. Therefore, we recommend applying spatial smoothing only after careful consideration.

First-Person Virtual Embodiment Modulates the Cortical Network that Encodes the Bodily Self and Its Surrounding Space during the Experience of Domestic Violence

Social aggression, such as domestic violence, has been associated with a reduced ability to take on others’ perspectives. In this naturalistic imaging study, we investigated whether training human participants to take on a first-person embodied perspective during the experience of domestic violence enhances the identification with the victim and elicits brain activity associated with the monitoring of the body and surrounding space and the experience of threat. We combined fMRI measurements with preceding virtual reality exposure from either first-person perspective (1PP) or third-person perspective (3PP) to manipulate whether the domestic abuse stimulus was perceived as directed to oneself or another. We found that 1PP exposure increased body ownership and identification with the virtual victim. Furthermore, when the stimulus was perceived as directed toward oneself, the brain network that encodes the bodily self and its surrounding space was more strongly synchronized across participants and connectivity increased from premotor cortex (PM) and intraparietal sulcus towards superior parietal lobe. Additionally, when the stimulus came near the body, brain activity in the amygdala (AMG) strongly synchronized across participants. Exposure to 3PP reduced synchronization of brain activity in the personal space network, increased modulation of visual areas and strengthened functional connectivity between PM, supramarginal gyrus and primary visual cortex. In conclusion, our results suggest that 1PP embodiment training enhances experience from the viewpoint of the virtual victim, which is accompanied by synchronization in the fronto-parietal network to predict actions toward the body and in the AMG to signal the proximity of the stimulus.

Dissociable neural systems for unconditioned acute and sustained fear

Fear protects organisms by increasing vigilance and preparedness, and by coordinating survival responses during life-threatening encounters. The fear circuit must thus operate on multiple timescales ranging from preparatory sustained alertness to acute fight-or-flight responses. Here we studied the brain basis of sustained and acute fear using naturalistic functional magnetic resonance imaging (fMRI) enabling analysis of different time-scales of fear responses. Subjects (N ​= ​37) watched feature-length horror movies while their hemodynamic brain activity was measured with fMRI. Time-variable intersubject correlation (ISC) was used to quantify the reliability of brain activity across participants, and seed-based phase synchronization was used for characterizing dynamic connectivity. Subjective ratings of fear were used to assess how synchronization and functional connectivity varied with emotional intensity. These data suggest that acute and sustained fear are supported by distinct neural pathways, with sustained fear amplifying mainly sensory responses, and acute fear increasing activity in brainstem, thalamus, amygdala and cingulate cortices. Sustained fear increased ISC in regions associated with acute fear, and also amplified functional connectivity within this network. The results were replicated in an independent experiment with a different subject sample and stimulus movie. The functional interplay between cortical networks involved in sustained anticipation of, and acute response to, threat involves a complex and dynamic interaction that depends on the proximity of threat, and the need to employ threat appraisals and vigilance for decision making and response selection.

Connectivity of sleep- and wake-promoting regions of the human hypothalamus observed during resting wakefulness

The hypothalamus is a central hub for regulating sleep-wake patterns, the circuitry of which has been investigated extensively in experimental animals. This work has identified a wake-promoting region in the posterior hypothalamus, with connections to other wake-promoting regions, and a sleep-promoting region in the anterior hypothalamus, with inhibitory projections to the posterior hypothalamus. It is unclear whether a similar organization exists in humans. Here, we use anatomical landmarks to identify homologous sleep- and wake-promoting regions of the human hypothalamus and investigate their functional relationships using resting-state functional connectivity magnetic resonance imaging in healthy awake participants. First, we identify a negative correlation (anticorrelation) between the anterior and posterior hypothalamus, two regions with opposing roles in sleep-wake regulation. Next, we show that hypothalamic connectivity predicts a pattern of regional sleep-wake changes previously observed in humans. Specifically, regions that are more positively correlated with the posterior hypothalamus and more negatively correlated with the anterior hypothalamus correspond to regions with the greatest change in cerebral blood flow between sleep-wake states. Taken together, these findings provide preliminary evidence relating a hypothalamic circuit investigated in animals to sleep-wake neuroimaging results in humans, with implications for our understanding of human sleep-wake regulation and the functional significance of anticorrelations.

Emotions amplify speaker-listener neural alignment

Individuals often align their emotional states during conversation. Here, we reveal how such emotional alignment is reflected in synchronization of brain activity across speakers and listeners. Two "speaker" subjects told emotional and neutral autobiographical stories while their hemodynamic brain activity was measured with functional magnetic resonance imaging (fMRI). The stories were recorded and played back to 16 "listener" subjects during fMRI. After scanning, both speakers and listeners rated the moment-to-moment valence and arousal of the stories. Time-varying similarity of the blood-oxygenation-level-dependent (BOLD) time series was quantified by intersubject phase synchronization (ISPS) between speaker-listener pairs. Telling and listening to the stories elicited similar emotions across speaker-listener pairs. Arousal was associated with increased speaker-listener neural synchronization in brain regions supporting attentional, auditory, somatosensory, and motor processing. Valence was associated with increased speaker-listener neural synchronization in brain regions involved in emotional processing, including amygdala, hippocampus, and temporal pole. Speaker-listener synchronization of subjective feelings of arousal was associated with increased neural synchronization in somatosensory and subcortical brain regions; synchronization of valence was associated with neural synchronization in parietal cortices and midline structures. We propose that emotion-dependent speaker-listener neural synchronization is associated with emotional contagion, thereby implying that listeners reproduce some aspects of the speaker's emotional state at the neural level.

Measuring shared responses across subjects using intersubject correlation

Our capacity to jointly represent information about the world underpins our social experience. By leveraging one individual's brain activity to model another's, we can measure shared information across brains-even in dynamic, naturalistic scenarios where an explicit response model may be unobtainable. Introducing experimental manipulations allows us to measure, for example, shared responses between speakers and listeners or between perception and recall. In this tutorial, we develop the logic of intersubject correlation (ISC) analysis and discuss the family of neuroscientific questions that stem from this approach. We also extend this logic to spatially distributed response patterns and functional network estimation. We provide a thorough and accessible treatment of methodological considerations specific to ISC analysis and outline best practices.

Altered Local and Large-Scale Dynamic Functional Connectivity Variability in Posttraumatic Stress Disorder: A Resting-State fMRI Study

Posttraumatic stress disorder (PTSD) is a psychiatric condition that can emerge after exposure to an exceedingly traumatic event. Previous neuroimaging studies have indicated that PTSD is characterized by aberrant resting-state functional connectivity (FC). However, few existing studies on PTSD have examined dynamic changes in resting-state FC related to network formation, interaction, and dissolution over time. In this study, we compared the dynamic resting-state local and large-scale FC between PTSD patients (n = 22) and healthy controls (HC; n = 22; conducted as standard deviation in resting-state local and large-scale FC over a series of sliding windows). Local dynamic FC was examined by calculating the dynamic regional homogeneity (dReHo), and large-scale dynamic FC (dFC) was investigated between regions with significant dReHo group differences. For the PTSD patients, we also investigated the relationship between symptom severity and dFC/dReHo. Our results showed that PTSD patients were characterized by I) increased dynamic (more variable) dReHo in left precuneus (PCu); II) increased dynamic (more variable) dFC between the left PCu and left insula; and III) decreased dFC between left PCu and left inferior parietal lobe (IPL), and decreased dFC between left PCu and right PCu. However, there is no significant correlation between the clinical indicators and dReHo/dFC after the family-wise-error (FWE) correction. These findings provided the initial evidence that PTSD is characterized by aberrant patterns of fluctuating communication within brain system such as the default mode network (DMN) and among different brain systems such as the salience network and the DMN.

Developmental changes in functional brain networks from birth through adolescence

Investigation of the brain connectome using functional magnetic resonance imaging (fMRI) and measures derived from graph theory analysis has emerged as a new approach to study brain development, cognitive function, and neurophysiological disorders. Here we use graph theory analysis to examine the influence of age, sex, and neurocognitive measures on developmental changes to the global and regional topology of functional brain networks derived from fMRI data recorded in 189 healthy subjects from the age of 0-18 years during rest. We observed that Global Efficiency and Rich-Club coefficient increased with age and Local Efficiency and Small-Worldness decreased with age, while Modularity at the global level showed an inverted U-shaped trajectory during development. Marginally significant differences were observed in Local Efficiency, Small-Worldness, and Modularity at a global level between boys and girls throughout development. We also examine the effects of neurocognitive measures in boys and girls globally and locally. Our results provide new insight to understand brain maturation of functional brain connectome and its relation to cognitive development from birth through adolescence.

Effective connectivity of the anterior hippocampus predicts recollection confidence during natural memory retrieval

Human interactions with the world are influenced by memories of recent events. This effect, often triggered by perceptual cues, occurs naturally and without conscious effort. However, the neuroscience of involuntary memory in a dynamic milieu has received much less attention than the mechanisms of voluntary retrieval with deliberate purpose. Here, we investigate the neural processes driven by naturalistic cues that relate to, and presumably trigger the retrieval of recent experiences. Viewing the continuation of recently viewed clips evokes greater bilateral activation in anterior hippocampus, precuneus and angular gyrus than naïve clips. While these regions manifest reciprocal connectivity, continued viewing specifically modulates the effective connectivity from the anterior hippocampus to the precuneus. The strength of this modulation predicts participants' confidence in later voluntary recall of news details. Our study reveals network mechanisms of dynamic, involuntary memory retrieval and its relevance to metacognition in a rich context resembling everyday life.

Impact of global signal regression on characterizing dynamic functional connectivity and brain states

Recently, resting-state functional magnetic resonance imaging (fMRI) studies have been extended to explore fluctuations in correlations over shorter timescales, referred to as dynamic functional connectivity (dFC). However, the impact of global signal regression (GSR) on dFC is not well established, despite the intensive investigations of the influence of GSR on static functional connectivity (sFC). This study aimed to examine the effect of GSR on the performance of the sliding-window correlation, a commonly used method for capturing functional connectivity (FC) dynamics based on resting-state fMRI and simultaneous electroencephalograph (EEG)-fMRI data. The results revealed that the impact of GSR on dFC was spatially heterogeneous, with some susceptible regions including the occipital cortex, sensorimotor area, precuneus, posterior insula and superior temporal gyrus, and that the impact was temporally modulated by the mean global signal (GS) magnitude across windows. Furthermore, GSR substantially changed the connectivity structures of the FC states responding to a high GS magnitude, as well as their temporal features, and even led to the emergence of new FC states. Conversely, those FC states marked by obvious anti-correlation structures associated with the default model network (DMN) were largely unaffected by GSR. Finally, we reported an association between the fluctuations in the windowed magnitude of GS and the time-varying EEG power within subjects, which implied changes in mental states underlying GS dynamics. Overall, this study suggested a potential neuropsychological basis, in addition to nuisance sources, for GS dynamics and highlighted the need for caution in applying GSR to sliding-window correlation analyses. At a minimum, the mental fluctuations of an individual subject, possibly related to ongoing vigilance, should be evaluated during the entire scan when the dynamics of FC is estimated.

Inter-subject synchrony as an index of functional specialization in early childhood

Early childhood is a time of significant change within multiple cognitive domains including social cognition, memory, executive function, and language; however, the corresponding neural changes remain poorly understood. This is likely due to the difficulty in acquiring artifact-free functional MRI data during complex task-based or unconstrained resting-state experiments in young children. In addition, task-based and resting state experiments may not capture dynamic real-world processing. Here we overcome both of these challenges through use of naturalistic viewing (i.e., passively watching a movie in the scanner) combined with inter-subject neural synchrony to examine functional specialization within 4- and 6-year old children. Using a novel and stringent crossed random effect statistical analysis, we find that children show more variable patterns of activation compared to adults, particularly within regions of the default mode network (DMN). In addition, we found partial evidence that child-to-adult synchrony increased as a function of age within a DMN region: the temporoparietal junction. Our results suggest age-related differences in functional brain organization within a cross-sectional sample during an ecologically valid context and demonstrate that neural synchrony during naturalistic viewing fMRI can be used to examine functional specialization during early childhood - a time when neural and cognitive systems are in flux.

Effect of Spatial Smoothing on Task fMRI ICA and Functional Connectivity

Spatial smoothing is a widely used preprocessing step in functional magnetic resonance imaging (fMRI) data analysis. In this work, we report on the spatial smoothing effect on task-evoked fMRI brain functional mapping and functional connectivity. Initially, we decomposed the task fMRI data into a collection of components or networks by independent component analysis (ICA). The designed task paradigm helps identify task-modulated ICA components (highly correlated with the task stimuli). For the ICA-extracted primary task component, we then measured the task activation volume at the task response foci. We used the task timecourse (designed) as a reference to order the ICA components according to the task correlations of the ICA timecourses. With the re-ordered ICA components, we calculated the inter-component function connectivity (FC) matrix (correlations among the ICA timecourses). By repeating the spatial smoothing of fMRI data with a Gaussian smoothing kernel with a full width at half maximum (FWHM) of {1, 3, 6, 9, 12, 15, 20, 25, 30, 35} mm, we measured the spatial smoothing effects. Our results show spatial smoothing reveals the following effects: (1) It decreases the task extraction performance of single-subject ICA more than that of multi-subject ICA; (2) It increases the task volume of multi-subject ICA more than that of single-subject ICA; (3) It strengthens the functional connectivity of single-subject ICA more than that of multi-subject ICA; and (4) It impacts the positive-negative imbalance of single-subject ICA more than that of multi-subject ICA. Our experimental results suggest a 2~3 voxel FWHM spatial smoothing for single-subject ICA in achieving an optimal balance of functional connectivity, and a wide range (2~5 voxels) of FWHM for multi-subject ICA.

Listening in Naturalistic Scenes: What Can Functional Near-Infrared Spectroscopy and Intersubject Correlation Analysis Tell Us About the Underlying Brain Activity?

Listening to speech in the noisy conditions of everyday life can be effortful, reflecting the increased cognitive workload involved in extracting meaning from a degraded acoustic signal. Studying the underlying neural processes has the potential to provide mechanistic insight into why listening is effortful under certain conditions. In a move toward studying listening effort under ecologically relevant conditions, we used the silent and flexible neuroimaging technique functional near-infrared spectroscopy (fNIRS) to examine brain activity during attentive listening to speech in naturalistic scenes. Thirty normally hearing participants listened to a series of narratives continuously varying in acoustic difficulty while undergoing fNIRS imaging. Participants then listened to another set of closely matched narratives and rated perceived effort and intelligibility for each scene. As expected, self-reported effort generally increased with worsening signal-to-noise ratio. After controlling for better-ear signal-to-noise ratio, perceived effort was greater in scenes that contained competing speech than in those that did not, potentially reflecting an additional cognitive cost of overcoming informational masking. We analyzed the fNIRS data using intersubject correlation, a data-driven approach suitable for analyzing data collected under naturalistic conditions. Significant intersubject correlation was seen in the bilateral auditory cortices and in a range of channels across the prefrontal cortex. The involvement of prefrontal regions is consistent with the notion that higher order cognitive processes are engaged during attentive listening to speech in complex real-world conditions. However, further research is needed to elucidate the relationship between perceived listening effort and activity in these extended cortical networks.

Effects of spatial smoothing on functional brain networks

Graph-theoretical methods have rapidly become a standard tool in studies of the structure and function of the human brain. Whereas the structural connectome can be fairly straightforwardly mapped onto a complex network, there are more degrees of freedom in constructing networks that represent functional connections between brain areas. For functional magnetic resonance imaging (fMRI) data, such networks are typically built by aggregating the blood-oxygen-level dependent signal time series of voxels into larger entities (such as Regions of Interest in some brain atlas) and determining their connection strengths from some measure of time-series correlations. Although it is evident that the outcome must be affected by how the voxel-level time series are treated at the preprocessing stage, there is a lack of systematic studies of the effects of preprocessing on network structure. Here, we focus on the effects of spatial smoothing, a standard preprocessing method for fMRI. We apply various levels of spatial smoothing to resting-state fMRI data and measure the changes induced in functional networks. We show that the level of spatial smoothing clearly affects the degrees and other centrality measures of functional network nodes; these changes are non-uniform, systematic, and depend on the geometry of the brain. The composition of the largest connected network component is also affected in a way that artificially increases the similarity of the networks of different subjects. Our conclusion is that wherever possible, spatial smoothing should be avoided when preprocessing fMRI data for network analysis.

Functional brain segmentation using inter-subject correlation in fMRI

The human brain continuously processes massive amounts of rich sensory information. To better understand such highly complex brain processes, modern neuroimaging studies are increasingly utilizing experimental setups that better mimic daily-life situations. A new exploratory data-analysis approach, functional segmentation inter-subject correlation analysis (FuSeISC), was proposed to facilitate the analysis of functional magnetic resonance (fMRI) data sets collected in these experiments. The method provides a new type of functional segmentation of brain areas, not only characterizing areas that display similar processing across subjects but also areas in which processing across subjects is highly variable. FuSeISC was tested using fMRI data sets collected during traditional block-design stimuli (37 subjects) as well as naturalistic auditory narratives (19 subjects). The method identified spatially local and/or bilaterally symmetric clusters in several cortical areas, many of which are known to be processing the types of stimuli used in the experiments. The method is not only useful for spatial exploration of large fMRI data sets obtained using naturalistic stimuli, but also has other potential applications, such as generation of a functional brain atlases including both lower- and higher-order processing areas. Finally, as a part of FuSeISC, a criterion-based sparsification of the shared nearest-neighbor graph was proposed for detecting clusters in noisy data. In the tests with synthetic data, this technique was superior to well-known clustering methods, such as Ward's method, affinity propagation, and K-means ++. Hum Brain Mapp 38:2643-2665, 2017. © 2017 Wiley Periodicals, Inc.

Untangling the relatedness among correlations, part I: Nonparametric approaches to inter-subject correlation analysis at the group level

FMRI data acquisition under naturalistic and continuous stimuli (e.g., watching a video or listening to music) has become popular recently due to the fact that it entails less manipulation and more realistic/complex contexts involved in the task, compared to the conventional task-based experimental designs. The synchronization or response similarities among subjects are typically measured through inter-subject correlation (ISC) between any pair of subjects. At the group level, summarizing the collection of ISC values is complicated by their intercorrelations, which necessarily lead to the violation of independence assumed in typical parametric approaches such as Student's t-test. Nonparametric methods, such as bootstrapping and permutation testing, have previously been adopted for testing purposes by resampling the time series of each subject, but the quantitative validity of these specific approaches in terms of controllability of false positive rate (FPR) has never been explored before. Here we survey the methods of ISC group analysis that have been employed in the literature, and discuss the issues involved in those methods. We then propose less computationally intensive nonparametric methods that can be performed at the group level (for both one- and two-sample analyses), as compared to the popular method of circularly shifting the EPI time series at the individual level. As part of the new approaches, subject-wise (SW) resampling is adopted instead of element-wise (EW) resampling, so that exchangeability and independence assumptions are satisfied, and the patterned correlation structure among the ISC values can be more accurately captured. We examine the FPR controllability and power achievement of all the methods through simulations, as well as their performance when applied to a real experimental dataset.

How Many Is Enough? Effect of Sample Size in Inter-Subject Correlation Analysis of fMRI

Inter-subject correlation (ISC) is a widely used method for analyzing functional magnetic resonance imaging (fMRI) data acquired during naturalistic stimuli. A challenge in ISC analysis is to define the required sample size in the way that the results are reliable. We studied the effect of the sample size on the reliability of ISC analysis and additionally addressed the following question: How many subjects are needed for the ISC statistics to converge to the ISC statistics obtained using a large sample? The study was realized using a large block design data set of 130 subjects. We performed a split-half resampling based analysis repeatedly sampling two nonoverlapping subsets of 10–65 subjects and comparing the ISC maps between the independent subject sets. Our findings suggested that with 20 subjects, on average, the ISC statistics had converged close to a large sample ISC statistic with 130 subjects. However, the split-half reliability of unthresholded and thresholded ISC maps improved notably when the number of subjects was increased from 20 to 30 or more.