Effects of Tissue-Specific Functional Magnetic Resonance Imaging Signal Regression on Resting-State Functional Connectivity

Distinctive and Complementary Roles of Default Mode Network Subsystems in Semantic Cognition

The default mode network (DMN) typically deactivates to external tasks, yet supports semantic cognition. It comprises medial temporal (MT), core, and frontotemporal (FT) subsystems, but its functional organization is unclear: the requirement for perceptual coupling versus decoupling, input modality (visual/verbal), type of information (social/spatial), and control demands all potentially affect its recruitment. We examined the effect of these factors on activation and deactivation of DMN subsystems during semantic cognition, across four task-based human functional magnetic resonance imaging (fMRI) datasets, and localized these responses in whole-brain state space defined by gradients of intrinsic connectivity. FT showed activation consistent with a central role across domains, tasks, and modalities, although it was most responsive to abstract, verbal tasks; this subsystem uniquely showed more "tuned" states characterized by increases in both activation and deactivation when semantic retrieval demands were higher. MT also activated to both perceptually coupled (scenes) and decoupled (autobiographical memory) tasks and showed stronger responses to picture associations, consistent with a role in scene construction. Core DMN consistently showed deactivation, especially to externally oriented tasks. These diverse contributions of DMN subsystems to semantic cognition were related to their location on intrinsic connectivity gradients: activation was closer to the sensory-motor cortex than deactivation, particularly for FT and MT, while activation for core DMN was distant from both visual cortex and cognitive control. These results reveal distinctive yet complementary DMN responses: MT and FT support different memory-based representations that are accessed externally and internally, while deactivation in core DMN is associated with demanding, external semantic tasks.

Micapipe: A pipeline for multimodal neuroimaging and connectome analysis

Multimodal magnetic resonance imaging (MRI) has accelerated human neuroscience by fostering the analysis of brain microstructure, geometry, function, and connectivity across multiple scales and in living brains. The richness and complexity of multimodal neuroimaging, however, demands processing methods to integrate information across modalities and to consolidate findings across different spatial scales. Here, we present micapipe, an open processing pipeline for multimodal MRI datasets. Based on BIDS-conform input data, micapipe can generate i) structural connectomes derived from diffusion tractography, ii) functional connectomes derived from resting-state signal correlations, iii) geodesic distance matrices that quantify cortico-cortical proximity, and iv) microstructural profile covariance matrices that assess inter-regional similarity in cortical myelin proxies. The above matrices can be automatically generated across established 18 cortical parcellations (100-1000 parcels), in addition to subcortical and cerebellar parcellations, allowing researchers to replicate findings easily across different spatial scales. Results are represented on three different surface spaces (native, conte69, fsaverage5), and outputs are BIDS-conform. Processed outputs can be quality controlled at the individual and group level. micapipe was tested on several datasets and is available at https://github.com/MICA-MNI/micapipe, documented at https://micapipe.readthedocs.io/, and containerized as a BIDS App http://bids-apps.neuroimaging.io/apps/. We hope that micapipe will foster robust and integrative studies of human brain microstructure, morphology, function, cand connectivity.

NeoRS: A Neonatal Resting State fMRI Data Preprocessing Pipeline

Resting state functional MRI (rsfMRI) has been shown to be a promising tool to study intrinsic brain functional connectivity and assess its integrity in cerebral development. In neonates, where functional MRI is limited to very few paradigms, rsfMRI was shown to be a relevant tool to explore regional interactions of brain networks. However, to identify the resting state networks, data needs to be carefully processed to reduce artifacts compromising the interpretation of results. Because of the non-collaborative nature of the neonates, the differences in brain size and the reversed contrast compared to adults due to myelination, neonates can’t be processed with the existing adult pipelines, as they are not adapted. Therefore, we developed NeoRS, a rsfMRI pipeline for neonates. The pipeline relies on popular neuroimaging tools (FSL, AFNI, and SPM) and is optimized for the neonatal brain. The main processing steps include image registration to an atlas, skull stripping, tissue segmentation, slice timing and head motion correction and regression of confounds which compromise functional data interpretation. To address the specificity of neonatal brain imaging, particular attention was given to registration including neonatal atlas type and parameters, such as brain size variations, and contrast differences compared to adults. Furthermore, head motion was scrutinized, and motion management optimized, as it is a major issue when processing neonatal rsfMRI data. The pipeline includes quality control using visual assessment checkpoints. To assess the effectiveness of NeoRS processing steps we used the neonatal data from the Baby Connectome Project dataset including a total of 10 neonates. NeoRS was designed to work on both multi-band and single-band acquisitions and is applicable on smaller datasets. NeoRS also includes popular functional connectivity analysis features such as seed-to-seed or seed-to-voxel correlations. Language, default mode, dorsal attention, visual, ventral attention, motor and fronto-parietal networks were evaluated. Topology found the different analyzed networks were in agreement with previously published studies in the neonate. NeoRS is coded in Matlab and allows parallel computing to reduce computational times; it is open-source and available on GitHub (https://github.com/venguix/NeoRS). NeoRS allows robust image processing of the neonatal rsfMRI data that can be readily customized to different datasets.

Tonic pain alters functional connectivity of the descending pain modulatory network involving amygdala, periaqueductal gray, parabrachial nucleus and anterior cingulate cortex

INTRODUCTION

Resting state functional connectivity (FC) is widely used to assess functional brain alterations in patients with chronic pain. However, reports of FC accompanying tonic pain in pain-free persons are rare. A network we term the Descending Pain Modulatory Network (DPMN) is implicated in healthy and pathologic pain modulation. Here, we evaluate the effect of tonic pain on FC of specific nodes of this network: anterior cingulate cortex (ACC), amygdala (AMYG), periaqueductal gray (PAG), and parabrachial nuclei (PBN).

METHODS

In 50 pain-free participants (30F), we induced tonic pain using a capsaicin-heat pain model. functional MRI measured resting BOLD signal during pain-free rest with a 32°C thermode and then tonic pain where participants experienced a previously warm temperature combined with capsaicin. We evaluated FC from ACC, AMYG, PAG, and PBN with correlation of self-report pain intensity during both states. We hypothesized tonic pain would diminish FC dyads within the DPMN.

RESULTS

Of all hypothesized FC dyads, only PAG and subgenual ACC was weakly altered during pain (F=3.34; p=0.074; pain-free>pain d=0.25). After pain induction sACC-PAG FC became positively correlated with pain intensity (R=0.38; t=2.81; p=0.007). Right PBN-PAG FC during pain-free rest positively correlated with subsequently experienced pain (R=0.44; t=3.43; p=0.001). During pain, this connection's FC was diminished (paired t=-3.17; p=0.0026). In whole-brain analyses, during pain-free rest, FC between left AMYG and right superior parietal lobule and caudate nucleus were positively correlated with subsequent pain. During pain, FC between left AMYG and right inferior temporal gyrus negatively correlated with pain. Subsequent pain positively correlated with right AMYG FC with right claustrum; right primary visual cortex and right temporo-occipitoparietal junction Conclusion: We demonstrate sACC-PAG tonic pain FC positively correlates with experienced pain and resting right PBN-PAG FC correlates with subsequent pain and is diminished during tonic pain. Finally, we reveal PAG- and right AMYG-anchored networks which correlate with subsequently experienced pain intensity. Our findings suggest specific connectivity patterns within the DPMN at rest are associated with subsequently experienced pain and modulated by tonic pain. These nodes and their functional modulation may reveal new therapeutic targets for neuromodulation or biomarkers to guide interventions.

Structural Connectivity Gradients of the Temporal Lobe Serve as Multiscale Axes of Brain Organization and Cortical Evolution

The temporal lobe is implicated in higher cognitive processes and is one of the regions that underwent substantial reorganization during primate evolution. Its functions are instantiated, in part, by the complex layout of its structural connections. Here, we identified low-dimensional representations of structural connectivity variations in human temporal cortex and explored their microstructural underpinnings and associations to macroscale function. We identified three eigenmodes which described gradients in structural connectivity. These gradients reflected inter-regional variations in cortical microstructure derived from quantitative magnetic resonance imaging and postmortem histology. Gradient-informed models accurately predicted macroscale measures of temporal lobe function. Furthermore, the identified gradients aligned closely with established measures of functional reconfiguration and areal expansion between macaques and humans, highlighting their potential role in shaping temporal lobe function throughout primate evolution. Findings were replicated in several datasets. Our results provide robust evidence for three axes of structural connectivity in human temporal cortex with consistent microstructural underpinnings and contributions to large-scale brain network function.

Abnormal medial prefrontal cortex functional connectivity and its association with clinical symptoms in chronic low back pain

Accumulating evidence has shown that complicated brain systems are involved in the development and maintenance of chronic low back pain (cLBP), but the association between brain functional changes and clinical outcomes remains unclear. Here, we used resting-state functional magnetic resonance imaging (fMRI) and multivariate pattern analysis to identify abnormal functional connectivity (FC) between the default mode, sensorimotor, salience, and central executive brain networks in cLBP and tested whether abnormal FCs are related to pain and comorbid symptoms. Fifty cLBP patients and 44 matched healthy controls (HCs) underwent an fMRI scan, from which brain networks were identified by independent component analysis. Multivariate pattern analysis, graph theory approaches, and correlation analyses were applied to find abnormal FCs that were associated with clinical symptoms. Findings were validated on a second cohort of 30 cLBP patients and 30 matched HCs. Results showed that the medial prefrontal cortex/rostral anterior cingulate cortex had abnormal FCs with brain regions within the default mode network and with other brain networks in cLBP patients. These altered FCs were also correlated with pain duration, pain severity, and pain interference. Finally, we found that resting-state FC could discriminate cLBP patients from HCs with 91% accuracy in the first cohort and 78% accuracy in the validation cohort. Our findings suggest that the medial prefrontal cortex/rostral anterior cingulate cortex may be an important hub for linking the default mode network with the other 3 networks in cLBP patients. Elucidating the altered FCs and their association with clinical outcomes will enhance our understanding of the pathophysiology of cLBP and may facilitate the development of pain management approaches.

A decade of test-retest reliability of functional connectivity: A systematic review and meta-analysis

BACKGROUND

Once considered mere noise, fMRI-based functional connectivity has become a major neuroscience tool in part due to early studies demonstrating its reliability. These fundamental studies revealed only the tip of the iceberg; over the past decade, many test-retest reliability studies have continued to add nuance to our understanding of this complex topic. A summary of these diverse and at times contradictory perspectives is needed.

OBJECTIVES

We aimed to summarize the existing knowledge regarding test-retest reliability of functional connectivity at the most basic unit of analysis: the individual edge level. This entailed (1) a meta-analytic estimate of reliability and (2) a review of factors influencing reliability.

METHODS

A search of Scopus was conducted to identify studies that estimated edge-level test-retest reliability. To facilitate comparisons across studies, eligibility was restricted to studies measuring reliability via the intraclass correlation coefficient (ICC). The meta-analysis included a random effects pooled estimate of mean edge-level ICC, with studies nested within datasets. The review included a narrative summary of factors influencing edge-level ICC.

RESULTS

From an initial pool of 212 studies, 44 studies were identified for the qualitative review and 25 studies for quantitative meta-analysis. On average, individual edges exhibited a "poor" ICC of 0.29 (95% CI = 0.23 to 0.36). The most reliable measurements tended to involve: (1) stronger, within-network, cortical edges, (2) eyes open, awake, and active recordings, (3) more within-subject data, (4) shorter test-retest intervals, (5) no artifact correction (likely due in part to reliable artifact), and (6) full correlation-based connectivity with shrinkage.

CONCLUSION

This study represents the first meta-analysis and systematic review investigating test-retest reliability of edge-level functional connectivity. Key findings suggest there is room for improvement, but care should be taken to avoid promoting reliability at the expense of validity. By pooling existing knowledge regarding this key facet of accuracy, this study supports broader efforts to improve inferences in the field.

Microstructural and functional gradients are increasingly dissociated in transmodal cortices

While the role of cortical microstructure in organising neural function is well established, it remains unclear how structural constraints can give rise to more flexible elements of cognition. While nonhuman primate research has demonstrated a close structure-function correspondence, the relationship between microstructure and function remains poorly understood in humans, in part because of the reliance on post mortem analyses, which cannot be directly related to functional data. To overcome this barrier, we developed a novel approach to model the similarity of microstructural profiles sampled in the direction of cortical columns. Our approach was initially formulated based on an ultra-high-resolution 3D histological reconstruction of an entire human brain and then translated to myelin-sensitive magnetic resonance imaging (MRI) data in a large cohort of healthy adults. This novel method identified a system-level gradient of microstructural differentiation traversing from primary sensory to limbic regions that followed shifts in laminar differentiation and cytoarchitectural complexity. Importantly, while microstructural and functional gradients described a similar hierarchy, they became increasingly dissociated in transmodal default mode and fronto-parietal networks. Meta-analytic decoding of these topographic dissociations highlighted involvement in higher-level aspects of cognition, such as cognitive control and social cognition. Our findings demonstrate a relative decoupling of macroscale functional from microstructural gradients in transmodal regions, which likely contributes to the flexible role these regions play in human cognition.

Impact of Global Mean Normalization on Regional Glucose Metabolism in the Human Brain

Because the human brain consumes a disproportionate fraction of the resting body's energy, positron emission tomography (PET) measurements of absolute glucose metabolism (CMR_glc) can serve as disease biomarkers. Global mean normalization (GMN) of PET data reveals disease-based differences from healthy individuals as fractional changes across regions relative to a global mean. To assess the impact of GMN applied to metabolic data, we compared CMR_glc with and without GMN in healthy awake volunteers with eyes closed (i.e., control) against specific physiological/clinical states, including healthy/awake with eyes open, healthy/awake but congenitally blind, healthy/sedated with anesthetics, and patients with disorders of consciousness. Without GMN, global CMR_glc alterations compared to control were detected in all conditions except in congenitally blind where regional CMR_glc variations were detected in the visual cortex. However, GMN introduced regional and bidirectional CMR_glc changes at smaller fractions of the quantitative delocalized changes. While global information was lost with GMN, the quantitative approach (i.e., a validated method for quantitative baseline metabolic activity without GMN) not only preserved global CMR_glc alterations induced by opening eyes, sedation, and varying consciousness but also detected regional CMR_glc variations in the congenitally blind. These results caution the use of GMN upon PET-measured CMR_glc data in health and disease.

Amygdala hyper-connectivity in a mouse model of unpredictable early life stress

Childhood maltreatment is associated with a wide range of psychopathologies including anxiety that emerge in childhood and in many cases persist in adulthood. Increased amygdala activation in response to threat and abnormal amygdala connectivity with frontolimbic brain regions, such as the hippocampus and the prefrontal cortex, are some of the most consistent findings seen in individuals exposed to childhood maltreatment. The underlying mechanisms responsible for these changes are difficult to study in humans but can be elucidated using animal models of early-life stress. Such studies are especially powerful in the mouse where precise control of the genetic background and the stress paradigm can be coupled with resting-state fMRI (rsfMRI) to map abnormal connectivity in circuits that regulate anxiety. To address this issue we first compared the effects of two models of early-life stress, limited bedding (LB) and unpredictable postnatal stress (UPS), on anxiety-like behavior in juvenile and adult mice. We found that UPS, but not LB, causes a robust increase in anxiety in juvenile and adult male mice. Next, we used rsfMRI to compare frontolimbic connectivity in control and UPS adult male mice. We found increased amygdala-prefrontal cortex and amygdala-hippocampus connectivity in UPS. The strength of the amygdala-hippocampal and amygdala-prefrontal cortex connectivity was highly correlated with anxiety-like behavior in the open-field test and elevated plus maze. These findings are the first to link hyperconnectivity in frontolimbic circuits and increased anxiety in a mouse model of early-life stress, allowing for more mechanistic understanding of parallel findings in humans.

Spontaneous activity forms a foundation for odor-evoked activation maps in the rat olfactory bulb

Fluctuations in spontaneous activity have been observed by many neuroimaging techniques, but because these resting-state changes are not evoked by stimuli, it is difficult to determine how they relate to task-evoked activations. We conducted multi-modal neuroimaging scans of the rat olfactory bulb, both with and without odor, to examine interaction between spontaneous and evoked activities. Independent component analysis of spontaneous fluctuations revealed resting-state networks, and odor-evoked changes revealed activation maps. We constructed simulated activation maps using resting-state networks that were highly correlated to evoked activation maps. Simulated activation maps derived by intrinsic optical signal (IOS), which covers the dorsal portion of the glomerular sheet, significantly differentiated one odor's evoked activation map from the other two. To test the hypothesis that spontaneous activity of the entire glomerular sheet is relevant for representing odor-evoked activations, we used functional magnetic resonance imaging (fMRI) to map the entire glomerular sheet. In contrast to the IOS results, the fMRI-derived simulated activation maps significantly differentiated all three odors' evoked activation maps. Importantly, no evoked activation maps could be significantly differentiated using simulated activation maps produced using phase-randomized resting-state networks. Given that some highly organized resting-state networks did not correlate with any odors' evoked activation maps, we posit that these resting-state networks may characterize evoked activation maps associated with odors not studied. These results emphasize that fluctuations in spontaneous activity form a foundation for active processing, signifying the relevance of resting-state mapping to functional neuroimaging.

Detection of synchronous brain activity in white matter tracts at rest and under functional loading

Functional MRI based on blood oxygenation level-dependent (BOLD) contrast is well established as a neuroimaging technique for detecting neural activity in the cortex of the human brain. While detection and characterization of BOLD signals, as well as their electrophysiological and hemodynamic/metabolic origins, have been extensively studied in gray matter (GM), the detection and interpretation of BOLD signals in white matter (WM) remain controversial. We have previously observed that BOLD signals in a resting state reveal structure-specific anisotropic temporal correlations in WM and that external stimuli alter these correlations and permit visualization of task-specific fiber pathways, suggesting variations in WM BOLD signals are related to neural activity. In this study, we provide further strong evidence that BOLD signals in WM reflect neural activities both in a resting state and under functional loading. We demonstrate that BOLD signal waveforms in stimulus-relevant WM pathways are synchronous with the applied stimuli but with various degrees of time delay and that signals in WM pathways exhibit clear task specificity. Furthermore, resting-state signal fluctuations in WM tracts show significant correlations with specific parcellated GM volumes. These observations support the notion that neural activities are encoded in WM circuits similarly to cortical responses.